模拟滤波器matlab程序,巴特沃兹模拟滤波器的MATLAB实现

巴特沃兹模拟滤波器的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)]);