main-01-01

# - 입력 데이터 준비
# - 데이터 범위 산정
# - 데이터 출력
# - 연산
# - 시각화
clc; clear all; close all;

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

tspan = [0.5, 3.5];
nn = 100;
dt = (tspan(2) - tspan(1)) / nn;
tt = (0:nn-1) * dt + tspan(1);
t_ = tt(2:end);

pp_val = polyval(pp, tt);
dp_val = polyval(dp, tt);

disp([min(tt) max(tt)]);
disp([min(pp_val) max(pp_val)]);
disp([min(dp_val) max(dp_val)]);

dp_v_0 = diff(pp_val) / dt;
dp_v_1 = (pp_val(2:end) - pp_val(1:end-1)) / dt;

figure
hold on
plot(tt, pp_val);
plot(tt, dp_val);
plot(t_, dp_v_0);
plot(t_, dp_v_1);
legend("pp", "dp", "dp_v_0", "dp_v_1");