%%模擬濾波器
%程序作者：郭源蕊?學號：?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（