(pi*m); %理想脉冲响应 function [db,mag,pha,W]=freqz_m(b,a) ;%自定义函数%db 为相对振幅(dB) ; mag 为绝对振幅; pha 为相位响应; W频率样本点向量。[H,W]=freqz(b,a,1000, 'whole'); %b 和 a分别为 H(z) 分子和分母多项式系数; H=(H(1:501))'; W=(W(1:501))'; mag=abs(H); db=20*log10((mag+eps)/max(mag)); pha=angle(H); %。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。 Wp=2* 3000 /7000; Ap= 1;%滤波器性能指标 Wr=2* 2400 /7000; Ar= 50; Width= Wp-Wr; %过渡带宽度 Wc= (Wr+Wp)/2; %理想滤波器的截止频率 N1=ceil(6.6*pi/Width); %滤波器长度 N=N1+mod(N1+1,2); %确保 N为奇数 n=0:1:N-1; %步长 hd=ideallp(Wc,N); %理想单位脉冲响应 window=hamming(N)'; %海明窗 h=window.*hd; %截取得到实际单位脉冲响应[db,mag,pha,W]=freqz_m(h,1); %计算实际滤波器的幅频响应 subplot(2, 1 ,1);%分割窗口 stem(n,h); %绘出实际单位脉冲响应 title(' 实际脉冲响应 h(n)'); %加标题说明 grid; %添加网格 subplot(2, 1 ,2);%分割窗口 plot(W/pi,db); %绘出实际滤波器的幅频响应 title(' 幅度响应(db)'); %加标题说明 grid%添加网格
全文预览
