不同模拟滤波器设计matlab代码

%%模拟滤波器
%程序作者：郭源蕊 学号： 201614144052
%程序描述：利用matlab实现模拟滤波器的设计，低通、低通到高通、低通到
%带通，低通到带阻的滤波器设计和波形比较。

%巴特沃斯——低通滤波器设计例题6.2.1
wp=2*pi*5000;ws=2*pi*12000;Rp=2;As=30;
[Nwc]=buttord（wpwsRpAs‘s‘）%计算率波器的阶数和3dB截止频率
[BA]=butter（Nwc‘s‘）;%计算滤波器系统函数分子分母多项式
fk=0:800/512:20000;wk=2*pi*fk;
Hk=freqs（BAwk）;
%figure
subplot（421）
subplot（421）
plot（fk/100020*log10（abs（Hk）））;
grid onxlabel（‘频率（kHz）‘）ylabel（‘幅度（dB）‘）
title（‘巴特沃斯低通模拟滤波器‘）
axis（[020-355]）

%巴特沃斯——高通滤波器设计例题6.2.6
wp=1;ws=4;Rp=0.1;As=40;
[Nwc]=buttord（wpwsRpAs‘s‘）%计算率波器的阶数和3dB截止频率
[BA]=butter（Nwc‘s‘）;%计算滤波器系统函数分子分母多项式
wph=2*pi*4000;
[BHAH]=lp2hp（BAwph）;
fk=0:800/512:20000;wk=2*pi*fk;
Hk=freqs（BHAHwk）;
Hk1=freqs（BAwk）;
%figure
subplot（423）;
plot（fk/100020*log10（abs（Hk1）））;
grid onxlabel（‘频率（kHz）‘）ylabel（‘幅度（dB）‘）
title（‘巴特沃斯低通模拟滤波器（利用低通转高通）‘）
subplot（424）
plot（fk/100020*log10（abs（Hk）））;
grid onxlabel（‘频率（kHz）‘）ylabel（‘幅度（dB）‘）
title（‘巴特沃斯高通模拟滤波器（利用低通转高通）‘）
axis（[010-355]）

%巴特沃斯——带通滤波器设计例题6.2.7
wp=1;ws=1.963;Rp=1;As=20;
[Nwc]=buttord（wpwsRpAs‘s‘）%计算率波器的阶数和3dB截止频率
[BA]=butter（Nwc‘s‘）;%计算滤波器系统函数分子分母多项式
wph=2*pi*4000;
fk=0:800/512:20000;wk=2*pi*fk;
Hk1=freqs（