bài 2.6 tín hiệu và hệ thống

clear;
clf;
ym1=1;ym2=-2;
xm1=0;xm2=0;
b2=3;b1=-2;b0=1;
a1=0.9;a0=0.8;
x=(0:10)*0;y=x;
for n=0:10
x(n+1)=1;
yc=b2*x(n+1)+b1*xm1+b0*xm2-a1*ym1-a0*ym2;
ym2=ym1;ym1=yc;
xm2=xm1;xm1=x(n+1);
y(n+1)=yc;
end
n=0:1:10;
plot(n,y,'*',n,x,'o');
grid on;
ylabel('y(n)')
xlabel('n')
legend('output','input')
title('Response of a second-order system by iteration')