다음의 차분 방정식으로 기술되는 인과적 LTI 시스템에 대해,

a) 시스템 함수( H(z) )를 구해보자.

b) 단위 임펄스 응답( h(n) )를 구해보자.

c) 단위 계단 응답( v(n) ), 즉 단위 계단 ( u(n) )에 대한 응답을 구해보자.

d) 주파수 응답 함수 H(e^{j \omega})와 크기, 위상을  구간에서 그려보자.

ex4_13.cem

/*   Chapter 4: Example 4.13 */
 
/*b.*/
 b = [1,0,-1]; a = [1,0,-0.81];
 [R,p,C] = resid_z(b,a) 
/*d.*/
 
 w = [0:500:1]*pi/500;
 H = freq_z(b,a,w);
 magH = abs(H); phaH = angle(H);
 subplot(2,1,1); plot(w/pi,magH);
 xtitle("frequency in pi units"); ytitle("Magnitude");
 title("Magnitude Response");
 subplot(2,1,2); plot(w/pi,phaH/pi);
 xtitle("frequency in pi units"); ytitle("Phase in pi units");
 title("Phase Response");