main-01

# filename: main-01.m
# writer: won sunggyu
# date: 2025-03-18
# language: octave
# description: manipulation of polinomial and graphical representation

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

pp = poly([1, 2, 3]);
dp = polyder(pp);

tspan = [0.8, 3.2]; 
nn = 10000;
[tt, dt] = timeaxis(tspan, nn); # 1xN
t_ = tt(2:end); # 1x(N-1)

pp_val = polyval(pp, tt); # 1xN
dp_val = polyval(dp, tt); # 1xN

# printf("tt     min: %6.2f, max: %6.2f\n", min(tt), max(tt));
# printf("pp_val min: %6.2f, max: %6.2f\n", min(pp_val), max(pp_val));
# printf("dp_val min: %6.2f, max: %6.2f\n", min(dp_val), max(dp_val));

dp_v_0 = diff(pp_val) / dt; # 1x(N-1)
dp_v_1 = my_diff(pp_val) / dt; # 1x(N-1)

# Visualization
param_f = {"Size", [960, 720], "Name", "ECO-5-Ex-01"};
param_a = {"Xlabel", "X-axis", "Ylabel", "Y-axis", "XLim", [0.5, 3.5], "YLim", [-2, 4]};

figured(param_f);
subplots(param_a);
plotd(tt, pp_val, ";pp;");
plotd(tt, dp_val, ";dp;");
textd(0.05, 0.90, fmt("pp = [ {pp} ]"));
textd(0.05, 0.85, fmt("dp = [ {dp} ]"));

# Animation
figured(param_f);
ax = subplots(2, 1, param_a);
plotd(ax(1), tt, pp_val, ";pp;");
plotd(ax(2), tt, dp_val, ";dp;");
textd(ax(1), 0.05, 0.90, fmt("pp = [ {pp} ]"));
textd(ax(2), 0.05, 0.85, fmt("dp = [ {dp} ]"));

figured(param_f);
subplots(param_a);
h1 = plotd(tt(1:2), pp_val(1:2), ";pp;");
h2 = plotd(tt(1:2), dp_val(1:2), ";dp;");
anim = Animator();
for i = 2:20:length(tt)
    set(h1, "Xdata", tt(1:i), "Ydata", pp_val(1:i));
    set(h2, "Xdata", tt(1:i), "Ydata", dp_val(1:i));
    drawnow;
    anim.capture();
end
anim.close();