线性代数二三事

第4章:作用于内积空间的变换

4.3 矩阵的正交对角化

第4章 作用于内积空间的变换

回顾我们在第三章中提到的矩阵的对角化: 对于矩阵A\vec A而言, 其可以对角化的条件是存在对角矩阵P\vec P与可逆矩阵Q\vec Q, 使得

A=QPQ1\vec A = \vec Q\vec P\vec Q^{-1}

成立. 其中位于P\vec P主对角线上的元素即为A\vec A的特征值; Q\vec Q的列向量即为A\vec A的特征向量. 如果我们回顾上一节中提到的正交矩阵和酉矩阵的定义, 我们发现当Q\vec Q为正交矩阵时, Q=Q1\vec Q^\top = \vec Q^{-1}. 因此我们能否尝试构造出A=QPQ\vec A = \vec Q\vec P\vec Q^\top的形式呢?

定义 4.8

AMn(F)\vec A \in M_{n}(\mathbb{F}).

① 若F=R\mathbb{F}=\mathbb{R}, 且存在由Rn\mathbb{R}^n的一组单位正交基为列向量构成的矩阵Q\vec Q以及对角矩阵P\vec P, 使得A=QPQ\vec A = \vec Q\vec P\vec Q^\top, 我们称A\vec A正交对角化(Orthogonally Diagonalizable);

② 若F=C\mathbb{F}=\mathbb{C}, 且存在由Cn\mathbb{C}^n的一组单位正交基为列向量构成的矩阵Q\vec Q以及对角矩阵P\vec P, 使得A=QPQ\vec A = \vec Q\vec P\vec Q^* (Q=(Q)\vec Q^* = (\overline{\vec Q})^\top), 我们称A\vec A酉对角化(Unitarily Diagonalizable).

通过定义, 我们可以分析出使得矩阵A\vec A正交对角化的条件:

  • 首先, A\vec A必须可以对角化. 否则我们就无法构造出矩阵P\vec PQ\vec Q.
  • 其次, 由于Q\vec Q为正交(酉)矩阵, 这也就代表了Q\vec Q中的列向量彼此正交. 由于我们知道, 不同的特征向量所对应的特征空间满足直和的关系, 因此A\vec A的特征空间彼此正交, 即EigA(λi)EigA(λj),ijEig_{\vec A}(\lambda_i) \perp Eig_{\vec A}(\lambda_j), i\neq j;
  • 最后, 在一个特征空间中, 我们同样要满足其基底向量彼此正交, 因此我们需要通过Gram-Schmidt算法去找出特征空间的单位正交基.

例题 4.5

设线性变换T:R2R2\vec T: \mathbb{R}^2 \to \mathbb{R}^2满足

T(11)=(22);T(10)=(30).\vec T\begin{pmatrix} 1 \\ 1 \end{pmatrix} = \begin{pmatrix} 2 \\ 2 \end{pmatrix}; \quad \vec T\begin{pmatrix} 1 \\ 0 \end{pmatrix} = \begin{pmatrix} 3 \\ 0 \end{pmatrix}.

那么令E:={(10),(01)}E := \left\{ \begin{pmatrix} 1 \\ 0 \end{pmatrix}, \begin{pmatrix} 0 \\ 1 \end{pmatrix} \right\}, 判断[T]EE[\vec T]_E^E能否正交对角化.

解答 4.5

不难发现, A:=[T]EE=(3102)\vec A := [\vec T]_E^E=\begin{pmatrix} 3 & -1 \\ 0 & 2 \end{pmatrix}. 那么很显然A\vec A有两个不同的特征值: λ1=3;λ2=2\lambda_1=3;\lambda_2=2. 由此我们知道A\vec A可对角化. 对于λ1=3\lambda_1=3而言, 其特征向量v1\vec v_1满足

v1Ker(Aλ1I)=Ker(0101)=Span(10).\vec v_1 \in \Ker(\vec A-\lambda_1\vec I) = \Ker\begin{pmatrix} 0 &-1\\0 & -1\end{pmatrix} = \Span\begin{pmatrix} 1\\0 \end{pmatrix}.

对于λ2=2\lambda_2=2而言, 其特征向量v2\vec v_2满足

v2Ker(Aλ2I)=Ker(1100)=Span(11).\vec v_2 \in \Ker(\vec A-\lambda_2\vec I) = \Ker\begin{pmatrix} 1 & -1 \\ 0 & 0 \end{pmatrix} = \Span\begin{pmatrix} 1\\ 1 \end{pmatrix}.

我们发现, v1,v2\vec v_1 , \vec v_2所处的特征空间不正交. 因此我们无法构造出正交矩阵Q\vec Q, 因此[T]EE[T]_E^E不可正交对角化.

定理 4.8

AMn(F)\vec A \in M_n(\mathbb{F}).

① 若F=R\mathbb{F}=\mathbb{R}, 且A\vec A可正交对角化, 那么A=A\vec A=\vec A^\top (A\vec A为对称矩阵);

② 若F=C\mathbb{F}=\mathbb{C}, 且A\vec A可酉对角化, 那么A\vec A正规矩阵 (Normal Matrix), 即 AA=AA\vec A\vec A^*=\vec A^*\vec A.

证明

对于①而言, 我们知道此时存在对角矩阵P\vec P和正交矩阵Q\vec Q, 使得A=QPQ\vec A = \vec Q \vec P \vec Q^\top. 那么

A=(QPQ)=(Q)PQ=QPQ.\vec A^\top = (\vec Q\vec P\vec Q^\top)^\top = (\vec Q^\top)^\top\vec P^\top\vec Q^\top = \vec Q\vec P\vec Q^\top.

对于②而言, 类似地

AA=QPQ(QPQ)=QPQQPQ=QPPQ=QPPQ=QPQQPQ=AA.\vec A\vec A^* = \vec Q\vec P\vec Q^*(\vec Q\vec P\vec Q^*)^*=\vec Q\vec P\vec Q^*\vec Q\vec P^*\vec Q^* = \vec Q\vec P\vec P^*\vec Q^*=\vec Q\vec P^*\vec P\vec Q^*=\vec Q\vec P^*\vec Q^*\vec Q\vec P\vec Q^* = \vec A^*\vec A.

该定理说明,实矩阵若可正交对角化,则必须是对称矩阵; 反过来, 我们随后便会知道:当AMn(R)\vec A \in M_n(\mathbb{R})时任何对称矩阵均可正交对角化. 对于AMn(C)\vec A \in M_n(\mathbb{C})而言, 我们把正规矩阵主要分成以下几类:

• 若AMn(C)\vec A\in M_n(\mathbb{C})满足A=A\vec A = \vec A^*, 则A\vec A自伴矩阵 (Self-adjoint Matrices). 其中

x,Ay=(x)(Ay)=(xA)y=(Ax)y=Ax,y=Ax,y.\la \vec x,\vec A\vec y \ra = (\vec x)^\top (\overline{\vec A \vec y}) = (\vec x^\top \overline{\vec A})\overline{\vec y} = (\overline{\vec A}^\top\vec x)^\top\vec y = \la \vec A^*\vec x,\vec y\ra = \la \vec A\vec x,\vec y\ra.

若定义在内积空间VV上的线性变换T:VV\vec T: V \to V满足Tu,v=u,Tv\la \vec T\vec u, \vec v \ra = \la \vec u, \vec T\vec v \ra, T\vec T也被称作自伴算子 (Self-adjoint Operator).

定理 4.9

AMn(C)\vec A \in M_n( \mathbb{C})为自伴矩阵, 则A\vec A的所有特征值为实数.

证明

λ\lambdaA\vec A的一个特征值, v\vec v为与之对应的特征向量. 那么根据推论4.24.2, 当AMn(C)\vec A \in M_n(\mathbb{C})时我们可以得出相似的结果:

x,Ay=(x)(Ay)=(xA)y=(Ax)y=Ax,y.\la \vec x,\vec A\vec y \ra = (\vec x)^\top (\overline{\vec A \vec y}) = (\vec x^\top \overline{\vec A})\overline{\vec y} = (\overline{\vec A}^\top\vec x)^\top\vec y = \la \vec A^*\vec x,\vec y\ra.

λ\lambdaA\vec A的一个特征值, v\vec v为与之对应的特征向量, 那么我们有 Av,v=v,Av\la \vec A\vec v, \vec v \ra = \la \vec v, \vec A\vec v \ra, 因此由Av=λv\vec A\vec v = \lambda\vec v, 我们得到λv,v=v,λv\la \lambda\vec v,\vec v \ra = \la \vec v,\lambda\vec v \ra. 根据内积的性质, λv,v=λv2\la \lambda\vec v , \vec v \ra = \lambda \cdot|| \vec v ||^2; v,Av=λv2\la \vec v, \vec A\vec v \ra = \overline{\lambda} \cdot ||\vec v||^2. 即我们有λv2=λv2\lambda ||\vec v||^2 = \overline{\lambda} ||\vec v||^2. 由于v\vec v为非零向量, 则λR\lambda \in \mathbb{R}.

同理, 所有的自伴算子的特征值也均为实数.

• 若AMn(C)\vec A\in M_n(\mathbb{C})满足A=A\vec A = -\vec A^*, 则A\vec A斜自伴矩阵 (Skew-adjoint Matrices). 其中

x,Ay=(x)(Ay)=(xA)y=(Ax)y=Ax,y=Ax,y.\la \vec x,\vec A\vec y \ra = (\vec x)^\top (\overline{\vec A \vec y}) = (\vec x^\top \overline{\vec A})\overline{\vec y} = (\overline{\vec A}^\top\vec x)^\top\vec y = \la \vec A^*\vec x,\vec y\ra = -\la \vec A\vec x,\vec y\ra.

若定义在内积空间VV上的线性变换T:VV\vec T: V \to V满足Tu,v=u,Tv\la \vec T\vec u, \vec v \ra = -\la \vec u, \vec T\vec v \ra, T\vec T也被称作斜自伴算子 (Skew-adjoint Operator).

定理 4.10

AMn(C)\vec A \in M_n( \mathbb{C})为斜自伴矩阵, 则A\vec A的所有非零特征值为纯虚数.

证明

λ\lambdaA\vec A的一个特征值, v\vec v为与之对应的特征向量. 则 Av,v=λv,v=λv2.\la \vec A\vec v, \vec v \ra = \la \lambda\vec v, \vec v \ra = \lambda ||\vec v||^2. 同时 Av,v=v,Av=v,λv=λv2.\la \vec A\vec v, \vec v \ra = -\la \vec v, \vec A\vec v \ra = -\la \vec v, \lambda\vec v \ra = -\overline{\lambda} ||\vec v||^2.λ=λ\lambda = -\overline{\lambda}, 即Re(λ)=0\mathbf{Re}(\lambda)=0.

同理, 所有的斜自伴算子的非零特征值也均为纯虚数.

• 若AMn(C)\vec A\in M_n(\mathbb{C})满足A=A1\vec A^* = \vec A^{-1}, 则A\vec A酉矩阵 (Unitary Matrices). 其中

Ax,Ay=(Ax)(Ay)=A(xA)y=A(Ax)y=AAx,y=AAx,y=x,y.\la \vec A\vec x , \vec A\vec y \ra = (\vec A\vec x)^\top (\overline{\vec A\vec y}) = \vec A(\vec x^\top \overline{\vec A})\overline{\vec y} = \vec A(\vec A^*\vec x)^\top\overline{\vec y} = \vec A \la \vec A^*\vec x, \vec y \ra = \vec A\vec A^* \la \vec x,\vec y\ra = \la \vec x,\vec y \ra.

若定义在内积空间VV上的线性变换T:VV\vec T: V \to V满足Tu,Tv=u,v\la \vec T\vec u, \vec T\vec v \ra = \la \vec u, \vec v \ra, T\vec T也被称作酉算子 (Unitary Operator).

定理 4.11

AMn(C)\vec A \in M_n( \mathbb{C})为酉矩阵, 则对A\vec A中任意的特征值λ\lambda而言, λλ=1\lambda \overline{\lambda}=1.

证明

根据酉算子的定义, Av,Av=v,v.\la \vec A\vec v,\vec A\vec v \ra = \la \vec v, \vec v \ra. 此外Av,Av=λv,λv=λλv,v\la \vec A\vec v, \vec A\vec v \ra = \la \lambda\vec v, \lambda \vec v \ra = \lambda \cdot \overline{\lambda} \la \vec v, \vec v \ra. 因此原式得证.

同理, 所有的酉算子的特征值也均满足上述性质.

我们接下来将重点研究正规矩阵的性质:

定理 4.12

AMn(C)\vec A \in M_n(\mathbb{C})为正规矩阵, u\vec uA\vec A中与λ\lambda对应的特征值, 则u\vec u同为A\vec A^*中与λ\overline{\lambda}对应的特征值.

证明

可知Av,v=λv,v=λv,v=v,λv\la \vec A\vec v,\vec v \ra = \la \lambda\vec v, \vec v \ra = \lambda \la \vec v, \vec v \ra = \la \vec v, \overline{\lambda} \vec v \ra. 同时结合Av,v=v,Av\la \vec A\vec v, \vec v \ra = \la \vec v, \vec A^*\vec v \ra即可得证.

定理 4.13

AMn(C)\vec A \in M_n(\mathbb{C})为正规矩阵, u,v\vec u, \vec v分别为特征值λ,μ\lambda, \mu所对应的特征向量(λμ\lambda \neq \mu). 则u,v=0\la \vec u, \vec v \ra = 0.

证明

Au=λu,Av=μv\vec A\vec u = \lambda\vec u, \vec A\vec v=\mu\vec v. 那么Au,v=λu,v=λu,v\la \vec A\vec u, \vec v \ra = \la \lambda\vec u, \vec v \ra = \lambda \la \vec u, \vec v \ra. 同时

Au,v=u,Av=u,μv=μu,v.\la \vec A\vec u, \vec v \ra = \la \vec u, \vec A^*\vec v \ra = \la \vec u, \overline{\mu}\vec v \ra = \mu \la \vec u, \vec v \ra.

由于λμ\lambda \neq \mu, 则u,v=0\la \vec u, \vec v \ra = 0.

也就是说, 对于所有的正规矩阵或正规算子而言, 其特征空间两两彼此正交. 那么是不是就代表了这些矩阵可以正交对角化了呢? 我们现在提出本章最重要的定理:

定理 4.14

VV为定义在数域F\mathbb{F} (R\mathbb{R}C\mathbb{C})上的内积空间, 记线性变换T:VV\vec T: V \to V, 其矩阵形式为A\vec A. 假设T\vec T的特征多项式可在F\mathbb{F}中分解, 那么T\vec T可正交对角化(A\vec A可正交对角化)的充要条件是T\vec T为正规算子(A\vec A为正规矩阵).

该定理的证明我们在本书中忽略. 感兴趣的同学请在今后选修泛函分析或者高等代数课程. 我们从两个方面去理解谱定理:

• 第一, 如果F=R\mathbb{F}=\mathbb{R}, 那么A=A\vec A^*=\vec A^\top, 即AA=AA\vec A\vec A^\top = \vec A^\top\vec A, 即A=A\vec A= \vec A^\top, 因此A\vec A为对称矩阵. 也就是说A\vec A可以正交对角化的充要条件是A\vec A为对称矩阵. 此时, 存在矩阵Q\vec Q, 使得Q\vec Q的列向量为A\vec A的特征向量, 且它们构成VV的一组标准正交基. 记P\vec P为包含A\vec A的特征值的对角矩阵, 那么

A=QPQ.\vec A = \vec Q\vec P \vec Q^\top.

• 第二, 如果F=C\mathbb{F}=\mathbb{C}, 此时A\vec A为正规矩阵, 也就包括了我们先前提出的几种情况: 自伴矩阵, 斜自伴矩阵, 酉矩阵. 因此, A\vec A可以正交对角化的充要条件是A\vec A为正规矩阵. 此时, 存在矩阵Q\vec Q, 使得Q\vec Q的列向量为A\vec A的特征向量, 且它们构成VV的一组标准正交基. 记P\vec P为包含A\vec A的特征值的对角矩阵, 那么

A=QPQ.\vec A = \vec Q\vec P \vec Q^*.

例题 4.6

A=(0i0i00000)\vec A = \begin{pmatrix} 0 & i & 0 \\ -i & 0 & 0 \\ 0 & 0 & 0 \end{pmatrix}. 找出P,QM3(C)\vec P,\vec Q\in M_3(\mathbb{C}), 使得A=QPQ\vec A = \vec Q\vec P\vec Q^*.

解答 4.6

不难发现, A\vec A为正规矩阵, 因此A\vec A可正交对角化. A\vec A的特征多项式为

CA(x)=det(xi0ix000x)=x(x21)=x(x+1)(x1).C_{\vec A}(x) = \det\begin{pmatrix} -x & i & 0 \\ -i & -x & 0 \\ 0& 0 &-x \end{pmatrix} = -x(x^2-1) = -x(x+1)(x-1).

因此A\vec A有三个特征值λ1=0,λ2=1,λ3=1\lambda_1=0,\lambda_2=-1,\lambda_3=1. 三个特征值对应的特征空间分别为

EigA(λ1)=Ker(0i0i00000)=Span(001).\Eig_{\vec A}(\lambda_1) = \Ker\begin{pmatrix} 0 & i & 0 \\ -i & 0 & 0 \\ 0 & 0 & 0 \end{pmatrix} = \Span\begin{pmatrix} 0 \\ 0 \\ 1\end{pmatrix}.EigA(λ2)=Ker(1i0i10001)=Span(i10).\Eig_{\vec A}(\lambda_2) = \Ker\begin{pmatrix} 1 & i & 0 \\ -i & 1 & 0 \\ 0 & 0 & 1 \end{pmatrix} = \Span\begin{pmatrix} i \\ -1 \\ 0\end{pmatrix}.EigA(λ3)=Ker(1i0i10001)=Span(1i0).\Eig_{\vec A}(\lambda_3) = \Ker\begin{pmatrix} 1 & i & 0 \\ -i & 1 & 0 \\ 0 & 0 & 1 \end{pmatrix} = \Span\begin{pmatrix} 1 \\ i \\ 0\end{pmatrix}.

由于A\vec A的特征空间两两正交, 因此我们得到的三个向量即构成了C3\mathbb{C}^3的一组正交基. 通过单位化, 我们得到的单位正交基为

v1=(001);v2=(i2120),v3=(12i20).\vec v_1=\begin{pmatrix} 0 \\ 0 \\ 1\end{pmatrix}; \vec v_2= \begin{pmatrix} \frac{i}{\sqrt{2}} \\ -\frac{1}{\sqrt{2}} \\ 0\end{pmatrix}, \vec v_3=\begin{pmatrix} \frac{1}{\sqrt{2}} \\ \frac{i}{\sqrt{2}} \\ 0\end{pmatrix}.

因此我们得到

P=(000010001),Q=(0i212012i2100).\vec P =\begin{pmatrix} 0 & 0 & 0 \\ 0 & -1 & 0 \\ 0 & 0 & 1 \end{pmatrix}, \quad\vec Q = \begin{pmatrix} 0 & \frac{i}{\sqrt{2}} & \frac{1}{\sqrt{2}} \\ 0 & -\frac{1}{\sqrt{2}} & \frac{i}{\sqrt{2}} \\ 1 & 0 & 0 \end{pmatrix}.

在本节的最后, 我们浅谈一下对角化概念的扩充: 三角化 (Triangularization).

定义 4.9

V\vec V为维数为nn的线性空间, 设线性变换T:VV\vec T: \vec V \to \vec V. 若存在V\vec V的一组基底β\beta, 使得[T]β[\vec T]_\beta为上三角矩阵, 且位于主对角线上的元素为T\vec T的特征值, 则我们称T\vec T为可三角化的线性算子.

矩阵的三角化和线性变换的三角化具有相似的定义: 若A\vec A相似于以A\vec A特征值为主对角线元素的上三角矩阵, 我们则称矩阵A\vec A可三角化, 反之同样成立: 此时存在可逆矩阵Q\vec Q和上三角矩阵P\vec P, 使得A=QPQ1\vec A = \vec Q \vec P \vec Q^{-1}. 其中

P=(λ1a12a1nλ2a2nλn),\vec P = \begin{pmatrix} \lambda_1 & a_{12} & \cdots & a_{1n} \\ & \lambda_2 & \cdots & a_{2n} \\ & & \ddots & \\ & & & \lambda_n \end{pmatrix},

λ1,,λn\lambda_1,\cdots,\lambda_nA\vec A的特征值. 如果Q\vec Q为正交矩阵或酉矩阵, 原结论同样成立, 我们称A\vec A可正交三角化.

定义 4.10

A,PMn(F)\vec A,\vec P\in M_n(\mathbb{F}).

① 若F=R\mathbb{F}=\mathbb{R}, 且存在正交矩阵Q\vec Q, 使得A=QPQ\vec A = \vec Q\vec P \vec Q^\top, 我们称A\vec A正交相似于P\vec P;

② 若F=C\mathbb{F}=\mathbb{C}, 且存在酉矩阵Q\vec Q, 使得A=QPQ\vec A = \vec Q\vec P \vec Q^*, 我们称A\vec A酉相似于P\vec P.

下面的定理告诉了我们矩阵A\vec A在什么情况下可以正交三角化:

定理 4.15

V\vec V为维数为nn, 定义在数域F\mathbb{F} (R\mathbb{R}C\mathbb{C})上的内积空间. 则以下两个推论等价:

① 线性变换T:VV\vec T: V \to V的特征多项式可在F\mathbb{F}中分解, 即

CT(x)=(xλ1)a1(xλk)ak,a1++ak=n;λiF;C_{\vec T}(x) = (x-\lambda_1)^{a_1}\cdots(x-\lambda_k)^{a_k},\quad a_1+\cdots+a_k=n; \lambda_i \in \mathbb{F};

② 存在一组VV的单位正交基β\beta, 使得[T]β[\vec T]_\beta为上三角矩阵, 且[T]β[\vec T]_\beta的主对角线上的元素为T\vec T的特征值.

证明

\to ①: 假设存在满足条件的单位正交基β\beta. 那么设T\vec T在任意基底下的矩阵形式为A\vec A, 则存在由β\beta中的元素为列向量的矩阵Q\vec Q, 使得A=Q[T]βQ\vec A = \vec Q [\vec T]_\beta \vec Q^*, 即A,[T]β\vec A, [\vec T]_\beta相似, 又因为相似矩阵的特征值相同, 故①得证.

\to ②. 我们将采用数学归纳法. 假设①成立. 设dim(V)=n\dim(V)=n. 首先若n=1n=1, 那么结果显然成立. 我们假设结论②对所有dim(V)<n\dim(V)<n的正整数成立, 则我们需验证结论②对dim(V)=n\dim(V)=n时成立. 设λ1\lambda_1T\vec T的一个特征值, v1\vec v_1为与之对应的特征向量, 且不失一般性地讲令v1=1||\vec v_1||=1. 利用Gram-Schmidt算法, 以v1\vec v_1为基础构造出V\vec V的一组单位正交基β={v1,,vn}\beta=\{\vec v_1,\cdots,\vec v_n\}. 令P=(v1vn)\vec P=\begin{pmatrix} \vec v_1 & \cdots & \vec v_n \end{pmatrix}, 则

A=P[T]βP1.\vec A = \vec P [\vec T]_\beta\vec P^{-1}.

其中, 我们可以将[T]β[\vec T]_\beta写作

[T]β=([T(v1)]β[T(vn)]β)=([λ1v1]β[T(vn)]β)=(λ1u0n1A1),[\vec T]_\beta = \begin{pmatrix} [\vec T(\vec v_1)]_\beta & \cdots & [\vec T(\vec v_n)]_\beta \end{pmatrix} = \begin{pmatrix} [\lambda_1\vec v_1]_\beta & \cdots & [\vec T(\vec v_n)]_\beta \end{pmatrix} = \begin{pmatrix} \lambda_1 & \vec u^\top \\ \vec 0_{n-1} & \vec A_1 \end{pmatrix},

0n1\vec 0_{n-1}代表维数为n1n-1的零向量; u\vec u为一未知向量, A1\vec A_1为未知(n1)×(n1)(n-1)\times(n-1)矩阵. 由于A,[T]β\vec A, [\vec T]_\beta相似, 因此二者有相同的特征多项式, 因此[T]β[\vec T]_\beta的特征多项式可在F\mathbb{F}中分解. 那么根据数学归纳法原理, 矩阵A1\vec A_1可正交三角化. 即存在正交矩阵(或酉矩阵)QMn1(F)\vec Q\in M_{n-1}(\mathbb{F}), 由A1\vec A_1的特征值为主对角线元素的上三角矩阵A2\vec A_2, 使得

A1=QA2Q1.\vec A_1 = \vec Q \vec A_2\vec Q^{-1}.

我们现在假设

R=(10n10n1Q),\vec R=\begin{pmatrix} 1 & \vec 0_{n-1}^\top \\ \vec 0_{n-1} & \vec Q \end{pmatrix},

那么不难发现矩阵R\vec R可逆. 除此之外, 注意到

R1[T]βR=(10n10n1Q1)(λ0u0n1A1)(10n10n1Q)=(λ0u0n1A2).\vec R^{-1} [\vec T]_\beta \vec R = \begin{pmatrix} 1 & \vec 0_{n-1}^\top \\ \vec 0_{n-1} & \vec Q^{-1} \end{pmatrix} \begin{pmatrix} \lambda_0 & \vec u^\top \\ \vec 0_{n-1} & \vec A_1 \end{pmatrix} \begin{pmatrix} 1 & \vec 0_{n-1}^\top \\ \vec 0_{n-1} & \vec Q \end{pmatrix} = \begin{pmatrix} \lambda_0 & \vec u^\top \\ \vec 0_{n-1} & \vec A_2 \end{pmatrix}.

因此R1[T]βR\vec R^{-1} [\vec T]_\beta \vec R为上三角矩阵. 同时由 R1[T]βR=R1P1APR=(PR)1A(PR)\vec R^{-1}[\vec T]_\beta\vec R = \vec R^{-1}\vec P^{-1} \vec A \vec P\vec R = (\vec P \vec R)^{-1} \vec A (\vec P \vec R) 可知, (PR)1A(PR)(\vec P \vec R)^{-1} \vec A (\vec P \vec R)为上三角矩阵. 由于P,Q\vec P,\vec Q均为正交矩阵(或酉矩阵), 则PQ\vec P\vec Q同为正交矩阵(或酉矩阵). 因此根据定义T\vec T可正交三角化.

推论 4.3

对任意的AMn(C)\vec A \in M_n(\mathbb{C})而言, A\vec A单一相似(Unitarily Similar) 于一个上三角矩阵P\vec P. 即存在一组Cn\mathbb{C}^n的单位正交基{v1,,vn}\{\vec v_1,\cdots,\vec v_n\}, 设矩阵Q=(v1vn)\vec Q = \begin{pmatrix} \vec v_1 & \cdots & \vec v_n \end{pmatrix}, 有

A=QPQ.\vec A = \vec Q \vec P \vec Q^*.

其中P\vec P位于主对角线上的元素即为A\vec A的特征值.

{4.3 练习}

1. 对于下面定义在内积空间VV上的线性变换T:VV\vec T: V \to V, 判断T\vec T是否为: 正规算子; 自伴算子; 酉算子 (当数域为R\mathbb{R}时, 我们把A\vec A^*当成A\vec A^\top处理即可).

(i) V=R2,T(ab)=(2a2b2a+5b)V=\mathbb{R}^2, \vec T\begin{pmatrix} a \\ b \end{pmatrix} = \begin{pmatrix} 2a-2b \\ -2a+5b \end{pmatrix};

(ii) V=R3,T(abc)=(a+b5b4a2b+5c)V =\mathbb{R}^3, \vec T\begin{pmatrix} a \\ b \\ c \end{pmatrix} = \begin{pmatrix} -a+b \\ 5b \\ 4a-2b+5c \end{pmatrix};

(iii) V=C2,T(ab)=(2a+iba+2b)\vec V=\mathbb{C}^2, \vec T\begin{pmatrix} a \\ b \end{pmatrix} = \begin{pmatrix} 2a+ib \\ a+2b \end{pmatrix};

(iv) V=P2(R),T(f(x))=f(x)\vec V = P_2(\mathbb{R}), \vec T(f(x)) = f'(x). 其中f(x),g(x)=01f(x)g(x)dx\la f(x), g(x) \ra = \int_0^1 f(x)\cdot g(x) dx;

(v) V=M2(R),T(A)=AV = M_2(\mathbb{R}), \vec T(\vec A) = \vec A^\top.

2.A=(111111111)\vec A = \begin{pmatrix} 1 & 1 & 1 \\ 1 & 1 & 1 \\ 1 & 1 & 1 \end{pmatrix}, 求出矩阵Q,PM3(R)\vec Q, \vec P\in M_3(\mathbb{R}), 使得A=QPQ\vec A = \vec Q\vec P\vec Q^\top.

3.A=(022201210)\vec A = \begin{pmatrix} 0 & -2 & 2 \\ 2 & 0 & -1 \\ -2 & 1 & 0 \end{pmatrix}, 求出矩阵Q,PM3(C)\vec Q, \vec P\in M_3(\mathbb{C}), 使得A=QPQ\vec A = \vec Q\vec P\vec Q^*.

4. 设定义在维数有限的内积空间VV上的线性变换T:VV\vec T:V \to V为正规算子(即TT=TT\vec T^* \vec T = \vec T \vec T^*, 或[T]β[\vec T]_\beta为正规矩阵). 证明:

(i) 对任意的vV,Tv=Tv\vec v \in V, || \vec T \vec v|| = ||\vec T^* \vec v ||;

(ii) 对任意的λF\lambda \in \mathbb{F}, TλI\vec T - \lambda \vec I为正则算子.

5.VV为定义在C\mathbb{C}上的内积空间, T:VV\vec T: V \to V为线性变换. 定义T1,T2:VV\vec T_1,\vec T_2: V\to V, 使得

T1=12(T+T);T2=12i(TT).\vec T_1 = \frac{1}{2}(\vec T+\vec T^*);\quad\vec T_2=\frac{1}{2i}(\vec T-\vec T^*).

证明: T1,T2\vec T_1,\vec T_2为自伴算子.

6. [简化版谱定理]VV为定义在F\mathbb{F}上的内积空间, T:VV\vec T: V \to V为线性变换. 记T\vec T的矩阵形式为A\vec A. 证明: T\vec T为自伴算子的充要条件是矩阵A\vec A正交相似(或单一相似)于对角矩阵. 即存在正交矩阵(或酉矩阵) Q\vec Q与对角矩阵P\vec P, 使得A=QPQ1\vec A = \vec Q \vec P \vec Q^{-1}成立. 提示: 用Schur定理. 同时该定理还告诉我们如下性质:

  • F=R\mathbb{F} = \mathbb{R}, 那么A\vec A正交相似于对角矩阵的充要条件是A=A\vec A = \vec A^\top;
  • F=C\mathbb{F} = \mathbb{C}, 那么A\vec A单一相似于对角矩阵的充要条件是A=A\vec A = \vec A^*.