(try2).m

FIR数字滤波器的最优设计(基于MATLAB软件的设计)

Wp=input('请输入FIR滤波器的通带截止频率(单位：弧度)Wp=');
Ws=input('请输入FIR滤波器的阻带截止频率(单位：弧度)Ws=');
delta1=input('请输入FIR滤波器的通带偏差:');
delta2=input('请输入FIR滤波器的阻带偏差:');
Ap=20*log10(2/(1-delta1)-1)
As=-20*log10(delta2/(delta1+1))
delta_w=2*pi/1000;
Wsi=Ws/delta_w+1;
Wpi=Wp/delta_w;
f=[Wp/pi,Ws/pi];
m=[1,0];
dev=[delta1,delta2];
[n,fo,ao,weights]=remezord(f,m,dev)
h=remez(n,fo,ao,weights);
[db,mag,pha,w]=freqz_m(h,[1]);
Asd=-max(db(Wsi:1:501));
while Asd<As
n=n+1;
h=remez(n,fo,ao,weights);
[db,mag,pha,w]=freqz_m(h,[1]);
Asd=-max(db(Wsi:1:501));
end
M=n;
h=remez(n,fo,ao,weights);
[db,mag,pha,w]=freqz_m(h,[1]);
Apd=-(min(db(1:1:Wp/delta_w+1)));
Asd=-round(max(db(Wsi:1:501)));
subplot(1,1,1);
subplot(2,1,1);
stem([0:1:M],h);
title('实际的脉冲响应')
axis([0,M,-0.1,0.5]);
xlabel('n');ylabel('h(n)');
subplot(2,1,2);
plot(w/pi,db);
title('幅度响应');
axis([0 1 -100 10]);
grid;
xlabel('以pi为单位的频率');
ylabel('分贝数／dB');