Octave Encyclopedia

main-01

manipulation of polinomial and graphical representation

ex-01/main-01.m

코드 인덱스로 돌아가기

카테고리

Exercise 01

학습 소스 그룹

코드 길이

60

lines

작성자

won sunggyu

2025-03-18

패키지

none

pkg load 기준

전체 코드

전체 코드를 복사해서 Octave에서 바로 실행할 수 있습니다.

# 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();

코드 해설

목적

manipulation of polinomial and graphical representation

입력

  • 스크립트 상단에서 정의한 파라미터/입력 데이터를 사용합니다.

출력

  • 콘솔 텍스트 출력

실행 흐름

  1. Visualization

핵심 함수/주제

plotdfmtttaxprintftextdfiguredmax

실습 과제

  • 질량/감쇠/강성 또는 전달함수 계수를 바꿔 응답 변화를 확인해보세요.
  • 축 범위와 라벨을 바꿔 그래프 해석성이 어떻게 달라지는지 확인해보세요.
  • 핵심 함수 plotd의 인자를 한 가지 바꿔 결과 변화를 기록해보세요.

학습 팁

같은 카테고리의 다른 코드