%-----------------------------------------------------------------
% exa020802, for example2.8.1 ~ 2.8.3  
% to test impz.m,filter and freqz
%-----------------------------------------------------------------
clear;
% 計算濾波器的沖擊響應
b=[.001836,.007344,.011016,.007374,.001836];
a=[1,-3.0544,3.8291,-2.2925,.55075];
[h,t]=impz(b,a,40);
subplot(221)
stem(t,h,'.');grid on;
ylabel('h(n)')
xlabel('n')

% 計算濾波器的階躍響應
x=ones(100);t=1:100;
y=filter(b,a,x);
subplot(222)
plot(t,x,'g.',t,y,'k-');grid on;
ylabel('x(n)  and  y(n)')
xlabel('n')

% 計算濾波器的幅頻響應和相頻響應(解卷繞表示)
[H,w]=freqz(b,a,256,'whole',1);      %求頻率響應；
Hr=abs(H);                           %求幅頻；
Hphase=angle(H);                     %求相頻；
Hphase=unwrap(Hphase);               %解卷繞；
subplot(223)
plot(w,Hr);grid on;
ylabel('幅頻響應')
xlabel('歸一化頻率')
subplot(224)
plot(w,Hphase);grid on;
ylabel('相頻響應')
xlabel('歸一化頻率')