计算方法(四)特征值问题的数值解

如果不是我的理解有问题,或者因为我没做题的话,那这章就是没什么知识点。

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.补充

考前写的,填上来。
在这里插入图片描述

  • 3
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值