demo-01

# filename: demo-01.m
# writer: won sunggyu
# date: 2025-04-22
# language: octave
# description: 배열 생성, 기본 연산, 기본 수학 함수 익히기

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

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

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

dx = 0.01;
nx = 1000;
x = 0:dx:dx*(nx-1); # x축의 범위
pp = poly([1, 4, 7]); # 다항식

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

r_pow1 = x .^ 2; # 지수 연산자
r_pow2 = power(x, 0.5); # 지수 함수

r_sin = sin(x); # 삼각 함수
r_cos = cos(x); # 삼각 함수
r_exp = exp(x); # 지수 함수
r_ply = polyval(pp, x); # 다항 함수

r_add = r_sin + r_cos; # 배열 연산 (element-wise)
r_sub = r_sin - r_cos; # 배열 연산 (element-wise)
r_mul = r_sin .* r_ply; # 배열 연산 (element-wise)
r_div = r_sin ./ r_exp; # 배열 연산 (element-wise)

r_log = log(x);
r_log10 = log10(x);
r_tanh = tanh(x);
r_coth = coth(x);

#------------------------------------------------------------------------------
# 그래프 그리기
#------------------------------------------------------------------------------
# graphics_toolkit("qt");  % 툴킷 변경
# graphics_toolkit("fltk");  % 툴킷 변경
graphics_toolkit("gnuplot");  % 툴킷 변경


# 그래프
fig = figured("Size", [1440, 960], "Move", [0, 0], "Name", mfilename);
# fig = figured("Size", [1440, 960], "Move", [0, 0], "Name", mfilename, "visible", "off");
ax = subplots(4, 3);

plot(ax(1, 1), x, r_sin);
plot(ax(2, 1), x, r_cos);
plot(ax(3, 1), x, r_exp);
plot(ax(4, 1), x, r_ply);

plot(ax(1, 2), x, r_add);
plot(ax(2, 2), x, r_sub);
plot(ax(3, 2), x, r_mul);
plot(ax(4, 2), x, r_div);

plot(ax(1, 3), x, r_log);
plot(ax(2, 3), x, r_log10);
plot(ax(3, 3), x, r_tanh);
plot(ax(4, 3), x, r_coth);



print("myplot.png", "-dpng");