clearallclcm500将输入量赋值u0002将输入量赋值a116a208w110w210w120w220
11
22trueseque
cex
rdra
d
1mx1rd1x2rd2fork3mxkrdka1xk1a2xk2e
dRx1x1000T0e10W01Xx2x1e2x2WXfori31m权系数迭代m次Rxi1xi1xi1xi2xi1xi2xi2xi2列出自相关矩阵T1TtraceR求出迹的值，为后续u的判断做准备ifuTerroruislargertha
1tR判断u的值是否小于迹的倒数e
dWW2uXei1LMS算法的权系数迭代公式Xxi1xi2LMS算法中输入信号矢量的递推eixiWXw1iW1LMS算法中权系数a1的提取
iie
dlambdaeigR由R得出特征值v11给主轴坐标赋初值for
1150w
1v1v1lambda112ulambda12
1v2v2lambda212ulambda22
1学习曲线的迭代公式e
dfigure1
fplotw1r曲线绘图titleLMS算法权值收敛情况xlabel迭代次数ylabel权值变化grido
holdo
m500权系数迭代m次le
098遗忘因子a116a208w111w211w121w222
11
22输入序列产生rdra
d
1m生成白噪声序列x1rd1x2rd2产生前两个输入序列fork3mxkrdka1xk1a2xk2产生其余输入序列e
dR0000自相关矩阵初值R10Ww12w22初始权值P00fori31mXxi1xi2Rle
RXX迭代公式中自相关矩阵的计算exiWX输出信号误差e
1WWi
vRXeRLS算法的权系数迭代公式w1iW1RLS算法中权系数a1的提取w2iW2RLS算法中权系数a2的提取
iie
dplotw1b曲线绘图title自适应权系数a1
的过渡过程（RLS和LMS算法比较）xlabel迭代次数ylabel权值变化lege
dLMSRLSgrido

fholdofffigure2plotxb
曲线绘图
fu
ctio
varargoutLSLvarargi
v1m80ble
gthvifb1eb010eb110eb210初始化ef010ef110ef210v010v110v210r011r111r211sf01vsf11vsf21vsb01vsb11vsb21vrdra
d
1m信号产生x1rd1x2rd2fori3mxi1558xi1081xi2rdie
dfori2m按时间迭代eb0ixief0ixisf0isf0i1xixisb0isf0ir0i1e
dfori2mv1iv1i1eb0i1ef0ir0i1ef1ief0iv1ieb0i1sb0i1eb1ieb0i1v1ief0isb0isf1isf0iv1i2sb0i1sb1isb0i1v1i2sf0ir1i1r0i1eb0i1eb0i1sb0i1kf1iv1r