关于递推最小二乘算法的程序，请教各位高手了，先说明如下：

给出m个n维向量q(k)和m个实数y(k),给出m*m阶矩阵p(0)(为初始的),向量w(0)(为初始的),常数a,最大计算次数K'

先计算g(k)=(p(k-1)q(k))/(a+q(k)p(k-1)q(k))(式1)以及p(k)=(1/a)*(p(k-1)-g(k)q(k)p(k-1))(式2),

再按照下式计算w(k')=w(k'-1)+g(k')[y(k)-q(k')w(k'-1)],(式3)

计算J(k')=J(k'-1)+(y(k)-q(k')w(k'-1)),判断|J(K')-J(K'-1)|是否小于给定某一值b

若是则取w(k*)=w(k'),不是则判断k'是否小于K',若是k'=k'+1,再计算 w(k'),即从式3计算,不是,则再计算g(k),p(k)即从(式1)(式2),

说明一下:其中涉及到矩阵运算的,没有标出转置,见谅