clear?all;close?all;clc
%斷點(diǎn)位置：14、22、41、61行
%讀入無噪聲音樂文件y;采樣頻率為fs
[yfs]=audioread（‘Clean.wav‘）;
x=y（:1）‘;
%?sound（xfs）;%播放音樂文件
N=length（x）;%音樂文件的采樣點(diǎn)數(shù)總和（采樣長(zhǎng)度）
fft_x=fft（xN）;%音樂信號(hào)的傅里葉變換
w=fs/N*（0:N-1）;%頻譜空間的譜線間隔矩陣，即采樣頻率間隔，fs為歌曲的采樣頻率，它最高
figure（1）subplot（211）;plot（x）title（‘Clean?Music?Signal‘）;
subplot（212）plot（wabs（fft_x））title（‘Frequency?domain?of?Clean?Music?Signal?‘）;%可以看出，音樂文件最高頻率是大概8000Hz

%讀入有噪聲音樂文件
[ynfs]=audioread（‘Corrupt.wav‘）;
xn=yn（:1）‘;
fft_xn=fft（xnN）;%音樂信號(hào)的傅里葉變換
?sound（xnfs）;
figure（2）subplot（211）;plot（xn）title（‘Noise?Music?Signal‘）;
subplot（212）plot（wabs（fft_xn））title（‘Frequency?domain?of?Noise?Music?Signal‘）;%對(duì)比figure（1）中的波形可以看出在750和1750Hz處各存在一個(gè)噪聲，因此需要設(shè)計(jì)帶阻濾波器

%%?band-pass?filter
a=0.16;
b=0.22;
wp=[a-0.05b+0.05];
ws=[ab];
Rp=3;
As=13;
[Nbwc]=buttord（wpwsRpAs）;
[BSbASb]=butter（Nbwc‘stop‘）;
fk=0:1/8000:1;
Hk=freqz（BSbASbfk*pi）;
figure（3）
plot（fk20*log（abs（Hk）））;
%對(duì)噪聲信號(hào)濾波
F_y1=filter（BSbASbxn）;
fft_F_y1=fft（F_y1N）;%濾波后音樂信號(hào)傅里葉變換
?sound（F_y1fs）;
figure（4）subplot（211）plot（F_y1）title（‘1st?Music?signal?after?filtering‘）
subplot（212）plot（wabs（fft_F_y1））title（‘Frequency?domain?of?Music?signal?after?1st?filtering‘）;

c=0.4;
d=0.46;
wp=[c-0.05d+0.05];
ws=[cd];
Rp=3;
As=13;
[Nbwc]=buttord（wpwsRpAs）;
[BSbASb]=butter（Nbwc‘stop‘）;
fk=0:1/8000:1;
Hk=freqz（BSbASbfk*pi）;
figure（5）
plot（fk20*log（abs（Hk）））;
%對(duì)噪聲信號(hào)濾波
F_y2=filter（BSbASbF_y1）;
fft_F_y2=fft（F_y2N）;%濾波后音樂信號(hào)傅里葉變換
sound（F_y2fs）;
figure（6）subplot（211）plot（F_y2）title（‘2nd?Music?signal?after?filtering‘）
subplot（212）plot（wabs（fft_F_y2））title（‘Frequency?domain?of?Music?signal?after?2nd?filtering‘）;

%求平均值方法濾除白噪聲
music=audioread（‘Clean.wav‘）;
music=music（:1）‘;
sum=F_y2;
for?ii=1:1999
????r=rand（1N）;
????x2=r-0.5;
????yn=music+x2;
????sum=sum+yn;
end
sum=sum/2000;
sound（sumfs）;
fft_sum=fft（sumN）;
figure（7）subplot（211）plot（sum）title（‘Music?signal?after?averaging‘）;
subplot（212）plot（wabs（fft_sum））title（‘Frequency?domain?of?Music?signal?after?averaging‘）;

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----

?????文件???????2304??2018-11-20?22:33??BS.m

?????文件??????50044??2018-11-15?13:05??Clean.wav

?????文件?????100044??2018-11-15?13:05??Corrupt.wav

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

???????????????152392????????????????????3