[Octave] HW#1

# 옥타브 Data Plot 숙제

# 1) 다항식 정의
roots = [0,1,3,5];    # 4차다항식 근
poly_coef = poly(roots);


# 2,3) 정의역 및 공역 설정, 다항식 그리기
x_range = [-1,6];
y_range = [-20,100];
xdata_interval = 0.1;

xlim(x_range);
ylim(y_range);

x = x_range(1) : xdata_interval : x_range(2);   # x 데이터 정의
y = polyval(poly_coef, x);  # 다항식의 x에 매칭되는 y값

figure(1);
plot(x, y);
title('Original Func','Fontsize',14);

# 4) 함수의 미분 그리기
der_coef = polyder(poly_coef);
y_der = polyval(der_coef, x);
figure(2);
plot(x, y_der);
title('Polyder','Fontsize',14);

# 5) 차분 정의에 의한 차분곡선
y_shifted = [y(2:length(y)), 0];    # 기존 y값을 1칸 옮긴 후, 0추가
y_hand_diff = (y_shifted-y)(1:length(y_shifted-y)-1);    # 끝 부분 데이터 제거
figure(3);
plot(x(1:length(x)-1), y_hand_diff(1:length(y_hand_diff))/xdata_interval);
title('Hand Diff','Fontsize',14);

# 6) 차분 함수를 사용한 차분곡선
y_diff = diff(y);   # 데이터 소실(x 길이 감소)
figure(4);
plot(x(1:length(x)-1), y_diff/xdata_interval);
title('Diff Func','Fontsize',14);

# 7) 한번에 그리기
figure(5);
set(gcf, "Name", "Total");
hold on;
plot(x,y,';Original Func;');
plot(x,y_der,';Polyder;');
plot(x(1:length(x)-1),y_hand_diff(1:length(y_hand_diff))/xdata_interval, ';Hand Diff;')
plot(x(1:length(x)-1),y_diff/xdata_interval, ';Diff Func;');
title('Total','Fontsize',14);
hold off;
grid on;