clear
close?all
clc
a1=-1.6;???
a2=0.8;?
[xfs1]=audioread（‘handel.wav‘）;
[dnfs2]=audioread（‘handel_echo.wav‘）;
figure;
subplot（211）;
plot（x）;
title（‘spk參考信號(hào)‘）;
subplot（212）;
plot（dn）;
title（‘MICin-期望信號(hào)‘）;

n=max（size（x））;??%?信號(hào)中的時(shí)間點(diǎn)個(gè)數(shù)???
P=5;???%?LMS算法重復(fù)運(yùn)算5次，用于評(píng)估五次運(yùn)算產(chǎn)生的誤差的平均水平
e=zeros（1n）;???%?用于存放誤差
ep=zeros（1n）;??%?用于存放五次運(yùn)算累積的誤差
ee=zeros（1n）;??%?用于存放平方差
%算法
h=waitbar（0‘計(jì)算進(jìn)度‘）;
steps?=?P;
for?p=1:P?
L=100;?????%?濾波器階數(shù)，考慮到兩個(gè)信號(hào)之間沒有延時(shí)，且這里只是用來分離出輸入的x信號(hào)，而且誤差不做要求，因此不需要設(shè)置太高的階數(shù)?
u=0.0022;???%?增益常數(shù)?
wL=zeros（Ln）;??%?產(chǎn)生一個(gè)權(quán)向量矩陣
for?i=（L+1）:n???%?計(jì)算權(quán)向量矩陣中第三組權(quán)向量到最后一組權(quán)向量相關(guān)的變量，根據(jù)x和e=x-y求下一個(gè)y沒有d（n）
????X=x（i-L:1:（i-1））;???%?更新濾波器的參考矢量X（n）
????y（i）=X‘*wL（:i）;????%?根據(jù)x計(jì)算i時(shí)刻輸出信號(hào)?
????e（i）=dn（i）-y（i）;?????%?計(jì)算i時(shí)刻誤差信號(hào)?
????wL（:（i+1））=wL（:i）+2*u*e（i）*X;?????%?更新i時(shí)刻濾波器的權(quán)向量?
????ee（i）=e（i）^2;???%?更新平方差
end?
ep=ep+ee;??%?平方差累積
waitbar（p/steps）;
end
close（h）;
eq=ep/P;????%?五十次重復(fù)計(jì)算平方差求均值
a1L=-wL（21:n）;???%?a1在LMS算法下值的變化，wL矩陣中第一行的1到n個(gè)數(shù)，權(quán)向量矩陣第一行?
a2L=-wL（11:n）;???%?a2在LMS算法下值的變化?，wL矩陣中第二行的1到n個(gè)數(shù)，權(quán)向量矩陣第二行?

%畫圖
figure;
subplot（321）;
plot（x）;?
title（‘需要采集的聲音Voice‘）;????%?根據(jù)w產(chǎn)生的隨機(jī)信號(hào)x

subplot（322）;
plot（dn）;
title（‘被回聲污染的了Voice，就是MICin‘）;???%?五十次運(yùn)算誤差求均值

subplot（323）;
plot（y）;?
title（‘合成回聲信號(hào)‘）;????%?根據(jù)w產(chǎn)生的隨機(jī)信號(hào)x

subplot（324）;
plot（e）;
title（‘誤差‘）;???%?五十次運(yùn)算誤差求均值


subplot（325）;?
plot（a1L‘r-‘）;?%?權(quán)向量矩陣第一行
hold?on;
plot（a2L‘k-‘）;?%?權(quán)向量矩陣第二行
title（‘權(quán)向量參數(shù)的學(xué)習(xí)過程‘）;??
legend（‘a(chǎn)1L‘‘a(chǎn)2L‘）;

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2018-12-21?13:27??demo4\
?????文件????????1862??2018-12-21?12:53??demo4\AdaFilter.m
?????文件?????????427??2018-12-21?13:28??demo4\AdaFilter2.m
?????文件?????????325??2018-05-13?22:22??demo4\add_noise.m
?????文件??????146270??2018-12-21?12:01??demo4\handel.wav
?????文件??????146270??2018-12-21?12:01??demo4\handel_echo.wav
?????文件????????2491??2018-12-21?11:53??demo4\kalman.m
?????文件?????????729??2018-05-13?20:52??demo4\Kframe.m
?????文件???????45100??2018-12-21?13:26??demo4\signalnoisy.wav
?????文件??????107294??2000-03-28?23:16??demo4\six.wav
?????文件?????????253??2018-12-21?11:22??demo4\SNR_singlech.m
?????文件???????45102??2005-07-12?08:58??demo4\sp01.wav
?????文件????????1978??2018-05-13?22:22??demo4\spectruesub.m
?????文件????????2849??2018-12-21?11:28??demo4\SpeechEnhancement.asv
?????文件????????2706??2018-12-21?11:53??demo4\SpeechEnhancement.m
?????文件????????1002??2018-05-13?22:23??demo4\stationary_noise_evaluate.m
?????文件????????1874??2018-05-13?22:23??demo4\wienerfilter.m