%---------------LMS法------------------%
clc;
clear?all;
close?all;
%%?參數設置
t=0:1:255;
DataLen=length（t）;%數據長度
sigma=0.2;%噪聲方差
order=50;%濾波器階數??
f=0.05;%正弦信號頻率
%%?開始
SinData=sin（2*pi*f*t）;%產生正弦信號
Data=SinData+sigma*randn（1DataLen）;%加入噪聲
d=SinData;
%----------自相關矩陣??
R=zeros（orderorder）;
for?ii=1:order
????for?jj=ii:1:DataLen
????????R（1ii）=R（1ii）+Data（1jj）*Data（1jj-ii+1）;
????end?
????R（1ii）=R（1ii）/（DataLen-ii+1）;
end??
for?ii=2:order??
????for?jj=2:order??
????????R（iijj）=R（ii-1jj-1）;??
????end??
end??
R_T=R‘;??
for?ii=1:order??
????R_T（iiii）=0;??
end??
R=R+R_T;??
inv_R=inv（R）;
%----------求最大特征值
[u?s?v]=svd（R）;
lamda_max=s（11）;
miu=1/lamda_max-0.05;
%----------輸出數據及誤差
W=zeros（order1）;
y=zeros（1DataLen）;
for?ii=order:DataLen
????xn=Data（ii:-1:ii-order+1）;
????y（ii