%自適應算法
function?[xenhancedsnr]=adaptive（filenamemnlevel）

%?if?nargin<2
%????fprintf（‘Usage:?adaptive（speech_clean1.wavnoisy_ad.wavenhanced_ad.wav）?\n\n‘）;
%????return;
%?end

[speech?Fs?nbits]=?wavread（filename）;?%?讀取清晰語音文件;

L=40;???%?濾波器長度

len=length（speech）;
n=1:len;

noise=nlevel*randn（len1）;??%加性噪聲
e=speech+noise;??%?產生帶噪語音


%?帶噪信號延遲作為參考信號輸入
for?t=2:len
????x（t:）=e（t-1）;
end
x（1）=0;

ws=zeros（Llen-L）;??%?存儲濾波器權重的矩陣，列為權重
errs=zeros（len1）;??%?誤差信號即純凈語音的估計
mus=m;

%?自適應濾波
for?i=L+1:len-L
j=i-L:1:i-1;
x1=x（j）;
y（i）=ws（:i）‘*x1;
errs（i）=e（i）-y（i）;
ws（:i+1）=ws（:i）+mus*x1*errs（i）;
end

enhanced=y*（max（e）-min（e））/（max（y）-min（y））;
SNR1?=?10*log10（var（speech‘）/var（e‘-speech‘））;???????%加噪語音信噪比
SNR2?=?10*log10（var（speech‘）/var（enhanced‘-speech（1:length（enhanced）））‘）;????%增強語音信噪比

%disp（‘自適應濾波信噪比改善量（dB）：‘）
snr=SNR2-SNR1;
x=e;
%?wavwrite（eFsnbitsoutfile1）;
%?wavwrite（errsFsnbitsoutfile2）;
%?figure（‘Name‘‘自適應濾波‘‘Numbertitle‘‘off‘）
%?subplot（311）;plot（s）title（‘純凈語音波形‘）;
%?subplot（312）;plot（e）title（‘帶噪語音波形‘）;
%?subplot（313）;plot（errs）title（‘增強語音波形‘）;

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

?????文件???????1270??2015-07-02?15:36??zengqiang1\adaptive.m

?????文件??????45022??2015-07-02?09:39??zengqiang1\enhanced_a.wav

?????文件??????48044??2015-07-02?01:56??zengqiang1\enhanced_p.wav

?????文件??????48044??2015-07-02?09:11??zengqiang1\enhanced_w.wav

?????文件??????48044??2015-06-30?18:21??zengqiang1\luyin.wav

?????文件???????8270??2015-07-01?16:44??zengqiang1\modaldlg.fig

?????文件???????6999??2015-07-01?16:44??zengqiang1\modaldlg.m

?????文件??????45102??2015-07-02?09:39??zengqiang1\noisy_a.wav

?????文件??????48044??2015-07-02?01:56??zengqiang1\noisy_p.wav

?????文件??????48044??2015-07-02?09:11??zengqiang1\noisy_w.wav

?????文件???????2889??2015-07-02?11:08??zengqiang1\pujian.m

?????文件????????166??2015-07-01?16:44??zengqiang1\readme.txt

?????文件??????16044??2015-06-30?18:21??zengqiang1\sunday_2.wav

?????文件??????25808??2015-07-02?15:55??zengqiang1\untitled.asv

?????文件??????16185??2015-07-02?15:27??zengqiang1\untitled.fig

?????文件??????25808??2015-07-02?01:10??zengqiang1\untitled.m

?????文件???????3971??2015-07-02?11:08??zengqiang1\wiener_as.m

?????文件????????940??2015-07-01?16:44??zengqiang1\ypt.m

?????目錄??????????0??2015-07-06?01:58??zengqiang1

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

???????????????438694????????????????????19