三组分层:A 组概念判断(快速自检)、B 组推理证明(考试大题)、C 组计算实操。全部附解答 (默认折叠),建议先独立完成再展开对照。所有题目只使用 3.1 、3.2 、3.4 、3.5 的概念与第一、二章先修知识。
本页 A、B、C 三组 61 题:A 组(20 题)、B 组(17 题)、C 组(24 题)。
组 定位 建议时长 卡壳时怎么办
A 组 (20 题)概念判断,口答即可 每题 1–2 分钟 直接看解答,回到对应小节的「即时自测」重读
B 组 (17 题)推理证明,需完整写出 每题 8–12 分钟 先看提示,再用 10 分钟;仍无思路才展开解答
C 组 (24 题)计算实操,注意验算 每题 5–8 分钟 看提示定位方法,算完一定代回验证
关于「提示」层
B5、B7、B8、B11、B15、C4、C13、C14、C17 这 9 题提供了折叠的「提示」——它只给从哪里下手 ,不给完整过程。先看提示、自己推完,再对照解答,学习效果远好于直接读解答。
关于难度分布
A 组覆盖全部「辨析点」:相抵分解不唯一 / 标准形唯一、相抵与相似的分野、单侧逆的存在判别与「方阵前提」、rank ( A ) = rank ( A H A ) \operatorname{rank}(A)=\operatorname{rank}(A^{\mathrm H}A) rank ( A ) = rank ( A H A ) 的四个等式、满秩分解的形状与读法、两种分解各自「不唯一的程度」、A X A = A AXA=A A X A = A 的解(A1–A12);QR 分解的存在 / 唯一 / 正线、列正交规范与酉矩阵的分野、QR 算法的相似性与收敛条件(A13–A16);Schur 引理的适用范围、实数域边界、谱映射与特征向量继承的单向性、最小多项式的四条事实、逆矩阵公式(A17–A20)。
B 组以相抵 ⟺ \iff ⟺ 同秩 、单侧逆的构造与判别 、满秩分解的两条主干证明 (rank C = r \operatorname{rank}C=r rank C = r 、定理 3.2.2 )与 A X A = A ⟺ C X B = I r AXA=A\iff CXB=I_r A X A = A ⟺ CXB = I r 为重心(B1–B10,B9 把 3.1 的结论接回第一章的正定矩阵);B11–B13 补 QR 的两条主干(对角酉阵刻画、迭代相似链与等模卡死),B14–B17 补 Schur / 谱映射 / 最小多项式的证明重做(B15 给出不依赖定理 3.5.1 的直接构造路线)。
C 组训练六件基本功:相抵分解(含分块矩阵法求 P − 1 P^{-1} P − 1 、Q − 1 Q^{-1} Q − 1 ) 、用相抵标准形解方程组 、单侧逆的构造与「全部右逆」的写法 、满秩分解的三种算法 (取列 / RREF / 从相抵分解读)、GS 手算 QR 与长方阵补齐 (C17、C18、C22)、H–C 与最小多项式驱动的三套计算 (C19–C21、C23、C24)。3.3 及以后素材到位后,本页按同样的三组结构续写。
A1 · 相抵分解的四个基本事实
判断对错并说明理由:
(a) 任何矩阵都有相抵分解;(b) 矩阵的相抵分解是唯一的;(c) 矩阵的相抵标准形是唯一的;(d) 两个同型矩阵相抵 ⟺ \iff ⟺ 它们的秩相等。
解答 · A1
(a) 对 (引理 3.1.1 ;零矩阵也成立:r = 0 r=0 r = 0 时 Γ = O \Gamma=O Γ = O ,取 P = Q = I P=Q=I P = Q = I )。(b) 错 ——分解式 A = P − 1 Γ Q A=P^{-1}\Gamma Q A = P − 1 Γ Q 中的 P P P 、Q Q Q 有无穷多种取法(对任意 ∗ * ∗ ,把 [ I r ∗ 0 I ] \begin{bmatrix}I_r&*\\0&I\end{bmatrix} [ I r 0 ∗ I ] 并进 P − 1 P^{-1} P − 1 、[ I r 0 ∗ I ] \begin{bmatrix}I_r&0\\*&I\end{bmatrix} [ I r ∗ 0 I ] 并进 Q Q Q 仍给出分解,见 3.1.1 )。(c) 对 ——Γ \Gamma Γ 只由秩 r r r 决定(秩在相抵下不变)。(d) 对 (本编补注 :同秩者都相抵于同一个 Γ \Gamma Γ ,再用相抵的传递性)。
要点 :「不唯一」说的是分解式 ,「唯一」说的是标准形 ——两句话各有对象,不能互换。
A2 · 相抵与相似
B = P − 1 A Q B=P^{-1}AQ B = P − 1 A Q 与 B = P − 1 A P B=P^{-1}AP B = P − 1 A P 分别刻画什么?它们之间是什么关系?何时两者重合?
解答 · A2
B = P − 1 A Q B=P^{-1}AQ B = P − 1 A Q (A , B A,B A , B 同型 )是相抵 :定理 2.3.5 的换基公式——同一个线性映射、出发端换基 Q Q Q 、到达端换基 P P P (P , Q P,Q P , Q 相互独立),或更一般地「A A A 经初等变换化为 B B B 」;
B = P − 1 A P B=P^{-1}AP B = P − 1 A P 是相似 (推论 2.3.4 ):同一个线性变换 (同一空间 V → V V\to V V → V )换一组基,两侧用同一个 P P P 。
关系:相似是相抵的特例 (取 Q = P Q=P Q = P ),故相似必相抵;反之不然(相抵允许 P , Q P,Q P , Q 独立,且不要求方阵——2 × 3 2\times3 2 × 3 矩阵也谈相抵,但只有同阶方阵才谈相似)。两者重合:当 A , B A,B A , B 同阶方阵且换基在出发端与到达端是同一件事时(Q = P Q=P Q = P )。
A3 · 单侧逆的存在性判断
判断下列矩阵是否有左逆、右逆(有则说明理由,没有也说明理由):
(a) A 1 = [ 1 0 0 0 1 0 ] A_1=\begin{bmatrix}1&0&0\\ 0&1&0\end{bmatrix} A 1 = [ 1 0 0 1 0 0 ] ;(b) A 2 = [ 1 0 0 1 0 0 ] A_2=\begin{bmatrix}1&0\\ 0&1\\ 0&0\end{bmatrix} A 2 = 1 0 0 0 1 0 ;(c) A 3 = [ 1 2 2 4 ] A_3=\begin{bmatrix}1&2\\ 2&4\end{bmatrix} A 3 = [ 1 2 2 4 ] ;(d) A 4 = [ 1 0 ] A_4=[1\ \ 0] A 4 = [ 1 0 ] 。
解答 · A3
由 定理 3.1.1 :右逆 ⟺ \iff ⟺ 行满秩;左逆 ⟺ \iff ⟺ 列满秩 。
(a) 2 × 3 2\times3 2 × 3 、秩 2 = m 2=m 2 = m :行满秩,有右逆 ;列不满秩(2 < 3 2<3 2 < 3 ),无左逆。
(b) 3 × 2 3\times2 3 × 2 、秩 2 = n 2=n 2 = n :列满秩,有左逆 ;行不满秩,无右逆。
(c) 秩 1 < min { 2 , 2 } 1<\min\{2,2\} 1 < min { 2 , 2 } :两者都没有 (顺带:方阵行满秩 ⟺ \iff ⟺ 可逆,故也可直接说「不可逆」)。
(d) 1 × 2 1\times2 1 × 2 、秩 1 = m 1=m 1 = m :行满秩,有右逆 (如 [ 1 0 ] \begin{bmatrix}1\\0\end{bmatrix} [ 1 0 ] );无左逆。
要点 :先看形状定方向(行满秩还是列满秩),再看秩是否达到该方向的上限。
A4 · 乘积为单位阵与「方阵前提」
(a) 设 A A A 是 n n n 阶方阵 ,若 A B = I n AB=I_n A B = I n ,B B B 是什么?唯一吗?(b) 设 A = [ 1 0 ] A=\begin{bmatrix}1&0\end{bmatrix} A = [ 1 0 ] ,A B = I 1 AB=I_1 A B = I 1 说明什么?B B B 唯一吗?(c) 把「乘积为单位阵的两矩阵必可逆」改成一个正确的命题。
解答 · A4
(a) A A A 可逆且 B = A − 1 B=A^{-1} B = A − 1 ,唯一 (这正是「方阵乘积为单位阵」的结论:A B = I ⇒ B = A − 1 ⇒ B A = I AB=I\Rightarrow B=A^{-1}\Rightarrow BA=I A B = I ⇒ B = A − 1 ⇒ B A = I )。
(b) 只说明 A A A 行满秩、B B B 是右逆 :设 B = ( b 1 , b 2 ) T B=(b_1,b_2)^{\mathrm T} B = ( b 1 , b 2 ) T ,则 A B = b 1 = 1 AB=b_1=1 A B = b 1 = 1 ,B = [ 1 t ] B=\begin{bmatrix}1\\ t\end{bmatrix} B = [ 1 t ] (t t t 任意)——不唯一 ;此时 B A = [ 1 0 t 0 ] ≠ I 2 BA=\begin{bmatrix}1&0\\ t&0\end{bmatrix}\neq I_2 B A = [ 1 t 0 0 ] = I 2 ,「逆矩阵」的语言完全失效。
(c) 正确说法:若两个方阵的乘积为单位阵,则这两个方阵都可逆 (且互为逆矩阵)。——「方阵」是关键词;长方阵情形只能谈单侧逆(定理 3.1.1 )。
A5 · 秩等式的四个判断
设 A ∈ C m × n A\in\mathbb C^{m\times n} A ∈ C m × n 。判断对错:
(a) rank ( A ) = rank ( A H A ) \operatorname{rank}(A)=\operatorname{rank}(A^{\mathrm H}A) rank ( A ) = rank ( A H A ) 恒成立;(b) A H A A^{\mathrm H}A A H A 可逆 ⟺ \iff ⟺ A A A 列满秩;(c) rank ( A H A ) = rank ( A A H ) \operatorname{rank}(A^{\mathrm H}A)=\operatorname{rank}(AA^{\mathrm H}) rank ( A H A ) = rank ( A A H ) 恒成立;(d) N ( A ) = N ( A H A ) N(A)=N(A^{\mathrm H}A) N ( A ) = N ( A H A ) 。
解答 · A5
四个全对 。
(a)(c)(d) 都是 命题 3.1.1 的内容:证明主线是 N ( A ) = N ( A H A ) N(A)=N(A^{\mathrm H}A) N ( A ) = N ( A H A ) (「⊆ \subseteq ⊆ 」两边左乘 A H A^{\mathrm H} A H ;「⊇ \supseteq ⊇ 」用内积正定性),再用亏加秩;对 A H A^{\mathrm H} A H 重复一遍得 rank ( A H ) = rank ( A A H ) \operatorname{rank}(A^{\mathrm H})=\operatorname{rank}(AA^{\mathrm H}) rank ( A H ) = rank ( A A H ) ,结合 rank A = rank A H \operatorname{rank}A=\operatorname{rank}A^{\mathrm H} rank A = rank A H 得四项等式链。
(b) 由 (a):A H A A^{\mathrm H}A A H A (n n n 阶)可逆 ⟺ \iff ⟺ rank ( A H A ) = n \operatorname{rank}(A^{\mathrm H}A)=n rank ( A H A ) = n ⟺ \iff ⟺ rank A = n \operatorname{rank}A=n rank A = n ⟺ \iff ⟺ A A A 列满秩。
注意 :结论依赖 A H A^{\mathrm H} A H (实 / 复内积空间);对一般数域 F F F 上没有共轭转置,此命题不适用。
A6 · 满秩分解的形状
设 A ∈ C 3 × 5 A\in\mathbb C^{3\times5} A ∈ C 3 × 5 ,rank A = 2 \operatorname{rank}A=2 rank A = 2 。若 A = B C A=BC A = BC 是满秩分解,B B B 、C C C 分别是几乘几、秩各是多少?B B B 的列与 C C C 的内容各有什么几何含义?
解答 · A6
B ∈ C 3 × 2 B\in\mathbb C^{3\times2} B ∈ C 3 × 2 ,rank B = 2 \operatorname{rank}B=2 rank B = 2 (列 满秩,2 2 2 列无关);C ∈ C 2 × 5 C\in\mathbb C^{2\times5} C ∈ C 2 × 5 ,rank C = 2 \operatorname{rank}C=2 rank C = 2 (行 满秩,2 2 2 行无关)。
几何含义:B B B 的 2 2 2 列是列空间 R ( A ) R(A) R ( A ) (dim = 2 \dim=2 dim = 2 )的一组基 ;C C C 的每一列 是 A A A 对应列在这组基下的坐标 (C = [ c 1 , ⋯ , c 5 ] C=[\boldsymbol{c}_1,\cdots,\boldsymbol{c}_5] C = [ c 1 , ⋯ , c 5 ] ,a i = B c i \boldsymbol{a}_i=B\boldsymbol{c}_i a i = B c i )。
形状记忆 :B B B 的行数 = m =m = m 、列数 = r =r = r ;C C C 的行数 = r =r = r 、列数 = n =n = n ——中间夹着 r r r 。
A7 · 满秩分解的「读法」
判断对错:
(a) R ( A ) = R ( B ) R(A)=R(B) R ( A ) = R ( B ) ;(b) N ( A ) = N ( C ) N(A)=N(C) N ( A ) = N ( C ) ;(c) A A A 的列之间的全部线性关系都记录在 C C C 里;(d) 满秩分解中的 B B B 必须由 A A A 的列向量组成。
解答 · A7
(a) 对 (定理 3.2.2 证明的第一步:由 C C C 行满秩、C C C 满射得 R ( A ) = R ( B C ) = R ( B ) R(A)=R(BC)=R(B) R ( A ) = R ( BC ) = R ( B ) )。
(b) 对 :C x = 0 ⇒ B C x = 0 ⇒ A x = 0 C\boldsymbol{x}=\boldsymbol{0}\Rightarrow BC\boldsymbol{x}=\boldsymbol{0}\Rightarrow A\boldsymbol{x}=\boldsymbol{0} C x = 0 ⇒ BC x = 0 ⇒ A x = 0 ;反之 A x = 0 ⇒ B ( C x ) = 0 ⇒ C x = 0 A\boldsymbol{x}=\boldsymbol{0}\Rightarrow B(C\boldsymbol{x})=\boldsymbol{0}\Rightarrow C\boldsymbol{x}=\boldsymbol{0} A x = 0 ⇒ B ( C x ) = 0 ⇒ C x = 0 (B B B 的列无关)。故两零空间相等——这正是 (c) 的严格含义:∑ k i a i = 0 ⟺ ∑ k i c i = 0 \sum k_i\boldsymbol{a}_i=\boldsymbol{0}\iff\sum k_i\boldsymbol{c}_i=\boldsymbol{0} ∑ k i a i = 0 ⟺ ∑ k i c i = 0 。
(c) 对 (即 (b) 的改述)。
(d) 错 :注 2 (注 2)的反例 A = I 2 A A=I_2A A = I 2 A :B = I 2 B=I_2 B = I 2 的任一列都不是 A A A 的列向量,但仍是合法的满秩分解。
A8 · 「不唯一」的两种程度
(a) 相抵分解不唯一到什么程度?(b) 满秩分解不唯一到什么程度?两句话的精确表述各是什么?
解答 · A8
(a) 相抵分解 A = P − 1 Γ Q A=P^{-1}\Gamma Q A = P − 1 Γ Q :P P P 、Q Q Q 各自独立 地有无穷多种取法(可被形如 [ I r ∗ 0 I ] \begin{bmatrix}I_r&*\\0&I\end{bmatrix} [ I r 0 ∗ I ] 、[ I r 0 ∗ I ] \begin{bmatrix}I_r&0\\*&I\end{bmatrix} [ I r ∗ 0 I ] 的可逆矩阵吸收);没有「相差一个矩阵」式的统一刻画,因为自由度分散在两侧。
(b) 满秩分解 A = B C A=BC A = BC :任意两种分解恰好相差一个可逆矩阵 D D D :
B 1 = B 2 D , C 1 = D − 1 C 2 B_1=B_2D,\qquad C_1=D^{-1}C_2 B 1 = B 2 D , C 1 = D − 1 C 2
(这就是 定理 3.2.2 。)
且 D D D 必可逆(是两组基之间的过渡矩阵)。反过来,任取可逆 D D D ,( B 1 D ) ( D − 1 C 1 ) = A (B_1D)(D^{-1}C_1)=A ( B 1 D ) ( D − 1 C 1 ) = A 仍给出满秩分解——一个分解即生成全部 (注 3 )。
要点 :满秩分解的不唯一被「一个 D D D 」精确刻画;相抵分解的不唯一更「松散」(两侧独立),但标准形唯一。
A9 · 矩阵方程 A X A = A AXA=A A X A = A
(a) 若 A A A 可逆,方程 A X A = A AXA=A A X A = A 的解是什么?唯一吗?(b) 一般情形(A ∈ C r m × n A\in\mathbb C^{m\times n}_r A ∈ C r m × n )课件是用什么办法把方程「降维」的?用一句话说清。
解答 · A9
(a) A A A 可逆时两边左乘 A − 1 A^{-1} A − 1 、右乘 A − 1 A^{-1} A − 1 得 X = A − 1 X=A^{-1} X = A − 1 ,唯一 (这也是 例 3.2.2 中「不能直接约掉 A A A 」这句话的反面:可逆时才能约)。
(b) 用满秩分解 A = B C A=BC A = BC :方程变成 B C X B C = B C BCXBC=BC BCXBC = BC ,其中 C X B CXB CXB 是 r r r 阶方阵——把 m × n m\times n m × n 的未知量 X X X 的问题「夹」成了 r r r 阶方阵上的条件 C X B = I r CXB=I_r CXB = I r ;再由 本编补注 ,二者等价 (课件写「充分条件」,其实是充要条件)。
A10 · 秩不等式的使用
(a) 判断:rank ( B C ) ≤ min { rank B , rank C } \operatorname{rank}(BC)\le\min\{\operatorname{rank}B,\operatorname{rank}C\} rank ( BC ) ≤ min { rank B , rank C } 。(b) 满秩分解中「rank C = r \operatorname{rank}C=r rank C = r 」的两步夹逼各是什么?(c) 为什么 rank ( B C ) ≤ rank C \operatorname{rank}(BC)\le\operatorname{rank}C rank ( BC ) ≤ rank C ?
解答 · A10
(a) 对 :R ( B C ) ⊆ R ( B ) R(BC)\subseteq R(B) R ( BC ) ⊆ R ( B ) 给出 rank ( B C ) ≤ rank B \operatorname{rank}(BC)\le\operatorname{rank}B rank ( BC ) ≤ rank B ;N ( C ) ⊆ N ( B C ) N(C)\subseteq N(BC) N ( C ) ⊆ N ( BC ) 加亏加秩给出 rank ( B C ) ≤ rank C \operatorname{rank}(BC)\le\operatorname{rank}C rank ( BC ) ≤ rank C 。
(b) 上界:C C C 只有 r r r 行,故 rank C ≤ r \operatorname{rank}C\le r rank C ≤ r ;下界:N ( C ) ⊆ N ( A ) N(C)\subseteq N(A) N ( C ) ⊆ N ( A ) (A = B C A=BC A = BC )加亏加秩得 rank C = n − dim N ( C ) ≥ n − dim N ( A ) = rank A = r \operatorname{rank}C=n-\dim N(C)\ge n-\dim N(A)=\operatorname{rank}A=r rank C = n − dim N ( C ) ≥ n − dim N ( A ) = rank A = r 。故 rank C = r \operatorname{rank}C=r rank C = r 。
(c) 若 C x = 0 C\boldsymbol{x}=\boldsymbol{0} C x = 0 则 B C x = 0 BC\boldsymbol{x}=\boldsymbol{0} BC x = 0 ,即 N ( C ) ⊆ N ( B C ) N(C)\subseteq N(BC) N ( C ) ⊆ N ( BC ) ;两个矩阵列数相同(都是 n n n ),亏加秩给出 rank ( B C ) = n − dim N ( B C ) ≤ n − dim N ( C ) = rank C \operatorname{rank}(BC)=n-\dim N(BC)\le n-\dim N(C)=\operatorname{rank}C rank ( BC ) = n − dim N ( BC ) ≤ n − dim N ( C ) = rank C 。
A11 · 相抵标准形的「几何」翻译
用基与线性映射的语言复述 引理 3.1.1 :Γ \Gamma Γ 的存在性等价于「总可以选到两组基,使映射的矩阵长成 Γ \Gamma Γ 的样子」。请说清这两组基是怎么选的。
解答 · A11
设 T T T 是「以 A A A 为矩阵」的映射(取定标准基后由 定理 2.3.1 唯一确定)。做法(引理 3.1.1 的本编证明):
出发端 :取零空间 N ( T ) N(T) N ( T ) 的一组基(n − r n-r n − r 个),用基扩充定理补成 F n F^n F n 的一组基 ξ 1 , ⋯ , ξ n − r , ξ n − r + 1 , ⋯ , ξ n \xi_1,\cdots,\xi_{n-r},\xi_{n-r+1},\cdots,\xi_n ξ 1 , ⋯ , ξ n − r , ξ n − r + 1 , ⋯ , ξ n ;
到达端 :证明后 r r r 个基向量的像 ξ n − r + 1 , ⋯ , ξ n \xi_{n-r+1},\cdots,\xi_n ξ n − r + 1 , ⋯ , ξ n 的像是 R ( T ) R(T) R ( T ) 的基(无关性由 ξ 1 , ⋯ , ξ n − r \xi_1,\cdots,\xi_{n-r} ξ 1 , ⋯ , ξ n − r 是零空间的基得出),把它们排在前面,扩充成 F m F^m F m 的基 η 1 , ⋯ , η m \eta_1,\cdots,\eta_m η 1 , ⋯ , η m ;
读出矩阵 :在这两组基下 T ( ξ 1 ) = ⋯ = T ( ξ n − r ) = 0 T(\xi_1)=\cdots=T(\xi_{n-r})=\boldsymbol{0} T ( ξ 1 ) = ⋯ = T ( ξ n − r ) = 0 、T ( ξ n − r + i ) = η i T(\xi_{n-r+i})=\eta_i T ( ξ n − r + i ) = η i ,故矩阵恰为 Γ \Gamma Γ ;换基公式(定理 2.3.5 )给出 A = P − 1 Γ Q A=P^{-1}\Gamma Q A = P − 1 Γ Q 。
一句话 :零空间的基排前面、值域的基接后面 ——Γ \Gamma Γ 就是「把 A A A 的活与死的部分分开摆好」的样子。
A12 · 纠错四句
改正下列说法:
(a)「相抵分解是唯一的」;(b)「行满秩矩阵有左逆」;(c)「A H A A^{\mathrm H}A A H A 与 A A A 的秩相等,因为它们是同一个二次型的矩阵」;(d)「满秩分解的 B B B 、C C C 可以随意取,只要 B C = A BC=A BC = A 」。
解答 · A12
(a) 应说:「相抵分解式 不唯一,相抵标准形 唯一」——见 A1。
(b) 应说:「行满秩矩阵有右 逆;列 满秩矩阵有左 逆」(定理 3.1.1 )。方向由 A B = I AB=I A B = I (右)与 B A = I BA=I B A = I (左)决定:行满秩 A A A 是「宽」的,可以右乘一个「高」的 B B B 凑出 I I I 。
(c) 二者秩相等对 ,但理由不对:理由不是「同一个二次型」(A H A A^{\mathrm H}A A H A 与 A A A 甚至不同型),而是 命题 3.1.1 的同解方程组证明(N ( A ) = N ( A H A ) N(A)=N(A^{\mathrm H}A) N ( A ) = N ( A H A ) ,再用亏加秩)。
(d) 不够:还要 B B B 列满秩 、C C C 行满秩 (定理 3.2.1 );否则如 [ 1 2 3 4 ] = [ 1 2 ] [ 1 2 ] \begin{bmatrix}1&2\\3&4\end{bmatrix}=\begin{bmatrix}1\\2\end{bmatrix}[1\ \ 2] [ 1 3 2 4 ] = [ 1 2 ] [ 1 2 ] 之类根本不成立,或退化成「秩不足」的分解(B B B 的列不满秩 ⇒ \Rightarrow ⇒ B B B 的列不是 R ( A ) R(A) R ( A ) 的基)。
A13 · QR 分解的四个基本事实
判断对错并说明理由:
(a) 任何满秩复方阵都有唯一的 QR 分解;(b) 不满秩的方阵一定没有 QR 分解;(c) 列满秩的长方阵 A ∈ C n m × n A\in\mathbb C^{m\times n}_n A ∈ C n m × n (n ≤ m n\le m n ≤ m )可分解为 A = U [ R 1 O ] A=U\begin{bmatrix}R_1\\O\end{bmatrix} A = U [ R 1 O ] ;(d) n n n 阶实满秩方阵在不要求「正线」时恰有 2 n 2^n 2 n 种正交三角分解。
解答 · A13
(a) 对 (定理 3.4.1 :满秩保证存在;「正线」约定保证唯一)。(b) 错 ——例 3.4.2 的 A = [ 1 1 0 0 ] A=\begin{bmatrix}1&1\\0&0\end{bmatrix} A = [ 1 0 1 0 ] 有 A = I 2 [ 1 1 0 0 ] A=I_2\begin{bmatrix}1&1\\0&0\end{bmatrix} A = I 2 [ 1 0 1 0 ] ,丢掉的只是唯一性 (R R R 不可逆)。(c) 对 (推论 3.4.1 :GS + + + 基扩充)。(d) 对 ——任意两种分解相差对角正交阵 D = diag ( ± 1 , ⋯ , ± 1 ) D=\operatorname{diag}(\pm1,\cdots,\pm1) D = diag ( ± 1 , ⋯ , ± 1 ) (2 n 2^n 2 n 个),见 3.4.1 的本编补注。
要点 :存在性的保底比想象中宽(不满秩、长方阵都常有分解),唯一性的门槛比想象中高(要满秩加 正线)。
A14 · 数清一个矩阵的全部 QR 分解
写出 A = [ 0 1 1 0 ] A=\begin{bmatrix}0&1\\ 1&0\end{bmatrix} A = [ 0 1 1 0 ] 的全部 正交三角分解(Q Q Q 正交、R R R 可逆上三角,不要求正线),指出其中正线的版本,并说明为什么没有第五种。
解答 · A14
A A A 本身是正交矩阵(置换阵,且 A 2 = I A^2=I A 2 = I ),两列 e 2 , e 1 \boldsymbol{e}_2,\boldsymbol{e}_1 e 2 , e 1 已标准正交。GS 直接给出正线版本 Q = A Q=A Q = A 、R = I R=I R = I (第一列 r 11 = ∥ e 2 ∥ = 1 r_{11}=\|\boldsymbol{e}_2\|=1 r 11 = ∥ e 2 ∥ = 1 ,r 12 = ( e 1 , e 2 ) = 0 r_{12}=(\boldsymbol{e}_1,\boldsymbol{e}_2)=0 r 12 = ( e 1 , e 2 ) = 0 ,r 22 = 1 r_{22}=1 r 22 = 1 )。由对角正交阵生成全部 2 2 = 4 2^2=4 2 2 = 4 种,形如 ( A ⋅ D , D ) \big(A\cdot D,\ D\big) ( A ⋅ D , D ) (R = I R=I R = I 被 D D D 左乘后仍是上三角):
( A , I ) , ( A diag ( − 1 , 1 ) , diag ( − 1 , 1 ) ) , ( A diag ( 1 , − 1 ) , diag ( 1 , − 1 ) ) , ( A diag ( − 1 , − 1 ) , diag ( − 1 , − 1 ) ) . (A,\ I),\qquad\big(A\operatorname{diag}(-1,1),\ \operatorname{diag}(-1,1)\big),\qquad\big(A\operatorname{diag}(1,-1),\ \operatorname{diag}(1,-1)\big),\qquad\big(A\operatorname{diag}(-1,-1),\ \operatorname{diag}(-1,-1)\big). ( A , I ) , ( A diag ( − 1 , 1 ) , diag ( − 1 , 1 ) ) , ( A diag ( 1 , − 1 ) , diag ( 1 , − 1 ) ) , ( A diag ( − 1 , − 1 ) , diag ( − 1 , − 1 ) ) .
(验算其一:A diag ( − 1 , 1 ) = [ 0 1 − 1 0 ] A\operatorname{diag}(-1,1)=\begin{bmatrix}0&1\\-1&0\end{bmatrix} A diag ( − 1 , 1 ) = [ 0 − 1 1 0 ] 正交 ✓,diag ( − 1 , 1 ) \operatorname{diag}(-1,1) diag ( − 1 , 1 ) 上三角 ✓,乘积 = [ 0 1 1 0 ] = A =\begin{bmatrix}0&1\\1&0\end{bmatrix}=A = [ 0 1 1 0 ] = A ✓。)正线版本 :只有 ( A , I ) (A,I) ( A , I ) (对角全 1 > 0 1>0 1 > 0 )。没有第五种 :若 A = Q ′ R ′ A=Q'R' A = Q ′ R ′ 是任一分解,则 W = R ′ R − 1 = Q ′ T Q W=R'R^{-1}=Q'^{\mathrm T}Q W = R ′ R − 1 = Q ′ T Q 既酉又上三角,故为对角阵且 ∣ w i i ∣ = 1 |w_{ii}|=1 ∣ w ii ∣ = 1 (3.4.1 的归纳),实数情形 w i = ± 1 w_i=\pm1 w i = ± 1 ——恰好穷尽。∎
A15 · 列正交规范矩阵辨析
设 Q ∈ C m × n Q\in\mathbb C^{m\times n} Q ∈ C m × n (m > n m>n m > n )满足 Q H Q = I n Q^{\mathrm H}Q=I_n Q H Q = I n 。判断:
(a) Q Q Q 是酉矩阵;(b) 对一切 x ∈ C n \boldsymbol{x}\in\mathbb C^n x ∈ C n 有 ∥ Q x ∥ = ∥ x ∥ \|Q\boldsymbol{x}\|=\|\boldsymbol{x}\| ∥ Q x ∥ = ∥ x ∥ ;(c) Q Q Q 列满秩,Q H Q^{\mathrm H} Q H 是它的一个左逆;(d) Q Q H = I m QQ^{\mathrm H}=I_m Q Q H = I m 。
解答 · A15
(a) 错 ——酉矩阵必须是方阵(正交矩阵与酉矩阵 2.5.1 );Q Q Q 只是列正交规范矩阵 (定义 3.4.2 )。(b) 对 :∥ Q x ∥ 2 = x H Q H Q x = ∥ x ∥ 2 \|Q\boldsymbol{x}\|^2=\boldsymbol{x}^{\mathrm H}Q^{\mathrm H}Q\boldsymbol{x}=\|\boldsymbol{x}\|^2 ∥ Q x ∥ 2 = x H Q H Q x = ∥ x ∥ 2 。(c) 对 :Q H Q = I n Q^{\mathrm H}Q=I_n Q H Q = I n 本身就是「列满秩 + + + 给出左逆」(定理 3.1.1 )。(d) 错 :Q Q H QQ^{\mathrm H} Q Q H 是把 C m \mathbb C^m C m 先投影到 n n n 维再放回的算子,与 I m I_m I m 差一个 m − n m-n m − n 维的正交补——例如 Q = [ 1 0 0 ] Q=\begin{bmatrix}1\\0\\0\end{bmatrix} Q = 1 0 0 时 Q Q H = diag ( 1 , 0 , 0 ) QQ^{\mathrm H}=\operatorname{diag}(1,0,0) Q Q H = diag ( 1 , 0 , 0 ) 。
要点 :Q H Q = I n Q^{\mathrm H}Q=I_n Q H Q = I n 与 Q Q H = I m QQ^{\mathrm H}=I_m Q Q H = I m 只在 m = n m=n m = n (方阵)时可能同时成立——形状就是命题的一部分 。
A16 · QR 算法四问
判断对错:
(a) 迭代的每个 A k A_k A k 都与 A A A 相似;(b)「各特征值模值互不相同」是 QR 算法收敛的必要条件;(c) 模互异时,收敛所得上三角阵的对角元按模从大到小排定;(d) 特征值等模时 QR 迭代一定不收敛。
解答 · A16
(a) 对 :A k + 1 = Q k H A k Q k A_{k+1}=Q_k^{\mathrm H}A_kQ_k A k + 1 = Q k H A k Q k ,链起来 A k = ( Q 1 ⋯ Q k − 1 ) H A ( Q 1 ⋯ Q k − 1 ) A_k=(Q_1\cdots Q_{k-1})^{\mathrm H}A(Q_1\cdots Q_{k-1}) A k = ( Q 1 ⋯ Q k − 1 ) H A ( Q 1 ⋯ Q k − 1 ) (QR 算法 )。(b) 错 ——是充分 条件,不是必要(例如 A A A 已经是对角元全正的上三角阵时迭代一步不动、天然「收敛」,无论对角元模是否互异)。(c) 对 (正文的本编注;例 3.4.4 中 4.7321 > 3 > 1.2679 4.7321>3>1.2679 4.7321 > 3 > 1.2679 即此序)。(d) 错 :等模时可能 卡死(旋转阵:A = A ⋅ I ⇒ A 2 = A A=A\cdot I\Rightarrow A_2=A A = A ⋅ I ⇒ A 2 = A ),但「一定不收敛」过强——见 (b) 的三角阵例子,特征值 λ 1 = λ 2 \lambda_1=\lambda_2 λ 1 = λ 2 (等模)迭代照样稳定。
要点 :背结论要背成「模互异 ⇒ \Rightarrow ⇒ 收敛」,别擅自加强成 ⟺ \Longleftrightarrow ⟺ 。
A17 · Schur 引理辨析
判断对错:
(a) 任意复方阵都酉相似于上三角阵;(b) 任意实方阵都正交相似于上三角阵;(c) 特征值全实的实方阵正交相似于上三角阵;(d) 任意复方阵都相似于对角阵。
解答 · A17
(a) 对 (定理 3.5.2 )。(b) 错 ——旋转阵 [ 0 − 1 1 0 ] \begin{bmatrix}0&-1\\1&0\end{bmatrix} [ 0 1 − 1 0 ] 特征值 ± i \pm\mathrm i ± i :正交相似保持实性,实上三角阵特征值全实,矛盾(3.5.1 的反例分析)。(c) 对 (定理 3.5.3 ,前提正是「特征值均为实数」)。(d) 错 (注 1:三角化人人有份,对角化不是——经典三矩阵:同一特征多项式,三种几何重数 2.4.2 )。
要点 :(a) 与 (d) 的差别就是「对角化」与「三角化」的差别;(b) 与 (c) 的差别就是复结论进入实数域时要补的门票。
A18 · 谱映射与特征向量继承
设 A A A 的特征值为 1 , − 1 1,-1 1 , − 1 (各一重,A A A 为二阶)。判断:
(a) A 2 + I A^2+I A 2 + I 的特征值为 2 , 2 2,2 2 , 2 ;(b) A A A 的特征向量都是 A 2 + I A^2+I A 2 + I 的特征向量;(c) A 2 + I A^2+I A 2 + I 的特征向量都是 A A A 的特征向量;(d) A 2026 A^{2026} A 2026 的特征值为 1 , 1 1,1 1 , 1 。
解答 · A18
(a) 对 (推论 3.5.1 :φ ( λ ) = λ 2 + 1 \varphi(\lambda)=\lambda^2+1 φ ( λ ) = λ 2 + 1 ,φ ( ± 1 ) = 2 \varphi(\pm1)=2 φ ( ± 1 ) = 2 )。(b) 对 (注 4 继承:A α = λ α ⇒ φ ( A ) α = φ ( λ ) α A\boldsymbol{\alpha}=\lambda\boldsymbol{\alpha}\Rightarrow\varphi(A)\boldsymbol{\alpha}=\varphi(\lambda)\boldsymbol{\alpha} A α = λ α ⇒ φ ( A ) α = φ ( λ ) α )。(c) 错 :A = diag ( 1 , − 1 ) A=\operatorname{diag}(1,-1) A = diag ( 1 , − 1 ) 时 A 2 + I = 2 I A^2+I=2I A 2 + I = 2 I ,一切非零向量都是它的特征向量,而 e 1 + e 2 \boldsymbol{e}_1+\boldsymbol{e}_2 e 1 + e 2 不是 A A A 的特征向量——φ \varphi φ 把两个特征值黏在一起了。(d) 对 :φ ( λ ) = λ 2026 \varphi(\lambda)=\lambda^{2026} φ ( λ ) = λ 2026 ,φ ( ± 1 ) = 1 \varphi(\pm1)=1 φ ( ± 1 ) = 1 (偶次幂把 ± 1 \pm1 ± 1 黏成 1 1 1 ,与 (c) 同源)。
A19 · 最小多项式四问
判断对错:
(a) 任何复方阵的最小多项式存在且唯一;(b) 最小多项式恒等于特征多项式;(c) m A ( λ ) m_A(\lambda) m A ( λ ) 整除 f A ( λ ) f_A(\lambda) f A ( λ ) ;(d) m A m_A m A 与 f A f_A f A 的根相同(不计重数)。
解答 · A19
(a) 对 (存在:f A f_A f A 本身零化(H–C),次数取最小;唯一:互相整除 + + + 同次首一,定理 3.5.5 (1))。(b) 错 :A = I 2 A=I_2 A = I 2 时 f A = ( λ − 1 ) 2 f_A=(\lambda-1)^2 f A = ( λ − 1 ) 2 而 m A = λ − 1 m_A=\lambda-1 m A = λ − 1 ;一般 m A = f A m_A=f_A m A = f A 当且仅当每个特征值的「最大块」都顶到代数重数(例 3.5.B 的情形)。(c) 对 (f A f_A f A 是零化多项式,m A m_A m A 整除一切零化多项式)。(d) 对 (定理 3.5.5 (2):正向由整除,反向由特征向量继承)。
A20 · 逆矩阵公式辨析
设 f A ( λ ) = λ n + a n − 1 λ n − 1 + ⋯ + a 1 λ + a 0 f_A(\lambda)=\lambda^n+a_{n-1}\lambda^{n-1}+\cdots+a_1\lambda+a_0 f A ( λ ) = λ n + a n − 1 λ n − 1 + ⋯ + a 1 λ + a 0 。判断:
(a) a 0 = 0 ⟺ A a_0=0\iff A a 0 = 0 ⟺ A 不可逆;(b) A A A 可逆时,A − 1 A^{-1} A − 1 可以写成 A A A 的多项式;(c) 推论 3.5.2 的公式是 A − 1 = 1 a 0 ( A n − 1 + a n − 1 A n − 2 + ⋯ + a 1 I ) A^{-1}=\frac{1}{a_0}(A^{n-1}+a_{n-1}A^{n-2}+\cdots+a_1I) A − 1 = a 0 1 ( A n − 1 + a n − 1 A n − 2 + ⋯ + a 1 I ) ;(d) 例 3.5.3 课件把答案 ( 1 , 3 ) (1,3) ( 1 , 3 ) 元印成 0 0 0 ,与它前面的特征多项式符号错无关,是独立的抄写错误。
解答 · A20
(a) 对 :a 0 = f A ( 0 ) = det ( − A ) = ( − 1 ) n det A a_0=f_A(0)=\det(-A)=(-1)^n\det A a 0 = f A ( 0 ) = det ( − A ) = ( − 1 ) n det A (推论 3.5.2 证明)。(b) 对 ——推论 3.5.2 的公式就是一个这样的多项式(这补上了谱映射讨论里「1 / λ 1/\lambda 1/ λ 不是多项式」的口子)。(c) 错 ——差一个负号:正确为 − 1 a 0 ( ⋯ ) -\frac{1}{a_0}(\cdots) − a 0 1 ( ⋯ ) (由 A ( ⋯ ) = − a 0 I A(\cdots)=-a_0I A ( ⋯ ) = − a 0 I 移项而来)。(d) 错 ——三处错误同源连锁:特征多项式符号错(( λ − 1 ) 3 (\lambda-1)^3 ( λ − 1 ) 3 印成 ( λ + 1 ) 3 (\lambda+1)^3 ( λ + 1 ) 3 )⇒ \Rightarrow ⇒ 公式行系数全错 ⇒ \Rightarrow ⇒ 答案矩阵错;即使按印错的公式行计算,答案也不是印出的那个矩阵(见 E3.5 的逐条核算)。
B1 · 相抵 ⟺ \iff ⟺ 同秩
证明:两个 m × n m\times n m × n 矩阵 A , B A,B A , B 相抵(存在可逆 P , Q P,Q P , Q 使 B = P − 1 A Q B=P^{-1}AQ B = P − 1 A Q )当且仅当 rank A = rank B \operatorname{rank}A=\operatorname{rank}B rank A = rank B 。
解答 · B1
必要性(⟹ \Longrightarrow ⟹ ) :初等变换不改变秩(本编提示:逐个验证三类行 / 列初等变换,或直接用「相抵 ⇒ \Rightarrow ⇒ 秩相等」的通用论证——P , Q P,Q P , Q 可逆给出 rank B = rank ( P − 1 A Q ) = rank A \operatorname{rank}B=\operatorname{rank}(P^{-1}AQ)=\operatorname{rank}A rank B = rank ( P − 1 A Q ) = rank A ,因为乘可逆矩阵不改变秩:rank ( P − 1 A Q ) = rank ( A Q ) = rank A \operatorname{rank}(P^{-1}AQ)=\operatorname{rank}(AQ)=\operatorname{rank}A rank ( P − 1 A Q ) = rank ( A Q ) = rank A ,每一步都用了「乘可逆矩阵不改变秩」。这里 rank ( X A ) = rank A \operatorname{rank}(XA)=\operatorname{rank}A rank ( X A ) = rank A 对可逆 X X X 成立,理由:R ( X A ) = X ( R ( A ) ) R(XA)=X(R(A)) R ( X A ) = X ( R ( A )) 且 X X X 是同构(命题 2.3.1 (5)),维数不变)。
充分性(⟸ \Longleftarrow ⟸ ) :设 rank A = rank B = r \operatorname{rank}A=\operatorname{rank}B=r rank A = rank B = r 。由 引理 3.1.1 ,A = P 1 − 1 Γ P 2 A=P_1^{-1}\Gamma P_2 A = P 1 − 1 Γ P 2 、B = Q 1 − 1 Γ Q 2 B=Q_1^{-1}\Gamma Q_2 B = Q 1 − 1 Γ Q 2 (同一个 Γ \Gamma Γ ,因为秩都是 r r r )。于是
B = Q 1 − 1 Γ Q 2 = Q 1 − 1 ( P 1 A P 2 − 1 ) Q 2 = ( Q 1 − 1 P 1 ) A ( P 2 − 1 Q 2 ) , B=Q_1^{-1}\Gamma Q_2=Q_1^{-1}(P_1AP_2^{-1})Q_2=(Q_1^{-1}P_1)\,A\,(P_2^{-1}Q_2), B = Q 1 − 1 Γ Q 2 = Q 1 − 1 ( P 1 A P 2 − 1 ) Q 2 = ( Q 1 − 1 P 1 ) A ( P 2 − 1 Q 2 ) ,
其中 Q 1 − 1 P 1 Q_1^{-1}P_1 Q 1 − 1 P 1 、P 2 − 1 Q 2 P_2^{-1}Q_2 P 2 − 1 Q 2 都可逆,故 A ≅ B A\cong B A ≅ B 。∎
B2 · 右逆的必要性(并说明左逆如何对称)
证明:若 A B = I m AB=I_m A B = I m ,则 A A A 行满秩。
解答 · B2
由 定理 2.3.6 ,rank ( A B ) = dim R ( A B ) \operatorname{rank}(AB)=\dim R(AB) rank ( A B ) = dim R ( A B ) ;而 R ( A B ) ⊆ R ( A ) R(AB)\subseteq R(A) R ( A B ) ⊆ R ( A ) ,故
rank ( A B ) ≤ rank A . \operatorname{rank}(AB)\le\operatorname{rank}A . rank ( A B ) ≤ rank A .
取 B = I m B=I_m B = I m 的情形即 rank ( A B ) = m \operatorname{rank}(AB)=m rank ( A B ) = m (A B AB A B 本身就是 I m I_m I m ,秩为 m m m ),于是 m ≤ rank A m\le\operatorname{rank}A m ≤ rank A ;另一方面 A A A 是 m m m 行矩阵,rank A ≤ m \operatorname{rank}A\le m rank A ≤ m 。两边夹得 rank A = m \operatorname{rank}A=m rank A = m ,即 A A A 行满秩。∎
左逆对称处理 :设 B A = I n BA=I_n B A = I n ,两边共轭转置得 A H B H = I n A^{\mathrm H}B^{\mathrm H}=I_n A H B H = I n ,即 A H A^{\mathrm H} A H 有右逆;刚刚证的结论给出 rank ( A H ) = n \operatorname{rank}(A^{\mathrm H})=n rank ( A H ) = n ;又 rank A = rank A H \operatorname{rank}A=\operatorname{rank}A^{\mathrm H} rank A = rank A H ,故 A A A 列满秩。(也可直接对 B A BA B A 用 rank ( B A ) ≤ rank A \operatorname{rank}(BA)\le\operatorname{rank}A rank ( B A ) ≤ rank A 。)
B3 · 列满秩时的左逆公式
证明:若 A ∈ C m × n A\in\mathbb C^{m\times n} A ∈ C m × n 列满秩,则 B = ( A H A ) − 1 A H B=(A^{\mathrm H}A)^{-1}A^{\mathrm H} B = ( A H A ) − 1 A H 是 A A A 的左逆。
解答 · B3
可逆性 :由 命题 3.1.1 ,rank ( A H A ) = rank A = n \operatorname{rank}(A^{\mathrm H}A)=\operatorname{rank}A=n rank ( A H A ) = rank A = n ,故 n n n 阶方阵 A H A A^{\mathrm H}A A H A 满秩、可逆。
验算 :
B A = ( A H A ) − 1 A H A = I n ✓ BA=(A^{\mathrm H}A)^{-1}A^{\mathrm H}A=I_n\ \checkmark B A = ( A H A ) − 1 A H A = I n ✓
要点 :这与 定理 3.1.1 「法 2」的右逆公式 A H ( A A H ) − 1 A^{\mathrm H}(AA^{\mathrm H})^{-1} A H ( A A H ) − 1 完全对称——把 A A A 换成 A H A^{\mathrm H} A H 即得。∎
B4 · 满秩分解中 rank C = r \operatorname{rank}C=r rank C = r
设 A = B C A=BC A = BC 是满秩分解(A ∈ C r m × n A\in\mathbb C^{m\times n}_r A ∈ C r m × n )。证明 rank C = r \operatorname{rank}C=r rank C = r (写出两步夹逼的完整理由)。
解答 · B4
上界 :C C C 是 r × n r\times n r × n 矩阵,其秩不超过行数:
rank C ≤ r . \operatorname{rank}C\le r . rank C ≤ r .
下界 :设 C x = 0 C\boldsymbol{x}=\boldsymbol{0} C x = 0 ,两端左乘 B B B 得 A x = B C x = 0 A\boldsymbol{x}=BC\boldsymbol{x}=\boldsymbol{0} A x = BC x = 0 ,故 N ( C ) ⊆ N ( A ) N(C)\subseteq N(A) N ( C ) ⊆ N ( A ) 。由亏加秩定理(定理 2.2.5 ,两矩阵列数同为 n n n ):
rank C = n − dim N ( C ) ≥ n − dim N ( A ) = rank A = r . \operatorname{rank}C=n-\dim N(C)\ \ge\ n-\dim N(A)=\operatorname{rank}A=r . rank C = n − dim N ( C ) ≥ n − dim N ( A ) = rank A = r .
结论 :r ≤ rank C ≤ r r\le\operatorname{rank}C\le r r ≤ rank C ≤ r ,故 rank C = r \operatorname{rank}C=r rank C = r 。∎
(B B B 的满秩是构造的一部分:B B B 的列取成 R ( A ) R(A) R ( A ) 的基,故 rank B = r \operatorname{rank}B=r rank B = r 显然。)
B5 · 满秩分解的唯一性(定理 3.2.2)
设 A = B 1 C 1 = B 2 C 2 A=B_1C_1=B_2C_2 A = B 1 C 1 = B 2 C 2 是 A ∈ C r m × n A\in\mathbb C^{m\times n}_r A ∈ C r m × n 的两种满秩分解。证明:存在可逆矩阵 D ∈ C r × r D\in\mathbb C^{r\times r} D ∈ C r × r 使 B 1 = B 2 D B_1=B_2D B 1 = B 2 D 、C 1 = D − 1 C 2 C_1=D^{-1}C_2 C 1 = D − 1 C 2 。
提示 · B5
先证两种分解的列空间相同 (都要等于 R ( A ) R(A) R ( A ) ,其中「R ( B 1 ) ⊆ R ( A ) R(B_1)\subseteq R(A) R ( B 1 ) ⊆ R ( A ) 」这一步需要 C 1 C_1 C 1 满射);然后把 B 1 B_1 B 1 、B 2 B_2 B 2 都看作 R ( A ) R(A) R ( A ) 的基,用过渡矩阵取 D D D ;最后用 B 2 B_2 B 2 的列无关把 C C C 的关系「约」出来。
解答(本编补,与 [3.2.3](sec32.md#sec-3.2.3) 正文一致)
第一步(列空间相等) :C 1 C_1 C 1 行满秩,故有右逆(定理 3.1.1 ),即 C 1 C_1 C 1 满射。任取 y = B 1 u ∈ R ( B 1 ) \boldsymbol{y}=B_1\boldsymbol{u}\in R(B_1) y = B 1 u ∈ R ( B 1 ) ,取 x \boldsymbol{x} x 使 C 1 x = u C_1\boldsymbol{x}=\boldsymbol{u} C 1 x = u ,则 y = B 1 C 1 x = A x ∈ R ( A ) \boldsymbol{y}=B_1C_1\boldsymbol{x}=A\boldsymbol{x}\in R(A) y = B 1 C 1 x = A x ∈ R ( A ) ——故 R ( B 1 ) ⊆ R ( A ) R(B_1)\subseteq R(A) R ( B 1 ) ⊆ R ( A ) ;反向 R ( A ) ⊆ R ( B 1 ) R(A)\subseteq R(B_1) R ( A ) ⊆ R ( B 1 ) 显然。于是 R ( B 1 ) = R ( A ) R(B_1)=R(A) R ( B 1 ) = R ( A ) ;同理 R ( B 2 ) = R ( A ) R(B_2)=R(A) R ( B 2 ) = R ( A ) 。
第二步(过渡矩阵) :B 1 B_1 B 1 、B 2 B_2 B 2 的列都是 r r r 维子空间 R ( A ) R(A) R ( A ) 的基。由过渡矩阵的存在唯一性(定义 1.3.5 、命题 1.3.1 ),存在可逆 D D D 使 B 1 = B 2 D B_1=B_2D B 1 = B 2 D 。
第三步(C C C 的关系) :B 2 C 2 = A = B 1 C 1 = B 2 D C 1 B_2C_2=A=B_1C_1=B_2DC_1 B 2 C 2 = A = B 1 C 1 = B 2 D C 1 ,移项得 B 2 ( C 2 − D C 1 ) = O B_2(C_2-DC_1)=O B 2 ( C 2 − D C 1 ) = O 。B 2 B_2 B 2 列满秩即列无关,B 2 X = O B_2X=O B 2 X = O 只有零解(按列用 定理 1.3.2 ),故 C 2 = D C 1 C_2=DC_1 C 2 = D C 1 ,即 C 1 = D − 1 C 2 C_1=D^{-1}C_2 C 1 = D − 1 C 2 。∎
B6 · 给定 B B B 后 C C C 唯一
设 A = B C A=BC A = BC 是满秩分解。证明:C = ( B H B ) − 1 B H A C=(B^{\mathrm H}B)^{-1}B^{\mathrm H}A C = ( B H B ) − 1 B H A (即 C C C 由 B B B 唯一确定)。
解答 · B6
B B B 列满秩 ⇒ \Rightarrow ⇒ B H B B^{\mathrm H}B B H B 可逆(命题 3.1.1 )。两边左乘 ( B H B ) − 1 B H (B^{\mathrm H}B)^{-1}B^{\mathrm H} ( B H B ) − 1 B H :
( B H B ) − 1 B H A = ( B H B ) − 1 B H B C = ( B H B ) − 1 ( B H B ) C = C ✓ (B^{\mathrm H}B)^{-1}B^{\mathrm H}A=(B^{\mathrm H}B)^{-1}B^{\mathrm H}BC=(B^{\mathrm H}B)^{-1}(B^{\mathrm H}B)C=C\ \checkmark ( B H B ) − 1 B H A = ( B H B ) − 1 B H BC = ( B H B ) − 1 ( B H B ) C = C ✓
∎
对偶结论 :给定 C C C 时 B = A C H ( C C H ) − 1 B=AC^{\mathrm H}(CC^{\mathrm H})^{-1} B = A C H ( C C H ) − 1 (同样的一句话验算)——见 3.2.2 的本编补注。
B7 · R ( A ) = R ( B ) R(A)=R(B) R ( A ) = R ( B )
证明:满秩分解 A = B C A=BC A = BC 中 R ( A ) = R ( B ) R(A)=R(B) R ( A ) = R ( B ) 。
提示 · B7
两个包含方向难度不同:一个方向连乘法都不用(B B B 的列组合);另一个方向要「解 C x = u C\boldsymbol{x}=\boldsymbol{u} C x = u 」,这正是 C C C 行满秩的用处。
解答 · B7
R ( A ) ⊆ R ( B ) R(A)\subseteq R(B) R ( A ) ⊆ R ( B ) :任取 y = A x = B ( C x ) \boldsymbol{y}=A\boldsymbol{x}=B(C\boldsymbol{x}) y = A x = B ( C x ) ,它是 B B B 的列的线性组合,故 y ∈ R ( B ) \boldsymbol{y}\in R(B) y ∈ R ( B ) 。(不需要任何前提。)
R ( B ) ⊆ R ( A ) R(B)\subseteq R(A) R ( B ) ⊆ R ( A ) :任取 y = B u ∈ R ( B ) \boldsymbol{y}=B\boldsymbol{u}\in R(B) y = B u ∈ R ( B ) 。因 C C C 行满秩,由 定理 3.1.1 它有右逆,即 C C C 作为映射 C n → C r \mathbb C^n\to\mathbb C^r C n → C r 是满射,故存在 x \boldsymbol{x} x 使 C x = u C\boldsymbol{x}=\boldsymbol{u} C x = u 。于是
y = B u = B C x = A x ∈ R ( A ) ✓ \boldsymbol{y}=B\boldsymbol{u}=BC\boldsymbol{x}=A\boldsymbol{x}\in R(A)\ \checkmark y = B u = BC x = A x ∈ R ( A ) ✓
故 R ( A ) = R ( B ) R(A)=R(B) R ( A ) = R ( B ) 。∎
推论 :dim R ( A ) = dim R ( B ) = r \dim R(A)=\dim R(B)=r dim R ( A ) = dim R ( B ) = r ——满秩分解把列空间「显式地」交给了 B B B 。∎
B8 · A X A = A AXA=A A X A = A 与 C X B = I r CXB=I_r CXB = I r 等价
设 A = B C ∈ C r m × n A=BC\in\mathbb C^{m\times n}_r A = BC ∈ C r m × n 是满秩分解。证明:对方阵 X ∈ C n × m X\in\mathbb C^{n\times m} X ∈ C n × m ,
A X A = A ⟺ C X B = I r . AXA=A\iff CXB=I_r . A X A = A ⟺ CXB = I r .
提示 · B8
「⟸ \Longleftarrow ⟸ 」直接把 C X B CXB CXB 换成 I r I_r I r ;「⟹ \Longrightarrow ⟹ 」分别左乘 B B B 的左逆、右乘 C C C 的右逆——两条单侧逆的存在性由 定理 3.1.1 与 命题 3.1.1 保证。
解答 · B8
(⟸ \Longleftarrow ⟸ ) 若 C X B = I r CXB=I_r CXB = I r ,则
A X A = B C X B C = B ( C X B ) C = B C = A ✓ AXA=BCXBC=B\,(CXB)\,C=BC=A\ \checkmark A X A = BCXBC = B ( CXB ) C = BC = A ✓
(⟹ \Longrightarrow ⟹ ) 设 B C X B C = B C BCXBC=BC BCXBC = BC 。先左乘 B B B 的左逆 B L = ( B H B ) − 1 B H B^{\mathrm L}=(B^{\mathrm H}B)^{-1}B^{\mathrm H} B L = ( B H B ) − 1 B H :
B L B C X B C = B L B C ⟹ C X B C = C ; B^{\mathrm L}BCXBC=B^{\mathrm L}BC\ \Longrightarrow\ CXBC=C ; B L BCXBC = B L BC ⟹ CXBC = C ;
再右乘 C C C 的右逆 C R = C H ( C C H ) − 1 C^{\mathrm R}=C^{\mathrm H}(CC^{\mathrm H})^{-1} C R = C H ( C C H ) − 1 :
C X B C C R = C C R ⟹ C X B = I r ✓ CXBCC^{\mathrm R}=CC^{\mathrm R}\ \Longrightarrow\ CXB=I_r\ \checkmark CXBC C R = C C R ⟹ CXB = I r ✓
两个单侧逆都存在(B H B B^{\mathrm H}B B H B 、C C H CC^{\mathrm H} C C H 可逆,命题 3.1.1 )。∎
B9 · 行满秩 ⟹ \Longrightarrow ⟹ A A H AA^{\mathrm H} A A H 正定
证明:若 A ∈ C m × n A\in\mathbb C^{m\times n} A ∈ C m × n 行满秩,则 A A H AA^{\mathrm H} A A H 是正定 Hermite 矩阵(用第一章的正定定义 定义 1.4.3 与 命题 3.1.1 各自完成一部分)。
解答 · B9
Hermite 性 :直接算 ( A A H ) H = A ( A H ) H = A A H (AA^{\mathrm H})^{\mathrm H}=A\,(A^{\mathrm H})^{\mathrm H}=AA^{\mathrm H} ( A A H ) H = A ( A H ) H = A A H ✓(共轭转置取两次回到自身)。
正定性 :对任意非零 x ∈ C m \boldsymbol{x}\in\mathbb C^m x ∈ C m ,
x H ( A A H ) x = ( A H x ) H ( A H x ) = ∥ A H x ∥ 2 ≥ 0. \boldsymbol{x}^{\mathrm H}(AA^{\mathrm H})\boldsymbol{x}=(A^{\mathrm H}\boldsymbol{x})^{\mathrm H}(A^{\mathrm H}\boldsymbol{x})=\|A^{\mathrm H}\boldsymbol{x}\|^2\ge0 . x H ( A A H ) x = ( A H x ) H ( A H x ) = ∥ A H x ∥ 2 ≥ 0.
等号成立 ⟺ \iff ⟺ A H x = 0 A^{\mathrm H}\boldsymbol{x}=\boldsymbol{0} A H x = 0 。而 A A A 行满秩 ⟺ \iff ⟺ rank ( A H ) = m \operatorname{rank}(A^{\mathrm H})=m rank ( A H ) = m ⟺ \iff ⟺ A H A^{\mathrm H} A H 列满秩 ⟺ \iff ⟺ A H A^{\mathrm H} A H 单射(推论 2.2.2 或 定理 2.2.2 ),故 A H x = 0 ⟺ x = 0 A^{\mathrm H}\boldsymbol{x}=\boldsymbol{0}\iff\boldsymbol{x}=\boldsymbol{0} A H x = 0 ⟺ x = 0 。于是对非零 x \boldsymbol{x} x 恒有 x H A A H x > 0 \boldsymbol{x}^{\mathrm H}AA^{\mathrm H}\boldsymbol{x}\gt0 x H A A H x > 0 ,A A H AA^{\mathrm H} A A H 正定。∎
对照 :3.1.4 用 命题 3.1.1 说明「A A A 行满秩 ⇒ \Rightarrow ⇒ A A H AA^{\mathrm H} A A H 可逆」;本题进一步指出可逆的类型 ——它是正定 Hermite 矩阵(可逆只是正定的一个推论)。
B10 · 相抵不唯一的「吸收」公式
验证 3.1.1 的不唯一性构造:对 Γ = [ I r 0 0 0 ] \Gamma=\begin{bmatrix}I_r&0\\0&0\end{bmatrix} Γ = [ I r 0 0 0 ] 与任意适当阶数矩阵 ∗ * ∗ ,证明
[ I r ∗ 0 I ] Γ = Γ , Γ [ I r 0 ∗ I ] = Γ , \begin{bmatrix}I_r&*\\0&I\end{bmatrix}\Gamma=\Gamma,\qquad \Gamma\begin{bmatrix}I_r&0\\*&I\end{bmatrix}=\Gamma, [ I r 0 ∗ I ] Γ = Γ , Γ [ I r ∗ 0 I ] = Γ ,
并说明这两个三角块都可逆。
解答 · B10
第一式 (m × m m\times m m × m 的块左乘):记 Γ = [ I r 0 0 0 ] \Gamma=\begin{bmatrix}I_r&0\\0&0\end{bmatrix} Γ = [ I r 0 0 0 ] (分块尺寸 r + ( m − r ) r+(m-r) r + ( m − r ) 、r + ( n − r ) r+(n-r) r + ( n − r ) )。按分块乘法:
[ I r ∗ 0 I ] [ I r 0 0 0 ] = [ I r ⋅ I r + ∗ ⋅ 0 I r ⋅ 0 + ∗ ⋅ 0 0 ⋅ I r + I ⋅ 0 0 ⋅ 0 + I ⋅ 0 ] = [ I r 0 0 0 ] = Γ ✓ \begin{bmatrix}I_r&*\\0&I\end{bmatrix}\begin{bmatrix}I_r&0\\0&0\end{bmatrix}=\begin{bmatrix}I_r\cdot I_r+*\cdot 0&\ I_r\cdot 0+*\cdot 0\\ 0\cdot I_r+I\cdot 0&\ 0\cdot0+I\cdot0\end{bmatrix}=\begin{bmatrix}I_r&0\\0&0\end{bmatrix}=\Gamma\ \checkmark [ I r 0 ∗ I ] [ I r 0 0 0 ] = [ I r ⋅ I r + ∗ ⋅ 0 0 ⋅ I r + I ⋅ 0 I r ⋅ 0 + ∗ ⋅ 0 0 ⋅ 0 + I ⋅ 0 ] = [ I r 0 0 0 ] = Γ ✓
(关键:Γ \Gamma Γ 的后 n − r n-r n − r 列全零,故 ∗ * ∗ 乘上去没有贡献。)
第二式 (n × n n\times n n × n 的块右乘):Γ \Gamma Γ 的后 m − r m-r m − r 行 全零,故
[ I r 0 0 0 ] [ I r 0 ∗ I ] = [ I r 0 0 0 ] = Γ ✓ \begin{bmatrix}I_r&0\\0&0\end{bmatrix}\begin{bmatrix}I_r&0\\*&I\end{bmatrix}=\begin{bmatrix}I_r&0\\0&0\end{bmatrix}=\Gamma\ \checkmark [ I r 0 0 0 ] [ I r ∗ 0 I ] = [ I r 0 0 0 ] = Γ ✓
可逆性 :两块都是「单位块 + 严格三角块」,行列式为 1 ⋅ det I = 1 ≠ 0 1\cdot\det I=1\neq0 1 ⋅ det I = 1 = 0 ,可逆。
用法 :由两式得 P − 1 Γ Q = ( P − 1 [ I r ∗ 0 I ] ) Γ ( [ I r 0 ∗ I ] Q ) P^{-1}\Gamma Q=\big(P^{-1}\begin{bmatrix}I_r&*\\0&I\end{bmatrix}\big)\Gamma\big(\begin{bmatrix}I_r&0\\*&I\end{bmatrix}Q\big) P − 1 Γ Q = ( P − 1 [ I r 0 ∗ I ] ) Γ ( [ I r ∗ 0 I ] Q ) ,右边仍是合法的相抵分解式——取不同的 ∗ * ∗ 即得不同的 P P P 、Q Q Q 。∎
B11 · 两种 QR 分解恰好相差一个对角酉阵
设 A A A 为 n n n 阶满秩复方阵,A = Q 1 R 1 = Q 2 R 2 A=Q_1R_1=Q_2R_2 A = Q 1 R 1 = Q 2 R 2 是两种 QR 分解(Q i Q_i Q i 酉、R i R_i R i 可逆上三角,不要求正线)。证明:存在对角酉阵 W = diag ( w 1 , ⋯ , w n ) W=\operatorname{diag}(w_1,\cdots,w_n) W = diag ( w 1 , ⋯ , w n ) (∣ w i ∣ = 1 |w_i|=1 ∣ w i ∣ = 1 )使 Q 2 = Q 1 W Q_2=Q_1W Q 2 = Q 1 W 、R 2 = W H R 1 R_2=W^{\mathrm H}R_1 R 2 = W H R 1 。
提示 · B11
考虑 W = R 2 R 1 − 1 = Q 2 H Q 1 W=R_2R_1^{-1}=Q_2^{\mathrm H}Q_1 W = R 2 R 1 − 1 = Q 2 H Q 1 ,证它「既酉又上三角」;再对列归纳证上三角酉阵必为对角。
解答 · B11
W W W 良定且两种身份 :R 1 R_1 R 1 可逆(A A A 满秩 ⇒ rank R 1 = n \Rightarrow\operatorname{rank}R_1=n ⇒ rank R 1 = n ,上三角满秩 ⟺ \iff ⟺ 对角全非零)。由 Q 1 R 1 = Q 2 R 2 Q_1R_1=Q_2R_2 Q 1 R 1 = Q 2 R 2 左乘 Q 1 H Q_1^{\mathrm H} Q 1 H 、右乘 R 1 − 1 R_1^{-1} R 1 − 1 :
W : = Q 1 H Q 2 = R 1 R 2 − 1 . W:=Q_1^{\mathrm H}Q_2=R_1R_2^{-1} . W := Q 1 H Q 2 = R 1 R 2 − 1 .
左边是酉阵之积,酉(正交 / 酉矩阵的运算性质 2.5.2 );右边是上三角阵之积再求逆,上三角。上三角且酉 ⇒ \Rightarrow ⇒ 对角、∣ w i i ∣ = 1 |w_{ii}|=1 ∣ w ii ∣ = 1 (对列归纳:第 1 1 1 列只有 w 11 w_{11} w 11 非零而范数为 1 1 1 ;设前 j − 1 j-1 j − 1 列已是 w i i e i w_{ii}\boldsymbol{e}_i w ii e i ,第 j j j 列与它们的内积为 w i i ‾ w i j = 0 \overline{w_{ii}}w_{ij}=0 w ii w ij = 0 ,逐个逼出 w i j = 0 w_{ij}=0 w ij = 0 (i < j i<j i < j ),范数给出 ∣ w j j ∣ = 1 |w_{jj}|=1 ∣ w jj ∣ = 1 )。于是 W W W 是对角酉阵。
收尾 :W = Q 1 H Q 2 ⇒ Q 2 = Q 1 W W=Q_1^{\mathrm H}Q_2\Rightarrow Q_2=Q_1W W = Q 1 H Q 2 ⇒ Q 2 = Q 1 W (两边左乘 Q 1 Q_1 Q 1 );W = R 1 R 2 − 1 ⇒ R 2 = W − 1 R 1 = W H R 1 W=R_1R_2^{-1}\Rightarrow R_2=W^{-1}R_1=W^{\mathrm H}R_1 W = R 1 R 2 − 1 ⇒ R 2 = W − 1 R 1 = W H R 1 (酉阵 W − 1 = W H W^{-1}=W^{\mathrm H} W − 1 = W H ),其第 i i i 个对角元为 w i ‾ r i i ( 1 ) \overline{w_i}\,r_{ii}^{(1)} w i r ii ( 1 ) 。∎
推论 :实数情形 w i = ± 1 w_i=\pm1 w i = ± 1 ,恰 2 n 2^n 2 n 种;再要求两个分解都正线(r i i ( 1 ) > 0 r_{ii}^{(1)}>0 r ii ( 1 ) > 0 、w i ‾ r i i ( 1 ) > 0 \overline{w_i}r_{ii}^{(1)}>0 w i r ii ( 1 ) > 0 )则逼出 w i = 1 w_i=1 w i = 1 ——唯一性(定理 3.4.1 )。
B12 · 上三角酉阵必为对角阵
证明:若 W = ( w i j ) W=(w_{ij}) W = ( w ij ) 既是 n n n 阶上三角阵(w i j = 0 , i > j w_{ij}=0,\ i>j w ij = 0 , i > j )又是酉矩阵,则 w i j = 0 w_{ij}=0 w ij = 0 (i < j i<j i < j )且 ∣ w i i ∣ = 1 |w_{ii}|=1 ∣ w ii ∣ = 1 。
解答 · B12
对列 j j j 归纳,断言「第 j j j 列 = w j j e j =w_{jj}\boldsymbol{e}_j = w jj e j 且 ∣ w j j ∣ = 1 |w_{jj}|=1 ∣ w jj ∣ = 1 」。
j = 1 j=1 j = 1 :上三角使第一列只有 w 11 w_{11} w 11 可能非零;酉矩阵各列范数为 1 1 1 (W H W = I W^{\mathrm H}W=I W H W = I 的 ( 1 , 1 ) (1,1) ( 1 , 1 ) 元),故 ∣ w 11 ∣ = 1 |w_{11}|=1 ∣ w 11 ∣ = 1 ,第 1 1 1 列 = w 11 e 1 =w_{11}\boldsymbol{e}_1 = w 11 e 1 。
j ⇒ j + 1 j\Rightarrow j+1 j ⇒ j + 1 :第 j + 1 j+1 j + 1 列的非零分量只在行 1 ⋯ j + 1 1\cdots j+1 1 ⋯ j + 1 。对每个 i ≤ j i\le j i ≤ j ,酉矩阵不同列正交(W H W = I W^{\mathrm H}W=I W H W = I 的 ( i , j + 1 ) (i,j+1) ( i , j + 1 ) 元):
0 = ∑ k ≤ i w k i ‾ w k , j + 1 = w i i ‾ w i , j + 1 ( 归纳假设:第 i 列 = w i i e i , k > i 项为零 ) , 0=\sum_{k\le i}\overline{w_{ki}}\,w_{k,j+1}=\overline{w_{ii}}\,w_{i,j+1}\quad(\text{归纳假设:第 }i\text{ 列 }=w_{ii}\boldsymbol{e}_i,\ k>i\text{ 项为零}), 0 = k ≤ i ∑ w ki w k , j + 1 = w ii w i , j + 1 ( 归纳假设:第 i 列 = w ii e i , k > i 项为零 ) ,
而 w i i ‾ ≠ 0 \overline{w_{ii}}\neq0 w ii = 0 ,故 w i , j + 1 = 0 w_{i,j+1}=0 w i , j + 1 = 0 。于是第 j + 1 j+1 j + 1 列只剩 w j + 1 , j + 1 w_{j+1,j+1} w j + 1 , j + 1 ,范数 1 1 1 给出 ∣ w j + 1 , j + 1 ∣ = 1 |w_{j+1,j+1}|=1 ∣ w j + 1 , j + 1 ∣ = 1 。∎
(这条小事实是 定理 3.4.1 唯一性与 B11 的共同引擎;它对「行满 + + + 酉」版本(下三角)同样成立,对 A H A^{\mathrm H} A H 说一遍即可。)
B13 · QR 迭代的相似链与「卡死」分析
(a) 证明 QR 迭代满足 A k = ( Q 1 Q 2 ⋯ Q k − 1 ) H A ( Q 1 Q 2 ⋯ Q k − 1 ) A_k=(Q_1Q_2\cdots Q_{k-1})^{\mathrm H}A(Q_1Q_2\cdots Q_{k-1}) A k = ( Q 1 Q 2 ⋯ Q k − 1 ) H A ( Q 1 Q 2 ⋯ Q k − 1 ) ,从而每个 A k A_k A k 与 A A A 相似;
(b) 证明:旋转阵 A = [ 0 − 1 1 0 ] A=\begin{bmatrix}0&-1\\1&0\end{bmatrix} A = [ 0 1 − 1 0 ] 的迭代序列恒为 A k = A A_k=A A k = A ,永不收敛到上三角阵。
解答 · B13
(a) 归纳。k = 1 k=1 k = 1 :空积为 I I I ,A 1 = A A_1=A A 1 = A ✓。设 A k = W k − 1 H A W k − 1 A_k=W_{k-1}^{\mathrm H}AW_{k-1} A k = W k − 1 H A W k − 1 ,W k − 1 = Q 1 ⋯ Q k − 1 W_{k-1}=Q_1\cdots Q_{k-1} W k − 1 = Q 1 ⋯ Q k − 1 。由迭代式与 R k = Q k H A k R_k=Q_k^{\mathrm H}A_k R k = Q k H A k :
A k + 1 = R k Q k = Q k H A k Q k = Q k H W k − 1 H A W k − 1 Q k = ( W k − 1 Q k ) H A ( W k − 1 Q k ) = W k H A W k , A_{k+1}=R_kQ_k=Q_k^{\mathrm H}A_kQ_k=Q_k^{\mathrm H}W_{k-1}^{\mathrm H}AW_{k-1}Q_k=(W_{k-1}Q_k)^{\mathrm H}A(W_{k-1}Q_k)=W_k^{\mathrm H}AW_k, A k + 1 = R k Q k = Q k H A k Q k = Q k H W k − 1 H A W k − 1 Q k = ( W k − 1 Q k ) H A ( W k − 1 Q k ) = W k H A W k ,
其中 W k = W k − 1 Q k = Q 1 ⋯ Q k W_k=W_{k-1}Q_k=Q_1\cdots Q_k W k = W k − 1 Q k = Q 1 ⋯ Q k (酉阵之积仍酉)。故 A k ∼ A A_k\sim A A k ∼ A ,特征值全程不变(相似不变量 2.4.4 )。
(b) A A A 是正交矩阵,其正线 QR 分解是 A = A ⋅ I A=A\cdot I A = A ⋅ I :唯一性(定理 3.4.1 )保证 GS 只能给出这一组——事实上 GS 对 A A A 的两列 e 2 , e 1 \boldsymbol{e}_2,\boldsymbol{e}_1 e 2 , e 1 原样通过(已标准正交),Q 1 = A Q_1=A Q 1 = A 、R 1 = I R_1=I R 1 = I 。于是
A 2 = R 1 Q 1 = I ⋅ A = A , A_2=R_1Q_1=I\cdot A=A, A 2 = R 1 Q 1 = I ⋅ A = A ,
归纳得 A k = A A_k=A A k = A 对一切 k k k 成立;而 A A A 不是上三角阵,序列不收敛于上三角阵。根源:特征值 ± i \pm\mathrm i ± i 等模 ,违反收敛的充分条件(QR 算法 )。∎
B14 · 重写定理 3.5.1 的归纳证明(不看正文)
不翻正文,独立写出 Schur 引理(复数域版本)的完整证明:任意 n n n 阶复方阵相似于上三角阵,且对角元为其特征值。
解答 · B14
对 n n n 归纳 。n = 1 n=1 n = 1 平凡。设对 n − 1 n-1 n − 1 成立。A A A 的特征多项式在 C \mathbb C C 上有根 λ 1 \lambda_1 λ 1 (代数基本定理),取特征向量 x 1 ≠ 0 \boldsymbol{x}_1\neq\boldsymbol{0} x 1 = 0 (A x 1 = λ 1 x 1 A\boldsymbol{x}_1=\lambda_1\boldsymbol{x}_1 A x 1 = λ 1 x 1 ),扩充为 C n \mathbb C^n C n 的基(推论 1.3.1 )并拼成可逆阵 P 1 = [ x 1 , ⋯ ] P_1=[\boldsymbol{x}_1,\cdots] P 1 = [ x 1 , ⋯ ] 。则 P 1 − 1 A P 1 P_1^{-1}AP_1 P 1 − 1 A P 1 的第一列为 P 1 − 1 A x 1 = λ 1 e 1 P_1^{-1}A\boldsymbol{x}_1=\lambda_1\boldsymbol{e}_1 P 1 − 1 A x 1 = λ 1 e 1 :
P 1 − 1 A P 1 = [ λ 1 β 0 A 1 ] . P_1^{-1}AP_1=\begin{bmatrix}\lambda_1&\boldsymbol{\beta}\\\boldsymbol{0}&A_1\end{bmatrix}. P 1 − 1 A P 1 = [ λ 1 0 β A 1 ] .
对 A 1 A_1 A 1 (n − 1 n-1 n − 1 阶)用归纳假设:存在可逆 Q Q Q 使 Q − 1 A 1 Q = Λ 1 Q^{-1}A_1Q=\Lambda_1 Q − 1 A 1 Q = Λ 1 上三角。令 P = P 1 diag ( 1 , Q ) P=P_1\operatorname{diag}(1,Q) P = P 1 diag ( 1 , Q ) ,则
P − 1 A P = [ λ 1 β Q 0 Λ 1 ] P^{-1}AP=\begin{bmatrix}\lambda_1&\boldsymbol{\beta}Q\\\boldsymbol{0}&\Lambda_1\end{bmatrix} P − 1 A P = [ λ 1 0 β Q Λ 1 ]
上三角。对角元:分块三角行列式给 f A ( λ ) = ( λ − λ 1 ) f A 1 ( λ ) f_A(\lambda)=(\lambda-\lambda_1)f_{A_1}(\lambda) f A ( λ ) = ( λ − λ 1 ) f A 1 ( λ ) ,故对角元全体(计重数)= λ 1 =\lambda_1 = λ 1 加 A 1 A_1 A 1 的特征值 = A =A = A 的特征值。∎
自查两处 :① 第一步为什么必能取到特征向量?(复数域上特征多项式必有根——这是「复」字在定理里干的活。)② 收尾为什么对角元恰是 A A A 的特征值?(相似不改特征多项式 + + + 分块三角行列式。)
B15 · 直接构造酉矩阵(不走定理 3.5.1 的路)
不经过定理 3.5.1,直接用「单位特征向量 + + + Gram–Schmidt 扩充」归纳证明定理 3.5.2:任意复方阵酉相似于上三角阵。
提示 · B15
取单位 特征向量 u 1 \boldsymbol{u}_1 u 1 ,用基扩充 + + + GS 把它扩成标准正交基;[ u 1 , ⋯ ] [\boldsymbol{u}_1,\cdots] [ u 1 , ⋯ ] 拼成的矩阵是酉矩阵,算出它的第一列象。
解答 · B15
对 n n n 归纳 。n = 1 n=1 n = 1 平凡。设对 n − 1 n-1 n − 1 成立。取 A A A 的特征值 λ 1 \lambda_1 λ 1 与单位 特征向量 u 1 \boldsymbol{u}_1 u 1 (∥ u 1 ∥ = 1 \|\boldsymbol{u}_1\|=1 ∥ u 1 ∥ = 1 ;特征向量除以范数仍是特征向量)。把 u 1 \boldsymbol{u}_1 u 1 扩充为 C n \mathbb C^n C n 的基(推论 1.3.1 ),再 GS 并单位化(Gram–Schmidt 正交化 ;「张成不变」保证 u 1 \boldsymbol{u}_1 u 1 原样保留),得标准正交基 u 1 , u 2 , ⋯ , u n \boldsymbol{u}_1,\boldsymbol{u}_2,\cdots,\boldsymbol{u}_n u 1 , u 2 , ⋯ , u n 。令 U 1 = [ u 1 , ⋯ , u n ] U_1=[\boldsymbol{u}_1,\cdots,\boldsymbol{u}_n] U 1 = [ u 1 , ⋯ , u n ] :列标准正交的方阵,酉(列向量刻画 2.5.1 )。计算 U 1 H A U 1 U_1^{\mathrm H}AU_1 U 1 H A U 1 的第一列:
U 1 H A u 1 = λ 1 U 1 H u 1 = λ 1 e 1 , U_1^{\mathrm H}A\boldsymbol{u}_1=\lambda_1U_1^{\mathrm H}\boldsymbol{u}_1=\lambda_1\boldsymbol{e}_1, U 1 H A u 1 = λ 1 U 1 H u 1 = λ 1 e 1 ,
故 U 1 H A U 1 = [ λ 1 β 0 A 1 ] U_1^{\mathrm H}AU_1=\begin{bmatrix}\lambda_1&\boldsymbol{\beta}\\\boldsymbol{0}&A_1\end{bmatrix} U 1 H A U 1 = [ λ 1 0 β A 1 ] 。对 A 1 A_1 A 1 用归纳假设:存在酉阵 V V V 使 V H A 1 V = Λ 1 V^{\mathrm H}A_1V=\Lambda_1 V H A 1 V = Λ 1 上三角。令 U = U 1 diag ( 1 , V ) U=U_1\operatorname{diag}(1,V) U = U 1 diag ( 1 , V ) (酉阵之积,酉),则
U H A U = diag ( 1 , V H ) [ λ 1 β 0 A 1 ] diag ( 1 , V ) = [ λ 1 β V 0 Λ 1 ] U^{\mathrm H}AU=\operatorname{diag}(1,V^{\mathrm H})\begin{bmatrix}\lambda_1&\boldsymbol{\beta}\\\boldsymbol{0}&A_1\end{bmatrix}\operatorname{diag}(1,V)=\begin{bmatrix}\lambda_1&\boldsymbol{\beta}V\\\boldsymbol{0}&\Lambda_1\end{bmatrix} U H A U = diag ( 1 , V H ) [ λ 1 0 β A 1 ] diag ( 1 , V ) = [ λ 1 0 β V Λ 1 ]
上三角。∎
与正文路线的对照 :正文(定理 3.5.2 )「先相似再 QR 升级」是一步到位的代数戏法;本题是「边归纳边正交化」的施工图。两条路在二阶情形的差别值得一眼:前者需要先有可逆 P P P ,后者从第一个特征向量起就全程保范。
B16 · 谱映射的一半与反例的另一半
(a) 不用 Schur 分解,证明:若 α \boldsymbol{\alpha} α 是 A A A 的属于 λ \lambda λ 的特征向量,则 α \boldsymbol{\alpha} α 是 φ ( A ) \varphi(A) φ ( A ) 的属于 φ ( λ ) \varphi(\lambda) φ ( λ ) 的特征向量(注 4),且 φ ( λ 1 ) , ⋯ , φ ( λ n ) \varphi(\lambda_1),\cdots,\varphi(\lambda_n) φ ( λ 1 ) , ⋯ , φ ( λ n ) (计重数)都是 φ ( A ) \varphi(A) φ ( A ) 的特征值;
(b) 举反例说明 φ ( A ) \varphi(A) φ ( A ) 的特征向量不必是 A A A 的特征向量。
解答 · B16
(a) 由 A α = λ α A\boldsymbol{\alpha}=\lambda\boldsymbol{\alpha} A α = λ α 归纳得 A k α = λ k α A^k\boldsymbol{\alpha}=\lambda^k\boldsymbol{\alpha} A k α = λ k α (k ≥ 0 k\ge0 k ≥ 0 ,含 A 0 = I A^0=I A 0 = I )。线性组合:
φ ( A ) α = ∑ i a i A i α = ( ∑ i a i λ i ) α = φ ( λ ) α , \varphi(A)\boldsymbol{\alpha}=\sum_ia_iA^i\boldsymbol{\alpha}=\Big(\sum_ia_i\lambda^i\Big)\boldsymbol{\alpha}=\varphi(\lambda)\boldsymbol{\alpha}, φ ( A ) α = i ∑ a i A i α = ( i ∑ a i λ i ) α = φ ( λ ) α ,
即 α \boldsymbol{\alpha} α 是 φ ( A ) \varphi(A) φ ( A ) 的特征向量(特征值 φ ( λ ) \varphi(\lambda) φ ( λ ) )。于是对每个特征值 λ i \lambda_i λ i 取一个特征向量 α i \boldsymbol{\alpha}_i α i ,φ ( λ i ) \varphi(\lambda_i) φ ( λ i ) 就有证人 地成为 φ ( A ) \varphi(A) φ ( A ) 的特征值;重数叙述「计重数」的完整版需数齐 n n n 个,那一步用三角化最省(推论 3.5.1 的证明)。
(b) A = diag ( 1 , − 1 ) A=\operatorname{diag}(1,-1) A = diag ( 1 , − 1 ) 、φ ( λ ) = λ 2 \varphi(\lambda)=\lambda^2 φ ( λ ) = λ 2 :φ ( A ) = A 2 = I \varphi(A)=A^2=I φ ( A ) = A 2 = I ,任何非零向量都是 I I I 的特征向量(特征值 1 1 1 ),但 e 1 + e 2 \boldsymbol{e}_1+\boldsymbol{e}_2 e 1 + e 2 不是 A A A 的特征向量。∎
要点 :继承是单向的——φ \varphi φ 可以把不同的特征值黏到同一个值上,信息只减不增(这正是最小多项式要「计账」的对象)。
B17 · 重写定理 3.5.5(不看正文)
独立证明最小多项式的两条性质:(1) m A m_A m A 唯一、整除 A A A 的任一零化多项式(特别地 m A ∣ f A m_A\mid f_A m A ∣ f A );(2) f A f_A f A 与 m A m_A m A 同根(不计重数)。
解答 · B17
(1) 存在性 (备查):f A f_A f A 零化 A A A (H–C,定理 3.5.4 ),故零化多项式之族非空,取次数最小者再除以首项系数(仍零化、次数不变)得首一的最小者 m A m_A m A 。
整除 :设 g g g 任一零化多项式。带余除法(定理 2.1.3 ):g = q m A + r g=qm_A+r g = q m A + r ,deg r < deg m A \deg r<\deg m_A deg r < deg m A 或 r = 0 r=0 r = 0 。代入矩阵并注意同一矩阵的多项式可交换:
O = g ( A ) = q ( A ) m A ( A ) ⏟ = O + r ( A ) ⟹ r ( A ) = O . O=g(A)=q(A)\underbrace{m_A(A)}_{=O}+r(A)\ \Longrightarrow\ r(A)=O . O = g ( A ) = q ( A ) = O m A ( A ) + r ( A ) ⟹ r ( A ) = O .
若 r ≠ 0 r\neq0 r = 0 ,则 r r r (除以首项系数后首一)是次数更低 的零化多项式,与最小性矛盾。故 r = 0 r=0 r = 0 ,m A ∣ g m_A\mid g m A ∣ g ;取 g = f A g=f_A g = f A 得 m A ∣ f A m_A\mid f_A m A ∣ f A 。
唯一 :若 m 1 , m 2 m_1,m_2 m 1 , m 2 都是最小的首一零化多项式,则 m 1 ∣ m 2 m_1\mid m_2 m 1 ∣ m 2 且 m 2 ∣ m 1 m_2\mid m_1 m 2 ∣ m 1 ;两者同次(都最小)且首一,故 m 1 = m 2 m_1=m_2 m 1 = m 2 (互相整除的同次首一多项式只能相差 1 1 1 倍)。
(2) 同根 :「m A m_A m A 的根 ⊆ f A \subseteq f_A ⊆ f A 的根」:m A ∣ f A m_A\mid f_A m A ∣ f A ,根只会少不会多。「f A f_A f A 的根 ⊆ m A \subseteq m_A ⊆ m A 的根」:设 λ 0 \lambda_0 λ 0 是 A A A 的特征值,取特征向量 α ≠ 0 \boldsymbol{\alpha}\neq\boldsymbol{0} α = 0 。由注 4(B16 的 (a)):
0 = m A ( A ) α = m A ( λ 0 ) α ⟹ m A ( λ 0 ) = 0. ■ \boldsymbol{0}=m_A(A)\boldsymbol{\alpha}=m_A(\lambda_0)\boldsymbol{\alpha}\ \Longrightarrow\ m_A(\lambda_0)=0 .\qquad\blacksquare 0 = m A ( A ) α = m A ( λ 0 ) α ⟹ m A ( λ 0 ) = 0. ■
一处易错 :不能用「f A ( A ) = O ⇒ f A = m A f_A(A)=O\Rightarrow f_A=m_A f A ( A ) = O ⇒ f A = m A 」——零化多项式有无数个(注 5),f A f_A f A 只是其中之一;同根但重数可以差(I 2 I_2 I 2 :f A = ( λ − 1 ) 2 f_A=(\lambda-1)^2 f A = ( λ − 1 ) 2 ,m A = λ − 1 m_A=\lambda-1 m A = λ − 1 )。
C1 · 相抵分解全流程(2 × 3 2\times3 2 × 3 、秩 1 1 1 )
求 A = [ 1 2 3 2 4 6 ] A=\begin{bmatrix}1&2&3\\ 2&4&6\end{bmatrix} A = [ 1 2 2 4 3 6 ] 的相抵分解式:写出 Γ \Gamma Γ 、E E E 、F F F (即 P P P 与 Q − 1 Q^{-1} Q − 1 ),再写出 P − 1 P^{-1} P − 1 、Q Q Q 并回代验证。
解答 · C1
第一步(定秩) :第二行 = 2 × =2\times = 2 × 第一行,故 r = 1 r=1 r = 1 ,Γ = [ 1 0 0 0 0 0 ] \Gamma=\begin{bmatrix}1&0&0\\ 0&0&0\end{bmatrix} Γ = [ 1 0 0 0 0 0 ] 。
第二步(建分块矩阵并消元) :挂上 [ A I 2 I 3 ∗ ] \begin{bmatrix}A&I_2\\ I_3&*\end{bmatrix} [ A I 3 I 2 ∗ ] ,只动前 2 2 2 行、前 3 3 3 列。
行变换(记总积 E E E ):r 2 ← r 2 − 2 r 1 r_2\leftarrow r_2-2r_1 r 2 ← r 2 − 2 r 1 ,即 E = [ 1 0 − 2 1 ] E=\begin{bmatrix}1&0\\ -2&1\end{bmatrix} E = [ 1 − 2 0 1 ] ;
列变换(记总积 F F F ):c 2 ← c 2 − 2 c 1 c_2\leftarrow c_2-2c_1 c 2 ← c 2 − 2 c 1 、c 3 ← c 3 − 3 c 1 c_3\leftarrow c_3-3c_1 c 3 ← c 3 − 3 c 1 ,即 F = [ 1 − 2 − 3 0 1 0 0 0 1 ] F=\begin{bmatrix}1&-2&-3\\ 0&1&0\\ 0&0&1\end{bmatrix} F = 1 0 0 − 2 1 0 − 3 0 1 。
第三步(回代验算,用 E A F = Γ EAF=\Gamma E A F = Γ ) :
E A F = [ 1 0 − 2 1 ] [ 1 2 3 2 4 6 ] [ 1 − 2 − 3 0 1 0 0 0 1 ] = [ 1 2 3 0 0 0 ] [ 1 − 2 − 3 0 1 0 0 0 1 ] = [ 1 0 0 0 0 0 ] = Γ ✓ EAF=\begin{bmatrix}1&0\\-2&1\end{bmatrix}\begin{bmatrix}1&2&3\\2&4&6\end{bmatrix}\begin{bmatrix}1&-2&-3\\0&1&0\\0&0&1\end{bmatrix} =\begin{bmatrix}1&2&3\\0&0&0\end{bmatrix}\begin{bmatrix}1&-2&-3\\0&1&0\\0&0&1\end{bmatrix}=\begin{bmatrix}1&0&0\\0&0&0\end{bmatrix}=\Gamma\ \checkmark E A F = [ 1 − 2 0 1 ] [ 1 2 2 4 3 6 ] 1 0 0 − 2 1 0 − 3 0 1 = [ 1 0 2 0 3 0 ] 1 0 0 − 2 1 0 − 3 0 1 = [ 1 0 0 0 0 0 ] = Γ ✓
第四步(写成 A = P − 1 Γ Q A=P^{-1}\Gamma Q A = P − 1 Γ Q ) :由 E A F = Γ EAF=\Gamma E A F = Γ 得 A = E − 1 Γ F − 1 A=E^{-1}\Gamma F^{-1} A = E − 1 Γ F − 1 ,与 A = P − 1 Γ Q A=P^{-1}\Gamma Q A = P − 1 Γ Q 对照(相抵变换求 P P P 与 Q − 1 Q^{-1} Q − 1 ):
P − 1 = E − 1 = [ 1 0 2 1 ] , Q = F − 1 = [ 1 2 3 0 1 0 0 0 1 ] , P = E , Q − 1 = F \boxed{\,P^{-1}=E^{-1}=\begin{bmatrix}1&0\\ 2&1\end{bmatrix},\qquad Q=F^{-1}=\begin{bmatrix}1&2&3\\ 0&1&0\\ 0&0&1\end{bmatrix},\qquad P=E,\qquad Q^{-1}=F\,} P − 1 = E − 1 = [ 1 2 0 1 ] , Q = F − 1 = 1 0 0 2 1 0 3 0 1 , P = E , Q − 1 = F
最终验证 :
P − 1 Γ Q = [ 1 0 2 1 ] [ 1 0 0 0 0 0 ] [ 1 2 3 0 1 0 0 0 1 ] = [ 1 0 0 2 0 0 ] [ 1 2 3 0 1 0 0 0 1 ] = [ 1 2 3 2 4 6 ] = A ✓ P^{-1}\Gamma Q=\begin{bmatrix}1&0\\2&1\end{bmatrix}\begin{bmatrix}1&0&0\\0&0&0\end{bmatrix}\begin{bmatrix}1&2&3\\0&1&0\\0&0&1\end{bmatrix} =\begin{bmatrix}1&0&0\\2&0&0\end{bmatrix}\begin{bmatrix}1&2&3\\0&1&0\\0&0&1\end{bmatrix}=\begin{bmatrix}1&2&3\\2&4&6\end{bmatrix}=A\ \checkmark P − 1 Γ Q = [ 1 2 0 1 ] [ 1 0 0 0 0 0 ] 1 0 0 2 1 0 3 0 1 = [ 1 2 0 0 0 0 ] 1 0 0 2 1 0 3 0 1 = [ 1 2 2 4 3 6 ] = A ✓
要点 :读出的两块是「变换的总积」E E E 、F F F ,满足 E A F = Γ EAF=\Gamma E A F = Γ ;写分解式时要取逆 (P − 1 = E − 1 P^{-1}=E^{-1} P − 1 = E − 1 、Q = F − 1 Q=F^{-1} Q = F − 1 )。用 E A F = Γ EAF=\Gamma E A F = Γ 验算比用 P − 1 Γ Q P^{-1}\Gamma Q P − 1 Γ Q 少两步矩阵乘法,也不容易把栏目记错。∎
C2 · 用相抵标准形解方程组
用 C1 的分解解 A x = b A\boldsymbol{x}=\boldsymbol{b} A x = b (A = [ 1 2 3 2 4 6 ] A=\begin{bmatrix}1&2&3\\2&4&6\end{bmatrix} A = [ 1 2 2 4 3 6 ] ,b = ( 1 , 2 ) T \boldsymbol{b}=(1,2)^{\mathrm T} b = ( 1 , 2 ) T ),写出通解;再判断 b ′ = ( 1 , 3 ) T \boldsymbol{b}'=(1,3)^{\mathrm T} b ′ = ( 1 , 3 ) T 是否有解。
解答 · C2
取 E = [ 1 0 − 2 1 ] E=\begin{bmatrix}1&0\\-2&1\end{bmatrix} E = [ 1 − 2 0 1 ] (E A F = Γ EAF=\Gamma E A F = Γ )。方程 Γ y = E b \Gamma\boldsymbol{y}=E\boldsymbol{b} Γ y = E b ,x = F y \boldsymbol{x}=F\boldsymbol{y} x = F y (因为 A = E − 1 Γ F − 1 A=E^{-1}\Gamma F^{-1} A = E − 1 Γ F − 1 ,故 Γ ( F − 1 x ) = E b \Gamma(F^{-1}\boldsymbol{x})=E\boldsymbol{b} Γ ( F − 1 x ) = E b ,记 y = F − 1 x \boldsymbol{y}=F^{-1}\boldsymbol{x} y = F − 1 x )。
(a) E b = [ 1 0 − 2 1 ] [ 1 2 ] = [ 1 0 ] E\boldsymbol{b}=\begin{bmatrix}1&0\\-2&1\end{bmatrix}\begin{bmatrix}1\\2\end{bmatrix}=\begin{bmatrix}1\\0\end{bmatrix} E b = [ 1 − 2 0 1 ] [ 1 2 ] = [ 1 0 ] ,后 m − r = 1 m-r=1 m − r = 1 个分量为 0 0 0 ✓ 有解。由 Γ y = ( y 1 , 0 ) T \Gamma\boldsymbol{y}=(y_1,0)^{\mathrm T} Γ y = ( y 1 , 0 ) T 得 y 1 = 1 y_1=1 y 1 = 1 、y 2 , y 3 y_2,y_3 y 2 , y 3 任意(记 y 2 = s , y 3 = t y_2=s,y_3=t y 2 = s , y 3 = t ):
x = F y = [ 1 − 2 − 3 0 1 0 0 0 1 ] [ 1 s t ] = [ 1 − 2 s − 3 t s t ] . \boldsymbol{x}=F\boldsymbol{y}=\begin{bmatrix}1&-2&-3\\0&1&0\\0&0&1\end{bmatrix}\begin{bmatrix}1\\s\\t\end{bmatrix}=\begin{bmatrix}1-2s-3t\\ s\\ t\end{bmatrix}. x = F y = 1 0 0 − 2 1 0 − 3 0 1 1 s t = 1 − 2 s − 3 t s t .
验算 :A x = ( 1 − 2 s − 3 t + 2 s + 3 t , 2 − 4 s − 6 t + 4 s + 6 t ) T = ( 1 , 2 ) T = b A\boldsymbol{x}=(1-2s-3t+2s+3t,\ 2-4s-6t+4s+6t)^{\mathrm T}=(1,2)^{\mathrm T}=\boldsymbol{b} A x = ( 1 − 2 s − 3 t + 2 s + 3 t , 2 − 4 s − 6 t + 4 s + 6 t ) T = ( 1 , 2 ) T = b ✓;自由参数 2 = n − r 2=n-r 2 = n − r ✓。
(b) E b ′ = [ 1 − 2 + 3 ] = [ 1 1 ] E\boldsymbol{b}'=\begin{bmatrix}1\\-2+3\end{bmatrix}=\begin{bmatrix}1\\1\end{bmatrix} E b ′ = [ 1 − 2 + 3 ] = [ 1 1 ] ,第二个分量 ≠ 0 \neq0 = 0 ,无解 ✓(与经典判据一致:rank ( A ∣ b ′ ) = 2 > rank A = 1 \operatorname{rank}(A|\boldsymbol{b}')=2>\operatorname{rank}A=1 rank ( A ∣ b ′ ) = 2 > rank A = 1 ,因为 b ′ \boldsymbol{b}' b ′ 不满足「第二分量 = 2 × =2\times = 2 × 第一分量」)。
C3 · 相抵分解重做(例 3.1.A)
求 A = [ 1 2 3 2 5 7 3 7 10 ] A=\begin{bmatrix}1&2&3\\ 2&5&7\\ 3&7&10\end{bmatrix} A = 1 2 3 2 5 7 3 7 10 的相抵分解式(Γ \Gamma Γ 、E E E 、F F F ,并写出 A = P − 1 Γ Q A=P^{-1}\Gamma Q A = P − 1 Γ Q )。
解答 · C3
过程见 例 3.1.A (分行、列两步;r = 2 r=2 r = 2 ):
Γ = [ 1 0 0 0 1 0 0 0 0 ] , E = [ 1 0 0 − 2 1 0 − 1 − 1 1 ] , F = [ 1 − 2 − 1 0 1 − 1 0 0 1 ] , E A F = Γ ✓ \Gamma=\begin{bmatrix}1&0&0\\ 0&1&0\\ 0&0&0\end{bmatrix},\qquad E=\begin{bmatrix}1&0&0\\ -2&1&0\\ -1&-1&1\end{bmatrix},\qquad F=\begin{bmatrix}1&-2&-1\\ 0&1&-1\\ 0&0&1\end{bmatrix},\qquad EAF=\Gamma\ \checkmark Γ = 1 0 0 0 1 0 0 0 0 , E = 1 − 2 − 1 0 1 − 1 0 0 1 , F = 1 0 0 − 2 1 0 − 1 − 1 1 , E A F = Γ ✓
A = E − 1 Γ F − 1 , E − 1 = [ 1 0 0 2 1 0 3 1 1 ] , F − 1 = [ 1 2 3 0 1 1 0 0 1 ] . A=E^{-1}\Gamma F^{-1},\qquad E^{-1}=\begin{bmatrix}1&0&0\\ 2&1&0\\ 3&1&1\end{bmatrix},\qquad F^{-1}=\begin{bmatrix}1&2&3\\ 0&1&1\\ 0&0&1\end{bmatrix}. A = E − 1 Γ F − 1 , E − 1 = 1 2 3 0 1 1 0 0 1 , F − 1 = 1 0 0 2 1 0 3 1 1 .
验算 :E − 1 Γ F − 1 = [ 1 0 0 2 1 0 3 1 0 ] [ 1 2 3 0 1 1 0 0 1 ] = [ 1 2 3 2 5 7 3 7 10 ] = A E^{-1}\Gamma F^{-1}=\begin{bmatrix}1&0&0\\2&1&0\\3&1&0\end{bmatrix}\begin{bmatrix}1&2&3\\0&1&1\\0&0&1\end{bmatrix}=\begin{bmatrix}1&2&3\\2&5&7\\3&7&10\end{bmatrix}=A E − 1 Γ F − 1 = 1 2 3 0 1 1 0 0 0 1 0 0 2 1 0 3 1 1 = 1 2 3 2 5 7 3 7 10 = A ✓。∎
C4 · 全部右逆
求 A = [ 1 2 0 0 1 1 ] A=\begin{bmatrix}1&2&0\\ 0&1&1\end{bmatrix} A = [ 1 0 2 1 0 1 ] 的全部 右逆。
提示 · C4
先用公式 B 0 = A H ( A A H ) − 1 B_0=A^{\mathrm H}(AA^{\mathrm H})^{-1} B 0 = A H ( A A H ) − 1 造一个,再考虑「差」:A ( B − B 0 ) = O A(B-B_0)=O A ( B − B 0 ) = O ,即 B − B 0 B-B_0 B − B 0 的每一列都属于 N ( A ) N(A) N ( A ) ——把 N ( A ) N(A) N ( A ) 的基求出来后加上参数即可。
解答 · C4
基本解 :A A H = [ 5 2 2 2 ] AA^{\mathrm H}=\begin{bmatrix}5&2\\2&2\end{bmatrix} A A H = [ 5 2 2 2 ] ,( A A H ) − 1 = 1 6 [ 2 − 2 − 2 5 ] (AA^{\mathrm H})^{-1}=\dfrac16\begin{bmatrix}2&-2\\-2&5\end{bmatrix} ( A A H ) − 1 = 6 1 [ 2 − 2 − 2 5 ] ,
B 0 = A H ( A A H ) − 1 = 1 6 [ 2 − 2 2 1 − 2 5 ] , A B 0 = I 2 ✓ B_0=A^{\mathrm H}(AA^{\mathrm H})^{-1}=\dfrac16\begin{bmatrix}2&-2\\ 2&1\\ -2&5\end{bmatrix},\qquad AB_0=I_2\ \checkmark B 0 = A H ( A A H ) − 1 = 6 1 2 2 − 2 − 2 1 5 , A B 0 = I 2 ✓
通式 :设 B = B 0 + Z B=B_0+Z B = B 0 + Z ,则 A B = B 0 AB=B_0 A B = B 0 与 A = I A=I A = I 的关系为 A B = A ( B 0 + Z ) = I 2 + A Z AB=A(B_0+Z)=I_2+AZ A B = A ( B 0 + Z ) = I 2 + A Z ;要 A B = I 2 AB=I_2 A B = I 2 只需 A Z = O AZ=O A Z = O ,即 Z Z Z 的每一列 属于 N ( A ) N(A) N ( A ) 。解 A x = 0 A\boldsymbol{x}=\boldsymbol{0} A x = 0 :x 1 + 2 x 2 = 0 x_1+2x_2=0 x 1 + 2 x 2 = 0 、x 2 + x 3 = 0 x_2+x_3=0 x 2 + x 3 = 0 ,取 x 3 = 1 x_3=1 x 3 = 1 得 n = ( 2 , − 1 , 1 ) T \boldsymbol{n}=(2,-1,1)^{\mathrm T} n = ( 2 , − 1 , 1 ) T ,N ( A ) = span { n } N(A)=\operatorname{span}\{\boldsymbol{n}\} N ( A ) = span { n } (dim = n − r = 1 \dim=n-r=1 dim = n − r = 1 ✓)。
于是全部右逆为
B = 1 6 [ 2 − 2 2 1 − 2 5 ] + [ 2 s 2 t − s − t s t ] , s , t ∈ C B=\dfrac16\begin{bmatrix}2&-2\\ 2&1\\ -2&5\end{bmatrix}+\begin{bmatrix}2s&2t\\ -s&-t\\ s&t\end{bmatrix},\qquad s,t\in\mathbb C B = 6 1 2 2 − 2 − 2 1 5 + 2 s − s s 2 t − t t , s , t ∈ C
(每个参数配一列)。验算 :A [ 2 s 2 t − s − t s t ] = [ 2 s − 2 s 2 t − 2 t − s + s − t + t ] = O A\begin{bmatrix}2s&2t\\-s&-t\\ s&t\end{bmatrix}=\begin{bmatrix}2s-2s&2t-2t\\ -s+s&-t+t\end{bmatrix}=O A 2 s − s s 2 t − t t = [ 2 s − 2 s − s + s 2 t − 2 t − t + t ] = O ✓,故 A B = I 2 AB=I_2 A B = I 2 ✓。两个自由参数 = = = n − r n-r n − r (每列一个)× \times × m m m (m = 2 m=2 m = 2 列)÷ \div ÷ ⋯——不背计数公式,只记住「多余的自由度来自 N ( A ) N(A) N ( A ) 」。
C5 · 判断并构造单侧逆
判断下列矩阵是否有右逆 / 左逆;对存在的,用公式构造一个:
(a) A 1 = [ 1 0 0 0 1 0 ] A_1=\begin{bmatrix}1&0&0\\ 0&1&0\end{bmatrix} A 1 = [ 1 0 0 1 0 0 ] ;(b) A 2 = [ 1 0 0 1 0 0 ] A_2=\begin{bmatrix}1&0\\ 0&1\\ 0&0\end{bmatrix} A 2 = 1 0 0 0 1 0 ;(c) A 3 = [ 1 2 2 4 ] A_3=\begin{bmatrix}1&2\\ 2&4\end{bmatrix} A 3 = [ 1 2 2 4 ] 。
解答 · C5
(a) 行满秩(r = 2 = m r=2=m r = 2 = m ):有右逆。A 1 A 1 H = I 2 A_1A_1^{\mathrm H}=I_2 A 1 A 1 H = I 2 (A 1 A_1 A 1 的行已标准正交),故 B = A 1 H ( A 1 A 1 H ) − 1 = A 1 H = [ 1 0 0 1 0 0 ] B=A_1^{\mathrm H}(A_1A_1^{\mathrm H})^{-1}=A_1^{\mathrm H}=\begin{bmatrix}1&0\\0&1\\0&0\end{bmatrix} B = A 1 H ( A 1 A 1 H ) − 1 = A 1 H = 1 0 0 0 1 0 ;验算 A 1 B = I 2 A_1B=I_2 A 1 B = I 2 ✓。
(b) 列满秩(r = 2 = n r=2=n r = 2 = n ):有左逆。A 2 H A 2 = I 2 A_2^{\mathrm H}A_2=I_2 A 2 H A 2 = I 2 ,故 B = ( A 2 H A 2 ) − 1 A 2 H = A 2 H = [ 1 0 0 ; 0 1 0 ] B=(A_2^{\mathrm H}A_2)^{-1}A_2^{\mathrm H}=A_2^{\mathrm H}=[1\ 0\ 0;\ 0\ 1\ 0] B = ( A 2 H A 2 ) − 1 A 2 H = A 2 H = [ 1 0 0 ; 0 1 0 ] ,即 [ 1 0 0 0 1 0 ] \begin{bmatrix}1&0&0\\0&1&0\end{bmatrix} [ 1 0 0 1 0 0 ] ;验算 B A 2 = I 2 BA_2=I_2 B A 2 = I 2 ✓。
(c) 秩 1 < 2 1<2 1 < 2 :两者都没有(若硬算会发现 A 3 H A 3 = [ 5 10 10 20 ] A_3^{\mathrm H}A_3=\begin{bmatrix}5&10\\10&20\end{bmatrix} A 3 H A 3 = [ 5 10 10 20 ] 不可逆,公式失效——这与 定理 3.1.1 的判断一致)。
要点 :(a)(b) 互为「转置镜像」,两个公式也互为镜像。
C6 · 右逆公式实操
设 A = [ 1 1 0 0 1 1 ] A=\begin{bmatrix}1&1&0\\ 0&1&1\end{bmatrix} A = [ 1 0 1 1 0 1 ] :判断 A A A 是否有右逆;用公式构造一个右逆并验算;说明为什么不用相抵分解也能得到答案。
解答 · C6
判断 :2 × 3 2\times3 2 × 3 、两行不成比例 ⇒ \Rightarrow ⇒ rank A = 2 = m \operatorname{rank}A=2=m rank A = 2 = m ,行满秩,有右逆 。
构造 :
A A H = [ 1 1 0 0 1 1 ] [ 1 0 1 1 0 1 ] = [ 2 1 1 2 ] , ( A A H ) − 1 = 1 3 [ 2 − 1 − 1 2 ] , AA^{\mathrm H}=\begin{bmatrix}1&1&0\\0&1&1\end{bmatrix}\begin{bmatrix}1&0\\1&1\\0&1\end{bmatrix}=\begin{bmatrix}2&1\\1&2\end{bmatrix},\qquad (AA^{\mathrm H})^{-1}=\frac13\begin{bmatrix}2&-1\\-1&2\end{bmatrix}, A A H = [ 1 0 1 1 0 1 ] 1 1 0 0 1 1 = [ 2 1 1 2 ] , ( A A H ) − 1 = 3 1 [ 2 − 1 − 1 2 ] ,
B = A H ( A A H ) − 1 = 1 3 [ 1 0 1 1 0 1 ] [ 2 − 1 − 1 2 ] = 1 3 [ 2 − 1 1 1 − 1 2 ] . B=A^{\mathrm H}(AA^{\mathrm H})^{-1}=\frac13\begin{bmatrix}1&0\\1&1\\0&1\end{bmatrix}\begin{bmatrix}2&-1\\-1&2\end{bmatrix}=\frac13\begin{bmatrix}2&-1\\1&1\\-1&2\end{bmatrix}. B = A H ( A A H ) − 1 = 3 1 1 1 0 0 1 1 [ 2 − 1 − 1 2 ] = 3 1 2 1 − 1 − 1 1 2 .
验算 :A B = 1 3 [ 1 1 0 0 1 1 ] [ 2 − 1 1 1 − 1 2 ] = 1 3 [ 3 0 0 3 ] = I 2 AB=\dfrac13\begin{bmatrix}1&1&0\\0&1&1\end{bmatrix}\begin{bmatrix}2&-1\\1&1\\-1&2\end{bmatrix}=\dfrac13\begin{bmatrix}3&0\\0&3\end{bmatrix}=I_2 A B = 3 1 [ 1 0 1 1 0 1 ] 2 1 − 1 − 1 1 2 = 3 1 [ 3 0 0 3 ] = I 2 ✓。
为什么不用相抵分解 :法 2 的公式依赖 命题 3.1.1 (rank ( A A H ) = rank A = m \operatorname{rank}(AA^{\mathrm H})=\operatorname{rank}A=m rank ( A A H ) = rank A = m ,故 A A H AA^{\mathrm H} A A H 可逆),而这一命题的证明只用方程组同解 + 内积正定性 ,与相抵分解无关——所以公式法比相抵法短得多。
C7 · 满秩分解(取 A A A 的列)
求 A = [ 1 2 3 0 1 1 ] A=\begin{bmatrix}1&2&3\\ 0&1&1\end{bmatrix} A = [ 1 0 2 1 3 1 ] 的满秩分解,并写出 R ( A ) R(A) R ( A ) 的一组基。
解答 · C7
定基 :a 1 = ( 1 , 0 ) T \boldsymbol{a}_1=(1,0)^{\mathrm T} a 1 = ( 1 , 0 ) T 、a 2 = ( 2 , 1 ) T \boldsymbol{a}_2=(2,1)^{\mathrm T} a 2 = ( 2 , 1 ) T 、a 3 = ( 3 , 1 ) T \boldsymbol{a}_3=(3,1)^{\mathrm T} a 3 = ( 3 , 1 ) T 。a 1 , a 2 \boldsymbol{a}_1,\boldsymbol{a}_2 a 1 , a 2 无关;a 3 = a 1 + a 2 \boldsymbol{a}_3=\boldsymbol{a}_1+\boldsymbol{a}_2 a 3 = a 1 + a 2 ✓。r = rank A = 2 r=\operatorname{rank}A=2 r = rank A = 2 ,故 a 1 , a 2 \boldsymbol{a}_1,\boldsymbol{a}_2 a 1 , a 2 是 R ( A ) R(A) R ( A ) 的基。
满秩分解 :
B = [ a 1 , a 2 ] = [ 1 2 0 1 ] , C = [ 1 0 1 0 1 1 ] ( 坐标: a 3 = a 1 + a 2 ) , B=[\boldsymbol{a}_1,\boldsymbol{a}_2]=\begin{bmatrix}1&2\\ 0&1\end{bmatrix},\qquad C=\begin{bmatrix}1&0&1\\ 0&1&1\end{bmatrix}\quad(\text{坐标:}\boldsymbol{a}_3=\boldsymbol{a}_1+\boldsymbol{a}_2), B = [ a 1 , a 2 ] = [ 1 0 2 1 ] , C = [ 1 0 0 1 1 1 ] ( 坐标: a 3 = a 1 + a 2 ) ,
B C = [ 1 2 0 1 ] [ 1 0 1 0 1 1 ] = [ 1 2 3 0 1 1 ] = A ✓ BC=\begin{bmatrix}1&2\\ 0&1\end{bmatrix}\begin{bmatrix}1&0&1\\ 0&1&1\end{bmatrix}=\begin{bmatrix}1&2&3\\ 0&1&1\end{bmatrix}=A\ \checkmark BC = [ 1 0 2 1 ] [ 1 0 0 1 1 1 ] = [ 1 0 2 1 3 1 ] = A ✓
R ( A ) = span { ( 1 , 0 ) T , ( 2 , 1 ) T } R(A)=\operatorname{span}\{(1,0)^{\mathrm T},(2,1)^{\mathrm T}\} R ( A ) = span {( 1 , 0 ) T , ( 2 , 1 ) T } 。∎
C8 · 满秩分解(RREF 法)
用 用 RREF 批量求满秩分解 求 A = [ 1 1 2 1 2 3 5 2 1 2 3 1 ] A=\begin{bmatrix}1&1&2&1\\ 2&3&5&2\\ 1&2&3&1\end{bmatrix} A = 1 2 1 1 3 2 2 5 3 1 2 1 的满秩分解,并读出 A A A 的列之间的线性关系。
解答 · C8
化 RREF :r 2 − 2 r 1 = ( 0 , 1 , 1 , 0 ) r_2-2r_1=(0,1,1,0) r 2 − 2 r 1 = ( 0 , 1 , 1 , 0 ) ;r 3 − r 2 = ( 0 , 0 , 0 , 0 ) r_3-r_2=(0,0,0,0) r 3 − r 2 = ( 0 , 0 , 0 , 0 ) ;回代 r 1 − r 2 = ( 1 , 0 , 1 , 1 ) r_1-r_2=(1,0,1,1) r 1 − r 2 = ( 1 , 0 , 1 , 1 ) 。得
R = [ 1 0 1 1 0 1 1 0 0 0 0 0 ] , r = 2 , 主元列 = 1 , 2. R=\begin{bmatrix}1&0&1&1\\ 0&1&1&0\\ 0&0&0&0\end{bmatrix},\qquad r=2,\ \text{主元列}=1,2 . R = 1 0 0 0 1 0 1 1 0 1 0 0 , r = 2 , 主元列 = 1 , 2.
取 B B B 、C C C :B = A B=A B = A 的第 1、2 列 = [ 1 1 2 3 1 2 ] =\begin{bmatrix}1&1\\2&3\\1&2\end{bmatrix} = 1 2 1 1 3 2 ,C = R C=R C = R 的前两行 = [ 1 0 1 1 0 1 1 0 ] =\begin{bmatrix}1&0&1&1\\ 0&1&1&0\end{bmatrix} = [ 1 0 0 1 1 1 1 0 ] 。
验算 :[ 1 1 2 3 1 2 ] [ 1 0 1 1 0 1 1 0 ] = [ 1 1 2 1 2 3 5 2 1 2 3 1 ] = A \begin{bmatrix}1&1\\2&3\\1&2\end{bmatrix}\begin{bmatrix}1&0&1&1\\0&1&1&0\end{bmatrix}=\begin{bmatrix}1&1&2&1\\2&3&5&2\\1&2&3&1\end{bmatrix}=A 1 2 1 1 3 2 [ 1 0 0 1 1 1 1 0 ] = 1 2 1 1 3 2 2 5 3 1 2 1 = A ✓。
列关系 :C C C 的第三列 ( 1 , 1 ) T (1,1)^{\mathrm T} ( 1 , 1 ) T ⟹ a 3 = a 1 + a 2 \boldsymbol{a}_3=\boldsymbol{a}_1+\boldsymbol{a}_2 a 3 = a 1 + a 2 ;第四列 ( 1 , 0 ) T (1,0)^{\mathrm T} ( 1 , 0 ) T ⟹ a 4 = a 1 \boldsymbol{a}_4=\boldsymbol{a}_1 a 4 = a 1 。代入验算:a 3 = ( 2 , 5 , 3 ) \boldsymbol{a}_3=(2,5,3) a 3 = ( 2 , 5 , 3 ) 而 a 1 + a 2 = ( 1 , 1 ) + ( 1 , 2 ) = ( 2 , 3 , 1 ) \boldsymbol{a}_1+\boldsymbol{a}_2=(1,1)+(1,2)=(2,3,1) a 1 + a 2 = ( 1 , 1 ) + ( 1 , 2 ) = ( 2 , 3 , 1 ) ——不对 !查 A A A :a 3 = ( 2 , 5 , 3 ) T \boldsymbol{a}_3=(2,5,3)^{\mathrm T} a 3 = ( 2 , 5 , 3 ) T ,a 1 + a 2 = ( 1 , 1 , 1 ) + ( 2 , 3 , 2 ) = ( 3 , 4 , 3 ) \boldsymbol{a}_1+\boldsymbol{a}_2=(1,1,1)+(2,3,2)=(3,4,3) a 1 + a 2 = ( 1 , 1 , 1 ) + ( 2 , 3 , 2 ) = ( 3 , 4 , 3 ) ,都不等于。重新核对 R R R 的行 r 1 = ( 1 , 0 , 1 , 1 ) r_1=(1,0,1,1) r 1 = ( 1 , 0 , 1 , 1 ) 、r 2 = ( 0 , 1 , 1 , 0 ) r_2=(0,1,1,0) r 2 = ( 0 , 1 , 1 , 0 ) ,且 B B B 取的是 A A A 的第 1、2 列 :a 1 = ( 1 , 2 , 1 ) \boldsymbol{a}_1=(1,2,1) a 1 = ( 1 , 2 , 1 ) 、a 2 = ( 1 , 3 , 2 ) \boldsymbol{a}_2=(1,3,2) a 2 = ( 1 , 3 , 2 ) 。则列关系应为 a 3 = 1 ⋅ a 1 + 1 ⋅ a 2 = ( 2 , 5 , 3 ) \boldsymbol{a}_3=1\cdot\boldsymbol{a}_1+1\cdot\boldsymbol{a}_2=(2,5,3) a 3 = 1 ⋅ a 1 + 1 ⋅ a 2 = ( 2 , 5 , 3 ) ✓(正是 A A A 的第三列);a 4 = 1 ⋅ a 1 + 0 ⋅ a 2 = ( 1 , 2 , 1 ) \boldsymbol{a}_4=1\cdot\boldsymbol{a}_1+0\cdot\boldsymbol{a}_2=(1,2,1) a 4 = 1 ⋅ a 1 + 0 ⋅ a 2 = ( 1 , 2 , 1 ) ✓。上面「不对」的检查把 B B B 的列抄着用了 A A A 的行,属笔误——教训:B B B 是 A A A 的列拼成的矩阵,核对时要用列向量 。∎
C9 · 两种满秩分解与过渡矩阵
对 A = [ 1 2 3 0 1 1 ] A=\begin{bmatrix}1&2&3\\ 0&1&1\end{bmatrix} A = [ 1 0 2 1 3 1 ] :写出两种不同的满秩分解 A = B 1 C 1 = B 2 C 2 A=B_1C_1=B_2C_2 A = B 1 C 1 = B 2 C 2 ,求出 定理 3.2.2 中的 D D D ,并验证 B 1 = B 2 D B_1=B_2D B 1 = B 2 D 、C 1 = D − 1 C 2 C_1=D^{-1}C_2 C 1 = D − 1 C 2 。
解答 · C9
分解一 (取第 1、2 列):B 1 = [ 1 2 0 1 ] B_1=\begin{bmatrix}1&2\\ 0&1\end{bmatrix} B 1 = [ 1 0 2 1 ] ,C 1 = [ 1 0 1 0 1 1 ] C_1=\begin{bmatrix}1&0&1\\ 0&1&1\end{bmatrix} C 1 = [ 1 0 0 1 1 1 ] (B 1 B_1 B 1 可逆,C 1 = B 1 − 1 A C_1=B_1^{-1}A C 1 = B 1 − 1 A )。
分解二 (取第 1、3 列):B 2 = [ 1 3 0 1 ] B_2=\begin{bmatrix}1&3\\ 0&1\end{bmatrix} B 2 = [ 1 0 3 1 ] ,C 2 = [ 1 − 1 0 0 1 1 ] C_2=\begin{bmatrix}1&-1&0\\ 0&1&1\end{bmatrix} C 2 = [ 1 0 − 1 1 0 1 ] 。验算:[ 1 3 0 1 ] [ 1 − 1 0 0 1 1 ] = [ 1 2 3 0 1 1 ] = A \begin{bmatrix}1&3\\0&1\end{bmatrix}\begin{bmatrix}1&-1&0\\0&1&1\end{bmatrix}=\begin{bmatrix}1&2&3\\0&1&1\end{bmatrix}=A [ 1 0 3 1 ] [ 1 0 − 1 1 0 1 ] = [ 1 0 2 1 3 1 ] = A ✓。
求 D D D :由 B 1 = B 2 D B_1=B_2D B 1 = B 2 D 得
D = B 2 − 1 B 1 = [ 1 − 3 0 1 ] [ 1 2 0 1 ] = [ 1 − 1 0 1 ] , det D = 1 ≠ 0 ✓ D=B_2^{-1}B_1=\begin{bmatrix}1&-3\\ 0&1\end{bmatrix}\begin{bmatrix}1&2\\ 0&1\end{bmatrix}=\begin{bmatrix}1&-1\\ 0&1\end{bmatrix},\qquad \det D=1\neq0\ \checkmark D = B 2 − 1 B 1 = [ 1 0 − 3 1 ] [ 1 0 2 1 ] = [ 1 0 − 1 1 ] , det D = 1 = 0 ✓
验算 C 1 = D − 1 C 2 C_1=D^{-1}C_2 C 1 = D − 1 C 2 :D − 1 = [ 1 1 0 1 ] D^{-1}=\begin{bmatrix}1&1\\ 0&1\end{bmatrix} D − 1 = [ 1 0 1 1 ] ,
D − 1 C 2 = [ 1 1 0 1 ] [ 1 − 1 0 0 1 1 ] = [ 1 0 1 0 1 1 ] = C 1 ✓ D^{-1}C_2=\begin{bmatrix}1&1\\ 0&1\end{bmatrix}\begin{bmatrix}1&-1&0\\ 0&1&1\end{bmatrix}=\begin{bmatrix}1&0&1\\ 0&1&1\end{bmatrix}=C_1\ \checkmark D − 1 C 2 = [ 1 0 1 1 ] [ 1 0 − 1 1 0 1 ] = [ 1 0 0 1 1 1 ] = C 1 ✓
读 D D D :D D D 是「由 B 2 B_2 B 2 的列到 B 1 B_1 B 1 的列」的过渡矩阵——即 ( [ a 1 , a 3 ] ) = ( [ a 1 , a 2 ] ) D ([\boldsymbol{a}_1,\boldsymbol{a}_3])=([\boldsymbol{a}_1,\boldsymbol{a}_2])D ([ a 1 , a 3 ]) = ([ a 1 , a 2 ]) D ,与 注 3 一致。∎
C10 · 相抵分解(3 × 2 3\times2 3 × 2 、秩 1 1 1 )
求 A = [ 1 − 1 0 0 2 − 2 ] A=\begin{bmatrix}1&-1\\ 0&0\\ 2&-2\end{bmatrix} A = 1 0 2 − 1 0 − 2 的相抵分解式。
解答 · C10
秩 :r = 1 r=1 r = 1 (三行都与 ( 1 , − 1 ) (1,-1) ( 1 , − 1 ) 成比例;实际 r 3 = 2 r 1 \boldsymbol{r}_3=2\boldsymbol{r}_1 r 3 = 2 r 1 、r 2 = 0 \boldsymbol{r}_2=\boldsymbol{0} r 2 = 0 )。
消元 :r 3 ← r 3 − 2 r 1 r_3\leftarrow r_3-2r_1 r 3 ← r 3 − 2 r 1 ;c 2 ← c 2 + c 1 c_2\leftarrow c_2+c_1 c 2 ← c 2 + c 1 。得 E = [ 1 0 0 0 1 0 − 2 0 1 ] E=\begin{bmatrix}1&0&0\\ 0&1&0\\ -2&0&1\end{bmatrix} E = 1 0 − 2 0 1 0 0 0 1 、F = [ 1 1 0 1 ] F=\begin{bmatrix}1&1\\ 0&1\end{bmatrix} F = [ 1 0 1 1 ] ,E A F = [ 1 0 0 0 0 0 ] = Γ EAF=\begin{bmatrix}1&0\\0&0\\0&0\end{bmatrix}=\Gamma E A F = 1 0 0 0 0 0 = Γ ✓。
分解式 :A = E − 1 Γ F − 1 A=E^{-1}\Gamma F^{-1} A = E − 1 Γ F − 1 ,
E − 1 = [ 1 0 0 0 1 0 2 0 1 ] , F − 1 = [ 1 − 1 0 1 ] , E − 1 Γ F − 1 = [ 1 − 1 0 0 2 − 2 ] = A ✓ E^{-1}=\begin{bmatrix}1&0&0\\ 0&1&0\\ 2&0&1\end{bmatrix},\qquad F^{-1}=\begin{bmatrix}1&-1\\ 0&1\end{bmatrix},\qquad E^{-1}\Gamma F^{-1}=\begin{bmatrix}1&-1\\ 0&0\\ 2&-2\end{bmatrix}=A\ \checkmark E − 1 = 1 0 2 0 1 0 0 0 1 , F − 1 = [ 1 0 − 1 1 ] , E − 1 Γ F − 1 = 1 0 2 − 1 0 − 2 = A ✓
∎
C11 · 四个秩的核对
对 A = [ 1 2 2 4 3 6 ] A=\begin{bmatrix}1&2\\ 2&4\\ 3&6\end{bmatrix} A = 1 2 3 2 4 6 :计算 A H A A^{\mathrm H}A A H A 、A A H AA^{\mathrm H} A A H ,并核对 命题 3.1.1 的四个秩。
解答 · C11
A H A = [ 1 2 3 2 4 6 ] [ 1 2 2 4 3 6 ] = [ 14 28 28 56 ] , A A H = [ 1 2 2 4 3 6 ] [ 1 2 3 2 4 6 ] = [ 5 10 15 10 20 30 15 30 45 ] . A^{\mathrm H}A=\begin{bmatrix}1&2&3\\2&4&6\end{bmatrix}\begin{bmatrix}1&2\\2&4\\3&6\end{bmatrix}=\begin{bmatrix}14&28\\ 28&56\end{bmatrix},\qquad AA^{\mathrm H}=\begin{bmatrix}1&2\\2&4\\3&6\end{bmatrix}\begin{bmatrix}1&2&3\\2&4&6\end{bmatrix}=\begin{bmatrix}5&10&15\\ 10&20&30\\ 15&30&45\end{bmatrix}. A H A = [ 1 2 2 4 3 6 ] 1 2 3 2 4 6 = [ 14 28 28 56 ] , A A H = 1 2 3 2 4 6 [ 1 2 2 4 3 6 ] = 5 10 15 10 20 30 15 30 45 .
秩 :rank A = 1 \operatorname{rank}A=1 rank A = 1 (列成比例);A H A A^{\mathrm H}A A H A 第二行 = 2 × =2\times = 2 × 第一行 ⇒ \Rightarrow ⇒ 秩 1 1 1 ;A A H AA^{\mathrm H} A A H 的各行成比例 ⇒ \Rightarrow ⇒ 秩 1 1 1 ;rank ( A H ) = 1 \operatorname{rank}(A^{\mathrm H})=1 rank ( A H ) = 1 。四项相等 ✓:
rank ( A ) = rank ( A H A ) = rank ( A A H ) = rank ( A H ) = 1. \operatorname{rank}(A)=\operatorname{rank}(A^{\mathrm H}A)=\operatorname{rank}(AA^{\mathrm H})=\operatorname{rank}(A^{\mathrm H})=1 . rank ( A ) = rank ( A H A ) = rank ( A A H ) = rank ( A H ) = 1.
要点 :A H A A^{\mathrm H}A A H A 、A A H AA^{\mathrm H} A A H 与 A A A 同秩,但尺寸不同 (A H A A^{\mathrm H}A A H A 是 2 × 2 2\times2 2 × 2 、A A H AA^{\mathrm H} A A H 是 3 × 3 3\times3 3 × 3 )——秩相同、零空间不同(N ( A H A ) = N ( A ) N(A^{\mathrm H}A)=N(A) N ( A H A ) = N ( A ) ,而 N ( A A H ) = N ( A H ) N(AA^{\mathrm H})=N(A^{\mathrm H}) N ( A A H ) = N ( A H ) )。∎
C12 · 满秩分解与四个子空间
对 A = [ 1 0 1 0 1 1 1 1 2 ] A=\begin{bmatrix}1&0&1\\ 0&1&1\\ 1&1&2\end{bmatrix} A = 1 0 1 0 1 1 1 1 2 :求满秩分解;由此写出 R ( A ) R(A) R ( A ) 的基、rank A \operatorname{rank}A rank A 与 N ( A ) N(A) N ( A ) 的基。
解答 · C12
满秩分解 :a 3 = a 1 + a 2 \boldsymbol{a}_3=\boldsymbol{a}_1+\boldsymbol{a}_2 a 3 = a 1 + a 2 ,r = 2 r=2 r = 2 ,取第 1、2 列:
B = [ 1 0 0 1 1 1 ] , C = [ 1 0 1 0 1 1 ] , B C = [ 1 0 1 0 1 1 1 1 2 ] = A ✓ B=\begin{bmatrix}1&0\\ 0&1\\ 1&1\end{bmatrix},\qquad C=\begin{bmatrix}1&0&1\\ 0&1&1\end{bmatrix},\qquad BC=\begin{bmatrix}1&0&1\\0&1&1\\1&1&2\end{bmatrix}=A\ \checkmark B = 1 0 1 0 1 1 , C = [ 1 0 0 1 1 1 ] , BC = 1 0 1 0 1 1 1 1 2 = A ✓
R ( A ) R(A) R ( A ) 的基 :B B B 的两列,R ( A ) = span { ( 1 , 0 , 1 ) T , ( 0 , 1 , 1 ) T } R(A)=\operatorname{span}\{(1,0,1)^{\mathrm T},(0,1,1)^{\mathrm T}\} R ( A ) = span {( 1 , 0 , 1 ) T , ( 0 , 1 , 1 ) T } ,rank A = 2 \operatorname{rank}A=2 rank A = 2 。
N ( A ) N(A) N ( A ) 的基 :解 A x = 0 A\boldsymbol{x}=\boldsymbol{0} A x = 0 :x 1 + x 3 = 0 x_1+x_3=0 x 1 + x 3 = 0 、x 2 + x 3 = 0 x_2+x_3=0 x 2 + x 3 = 0 ,取 x 3 = 1 x_3=1 x 3 = 1 :n = ( − 1 , − 1 , 1 ) T \boldsymbol{n}=(-1,-1,1)^{\mathrm T} n = ( − 1 , − 1 , 1 ) T 。验算 A n = ( − 1 + 1 , − 1 + 1 , − 1 − 1 + 2 ) T = 0 A\boldsymbol{n}=(-1+1,-1+1,-1-1+2)^{\mathrm T}=\boldsymbol{0} A n = ( − 1 + 1 , − 1 + 1 , − 1 − 1 + 2 ) T = 0 ✓。N ( A ) = span { n } N(A)=\operatorname{span}\{\boldsymbol{n}\} N ( A ) = span { n } ,dim = 1 = n − r \dim=1=n-r dim = 1 = n − r ✓。
一句话 :满秩分解给出「活的」部分(R ( A ) R(A) R ( A ) 的基在 B B B 里),而「死的」部分(N ( A ) N(A) N ( A ) )要另解方程——但两者的维数由 r r r 一分为二:2 + 1 = 3 2+1=3 2 + 1 = 3 ✓。
C13 · 右逆与全部右逆(2 × 3 2\times3 2 × 3 )
设 A = [ 1 2 3 0 1 1 ] A=\begin{bmatrix}1&2&3\\ 0&1&1\end{bmatrix} A = [ 1 0 2 1 3 1 ] :用公式构造一个右逆,并写出全部右逆。
提示 · C13
先算 A A H AA^{\mathrm H} A A H 与其逆;「全部右逆」的差集是「每列属于 N ( A ) N(A) N ( A ) 」的矩阵——N ( A ) N(A) N ( A ) 的基是 ( − 1 , − 1 , 1 ) T (-1,-1,1)^{\mathrm T} ( − 1 , − 1 , 1 ) T 。
解答 · C13
基本解 :A A H = [ 14 5 5 2 ] AA^{\mathrm H}=\begin{bmatrix}14&5\\5&2\end{bmatrix} A A H = [ 14 5 5 2 ] ,det = 3 \det=3 det = 3 ,( A A H ) − 1 = 1 3 [ 2 − 5 − 5 14 ] (AA^{\mathrm H})^{-1}=\dfrac13\begin{bmatrix}2&-5\\-5&14\end{bmatrix} ( A A H ) − 1 = 3 1 [ 2 − 5 − 5 14 ] ,
B 0 = A H ( A A H ) − 1 = 1 3 [ 1 0 2 1 3 1 ] [ 2 − 5 − 5 14 ] = 1 3 [ 2 − 5 − 1 4 1 − 1 ] , A B 0 = I 2 ✓ B_0=A^{\mathrm H}(AA^{\mathrm H})^{-1}=\frac13\begin{bmatrix}1&0\\2&1\\3&1\end{bmatrix}\begin{bmatrix}2&-5\\-5&14\end{bmatrix}=\frac13\begin{bmatrix}2&-5\\ -1&4\\ 1&-1\end{bmatrix},\qquad AB_0=I_2\ \checkmark B 0 = A H ( A A H ) − 1 = 3 1 1 2 3 0 1 1 [ 2 − 5 − 5 14 ] = 3 1 2 − 1 1 − 5 4 − 1 , A B 0 = I 2 ✓
N ( A ) N(A) N ( A ) :x 1 + 2 x 2 + 3 x 3 = 0 x_1+2x_2+3x_3=0 x 1 + 2 x 2 + 3 x 3 = 0 、x 2 + x 3 = 0 x_2+x_3=0 x 2 + x 3 = 0 得 n = ( − 1 , − 1 , 1 ) T \boldsymbol{n}=(-1,-1,1)^{\mathrm T} n = ( − 1 , − 1 , 1 ) T (验算 A n = ( − 1 − 2 + 3 , − 1 + 1 ) T = 0 A\boldsymbol{n}=(-1-2+3,-1+1)^{\mathrm T}=\boldsymbol{0} A n = ( − 1 − 2 + 3 , − 1 + 1 ) T = 0 ✓)。
全部右逆 :
B = B 0 + [ − s − t − s − t s t ] , s , t ∈ C B=B_0+\begin{bmatrix}-s&-t\\ -s&-t\\ s&t\end{bmatrix},\qquad s,t\in\mathbb C B = B 0 + − s − s s − t − t t , s , t ∈ C
验算 :A [ − s − t − s − t s t ] = [ − s − 2 s + 3 s − t − 2 t + 3 t − s + s − t + t ] = O A\begin{bmatrix}-s&-t\\-s&-t\\s&t\end{bmatrix}=\begin{bmatrix}-s-2s+3s&-t-2t+3t\\ -s+s&-t+t\end{bmatrix}=O A − s − s s − t − t t = [ − s − 2 s + 3 s − s + s − t − 2 t + 3 t − t + t ] = O ✓。∎
C14 · A X A = A AXA=A A X A = A 的一个解
设 A = [ 1 1 1 1 ] A=\begin{bmatrix}1&1\\ 1&1\end{bmatrix} A = [ 1 1 1 1 ] :用满秩分解求出 A X A = A AXA=A A X A = A 的一个解 X X X 并验证。
提示 · C14
r = 1 r=1 r = 1 :取 B = [ 1 1 ] B=\begin{bmatrix}1\\1\end{bmatrix} B = [ 1 1 ] 、C = [ 1 1 ] C=[1\ \ 1] C = [ 1 1 ] ,代 X = C H ( C C H ) − 1 ( B H B ) − 1 B H X=C^{\mathrm H}(CC^{\mathrm H})^{-1}(B^{\mathrm H}B)^{-1}B^{\mathrm H} X = C H ( C C H ) − 1 ( B H B ) − 1 B H 再验算 A X A = A AXA=A A X A = A 。
解答 · C14
满秩分解 :B = [ 1 1 ] B=\begin{bmatrix}1\\1\end{bmatrix} B = [ 1 1 ] 、C = [ 1 1 ] C=[1\ \ 1] C = [ 1 1 ] ,B C = A BC=A BC = A ✓。
构造 X X X :
X = C H ( C C H ) − 1 ( B H B ) − 1 B H = [ 1 1 ] ⋅ 1 2 ⋅ 1 2 ⋅ [ 1 1 ] = 1 4 [ 1 1 1 1 ] . X=C^{\mathrm H}(CC^{\mathrm H})^{-1}(B^{\mathrm H}B)^{-1}B^{\mathrm H} =\begin{bmatrix}1\\1\end{bmatrix}\cdot\frac12\cdot\frac12\cdot[1\ \ 1] =\frac14\begin{bmatrix}1&1\\ 1&1\end{bmatrix}. X = C H ( C C H ) − 1 ( B H B ) − 1 B H = [ 1 1 ] ⋅ 2 1 ⋅ 2 1 ⋅ [ 1 1 ] = 4 1 [ 1 1 1 1 ] .
验算 :
A X A = 1 4 [ 1 1 1 1 ] [ 1 1 1 1 ] [ 1 1 1 1 ] = 1 4 [ 2 2 2 2 ] [ 1 1 1 1 ] = 1 4 [ 4 4 4 4 ] = [ 1 1 1 1 ] = A ✓ AXA=\frac14\begin{bmatrix}1&1\\1&1\end{bmatrix}\begin{bmatrix}1&1\\1&1\end{bmatrix}\begin{bmatrix}1&1\\1&1\end{bmatrix} =\frac14\begin{bmatrix}2&2\\2&2\end{bmatrix}\begin{bmatrix}1&1\\1&1\end{bmatrix}=\frac14\begin{bmatrix}4&4\\4&4\end{bmatrix}=\begin{bmatrix}1&1\\1&1\end{bmatrix}=A\ \checkmark A X A = 4 1 [ 1 1 1 1 ] [ 1 1 1 1 ] [ 1 1 1 1 ] = 4 1 [ 2 2 2 2 ] [ 1 1 1 1 ] = 4 1 [ 4 4 4 4 ] = [ 1 1 1 1 ] = A ✓
补充核对 (用 本编补注 ):C X B = [ 1 1 ] ⋅ 1 4 [ 1 1 1 1 ] ⋅ [ 1 1 ] = [ 1 1 ] ⋅ 1 4 [ 2 2 ] = 1 = I 1 CXB=[1\ \ 1]\cdot\frac14\begin{bmatrix}1&1\\1&1\end{bmatrix}\cdot\begin{bmatrix}1\\1\end{bmatrix}=[1\ \ 1]\cdot\frac14\begin{bmatrix}2\\2\end{bmatrix}=1=I_1 CXB = [ 1 1 ] ⋅ 4 1 [ 1 1 1 1 ] ⋅ [ 1 1 ] = [ 1 1 ] ⋅ 4 1 [ 2 2 ] = 1 = I 1 ✓。
观察 :解不唯一 ——所有元素之和为 1 1 1 的 2 × 2 2\times2 2 × 2 矩阵都满足 A X A = A AXA=A A X A = A (因为 A X A = ( ∑ i j x i j ) A AXA=(\sum_{ij}x_{ij})\,A A X A = ( ∑ ij x ij ) A )。∎
C15 · 判断是否为满秩分解
指出下列各式是否为满秩分解(说明理由):
(a) [ 1 2 3 4 ] = [ 1 2 ] [ 1 2 ] \begin{bmatrix}1&2\\ 3&4\end{bmatrix}=\begin{bmatrix}1\\2\end{bmatrix}[1\ \ 2] [ 1 3 2 4 ] = [ 1 2 ] [ 1 2 ] ;(b) [ 1 2 3 4 ] = I 2 [ 1 2 3 4 ] \begin{bmatrix}1&2\\ 3&4\end{bmatrix}=I_2\begin{bmatrix}1&2\\ 3&4\end{bmatrix} [ 1 3 2 4 ] = I 2 [ 1 3 2 4 ] ;(c) [ 2 4 1 2 ] = [ 2 1 ] [ 1 2 ] \begin{bmatrix}2&4\\ 1&2\end{bmatrix}=\begin{bmatrix}2\\1\end{bmatrix}[1\ \ 2] [ 2 1 4 2 ] = [ 2 1 ] [ 1 2 ] ;(d) [ 2 4 1 2 ] = [ 1 1 ] [ 2 4 ] \begin{bmatrix}2&4\\ 1&2\end{bmatrix}=\begin{bmatrix}1\\1\end{bmatrix}[2\ \ 4] [ 2 1 4 2 ] = [ 1 1 ] [ 2 4 ] 。
解答 · C15
(a) 不是 :右边乘积 = [ 1 2 2 4 ] ≠ A =\begin{bmatrix}1&2\\2&4\end{bmatrix}\neq A = [ 1 2 2 4 ] = A (且 A A A 可逆、秩 2 2 2 ,而右边秩 ≤ 1 \le1 ≤ 1 ——秩都不等,必不是)。
(b) 是 :I 2 I_2 I 2 列满秩(秩 2 2 2 )、A A A 行满秩(2 × 2 2\times2 2 × 2 秩 2 2 2 ),而 A = I 2 A A=I_2A A = I 2 A ——注意这正是 注 2 (注 2)型的「平凡但合法」的分解。
(c) 是 :B = [ 2 1 ] B=\begin{bmatrix}2\\1\end{bmatrix} B = [ 2 1 ] 列满秩(1 1 1 列,非零)、C = [ 1 2 ] C=[1\ \ 2] C = [ 1 2 ] 行满秩(1 1 1 行,非零),乘积 = [ 2 4 1 2 ] = A =\begin{bmatrix}2&4\\1&2\end{bmatrix}=A = [ 2 1 4 2 ] = A ✓;r = 1 r=1 r = 1 ✓。
(d) 不是 :乘积 = [ 2 4 2 4 ] ≠ [ 2 4 1 2 ] =\begin{bmatrix}2&4\\2&4\end{bmatrix}\neq\begin{bmatrix}2&4\\1&2\end{bmatrix} = [ 2 2 4 4 ] = [ 2 1 4 2 ] (B B B 的列不在 R ( A ) R(A) R ( A ) 的方向上:R ( A ) = span { ( 2 , 1 ) T } R(A)=\operatorname{span}\{(2,1)^{\mathrm T}\} R ( A ) = span {( 2 , 1 ) T } ,而 ( 1 , 1 ) T (1,1)^{\mathrm T} ( 1 , 1 ) T 不属于它)。
要点 :判断满秩分解要过三关——乘积等于 A A A 、B B B 列满秩、C C C 行满秩;且 B B B 的列必须落在 R ( A ) R(A) R ( A ) 里((d) 就是这一关不过)。
C16 · 从相抵分解读出满秩分解
设 A A A 如 C3 (3 × 3 3\times3 3 × 3 、秩 2 2 2 )。用相抵分解的「方法 1」(定理 3.2.1 )从 E − 1 E^{-1} E − 1 、F − 1 F^{-1} F − 1 直接读出 A A A 的一个满秩分解,并验算。
解答 · C16
方法 1 的配方:A = P − 1 Γ Q = P − 1 [ I r 0 ] [ I r , 0 ] Q A=P^{-1}\Gamma Q=P^{-1}\begin{bmatrix}I_r\\0\end{bmatrix}[I_r,\ 0]Q A = P − 1 Γ Q = P − 1 [ I r 0 ] [ I r , 0 ] Q ,故
B = P − 1 [ I r 0 ] = ( P − 1 的前 r 列) , C = [ I r , 0 ] Q = ( Q 的前 r 行) . B=P^{-1}\begin{bmatrix}I_r\\0\end{bmatrix}=\text{(}P^{-1}\ \text{的前 } r \text{ 列)},\qquad C=[I_r,\ 0]Q=\text{(}Q\ \text{的前 } r \text{ 行)}. B = P − 1 [ I r 0 ] = ( P − 1 的前 r 列) , C = [ I r , 0 ] Q = ( Q 的前 r 行) .
代入 C3 的结果 P − 1 = E − 1 = [ 1 0 0 2 1 0 3 1 1 ] P^{-1}=E^{-1}=\begin{bmatrix}1&0&0\\2&1&0\\3&1&1\end{bmatrix} P − 1 = E − 1 = 1 2 3 0 1 1 0 0 1 、Q = F − 1 = [ 1 2 3 0 1 1 0 0 1 ] Q=F^{-1}=\begin{bmatrix}1&2&3\\0&1&1\\0&0&1\end{bmatrix} Q = F − 1 = 1 0 0 2 1 0 3 1 1 ,取前 r = 2 r=2 r = 2 列 / 行:
B = [ 1 0 2 1 3 1 ] , C = [ 1 2 3 0 1 1 ] . B=\begin{bmatrix}1&0\\ 2&1\\ 3&1\end{bmatrix},\qquad C=\begin{bmatrix}1&2&3\\ 0&1&1\end{bmatrix}. B = 1 2 3 0 1 1 , C = [ 1 0 2 1 3 1 ] .
验算 :
B C = [ 1 0 2 1 3 1 ] [ 1 2 3 0 1 1 ] = [ 1 2 3 2 5 7 3 7 10 ] = A ✓ BC=\begin{bmatrix}1&0\\2&1\\3&1\end{bmatrix}\begin{bmatrix}1&2&3\\0&1&1\end{bmatrix}=\begin{bmatrix}1&2&3\\2&5&7\\3&7&10\end{bmatrix}=A\ \checkmark BC = 1 2 3 0 1 1 [ 1 0 2 1 3 1 ] = 1 2 3 2 5 7 3 7 10 = A ✓
且 B B B 列满秩(两列无关)、C C C 行满秩(两行无关)✓。
读信息 :R ( A ) = span { ( 1 , 2 , 3 ) T , ( 0 , 1 , 1 ) T } R(A)=\operatorname{span}\{(1,2,3)^{\mathrm T},(0,1,1)^{\mathrm T}\} R ( A ) = span {( 1 , 2 , 3 ) T , ( 0 , 1 , 1 ) T } ——注意这不是 A A A 的前两列 (A A A 的前两列是 ( 1 , 2 , 3 ) T , ( 2 , 5 , 7 ) T (1,2,3)^{\mathrm T},(2,5,7)^{\mathrm T} ( 1 , 2 , 3 ) T , ( 2 , 5 , 7 ) T ):同一组基可以由不同的向量担任,正印证了注 2(注 2 )「基不一定从 A A A 的列中选取」。∎
C17 · GS 手算 QR 分解(2 × 2 2\times2 2 × 2 )
求 A = [ 2 1 1 2 ] A=\begin{bmatrix}2&1\\ 1&2\end{bmatrix} A = [ 2 1 1 2 ] 的 QR 分解(Q Q Q 正交、R R R 正线上三角),并回代验证。
提示 · C17
照 定理 3.4.1 存在性证明的流程:先正交化得 y k \boldsymbol{y}_k y k ,对角元取 ∥ y k ∥ \|\boldsymbol{y}_k\| ∥ y k ∥ (残差的范数,不是 ∥ a k ∥ \|\boldsymbol{a}_k\| ∥ a k ∥ ),非对角元取 ( a k , z i ) (\boldsymbol{a}_k,\boldsymbol{z}_i) ( a k , z i ) 。
解答 · C17
det A = 3 ≠ 0 \det A=3\neq0 det A = 3 = 0 ,满秩 ✓。a 1 = ( 2 , 1 ) T \boldsymbol{a}_1=(2,1)^{\mathrm T} a 1 = ( 2 , 1 ) T 、a 2 = ( 1 , 2 ) T \boldsymbol{a}_2=(1,2)^{\mathrm T} a 2 = ( 1 , 2 ) T 。
GS :∥ y 1 ∥ = ∥ a 1 ∥ = 5 \|\boldsymbol{y}_1\|=\|\boldsymbol{a}_1\|=\sqrt5 ∥ y 1 ∥ = ∥ a 1 ∥ = 5 ,z 1 = 1 5 ( 2 , 1 ) T \boldsymbol{z}_1=\tfrac{1}{\sqrt5}(2,1)^{\mathrm T} z 1 = 5 1 ( 2 , 1 ) T ;( a 2 , z 1 ) = 4 5 (\boldsymbol{a}_2,\boldsymbol{z}_1)=\tfrac{4}{\sqrt5} ( a 2 , z 1 ) = 5 4 ,
y 2 = a 2 − 4 5 z 1 = ( 1 , 2 ) T − 4 5 ( 2 , 1 ) T = ( − 3 5 , 6 5 ) T , ∥ y 2 ∥ = 3 5 , z 2 = 1 5 ( − 1 , 2 ) T . \boldsymbol{y}_2=\boldsymbol{a}_2-\tfrac{4}{\sqrt5}\boldsymbol{z}_1=(1,2)^{\mathrm T}-\tfrac45(2,1)^{\mathrm T}=\left(-\tfrac35,\ \tfrac65\right)^{\mathrm T},\qquad \|\boldsymbol{y}_2\|=\frac{3}{\sqrt5},\quad \boldsymbol{z}_2=\frac{1}{\sqrt5}(-1,2)^{\mathrm T}. y 2 = a 2 − 5 4 z 1 = ( 1 , 2 ) T − 5 4 ( 2 , 1 ) T = ( − 5 3 , 5 6 ) T , ∥ y 2 ∥ = 5 3 , z 2 = 5 1 ( − 1 , 2 ) T .
拼装 :
Q = 1 5 [ 2 − 1 1 2 ] , R = [ 5 4 5 0 3 5 ] . Q=\frac{1}{\sqrt5}\begin{bmatrix}2&-1\\ 1&2\end{bmatrix},\qquad R=\begin{bmatrix}\sqrt5&\tfrac{4}{\sqrt5}\\ 0&\tfrac{3}{\sqrt5}\end{bmatrix}. Q = 5 1 [ 2 1 − 1 2 ] , R = [ 5 0 5 4 5 3 ] .
验算 :Q R = [ 2 8 5 − 3 5 1 4 5 + 6 5 ] = [ 2 1 1 2 ] = A QR=\begin{bmatrix}2&\tfrac85-\tfrac35\\ 1&\tfrac45+\tfrac65\end{bmatrix}=\begin{bmatrix}2&1\\ 1&2\end{bmatrix}=A QR = [ 2 1 5 8 − 5 3 5 4 + 5 6 ] = [ 2 1 1 2 ] = A ✓;Q T Q = I Q^{\mathrm T}Q=I Q T Q = I ✓、det Q = 4 + 1 5 = 1 \det Q=\tfrac{4+1}{5}=1 det Q = 5 4 + 1 = 1 (这又是一个 Givens 旋转,Q = [ cos θ − sin θ sin θ cos θ ] Q=\begin{bmatrix}\cos\theta&-\sin\theta\\\sin\theta&\cos\theta\end{bmatrix} Q = [ cos θ sin θ − sin θ cos θ ] ,cos θ = 2 5 \cos\theta=\tfrac2{\sqrt5} cos θ = 5 2 ——对称矩阵的 QR 因子常常「长得好看」)。∎
C18 · 长方阵的完全 QR 分解(推论 3.4.1 实操)
设 B = [ 1 0 0 1 1 1 ] B=\begin{bmatrix}1&0\\ 0&1\\ 1&1\end{bmatrix} B = 1 0 1 0 1 1 。(a) 求 B = Q 1 R 1 B=Q_1R_1 B = Q 1 R 1 (Q 1 Q_1 Q 1 列正交规范、R 1 R_1 R 1 正线上三角);(b) 补齐 C 3 \mathbb C^3 C 3 的标准正交基,写出 B = Q R B=QR B = QR (Q Q Q 酉、R = [ R 1 O ] R=\begin{bmatrix}R_1\\O\end{bmatrix} R = [ R 1 O ] ),并验证 Q H Q = I 3 Q^{\mathrm H}Q=I_3 Q H Q = I 3 。
解答 · C18
(a) a 1 = ( 1 , 0 , 1 ) T \boldsymbol{a}_1=(1,0,1)^{\mathrm T} a 1 = ( 1 , 0 , 1 ) T :∥ a 1 ∥ = 2 \|\boldsymbol{a}_1\|=\sqrt2 ∥ a 1 ∥ = 2 ,z 1 = 1 2 ( 1 , 0 , 1 ) T \boldsymbol{z}_1=\tfrac{1}{\sqrt2}(1,0,1)^{\mathrm T} z 1 = 2 1 ( 1 , 0 , 1 ) T 。( a 2 , z 1 ) = 1 2 (\boldsymbol{a}_2,\boldsymbol{z}_1)=\tfrac{1}{\sqrt2} ( a 2 , z 1 ) = 2 1 (a 2 = ( 0 , 1 , 1 ) T \boldsymbol{a}_2=(0,1,1)^{\mathrm T} a 2 = ( 0 , 1 , 1 ) T ):
y 2 = a 2 − 1 2 ( 1 , 0 , 1 ) T = ( − 1 2 , 1 , 1 2 ) T , ∥ y 2 ∥ = 6 2 , z 2 = 1 6 ( − 1 , 2 , 1 ) T . \boldsymbol{y}_2=\boldsymbol{a}_2-\tfrac12(1,0,1)^{\mathrm T}=\left(-\tfrac12,\ 1,\ \tfrac12\right)^{\mathrm T},\qquad \|\boldsymbol{y}_2\|=\frac{\sqrt6}{2},\quad \boldsymbol{z}_2=\frac{1}{\sqrt6}(-1,2,1)^{\mathrm T}. y 2 = a 2 − 2 1 ( 1 , 0 , 1 ) T = ( − 2 1 , 1 , 2 1 ) T , ∥ y 2 ∥ = 2 6 , z 2 = 6 1 ( − 1 , 2 , 1 ) T .
Q 1 = [ 1 2 − 1 6 0 2 6 1 2 1 6 ] , R 1 = [ 2 1 2 0 6 2 ] , Q 1 H Q 1 = I 2 ✓ Q_1=\begin{bmatrix}\tfrac{1}{\sqrt2}&-\tfrac{1}{\sqrt6}\\ 0&\tfrac{2}{\sqrt6}\\ \tfrac{1}{\sqrt2}&\tfrac{1}{\sqrt6}\end{bmatrix},\qquad R_1=\begin{bmatrix}\sqrt2&\tfrac{1}{\sqrt2}\\ 0&\tfrac{\sqrt6}{2}\end{bmatrix},\qquad Q_1^{\mathrm H}Q_1=I_2\ \checkmark Q 1 = 2 1 0 2 1 − 6 1 6 2 6 1 , R 1 = [ 2 0 2 1 2 6 ] , Q 1 H Q 1 = I 2 ✓
验算 Q 1 R 1 Q_1R_1 Q 1 R 1 :第一列 2 z 1 = ( 1 , 0 , 1 ) T \sqrt2\boldsymbol{z}_1=(1,0,1)^{\mathrm T} 2 z 1 = ( 1 , 0 , 1 ) T ✓;第二列 1 2 z 1 + 6 2 z 2 = ( 1 2 , 0 , 1 2 ) T + ( − 1 2 , 1 , 1 2 ) T = ( 0 , 1 , 1 ) T \tfrac{1}{\sqrt2}\boldsymbol{z}_1+\tfrac{\sqrt6}2\boldsymbol{z}_2=(\tfrac12,0,\tfrac12)^{\mathrm T}+(-\tfrac12,1,\tfrac12)^{\mathrm T}=(0,1,1)^{\mathrm T} 2 1 z 1 + 2 6 z 2 = ( 2 1 , 0 , 2 1 ) T + ( − 2 1 , 1 , 2 1 ) T = ( 0 , 1 , 1 ) T ✓。
(b) 取 z 3 \boldsymbol{z}_3 z 3 同时垂直于 z 1 ∥ ( 1 , 0 , 1 ) \boldsymbol{z}_1\parallel(1,0,1) z 1 ∥ ( 1 , 0 , 1 ) 与 z 2 ∥ ( − 1 , 2 , 1 ) \boldsymbol{z}_2\parallel(-1,2,1) z 2 ∥ ( − 1 , 2 , 1 ) :叉积 ( 1 , 0 , 1 ) × ( − 1 , 2 , 1 ) = ( − 2 , − 2 , 2 ) (1,0,1)\times(-1,2,1)=(-2,-2,2) ( 1 , 0 , 1 ) × ( − 1 , 2 , 1 ) = ( − 2 , − 2 , 2 ) ,单位化得 z 3 = 1 3 ( 1 , 1 , − 1 ) T \boldsymbol{z}_3=\tfrac{1}{\sqrt3}(1,1,-1)^{\mathrm T} z 3 = 3 1 ( 1 , 1 , − 1 ) T (两两内积为零 ✓)。于是
Q = [ 1 2 − 1 6 1 3 0 2 6 1 3 1 2 1 6 − 1 3 ] ( 酉 ) , B = Q [ 2 1 2 0 6 2 0 0 ] ✓ Q=\begin{bmatrix}\tfrac{1}{\sqrt2}&-\tfrac{1}{\sqrt6}&\tfrac{1}{\sqrt3}\\ 0&\tfrac{2}{\sqrt6}&\tfrac{1}{\sqrt3}\\ \tfrac{1}{\sqrt2}&\tfrac{1}{\sqrt6}&-\tfrac{1}{\sqrt3}\end{bmatrix}\ (\text{酉}),\qquad B=Q\begin{bmatrix}\sqrt2&\tfrac{1}{\sqrt2}\\ 0&\tfrac{\sqrt6}{2}\\ 0&0\end{bmatrix}\ \checkmark Q = 2 1 0 2 1 − 6 1 6 2 6 1 3 1 3 1 − 3 1 ( 酉 ) , B = Q 2 0 0 2 1 2 6 0 ✓
Q H Q = I 3 Q^{\mathrm H}Q=I_3 Q H Q = I 3 逐列验算 ✓。∎
C19 · 高次幂的两条路
设 A = [ 4 − 1 2 1 ] A=\begin{bmatrix}4&-1\\ 2&1\end{bmatrix} A = [ 4 2 − 1 1 ] 。(a) 求 f A ( λ ) f_A(\lambda) f A ( λ ) ;(b) 用 H–C 定理把 A 10 A^{10} A 10 写成 A A A 的一次 多项式并算出显式结果;(c) 该矩阵可对角化吗?用对角化路复核 (b) 的系数。
解答 · C19
(a) tr A = 5 \operatorname{tr}A=5 tr A = 5 、det A = 6 \det A=6 det A = 6 ,f A ( λ ) = λ 2 − 5 λ + 6 = ( λ − 2 ) ( λ − 3 ) f_A(\lambda)=\lambda^2-5\lambda+6=(\lambda-2)(\lambda-3) f A ( λ ) = λ 2 − 5 λ + 6 = ( λ − 2 ) ( λ − 3 ) 。
(b) H–C:A 2 = 5 A − 6 I A^2=5A-6I A 2 = 5 A − 6 I ,故 A 10 = α A + β I A^{10}=\alpha A+\beta I A 10 = α A + β I 。在特征值处定系数(h ( λ ) = λ 10 − α λ − β h(\lambda)=\lambda^{10}-\alpha\lambda-\beta h ( λ ) = λ 10 − α λ − β 在 2 , 3 2,3 2 , 3 处为零 ⇒ f A ∣ h ⇒ h ( A ) = O \Rightarrow f_A\mid h\Rightarrow h(A)=O ⇒ f A ∣ h ⇒ h ( A ) = O ):
2 10 = 2 α + β , 3 10 = 3 α + β ⟹ α = 3 10 − 2 10 = 58025 , β = 2 11 − 3 10 = − 115026. 2^{10}=2\alpha+\beta,\qquad 3^{10}=3\alpha+\beta\ \Longrightarrow\ \alpha=3^{10}-2^{10}=58025,\quad \beta=2^{11}-3^{10}=-115026 . 2 10 = 2 α + β , 3 10 = 3 α + β ⟹ α = 3 10 − 2 10 = 58025 , β = 2 11 − 3 10 = − 115026.
A 10 = 58025 A − 115026 I = [ 4 × 58025 − 115026 − 58025 2 × 58025 58025 − 115026 ] = [ 117074 − 58025 116050 − 57001 ] . A^{10}=58025\,A-115026\,I=\begin{bmatrix}4\times58025-115026&-58025\\ 2\times58025&58025-115026\end{bmatrix} =\begin{bmatrix}117074&-58025\\ 116050&-57001\end{bmatrix}. A 10 = 58025 A − 115026 I = [ 4 × 58025 − 115026 2 × 58025 − 58025 58025 − 115026 ] = [ 117074 116050 − 58025 − 57001 ] .
(本编数值复核:与双精度 A 10 A^{10} A 10 逐位一致 ✓。)
(c) f A f_A f A 无重根 ⇒ \Rightarrow ⇒ 可对角化(不同特征值的特征向量线性无关 2.4.5 ):λ = 2 \lambda=2 λ = 2 时 ( A − 2 I ) x = 0 (A-2I)\boldsymbol{x}=\boldsymbol{0} ( A − 2 I ) x = 0 给 p 1 = ( 1 , 2 ) T \boldsymbol{p}_1=(1,2)^{\mathrm T} p 1 = ( 1 , 2 ) T ;λ = 3 \lambda=3 λ = 3 给 p 2 = ( 1 , 1 ) T \boldsymbol{p}_2=(1,1)^{\mathrm T} p 2 = ( 1 , 1 ) T 。P = [ p 1 , p 2 ] P=[\boldsymbol{p}_1,\boldsymbol{p}_2] P = [ p 1 , p 2 ] ,P − 1 A P = diag ( 2 , 3 ) P^{-1}AP=\operatorname{diag}(2,3) P − 1 A P = diag ( 2 , 3 ) ,A 10 = P diag ( 2 10 , 3 10 ) P − 1 A^{10}=P\operatorname{diag}(2^{10},3^{10})P^{-1} A 10 = P diag ( 2 10 , 3 10 ) P − 1 。系数复核:α , β \alpha,\beta α , β 满足的方程组正是「λ m = α λ + β \lambda^m=\alpha\lambda+\beta λ m = α λ + β 在 λ = 2 , 3 \lambda=2,3 λ = 2 , 3 处成立」,两条路的定系数步骤相同——对角化路把 P P P 显式算出来,H–C 路根本不需要 P P P (且对不可对角化的矩阵照用,见 例 3.5.2 )。∎
C20 · 带余除法与最小多项式算 A 4 A^4 A 4
设 A = [ 1 1 0 0 1 0 0 0 − 1 ] A=\begin{bmatrix}1&1&0\\ 0&1&0\\ 0&0&-1\end{bmatrix} A = 1 0 0 1 1 0 0 0 − 1 。(a) 求 f A f_A f A 与 m A m_A m A ;(b) 用「除以 m A m_A m A 」计算 A 4 A^4 A 4 ;(c) 直接计算 A 4 A^4 A 4 核对。
解答 · C20
(a) 分块三角:f A ( λ ) = ( λ − 1 ) 2 ( λ + 1 ) f_A(\lambda)=(\lambda-1)^2(\lambda+1) f A ( λ ) = ( λ − 1 ) 2 ( λ + 1 ) 。候选 ( λ − 1 ) ( λ + 1 ) (\lambda-1)(\lambda+1) ( λ − 1 ) ( λ + 1 ) 、( λ − 1 ) 2 ( λ + 1 ) (\lambda-1)^2(\lambda+1) ( λ − 1 ) 2 ( λ + 1 ) 。试低次:( A − I ) ( A + I ) = [ 0 2 0 0 0 0 0 0 0 ] ≠ O (A-I)(A+I)=\begin{bmatrix}0&2&0\\0&0&0\\0&0&0\end{bmatrix}\neq O ( A − I ) ( A + I ) = 0 0 0 2 0 0 0 0 0 = O ,落选;故 m A = ( λ − 1 ) 2 ( λ + 1 ) = f A m_A=(\lambda-1)^2(\lambda+1)=f_A m A = ( λ − 1 ) 2 ( λ + 1 ) = f A 。
(b) λ 4 = m A ⋅ ( λ + 1 ) + ( 2 λ 2 − 1 ) \lambda^4=m_A\cdot(\lambda+1)+\big(2\lambda^2-1\big) λ 4 = m A ⋅ ( λ + 1 ) + ( 2 λ 2 − 1 ) (验算:( λ − 1 ) 2 ( λ + 1 ) 2 + 2 λ 2 − 1 = ( λ 2 − 1 ) 2 + 2 λ 2 − 1 = λ 4 (\lambda-1)^2(\lambda+1)^2+2\lambda^2-1=(\lambda^2-1)^2+2\lambda^2-1=\lambda^4 ( λ − 1 ) 2 ( λ + 1 ) 2 + 2 λ 2 − 1 = ( λ 2 − 1 ) 2 + 2 λ 2 − 1 = λ 4 ✓)。故
A 4 = m A ( A ) ( A + I ) + 2 A 2 − I = 2 A 2 − I . A^4=m_A(A)(A+I)+2A^2-I=2A^2-I . A 4 = m A ( A ) ( A + I ) + 2 A 2 − I = 2 A 2 − I .
(c) A 2 = [ 1 2 0 0 1 0 0 0 1 ] A^2=\begin{bmatrix}1&2&0\\0&1&0\\0&0&1\end{bmatrix} A 2 = 1 0 0 2 1 0 0 0 1 (( − 1 ) 2 = 1 (-1)^2=1 ( − 1 ) 2 = 1 ),故 A 4 = 2 A 2 − I = [ 1 4 0 0 1 0 0 0 1 ] A^4=2A^2-I=\begin{bmatrix}1&4&0\\0&1&0\\0&0&1\end{bmatrix} A 4 = 2 A 2 − I = 1 0 0 4 1 0 0 0 1 。直接核对:A 4 = ( A 2 ) 2 A^4=(A^2)^2 A 4 = ( A 2 ) 2 ,上三角块 [ 1 2 0 1 ] 2 = [ 1 4 0 1 ] \begin{bmatrix}1&2\\0&1\end{bmatrix}^2=\begin{bmatrix}1&4\\0&1\end{bmatrix} [ 1 0 2 1 ] 2 = [ 1 0 4 1 ] ✓,右下 ( − 1 ) 4 = 1 (−1)^4=1 ( − 1 ) 4 = 1 ✓。∎
C21 · 候选法求最小多项式(低次落选型)
求 A = [ 2 1 0 0 2 1 0 0 4 ] A=\begin{bmatrix}2&1&0\\ 0&2&1\\ 0&0&4\end{bmatrix} A = 2 0 0 1 2 0 0 1 4 的最小多项式。
解答 · C21
f A ( λ ) = ( λ − 2 ) 2 ( λ − 4 ) f_A(\lambda)=(\lambda-2)^2(\lambda-4) f A ( λ ) = ( λ − 2 ) 2 ( λ − 4 ) (上三角,对角 2 , 2 , 4 2,2,4 2 , 2 , 4 )。候选:( λ − 2 ) ( λ − 4 ) (\lambda-2)(\lambda-4) ( λ − 2 ) ( λ − 4 ) 与 ( λ − 2 ) 2 ( λ − 4 ) (\lambda-2)^2(\lambda-4) ( λ − 2 ) 2 ( λ − 4 ) 。
试低次 :A − 2 I = [ 0 1 0 0 0 1 0 0 2 ] A-2I=\begin{bmatrix}0&1&0\\0&0&1\\0&0&2\end{bmatrix} A − 2 I = 0 0 0 1 0 0 0 1 2 ,A − 4 I = [ − 2 1 0 0 − 2 1 0 0 0 ] A-4I=\begin{bmatrix}-2&1&0\\0&-2&1\\0&0&0\end{bmatrix} A − 4 I = − 2 0 0 1 − 2 0 0 1 0 。乘积第一行 = A − 2 I =A-2I = A − 2 I 第一行 ( 0 , 1 , 0 ) (0,1,0) ( 0 , 1 , 0 ) 乘开 = A − 4 I =A-4I = A − 4 I 的第二行 = ( 0 , − 2 , 1 ) ≠ 0 =(0,-2,1)\neq\boldsymbol{0} = ( 0 , − 2 , 1 ) = 0 ——落选。
试高次 :( A − 2 I ) 2 = [ 0 0 1 0 0 2 0 0 4 ] (A-2I)^2=\begin{bmatrix}0&0&1\\0&0&2\\0&0&4\end{bmatrix} ( A − 2 I ) 2 = 0 0 0 0 0 0 1 2 4 ,( A − 2 I ) 2 ( A − 4 I ) (A-2I)^2(A-4I) ( A − 2 I ) 2 ( A − 4 I ) 的每一行都是 ( 0 , 0 , ∗ ) (0,0,*) ( 0 , 0 , ∗ ) 乘 A − 4 I A-4I A − 4 I (第三行为零行),故 = O =O = O ✓。
m A ( λ ) = ( λ − 2 ) 2 ( λ − 4 ) . m_A(\lambda)=(\lambda-2)^2(\lambda-4). m A ( λ ) = ( λ − 2 ) 2 ( λ − 4 ) .
读法 :λ = 2 \lambda=2 λ = 2 的特征子空间 N ( A − 2 I ) N(A-2I) N ( A − 2 I ) 由 e 1 \boldsymbol{e}_1 e 1 张成(解方程得 x 2 = x 3 = 0 x_2=x_3=0 x 2 = x 3 = 0 ),几何重数 1 < 1< 1 < 代数重数 2 2 2 (几何重数的计算 2.4.6 )——最小多项式里 ( λ − 2 ) (\lambda-2) ( λ − 2 ) 必须带平方,与 例 3.5.B 同型。∎
C22 · QR 迭代手算三轮(分数精确)
对 A = [ 2 1 1 2 ] A=\begin{bmatrix}2&1\\ 1&2\end{bmatrix} A = [ 2 1 1 2 ] (即 C17 )执行 QR 算法三轮:写出 A 2 A_2 A 2 、A 3 A_3 A 3 (精确分数)与 A 4 A_4 A 4 (四位小数),观察对角元向何值收敛、非对角元衰减多快。
解答 · C22
第一轮 :A 1 = Q 1 R 1 A_1=Q_1R_1 A 1 = Q 1 R 1 用 C17 的结果,A 2 = R 1 Q 1 A_2=R_1Q_1 A 2 = R 1 Q 1 :
A 2 = [ 5 4 5 0 3 5 ] ⋅ 1 5 [ 2 − 1 1 2 ] = [ 14 5 3 5 3 5 6 5 ] = [ 2.8 0.6 0.6 1.2 ] . A_2=\begin{bmatrix}\sqrt5&\tfrac{4}{\sqrt5}\\0&\tfrac{3}{\sqrt5}\end{bmatrix}\cdot\frac{1}{\sqrt5}\begin{bmatrix}2&-1\\1&2\end{bmatrix} =\begin{bmatrix}\tfrac{14}{5}&\tfrac35\\ \tfrac35&\tfrac65\end{bmatrix} =\begin{bmatrix}2.8&0.6\\ 0.6&1.2\end{bmatrix}. A 2 = [ 5 0 5 4 5 3 ] ⋅ 5 1 [ 2 1 − 1 2 ] = [ 5 14 5 3 5 3 5 6 ] = [ 2.8 0.6 0.6 1.2 ] .
第二轮 :对 A 2 A_2 A 2 再 GS(z 1 = 1 205 ( 14 , 3 ) T \boldsymbol{z}_1=\tfrac1{\sqrt{205}}(14,3)^{\mathrm T} z 1 = 205 1 ( 14 , 3 ) T ,( a 2 , z 1 ) = 12 205 (\boldsymbol{a}_2,\boldsymbol{z}_1)=\tfrac{12}{\sqrt{205}} ( a 2 , z 1 ) = 205 12 ,y 2 ∥ ( − 3 , 14 ) T \boldsymbol{y}_2\parallel(-3,14)^{\mathrm T} y 2 ∥ ( − 3 , 14 ) T 、∥ y 2 ∥ = 3 205 41 \|\boldsymbol{y}_2\|=\tfrac{3\sqrt{205}}{41} ∥ y 2 ∥ = 41 3 205 ):
Q 2 = 1 205 [ 14 − 3 3 14 ] , R 2 = [ 205 5 12 205 0 3 205 41 ] , A 3 = R 2 Q 2 = [ 122 41 9 41 9 41 42 41 ] ≈ [ 2.9756 0.2195 0.2195 1.0244 ] . Q_2=\frac{1}{\sqrt{205}}\begin{bmatrix}14&-3\\ 3&14\end{bmatrix},\qquad R_2=\begin{bmatrix}\tfrac{\sqrt{205}}{5}&\tfrac{12}{\sqrt{205}}\\ 0&\tfrac{3\sqrt{205}}{41}\end{bmatrix},\qquad A_3=R_2Q_2=\begin{bmatrix}\tfrac{122}{41}&\tfrac9{41}\\ \tfrac9{41}&\tfrac{42}{41}\end{bmatrix} \approx\begin{bmatrix}2.9756&0.2195\\ 0.2195&1.0244\end{bmatrix}. Q 2 = 205 1 [ 14 3 − 3 14 ] , R 2 = [ 5 205 0 205 12 41 3 205 ] , A 3 = R 2 Q 2 = [ 41 122 41 9 41 9 41 42 ] ≈ [ 2.9756 0.2195 0.2195 1.0244 ] .
第三轮 (数值):A 4 ≈ [ 2.9973 − 0.0740 − 0.0740 1.0027 ] A_4\approx\begin{bmatrix}2.9973&-0.0740\\ -0.0740&1.0027\end{bmatrix} A 4 ≈ [ 2.9973 − 0.0740 − 0.0740 1.0027 ] 。
观察 :特征值精确解 3 , 1 3,1 3 , 1 (f A = ( λ − 3 ) ( λ − 1 ) f_A=(\lambda-3)(\lambda-1) f A = ( λ − 3 ) ( λ − 1 ) ,模不同 ⇒ \Rightarrow ⇒ 收敛条件满足)。对角元 2.8 → 2.976 → 2.997 → 3 2.8\to2.976\to2.997\to3 2.8 → 2.976 → 2.997 → 3 、1.2 → 1.024 → 1.003 → 1 1.2\to1.024\to1.003\to1 1.2 → 1.024 → 1.003 → 1 (按模从大到小排定 ✓);非对角元 0.6 → 0.22 → 0.074 0.6\to0.22\to0.074 0.6 → 0.22 → 0.074 ——每轮约缩小到三分之一(收敛率由相邻特征值模比 ∣ λ 2 ∣ ∣ λ 1 ∣ = 1 3 \tfrac{|\lambda_2|}{|\lambda_1|}=\tfrac13 ∣ λ 1 ∣ ∣ λ 2 ∣ = 3 1 决定,本编注)。∎
C23 · 用推论 3.5.2 求逆
设 A = [ 1 2 0 0 1 0 0 0 2 ] A=\begin{bmatrix}1&2&0\\ 0&1&0\\ 0&0&2\end{bmatrix} A = 1 0 0 2 1 0 0 0 2 。(a) 求 f A ( λ ) f_A(\lambda) f A ( λ ) ;(b) 用逆矩阵公式求 A − 1 A^{-1} A − 1 ;(c) 回代验证。
解答 · C23
(a) 上三角:f A ( λ ) = ( λ − 1 ) 2 ( λ − 2 ) = λ 3 − 4 λ 2 + 5 λ − 2 f_A(\lambda)=(\lambda-1)^2(\lambda-2)=\lambda^3-4\lambda^2+5\lambda-2 f A ( λ ) = ( λ − 1 ) 2 ( λ − 2 ) = λ 3 − 4 λ 2 + 5 λ − 2 。
(b) n = 3 n=3 n = 3 ,a 2 = − 4 a_2=-4 a 2 = − 4 、a 1 = 5 a_1=5 a 1 = 5 、a 0 = − 2 ≠ 0 a_0=-2\neq0 a 0 = − 2 = 0 (det A = 2 \det A=2 det A = 2 ,( − 1 ) 3 det A = − 2 (-1)^3\det A=-2 ( − 1 ) 3 det A = − 2 ✓):
A − 1 = − 1 a 0 ( A 2 + a 2 A + a 1 I ) = 1 2 ( A 2 − 4 A + 5 I ) . A^{-1}=-\frac{1}{a_0}\big(A^2+a_2A+a_1I\big)=\frac12\big(A^2-4A+5I\big). A − 1 = − a 0 1 ( A 2 + a 2 A + a 1 I ) = 2 1 ( A 2 − 4 A + 5 I ) .
A 2 = [ 1 4 0 0 1 0 0 0 4 ] A^2=\begin{bmatrix}1&4&0\\0&1&0\\0&0&4\end{bmatrix} A 2 = 1 0 0 4 1 0 0 0 4 ,故
A 2 − 4 A + 5 I = [ 2 − 4 0 0 2 0 0 0 1 ] , A − 1 = [ 1 − 2 0 0 1 0 0 0 1 2 ] . A^2-4A+5I=\begin{bmatrix}2&-4&0\\0&2&0\\0&0&1\end{bmatrix},\qquad A^{-1}=\begin{bmatrix}1&-2&0\\ 0&1&0\\ 0&0&\tfrac12\end{bmatrix}. A 2 − 4 A + 5 I = 2 0 0 − 4 2 0 0 0 1 , A − 1 = 1 0 0 − 2 1 0 0 0 2 1 .
(c) A A − 1 = [ 1 2 0 0 1 0 0 0 2 ] [ 1 − 2 0 0 1 0 0 0 1 2 ] = [ 1 − 2 + 2 0 0 1 0 0 0 1 ] = I AA^{-1}=\begin{bmatrix}1&2&0\\0&1&0\\0&0&2\end{bmatrix}\begin{bmatrix}1&-2&0\\0&1&0\\0&0&\frac12\end{bmatrix}=\begin{bmatrix}1&-2+2&0\\0&1&0\\0&0&1\end{bmatrix}=I A A − 1 = 1 0 0 2 1 0 0 0 2 1 0 0 − 2 1 0 0 0 2 1 = 1 0 0 − 2 + 2 1 0 0 0 1 = I ✓。∎
C24 · 最小多项式与矩阵多项式综合
设 A = [ 4 1 0 0 4 0 0 0 2 ] A=\begin{bmatrix}4&1&0\\ 0&4&0\\ 0&0&2\end{bmatrix} A = 4 0 0 1 4 0 0 0 2 。(a) 求 f A f_A f A 与 m A m_A m A ;(b) 用 m A m_A m A 计算 A 3 A^3 A 3 ;(c) 直接乘 A 3 A^3 A 3 核对。
解答 · C24
(a) f A ( λ ) = ( λ − 4 ) 2 ( λ − 2 ) f_A(\lambda)=(\lambda-4)^2(\lambda-2) f A ( λ ) = ( λ − 4 ) 2 ( λ − 2 ) 。试低次 ( λ − 4 ) ( λ − 2 ) (\lambda-4)(\lambda-2) ( λ − 4 ) ( λ − 2 ) :( A − 4 I ) ( A − 2 I ) (A-4I)(A-2I) ( A − 4 I ) ( A − 2 I ) 的第一行 = A − 4 I =A-4I = A − 4 I 第一行 ( 0 , 1 , 0 ) (0,1,0) ( 0 , 1 , 0 ) 乘开 = A − 2 I =A-2I = A − 2 I 第二行 = ( 0 , 2 , 0 ) ≠ 0 =(0,2,0)\neq\boldsymbol{0} = ( 0 , 2 , 0 ) = 0 ,落选(λ = 4 \lambda=4 λ = 4 的几何重数 1 < 2 1<2 1 < 2 ,同 例 3.5.B 型)。故 m A = ( λ − 4 ) 2 ( λ − 2 ) = f A m_A=(\lambda-4)^2(\lambda-2)=f_A m A = ( λ − 4 ) 2 ( λ − 2 ) = f A 。
(b) λ 3 = m A ⋅ 1 + ( 10 λ 2 − 32 λ + 32 ) \lambda^3=m_A\cdot1+\big(10\lambda^2-32\lambda+32\big) λ 3 = m A ⋅ 1 + ( 10 λ 2 − 32 λ + 32 ) (验算:( λ − 4 ) 2 ( λ − 2 ) = λ 3 − 10 λ 2 + 32 λ − 32 (\lambda-4)^2(\lambda-2)=\lambda^3-10\lambda^2+32\lambda-32 ( λ − 4 ) 2 ( λ − 2 ) = λ 3 − 10 λ 2 + 32 λ − 32 ,λ 3 − ( λ 3 − 10 λ 2 + 32 λ − 32 ) = 10 λ 2 − 32 λ + 32 \lambda^3-(\lambda^3-10\lambda^2+32\lambda-32)=10\lambda^2-32\lambda+32 λ 3 − ( λ 3 − 10 λ 2 + 32 λ − 32 ) = 10 λ 2 − 32 λ + 32 ✓)。故
A 3 = 10 A 2 − 32 A + 32 I . A^3=10A^2-32A+32I . A 3 = 10 A 2 − 32 A + 32 I .
(c) A 2 = [ 16 8 0 0 16 0 0 0 4 ] A^2=\begin{bmatrix}16&8&0\\0&16&0\\0&0&4\end{bmatrix} A 2 = 16 0 0 8 16 0 0 0 4 ,10 A 2 − 32 A + 32 I = [ 160 − 128 + 32 80 − 32 0 0 160 − 128 + 32 0 0 0 40 − 64 + 32 ] = [ 64 48 0 0 64 0 0 0 8 ] 10A^2-32A+32I=\begin{bmatrix}160-128+32&80-32&0\\0&160-128+32&0\\0&0&40-64+32\end{bmatrix}=\begin{bmatrix}64&48&0\\0&64&0\\0&0&8\end{bmatrix} 10 A 2 − 32 A + 32 I = 160 − 128 + 32 0 0 80 − 32 160 − 128 + 32 0 0 0 40 − 64 + 32 = 64 0 0 48 64 0 0 0 8 。直接核对:[ 4 1 0 4 ] 3 = [ 64 48 0 64 ] \begin{bmatrix}4&1\\0&4\end{bmatrix}^3=\begin{bmatrix}64&48\\0&64\end{bmatrix} [ 4 0 1 4 ] 3 = [ 64 0 48 64 ] (超对角元 3 × 4 2 = 48 3\times4^2=48 3 × 4 2 = 48 ✓)、2 3 = 8 2^3=8 2 3 = 8 ✓。∎
依据 本课程第三章课件 3.1 相抵分解(15 页)、3.2 满秩分解(9 页)、3.4 QR 分解(16 页)、3.5 Schur 分解(22 页,后两份存于 _sources/matrix/)整理 · 四节均无录播 · 整理日期 2026-10-15 · 新增例题与全部计算题均经回代验算(C19、C22 另做双精度复核)· 3.3 及以后素材到位后本页续写。