%《數字信號處理教程——MATLAB釋義與實現》第七章例7.5.2程序hc752a
% （逐步逼近）用等波動法設計低通濾波器2
% 電子工業出版社出版  陳懷琛編著 2004年9月
%
N=11;L=floor((N-1)/2);              % 濾波器及方程階數
w0=[0, 0.65, 1, 1.5, 1.8, 2.4,  3.14]; % 初設頻率向量
D=[1,1,1,0,0,0,0];                  % 預期幅頻特性
P=cos(w0'*[0:L]);                   % 符幅特性諧波向量
Wt=[0.5,0.5,0.5,1,1,1,1];           % 權重函數
WT=(-1).^[0:6]./Wt;                 % 求系數中的WT項
Pa=[P,WT'];                          % 最小最大方程系數向量
d=Pa\ D'                            % 求解最小最大方程
h=[flipud(d(2:6)/2)',d(1),d(2:6)'/2]  % 按類型I的FIR濾波器由d求h
[H,w]=freqz(h,1);				% 由脈沖響應求頻率特性
plot(w,abs(H)),hold on					% 畫出實際的幅頻特性
H0=freqz(h,1,w0);abs(H0)-D		% 求給定頻點上的幅特性，并求點上的誤差
plot(w0,abs(H0),'.')   		% 畫出設定頻點上的幅特性