最陡下降法矢量迭代公式:
w(n+1)=w(n)+μ(−∇(n))
其中梯度
∇=Rw(n)−P代入上式:
w(n+1)=w(n)−2μ(Rw(n)−P)
而梯度为零的地方为为维纳解
w∗=P/R
所以:
w(n+1)=w(n)−2μR(w(n)−w∗)=[I−2μR]w(n)+2μRw∗
由于上式
w(n)的系数矩阵不是对角阵(不是对角阵的话求解矩阵方程计算量会很大),所以这里我把
w坐标系通过平移和旋转变成主轴坐标
v′,所以:
v′(n+1)=[I−2μΛ]v′(n+)
重点来了,通过标量迭代推导:
v′(n+1)=[I−2μΛ]nv′(0)
通过观察上式,系数矩阵是一个
n次幂阵,为了使算法稳定,而不会因为迭代次数爆炸,所以:
n→∞lim[I−2μΛ]n=0
或者写成标量表达:
n→∞lim[1−2μλk]n=0
所以:
∣1−2μλk∣<1,k=0,1,…
终于推导出来了
0<μ<λmax−1
上式中,
λmax是
R最大特征值。