matlab巴特沃兹滤波器,巴特沃兹模拟滤波器的MATLAB实现.docx

巴特沃兹模拟滤波器的MATLAB实现

设计一个模拟巴特沃特低通滤波器，它在30rad/s处具有1dB或更好的波动，在50rad/s处具有至少30dB的衰减。求出级联形式的系统函数，画出滤波器的幅度响应、对数幅度响应、相位响应和脉冲响应图。

解：程序清单如下

Wp=30;Ws=50;Rp=1;As=30; %技术指标

Ripple=10^(-Rp/20);

Attn=10^(-As/20);

[b,a]=afd_butt(Wp,Ws,Rp,As) %巴特沃兹低通滤波器

[C,B,A]=sdir2cas(b,a) %计算二阶节系数,级联型实现

[db,mag,pha,w]=freqs_m(b,a,50); %计算幅频响应

[ha,x,t]=impulse(b,a); %计算模拟滤波器的单位脉冲响应

figure(1);clf;

subplot(2,2,1);plot(w,mag);title('Magnitude Response');

xlabel('Analog frequency in rad/s');

ylabel('H');

%

axis([0,50,0,1.1])

set(gca,'XTickMode','manual','XTick',[0,30,40,50]);

set(gca,'YTickMode','manual','YTick',[0,Attn,Ripple,1]);

grid

subplot(2,2,2);plot(w,db);title('Magnitude in dB');

xlabel('Analog frequency in rad/s');

ylabel('decibels');

%

axis([0,50,-40,5])

set(gca,'XTickMode','manual','XTick',[0,30,40,50]);

set(gca,'YTickMode','manual','YTick',[-40,-As,-Rp,0]);

grid

subplot(2,2,3);plot(w,pha/pi);title('Phase Response');

xlabel('Analog frequency in rad/s');

ylabel('radians');

%

axis([0,50,-1.1,1.1])

set(gca,'XTickMode','manual','XTick',[0,30,40,50]);

set(gca,'YTickMode','manual','YTick',[-1,-0.5,0,0.5,1]);

grid

subplot(2,2,4);plot(t,ha);title('Impulse Response');

xlabel('time in seconds');

ylabel('ha(t)');

axis([0,max(t)+0.05,min(ha),max(ha)+0.025]);

set(gca,'XTickMode','manual','XTick',[0,0.1,max(t)]);

set(gca,'YTickMode','manual','YTick',[0,0.1,max(ha)]);

grid

%巴特沃兹模拟滤波器的设计子程序

function[b,a]=afd_butt(Wp,Ws,Rp,As);

if Wp<=0

error('Passband edge must be larger than 0')

end

if Ws<=Wp

error('Stopband edge must be larger than Passed edge')

end

if (Rp<=0)|(As<0)

error('PB ripple and /0r SB attenuation must be larger than 0')

end

N=ceil((log10((10^(Rp/10)-1)/(10^(As/10)-1)))/(2*log10(Wp/Ws)));

OmegaC=Wp/((10^(Rp/10)-1)^(1/(2*N)));

[b,a]=u_buttap(N,OmegaC);

%设计非归一化巴特沃兹模拟低通滤波器原型子程序

function [b,a]=u_buttap(N,OmegaC);

[z,p,k]=buttap(N);

p=p*OmegaC;

k=k*OmegaC^N;

B=real(poly(z));

b0=k;

b=k*B;

a=real(poly(p));

%计算系统函数的幅度响应和相位响应子程序

fun