矩阵特征值问题的数值方法

合集下载

矩阵特征值问题的数值方法.

矩阵特征值问题的数值方法.

矩阵特征值问题的数值方法矩阵特征值设A 是n 阶矩阵,x 是非零列向量. 如果有数λ 存在,满足那么,称x 是矩阵A 关于特征值λ的特征向量. 很显然一般地有主特征值的乘幂迭代法设n 阶矩阵A 的n 个特征值按模从大到小排序为:n 其对应的n 个线性无关的特征向量分别为:设是任意一个非零的n 维向量,则:假设,构造一个向量序列:则:或者:当时:如果是矩阵A 的关于特征值的一个特征向量,特征值个特征那么对于任意一个给定的,也是特征值的特征向量。

所以,是对主特征值对应的特征向量的近似。

如果则会变得很大或者如果,则会变得很大,或者如果,则会变得非常小,在实际计算中,为避免这种情况的出现需对做归一化处理况的出现,需对做归一化处理:由:左乘得:所以主特征值的近似值所以主特征值的近似值:残余误差向量定义为:当迭代次数充分大时,残余误差将充分小。

逆乘幂法:类似地,也可以求模最小特征值和对应的特征向量特征向量。

上述问题的主特征值问题就是矩阵A 的模最小特征值问题。

结果,逆乘幂法的迭代公式为:在实际应用中,无需计算逆矩阵,但需求解线性系统实对称矩阵的基本定理:对实对称矩阵A ,一定存在一个正交相似变换使得为对角矩阵且其对角矩阵P ,使得:为对角矩阵,且其对角的特征值元素为矩阵A 的特征值。

相似变换:相似变换保持矩阵特征值(但不是特征向量)不变不变。

(证明略)正交相似变换:中。

正交相似变换的例子—坐标旋转:叫旋转矩阵。

容易验证:。

适当选择旋转角,可消去xy 项—得到对角阵D 。

矩阵特征值问题的数值方法实对称矩阵的基本定理再看下面的例子:令:O 平面的坐标旋转变换适当同样地有:。

则是在x-O-z 平面的坐标旋转变换。

适当x z —D 。

选择旋转角可消去z 项得到对角阵实对称矩阵的Jacobi 方法:全部特征值和特征向量根据实对称矩阵的基本定理,求得矩阵A 的全部特征值的关键是找到正交相似变换矩阵P 使部特征值的关键,是找到正交相似变换矩阵P ,使得为对角阵。

特征值和特征向量计算的数值方法

特征值和特征向量计算的数值方法

特征值和特征向量计算的数值方法在数学和计算机科学领域中,特征值和特征向量是非常重要的概念。

特征值和特征向量的计算有许多不同的数值方法,本文将介绍其中一些常见的数值方法,并分析它们的优劣和适用范围。

一、特征值和特征向量的定义在矩阵理论中,给定一个n×n的矩阵A,如果存在一个非零向量v和一个标量λ,使得Av=λv,那么称v为矩阵A的特征向量,λ为矩阵A的特征值。

特征值和特征向量的计算可以帮助我们理解矩阵的性质以及解决一些实际问题。

二、幂法幂法是计算特征值和特征向量的常用数值方法之一。

幂法的基本思想是通过多次迭代,逐渐逼近矩阵的特征值和特征向量。

具体操作如下:1. 初始化一个非零向量b0;2. 进行迭代计算:bi+1 = A * bi / ||A * bi||;3. 取出近似特征向量的最后一列:v = bn;4. 进行迭代计算特征值:λ = (Av)T * v / (vT * v)。

幂法的主要优点是简单易懂,且只需要进行矩阵向量乘法和内积计算。

然而,幂法仅能求取具有最大特征值的特征向量,而且对于存在多个特征值相等的情况并不适用。

三、反幂法反幂法是幂法的一种改进方法,用于求取矩阵A的最小特征值和对应的特征向量。

反幂法的基本步骤如下:1. 初始化一个非零向量b0;2. 进行迭代计算:bi+1 = (A - μI)^-1 * bi / ||(A - μI)^-1 * bi||;3. 取出近似特征向量的最后一列:v = bn;4. 进行迭代计算特征值:λ = (Av)T * v / (vT * v)。

反幂法的改进之处在于引入了矩阵的逆运算,通过使用矩阵A减去一个合适的常数μ乘以单位矩阵来实现。

反幂法适用于矩阵A的特征值接近于μ的情况。

四、QR方法QR方法也是一种常用的特征值计算方法,它适用于求解所有特征值以及对应的特征向量。

QR方法的基本思想是将一个矩阵分解为正交矩阵Q和上三角矩阵R的乘积,然后迭代地将矩阵A转化为更接近上三角形的形式。

数值分析-第7章 矩阵特征值问题的数值解法n

数值分析-第7章 矩阵特征值问题的数值解法n
(-3.406542,-2.460280, 6.920561) (-2.832406, -2.028615, 6.210333)
7
9 11 12
6.104716
6.026349 6.006637 6.003327
(-0.450275, -0.322058, 1.0)
(-0.445914, -0.318617, 1.0) (-0.444814, -0.31775, 1.0) (-0.444630, -0.317606, 1.0)
其中i为A的特征值,P的各列为相应于i的特征向量。
P -1 AP D
2

n
2
定理7.1.3 ARnn,1, …, n为A的特征值,则
(1)A的迹数等于特征值之和,即 tr ( A) aii i
i 1 i 1
n
n
(2)A的行列式值等于全体特征值之积,即
1 xi(k +1) / xi(k )
i 1,2,, n
可见,当k充分大时, ( k ) 近似于主特征值, ( k +1) 与x ( k )的对应非零分量的比值 x x 近似于主特征值。
在实际计算中需要对计算结果进行规 , 范化。因为当 1 1时,x (k ) 趋于零, 当1 1时, x ( k )的非零分量趋于无穷。 从而计算时会出现下溢 或上溢。
特征值的范围. 解 我们先分别求出各个圆盘区域。 D1 = {z:|z – 1|£0.6};D2 = {z:|z – 3|£0.8} D3 = {z:|z + 1|£1.8};D4 = {z:|z + 4|£0.6}. 易见D2和D4为 弧立圆盘分别 包含A的两个实 特征值.

求矩阵特征值的方法

求矩阵特征值的方法

求矩阵特征值的方法矩阵特征值是矩阵理论中的一个重要概念,它在许多领域中都有着广泛的应用,如物理学、工程学、计算机科学等。

求矩阵特征值的方法有多种,下面将介绍其中的三种常用方法。

一、特征多项式法特征多项式法是求矩阵特征值的一种常用方法。

它的基本思想是将矩阵A与一个未知数λ相乘,得到一个新的矩阵B=A-λI,其中I为单位矩阵。

然后求解矩阵B的行列式,得到一个关于λ的多项式,称为特征多项式。

矩阵A的特征值就是使特征多项式等于零的λ值。

具体步骤如下:1. 构造矩阵B=A-λI。

2. 求解矩阵B的行列式det(B)。

3. 解特征多项式det(B)=0,得到矩阵A的特征值λ。

二、幂法幂法是求矩阵特征值的一种迭代方法。

它的基本思想是从一个任意的非零向量开始,不断地将其乘以矩阵A,直到向量的方向趋于特征向量的方向,同时向量的模长趋于特征值的绝对值。

具体步骤如下:1. 选择一个任意的非零向量x0。

2. 迭代计算xn+1=Axn/||Axn||,其中||Axn||为Axn的模长。

3. 当xn+1与xn的差值小于某个预设的精度时,停止迭代,此时xn 的模长即为矩阵A的最大特征值,xn/||xn||即为对应的特征向量。

三、QR分解法QR分解法是求矩阵特征值的一种数值方法。

它的基本思想是将矩阵A 分解为QR,其中Q为正交矩阵,R为上三角矩阵。

然后对R进行迭代,得到一个对角矩阵,对角线上的元素即为矩阵A的特征值。

具体步骤如下:1. 对矩阵A进行QR分解,得到A=QR。

2. 对R进行迭代,得到一个对角矩阵D,对角线上的元素即为矩阵A的特征值。

以上三种方法都有其优缺点,具体选择哪种方法取决于实际应用场景和计算需求。

在实际应用中,还可以结合多种方法进行求解,以提高计算精度和效率。

矩阵特征值的数值解法

矩阵特征值的数值解法

矩阵特征值的数值解法矩阵的特征值是在矩阵与其特征向量之间的关系中的数值解。

特征值在各个领域中都有广泛应用,包括物理、工程、金融等。

在解决实际问题时,我们经常需要计算矩阵的特征值,因此研究如何求解矩阵特征值的数值方法是非常重要的。

1. 幂迭代法(Power Iteration)幂迭代法是求解矩阵特征值的一种简单而常用的数值方法。

它的基本思想是通过不断迭代矩阵与向量的乘积,使得向量趋近于该矩阵的一个特征向量。

具体步骤如下:(1)初始化一个非零的初始向量x。

(2)进行迭代计算,即$x^{(k+1)}=Ax^{(k)}/,Ax^{(k)},$。

(3)当向量x的相对误差小于一些预设的精度要求时,停止迭代,此时的x即为矩阵A的一个特征向量。

(4)将x带入特征值的定义式$\frac{Ax}{x}$,计算出特征值。

幂迭代法的优点是简单易实现,计算速度较快,缺点是只能求解特征值模最大的特征向量,而且对于存在特征值模相近的情况,容易收敛到错误的特征值上。

2. QR迭代法(QR Iteration)QR迭代法是一种较为稳定的求解矩阵特征值的数值方法。

它的基本思想是通过不断进行QR分解,使得矩阵的特征值逐渐收敛。

具体步骤如下:(1)将矩阵A进行QR分解,得到正交矩阵Q和上三角矩阵R,令$A_1=RQ$。

(2)将$A_1$再次进行QR分解,得到新的矩阵$A_2=R_1Q_1$。

(3)重复步骤(2),直到得到收敛的矩阵$A_k$,此时$A_k$的对角线上的元素即为矩阵A的特征值。

QR迭代法的优点是对于特征值模相近的情况仍然能够收敛到正确的特征值上。

缺点是每次QR分解都需要消耗大量的计算量,迭代次数较多时计算速度较慢。

3. Jacobi迭代法(Jacobi's Method)Jacobi迭代法是一种通过对称矩阵的对角线元素进行迭代操作,逐步将非对角元素变为零的求解特征值的方法。

具体步骤如下:(1)初始化一个对称矩阵A。

矩阵特征值和特征向量的数值解法

矩阵特征值和特征向量的数值解法
(k )
的常用方法是迭代每一步对向量 u
规范化。引入函数 max( u
(k )
) ,它表示取
向 量 u (k ) 中 按模 最大 的分 量,例 如, u (k ) =(2,-5,4)T,则 max( u (k ) )=-5,这 样
u (k ) 的最大分量为 1,即完成了规范化。 (k ) max (u )
7.1 幂法
7.1.1 幂法原理及实用幂法 幂法主要用于求矩阵按模最大的特征值和相应的特征向量。设矩阵,2,..., n) 满足:
| λ1 |>| λ 2 |≥| λ3 |≥ ... ≥| λ n | (7.1.1)
相应的 n 个特征向量 xi (i = 1,2,..., n) 线性无关。上述假设表明, λ1 为非零单 实根, x1 为实特征向量。
k →∞
k →∞
lim v ( k ) =
x1 max( x1 )
事实上,由式(7.1.5)知
v
(k )
=
Ak u ( 0 )
∏m
i =0
k
i
算法 7.1.1 实用幂法 (1) 输入: aij (i, j = 1,2, L n), ui (i = 1,2, L), ε ; (2) k = 1; m0 = max(ui );
7.1 幂法
幂法基本原理是:任取非零实向量 u
(0)
,做迭代
u ( k ) = Au ( k −1) = Ak u ( 0 ) (k = 1,2,...)

( 7 .1 . 2 )
λ1 = lim
这里 u j 表示向量 u
(k ) (k )
u (jk +1) u (jk )
k →∞

矩阵的特征值与特征向量的数值解法

矩阵的特征值与特征向量的数值解法

第八章 矩阵地特征值与特征向量地数值解法某些工程计算涉及到矩阵地特征值与特征向量地求解 .如果从原始矩阵出发,先求 出特征多项式,再求特征多项式地根,在理论上是无可非议地•但一般不用这种方 法,因为了这种算法往往不稳定•常用地方法是迭代法或变换法•本章介绍求解特 征值与特征向量地一些方法•§ 1乘幂法乘幕法是通过求矩阵地特征向量来求特征值地一种 迭代法,它适用于求矩阵 地按模最大地特征值 及对应地特征向量.b5E2RGbCAP 定理8 • 1设矩阵Ai x n 有n 个线性无关地特征向量 X<i=1,2,…,n ),其对应地特征 值入 i (i =1,2,…,n> 满足 plEanqFDPw|入1|>|入2|三…三|入n |则对任何n 维非零初始向量 乙,构造Z k = AZ k-1(k=1,2.其中(Z k >j 表示向量Z<地第j 个分量. 证明:只就入i 是实数地情况证明如下 因为A 有n 个线性无关地特征向量X,<i = 1,2,用X<i = 1,2, …,n )线性表示,即Z 0=a 1X 1 + 用A 构造向量序列{Z k }其中由矩阵特征值定义知 AXm i X(i=1,2,…,n>,故Z k 二A k Z^ :1A k X^ : 2A k X 2nA kX n 「T ;X1 *〉2';X2- :'n'n Xn同理有li m (ZQ j_______________ <22?=■ 1<8 • 1) Z 1 二 AZ 0,乙二 AZ= A^Z。

,川,Zk-AZ kj-A Zo(8・2>- k' nkTX ii zz2-nJ 2-7k -AZk」=人X ii =2<A1」<8.3)<8.4 ),设a 1工0,并且注意到…,n )所以任何非零向量Z o 都可 a 2茨 + …+a nX <a 1 工 0) DXDiTa9E3d将<8.3 )与<8.4 )所得乙及Z k-1地第j个分量相除| 入i|<| 入…,n> 得RTCrpUDGiT1|(i=1,2,定理8 • 1地证明过程实际上是给出了矩阵地按模最大特征值地计算方法:1) 先任取一非零向量Z 0, 一般可取Z o =(1,1,1> T; 2) 按<8.2 )式计算 乙=AZ -i (k=1,2,…>;3)当K足够大时,即可求出詔;=6为了减少"1对于所选地第j个分量地依赖性,还可用各个分量比地平均值来代替,即关于对应于入1地特征向量地计算:由<8.1 )知,当k 充分大时,Z k =入1Z k-1,又由迭代式 Z k = AZ k-1,可知AZ k-1 =入1Z k-1故 由特征值定义知 Z k-1即为入1对应地特征向量,或Z k =入1Z k-1为入1对应地特征向 量.5PCzVD7HxA这种求矩阵地按模最大特征值及其对应特征向量地方法称为 乘幕法. 应用乘幕法计算A 地按模最大特征值入1和对应特征向量时,由<8.3)易知Z k = *-n厲入+送码J y1X ii 2当|入1|>1或|入1|<1时,Z k 中不为零地分量将会随 K 地增大而无限增大,或随K 地 「 ------------ 增大而趋于零,用计算机计算就会出现“上溢”或“下溢” .为了克服这个缺点,一」无 穷 常将迭代向量 乙先规范化,然后再计算,具体做法是:jLBHrnAILg 一,一用max (Z>S 示向量Z k 地绝对值最大地分量,任取一初始向量Z o =a 1X 1+ a 汎+…+ a n X^V a 1工0)构造与<8.2 )对应地向量序列.xHAQX74J0XAZ o由<8.3)可知Yk = maZk A kZ o max A kZ o max n:X 亠1 1 j ii =2X inM • r ii -2X i丿丿(k tmax X i<8.7J 二 AYA 2Z omax AZ0J 'max 乙max AZ oA 2Z 。

矩阵特征值问题的数值计算

矩阵特征值问题的数值计算

矩阵特征值问题的计算方法特征值问题:A V=λV¾直接计算:A的阶数较小,且特征值分离得较好 特征值:det(λI-A)=0,特征向量:(λI-A)V=0¾迭代法:幂法与反幂法¾变换法:雅可比方法与QR方法内容:一、 特征值的估计及其误差问题二、 幂法与反幂法三、 雅可比方法四、 QR方法一、 特征值的估计及其误差问题 (一)特征值的估计结论 1.1:n 阶矩阵()ij n n A a ×=的任何一个特征值必属于复平面上的n 个圆盘:1,||||,1,2,ni ii ij j j i D z z a a i n =≠⎧⎫⎪⎪=−≤=⎨⎬⎪⎪⎩⎭∑"(10.1) 的并集。

结论1.2:若(10.1)中的m个圆盘形成一个连通区域D,且D与其余的n-m个圆盘不相连,则D中恰有A的m个特征值。

(二)特征值的误差问题结论1.3:对于n 阶矩阵()ij n n A a ×=,若存在n 阶非奇异矩阵H ,使得11(,,)n H AH diag λλ−=Λ=", (10.2)则11min ||||||||||||||i p p p i nH H A λλ−≤≤−≤∆ (10.3)其中λ是A A +∆的一个特征值,而(1,,)i i n λ="是A 的特征值,1,2,p =∞。

结论1.4:若n 阶矩阵A 是实对称的,则1min ||||||i p i nA λλ≤≤−≤∆。

(10.4)注:(10.4)表明,当A 是实对称时,由矩阵的微小误差所引起的特征值摄动也是微小的。

但是对于非对称矩阵而言,特别是对条件数很大的矩阵,情况未必如此。

二、 幂法与反幂法(一) 幂法:求实矩阵按模最大的特征值与特征向量假设n 阶实矩阵A 具有n 个线性无关的特征向量,1,iV i n =",则对于任意的0nX R ∈,有 01ni ii X a V ==∑,从而有01111112((/))n nk k k i i i i ii i nk k i i i i A X a A V a V a V a V λλλλ======+∑∑∑.若A 的特征值分布如下:123||||||||n λλλλ>≥≥≥",则有01111()k kk A X a V λλ→∞⎯⎯⎯→为对应的特征向量须注意的是,若1||1λ<,则10kλ→,出现“下溢”,若1||1λ>,则1kλ→∞,出现“上溢”,为避免这些现象的发生,须对0kA X 进行规范化。

方法求矩阵全部特征值

方法求矩阵全部特征值

方法求矩阵全部特征值求解矩阵的全部特征值是一个重要的问题,它在线性代数和数值计算中都有广泛的应用。

在本文中,将详细介绍几种方法来求解矩阵的全部特征值,包括特征值分解方法、幂迭代方法、QR方法、Jacobi方法和带位移的QR方法。

特征值是一个矩阵的最重要的性质之一,它描述了矩阵的行为和性质。

特征值可以用于计算矩阵的条件数、正交变换、矩阵的相似性和对角化等。

求解矩阵的全部特征值可以通过特征值分解来实现。

特征值分解是将一个矩阵分解成一个对角矩阵和一个特征向量矩阵的乘积。

根据特征值分解的定义,可以得到以下公式:A=QΛQ^(-1)其中A是一个n×n的矩阵,Λ是一个对角矩阵,Q是一个n×n的正交矩阵,^(-1)表示矩阵的逆。

通过特征值分解,可以求解矩阵的全部特征值和对应的特征向量。

特征值分解的方法有很多种,比如QR方法、Jacobi方法和带位移的QR方法。

幂迭代是一种求解矩阵最大特征值和对应特征向量的迭代方法。

幂迭代的基本思想是通过不断迭代矩阵的幂次来逼近最大特征值和对应特征向量。

幂迭代的过程可以通过以下公式表示:x(k+1)=Ax(k)其中x(k)表示第k次迭代的特征向量,A表示待求解的矩阵。

幂迭代的收敛性取决于一个非零初始向量的选择和特征值的大小。

当初始向量与最大特征值对应的特征向量接近时,幂迭代可以得到最大特征值的逼近值。

通过迭代可以不断逼近最大特征值,同时得到对应特征向量。

QR方法是一种求解实对称矩阵全部特征值的方法,它通过迭代将矩阵变换为上三角矩阵。

QR方法的基本步骤包括QR分解、矩阵相似变换和迭代。

在每一次迭代中,矩阵A都被变换为一个上三角矩阵R,并且特征值逐步靠近对角线的元素。

Jacobi方法是一种通过旋转矩阵将矩阵对角化的方法。

Jacobi方法的基本思想是通过多次相似变换将矩阵的非对角元素逐步置为零,使得矩阵对角化。

Jacobi方法的关键步骤是选择旋转角度和旋转矩阵,通过旋转操作将非对角元素置为零。

线性方程组与矩阵特征值求解的数值方法

线性方程组与矩阵特征值求解的数值方法

(m-k)次乘法运算
A(k与1) A前(k)k行元素相同, A(k的1)左上角k阶阵
a (1) 11
A(k) 11
a1(kk )
为上三角阵。
a
(k kk
)
第k步约化公式:
Lk A(k ) A(k1)
Lk
b
(k
)
b(k 1)
(3)继续上述约化过程,且设a
(k kk
)
0(k
1,2,
, s),
(2.1)
A~
1 2
4 5
7 8
1
1
( (
E2 E3
) )
2( 3(
E1 E1
) )
E E
2 3
1 0
4 3
7 6
1(E3 ) 2(E2 ) E3 1
1
0
4 3
7 6
11,
3
6
11
1
0
6
10
2
0
0
2
0
(2)回代求解,得:x1
1, 3
x
2
1, 3
x
3
0。
结论:
整个计算过程可分为两部分:(1)消元:把原 方程组转化为系数矩阵为上三角矩阵的方程组; (2)回代:由系数矩阵为上三角矩阵的方程组求解
0, 计算乘数
m ik
a(k) ik
a
(k kk
)
,(i
(m-k)次除法运算 k 1, , m),
k)
)进行行初等变换,使
A( k )第k列
a(k kk
) 以下元素约为零,
即 ri mikrk ri (i k 1,L , m) ,得到与原方程组等价的方程组 A(k 1) x b(k 1)

线性方程组与矩阵特征值求解的数值方法

线性方程组与矩阵特征值求解的数值方法

线性方程组与矩阵特征值求解的数值方法线性方程组与矩阵特征值求解是线性代数中的两个重要问题。

线性方程组解决了形如Ax=b的方程组,其中A为一个m×n的矩阵,b为一个m 维的向量,求解x使得该方程组成立。

矩阵特征值求解是求解形如Ax=λx的特征值和特征向量问题,其中A为一个n×n的矩阵,λ为特征值,x为特征向量。

这两个问题在实际应用中有广泛的应用,如计算机图形学、仿真和优化等领域。

本文将介绍线性方程组和矩阵特征值求解的数值方法。

一、线性方程组的求解方法1.1直接法直接法是指通过一系列的代数运算和变换直接求解线性方程组的解。

经典的直接法有高斯消元法、LU分解法和Cholesky分解法等。

这些方法的时间复杂度通常为O(n^3)。

直接法的优点是解的精度高,稳定性好,适用于小规模的问题。

1.2迭代法迭代法是指通过迭代计算逼近线性方程组的解。

迭代法的基本思想是将原方程组转化为递推的形式,并选择一个初始解,通过递推计算得到趋于或精确的解。

常用的迭代法有Jacobi迭代法、Gauss-Seidel迭代法和SOR迭代法等。

这些方法的时间复杂度通常为O(n^2)。

迭代法的优点是适用于大规模问题,但收敛速度慢,精度较差。

二、矩阵特征值求解方法2.1幂法幂法是求解特征值最大的特征值与对应特征向量的方法。

假设有一个n×n的矩阵A,选择一个初始向量x(0),通过迭代计算x(k)=Ax(k-1)/,Ax(k-1),其中,·,表示向量的范数,直到收敛为止。

最后得到的x为特征向量,特征值为λ=(Ax·x)/(x·x)。

幂法的收敛速度较慢,但适用于特征值分布差异较大的情况。

2.2反幂法反幂法是求解特征值最小的特征值与对应特征向量的方法。

和幂法类似,反幂法选择一个初始向量x(0),通过迭代计算x(k)=(A-λI)^-1x(k-1)/,(A-λI)^-1x(k-1),其中I为单位矩阵,λ为近似的特征值,直到收敛为止。

第8章 特征值问题的计算方法

第8章 特征值问题的计算方法

第8章 特征值问题的计算方法本章中讨论求n 阶实矩阵的特征值的数值方法。

8.1 基本概念与性质设A 是n 阶方阵,若数λ和非0向量x 满足:x Ax λ=则λ称为A 的特征值,x 称为A 的对应于λ的特征向量。

A 的特征值的全体()A λ称为A 的谱集。

n 次多项式方程()0det =−A I λ称为A 的特征方程,()A I −λdet 称为A 的特征多项式。

8.2 幂法矩阵的模最大的特征值称为主特征值。

幂法可用于求矩阵的主特征值及其相应的特征向量。

设n 阶方阵A 有有n 个线性无关的特征向量。

设j j j x Ax λ=,j=1..n,其中j λ是A 的特征值,设A 的主特征值1λ是实数且是单重,n λλλ≥≥>L 21.特征向量乘以非0常数仍然是特征向量,故可增加约束,只求范数为1的向量。

设v 0是任意一个非0向量,则v 0可惟一地表示成n 个特征向量的线性组合,设∑==ni i i x v 10α,假设01≠α,令01v A Av v k k k ==−,则111211111~x x x x v k n i i ki i k ni i k i i k αλλλααλλα⎟⎟⎠⎞⎜⎜⎝⎛⎟⎟⎠⎞⎜⎜⎝⎛+==∑∑==,∞→k , 当k>>1时,11−≈k k v v λ,11λ→−k k v v ,1x v v k k →。

为避免计算机出现上溢或下溢现象,在每步计算中将v k 规格化。

111−−≈=k k k v Av u λ,k k k m u v =,, k=1,2,…… 则 1x v k →, ()()()111111,,,−−−−−≈=k k k k k k v v v Av v u λ())1111,,−−−≈k k k k v v v u λ若取2kk u m =(k=0,1,2,…),则()11,−≈k k v u λ,简化了运算。

算法8.1功能:用幂法求矩阵主特征值。

矩阵特征值的求法举例

矩阵特征值的求法举例

矩阵特征值的求法举例矩阵是线性代数中的重要概念,它在科学计算、工程领域以及图像处理等领域都有着广泛的应用。

而在矩阵中,特征值是一个非常重要的概念,它不仅能够描述矩阵的性质,还能够在很多实际问题中起到关键作用。

那么,特征值又是如何求解的呢?本文将通过几个具体的例子来说明矩阵特征值的求法。

一、矩阵特征值的定义我们来介绍一下矩阵的特征值是什么。

对于一个n阶矩阵A(n*n),如果存在一个数λ和一个非零向量v,使得Av=λv,那么我们称λ是矩阵A的特征值,v是对应的特征向量。

特征值和特征向量的求解对于矩阵的性质和应用有着非常重要的作用。

下面我们就通过具体的例子来说明矩阵特征值的求法。

二、特征值的求法1. 对角矩阵的特征值我们来看一个简单的例子,对于一个对角矩阵,特征值的求法非常简单。

对于一个对角矩阵D,我们有D=diag{d1, d2, …, dn},其中对角线元素为d1, d2, …, dn。

那么,对角矩阵的特征值为其对角线元素,即λ1=d1, λ2=d2, …, λn=dn。

特征向量可以取对应的单位向量,如e1=[1, 0, 0, …, 0],e2=[0, 1, 0, …, 0],以此类推。

对于一个2*2的对角矩阵A= [3, 0; 0, 5],其特征值为λ1=3, λ2=5,对应的特征向量可以分别取为v1=[1, 0]和v2=[0, 1]。

接下来,我们来看一个稍复杂一点的例子,对于一个3*3的矩阵,特征值的求法比较繁琐,通常采用特征多项式的方法进行求解。

假设矩阵A= [a, b, c; d, e, f; g, h, i],我们可以先求解其特征多项式:|A-λI| = det|a-λ, b, c; d, e-λ, f; g, h, i-λ|简化上式得到:(a-λ)(e-λ)(i-λ) + (b*d*λ + c*f*λ + a*e*λ) - (a*f*λ + c*d*λ + b*i) = 0然后,我们解出多项式的根,即为矩阵A的特征值。

第七章 求矩阵特征值的数值方法

第七章 求矩阵特征值的数值方法
T
4
定理 2.关于正交矩阵有如下结论
(1)单位矩阵是正交矩阵。 (2)若 A 为正交矩阵,则 A
1
AT 。
(3)若 A 为正交矩阵,则 A 也是正交矩阵。 (4)若 A 为正交矩阵,则 A 1 或 A 1 。 (5)若 A, B 为同阶正交矩阵,则 AB 与 BA 也是正交矩阵。
5
n
k
在 a1 a2 0 时,这是一个主特征向量。
18
② 当 1 2 时, 2 是 A 2 的主特征值,故
1,2 lim max( A2 zk ) .
k
计算时,迭代可能发散,但是
A2 zk Ayk 1 mk 1 Azk 1 mk 1mk 2 zk 2 max( A2 zk ) mk 1mk 2 , mk 1mk 2 1 k 。
17
j a1 x1 a2 x2 a j x j 1 j 3 zk k n j max a1 x1 a2 x2 a j x j 1 j 3 a1 x1 a2 x2 k max(a1 x1 a2 x2 )
P A x am Am x am1 Am1 x ...... a1 Ax a0 x
am x am1
m m 1
x ...... a1 x a0 x P x
7
即, P 是 P A 的一个特征值
定理 5. 为方阵 A 的特征值,若 x1 , x2 都是属于 的特 征向量,则
第七章 求矩阵特征值的数值方法
1
1、预备知识
定义 1 设 A (aij ) nn , 如果 AT A ,则称 A 为对称矩阵。

数值分析——矩阵特征值问题计算

数值分析——矩阵特征值问题计算

17
vk 1k a1x1
即为矩阵 A 的对应特征值 1 的近似特征向量。

vk 1 Avk 1k 1a1x1 1vk
用 (vk)i 表示 vk 的第 i 个分量,则当k充分大时,有
vk1 i
vk
i
1
即为主特征值的近似值。
18
定理 设 A Rnn 有 n 个线性无关的特征向量,
主特征值 1 满足
a11k x1 a2k2 x2 ankn xn
1k
a1
x1
a2
2 1
k
x2
an
n 1k Fra bibliotek xn 1k a1x1 k
16
其中
k
a2
2 1
k
x2
an
n 1
k
xn
由假设条件 从而
j 1
1 j 2, ,n, 所以
lim
k
vk
1k
a1x1
lim
k
k
0
所以当k充分大时,有
vk 1k a1x1
9
0.5, 1, 0.8611
3 5.7222, 11.4444, 8.361 4 5.4621, 10.9223, 8.2306 5 5.5075, 11.0142, 8.2576
11.4444 10.9223 11.0142
0.5, 1, 0.7360 0.5, 1, 0.7536 0.5, 1, 0.7494
n
( Ax, x) (x, x)
1
(2)
n
min x0
( Ax, x) (x, x)
(3)
1
max x0
( Ax, x) (x, x)
  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

H
x x xCni1且x0
H
由定理9.2.3,xCmnia1且xx0
xH Ax xH x
i
又由定理9.2.2,对任意x≠0,有
1
max
xCni1且x 0
xH Ex xH x
n
从而有 i i 1
另一方面, A=(A+E)-E. 记 1 2 E的特征值,那么, i ni1
为n 矩阵-
重复上面的过程,可得 i i 1
第9章 矩阵特征值问题的数值 方法
9.1 特征值与特征向量 9.2 Hermite矩阵特征值问题 9.3 Jacobi方法 9.4 对分法 9.5 乘幂法 9.6 反幂法 9.7 QR方法
9.1 特征值与特征向量
设A是n阶矩阵,x是非零列向量. 如果有
数λ存在,满足
, (1)
那么,称x是矩阵A关于特征值λ的特征向 量.
f () 0的根. 反之,如果λ是
的根,
f () 0
那么齐次方程组(2)有非零解向量x,使(1)式
成立. 从而,λ是A的一个特征值.
A的特征值也称为A的特征根.
矩阵特征值和特征向量有如下主要性质:
定理9.1.1 n阶矩阵A是降秩矩阵的充分必要 条件是A有零特征值.
定理9.1.2 设矩阵A与矩阵B相似,那么它们 有相同的特征值.
从而有 i i n
定理9.2.5通常又称为Hermite矩阵特征值 的扰动定理
பைடு நூலகம்
定理9.2.6 设矩阵A和A′=A+E都是n阶Hermite矩
阵,其特征值分别为 1 2 n 和1 2 ,n 那么 i E 2 2 i E 2
这个定理表明,扰动矩阵E使A的特征值的变化
不会超过 ‖E‖2. 一般‖E‖2 矩阵特征值是良态的.
定理9.1.3 n阶矩阵A与AT有相同的特征值.
定理9.1.4 设λi≠λj是n阶矩阵A的两个互异特 征值,x、y分别是其相应的右特征向 量和左特征向量,那么,xTy=0 .
9.2 Hermite矩阵特征值问题
• 设A为n阶矩阵,其共轭转置矩阵记为AH. 如 果A=AH,那么,A称为Hermite矩阵.
个特征值,那么
A
2
max
1in
i
n
A F
i2
i 1
证:
由 A 2 ( AH A) ( A2 ) (( A))2 2
因此
A
2
max
1in
i n
又由 A 2 tr( AH A) tr( A2 ) F
i2
n
i 1
得 A F
i2
i 1
设x是一个非零向量,A是Hermite矩阵,
称 xH Ax 为矩阵A关于向量x的Rayleigh商, 记为xRH(xx).
矩阵特征值问题的性态是很复杂的,通常 分别就单个特征值或整体特征值给出状态数进 行分 析. 对于Hermite矩阵,由于其特征值问题 的特殊性质,其特征值都是良态的.下面先证明 Hermite矩阵特征值的扰动定理.
定理9.2.5 设矩阵A,E,A+E都是n阶Hermite 矩阵,其特征值分别为1 2 n 1 2 n 1 2 n 那么, i n i i 1
9.2.1 Hermite矩阵的有关性质
设 1, 2 ,..., n 是Hermite矩阵A的n个特征
值. 有以下性质:
• 1, 2 ,..., n全是实数.
• 1, 2 ,..., n有相应的n个线性无关的特征
向量,它们可以化为一组标准酉交的特征
向量组 u1,u2 ,..., un,即 uiHu j ij
特征向量为u1, u2 ,..., un. 用Ck表示酉空间
Cn中任意的k维子空间,那么
k
max min R(x) Ck xCk且x0

k
max min R(x) Cnk1 xCnk1且x0
9.2.3 Hermite矩阵特征值问题的性态
矩阵特征值问题与求解线性方程组问题一 样,都存在当矩阵A的原始 数据有小变化(小扰 动)时,引起特征值问题的变化有大有小的问题, 如果引起的变化小,称 该特征值问题是良态的. 反之,称为病态的.
如果把(1)式右端写为 x,那么(1)式又可写
为:
(I A)x 0 (2)

即| I A | 0
f () | I A | n an1n1 ... a1 a0
它是关于参数λ的n次多项式,称为矩阵A的特 征多项式, 其中a0=(-1)n|A|.
显然,当λ是A的一个特征值时,它必然

• u1, u2 ,..., un是酉空间中的一组标准酉交基.
• 记U=(u1, u2 ,..., un ),它是一个酉阵,即
UHU=UUH=I,那么
1
UH
AU
D
n
即A与以1, 2 ,..., n为对角元的对角阵相似.
• A为正定矩阵的充分必要条件是1, 2 ,..., n
全为正数.
定理9.2.1 设 1, 2 ,..., 是nHermite矩阵A的n
定理9.2.2 如果A的n个特征值为1 2 ... n
其相应的标准酉交的特征向量为 u1,u2 ,..., un
那么有 1 R(x) n
定理9.2.3 设A是Hermite矩阵 ,那么
k
min
xCk且x0
R(
x)或k
min R(x)
xCnk1且x0
9.2.2 极值定理
定理9.2.4(极值定理) 设Hermite矩阵的n个特 征值为 1 2 ... n,其相应的标准酉交
设A
a11 a21
a12 a22
是二阶实对称矩阵,即a21=a12,
其特征值为λ1,λ2. 令
使得
RT
AR
1
1

cos sin
R
sin
cos
B
RT
AR
b11 b21
证 设矩阵A关于特征值λ1,λ2,…,λn 的标准
酉交特征向量为u1,u2,…,un,
Cni1
是由ui,ui+1,…,un生成的n-i+1维子空间.
对 Cni1中任意非零向量x,由极值定理,有
i
max
xCni 1且x 0
xH (A E)x xH x
xH Ax
xH Ex
max
max
x x xCni1且x0
Hermite
9.3 Jacobi方法
理论上,实对称矩阵A正交相似于以A的特征 值为对角元 的 对角阵. 问题是如何构造这 样的正交矩阵呢? Jacobi方法就是通过构造 特殊的正交矩阵 序列,通过相似变换使A 的非对角线元素逐次零化来实现对角化的.
9.3.1 平面旋转矩阵与相似约化
先看一个简单的例子.
相关文档
最新文档