每周一题（五）滤波器设计

创

wp = 0.2*pi;                          %  数字滤波器指标
Rp = 1;
As  = 15;
Fs = 100;
%  分贝形式的指标转化为绝对坐标下的指标（绘图标刻度用）
Ripp le = 1/(10^(Rp/20));
Attn =  1/(10^(As/20));
%  模拟原型指标的频率逆映射
T  = 1/Fs;
OmegaP = (2/T)*tan(wp/2);
OmegaS = (2/T)*tan(ws/2);
%  模拟巴特沃思原型滤波器计算
disp('模拟巴特沃思原型滤波器阶数 N 和 3dB 截止频率 OmegaC 分别为： ')


22

[N,OmegaC]  =buttord(OmegaP,OmegaS,Rp,As,'s')
[z0,p0,k0] = buttap( N);                       %  归一化巴特沃思原型设计
ba0 = real(po ly( z0));ba0  = k0*ba0 ;           %  由零点计算分子系数归一化向量
aa0 = real(po ly(p0));                          %  由极点计算分母系数归一化向量
d isp('分子系数向量 bd，分母系数向量 ad 分别为：')
[bd,ad] =  bilinear(ba0,aa0,Fs/OmegaC)        %  双线性变换
d isp('二阶环节级联结构的分子分母系数矩阵 sos,整体增益 G 分别为：')
[sos,G] = tf2sos(bd, ad)                       %  变为二阶环节级联结构


[H,w] =  freqz(bd,ad,1000);                    %  计算数字系统频率响应
mag = abs( H);
db = 20*log10(( mag+eps) /max(mag));
pha = angle(H);
grd = grpdelay(bd,ad,w);
%  绘图
subplot(2,1,1);zplane(bd,ad) ; title( '零极点分布')              %绘制零极点分布


原
xlabel( ''); ylabel( '|H|');axis([0,1,0,1.1]);
set(gca,'XTickMode','manual','XTick',[0,wp/pi,ws/pi,1]);
set(gca,'YTick mode', 'manual','YTick',[ 0,Attn,Ripple,1]);grid


subplot(4,2,7); plot(w/pi,db); title('模值(dB)');               % 绘制幅度响应（分贝）
xlabel( '频率：（单位：pi）');  ylabel( '分贝'); axis([0,1,-20,5]);
set(gca,'XTickMode','manual','XTick',[0,wp/pi,ws/pi,1]);
set(gca,'YTick mode', 'manual','YTick',[ -20,-As,-Rp,0]); gr id



subplot(4,2,6); plot(w/pi,pha/p i); title('相位响应')             %绘制相位响应
xlabel( ''); ylabel( '单位:pi') ; axis([0,1,-1,1]);
set(gca,'XTickMode','manual','XTick',[0,wp/pi,ws/pi,1]);

23

set(gca,'YTick mo de', 'manual','YTick',[ -1,0,1]); grid


subplot(4,2,8); plot(w/pi,grd); title('群延迟')                 %绘制群延迟
xlabel( '频率（单位：pi）');  ylabel( '样本'); axis([0,1,0,8])
set(gca,'XTickMode','manual','XTick',[0,wp/pi,ws/pi,1]);
set(gca,'YTick mode', 'manual','YTick',[ 0:2:8]); grid


set(gcf,'co lor','w')                                           %设置绘图背景为白色

肯定有错误，求大神帮我改正，有不一样的最好了。