clear?all?????????????????%釋放系統(tǒng)存儲(chǔ)空間
clf???????????????????????%清除運(yùn)行中的圖
[sfs]=audioread（‘guxiang.wav?‘）;%讀取音頻文件
d=s;??????????????????????%期望得到信號(hào)
x=awgn（s1‘measured‘）;?%音頻加入高斯噪聲
N=length（x）;??????????????%取讀采樣的長(zhǎng)度
ns=x-s;???????????????????%高斯白噪聲信號(hào)
u=x;??????????????????????%需要濾波的信號(hào)
%%%%%%%%%%利用LMS算法進(jìn)行自適應(yīng)濾波%%%%%%%%%%%%
mu=0.05;??????????????????%收斂步長(zhǎng)
M=100;?????????????????????%32階濾波器
w=zeros（M1）;?????????????%初始化
u=u（:）;
d=d（:）;
e=zeros（1N）;
for?n=M:N?????????????????%開(kāi)始迭代
????uvec=u（n:-1:n-M+1）;???%
????e（n）=d（n）-w‘*uvec;????%y（n）=w‘*