ms24_6.m

《数字信号处理实验(MATLAB版)》这

%FIR高通滤波器
fs=300;fp=450;Fs=2000;  %输入设计指标
ws=fs/(Fs/2)*pi;wp=fp/(Fs/2)*pi; %计算归一化角频率
deltaw=wp-ws;           %计算过渡带的宽度
N0=ceil(6.1*pi/deltaw); %按三角形窗计算滤波器长度N0
N=N0+mod(N0+1,2)        %为实现FIR类型I偶对称滤波器,应确保N为奇数
windows=triang(N);      %使用三角形窗
wc=(ws+wp)/2/pi;        %截止频率取归一化通阻带频率的平均值
b=fir1(N-1,wc,'high',windows);
[db,mag,pha,grd,w]=freqz_m(b,1); %求解频率特性
n=0:N-1;dw=2*pi/1000;   %dw为频率分辨率，将0～2π分为1000份
Rp=-(min(db(wp/dw+1:501)))         %检验通带波动
As=-round(max(db(1:ws/dw+1)))  %检验最小阻带衰减
%
subplot(2,2,1),stem(n,b,'filled','k');
axis([0,N,1.1*min(b),1.1*max(b)]);title('实际脉冲响应');
xlabel('n');ylabel('h(n)');
subplot(2,2,2),stem(n,windows,'filled','k');
axis([0,N,0,1.1]);title('窗函数特性');
xlabel('n');ylabel('wd(n)');
subplot(2,2,3),plot(w/2/pi*Fs,db,'k');
axis([0,Fs/2,-40,2]);title('幅频响应');
xlabel('f(Hz)');ylabel('H(e^{j\omega})');
set(gca,'XTickMode','manual','XTick',[0,fs,fp,Fs/2]);
set(gca,'YTickMode','manual','YTick',[-20,-3,0]);grid
subplot(2,2,4),plot(w/2/pi*Fs,pha,'k');
axis([0,Fs/2,-4,4]);
title('相频响应');
xlabel('f(Hz)');ylabel('\phi(\omega)');
set(gca,'XTickMode','manual','XTick',[0,fs,fp,Fs/2]);
set(gca,'YTickMode','manual','YTick',[-pi,0,pi]);grid