demo-05

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

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

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

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

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

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

# mat 파일로 저장
filename = "demo-05.mat";
save(filename, "x", "y1", "y2", "-7"); # save variables to a .mat file

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

# mat 파일에서 읽기
data = load(filename); # load variables from the .mat file

# fieldnames(data)  # 구조체 data의 필드 이름을 확인
# {
#   [1,1] = x
#   [2,1] = y1
#   [3,1] = y2
# }

x = data.x;
y1 = data.y1;
y2 = data.y2;
clear("data");

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

# 그래프
fig = 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");