다음 차분방정식에 대해

a)  에서 임펄스 응답을 구해보자.

ex2_11.cem

a=[1,-1,0.9];b=1;
% Part a) 
 x = implseq(0,-20,120); n=[-20:120];
 h = filter(b,a,x);
 subplot(2,1,1);stem(n,h,"-");
 title("Impulse Response"); xtitle("n"); ytitle("h(n)");

b)  에서 단위 계산 응답을 구해보자.

ex2_11.cem

% Part b) 
 x = stepseq(0,-20,120);
 s = filter(b,a,x);
 subplot(2,1,2); stem(n,s, "-") title("Step Response"); xtitle("n"); ytitle("s(n)");

c) 으로 특징된 이 시스템을 안정할까?

ex2_11.cem

% Part c)
sum(abs(h))z=roots(a);
magz=abs(z)