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參考信號‘）;
subplot（212）;
plot（dn）;
title（‘MICin-期望信號‘）;

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

%畫圖
figure;
subplot（321）;
plot（x）;?
title（‘需要采集的聲音Voice‘）;????%?根據w產生的隨機信號x

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

subplot（323）;
plot（y）;?
title（‘合成回聲信號‘）;????%?根據w產生的隨機信號x

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


subplot（325）;?
plot（a1L‘r-‘）;?%?權向量矩陣第一行
hold?on;
plot（a2L‘k-‘）;?%?權向量矩陣第二行
title（‘權向量參數的學習過程‘）;??
legend（‘a1L‘‘a2L‘）;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????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