如果不是我的理解有问题,或者因为我没做题的话,那这章就是没什么知识点。
1.乘幂法
两步一直重复:
y
(
k
)
^{(k)}
(k)=x
(
k
)
^{(k)}
(k)/(max(x
(
k
)
^{(k)}
(k)))
x
(
k
+
1
)
^{(k+1)}
(k+1)=A*y
(
k
)
^{(k)}
(k)
max(x
(
k
)
^{(k)}
(k))就是取x
(
k
)
^{(k)}
(k)中最大的分量。
看一下P
82
_{82}
82例4.1就好
有个重点:y才是特征向量
2. 反幂法
在上面基础上改一点:
y
(
k
)
^{(k)}
(k)=x
(
k
)
^{(k)}
(k)/(max(x
(
k
)
^{(k)}
(k)))
A*x
(
k
+
1
)
^{(k+1)}
(k+1)=y
(
k
)
^{(k)}
(k)
3.雅各比旋转
直接看P 91 _{91} 91例4.5吧
找到非主对角线元素中绝对值最大的元素a i j _{ij} ij,通过tan2θ=2a i j _{ij} ij/(a j j _{jj} jj-a i i _{ii} ii),计算出θ的cos值和sin值。
先列一个主对角线都是1,其余全是0的单位矩阵,然后将矩阵中的这四个点换掉就得到P的结果:
对于第i行、第i列、第j行、第j列形成的正方形的四个角(即a
i
i
_{ii}
ii、a
i
j
_{ij}
ij、a
j
i
_{ji}
ji、a
j
j
_{jj}
jj四个点),用:
这个方式替换掉,其他位置的元素不变。
已知P的计算方式,就可以循环利用下面这个公式迭代:
A k + 1 _{k+1} k+1=P k T ^{T}_{k} kT*A k _{k} k*P k _{k} k
最后得到只有主对角线有数字的矩阵A k _{k} k,主对角线上数字放在一起变成的向量就是特征根λ,S k _{k} k=P 1 _{1} 1*P 2 _{2} 2*P 3 _{3} 3……*P k _{k} k,其中S k _{k} k的每列都是一个特征向量。
注意:书中例题在P 92 _{92} 92最后求特征向量时公式写错了,其中每一个P都不带转置,即我上面所写,公式在P 90 _{90} 90。
4.补充
考前写的,填上来。