%《數字信號處理教程——MATLAB釋義與實現》第七章例7.2.2程序hc722
% 線性相位FIR濾波器的類型及符幅特性
% 電子工業出版社出版  陳懷琛編著 2004年9月
%
clear
h1=[3,-1,2,-3,5,-3,2,-1,3];             % 輸入一個類型I的線性相位濾波器系數
[A1,w1,typea,tao1]=amplres(h1);         % 用amplres函數求符幅特性并判別其類型
subplot(2,2,1),plot(w1,A1),grid on      % 畫出符幅特性
title('線性相位濾波器類型I的符幅特性')
h2= [3,-1,2,-3,-3,2,-1,3];              % 輸入一個類型II的線性相位濾波器系數
[A2,w2,typeb,tao2]=amplres(h2);         % 用amplres函數求符幅特性并判別其類型
subplot(2,2,2),plot(w2,A2),grid on      % 畫出符幅特性
title('線性相位濾波器類型II的符幅特性')
h3=[3,-1,2,-3,0, 3, -2, 1,-3];          % 輸入一個類型III的線性相位濾波器系數
[A3,w3,typec,tao3]=amplres(h3);         % 用amplres函數求符幅特性并判別其類型
subplot(2,2,3),plot(w3,A3),grid on      % 畫出符幅特性
title('線性相位濾波器類型III的符幅特性')
h4=[3,-1,2,-3,03, -2, 1,-3];            % 輸入一個類型IV的線性相位濾波器系數
[A4,w4,typed,tao4]=amplres(h4);         % 用amplres函數求符幅特性并判別其類型
subplot(2,2,4),plot(w4,A4),grid on      % 畫出符幅特性
title('線性相位濾波器類型IV的符幅特性')
set(gcf,'color','w');                   % 置圖形背景色為白