基于matlab的FIR低通,高通,带通,带阻滤波器设计

ning窗函数的单位脉冲响应’);Рfigure(4);Рplot(w/pi,20*log(abs(h2)/abs(h2(1))));Рgrid;Рxlabel('w/pi');Рylabel('幅度(dB)');Р4.带阻滤波器设计Р% 采用Hamming窗设计一个带阻FIR滤波器Р% 要求:Р% 阻带:0.35pi~0.65pi,带内最小衰减Rs=50dB;Р% 通带:0~0.2pi和0.8pi~pi,带内最大衰减:Rp=1dB。РclcРclearРWpl=0.2*pi; %根据阻带:0.35pi~0.65pi,通带:0~0.2pi和0.8pi~pi,РWph=0.8*pi; %确定两个通带截止频率和两个阻带截止频率。РWsl=0.35*pi;РWsh=0.65*pi;Рdert_w=min((Wsl-Wpl),(Wph-Wsh));РN=ceil(6.6*pi/dert_w); %根据过度带宽确定NРn=0:1:N-1;РWcl=(Wsl+Wpl)/2; %低通中心频率РWch=(Wsh+Wph)/2; %高通中心频率Рhd=ideal(pi,N)-ideal(Wch,N)+ideal(Wcl,N); %带通滤波器的原型高通+低通РB=hamming(N)'; %海明窗Рh=hd.*B; %加窗后Р[H,m]=freqz(h,[1],1024,'whole');%获取频率响应Рmag=abs(H); %幅值Рdb=20*log10((mag+eps)/max(mag));%分贝数Рpha=angle(H); %相位Рw=m/pi;Р%绘图Рsubplot(2,2,1);Рstem(n,hd);Рtitle('理想时域波形');Рxlabel('n');Рylabel('sa(n)');Рsubplot(2,2,2);Рstem(n,B);Рtitle('海明窗');