demo-04

# filename: demo-04.m
# writer: won sunggyu
# date: 2025-04-22
# language: octave
# description: csv 파일 읽기 및 쓰기

#------------------------------------------------------------------------------
# 초기화
#------------------------------------------------------------------------------

run("startup.m");
printf(fmt("{mfilename}\n", "#FF5733"));

#------------------------------------------------------------------------------
# 데이터 준비
#------------------------------------------------------------------------------

# 1D 진동 응답식
dx = 0.01;
nx = 1000;
x = 0:dx:dx*(nx-1);

fn = 5; # frequency
ze = 0.01; # damping ratio
wn = 2*pi*fn; # natural frequency
wd = sqrt(1-ze^2)*wn; # damped frequency

r_sin = sin(wd*x);
r_cos = cos(wd*x);
r_exp = exp(-ze*wn*x);

y1 = r_sin .* r_exp;  # damped sine wave
y2 = r_cos .* r_exp;  # damped cosine wave

#------------------------------------------------------------------------------
# 데이터 연산
#------------------------------------------------------------------------------

# csv 파일로 저장
filename = "demo-04.csv";
csvwrite(filename, [x', y1', y2']);

# 변수 삭제
clear("x", "y1", "y2");

# csv 파일에서 읽기
data = csvread(filename);

x = data(:, 1);
y1 = data(:, 2);
y2 = data(:, 3);

#------------------------------------------------------------------------------
# 그래프 그리기
#------------------------------------------------------------------------------

# 그래프
figured("Size", [1280, 720], "Move", [0, 0], "Name", mfilename);
ax = subplots(2, 1, "Xlabel", "Time (s)", "Ylabel", "Amplitude");

plot(ax(1, 1), x, y1);
plot(ax(2, 1), x, y2);

title(ax(1, 1), "Damped Sine Wave");
title(ax(2, 1), "Damped Cosine Wave");