有限实矩阵乘法法则的证明
设 m , n , p , q m,n,p,q m , n , p , q 为正整数,A ∈ R m × n A\in\mathbb R^{m\times n} A ∈ R m × n 、B ∈ R n × p B\in\mathbb R^{n\times p} B ∈ R n × p 、C ∈ R p × q C\in\mathbb R^{p\times q} C ∈ R p × q 。以下只用行乘列定义与实数的有限求和。对任意输出位置 ( i , ℓ ) (i,\ell) ( i , ℓ ) ,
( ( A B ) C ) i ℓ = ∑ k = 1 p ( ∑ j = 1 n a i j b j k ) c k ℓ = ∑ j = 1 n a i j ( ∑ k = 1 p b j k c k ℓ ) = ( A ( B C ) ) i ℓ . \begin{aligned}
((AB)C)_{i\ell}
&=\sum_{k=1}^{p}\left(\sum_{j=1}^{n}a_{ij}b_{jk}\right)c_{k\ell}\\
&=\sum_{j=1}^{n}a_{ij}\left(\sum_{k=1}^{p}b_{jk}c_{k\ell}\right)
=(A(BC))_{i\ell}.
\end{aligned} (( A B ) C ) i ℓ = k = 1 ∑ p ( j = 1 ∑ n a ij b jk ) c k ℓ = j = 1 ∑ n a ij ( k = 1 ∑ p b jk c k ℓ ) = ( A ( BC ) ) i ℓ .
中间等号只交换两个有限和的求和次序,并使用实数分配律。两边同为 m × q m\times q m × q 矩阵,每个位置相等,便证明结合律;没有调换矩阵因子的次序。
对 D ∈ R n × p D\in\mathbb R^{n\times p} D ∈ R n × p ,同一定义给出
[ A ( B + D ) ] i k = ∑ j = 1 n a i j ( b j k + d j k ) = ( A B ) i k + ( A D ) i k . [A(B+D)]_{ik}=\sum_{j=1}^{n}a_{ij}(b_{jk}+d_{jk})
=(AB)_{ik}+(AD)_{ik}. [ A ( B + D ) ] ik = j = 1 ∑ n a ij ( b jk + d jk ) = ( A B ) ik + ( A D ) ik .
对 E ∈ R m × n E\in\mathbb R^{m\times n} E ∈ R m × n ,则有
[ ( A + E ) B ] i k = ∑ j = 1 n ( a i j + e i j ) b j k = ( A B ) i k + ( E B ) i k . [(A+E)B]_{ik}=\sum_{j=1}^{n}(a_{ij}+e_{ij})b_{jk}
=(AB)_{ik}+(EB)_{ik}. [( A + E ) B ] ik = j = 1 ∑ n ( a ij + e ij ) b jk = ( A B ) ik + ( EB ) ik .
所以在所写的相容大小下,两条分配律均成立。对任意实标量 c c c ,
[ ( c A ) B ] i k = ∑ j ( c a i j ) b j k = c ∑ j a i j b j k = ∑ j a i j ( c b j k ) = [ A ( c B ) ] i k , [(cA)B]_{ik}=\sum_j(ca_{ij})b_{jk}
=c\sum_j a_{ij}b_{jk}
=\sum_j a_{ij}(cb_{jk})=[A(cB)]_{ik}, [( c A ) B ] ik = j ∑ ( c a ij ) b jk = c j ∑ a ij b jk = j ∑ a ij ( c b jk ) = [ A ( c B ) ] ik ,
逐项相等便证明 ( c A ) B = c ( A B ) = A ( c B ) (cA)B=c(AB)=A(cB) ( c A ) B = c ( A B ) = A ( c B ) 。特别地,对 x , y ∈ R n x,y\in\mathbb R^n x , y ∈ R n 与实数 s , t s,t s , t ,
[ A ( s x + t y ) ] i = ∑ j a i j ( s x j + t y j ) = s ( A x ) i + t ( A y ) i . [A(sx+ty)]_i=\sum_j a_{ij}(sx_j+ty_j)
=s(Ax)_i+t(Ay)_i. [ A ( s x + t y ) ] i = j ∑ a ij ( s x j + t y j ) = s ( A x ) i + t ( A y ) i .
这证明了比较同一线性系统的解时所用的矩阵乘向量线性性质。所有恒等式都只针对大小相容的有限实矩阵,不声称任意矩阵因子可交换或可消去。
向量空间是在域 F F F 上配备加法和数乘,并满足通常的阿贝尔群及分配律公理的集合。列向量只是一个模型。次数至多为 d d d 的多项式空间 P d ( F ) P_d(F) P d ( F ) 和空间之间的线性映射也是模型。零多项式属于 P d ( F ) P_d(F) P d ( F ) ,虽然其次数通常不定义。
子空间 包含零,且对加法和数乘封闭。若每个向量都是某列表中向量的线性组合,该列表张成空间;若只有零系数组合能得到零,它线性无关。基 同时满足两者。有限生成列表可通过删除相关项缩减成基;线性无关列表可通过添加其张成空间外的向量扩充成基,直到张成给定的有限生成列表。
交换引理。 若 r r r 个向量的线性无关列表位于 s s s 个向量的张成空间中,则 r ≤ s r\le s r ≤ s 。理由是:用生成列表表示第一个无关向量,利用某个非零系数,以它替换一个生成向量,再重复。每一步的线性无关性保证下一个无关向量对尚未被替换的某个生成向量有非零系数,因此至少需要 r r r 个位置。对两个基双向比较得到长度相等;共同长度为 dim V \dim V dim V 。
因此在 n n n 维空间中,任意 n n n 个线性无关向量构成基,任意 n n n 个张成向量也构成基。只有在证明线性无关性或张成性之后,计数列表才有用。
不预设计数结论的有限基存在与扩充
设实向量空间 V V V 已给定有限生成列表 t 1 , … , t N t_1,\ldots,t_N t 1 , … , t N 。先证明任何线性无关列表 u 1 , … , u k u_1,\ldots,u_k u 1 , … , u k 的长度不超过 N N N 。插入 j j j 个向量后,维持一个由 u 1 , … , u j u_1,\ldots,u_j u 1 , … , u j 与 N − j N-j N − j 个剩余旧向量组成的生成列表。将 u j + 1 u_{j+1} u j + 1 按此列表展开,至少一个剩余旧向量的系数非零;否则它在先前 u i u_i u i 的张成空间中,与无关性矛盾。解出这个旧向量,再以 u j + 1 u_{j+1} u j + 1 取代它。新列表仍张成同一空间,而且保留全部先前插入向量。若 k > N k>N k > N ,则 N N N 步后前 N N N 个向量已张成 V V V ,与 u N + 1 u_{N+1} u N + 1 的无关性矛盾。此即所需上界,也包括 N = 0 N=0 N = 0 。
每个子空间 W W W 都有有限基。从空的无关列表开始;若当前张成空间尚不是 W W W ,就在 W W W 中取一个不在该张成空间内的向量。若新向量在某关系中的系数非零,就能由旧向量表示它,矛盾,故新列表仍无关。已证上界允许至多添加 N N N 次,所以过程必停止于 W W W 的无关生成列表。若 W = { 0 } W=\{0\} W = { 0 } ,立即得到空基。改从 W W W 内任意给定无关列表开始,同样证明该列表可扩充为 W W W 的基;只需有界有限次选择。最后,若两组基的长度为 p , q p,q p , q ,分别把它们用作上界证明中的有限生成列表,得到 p ≤ q p\le q p ≤ q 与 q ≤ p q\le p q ≤ p 。故基的长度有良定义,称为维数。这些论证补足秩–零度定理所用的基存在、扩充与维数事实。
映射 T : V → W T:V\to W T : V → W 若满足 T ( α v + β u ) = α T v + β T u T(\alpha v+\beta u)=\alpha Tv+\beta Tu T ( αv + β u ) = α T v + βT u ,则为线性。其核与像都是子空间。选取 ker T \ker T ker T 的基 k 1 , … , k r k_1,\ldots,k_r k 1 , … , k r ,并以 v 1 , … , v s v_1,\ldots,v_s v 1 , … , v s 扩充为基,则 T v 1 , … , T v s Tv_1,\ldots,Tv_s T v 1 , … , T v s 构成 im T \operatorname{im}T im T 的基:展开任意输入可证明张成性;若这些像线性相关,相应输入组合落入核,而扩充基的线性无关性迫使所有系数为零。因此
dim V = dim ker T + dim im T . \dim V=\dim\ker T+\dim\operatorname{im}T. dim V = dim ker T + dim im T .
这就是秩–零度定理 。有限维且维数相等时,单射与满射等价;无限维空间中该等价可能失效。
本节默认标量域是实数域 R \mathbb{R} R :除非特别说明,矩阵、标量和向量都
含有实数元素。所有论证都是代数论证;把 R \mathbb{R} R 换成复数域
C \mathbb{C} C 后,同样的推理仍然成立,但可用的特征值可能更多。
单个特征向量已经有用;若一个方阵拥有足够多线性无关的特征向量,就可
以找到一个坐标系,使矩阵的作用变成对角的。这里的“足够多”有精确定义:
对于 n × n n\times n n × n 矩阵,需要 n n n 个线性无关的特征向量,它们才会形成
R n \mathbb{R}^n R n 的一组基底。此后矩阵的幂与逆矩阵可以转化为对角线上逐项
计算。
相似就是换坐标
坐标转换矩阵的构造见基底扩充与换基 。
设 S S S 可逆,令 x = S y x=Sy x = S y ,其中 y y y 是向量在由 S S S 的列组成的新基底下的坐标。
记 B = S − 1 A S B=S^{-1}AS B = S − 1 A S ,则
A x = A S y = S B y . Ax=ASy=SBy. A x = A S y = SB y .
所以 A A A 作用于原坐标中的 x x x ,正好对应 B B B 作用于新坐标中的 y y y 。
这也说明为什么等式左边必须有 S − 1 S^{-1} S − 1 :先由 y y y 转回 x = S y x=Sy x = S y ,再应用
A A A ,最后用 S − 1 S^{-1} S − 1 把输出转回新坐标,得到的正是 S − 1 A S y S^{-1}ASy S − 1 A S y 。矩阵
S S S 必须可逆,否则 y ↦ S y y\mapsto Sy y ↦ S y 不能提供坐标系统;相似允许改变坐标
标签,却不允许没有依据地把一个矩阵替换成另一个矩阵。
定义
相似 同一标量域上的两个 n × n n\times n n × n 矩阵 A A A 与 B B B ,若存在可逆矩阵 S S S 使
S − 1 A S = B , S^{-1}AS=B, S − 1 A S = B , 则称 A A A 与 B B B 相似 。
定义
对角化与可对角化 设 A A A 是 n × n n\times n n × n 矩阵。若存在可逆矩阵 S S S 与标量
λ 1 , … , λ n \lambda_1,\ldots,\lambda_n λ 1 , … , λ n ,使
S − 1 A S = D = diag ( λ 1 , … , λ n ) , S^{-1}AS=D=\operatorname{diag}(\lambda_1,\ldots,\lambda_n), S − 1 A S = D = diag ( λ 1 , … , λ n ) , 则称这个等式是 A A A 的一个对角化 ,并称 A A A 可对角化 。
对角化是目标矩阵为对角矩阵的特殊相似关系。对角线上的数值只有在与
S S S 的列保持同一次序时才有意义。
特征向量字典
定理
一个对角化等式的刻画 设 A A A 是 n × n n\times n n × n 矩阵,且
S = [ v 1 v 2 ⋯ v n ] S=[v_1\ v_2\ \cdots\ v_n] S = [ v 1 v 2 ⋯ v n ] 是可逆矩阵。令 D = diag ( λ 1 , … , λ n ) D=\operatorname{diag}(\lambda_1,\ldots,\lambda_n) D = diag ( λ 1 , … , λ n ) 。以下
叙述等价:
对每个 j = 1 , … , n j=1,\ldots,n j = 1 , … , n ,都有 A v j = λ j v j Av_j=\lambda_jv_j A v j = λ j v j ;
A S = S D AS=SD A S = S D ;
S − 1 A S = D S^{-1}AS=D S − 1 A S = D ;
A = S D S − 1 A=SDS^{-1} A = S D S − 1 。
证明
四种等价形式及其依赖关系 A S AS A S 的第 j j j 列是 A v j Av_j A v j ,S D SD S D 的第 j j j 列是 λ j v j \lambda_jv_j λ j v j ,所以所有
列方程合在一起恰好就是 A S = S D AS=SD A S = S D 。由于 S S S 已经声明可逆,左乘 S − 1 S^{-1} S − 1
得到 S − 1 A S = D S^{-1}AS=D S − 1 A S = D ;另一种做法是直接对原来的 A S = S D AS=SD A S = S D 右乘 S − 1 S^{-1} S − 1 ,得到 A = S D S − 1 A=SDS^{-1} A = S D S − 1 。反过来,对这两
个等式作相应的乘法即可恢复 A S = S D AS=SD A S = S D ,再逐列恢复特征向量方程。可逆性
必须在约去 S − 1 S S^{-1}S S − 1 S 之前说明;可逆矩阵的列也自动非零。
这条字典能防止次序错误。如果 S S S 的第一列对应 λ 1 \lambda_1 λ 1 、第二列对应
λ 2 \lambda_2 λ 2 ,则对角线必须按 λ 1 , λ 2 \lambda_1,\lambda_2 λ 1 , λ 2 排列。若用置换矩阵
Q Q Q 改变列的次序,就必须同时改变对角线的次序。
同样,若 V = [ v 1 ⋯ v p ] V=[v_1\ \cdots\ v_p] V = [ v 1 ⋯ v p ] 是 n × p n\times p n × p 矩阵,且 D = diag ( μ 1 , … , μ p ) D=\operatorname{diag}(\mu_1,\ldots,\mu_p) D = diag ( μ 1 , … , μ p ) ,则 A V = V D AV=VD A V = V D 等价于每一列都满足 A v j = μ j v j Av_j=\mu_jv_j A v j = μ j v j 。当 p < n p<n p < n 时,这仍能记录一部分特征向量的方程,却不是对角化:非方阵 V V V 没有双边逆,线性相关的列也不能组成完整坐标系。只有 p = n p=n p = n 且各列线性无关时,这个等式才给出对角化。
定理
可对角化的精确判准 实数域上的 n × n n\times n n × n 矩阵 A A A 可对角化,当且仅当它有 n n n 个线性无关
的实特征向量。等价地,这些特征向量可以组成 R n \mathbb{R}^n R n 的一组基底。
必要性来自对角化矩阵 S S S 的 n n n 个列:S S S 可逆,所以列向量线性无关,
特征向量字典又说明它们都是特征向量。充分性则是把 n n n 个线性无关的
特征向量作为列组成 S S S ;线性无关保证 S S S 可逆,字典随即给出对角化。
仅仅列出 n n n 个特征向量并不够,因为重复或线性相关的向量不能组成坐标
基底。
对角化例子
例题
由一个 2 乘 2 矩阵建立 S S S 和 D D D 令
A = [ 2 1 1 2 ] . A=\begin{bmatrix}2&1\\1&2\end{bmatrix}. A = [ 2 1 1 2 ] . 特征多项式为
det ( A − λ I ) = ( 2 − λ ) 2 − 1 = λ 2 − 4 λ + 3 = ( λ − 1 ) ( λ − 3 ) . \det(A-\lambda I)=(2-\lambda)^2-1
=\lambda^2-4\lambda+3=(\lambda-1)(\lambda-3). det ( A − λ I ) = ( 2 − λ ) 2 − 1 = λ 2 − 4 λ + 3 = ( λ − 1 ) ( λ − 3 ) . 对 λ = 3 \lambda=3 λ = 3 ,方程 ( A − 3 I ) x = 0 (A-3I)x=0 ( A − 3 I ) x = 0 给出 x 2 = x 1 x_2=x_1 x 2 = x 1 ,所以取
v 1 = [ 1 1 ] v_1=\begin{bmatrix}1\\1\end{bmatrix} v 1 = [ 1 1 ] 。对 λ = 1 \lambda=1 λ = 1 ,方程
( A − I ) x = 0 (A-I)x=0 ( A − I ) x = 0 给出 x 2 = − x 1 x_2=-x_1 x 2 = − x 1 ,所以取
v 2 = [ 1 − 1 ] v_2=\begin{bmatrix}1\\-1\end{bmatrix} v 2 = [ 1 − 1 ] 。按此顺序组成
S = [ 1 1 1 − 1 ] , D = [ 3 0 0 1 ] . S=\begin{bmatrix}1&1\\1&-1\end{bmatrix},\qquad
D=\begin{bmatrix}3&0\\0&1\end{bmatrix}. S = [ 1 1 1 − 1 ] , D = [ 3 0 0 1 ] . det ( S ) = − 2 \det(S)=-2 det ( S ) = − 2 ,并且
S − 1 = 1 2 [ 1 1 1 − 1 ] = 1 2 S . S^{-1}=\frac12\begin{bmatrix}1&1\\1&-1\end{bmatrix}=\frac12S. S − 1 = 2 1 [ 1 1 1 − 1 ] = 2 1 S . 两条列方程给出 A S = S D AS=SD A S = S D ,因而 S − 1 A S = D S^{-1}AS=D S − 1 A S = D 。再作坐标检查:取
x = [ 3 1 ] x=\begin{bmatrix}3\\1\end{bmatrix} x = [ 3 1 ] ,则
x = S y x=Sy x = S y ,其中 y = [ 2 1 ] y=\begin{bmatrix}2\\1\end{bmatrix} y = [ 2 1 ] ,所以
D y = [ 6 1 ] , S D y = [ 7 5 ] = A x . Dy=\begin{bmatrix}6\\1\end{bmatrix},\qquad
SDy=\begin{bmatrix}7\\5\end{bmatrix}=Ax. Dy = [ 6 1 ] , S Dy = [ 7 5 ] = A x . 在新坐标中作对角计算,和原矩阵直接作用的结果一致。
例题
一个可对角化的上三角矩阵 令
A = [ 1 1 1 0 2 2 0 0 3 ] . A=\begin{bmatrix}1&1&1\\0&2&2\\0&0&3\end{bmatrix}. A = 1 0 0 1 2 0 1 2 3 . 取以下特征向量
u 1 = [ 1 0 0 ] , u 2 = [ 1 1 0 ] , u 3 = [ 3 4 2 ] , u_1=\begin{bmatrix}1\\0\\0\end{bmatrix},\quad
u_2=\begin{bmatrix}1\\1\\0\end{bmatrix},\quad
u_3=\begin{bmatrix}3\\4\\2\end{bmatrix}, u 1 = 1 0 0 , u 2 = 1 1 0 , u 3 = 3 4 2 , 它们对应的特征值依次为 1 , 2 , 3 1,2,3 1 , 2 , 3 。令
U = [ u 1 u 2 u 3 ] = [ 1 1 3 0 1 4 0 0 2 ] , D = diag ( 1 , 2 , 3 ) . U=[u_1\ u_2\ u_3]=\begin{bmatrix}1&1&3\\0&1&4\\0&0&2\end{bmatrix},
\quad D=\operatorname{diag}(1,2,3). U = [ u 1 u 2 u 3 ] = 1 0 0 1 1 0 3 4 2 , D = diag ( 1 , 2 , 3 ) . det ( U ) = 2 \det(U)=2 det ( U ) = 2 ,所以 U U U 可逆,且
U − 1 = [ 1 − 1 1 / 2 0 1 − 2 0 0 1 / 2 ] . U^{-1}=\begin{bmatrix}1&-1&1/2\\0&1&-2\\0&0&1/2\end{bmatrix}. U − 1 = 1 0 0 − 1 1 0 1/2 − 2 1/2 . 直接相乘可得
A u 1 = [ 1 0 0 ] = 1 u 1 , A u 2 = [ 2 2 0 ] = 2 u 2 , A u 3 = [ 9 12 6 ] = 3 u 3 . Au_1=\begin{bmatrix}1\\0\\0\end{bmatrix}=1u_1,\quad
Au_2=\begin{bmatrix}2\\2\\0\end{bmatrix}=2u_2,\quad
Au_3=\begin{bmatrix}9\\12\\6\end{bmatrix} =3u_3. A u 1 = 1 0 0 = 1 u 1 , A u 2 = 2 2 0 = 2 u 2 , A u 3 = 9 12 6 = 3 u 3 . 所以 A U = U D AU=UD A U = U D ,并且 U − 1 A U = diag ( 1 , 2 , 3 ) U^{-1}AU=\operatorname{diag}(1,2,3) U − 1 A U = diag ( 1 , 2 , 3 ) 。真正保证可
对角化的是线性无关的特征向量,而不是上三角形状本身。
例题
重复特征值不一定给出完整特征基 考虑
J = [ 1 4 0 1 ] . J=\begin{bmatrix}1&4\\0&1\end{bmatrix}. J = [ 1 0 4 1 ] . det ( J − λ I ) = ( 1 − λ ) 2 \det(J-\lambda I)=(1-\lambda)^2 det ( J − λ I ) = ( 1 − λ ) 2 ,所以唯一特征值是 1 1 1 。而
( J − I ) x = [ 0 4 0 0 ] [ x 1 x 2 ] = 0 (J-I)x=\begin{bmatrix}0&4\\0&0\end{bmatrix}
\begin{bmatrix}x_1\\x_2\end{bmatrix}=0 ( J − I ) x = [ 0 0 4 0 ] [ x 1 x 2 ] = 0 迫使 x 2 = 0 x_2=0 x 2 = 0 ,x 1 x_1 x 1 任意。因此特征空间是
span { [ 1 0 ] } \operatorname{span}\{\begin{bmatrix}1\\0\end{bmatrix}\} span { [ 1 0 ] } ,维数只有一。
2 × 2 2\times2 2 × 2 矩阵需要两个线性无关的特征向量,所以 J J J 不可对角化。重复
特征值并不自动排除对角化,关键是对应特征空间是否提供足够多的方向。
置换与相似不变量
若 Q Q Q 是 n × n n\times n n × n 置换矩阵,且 A = S D S − 1 A=SDS^{-1} A = S D S − 1 ,则 S Q SQ SQ 仍可逆,并且
( S Q ) − 1 A ( S Q ) = Q − 1 S − 1 A S Q = Q − 1 D Q . (SQ)^{-1}A(SQ)=Q^{-1}S^{-1}ASQ=Q^{-1}DQ. ( SQ ) − 1 A ( SQ ) = Q − 1 S − 1 A SQ = Q − 1 D Q .
右侧仍是对角矩阵,只是对角元素按 Q Q Q 重排。因此对角化通常不唯一。
把任一特征向量列乘以非零标量也可以,只要相应的对角值和列次序仍然
同步。
定理
相似矩阵有相同的特征多项式 若 B = S − 1 A S B=S^{-1}AS B = S − 1 A S ,其中 S S S 可逆,则
det ( B − t I ) = det ( S − 1 ( A − t I ) S ) = det ( A − t I ) . \det(B-tI)=\det\bigl(S^{-1}(A-tI)S\bigr)=\det(A-tI). det ( B − t I ) = det ( S − 1 ( A − t I ) S ) = det ( A − t I ) . 所以相似矩阵有相同的特征多项式,也有相同的特征值(按代数重数计)。
这里使用 S − 1 ( A − t I ) S = S − 1 A S − t S − 1 I S S^{-1}(A-tI)S=S^{-1}AS-tS^{-1}IS S − 1 ( A − t I ) S = S − 1 A S − t S − 1 I S 、S − 1 I S = I S^{-1}IS=I S − 1 I S = I ,以及行列
式的乘法性质。相似并不表示两个矩阵的元素或特征向量相同;相应的特征
空间由坐标转换相互对应。
反例模式
相同特征资料不足以证明相似 令
N = [ 0 1 0 0 ] , 0 2 = [ 0 0 0 0 ] . N=\begin{bmatrix}0&1\\0&0\end{bmatrix},\qquad
0_2=\begin{bmatrix}0&0\\0&0\end{bmatrix}. N = [ 0 0 1 0 ] , 0 2 = [ 0 0 0 0 ] . 两者的特征多项式都是 t 2 t^2 t 2 ,但它们不相似。对任意可逆矩阵 S S S ,都有
S − 1 0 2 S = 0 2 S^{-1}0_2S=0_2 S − 1 0 2 S = 0 2 ,而 N ≠ 0 2 N\ne0_2 N = 0 2 。相同特征多项式是相似的必要条件,
不是充分条件。
幂、逆矩阵与转置
定理
由一次对角化得到幂、逆矩阵与转置 设 A = S D S − 1 A=SDS^{-1} A = S D S − 1 ,其中 S S S 可逆,且
D = diag ( λ 1 , … , λ n ) D=\operatorname{diag}(\lambda_1,\ldots,\lambda_n) D = diag ( λ 1 , … , λ n ) 。对每个正整数 m m m ,
A m = S D m S − 1 , D m = diag ( λ 1 m , … , λ n m ) . A^m=SD^mS^{-1},\qquad
D^m=\operatorname{diag}(\lambda_1^m,\ldots,\lambda_n^m). A m = S D m S − 1 , D m = diag ( λ 1 m , … , λ n m ) . 若 A A A 可逆,则所有 λ j \lambda_j λ j 非零,并且
A − 1 = S D − 1 S − 1 A^{-1}=SD^{-1}S^{-1} A − 1 = S D − 1 S − 1 。此外,A T A^T A T 也可对角化,并有相同特征值。
正整数公式从 m = 1 m=1 m = 1 开始;另行规定 A 0 = I n A^0=I_n A 0 = I n 、D 0 = I n D^0=I_n D 0 = I n 。负整数幂需要
可逆性:只有当每个对角元素非零时,D − 1 D^{-1} D − 1 和 A − 1 A^{-1} A − 1 才存在,才能
把公式延伸到 m < 0 m<0 m < 0 。
证明
幂与转置公式为什么成立 若已知 A m = S D m S − 1 A^m=SD^mS^{-1} A m = S D m S − 1 ,则
A m + 1 = ( S D m S − 1 ) ( S D S − 1 ) = S D m + 1 S − 1 , A^{m+1}=(SD^mS^{-1})(SDS^{-1})=SD^{m+1}S^{-1}, A m + 1 = ( S D m S − 1 ) ( S D S − 1 ) = S D m + 1 S − 1 , 其中约去了相邻的 S − 1 S S^{-1}S S − 1 S ,这给出归纳步骤。逆矩阵的因子次序反过来,
所以 ( S D S − 1 ) − 1 = S D − 1 S − 1 (SDS^{-1})^{-1}=SD^{-1}S^{-1} ( S D S − 1 ) − 1 = S D − 1 S − 1 。对转置,
A T = ( S − 1 ) T D S T . A^T=(S^{-1})^TDS^T. A T = ( S − 1 ) T D S T . 令 W = ( S − 1 ) T W=(S^{-1})^T W = ( S − 1 ) T ,则 W − 1 = S T W^{-1}=S^T W − 1 = S T ,从而
W − 1 A T W = D T = D W^{-1}A^TW=D^T=D W − 1 A T W = D T = D 。因此 W W W 的列是 A T A^T A T 的特征向量;原来 S S S 的列
未必是 A T A^T A T 的特征向量。
例题
通过对角化计算对称矩阵的所有正整数幂 令
B = [ 2 1 1 1 2 1 1 1 2 ] , S = [ 1 1 1 1 − 1 0 1 0 − 1 ] , S − 1 = 1 3 [ 1 1 1 1 − 2 1 1 1 − 2 ] , B=\begin{bmatrix}2&1&1\\1&2&1\\1&1&2\end{bmatrix},\quad
S=\begin{bmatrix}1&1&1\\1&-1&0\\1&0&-1\end{bmatrix},\quad
S^{-1}=\frac13\begin{bmatrix}1&1&1\\1&-2&1\\1&1&-2\end{bmatrix}, B = 2 1 1 1 2 1 1 1 2 , S = 1 1 1 1 − 1 0 1 0 − 1 , S − 1 = 3 1 1 1 1 1 − 2 1 1 1 − 2 , 并令 D = diag ( 4 , 1 , 1 ) D=\operatorname{diag}(4,1,1) D = diag ( 4 , 1 , 1 ) 。S S S 的三列是特征向量,所以
B = S D S − 1 B=SDS^{-1} B = S D S − 1 。对任意正整数 m m m ,
B m = 1 3 [ 4 m + 2 4 m − 1 4 m − 1 4 m − 1 4 m + 2 4 m − 1 4 m − 1 4 m − 1 4 m + 2 ] . B^m=\frac13\begin{bmatrix}
4^m+2&4^m-1&4^m-1\\
4^m-1&4^m+2&4^m-1\\
4^m-1&4^m-1&4^m+2
\end{bmatrix}. B m = 3 1 4 m + 2 4 m − 1 4 m − 1 4 m − 1 4 m + 2 4 m − 1 4 m − 1 4 m − 1 4 m + 2 . 当 m = 1 m=1 m = 1 时,对角元素为 2 2 2 、非对角元素为 1 1 1 ,恢复原矩阵;当 m = 2 m=2 m = 2
时,对角元素为 6 6 6 、非对角元素为 5 5 5 ,与直接相乘相同。之后计算别的
正整数幂时,不需要重新寻找特征向量。
复数扩展与共同特征基
在实数域上,
R = [ 1 − 1 1 1 ] R=\begin{bmatrix}1&-1\\1&1\end{bmatrix} R = [ 1 1 − 1 1 ]
的特征多项式是 ( 1 − λ ) 2 + 1 (1-\lambda)^2+1 ( 1 − λ ) 2 + 1 ,没有实特征值,所以不能在
R \mathbb{R} R 上对角化。在 C \mathbb{C} C 上可以取
U = [ i − i 1 1 ] , U − 1 = 1 2 [ − i 1 i 1 ] , D = diag ( 1 + i , 1 − i ) . U=\begin{bmatrix}i&-i\\1&1\end{bmatrix},\quad
U^{-1}=\frac12\begin{bmatrix}-i&1\\i&1\end{bmatrix},\quad
D=\operatorname{diag}(1+i,1-i). U = [ i 1 − i 1 ] , U − 1 = 2 1 [ − i i 1 1 ] , D = diag ( 1 + i , 1 − i ) .
此时 U − 1 R U = D U^{-1}RU=D U − 1 R U = D ,且 det ( U ) = 2 i ≠ 0 \det(U)=2i\ne0 det ( U ) = 2 i = 0 。这说明复数对角化是换了标量
域,并不与实数域上的结论矛盾。
定理
拥有共同特征基的和与积 若 A = S D S − 1 A=SDS^{-1} A = S D S − 1 、B = S E S − 1 B=SES^{-1} B = SE S − 1 使用同一个可逆矩阵 S S S ,且 D , E D,E D , E 都是
对角矩阵,则对任意标量 α , β \alpha,\beta α , β ,
α A + β B = S ( α D + β E ) S − 1 , A B = S D E S − 1 . \alpha A+\beta B=S(\alpha D+\beta E)S^{-1},\qquad
AB=SDE S^{-1}. α A + βB = S ( α D + βE ) S − 1 , A B = S D E S − 1 . 中间的矩阵仍为对角矩阵,所以这些矩阵由同一组基底同时对角化。
证明只需分配律和 S − 1 S = I S^{-1}S=I S − 1 S = I :
α S D S − 1 + β S E S − 1 = S ( α D + β E ) S − 1 \alpha SDS^{-1}+\beta SES^{-1}=S(\alpha D+\beta E)S^{-1} α S D S − 1 + βSE S − 1 = S ( α D + βE ) S − 1 ,而
( S D S − 1 ) ( S E S − 1 ) = S D E S − 1 (SDS^{-1})(SES^{-1})=SDE S^{-1} ( S D S − 1 ) ( SE S − 1 ) = S D E S − 1 。共同的 S S S 是必要条件;分别可对角化
并不能自动提供共同特征基。
重复特征值能说明什么
特征多项式记录代数重数,特征空间维数记录可用的线性无关特征向量数目,两者回答
不同问题。若重复特征值的特征空间维数不足,重新排列列向量也不能创造缺少的方向,
矩阵 J J J 正好显示了这种障碍。若特征空间有足够维数,重复的对角元素并不会造成
问题;对角矩阵本身也可以重复同一个特征值,同时仍有完整特征基。
相似关系保持特征多项式,是因为它只改变同一个线性作用的坐标表示。若
A v = λ v Av=\lambda v A v = λ v 且 B = S − 1 A S B=S^{-1}AS B = S − 1 A S ,则 B ( S − 1 v ) = λ ( S − 1 v ) B(S^{-1}v)=\lambda(S^{-1}v) B ( S − 1 v ) = λ ( S − 1 v ) ;由于 S − 1 S^{-1} S − 1
可逆,向量不会变成零。这说明特征值相同,但字面上的特征向量通常不同。也因此,
特征多项式只是必要的不变量,不能记录特征空间如何排列。
核对最后的证明依据
要判断给定数域上的可对角化性,先找出所有特征值,再把各特征空间的维数相加。
这个总数恰为 n n n ,当且仅当一个 n × n n\times n n × n 矩阵有特征基。把选出的有序基底
作为 S S S 的列,让 D D D 中的特征值使用同一次序,再验证 A S = S D AS=SD A S = S D 和 S S S 可逆。
这两项便能确认所写的对角化,不必重新计算一次 S − 1 S^{-1} S − 1 。
快速检查与练习
思考检查
设 S S S 可逆,且 S − 1 A S = D S^{-1}AS=D S − 1 A S = D 为对角矩阵。S S S 的各列与 A A A 的特征向量及 D D D 的对角元素有什么关系?
解答 · 答案 各列必须是 A A A 的线性无关特征向量,组成特征基;第 j j j 列对应 D D D 的第 j j j 个对角元素。
思考检查
一个可对角化的 3×3 矩阵,最多能有两个线性无关的特征向量吗?
解答 · 答案 不能。可对角化的 3 × 3 3\times3 3 × 3 矩阵最多的线性无关特征向量数目是三个,不可能只有两个。
思考检查
若 A A A 与对角矩阵 D D D 相似,它们有相同特征值吗?
解答 · 答案 有。相似矩阵有相同的特征多项式,因此有相同特征值(按代数重数计)。
思考检查
设 A A A 有特征向量 v 1 , v 2 , v 3 v_1,v_2,v_3 v 1 , v 2 , v 3 ,对应特征值为 2 , 5 , − 1 2,5,-1 2 , 5 , − 1 ,且三向量线性无关。3 × 3 3×3 3 × 3 矩阵 A A A 的对角化中出现什么对角矩阵?
解答 · 导引解答 若 S = [ v 1 v 2 v 3 ] S=[v_1\ v_2\ v_3] S = [ v 1 v 2 v 3 ] ,则
S − 1 A S = diag ( 2 , 5 , − 1 ) . S^{-1}AS=\operatorname{diag}(2,5,-1). S − 1 A S = diag ( 2 , 5 , − 1 ) .
思考检查
为什么 J = [ 1 4 0 1 ] J=\begin{bmatrix}1&4\\0&1\end{bmatrix} J = [ 1 0 4 1 ] 不可对角化?
解答 · 导引解答 唯一特征值是 1 1 1 ,而 ( J − I ) x = 0 (J-I)x=0 ( J − I ) x = 0 强制 x 2 = 0 x_2=0 x 2 = 0 ,所以特征空间只有一维;
2 × 2 2\times2 2 × 2 矩阵需要两个线性无关的特征向量。
思考检查
若 A = S D S − 1 A=SDS^{-1} A = S D S − 1 ,且 D = diag ( 2 , − 1 ) D=\operatorname{diag}(2,-1) D = diag ( 2 , − 1 ) ,则 A 3 A^3 A 3 是什么?
解答 · 导引解答 A 3 = S D 3 S − 1 = S diag ( 8 , − 1 ) S − 1 . A^3=SD^3S^{-1}=S\operatorname{diag}(8,-1)S^{-1}. A 3 = S D 3 S − 1 = S diag ( 8 , − 1 ) S − 1 .
相关笔记
若还不熟悉齐次方程组的观点,请先阅读
8.1 特征值、特征向量与特征空间 。
接着阅读
8.3 特征多项式与对角化判别 ,
学习判断特征向量是否足够的多项式工具。本节的基底语言也依赖
6.5 基底与维数 。