矩阵对角化

一、什么样的n阶矩阵才能对角化?

先说结论 {a.特征值无重根b.有重根{b.1.厄米特矩阵(Hermitian),特别的对称实数矩阵b.2.非厄矩阵{b.2.1.几何重数=代数重数b.2.2.几何重数≠代数重数

定理1:

λ1,λ2,,λm\lambda_1,\lambda_2,\dots,\lambda_m 是方阵AA的m个特征值,p1,p2,,pmp_1,p_2,\dots,p_m依次是与之对应的特征向量,如果λ1,λ2,,λm\lambda_1,\lambda_2,\dots,\lambda_m各不相等,则p1,p2,,pmp_1,p_2,\dots,p_m线性无关。

定理2:

n阶矩阵AA与对角矩阵相似(即AA能对角化)的充分必要条件是AAnn个线性无关的特征向量。

推论:

如果n阶矩阵AAnn个特征值互不相等,则AA与对角矩阵相似。

二、根据以上定理和推理特征值无重根的a和b.2.1情况一定可以对角化。

三、有重根情况b.1情况

3.1 厄米特矩阵

厄米特矩阵是满足AH=AA^H=A的矩阵,即AHA^HAA的共轭转置矩阵。

3.2 酉矩阵

满足AHA=IA^HA=I的矩阵。酉矩阵,其列向量是规范正交的,其秩必为n。 A1=AHA^{-1}=A^H n阶n维线性无关向量必可以通过施密特正交化得到标准正交基矩阵,即是酉矩阵。

3.3 舒尔化矩阵

设n阶矩阵AA的特征值λk\lambda_kkk重根,对应mm个特征值,p1,p2,,pmp_1,p_2,\dots,p_m依次是与之对应的特征向量,组成一组基,把{p1,p2,,pm}\{p_1,p_2,\dots,p_m\}扩充至n个(n维的最大线性无关组是n个),并施密特正交化,得到P=(p1,p2,,pm,pm+1,,pn)P=(p_1,p_2,\dots,p_m,p_{m+1},\dots,p_n) AP=(λkp1,λkp2,,λkpm,Apm+1,Apm+2,,Apn)AP=(\lambda_kp_1,\lambda_kp_2,\dots,\lambda_kp_m,Ap_{m+1},Ap_{m+2},\dots,Ap_n)

PHAP=(p1Hp2Hp3HpmHpm+1Hpn)(λkp1,λkp2,,λkpm,Apm+1,Apm+2,,Apn) P^HAP= \begin{pmatrix} p_1^H \\ p_2^H\\ p_3^H\\ \vdots\\ p_m^H\\ p_{m+1}^H\\ \vdots\\ p_{n} \end{pmatrix}(\lambda_kp_1,\lambda_kp_2,\dots,\lambda_kp_m,Ap_{m+1},Ap_{m+2},\dots,Ap_n)

A=(aij),pj=(qij)A=(a_{ij}), p_j=(q_{ij}) 因为PHP^H的正交性,所以,

PHAP=(λkp1HApm+1p1HApnλkp2HApm+1p2HApnλkpmHApm+1pmHApn000pm+1HApm+1pm+1HApn000pnHApm+1pnHApn) \begin{aligned} P^HAP &=\begin{pmatrix} \lambda_{k}& & & & p^H_1Ap_{m+1}&\cdots&p^H_1Ap_{n} \\ & \lambda_{k} & & & p^H_2Ap_{m+1}&\cdots&p^H_2Ap_{n}\\ & & \ddots & & \vdots& &\vdots \\ & & & \lambda_{k} & p^H_mAp_{m+1}&\cdots&p^H_mAp_{n}\\ 0 & 0 & \cdots & 0 & p^H_{m+1}Ap_{m+1}&\cdots&p^H_{m+1}Ap_{n} \\\vdots& \vdots& & \vdots& \vdots& &\vdots \\0&0& \cdots & 0&p^H_nAp_{m+1}&\cdots&p^H_nAp_{n}\\ \end{pmatrix} \end{aligned} P为我们构造的酉矩阵,于是\RightarrowP1AP=PHAP=(λkp1HApm+1p1HApnλkp2HApm+1p2HApnλkpmHApm+1pmHApn000pm+1HApm+1pm+1HApn000pnHApm+1pnHApn)=B \begin{equation} \begin{aligned} P^{-1}AP&=P^HAP &=\begin{pmatrix} \lambda_{k}& & & & p^H_1Ap_{m+1}&\cdots&p^H_1Ap_{n} \\ & \lambda_{k} & & & p^H_2Ap_{m+1}&\cdots&p^H_2Ap_{n}\\ & & \ddots & & \vdots& &\vdots \\ & & & \lambda_{k} & p^H_mAp_{m+1}&\cdots&p^H_mAp_{n}\\ 0 & 0 & \cdots & 0 & p^H_{m+1}Ap_{m+1}&\cdots&p^H_{m+1}Ap_{n} \\\vdots& \vdots& & \vdots& \vdots& &\vdots \\0&0& \cdots & 0&p^H_nAp_{m+1}&\cdots&p^H_nAp_{n}\\ \end{pmatrix} \\ &=B \end{aligned} \end{equation} 为书写方便,将(1)式BB简写为块矩阵:

B=(λkλkλk000000)=(λkWλkOM)

其中MM矩阵可以重复上面的步骤,最终得到上三角矩阵。

AABB相似,根据有关定理,它们具有相同的特征值。 |AλE|=|BλE||A-\lambda E|=|B-\lambda E|

|BλE|=(λkλWλkλOMλE)=(λkλ)mg(λ)

又因为λk\lambda_{k}kk重根,所以|AλE|=(λkλ)kf(λ)|A-\lambda E|=(\lambda_{k}-\lambda)^kf(\lambda) ,所以,几何重数mm必小于代数重数kk , 即mkm\leq k

又因为PP为酉矩阵,所以 BH=(PHAP)H=(PH(AP))H=(AP)H(PH)H=(AP)HP=PHAHP \begin{aligned} B^H&=(P^HAP)^H \\ &=(P^H(AP))^H \\ &=(AP)^H(P^H)^H \\ &=(AP)^HP \\ &=P^HA^HP \\ \end{aligned}

如果AA为厄米特矩阵,则BH=BB^H=B,所以BB也为厄米特矩阵。所以BB必为对角矩阵。 所以,b.1情况一定可以对角化。

参考:
1. 厄米特矩阵
2. 舒尔定理
3. 几何重数小于等于代数重数

百度大模型

发表回复

您的邮箱地址不会被公开。 必填项已用 * 标注

腾讯云图