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