线性代数应用基础补充2

Posted myjs999

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了线性代数应用基础补充2相关的知识,希望对你有一定的参考价值。

Gershgorin圆盘

格史高林(Gershgorin)圆盘:对于复方阵 A A A的每一行,以 a i i a_ii aii为中心,其它元素的模的和为半径的圆称为圆盘 D i D_i Di

Gershgorin定理 A A A的所有特征值都在所有 D i D_i Di的并集里。

更强的版本:对于每一列,以 a i i a_ii aii为中心,其它元素的模的和为半径的列圆盘为 D i ∗ D^*_i Di。那么 A A A的所有特征值在所有 D i D_i Di的并集里,也在所有 D i ∗ D^*_i Di的并集里(原来是废话)。即 σ ( A ) ⊂ ( ⋃ D i ) ∩ ( ⋃ D i ∗ ) \\sigma(A)\\subset(\\bigcup D_i)\\cap(\\bigcup D_i^*) σ(A)(Di)(Di)

定理 称几个圆盘连成的一个极大连通块是一个连通区域,如果这个区域由 m m m个圆盘构成,那么这里面恰好包含 m m m个特征值。

豪斯霍尔德变换

设实向量 ∣ ∣ w ∣ ∣ = 1 ||\\bmw||=1 w=1,称矩阵 P = 1 − 2 w w T P=1-2\\bmww^T P=12wwT是一个豪斯霍尔德(Householder)矩阵。

将任意向量 v \\bmv v分解为沿 w \\bmw w方向 v 2 \\bmv_2 v2和垂直于 w \\bmw w方向 v 1 \\bmv_1 v1,则有 P v 1 = v 1 Pv_1=v_1 Pv1=v1 P v 2 = − v 2 Pv_2=-v_2 Pv2=v2,因此 P v = v 1 − v 2 P\\bmv=\\bmv_1-\\bmv_2 Pv=v1v2

性质:

  • P P P对称。
  • P P P是幺正矩阵。
  • P 2 = I P^2=I P2=I,即 P = P − 1 P=P^-1 P=P1
  • ∣ P ∣ = 1 |P|=1 P=1
  • P P P的特征值是 n − 1 n-1 n1 1 1 1 1 1 1 − 1 -1 1

定理 x ≠ y \\bmx\\ne \\bmy x=y,则 y = P x , ∃ P 是 H o u s e h o l d e r    ⟺    ∣ ∣ x ∣ ∣ 2 = ∣ ∣ y ∣ ∣ 2 \\bmy=P\\bmx,\\exists P\\rm是Householder\\iff||\\bmx||_2=||\\bmy||_2 y=Px,PHouseholderx2=y2也就是经过Householder变换的向量模长不变。两个模长相等的向量也一定可以由Householder变换互相得到。

Givens变换

吉芬斯(Givens)矩阵 J ( i , k , θ ) J(i,k,\\theta) J(i,k,θ) n n n阶单位矩阵把 i , i i,i i,i k , k k,k k,k两个位置改成 cos ⁡ θ \\cos \\theta cosθ,然后把张出去的右上角改成 sin ⁡ θ \\sin\\theta sinθ、左下角改成 − sin ⁡ θ -\\sin\\theta sinθ得到的矩阵。

吉芬斯矩阵是正交矩阵而且行列式是 1 1 1

对于向量 x i \\x_i\\ xi,如果想要把第 i i i k k k个元素变成 r ≡ x i 2 + x k 2 r\\equiv\\sqrtx_i^2+x_k^2 rxi2+xk2 0 0 0,只需要左乘一个吉芬斯矩阵,其 cos ⁡ θ = x i / r \\cos \\theta=x_i/r cosθ=xi/r sin ⁡ θ = x k / r \\sin \\theta = x_k/r sinθ=xk/r。都是零的特殊情况除外,此时 cos ⁡ θ = 1 \\cos\\theta=1 cosθ=1 sin ⁡ θ = 0 \\sin\\theta=0 sinθ=0

称可对角化矩阵为稀碎矩阵,因为他的所有特征向量都线性无关,Jordan标准型是对角矩阵。

幂迭代法

A A A是实稀碎矩阵,特征向量为 x i \\bmx_i xi,对应 λ i \\lambda_i λi。任意向量 x ( 0 ) = ∑ α i x i \\bmx^(0)=\\sum\\alpha_i\\bmx_i x(0)=αixi

α 1 ≠ 0 \\alpha_1\\ne 0 α1=0,且 ∣ λ 1 ∣ |\\lambda_1| λ1最大,那么 A k x ( 0 ) ( k → + ∞ ) A^k\\bmx^(0)(k\\to+\\infty) A<

以上是关于线性代数应用基础补充2的主要内容,如果未能解决你的问题,请参考以下文章

线性代数应用基础补充1

线性代数应用基础补充1

数学笔记17——定积分的应用2(体积)

hdu6354 /// 圆的相交

数值分析基础应用线性代数

数值分析基础应用线性代数