fs=22050;%采樣頻率
[Yfs]=wavread（‘F:\學校相關\信號處理要求\Matlab文件\姜杰.wav‘）;%根據采樣頻率將音頻采樣到y中
sound（Y22050）;?%播放信號

%加入噪音
t=（0:length（y）-1）/fs;
d=[0.005*cos（2*pi*5500*t）+0.005*cos（2*pi*5500*t）]‘;
y=Y+d;

%設置性能指標
Ts=1/fs;R1=10;
wp=2*pi*1100/fs;
ws=2*pi*1400/fs;
Rp=1;
Rl=30;
wp1=2/Ts*tan（wp/2）;?%將模擬指標轉換成數字指標
ws1=2/Ts*tan（ws/2）;?%將模擬指標轉換成數字指標

%選用buttord函數將參數傳入
[NWn]=buttord（wp1ws1RpRl‘s‘）;?%選擇濾波器的最小階數
[ZPK]=buttap（N）;?%創建butterworth模擬濾波器
[BapAap]=zp2tf（ZPK）;
[ba]=lp2lp（BapAapWn）;
[bzaz]=bilinear（bafs）;?%用雙線性變換法實現模擬濾波器到數字濾波器的轉換
[HW]=freqz（bzaz）;?%繪制頻率響應曲線

%將圖畫出來
figure（4）
plot（W*fs/（2*pi）abs（H））%頻率響應
grid
xlabel（‘頻率／Hz‘）
ylabel（‘頻率響應幅度‘）
title（‘IIR低通濾波器‘）

%進行濾波并畫圖比較
f1=filter（bzazy）;
figure（5）
subplot（211）
plot（y）?%畫出濾波前的時域圖
title（‘IIR低通濾波器濾波前的時域波形‘）;
subplot（212）
plot（f1）;?%畫出濾波后的時域圖
title（‘IIR低通濾波器濾波后的時域波形‘）;
F0=fft（f11024）;
f=fs*（0:511）/1024;
figure（6）
y2=fft（y1024）;
subplot（211）;
plot（fabs（y2（1:512）））;?%畫出濾波前的頻譜圖
title（‘IIR低通濾波器濾波前的頻譜‘）
xlabel（‘頻率/Hz‘）;
ylabel（‘幅值‘）;
subplot（212）
F1=plot（fabs（F0（1:512）））;?%畫出濾波后的頻譜圖
title（‘IIR低通濾波器濾波后的頻譜‘）?
xlabel（‘頻率/Hz‘）;?
ylabel（‘幅值‘）;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件??????16364??2018-01-03?16:32??openDISK\2012052117213095.gif

?????文件??????29109??2018-01-03?16:32??openDISK\2012052117224297.gif

?????文件???????4326??2018-01-03?11:23??openDISK\css\magic.css

?????文件???????5018??2017-12-28?16:38??openDISK\img\draw.png

?????文件???????3895??2018-01-02?21:58??openDISK\index.html

?????文件???????1669??2018-01-05?19:29??openDISK\js\F_低通.m

?????文件???????1794??2018-01-05?18:21??openDISK\js\F_帶通.m

?????文件???????1602??2018-01-05?19:06??openDISK\js\F_帶阻.m

?????文件???????1675??2018-01-05?18:11??openDISK\js\F_高通.m

?????文件???????1081??2018-01-05?19:18??openDISK\js\I_低通.m

?????文件???????1095??2018-01-05?19:27??openDISK\js\I_帶通.m

?????文件???????1062??2018-01-05?19:23??openDISK\js\I_高通.m

?????文件??????16622??2018-01-03?20:52??openDISK\js\magic.js

?????目錄??????????0??2017-12-28?15:46??openDISK\css

?????目錄??????????0??2017-12-28?16:41??openDISK\img

?????目錄??????????0??2018-01-05?19:23??openDISK\js

?????目錄??????????0??2018-01-03?16:32??openDISK

-----------?---------??----------?-----??----

????????????????85312????????????????????17