线性代数应用基础补充2
Posted myjs999
tags:
篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了线性代数应用基础补充2相关的知识,希望对你有一定的参考价值。
特征值集合可以多项式运算
若 A A A的特征值集合是 σ ( A ) \\sigma(A) σ(A),则 p ( A ) p(A) p(A)的 特征值集合是 p ( σ ( A ) ) p(\\sigma(A)) p(σ(A))。 p ( A ) = k 0 I + k 1 A + k 2 A 2 + ⋯ + k − 1 A − 1 + k − 2 ( A − 1 ) 2 + ⋯ p(A)=k_0I+k_1A+k_2A^2+\\cdots+k_-1A^-1+k_-2(A^-1)^2+\\cdots p(A)=k0I+k1A+k2A2+⋯+k−1A−1+k−2(A−1)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=1−2wwT是一个豪斯霍尔德(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=v1−v2。
性质:
- P P P对称。
- P P P是幺正矩阵。
- P 2 = I P^2=I P2=I,即 P = P − 1 P=P^-1 P=P−1。
- ∣ P ∣ = 1 |P|=1 ∣P∣=1
- P P P的特征值是 n − 1 n-1 n−1个 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,∃P是Householder⟺∣∣x∣∣2=∣∣y∣∣2也就是经过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
r≡xi2+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/ 以上是关于线性代数应用基础补充2的主要内容,如果未能解决你的问题,请参考以下文章