www.pudn.com > (try2).rar > (try2).m

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