对角化

相似矩阵

相似:A、B是n阶方阵,存在n阶可逆矩阵P,如果:P1AP=BP^{-1}AP=B,则:A、B相似(A ~ B)
  • 反身性,A ~ A,E1AE=BE^{-1}AE=B
  • 对称性,A ~ B -> B ~ A,同时左乘P以及右乘P1P^{-1},即:PP1APP1=PBP1PP^{-1}APP^{-1}=PBP^{-1},得PBP1=APBP^{-1}=A,转换下:(P1)1B(P1)=A(P^{-1})^{-1}B(P^{-1})=A
  • 传递性,A ~ B 并且 B ~ C -> A ~ C,P1AP=BP^{-1}AP=BQ1BP=CQ^{-1}BP=C,代入B:Q1P1APQ=CQ^{-1}P^{-1}APQ=C(PQ)1A(PQ)=C(PQ)^{-1}A(PQ)=C
性质1:如果A ~ B,那么A、B有相同的特征值并且|A|=|B|(对角线相乘)、tr(A)=tr(B)

证明:λEB=λEP1AP=λP1EPP1AP=p1×λA×P=λA |\lambda E - B| = |\lambda E - P^{-1}AP| = |\lambda P^{-1}EP - P^{-1}AP| = |p^{-1}| \times |\lambda - A| \times |P| = |\lambda - A|

性质2:如果A ~ B,那么A可逆的充要条件是B可逆,并且如果A、B都可逆,那么A1 A^{-1} ~ B1 B^{-1}

证明:根据性质1可知,如果A ~ B,那么|A|=|B|,①|A|、|B|均不等于0,则A、B可逆,②|A|、|B|均等于0,则A、B不可逆,即:A、B同时可逆或同时不可逆

B1=(P1AP)1=P1A1PB^{-1}=(P^{-1}AP)^{-1} = P^{-1}A^{-1}P,所以A1 A^{-1} ~ B1 B^{-1}

性质3:如果A ~ B,那么Am A^m ~ Bm B^m

证明:Bm=(P1AP)m=P1AP×P1AP×P1AP...=P1AmP B^m = (P^{-1}AP)^m = P^{-1}AP \times P^{-1}AP \times P^{-1}AP ... = P^{-1}A^mP

A ~ B 性质总结:
  • tr(A) = tr(B)
  • |A| = |B|
  • 均可逆(或均不可逆)
  • A1=B1A^{-1}=B^{-1}
  • Am=BmA^{m}=B^{m}
  • 特征值相同
  • r(A) = r(B),矩阵左乘可逆矩阵或右乘可逆矩阵或左右都乘可逆矩阵,矩阵的秩不变

计算题(相似求未知量)

A=(20001102a)B=(10003201b)A = \begin{pmatrix} 2&0&0\\ 0&1&1\\ 0&2&a \end{pmatrix} \quad B = \begin{pmatrix} 1&0&0\\ 0&3&2\\ 0&-1&b \end{pmatrix} \quad A ~ B,求:a、b的值

用tr(A) = tr(B) 和 |A| = |B|,解方程即可

对角形

P1APP^{-1}AP = ∧(对角形)

与对角形相似的条件

定理:A ~ ∧ <=> A有n个线性无关的特征向量

证明必要性:P1APP^{-1}AP = ∧ = (λ10000λn)P=(α1,...,αn)\begin{pmatrix} {\lambda_1}&0&0\\ {\vdots}&{\ddots}&{\vdots}\\ 0&0&{\lambda_n} \end{pmatrix} \quad P=(\alpha_1,...,\alpha_n)

P1APP^{-1}AP=∧,同时左乘P,得:AP=P∧,代入向量组,A(α1,...,αn)=(α1,...,αn)(λ10000λn)A(\alpha_1,...,\alpha_n)=(\alpha_1,...,\alpha_n)\begin{pmatrix} {\lambda_1}&0&0\\ {\vdots}&{\ddots}&{\vdots}\\ 0&0&{\lambda_n} \end{pmatrix} = Aα1,...,Aαn=λ1α1,...,λnαnA \alpha_1,...,A \alpha_n = \lambda_1 \alpha_1,...,\lambda_n \alpha_n(根据分块矩阵定义)

所以,Aα1=λ1α1,...,Aαn=λnαn A \alpha_1 = \lambda_1 \alpha_1,...,A \alpha_n = \lambda_n \alpha_n

因为P可逆,所以α1,...,αn!=0\alpha_1,...,\alpha_n != 0

所以,满足特征向量的性质:Aα=λαA \alpha = \lambda \alpha且α != 0

因为,r(P)=n,所以r(α1,...,αn)=n r(\alpha_1,...,\alpha_n) =n(矩阵的秩等于矩阵的列秩),所以α1,...,αn\alpha_1,...,\alpha_n线性无关

证明充分性:把必要性的步骤反过来证明即可

推论:如果A有n个互异的特征根(充分条件),那么A ~ ∧

证明:λ1,...,λn \lambda_1,...,\lambda_n 是互异的话,那么每一个特征值找一个特征向量出来(α1,...,αn \alpha_1,...,\alpha_n ),它们必定是线性无关的(条件强于上面的定理)

  • 如果矩阵所有的根都是单根(必定线性无关),一定相似于对角形
  • 如果矩阵有重根,需要找到n(n等于重根的个数)个线性无关的特征向量,否则不能不相似于对角形

定理:A ~ ∧ <=> 每个rir_i重的特征根,其基础解系有rir_i个解

计算题(对角形相似)

A=(321222361)A = \begin{pmatrix} 3&2&-1\\ -2&-2&2\\ 3&6&-1 \end{pmatrix},是否与对角形相似?如果相似,P=?,∧=?

λEA=λ3212λ+2236λ+1|\lambda E - A| = \begin{vmatrix} {\lambda-3}&-2&1\\ 2&{\lambda+2}&-2\\ -3&-6&{\lambda+1} \end{vmatrix},第3列加到第1列,λ2210λ+22λ26λ+1=(λ2)2(λ+4)=0\begin{vmatrix} {\lambda-2}&-2&1\\ 0&{\lambda+2}&-2\\ {\lambda-2}&-6&{\lambda+1} \end{vmatrix} = (\lambda-2)^2(\lambda+4)=0

λ1=4λ2=λ3=2\lambda_1=-4 \quad \lambda_2=\lambda_3=2

λ1=4\lambda_1=-4时,对(λ1EA)x=0(\lambda_1 E - A )x=0(721222363)>(10130123000)\begin{pmatrix} -7&-2&1\\ 2&-2&-2\\ -3&-6&-3 \end{pmatrix} -> \begin{pmatrix} 1&0&{-\frac{1}{3}}\\ 0&1&{\frac{2}{3}}\\ 0&0&0 \end{pmatrix}

{x1=13x3x2=23x3\begin{cases} x_1=\frac{1}{3}x_3\\ x_2=-\frac{2}{3}x_3 \end{cases}

x3=1x_3=1,得:(13231)=α1\begin{pmatrix} {\frac{1}{3}}\\ {-\frac{2}{3}}\\ 1 \end{pmatrix}=\alpha_1

λ2=λ3=4\lambda_2=\lambda_3=4时,得:(210)=α2(101)=α3\begin{pmatrix} -2\\ 1\\ 0 \end{pmatrix} = \alpha_2 \quad \begin{pmatrix} 1\\ 0\\ 1 \end{pmatrix} = \alpha_3

P=(α1,α2,α3)=(13212310101) P=(\alpha_1,\alpha_2,\alpha_3) = \begin{pmatrix} {\frac{1}{3}}&-2&1\\ {-\frac{2}{3}}&1&0\\ 1&0&1 \end{pmatrix}

∧ = (400020002)\begin{pmatrix} -4&0&0\\ 0&2&0\\ 0&0&2 \end{pmatrix}(把特征值依次带入即可)

注意:特征值和特征向量的顺序必须是对应的

计算题(相似反推矩阵)

3阶方阵A,λ1=1,λ2=2,λ3=3,α1=(1,1,1)T,α2=(1,2,3)T,α3=(1,3,6)T\lambda_1=1 , \lambda_2=2 , \lambda_3=3 , \alpha_1=(1,1,1)^T , \alpha_2=(1,2,3)^T , \alpha_3=(1,3,6)^T,求:A、ATA^T的特征值和特征向量

P=(α1,α2,α3)=(111123136)P1AP=(100020003)P=(\alpha_1,\alpha_2,\alpha_3)=\begin{pmatrix} 1&1&1\\ 1&2&3\\ 1&3&6 \end{pmatrix} \quad P^{-1}AP=\begin{pmatrix} 1&0&0\\ 0&2&0\\ 0&0&3 \end{pmatrix}=∧

同时左乘P,右乘P1P^{-1},得:A=P∧P1P^{-1},得:A=(010012397)A=\begin{pmatrix} 0&1&0\\ 0&-1&2\\ 3&-9&7 \end{pmatrix}

ATA^T的特征值不变,所以还是λ1=1,λ2=2,λ3=3\lambda_1=1 , \lambda_2=2 , \lambda_3=3

P1AP=P^{-1}AP=∧,(P1AP)T=PTAT(p1)T=(P^{-1}AP)^T=P^TA^T(p^{-1})^T=∧(对角形的转置等于它本身)

逆矩阵的转置等于转置的逆矩阵,所以PTAT(PT)1=P^TA^T(P^T)^{-1}=

整理为定理:[(PT)1]1AT(PT)1=[(P^T)^{-1}]^{-1}A^T(P^T)^{-1}=∧,求得:(PT)1=(331352121)(P^T)^{-1}=\begin{pmatrix} 3&-3&1\\ -3&5&-2\\ 1&-2&1 \end{pmatrix}

(PT)1(P^T)^{-1}3个列向量就是ATA^T的特征向量

计算题(利用对角形求高次幂)

A=(011202110)A=\begin{pmatrix} 0&1&-1\\ -2&0&2\\ -1&1&0 \end{pmatrix},求:A100A^{100}

先求出特征值,λ1=1,λ2=0,λ3=1\lambda_1=-1 , \lambda_2=0 , \lambda_3=1

再求出特征向量并构成矩阵P,P=(111014113)P=\begin{pmatrix} 1&1&1\\ 0&1&4\\ 1&1&3 \end{pmatrix}

P1AP=P^{-1}AP=∧,同时左乘P和右乘P1P^{-1},得:A=P∧P1P^{-1}

P∧P1P^{-1}P∧P1P^{-1}...P∧P1P^{-1} = P∧的100次方P1P^{-1}

计算出结果...

实对称矩阵的对角化(预热)

所有的实对称矩阵都能对角化(也就是说:所有的实对称矩阵都能相似于一个对角形)

内积:两个向量对应的分量相乘后累加,得到的值就是内积(内积是个数)

α=(a1a2a3)β=(b1b2b3)\alpha=\begin{pmatrix} a_1\\ a_2\\ a_3 \end{pmatrix} \quad \beta=\begin{pmatrix} b_1\\ b_2\\ b_3 \end{pmatrix}

(α,β) = a1b1+a2b2+a3b3=numa_1b_1 + a_2b_2 + a_3b_3 = num

我们习惯上用"行"乘"列"的方式来表示内积。如:α、β都是列向量,那么可以表示为:αTβ\alpha^T \beta
  • (α,α) >= 0,因为:a12+a22+a32a_1^2+a_2^2+a_3^2肯定大于等于0,所以:(α,α) = 0 <=> α = 0
  • 对称性:(α,β) = (β,α)
  • 齐次性:(kα,β)=k(α,β),(α,kβ)=k(α,β),(kα,kβ)=k2k^2(α,β)
  • (α+β,γ)=(α,γ)+(β,γ), (α,β+γ)=(α,β)+(α,γ)
向量的长度(或者说:范数、模),α=(α,α)||\alpha||=\sqrt{(\alpha,\alpha)},两边同时平方:(α,α)=α2(\alpha,\alpha)=||\alpha||^2

α=(1,1,5)\alpha=(-1,-1,5),那么α=(1)2+(1)2+(5)2=27||\alpha|| = \sqrt{(-1)^2+(-1)^2+(5)^2} = \sqrt{27},点到原点的距离

如果||α||=1,则α是单位向量

单位化(标准化),1αα=β\frac{1}{||\alpha||}\alpha=\beta,乘以倒数后的向量的长度为1

证明:β=(1αα,1αα)=1α2(α,α)=1α(α,α)=1αα=1||\beta|| = \sqrt{(\frac{1}{||\alpha||}\alpha,\frac{1}{||\alpha||}\alpha)}=\sqrt{\frac{1}{||\alpha||^2}(\alpha,\alpha)}=\frac{1}{||\alpha||}\sqrt{(\alpha,\alpha)}=\frac{1}{||\alpha||}||\alpha||=1

性质1:α>=0 ||\alpha||>=0,||α||=0 <=> α=0
性质2:||kα||=|k| * ||α||(kα的长度=k的绝对值乘以α的长度)

证明:kα=(kα,kα)=k2(α,α)=k×α||k\alpha||=\sqrt{(k\alpha,k\alpha)}=\sqrt{k^2(\alpha,\alpha)}=|k| \times ||\alpha||

柯西施瓦茨不等式,|(α,β)| <= ||α|| * ||β||

a1b1+a2b2+...+anbn<=a12+a22+...+an2×b12+b22+...+bn2|a_1b_1+a_2b_2+...+a_nb_n| <= \sqrt{a_1^2+a_2^2+...+a_n^2} \times \sqrt{b_1^2+b_2^2+...+b_n^2}

三角不等式,||α+β|| <= ||α|| + ||β||,两边之和大于第三边,在向量中如果同方向就可以等于

正交

正交(垂直),(α,β) = 0,那么 α ⊥ β,0向量和任何向量都正交,和自身正交的向量只能是0向量
正交向量组:给一组向量(不含0向量),两两都正交,α1,α2,...,αs\alpha_1,\alpha_2,...,\alpha_s中随便取2个向量都是正交的
标准正交向量组,每一个向量都是单位向量,{(αi,αi)=1(αi,αj)=0\begin{cases} (\alpha_i,\alpha_i)=1\\ (\alpha_i,\alpha_j)=0 \end{cases}
定理:如果α1,α2,...,αs\alpha_1,\alpha_2,...,\alpha_s是正交向量组,那么α1,α2,...,αs\alpha_1,\alpha_2,...,\alpha_s一定线性无关(线性相关须在同一直线上)

证明: 设:k1α1+k2α2+...+ksαs=0k_1\alpha_1+k_2\alpha_2+...+k_s\alpha_s=0

(α1,k1α1+k2α2+...+ksαs)=(α1,0)=0(\alpha_1,k_1\alpha_1+k_2\alpha_2+...+k_s\alpha_s)=(\alpha_1,0)=0

k1(α1,α1)+k2(α1,α2)+...+ks(α1,αs)=0k_1(\alpha_1,\alpha_1) + k_2(\alpha_1,\alpha_2) + ... + k_s(\alpha_1,\alpha_s) = 0

根据正交的性质,后面的都等于0,只剩:k1(α1,α1)=0k_1(\alpha_1,\alpha_1)=0,因为正交向量组不含0向量,所以(α1,α1)!=0(\alpha_1,\alpha_1) != 0,所以:k1=0k_1 = 0

用同样的方式,求出k2,...,ks=0k_2,...,k_s=0,所以线性无关

施密特正交化,给一组线性无关的α1,α2,...,αs\alpha_1,\alpha_2,...,\alpha_s,求一组与之等价的并且是正交的β1,β2,...,βs\beta_1,\beta_2,...,\beta_s

β1=α1\beta_1 = \alpha_1

β2=α2(α2,β1)(β1,β1)×β1\beta_2 = \alpha_2 - \frac{(\alpha_2,\beta_1)}{(\beta_1,\beta_1)} \times \beta_1

β3=α3(α3,β1)(β1,β1)×β1(α3,β2)(β2,β2)×β2\beta_3 = \alpha_3 - \frac{(\alpha_3,\beta_1)}{(\beta_1,\beta_1)} \times \beta_1 - \frac{(\alpha_3,\beta_2)}{(\beta_2,\beta_2)} \times \beta_2

...

计算题(施密特正交化)

线性无关的向量,α1=(1,1,1,1)Tα2=(1,2,3,4)Tα3=(1,2,2,3)T\alpha_1=(1,1,1,1)^T \quad \alpha_2=(1,-2,-3,-4)^T \quad \alpha_3=(1,2,2,3)^T

一般流程:α1,α2,...αs\alpha_1,\alpha_2,...\alpha_s 正交化-> β1,β2,...βs\beta_1,\beta_2,...\beta_s 单位化-> η1,η2,...ηs\eta_1,\eta_2,...\eta_s

β1=α1\beta_1 = \alpha_1

β2=α2(α2,β1)(β1,β1)β1\beta_2 = \alpha_2 - \frac{(\alpha_2,\beta_1)}{(\beta_1,\beta_1)}\beta_1

β3=α3(α3,β1)(β1,β1)β1(α3,β2)(β2,β2)β2\beta_3 = \alpha_3 - \frac{(\alpha_3,\beta_1)}{(\beta_1,\beta_1)}\beta_1 - \frac{(\alpha_3,\beta_2)}{(\beta_2,\beta_2)}\beta_2

η1=1β1β1η2=1β2β2η3=1β3β3\eta_1 = \frac{1}{||\beta_1||}\beta_1 \quad \eta_2 = \frac{1}{||\beta_2||}\beta_2 \quad \eta_3 = \frac{1}{||\beta_3||}\beta_3

正交矩阵

正交矩阵:A是n阶方阵,ATA=EA^TA=E
性质1:如果A是正交矩阵,那么|A| = 1或-1

证明:AT×A=1 |A^T| \times |A| = 1 ,则:A2=1 |A|^2=1

性质2:如果A是正交矩阵,那么A1=ATA^{-1}=A^TA1A^{-1}ATA^T均为正交矩阵

证明:因为ATA=EA^TA=E,所以A1=ATA^{-1}=A^T

(A1)TA1=(A1)TAT=(AA1)T=E(A^{-1})^TA^{-1}=(A^{-1})^TA^T=(AA^{-1})^T=E

(AT)TAT=AAT=AA1=E(A^T)^TA^T=AA^T=AA^{-1}=E

性质3:A、B正交,那么AB也正交

证明:(AB)TAB=BTATAB=B1A1AB=E(AB)^TAB=B^TA^TAB=B^{-1}A^{-1}AB=E

性质4:如果A是正交矩阵,α、β是n维列向量,那么(Aα,Aβ)=(α,β)

证明:(Aα)TAβ=αTATAβ=αTβ=(α,β)(A\alpha)^TA\beta=\alpha^TA^TA\beta=\alpha^T\beta=(\alpha,\beta)

定理:A是正交矩阵 <=> A的列(行)向量组是标准正交向量组

证明:设A的列向量为(α1,α2,...,αn)(\alpha_1,\alpha_2,...,\alpha_n)

向量可以看作是一个分块矩阵,分块矩阵转置:先把子块元素求转置,再对每个子块求转置

ATA=(α1,α2,...,αn)T(α1,α2,...,αn) A^TA = (\alpha_1,\alpha_2,...,\alpha_n)^T(\alpha_1,\alpha_2,...,\alpha_n)

=(α1Tα2TαnT)(α1,α2,...,αn) =\begin{pmatrix} {\alpha_1^T}\\ {\alpha_2^T}\\ {\vdots}\\ {\alpha_n^T} \end{pmatrix} (\alpha_1,\alpha_2,...,\alpha_n)

=(α1Tα1α1Tα2α1Tαnα2Tα1α2Tα2α2TαnαnTα1αnTα2αnTαn)=\begin{pmatrix} {\alpha_1^T\alpha_1}&{\alpha_1^T\alpha_2}&{\cdots}&{\alpha_1^T\alpha_n}\\ {\alpha_2^T\alpha_1}&{\alpha_2^T\alpha_2}&{\cdots}&{\alpha_2^T\alpha_n}\\ {\vdots}&{\vdots}&{\ddots}&{\vdots}\\ {\alpha_n^T\alpha_1}&{\alpha_n^T\alpha_2}&{\cdots}&{\alpha_n^T\alpha_n} \end{pmatrix}

=((α1α1)(α1α2)(α1αn)(α2α1)(α2α2)(α2αn)(αnα1)(αnα2)(αnαn))=\begin{pmatrix} {(\alpha_1\alpha_1)}&{(\alpha_1\alpha_2)}&{\cdots}&{(\alpha_1\alpha_n)}\\ {(\alpha_2\alpha_1)}&{(\alpha_2\alpha_2)}&{\cdots}&{(\alpha_2\alpha_n)}\\ {\vdots}&{\vdots}&{\ddots}&{\vdots}\\ {(\alpha_n\alpha_1)}&{(\alpha_n\alpha_2)}&{\cdots}&{(\alpha_n\alpha_n)} \end{pmatrix}

因为是正交矩阵,所以最终结果等于单位阵E,即:(100010001)\begin{pmatrix} 1&0&{\cdots}&0\\ 0&1&{\cdots}&0\\ {\vdots}&{\vdots}&{\ddots}&{\vdots}\\ 0&0&{\cdots}&1 \end{pmatrix}

(αi,αi)=α2=1(\alpha_i,\alpha_i)=||\alpha||^2=1(ai,aj)=0(a_i,a_j)=0

所以,满足必要性(标准正交向量组)

充分性的证明逻辑跟必要性类似

计算题(正交矩阵)

三阶非0矩阵A,aij=Aija_{ij}=A_{ij},求:①|A|=1,②A是正交矩阵

证明:A=(A11A21A31A12A22A32A13A23A33)=(a11a21a31a12a22a32a13a23a33)=ATA^*=\begin{pmatrix} {A_{11}}&{A_{21}}&{A_{31}}\\ {A_{12}}&{A_{22}}&{A_{32}}\\ {A_{13}}&{A_{23}}&{A_{33}} \end{pmatrix}=\begin{pmatrix} {a_{11}}&{a_{21}}&{a_{31}}\\ {a_{12}}&{a_{22}}&{a_{32}}\\ {a_{13}}&{a_{23}}&{a_{33}} \end{pmatrix}=A^T

ATA=AA=AEA^TA=A^* A=|A|E

两边同时取行列式

AT×A=A|A^T | \times |A| = | |A| |

  • AT=A |A^T|=|A|
  • 行列式提取k要提n次

所以,A2=A3 |A|^2 = |A|^3

以上式子可以写成:A2(A1)=0|A|^2(|A|-1)=0,有2个解:1和0

因为A是三阶非0矩阵,所以必定有元素!=0,假设:a11!=0a_{11} != 0

那么,按第一行展开后,得:A=a11A11+a12A12+a13A13=a112+a122+a132>0 |A|=a_{11}A_{11} + a_{12}A_{12} + a_{13}A_{13} = a_{11}^2+a_{12}^2+a_{13}^2 > 0aij=Aija_{ij}=A_{ij}

所以,|A| != 0,只有一个解:|A| = 1

实对称矩阵的对角化

实数并且是对称矩阵的对角化

n阶实对称矩阵的n个特征值都是实数并且特征向量都是实向量

定理1:实对称矩阵A的不同特征值对应的特征向量一定正交

证明:AT=A A^T=A

λ1,λ2(λ1!=λ2)\lambda_1,\lambda_2 \quad (\lambda_1 != \lambda_2)

Aα1=λ1α1Aα2=λ2α2A\alpha_1 = \lambda_1\alpha_1 \quad A\alpha_2 = \lambda_2\alpha_2

只要证出内积(α1,α2)=0 (\alpha_1,\alpha_2) = 0 ,就表示正交

(Aα1,α2)=(λ1α1,α2)=λ1(α1,α2)(A\alpha_1,\alpha_2)=(\lambda_1\alpha_1,\alpha_2)=\lambda_1(\alpha_1,\alpha_2)

(Aα1,α2)=(Aα1)Tα2=α1TATα2=α1TAα2=λ2α1Tα2=λ2(α1,α2)(A\alpha_1,\alpha_2)=(A\alpha_1)^T\alpha_2=\alpha_1^TA^T\alpha_2=\alpha_1^TA\alpha_2=\lambda_2\alpha_1^T\alpha2=\lambda_2(\alpha_1,\alpha_2)

所以,λ1(α1,α2)=λ2(α1,α2)\lambda_1(\alpha_1,\alpha_2) = \lambda_2(\alpha_1,\alpha_2)

所以,(λ1λ2)(α1,α2)=0(\lambda_1-\lambda_2)(\alpha_1,\alpha_2)=0

因为,λ1!=λ2\lambda_1 != \lambda_2

所以,(α1,α2)=0(\alpha_1,\alpha_2) = 0,所以特征向量一定正交

正交相似:如果A、B是同阶方阵,如果存在正交矩阵P,使得 P1AP=BP^{-1}AP=B
定理2:假设A是实对称矩阵,一定存在正交矩阵Q,使得:Q1AQ=Q^{-1}AQ=∧=(λ1000λ2000λn)\begin{pmatrix} {\lambda_1}&0&{\cdots}&0\\ 0&{\lambda_2}&{\cdots}&0\\ {\vdots}&{\vdots}&{\ddots}&{\vdots}\\ 0&0&{\cdots}&{\lambda_n} \end{pmatrix}
这个定理说明:n阶实对称矩阵A一定有n个线性无关的特征向量
推论:假设A是实对称矩阵,λi\lambda_i是它的一个rir_i重特征根,那么A对应于rir_i的对应λi\lambda_i的特征向量一定有rir_i个(实对称矩阵一定能对角化)

计算步骤总结

3阶实对称矩阵A,求正交矩阵Q,使得Q1AQ=Q^{-1}AQ=

证明:①求特征值 ②求特征向量 ③特征向量正交化、单位化 ④特征向量做成列构成Q、特征值构成∧

注意:对角形的λ\lambda顺序要跟特征向量Q的λ\lambda一致

情况1,设:λ1=a,λ2=b,λ3=c\lambda_1=a,\lambda_2=b,\lambda_3=c(三个根都是单根)

根据上面的定理1:α1,α2,α3\alpha_1,\alpha_2,\alpha_3一定正交,直接单位化即可:η1=1α1α1,η2=1α2α2,η3=1α3α3\eta_1=\frac{1}{||\alpha_1||}\alpha_1,\eta_2=\frac{1}{||\alpha_2||}\alpha_2,\eta_3=\frac{1}{||\alpha_3||}\alpha_3

Q=(η1,η2,η3)Q = (\eta_1,\eta_2,\eta_3)

Q1AQ=(λ1000λ2000λ3)Q^{-1}AQ=\begin{pmatrix} {\lambda_1}&0&0\\ 0&{\lambda_2}&0\\ 0&0&{\lambda_3} \end{pmatrix}

情况2,设:λ1=a,λ2=λ3=b\lambda_1=a,\lambda_2=\lambda_3=b(一个单根、一个二重根)

α2,α3\alpha_2,\alpha_3做施密特正交化,α1\alpha_1不用做施密特正交化,因为一定是正交的,但是为了符号统一,可以写成β1=α1\beta_1 = \alpha_1

再单位化后,得到η1,η2,η3\eta_1,\eta_2,\eta_3

情况3,设:λ1=λ2=λ3=a\lambda_1=\lambda_2=\lambda_3=a(一个三重根)

α1,α2,α3\alpha_1,\alpha_2,\alpha_3都得做施密特正交化,得到:β1,β2,β3\beta_1,\beta_2,\beta_3

再单位化后,得到η1,η2,η3\eta_1,\eta_2,\eta_3

results matching ""

    No results matching ""