main_chan_01

clc
clear

#1
A = [5 8 1;2 4 1;3 4 8];
P = poly(A);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

#2
x_lim = linspace(-10, 10, 100);
y_lim = [-3000:1000];

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

#3
f = polyval(P,x_lim);

figure(1)
plot(x_lim,f)
grid on
ylim([-3000 1000]);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

#4
k = polyder(P)
f_p = polyval(k, x_lim)

figure(2)
plot(x_lim, f_p)
grid on
ylim([-3000 1000]);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

#5
for i = 1:length(x_lim)-1
m_di_P(i) =  f(i+1) - f(i);
end

figure(3)
plot(x_lim(2:end), m_di_P)
grid on
ylim([-3000 1000]);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

#6
di_P = diff(f)

figure(4)
plot(x_lim(2:end), di_P)
grid on
ylim([-3000 1000]);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

#7
figure(5)
plot(x_lim,f); hold on
plot(x_lim, f_p);
plot(x_lim(2:end), m_di_P , 'o');
plot(x_lim(2:end), di_P);
hold off
grid on
ylim([-3000 1000]);