Keyboard shortcuts

Press ← or → to navigate between chapters

Press S or / to search in the book

Press ? to show this help

Press Esc to hide this help

习题集 · 第三章(3.1–3.2 与 3.4–3.5)

三组分层: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⁡(AHA)\operatorname{rank}(A)=\operatorname{rank}(A^{\mathrm H}A) 的四个等式、满秩分解的形状与读法、两种分解各自「不唯一的程度」、AXA=AAXA=A 的解(A1–A12);QR 分解的存在 / 唯一 / 正线、列正交规范与酉矩阵的分野、QR 算法的相似性与收敛条件(A13–A16);Schur 引理的适用范围、实数域边界、谱映射与特征向量继承的单向性、最小多项式的四条事实、逆矩阵公式(A17–A20)。

B 组以相抵   ⟺  \iff 同秩、单侧逆的构造与判别、满秩分解的两条主干证明(rank⁡C=r\operatorname{rank}C=r、定理 3.2.2)与 AXA=A  ⟺  CXB=IrAXA=A\iff CXB=I_r 为重心(B1–B10,B9 把 3.1 的结论接回第一章的正定矩阵);B11–B13 补 QR 的两条主干(对角酉阵刻画、迭代相似链与等模卡死),B14–B17 补 Schur / 谱映射 / 最小多项式的证明重做(B15 给出不依赖定理 3.5.1 的直接构造路线)。

C 组训练六件基本功:相抵分解(含分块矩阵法求 P−1P^{-1}、Q−1Q^{-1})、用相抵标准形解方程组、单侧逆的构造与「全部右逆」的写法、满秩分解的三种算法(取列 / RREF / 从相抵分解读)、GS 手算 QR 与长方阵补齐(C17、C18、C22)、H–C 与最小多项式驱动的三套计算(C19–C21、C23、C24)。3.3 及以后素材到位后,本页按同样的三组结构续写。

A 组 · 概念判断

A1 · 相抵分解的四个基本事实

判断对错并说明理由:

(a) 任何矩阵都有相抵分解;(b) 矩阵的相抵分解是唯一的;(c) 矩阵的相抵标准形是唯一的;(d) 两个同型矩阵相抵   ⟺  \iff 它们的秩相等。

解答 · A1

(a) 对(引理 3.1.1;零矩阵也成立:r=0r=0 时 Γ=O\Gamma=O,取 P=Q=IP=Q=I)。(b) 错——分解式 A=P−1ΓQA=P^{-1}\Gamma Q 中的 PP、QQ 有无穷多种取法(对任意 ∗*,把 [Ir∗0I]\begin{bmatrix}I_r&*\\0&I\end{bmatrix} 并进 P−1P^{-1}、[Ir0∗I]\begin{bmatrix}I_r&0\\*&I\end{bmatrix} 并进 QQ 仍给出分解,见 3.1.1)。(c) 对——Γ\Gamma 只由秩 rr 决定(秩在相抵下不变)。(d) 对(本编补注:同秩者都相抵于同一个 Γ\Gamma,再用相抵的传递性)。

要点:「不唯一」说的是分解式,「唯一」说的是标准形——两句话各有对象,不能互换。

A2 · 相抵与相似

B=P−1AQB=P^{-1}AQ 与 B=P−1APB=P^{-1}AP 分别刻画什么?它们之间是什么关系?何时两者重合?

解答 · A2
  • B=P−1AQB=P^{-1}AQ(A,BA,B 同型)是相抵:定理 2.3.5 的换基公式——同一个线性映射、出发端换基 QQ、到达端换基 PP(P,QP,Q 相互独立),或更一般地「AA 经初等变换化为 BB」;
  • B=P−1APB=P^{-1}AP 是相似(推论 2.3.4):同一个线性变换(同一空间 V→VV\to V)换一组基,两侧用同一个 PP。

关系:相似是相抵的特例(取 Q=PQ=P),故相似必相抵;反之不然(相抵允许 P,QP,Q 独立,且不要求方阵——2×32\times3 矩阵也谈相抵,但只有同阶方阵才谈相似)。两者重合:当 A,BA,B 同阶方阵且换基在出发端与到达端是同一件事时(Q=PQ=P)。

A3 · 单侧逆的存在性判断

判断下列矩阵是否有左逆、右逆(有则说明理由,没有也说明理由):

(a) A1=[100010]A_1=\begin{bmatrix}1&0&0\\ 0&1&0\end{bmatrix};(b) A2=[100100]A_2=\begin{bmatrix}1&0\\ 0&1\\ 0&0\end{bmatrix};(c) A3=[1224]A_3=\begin{bmatrix}1&2\\ 2&4\end{bmatrix};(d) A4=[1  0]A_4=[1\ \ 0]。

解答 · A3

由 定理 3.1.1:右逆   ⟺  \iff 行满秩;左逆   ⟺  \iff 列满秩。

(a) 2×32\times3、秩 2=m2=m:行满秩,有右逆;列不满秩(2<32&lt;3),无左逆。 (b) 3×23\times2、秩 2=n2=n:列满秩,有左逆;行不满秩,无右逆。 (c) 秩 1<min⁡{2,2}1&lt;\min\{2,2\}:两者都没有(顺带:方阵行满秩   ⟺  \iff 可逆,故也可直接说「不可逆」)。 (d) 1×21\times2、秩 1=m1=m:行满秩,有右逆(如 [10]\begin{bmatrix}1\\0\end{bmatrix});无左逆。

要点:先看形状定方向(行满秩还是列满秩),再看秩是否达到该方向的上限。

A4 · 乘积为单位阵与「方阵前提」

(a) 设 AA 是 nn 阶方阵,若 AB=InAB=I_n,BB 是什么?唯一吗?(b) 设 A=[10]A=\begin{bmatrix}1&amp;0\end{bmatrix},AB=I1AB=I_1 说明什么?BB 唯一吗?(c) 把「乘积为单位阵的两矩阵必可逆」改成一个正确的命题。

解答 · A4

(a) AA 可逆且 B=A−1B=A^{-1},唯一(这正是「方阵乘积为单位阵」的结论:AB=I⇒B=A−1⇒BA=IAB=I\Rightarrow B=A^{-1}\Rightarrow BA=I)。 (b) 只说明 AA 行满秩、BB 是右逆:设 B=(b1,b2)TB=(b_1,b_2)^{\mathrm T},则 AB=b1=1AB=b_1=1,B=[1t]B=\begin{bmatrix}1\\ t\end{bmatrix}(tt 任意)——不唯一;此时 BA=[10t0]≠I2BA=\begin{bmatrix}1&amp;0\\ t&amp;0\end{bmatrix}\neq I_2,「逆矩阵」的语言完全失效。 (c) 正确说法:若两个方阵的乘积为单位阵,则这两个方阵都可逆(且互为逆矩阵)。——「方阵」是关键词;长方阵情形只能谈单侧逆(定理 3.1.1)。

A5 · 秩等式的四个判断

设 A∈Cm×nA\in\mathbb C^{m\times n}。判断对错:

(a) rank⁡(A)=rank⁡(AHA)\operatorname{rank}(A)=\operatorname{rank}(A^{\mathrm H}A) 恒成立;(b) AHAA^{\mathrm H}A 可逆   ⟺  \iff AA 列满秩;(c) rank⁡(AHA)=rank⁡(AAH)\operatorname{rank}(A^{\mathrm H}A)=\operatorname{rank}(AA^{\mathrm H}) 恒成立;(d) N(A)=N(AHA)N(A)=N(A^{\mathrm H}A)。

解答 · A5

四个全对。

(a)(c)(d) 都是 命题 3.1.1 的内容:证明主线是 N(A)=N(AHA)N(A)=N(A^{\mathrm H}A)(「⊆\subseteq」两边左乘 AHA^{\mathrm H};「⊇\supseteq」用内积正定性),再用亏加秩;对 AHA^{\mathrm H} 重复一遍得 rank⁡(AH)=rank⁡(AAH)\operatorname{rank}(A^{\mathrm H})=\operatorname{rank}(AA^{\mathrm H}),结合 rank⁡A=rank⁡AH\operatorname{rank}A=\operatorname{rank}A^{\mathrm H} 得四项等式链。 (b) 由 (a):AHAA^{\mathrm H}A(nn 阶)可逆   ⟺  \iff rank⁡(AHA)=n\operatorname{rank}(A^{\mathrm H}A)=n   ⟺  \iff rank⁡A=n\operatorname{rank}A=n   ⟺  \iff AA 列满秩。

注意:结论依赖 AHA^{\mathrm H}(实 / 复内积空间);对一般数域 FF 上没有共轭转置,此命题不适用。

A6 · 满秩分解的形状

设 A∈C3×5A\in\mathbb C^{3\times5},rank⁡A=2\operatorname{rank}A=2。若 A=BCA=BC 是满秩分解,BB、CC 分别是几乘几、秩各是多少?BB 的列与 CC 的内容各有什么几何含义?

解答 · A6

B∈C3×2B\in\mathbb C^{3\times2},rank⁡B=2\operatorname{rank}B=2(列满秩,22 列无关);C∈C2×5C\in\mathbb C^{2\times5},rank⁡C=2\operatorname{rank}C=2(行满秩,22 行无关)。

几何含义:BB 的 22 列是列空间 R(A)R(A)(dim⁡=2\dim=2)的一组基;CC 的每一列是 AA 对应列在这组基下的坐标(C=[c1,⋯ ,c5]C=[\boldsymbol{c}_1,\cdots,\boldsymbol{c}_5],ai=Bci\boldsymbol{a}_i=B\boldsymbol{c}_i)。

形状记忆:BB 的行数 =m=m、列数 =r=r;CC 的行数 =r=r、列数 =n=n——中间夹着 rr。

A7 · 满秩分解的「读法」

判断对错:

(a) R(A)=R(B)R(A)=R(B);(b) N(A)=N(C)N(A)=N(C);(c) AA 的列之间的全部线性关系都记录在 CC 里;(d) 满秩分解中的 BB 必须由 AA 的列向量组成。

解答 · A7

(a) 对(定理 3.2.2 证明的第一步:由 CC 行满秩、CC 满射得 R(A)=R(BC)=R(B)R(A)=R(BC)=R(B))。 (b) 对:Cx=0⇒BCx=0⇒Ax=0C\boldsymbol{x}=\boldsymbol{0}\Rightarrow BC\boldsymbol{x}=\boldsymbol{0}\Rightarrow A\boldsymbol{x}=\boldsymbol{0};反之 Ax=0⇒B(Cx)=0⇒Cx=0A\boldsymbol{x}=\boldsymbol{0}\Rightarrow B(C\boldsymbol{x})=\boldsymbol{0}\Rightarrow C\boldsymbol{x}=\boldsymbol{0}(BB 的列无关)。故两零空间相等——这正是 (c) 的严格含义:∑kiai=0  ⟺  ∑kici=0\sum k_i\boldsymbol{a}_i=\boldsymbol{0}\iff\sum k_i\boldsymbol{c}_i=\boldsymbol{0}。 (c) 对(即 (b) 的改述)。 (d) 错:注 2(注 2)的反例 A=I2AA=I_2A:B=I2B=I_2 的任一列都不是 AA 的列向量,但仍是合法的满秩分解。

A8 · 「不唯一」的两种程度

(a) 相抵分解不唯一到什么程度?(b) 满秩分解不唯一到什么程度?两句话的精确表述各是什么?

解答 · A8

(a) 相抵分解 A=P−1ΓQA=P^{-1}\Gamma Q:PP、QQ 各自独立地有无穷多种取法(可被形如 [Ir∗0I]\begin{bmatrix}I_r&amp;*\\0&amp;I\end{bmatrix}、[Ir0∗I]\begin{bmatrix}I_r&amp;0\\*&amp;I\end{bmatrix} 的可逆矩阵吸收);没有「相差一个矩阵」式的统一刻画,因为自由度分散在两侧。 (b) 满秩分解 A=BCA=BC:任意两种分解恰好相差一个可逆矩阵 DD:

B1=B2D,C1=D−1C2B_1=B_2D,\qquad C_1=D^{-1}C_2

(这就是 定理 3.2.2。)

且 DD 必可逆(是两组基之间的过渡矩阵)。反过来,任取可逆 DD,(B1D)(D−1C1)=A(B_1D)(D^{-1}C_1)=A 仍给出满秩分解——一个分解即生成全部(注 3)。

要点:满秩分解的不唯一被「一个 DD」精确刻画;相抵分解的不唯一更「松散」(两侧独立),但标准形唯一。

A9 · 矩阵方程 AXA=AAXA=A

(a) 若 AA 可逆,方程 AXA=AAXA=A 的解是什么?唯一吗?(b) 一般情形(A∈Crm×nA\in\mathbb C^{m\times n}_r)课件是用什么办法把方程「降维」的?用一句话说清。

解答 · A9

(a) AA 可逆时两边左乘 A−1A^{-1}、右乘 A−1A^{-1} 得 X=A−1X=A^{-1},唯一(这也是 例 3.2.2 中「不能直接约掉 AA」这句话的反面:可逆时才能约)。 (b) 用满秩分解 A=BCA=BC:方程变成 BCXBC=BCBCXBC=BC,其中 CXBCXB 是 rr 阶方阵——把 m×nm\times n 的未知量 XX 的问题「夹」成了 rr 阶方阵上的条件 CXB=IrCXB=I_r;再由 本编补注,二者等价(课件写「充分条件」,其实是充要条件)。

A10 · 秩不等式的使用

(a) 判断:rank⁡(BC)≤min⁡{rank⁡B,rank⁡C}\operatorname{rank}(BC)\le\min\{\operatorname{rank}B,\operatorname{rank}C\}。(b) 满秩分解中「rank⁡C=r\operatorname{rank}C=r」的两步夹逼各是什么?(c) 为什么 rank⁡(BC)≤rank⁡C\operatorname{rank}(BC)\le\operatorname{rank}C?

解答 · A10

(a) 对:R(BC)⊆R(B)R(BC)\subseteq R(B) 给出 rank⁡(BC)≤rank⁡B\operatorname{rank}(BC)\le\operatorname{rank}B;N(C)⊆N(BC)N(C)\subseteq N(BC) 加亏加秩给出 rank⁡(BC)≤rank⁡C\operatorname{rank}(BC)\le\operatorname{rank}C。 (b) 上界:CC 只有 rr 行,故 rank⁡C≤r\operatorname{rank}C\le r;下界:N(C)⊆N(A)N(C)\subseteq N(A)(A=BCA=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=r\operatorname{rank}C=r。 (c) 若 Cx=0C\boldsymbol{x}=\boldsymbol{0} 则 BCx=0BC\boldsymbol{x}=\boldsymbol{0},即 N(C)⊆N(BC)N(C)\subseteq N(BC);两个矩阵列数相同(都是 nn),亏加秩给出 rank⁡(BC)=n−dim⁡N(BC)≤n−dim⁡N(C)=rank⁡C\operatorname{rank}(BC)=n-\dim N(BC)\le n-\dim N(C)=\operatorname{rank}C。

A11 · 相抵标准形的「几何」翻译

用基与线性映射的语言复述 引理 3.1.1:Γ\Gamma 的存在性等价于「总可以选到两组基,使映射的矩阵长成 Γ\Gamma 的样子」。请说清这两组基是怎么选的。

解答 · A11

设 TT 是「以 AA 为矩阵」的映射(取定标准基后由 定理 2.3.1 唯一确定)。做法(引理 3.1.1 的本编证明):

  1. 出发端:取零空间 N(T)N(T) 的一组基(n−rn-r 个),用基扩充定理补成 FnF^n 的一组基 ξ1,⋯ ,ξn−r,ξn−r+1,⋯ ,ξn\xi_1,\cdots,\xi_{n-r},\xi_{n-r+1},\cdots,\xi_n;
  2. 到达端:证明后 rr 个基向量的像 ξn−r+1,⋯ ,ξn\xi_{n-r+1},\cdots,\xi_n 的像是 R(T)R(T) 的基(无关性由 ξ1,⋯ ,ξn−r\xi_1,\cdots,\xi_{n-r} 是零空间的基得出),把它们排在前面,扩充成 FmF^m 的基 η1,⋯ ,ηm\eta_1,\cdots,\eta_m;
  3. 读出矩阵:在这两组基下 T(ξ1)=⋯=T(ξn−r)=0T(\xi_1)=\cdots=T(\xi_{n-r})=\boldsymbol{0}、T(ξn−r+i)=ηiT(\xi_{n-r+i})=\eta_i,故矩阵恰为 Γ\Gamma;换基公式(定理 2.3.5)给出 A=P−1ΓQA=P^{-1}\Gamma Q。

一句话:零空间的基排前面、值域的基接后面——Γ\Gamma 就是「把 AA 的活与死的部分分开摆好」的样子。

A12 · 纠错四句

改正下列说法:

(a)「相抵分解是唯一的」;(b)「行满秩矩阵有左逆」;(c)「AHAA^{\mathrm H}A 与 AA 的秩相等,因为它们是同一个二次型的矩阵」;(d)「满秩分解的 BB、CC 可以随意取,只要 BC=ABC=A」。

解答 · A12

(a) 应说:「相抵分解式不唯一,相抵标准形唯一」——见 A1。 (b) 应说:「行满秩矩阵有右逆;列满秩矩阵有左逆」(定理 3.1.1)。方向由 AB=IAB=I(右)与 BA=IBA=I(左)决定:行满秩 AA 是「宽」的,可以右乘一个「高」的 BB 凑出 II。 (c) 二者秩相等对,但理由不对:理由不是「同一个二次型」(AHAA^{\mathrm H}A 与 AA 甚至不同型),而是 命题 3.1.1 的同解方程组证明(N(A)=N(AHA)N(A)=N(A^{\mathrm H}A),再用亏加秩)。 (d) 不够:还要 BB 列满秩、CC 行满秩(定理 3.2.1);否则如 [1234]=[12][1  2]\begin{bmatrix}1&amp;2\\3&amp;4\end{bmatrix}=\begin{bmatrix}1\\2\end{bmatrix}[1\ \ 2] 之类根本不成立,或退化成「秩不足」的分解(BB 的列不满秩 ⇒\Rightarrow BB 的列不是 R(A)R(A) 的基)。

A13 · QR 分解的四个基本事实

判断对错并说明理由:

(a) 任何满秩复方阵都有唯一的 QR 分解;(b) 不满秩的方阵一定没有 QR 分解;(c) 列满秩的长方阵 A∈Cnm×nA\in\mathbb C^{m\times n}_n(n≤mn\le m)可分解为 A=U[R1O]A=U\begin{bmatrix}R_1\\O\end{bmatrix};(d) nn 阶实满秩方阵在不要求「正线」时恰有 2n2^n 种正交三角分解。

解答 · A13

(a) 对(定理 3.4.1:满秩保证存在;「正线」约定保证唯一)。(b) 错——例 3.4.2 的 A=[1100]A=\begin{bmatrix}1&amp;1\\0&amp;0\end{bmatrix} 有 A=I2[1100]A=I_2\begin{bmatrix}1&amp;1\\0&amp;0\end{bmatrix},丢掉的只是唯一性(RR 不可逆)。(c) 对(推论 3.4.1:GS ++ 基扩充)。(d) 对——任意两种分解相差对角正交阵 D=diag⁡(±1,⋯ ,±1)D=\operatorname{diag}(\pm1,\cdots,\pm1)(2n2^n 个),见 3.4.1 的本编补注。

要点:存在性的保底比想象中宽(不满秩、长方阵都常有分解),唯一性的门槛比想象中高(要满秩加正线)。

A14 · 数清一个矩阵的全部 QR 分解

写出 A=[0110]A=\begin{bmatrix}0&amp;1\\ 1&amp;0\end{bmatrix} 的全部正交三角分解(QQ 正交、RR 可逆上三角,不要求正线),指出其中正线的版本,并说明为什么没有第五种。

解答 · A14

AA 本身是正交矩阵(置换阵,且 A2=IA^2=I),两列 e2,e1\boldsymbol{e}_2,\boldsymbol{e}_1 已标准正交。GS 直接给出正线版本 Q=AQ=A、R=IR=I(第一列 r11=∥e2∥=1r_{11}=\|\boldsymbol{e}_2\|=1,r12=(e1,e2)=0r_{12}=(\boldsymbol{e}_1,\boldsymbol{e}_2)=0,r22=1r_{22}=1)。由对角正交阵生成全部 22=42^2=4 种,形如 (A⋅D, D)\big(A\cdot D,\ D\big)(R=IR=I 被 DD 左乘后仍是上三角):

(A, I),(Adiag⁡(−1,1), diag⁡(−1,1)),(Adiag⁡(1,−1), diag⁡(1,−1)),(Adiag⁡(−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).

(验算其一:Adiag⁡(−1,1)=[01−10]A\operatorname{diag}(-1,1)=\begin{bmatrix}0&amp;1\\-1&amp;0\end{bmatrix} 正交 ✓,diag⁡(−1,1)\operatorname{diag}(-1,1) 上三角 ✓,乘积 =[0110]=A=\begin{bmatrix}0&amp;1\\1&amp;0\end{bmatrix}=A ✓。)正线版本:只有 (A,I)(A,I)(对角全 1>01&gt;0)。没有第五种:若 A=Q′R′A=Q&#x27;R&#x27; 是任一分解,则 W=R′R−1=Q′TQW=R&#x27;R^{-1}=Q&#x27;^{\mathrm T}Q 既酉又上三角,故为对角阵且 ∣wii∣=1|w_{ii}|=1(3.4.1 的归纳),实数情形 wi=±1w_i=\pm1——恰好穷尽。∎

A15 · 列正交规范矩阵辨析

设 Q∈Cm×nQ\in\mathbb C^{m\times n}(m>nm&gt;n)满足 QHQ=InQ^{\mathrm H}Q=I_n。判断:

(a) QQ 是酉矩阵;(b) 对一切 x∈Cn\boldsymbol{x}\in\mathbb C^n 有 ∥Qx∥=∥x∥\|Q\boldsymbol{x}\|=\|\boldsymbol{x}\|;(c) QQ 列满秩,QHQ^{\mathrm H} 是它的一个左逆;(d) QQH=ImQQ^{\mathrm H}=I_m。

解答 · A15

(a) 错——酉矩阵必须是方阵(正交矩阵与酉矩阵 2.5.1);QQ 只是列正交规范矩阵(定义 3.4.2)。(b) 对:∥Qx∥2=xHQHQx=∥x∥2\|Q\boldsymbol{x}\|^2=\boldsymbol{x}^{\mathrm H}Q^{\mathrm H}Q\boldsymbol{x}=\|\boldsymbol{x}\|^2。(c) 对:QHQ=InQ^{\mathrm H}Q=I_n 本身就是「列满秩 ++ 给出左逆」(定理 3.1.1)。(d) 错:QQHQQ^{\mathrm H} 是把 Cm\mathbb C^m 先投影到 nn 维再放回的算子,与 ImI_m 差一个 m−nm-n 维的正交补——例如 Q=[100]Q=\begin{bmatrix}1\\0\\0\end{bmatrix} 时 QQH=diag⁡(1,0,0)QQ^{\mathrm H}=\operatorname{diag}(1,0,0)。

要点:QHQ=InQ^{\mathrm H}Q=I_n 与 QQH=ImQQ^{\mathrm H}=I_m 只在 m=nm=n(方阵)时可能同时成立——形状就是命题的一部分。

A16 · QR 算法四问

判断对错:

(a) 迭代的每个 AkA_k 都与 AA 相似;(b)「各特征值模值互不相同」是 QR 算法收敛的必要条件;(c) 模互异时,收敛所得上三角阵的对角元按模从大到小排定;(d) 特征值等模时 QR 迭代一定不收敛。

解答 · A16

(a) 对:Ak+1=QkHAkQkA_{k+1}=Q_k^{\mathrm H}A_kQ_k,链起来 Ak=(Q1⋯Qk−1)HA(Q1⋯Qk−1)A_k=(Q_1\cdots Q_{k-1})^{\mathrm H}A(Q_1\cdots Q_{k-1})(QR 算法)。(b) 错——是充分条件,不是必要(例如 AA 已经是对角元全正的上三角阵时迭代一步不动、天然「收敛」,无论对角元模是否互异)。(c) 对(正文的本编注;例 3.4.4 中 4.7321>3>1.26794.7321&gt;3&gt;1.2679 即此序)。(d) 错:等模时可能卡死(旋转阵:A=A⋅I⇒A2=AA=A\cdot I\Rightarrow A_2=A),但「一定不收敛」过强——见 (b) 的三角阵例子,特征值 λ1=λ2\lambda_1=\lambda_2(等模)迭代照样稳定。

要点:背结论要背成「模互异 ⇒\Rightarrow 收敛」,别擅自加强成 ⟺\Longleftrightarrow。

A17 · Schur 引理辨析

判断对错:

(a) 任意复方阵都酉相似于上三角阵;(b) 任意实方阵都正交相似于上三角阵;(c) 特征值全实的实方阵正交相似于上三角阵;(d) 任意复方阵都相似于对角阵。

解答 · A17

(a) 对(定理 3.5.2)。(b) 错——旋转阵 [0−110]\begin{bmatrix}0&amp;-1\\1&amp;0\end{bmatrix} 特征值 ±i\pm\mathrm i:正交相似保持实性,实上三角阵特征值全实,矛盾(3.5.1 的反例分析)。(c) 对(定理 3.5.3,前提正是「特征值均为实数」)。(d) 错(注 1:三角化人人有份,对角化不是——经典三矩阵:同一特征多项式,三种几何重数 2.4.2)。

要点:(a) 与 (d) 的差别就是「对角化」与「三角化」的差别;(b) 与 (c) 的差别就是复结论进入实数域时要补的门票。

A18 · 谱映射与特征向量继承

设 AA 的特征值为 1,−11,-1(各一重,AA 为二阶)。判断:

(a) A2+IA^2+I 的特征值为 2,22,2;(b) AA 的特征向量都是 A2+IA^2+I 的特征向量;(c) A2+IA^2+I 的特征向量都是 AA 的特征向量;(d) A2026A^{2026} 的特征值为 1,11,1。

解答 · A18

(a) 对(推论 3.5.1:φ(λ)=λ2+1\varphi(\lambda)=\lambda^2+1,φ(±1)=2\varphi(\pm1)=2)。(b) 对(注 4 继承:Aα=λα⇒φ(A)α=φ(λ)αA\boldsymbol{\alpha}=\lambda\boldsymbol{\alpha}\Rightarrow\varphi(A)\boldsymbol{\alpha}=\varphi(\lambda)\boldsymbol{\alpha})。(c) 错:A=diag⁡(1,−1)A=\operatorname{diag}(1,-1) 时 A2+I=2IA^2+I=2I,一切非零向量都是它的特征向量,而 e1+e2\boldsymbol{e}_1+\boldsymbol{e}_2 不是 AA 的特征向量——φ\varphi 把两个特征值黏在一起了。(d) 对:φ(λ)=λ2026\varphi(\lambda)=\lambda^{2026},φ(±1)=1\varphi(\pm1)=1(偶次幂把 ±1\pm1 黏成 11,与 (c) 同源)。

A19 · 最小多项式四问

判断对错:

(a) 任何复方阵的最小多项式存在且唯一;(b) 最小多项式恒等于特征多项式;(c) mA(λ)m_A(\lambda) 整除 fA(λ)f_A(\lambda);(d) mAm_A 与 fAf_A 的根相同(不计重数)。

解答 · A19

(a) 对(存在:fAf_A 本身零化(H–C),次数取最小;唯一:互相整除 ++ 同次首一,定理 3.5.5(1))。(b) 错:A=I2A=I_2 时 fA=(λ−1)2f_A=(\lambda-1)^2 而 mA=λ−1m_A=\lambda-1;一般 mA=fAm_A=f_A 当且仅当每个特征值的「最大块」都顶到代数重数(例 3.5.B 的情形)。(c) 对(fAf_A 是零化多项式,mAm_A 整除一切零化多项式)。(d) 对(定理 3.5.5(2):正向由整除,反向由特征向量继承)。

A20 · 逆矩阵公式辨析

设 fA(λ)=λn+an−1λn−1+⋯+a1λ+a0f_A(\lambda)=\lambda^n+a_{n-1}\lambda^{n-1}+\cdots+a_1\lambda+a_0。判断:

(a) a0=0  ⟺  Aa_0=0\iff A 不可逆;(b) AA 可逆时,A−1A^{-1} 可以写成 AA 的多项式;(c) 推论 3.5.2 的公式是 A−1=1a0(An−1+an−1An−2+⋯+a1I)A^{-1}=\frac{1}{a_0}(A^{n-1}+a_{n-1}A^{n-2}+\cdots+a_1I);(d) 例 3.5.3 课件把答案 (1,3)(1,3) 元印成 00,与它前面的特征多项式符号错无关,是独立的抄写错误。

解答 · A20

(a) 对:a0=fA(0)=det⁡(−A)=(−1)ndet⁡Aa_0=f_A(0)=\det(-A)=(-1)^n\det A(推论 3.5.2 证明)。(b) 对——推论 3.5.2 的公式就是一个这样的多项式(这补上了谱映射讨论里「1/λ1/\lambda 不是多项式」的口子)。(c) 错——差一个负号:正确为 −1a0(⋯ )-\frac{1}{a_0}(\cdots)(由 A(⋯ )=−a0IA(\cdots)=-a_0I 移项而来)。(d) 错——三处错误同源连锁:特征多项式符号错((λ−1)3(\lambda-1)^3 印成 (λ+1)3(\lambda+1)^3)⇒\Rightarrow 公式行系数全错 ⇒\Rightarrow 答案矩阵错;即使按印错的公式行计算,答案也不是印出的那个矩阵(见 E3.5 的逐条核算)。

B 组 · 推理证明

B1 · 相抵   ⟺  \iff 同秩

证明:两个 m×nm\times n 矩阵 A,BA,B 相抵(存在可逆 P,QP,Q 使 B=P−1AQB=P^{-1}AQ)当且仅当 rank⁡A=rank⁡B\operatorname{rank}A=\operatorname{rank}B。

解答 · B1

必要性(⟹\Longrightarrow):初等变换不改变秩(本编提示:逐个验证三类行 / 列初等变换,或直接用「相抵 ⇒\Rightarrow 秩相等」的通用论证——P,QP,Q 可逆给出 rank⁡B=rank⁡(P−1AQ)=rank⁡A\operatorname{rank}B=\operatorname{rank}(P^{-1}AQ)=\operatorname{rank}A,因为乘可逆矩阵不改变秩:rank⁡(P−1AQ)=rank⁡(AQ)=rank⁡A\operatorname{rank}(P^{-1}AQ)=\operatorname{rank}(AQ)=\operatorname{rank}A,每一步都用了「乘可逆矩阵不改变秩」。这里 rank⁡(XA)=rank⁡A\operatorname{rank}(XA)=\operatorname{rank}A 对可逆 XX 成立,理由:R(XA)=X(R(A))R(XA)=X(R(A)) 且 XX 是同构(命题 2.3.1(5)),维数不变)。

充分性(⟸\Longleftarrow):设 rank⁡A=rank⁡B=r\operatorname{rank}A=\operatorname{rank}B=r。由 引理 3.1.1,A=P1−1ΓP2A=P_1^{-1}\Gamma P_2、B=Q1−1ΓQ2B=Q_1^{-1}\Gamma Q_2(同一个 Γ\Gamma,因为秩都是 rr)。于是

B=Q1−1ΓQ2=Q1−1(P1AP2−1)Q2=(Q1−1P1) A (P2−1Q2),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),

其中 Q1−1P1Q_1^{-1}P_1、P2−1Q2P_2^{-1}Q_2 都可逆,故 A≅BA\cong B。∎

B2 · 右逆的必要性(并说明左逆如何对称)

证明:若 AB=ImAB=I_m,则 AA 行满秩。

解答 · B2

由 定理 2.3.6,rank⁡(AB)=dim⁡R(AB)\operatorname{rank}(AB)=\dim R(AB);而 R(AB)⊆R(A)R(AB)\subseteq R(A),故

rank⁡(AB)≤rank⁡A.\operatorname{rank}(AB)\le\operatorname{rank}A .

取 B=ImB=I_m 的情形即 rank⁡(AB)=m\operatorname{rank}(AB)=m(ABAB 本身就是 ImI_m,秩为 mm),于是 m≤rank⁡Am\le\operatorname{rank}A;另一方面 AA 是 mm 行矩阵,rank⁡A≤m\operatorname{rank}A\le m。两边夹得 rank⁡A=m\operatorname{rank}A=m,即 AA 行满秩。∎

左逆对称处理:设 BA=InBA=I_n,两边共轭转置得 AHBH=InA^{\mathrm H}B^{\mathrm H}=I_n,即 AHA^{\mathrm H} 有右逆;刚刚证的结论给出 rank⁡(AH)=n\operatorname{rank}(A^{\mathrm H})=n;又 rank⁡A=rank⁡AH\operatorname{rank}A=\operatorname{rank}A^{\mathrm H},故 AA 列满秩。(也可直接对 BABA 用 rank⁡(BA)≤rank⁡A\operatorname{rank}(BA)\le\operatorname{rank}A。)

B3 · 列满秩时的左逆公式

证明:若 A∈Cm×nA\in\mathbb C^{m\times n} 列满秩,则 B=(AHA)−1AHB=(A^{\mathrm H}A)^{-1}A^{\mathrm H} 是 AA 的左逆。

解答 · B3

可逆性:由 命题 3.1.1,rank⁡(AHA)=rank⁡A=n\operatorname{rank}(A^{\mathrm H}A)=\operatorname{rank}A=n,故 nn 阶方阵 AHAA^{\mathrm H}A 满秩、可逆。

验算:

BA=(AHA)−1AHA=In ✓BA=(A^{\mathrm H}A)^{-1}A^{\mathrm H}A=I_n\ \checkmark

要点:这与 定理 3.1.1「法 2」的右逆公式 AH(AAH)−1A^{\mathrm H}(AA^{\mathrm H})^{-1} 完全对称——把 AA 换成 AHA^{\mathrm H} 即得。∎

B4 · 满秩分解中 rank⁡C=r\operatorname{rank}C=r

设 A=BCA=BC 是满秩分解(A∈Crm×nA\in\mathbb C^{m\times n}_r)。证明 rank⁡C=r\operatorname{rank}C=r(写出两步夹逼的完整理由)。

解答 · B4

上界:CC 是 r×nr\times n 矩阵,其秩不超过行数:

rank⁡C≤r.\operatorname{rank}C\le r .

下界:设 Cx=0C\boldsymbol{x}=\boldsymbol{0},两端左乘 BB 得 Ax=BCx=0A\boldsymbol{x}=BC\boldsymbol{x}=\boldsymbol{0},故 N(C)⊆N(A)N(C)\subseteq N(A)。由亏加秩定理(定理 2.2.5,两矩阵列数同为 nn):

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 .

结论:r≤rank⁡C≤rr\le\operatorname{rank}C\le r,故 rank⁡C=r\operatorname{rank}C=r。∎

(BB 的满秩是构造的一部分:BB 的列取成 R(A)R(A) 的基,故 rank⁡B=r\operatorname{rank}B=r 显然。)

B5 · 满秩分解的唯一性(定理 3.2.2)

设 A=B1C1=B2C2A=B_1C_1=B_2C_2 是 A∈Crm×nA\in\mathbb C^{m\times n}_r 的两种满秩分解。证明:存在可逆矩阵 D∈Cr×rD\in\mathbb C^{r\times r} 使 B1=B2DB_1=B_2D、C1=D−1C2C_1=D^{-1}C_2。

提示 · B5

先证两种分解的列空间相同(都要等于 R(A)R(A),其中「R(B1)⊆R(A)R(B_1)\subseteq R(A)」这一步需要 C1C_1 满射);然后把 B1B_1、B2B_2 都看作 R(A)R(A) 的基,用过渡矩阵取 DD;最后用 B2B_2 的列无关把 CC 的关系「约」出来。

解答(本编补,与 [3.2.3](sec32.md#sec-3.2.3) 正文一致)

第一步(列空间相等):C1C_1 行满秩,故有右逆(定理 3.1.1),即 C1C_1 满射。任取 y=B1u∈R(B1)\boldsymbol{y}=B_1\boldsymbol{u}\in R(B_1),取 x\boldsymbol{x} 使 C1x=uC_1\boldsymbol{x}=\boldsymbol{u},则 y=B1C1x=Ax∈R(A)\boldsymbol{y}=B_1C_1\boldsymbol{x}=A\boldsymbol{x}\in R(A)——故 R(B1)⊆R(A)R(B_1)\subseteq R(A);反向 R(A)⊆R(B1)R(A)\subseteq R(B_1) 显然。于是 R(B1)=R(A)R(B_1)=R(A);同理 R(B2)=R(A)R(B_2)=R(A)。

第二步(过渡矩阵):B1B_1、B2B_2 的列都是 rr 维子空间 R(A)R(A) 的基。由过渡矩阵的存在唯一性(定义 1.3.5、命题 1.3.1),存在可逆 DD 使 B1=B2DB_1=B_2D。

第三步(CC 的关系):B2C2=A=B1C1=B2DC1B_2C_2=A=B_1C_1=B_2DC_1,移项得 B2(C2−DC1)=OB_2(C_2-DC_1)=O。B2B_2 列满秩即列无关,B2X=OB_2X=O 只有零解(按列用 定理 1.3.2),故 C2=DC1C_2=DC_1,即 C1=D−1C2C_1=D^{-1}C_2。∎

B6 · 给定 BB 后 CC 唯一

设 A=BCA=BC 是满秩分解。证明:C=(BHB)−1BHAC=(B^{\mathrm H}B)^{-1}B^{\mathrm H}A(即 CC 由 BB 唯一确定)。

解答 · B6

BB 列满秩 ⇒\Rightarrow BHBB^{\mathrm H}B 可逆(命题 3.1.1)。两边左乘 (BHB)−1BH(B^{\mathrm H}B)^{-1}B^{\mathrm H}:

(BHB)−1BHA=(BHB)−1BHBC=(BHB)−1(BHB)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

∎

对偶结论:给定 CC 时 B=ACH(CCH)−1B=AC^{\mathrm H}(CC^{\mathrm H})^{-1}(同样的一句话验算)——见 3.2.2 的本编补注。

B7 · R(A)=R(B)R(A)=R(B)

证明:满秩分解 A=BCA=BC 中 R(A)=R(B)R(A)=R(B)。

提示 · B7

两个包含方向难度不同:一个方向连乘法都不用(BB 的列组合);另一个方向要「解 Cx=uC\boldsymbol{x}=\boldsymbol{u}」,这正是 CC 行满秩的用处。

解答 · B7

R(A)⊆R(B)R(A)\subseteq R(B):任取 y=Ax=B(Cx)\boldsymbol{y}=A\boldsymbol{x}=B(C\boldsymbol{x}),它是 BB 的列的线性组合,故 y∈R(B)\boldsymbol{y}\in R(B)。(不需要任何前提。)

R(B)⊆R(A)R(B)\subseteq R(A):任取 y=Bu∈R(B)\boldsymbol{y}=B\boldsymbol{u}\in R(B)。因 CC 行满秩,由 定理 3.1.1 它有右逆,即 CC 作为映射 Cn→Cr\mathbb C^n\to\mathbb C^r 是满射,故存在 x\boldsymbol{x} 使 Cx=uC\boldsymbol{x}=\boldsymbol{u}。于是

y=Bu=BCx=Ax∈R(A) ✓\boldsymbol{y}=B\boldsymbol{u}=BC\boldsymbol{x}=A\boldsymbol{x}\in R(A)\ \checkmark

故 R(A)=R(B)R(A)=R(B)。∎

推论:dim⁡R(A)=dim⁡R(B)=r\dim R(A)=\dim R(B)=r——满秩分解把列空间「显式地」交给了 BB。∎

B8 · AXA=AAXA=A 与 CXB=IrCXB=I_r 等价

设 A=BC∈Crm×nA=BC\in\mathbb C^{m\times n}_r 是满秩分解。证明:对方阵 X∈Cn×mX\in\mathbb C^{n\times m},

AXA=A  ⟺  CXB=Ir.AXA=A\iff CXB=I_r .

提示 · B8

「⟸\Longleftarrow」直接把 CXBCXB 换成 IrI_r;「⟹\Longrightarrow」分别左乘 BB 的左逆、右乘 CC 的右逆——两条单侧逆的存在性由 定理 3.1.1 与 命题 3.1.1 保证。

解答 · B8

(⟸\Longleftarrow) 若 CXB=IrCXB=I_r,则

AXA=BCXBC=B (CXB) C=BC=A ✓AXA=BCXBC=B\,(CXB)\,C=BC=A\ \checkmark

(⟹\Longrightarrow) 设 BCXBC=BCBCXBC=BC。先左乘 BB 的左逆 BL=(BHB)−1BHB^{\mathrm L}=(B^{\mathrm H}B)^{-1}B^{\mathrm H}:

BLBCXBC=BLBC ⟹ CXBC=C;B^{\mathrm L}BCXBC=B^{\mathrm L}BC\ \Longrightarrow\ CXBC=C ;

再右乘 CC 的右逆 CR=CH(CCH)−1C^{\mathrm R}=C^{\mathrm H}(CC^{\mathrm H})^{-1}:

CXBCCR=CCR ⟹ CXB=Ir ✓CXBCC^{\mathrm R}=CC^{\mathrm R}\ \Longrightarrow\ CXB=I_r\ \checkmark

两个单侧逆都存在(BHBB^{\mathrm H}B、CCHCC^{\mathrm H} 可逆,命题 3.1.1)。∎

B9 · 行满秩 ⟹\Longrightarrow AAHAA^{\mathrm H} 正定

证明:若 A∈Cm×nA\in\mathbb C^{m\times n} 行满秩,则 AAHAA^{\mathrm H} 是正定 Hermite 矩阵(用第一章的正定定义 定义 1.4.3 与 命题 3.1.1 各自完成一部分)。

解答 · B9

Hermite 性:直接算 (AAH)H=A (AH)H=AAH(AA^{\mathrm H})^{\mathrm H}=A\,(A^{\mathrm H})^{\mathrm H}=AA^{\mathrm H} ✓(共轭转置取两次回到自身)。

正定性:对任意非零 x∈Cm\boldsymbol{x}\in\mathbb C^m,

xH(AAH)x=(AHx)H(AHx)=∥AHx∥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 .

等号成立   ⟺  \iff AHx=0A^{\mathrm H}\boldsymbol{x}=\boldsymbol{0}。而 AA 行满秩   ⟺  \iff rank⁡(AH)=m\operatorname{rank}(A^{\mathrm H})=m   ⟺  \iff AHA^{\mathrm H} 列满秩   ⟺  \iff AHA^{\mathrm H} 单射(推论 2.2.2 或 定理 2.2.2),故 AHx=0  ⟺  x=0A^{\mathrm H}\boldsymbol{x}=\boldsymbol{0}\iff\boldsymbol{x}=\boldsymbol{0}。于是对非零 x\boldsymbol{x} 恒有 xHAAHx>0\boldsymbol{x}^{\mathrm H}AA^{\mathrm H}\boldsymbol{x}\gt0,AAHAA^{\mathrm H} 正定。∎

对照:3.1.4 用 命题 3.1.1 说明「AA 行满秩 ⇒\Rightarrow AAHAA^{\mathrm H} 可逆」;本题进一步指出可逆的类型——它是正定 Hermite 矩阵(可逆只是正定的一个推论)。

B10 · 相抵不唯一的「吸收」公式

验证 3.1.1 的不唯一性构造:对 Γ=[Ir000]\Gamma=\begin{bmatrix}I_r&amp;0\\0&amp;0\end{bmatrix} 与任意适当阶数矩阵 ∗*,证明

[Ir∗0I]Γ=Γ,Γ[Ir0∗I]=Γ,\begin{bmatrix}I_r&amp;*\\0&amp;I\end{bmatrix}\Gamma=\Gamma,\qquad \Gamma\begin{bmatrix}I_r&amp;0\\*&amp;I\end{bmatrix}=\Gamma,

并说明这两个三角块都可逆。

解答 · B10

第一式(m×mm\times m 的块左乘):记 Γ=[Ir000]\Gamma=\begin{bmatrix}I_r&amp;0\\0&amp;0\end{bmatrix}(分块尺寸 r+(m−r)r+(m-r)、r+(n−r)r+(n-r))。按分块乘法:

[Ir∗0I][Ir000]=[Ir⋅Ir+∗⋅0 Ir⋅0+∗⋅00⋅Ir+I⋅0 0⋅0+I⋅0]=[Ir000]=Γ ✓\begin{bmatrix}I_r&amp;*\\0&amp;I\end{bmatrix}\begin{bmatrix}I_r&amp;0\\0&amp;0\end{bmatrix}=\begin{bmatrix}I_r\cdot I_r+*\cdot 0&amp;\ I_r\cdot 0+*\cdot 0\\ 0\cdot I_r+I\cdot 0&amp;\ 0\cdot0+I\cdot0\end{bmatrix}=\begin{bmatrix}I_r&amp;0\\0&amp;0\end{bmatrix}=\Gamma\ \checkmark

(关键:Γ\Gamma 的后 n−rn-r 列全零,故 ∗* 乘上去没有贡献。)

第二式(n×nn\times n 的块右乘):Γ\Gamma 的后 m−rm-r 行全零,故

[Ir000][Ir0∗I]=[Ir000]=Γ ✓\begin{bmatrix}I_r&amp;0\\0&amp;0\end{bmatrix}\begin{bmatrix}I_r&amp;0\\*&amp;I\end{bmatrix}=\begin{bmatrix}I_r&amp;0\\0&amp;0\end{bmatrix}=\Gamma\ \checkmark

可逆性:两块都是「单位块 + 严格三角块」,行列式为 1⋅det⁡I=1≠01\cdot\det I=1\neq0,可逆。

用法:由两式得 P−1ΓQ=(P−1[Ir∗0I])Γ([Ir0∗I]Q)P^{-1}\Gamma Q=\big(P^{-1}\begin{bmatrix}I_r&amp;*\\0&amp;I\end{bmatrix}\big)\Gamma\big(\begin{bmatrix}I_r&amp;0\\*&amp;I\end{bmatrix}Q\big),右边仍是合法的相抵分解式——取不同的 ∗* 即得不同的 PP、QQ。∎

B11 · 两种 QR 分解恰好相差一个对角酉阵

设 AA 为 nn 阶满秩复方阵,A=Q1R1=Q2R2A=Q_1R_1=Q_2R_2 是两种 QR 分解(QiQ_i 酉、RiR_i 可逆上三角,不要求正线)。证明:存在对角酉阵 W=diag⁡(w1,⋯ ,wn)W=\operatorname{diag}(w_1,\cdots,w_n)(∣wi∣=1|w_i|=1)使 Q2=Q1WQ_2=Q_1W、R2=WHR1R_2=W^{\mathrm H}R_1。

提示 · B11

考虑 W=R2R1−1=Q2HQ1W=R_2R_1^{-1}=Q_2^{\mathrm H}Q_1,证它「既酉又上三角」;再对列归纳证上三角酉阵必为对角。

解答 · B11

WW 良定且两种身份:R1R_1 可逆(AA 满秩 ⇒rank⁡R1=n\Rightarrow\operatorname{rank}R_1=n,上三角满秩   ⟺  \iff 对角全非零)。由 Q1R1=Q2R2Q_1R_1=Q_2R_2 左乘 Q1HQ_1^{\mathrm H}、右乘 R1−1R_1^{-1}:

W:=Q1HQ2=R1R2−1.W:=Q_1^{\mathrm H}Q_2=R_1R_2^{-1} .

左边是酉阵之积,酉(正交 / 酉矩阵的运算性质 2.5.2);右边是上三角阵之积再求逆,上三角。上三角且酉 ⇒\Rightarrow 对角、∣wii∣=1|w_{ii}|=1(对列归纳:第 11 列只有 w11w_{11} 非零而范数为 11;设前 j−1j-1 列已是 wiieiw_{ii}\boldsymbol{e}_i,第 jj 列与它们的内积为 wii‾wij=0\overline{w_{ii}}w_{ij}=0,逐个逼出 wij=0w_{ij}=0(i<ji&lt;j),范数给出 ∣wjj∣=1|w_{jj}|=1)。于是 WW 是对角酉阵。

收尾:W=Q1HQ2⇒Q2=Q1WW=Q_1^{\mathrm H}Q_2\Rightarrow Q_2=Q_1W(两边左乘 Q1Q_1);W=R1R2−1⇒R2=W−1R1=WHR1W=R_1R_2^{-1}\Rightarrow R_2=W^{-1}R_1=W^{\mathrm H}R_1(酉阵 W−1=WHW^{-1}=W^{\mathrm H}),其第 ii 个对角元为 wi‾ rii(1)\overline{w_i}\,r_{ii}^{(1)}。∎

推论:实数情形 wi=±1w_i=\pm1,恰 2n2^n 种;再要求两个分解都正线(rii(1)>0r_{ii}^{(1)}&gt;0、wi‾rii(1)>0\overline{w_i}r_{ii}^{(1)}&gt;0)则逼出 wi=1w_i=1——唯一性(定理 3.4.1)。

B12 · 上三角酉阵必为对角阵

证明:若 W=(wij)W=(w_{ij}) 既是 nn 阶上三角阵(wij=0, i>jw_{ij}=0,\ i&gt;j)又是酉矩阵,则 wij=0w_{ij}=0(i<ji&lt;j)且 ∣wii∣=1|w_{ii}|=1。

解答 · B12

对列 jj 归纳,断言「第 jj 列 =wjjej=w_{jj}\boldsymbol{e}_j 且 ∣wjj∣=1|w_{jj}|=1」。

j=1j=1:上三角使第一列只有 w11w_{11} 可能非零;酉矩阵各列范数为 11(WHW=IW^{\mathrm H}W=I 的 (1,1)(1,1) 元),故 ∣w11∣=1|w_{11}|=1,第 11 列 =w11e1=w_{11}\boldsymbol{e}_1。

j⇒j+1j\Rightarrow j+1:第 j+1j+1 列的非零分量只在行 1⋯j+11\cdots j+1。对每个 i≤ji\le j,酉矩阵不同列正交(WHW=IW^{\mathrm H}W=I 的 (i,j+1)(i,j+1) 元):

0=∑k≤iwki‾ wk,j+1=wii‾ wi,j+1(归纳假设:第 i 列 =wiiei, 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&gt;i\text{ 项为零}),

而 wii‾≠0\overline{w_{ii}}\neq0,故 wi,j+1=0w_{i,j+1}=0。于是第 j+1j+1 列只剩 wj+1,j+1w_{j+1,j+1},范数 11 给出 ∣wj+1,j+1∣=1|w_{j+1,j+1}|=1。∎

(这条小事实是 定理 3.4.1 唯一性与 B11 的共同引擎;它对「行满 ++ 酉」版本(下三角)同样成立,对 AHA^{\mathrm H} 说一遍即可。)

B13 · QR 迭代的相似链与「卡死」分析

(a) 证明 QR 迭代满足 Ak=(Q1Q2⋯Qk−1)HA(Q1Q2⋯Qk−1)A_k=(Q_1Q_2\cdots Q_{k-1})^{\mathrm H}A(Q_1Q_2\cdots Q_{k-1}),从而每个 AkA_k 与 AA 相似; (b) 证明:旋转阵 A=[0−110]A=\begin{bmatrix}0&amp;-1\\1&amp;0\end{bmatrix} 的迭代序列恒为 Ak=AA_k=A,永不收敛到上三角阵。

解答 · B13

(a) 归纳。k=1k=1:空积为 II,A1=AA_1=A ✓。设 Ak=Wk−1HAWk−1A_k=W_{k-1}^{\mathrm H}AW_{k-1},Wk−1=Q1⋯Qk−1W_{k-1}=Q_1\cdots Q_{k-1}。由迭代式与 Rk=QkHAkR_k=Q_k^{\mathrm H}A_k:

Ak+1=RkQk=QkHAkQk=QkHWk−1HAWk−1Qk=(Wk−1Qk)HA(Wk−1Qk)=WkHAWk,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,

其中 Wk=Wk−1Qk=Q1⋯QkW_k=W_{k-1}Q_k=Q_1\cdots Q_k(酉阵之积仍酉)。故 Ak∼AA_k\sim A,特征值全程不变(相似不变量 2.4.4)。

(b) AA 是正交矩阵,其正线 QR 分解是 A=A⋅IA=A\cdot I:唯一性(定理 3.4.1)保证 GS 只能给出这一组——事实上 GS 对 AA 的两列 e2,e1\boldsymbol{e}_2,\boldsymbol{e}_1 原样通过(已标准正交),Q1=AQ_1=A、R1=IR_1=I。于是

A2=R1Q1=I⋅A=A,A_2=R_1Q_1=I\cdot A=A,

归纳得 Ak=AA_k=A 对一切 kk 成立;而 AA 不是上三角阵,序列不收敛于上三角阵。根源:特征值 ±i\pm\mathrm i 等模,违反收敛的充分条件(QR 算法)。∎

B14 · 重写定理 3.5.1 的归纳证明(不看正文)

不翻正文,独立写出 Schur 引理(复数域版本)的完整证明:任意 nn 阶复方阵相似于上三角阵,且对角元为其特征值。

解答 · B14

对 nn 归纳。n=1n=1 平凡。设对 n−1n-1 成立。AA 的特征多项式在 C\mathbb C 上有根 λ1\lambda_1(代数基本定理),取特征向量 x1≠0\boldsymbol{x}_1\neq\boldsymbol{0}(Ax1=λ1x1A\boldsymbol{x}_1=\lambda_1\boldsymbol{x}_1),扩充为 Cn\mathbb C^n 的基(推论 1.3.1)并拼成可逆阵 P1=[x1,⋯ ]P_1=[\boldsymbol{x}_1,\cdots]。则 P1−1AP1P_1^{-1}AP_1 的第一列为 P1−1Ax1=λ1e1P_1^{-1}A\boldsymbol{x}_1=\lambda_1\boldsymbol{e}_1:

P1−1AP1=[λ1β0A1].P_1^{-1}AP_1=\begin{bmatrix}\lambda_1&amp;\boldsymbol{\beta}\\\boldsymbol{0}&amp;A_1\end{bmatrix}.

对 A1A_1(n−1n-1 阶)用归纳假设:存在可逆 QQ 使 Q−1A1Q=Λ1Q^{-1}A_1Q=\Lambda_1 上三角。令 P=P1diag⁡(1,Q)P=P_1\operatorname{diag}(1,Q),则

P−1AP=[λ1βQ0Λ1]P^{-1}AP=\begin{bmatrix}\lambda_1&amp;\boldsymbol{\beta}Q\\\boldsymbol{0}&amp;\Lambda_1\end{bmatrix}

上三角。对角元:分块三角行列式给 fA(λ)=(λ−λ1)fA1(λ)f_A(\lambda)=(\lambda-\lambda_1)f_{A_1}(\lambda),故对角元全体(计重数)=λ1=\lambda_1 加 A1A_1 的特征值 =A=A 的特征值。∎

自查两处:① 第一步为什么必能取到特征向量?(复数域上特征多项式必有根——这是「复」字在定理里干的活。)② 收尾为什么对角元恰是 AA 的特征值?(相似不改特征多项式 ++ 分块三角行列式。)

B15 · 直接构造酉矩阵(不走定理 3.5.1 的路)

不经过定理 3.5.1,直接用「单位特征向量 ++ Gram–Schmidt 扩充」归纳证明定理 3.5.2:任意复方阵酉相似于上三角阵。

提示 · B15

取单位特征向量 u1\boldsymbol{u}_1,用基扩充 ++ GS 把它扩成标准正交基;[u1,⋯ ][\boldsymbol{u}_1,\cdots] 拼成的矩阵是酉矩阵,算出它的第一列象。

解答 · B15

对 nn 归纳。n=1n=1 平凡。设对 n−1n-1 成立。取 AA 的特征值 λ1\lambda_1 与单位特征向量 u1\boldsymbol{u}_1(∥u1∥=1\|\boldsymbol{u}_1\|=1;特征向量除以范数仍是特征向量)。把 u1\boldsymbol{u}_1 扩充为 Cn\mathbb C^n 的基(推论 1.3.1),再 GS 并单位化(Gram–Schmidt 正交化;「张成不变」保证 u1\boldsymbol{u}_1 原样保留),得标准正交基 u1,u2,⋯ ,un\boldsymbol{u}_1,\boldsymbol{u}_2,\cdots,\boldsymbol{u}_n。令 U1=[u1,⋯ ,un]U_1=[\boldsymbol{u}_1,\cdots,\boldsymbol{u}_n]:列标准正交的方阵,酉(列向量刻画 2.5.1)。计算 U1HAU1U_1^{\mathrm H}AU_1 的第一列:

U1HAu1=λ1U1Hu1=λ1e1,U_1^{\mathrm H}A\boldsymbol{u}_1=\lambda_1U_1^{\mathrm H}\boldsymbol{u}_1=\lambda_1\boldsymbol{e}_1,

故 U1HAU1=[λ1β0A1]U_1^{\mathrm H}AU_1=\begin{bmatrix}\lambda_1&amp;\boldsymbol{\beta}\\\boldsymbol{0}&amp;A_1\end{bmatrix}。对 A1A_1 用归纳假设:存在酉阵 VV 使 VHA1V=Λ1V^{\mathrm H}A_1V=\Lambda_1 上三角。令 U=U1diag⁡(1,V)U=U_1\operatorname{diag}(1,V)(酉阵之积,酉),则

UHAU=diag⁡(1,VH)[λ1β0A1]diag⁡(1,V)=[λ1βV0Λ1]U^{\mathrm H}AU=\operatorname{diag}(1,V^{\mathrm H})\begin{bmatrix}\lambda_1&amp;\boldsymbol{\beta}\\\boldsymbol{0}&amp;A_1\end{bmatrix}\operatorname{diag}(1,V)=\begin{bmatrix}\lambda_1&amp;\boldsymbol{\beta}V\\\boldsymbol{0}&amp;\Lambda_1\end{bmatrix}

上三角。∎

与正文路线的对照:正文(定理 3.5.2)「先相似再 QR 升级」是一步到位的代数戏法;本题是「边归纳边正交化」的施工图。两条路在二阶情形的差别值得一眼:前者需要先有可逆 PP,后者从第一个特征向量起就全程保范。

B16 · 谱映射的一半与反例的另一半

(a) 不用 Schur 分解,证明:若 α\boldsymbol{\alpha} 是 AA 的属于 λ\lambda 的特征向量,则 α\boldsymbol{\alpha} 是 φ(A)\varphi(A) 的属于 φ(λ)\varphi(\lambda) 的特征向量(注 4),且 φ(λ1),⋯ ,φ(λn)\varphi(\lambda_1),\cdots,\varphi(\lambda_n)(计重数)都是 φ(A)\varphi(A) 的特征值; (b) 举反例说明 φ(A)\varphi(A) 的特征向量不必是 AA 的特征向量。

解答 · B16

(a) 由 Aα=λαA\boldsymbol{\alpha}=\lambda\boldsymbol{\alpha} 归纳得 Akα=λkαA^k\boldsymbol{\alpha}=\lambda^k\boldsymbol{\alpha}(k≥0k\ge0,含 A0=IA^0=I)。线性组合:

φ(A)α=∑iaiAiα=(∑iaiλi)α=φ(λ)α,\varphi(A)\boldsymbol{\alpha}=\sum_ia_iA^i\boldsymbol{\alpha}=\Big(\sum_ia_i\lambda^i\Big)\boldsymbol{\alpha}=\varphi(\lambda)\boldsymbol{\alpha},

即 α\boldsymbol{\alpha} 是 φ(A)\varphi(A) 的特征向量(特征值 φ(λ)\varphi(\lambda))。于是对每个特征值 λi\lambda_i 取一个特征向量 αi\boldsymbol{\alpha}_i,φ(λi)\varphi(\lambda_i) 就有证人地成为 φ(A)\varphi(A) 的特征值;重数叙述「计重数」的完整版需数齐 nn 个,那一步用三角化最省(推论 3.5.1 的证明)。

(b) A=diag⁡(1,−1)A=\operatorname{diag}(1,-1)、φ(λ)=λ2\varphi(\lambda)=\lambda^2:φ(A)=A2=I\varphi(A)=A^2=I,任何非零向量都是 II 的特征向量(特征值 11),但 e1+e2\boldsymbol{e}_1+\boldsymbol{e}_2 不是 AA 的特征向量。∎

要点:继承是单向的——φ\varphi 可以把不同的特征值黏到同一个值上,信息只减不增(这正是最小多项式要「计账」的对象)。

B17 · 重写定理 3.5.5(不看正文)

独立证明最小多项式的两条性质:(1) mAm_A 唯一、整除 AA 的任一零化多项式(特别地 mA∣fAm_A\mid f_A);(2) fAf_A 与 mAm_A 同根(不计重数)。

解答 · B17

(1) 存在性(备查):fAf_A 零化 AA(H–C,定理 3.5.4),故零化多项式之族非空,取次数最小者再除以首项系数(仍零化、次数不变)得首一的最小者 mAm_A。

整除:设 gg 任一零化多项式。带余除法(定理 2.1.3):g=qmA+rg=qm_A+r,deg⁡r<deg⁡mA\deg r&lt;\deg m_A 或 r=0r=0。代入矩阵并注意同一矩阵的多项式可交换:

O=g(A)=q(A)mA(A)⏟=O+r(A) ⟹ r(A)=O.O=g(A)=q(A)\underbrace{m_A(A)}_{=O}+r(A)\ \Longrightarrow\ r(A)=O .

若 r≠0r\neq0,则 rr(除以首项系数后首一)是次数更低的零化多项式,与最小性矛盾。故 r=0r=0,mA∣gm_A\mid g;取 g=fAg=f_A 得 mA∣fAm_A\mid f_A。

唯一:若 m1,m2m_1,m_2 都是最小的首一零化多项式,则 m1∣m2m_1\mid m_2 且 m2∣m1m_2\mid m_1;两者同次(都最小)且首一,故 m1=m2m_1=m_2(互相整除的同次首一多项式只能相差 11 倍)。

(2) 同根:「mAm_A 的根 ⊆fA\subseteq f_A 的根」:mA∣fAm_A\mid f_A,根只会少不会多。「fAf_A 的根 ⊆mA\subseteq m_A 的根」:设 λ0\lambda_0 是 AA 的特征值,取特征向量 α≠0\boldsymbol{\alpha}\neq\boldsymbol{0}。由注 4(B16 的 (a)):

0=mA(A)α=mA(λ0)α ⟹ mA(λ0)=0.■\boldsymbol{0}=m_A(A)\boldsymbol{\alpha}=m_A(\lambda_0)\boldsymbol{\alpha}\ \Longrightarrow\ m_A(\lambda_0)=0 .\qquad\blacksquare

一处易错:不能用「fA(A)=O⇒fA=mAf_A(A)=O\Rightarrow f_A=m_A」——零化多项式有无数个(注 5),fAf_A 只是其中之一;同根但重数可以差(I2I_2:fA=(λ−1)2f_A=(\lambda-1)^2,mA=λ−1m_A=\lambda-1)。

C 组 · 计算实操

C1 · 相抵分解全流程(2×32\times3、秩 11)

求 A=[123246]A=\begin{bmatrix}1&amp;2&amp;3\\ 2&amp;4&amp;6\end{bmatrix} 的相抵分解式:写出 Γ\Gamma、EE、FF(即 PP 与 Q−1Q^{-1}),再写出 P−1P^{-1}、QQ 并回代验证。

解答 · C1

第一步(定秩):第二行 =2×=2\times 第一行,故 r=1r=1,Γ=[100000]\Gamma=\begin{bmatrix}1&amp;0&amp;0\\ 0&amp;0&amp;0\end{bmatrix}。

第二步(建分块矩阵并消元):挂上 [AI2I3∗]\begin{bmatrix}A&amp;I_2\\ I_3&amp;*\end{bmatrix},只动前 22 行、前 33 列。

  • 行变换(记总积 EE):r2←r2−2r1r_2\leftarrow r_2-2r_1,即 E=[10−21]E=\begin{bmatrix}1&amp;0\\ -2&amp;1\end{bmatrix};
  • 列变换(记总积 FF):c2←c2−2c1c_2\leftarrow c_2-2c_1、c3←c3−3c1c_3\leftarrow c_3-3c_1,即 F=[1−2−3010001]F=\begin{bmatrix}1&amp;-2&amp;-3\\ 0&amp;1&amp;0\\ 0&amp;0&amp;1\end{bmatrix}。

第三步(回代验算,用 EAF=ΓEAF=\Gamma):

EAF=[10−21][123246][1−2−3010001]=[123000][1−2−3010001]=[100000]=Γ ✓EAF=\begin{bmatrix}1&amp;0\\-2&amp;1\end{bmatrix}\begin{bmatrix}1&amp;2&amp;3\\2&amp;4&amp;6\end{bmatrix}\begin{bmatrix}1&amp;-2&amp;-3\\0&amp;1&amp;0\\0&amp;0&amp;1\end{bmatrix} =\begin{bmatrix}1&amp;2&amp;3\\0&amp;0&amp;0\end{bmatrix}\begin{bmatrix}1&amp;-2&amp;-3\\0&amp;1&amp;0\\0&amp;0&amp;1\end{bmatrix}=\begin{bmatrix}1&amp;0&amp;0\\0&amp;0&amp;0\end{bmatrix}=\Gamma\ \checkmark

第四步(写成 A=P−1ΓQA=P^{-1}\Gamma Q):由 EAF=ΓEAF=\Gamma 得 A=E−1ΓF−1A=E^{-1}\Gamma F^{-1},与 A=P−1ΓQA=P^{-1}\Gamma Q 对照(相抵变换求 PP 与 Q−1Q^{-1}):

 P−1=E−1=[1021],Q=F−1=[123010001],P=E,Q−1=F \boxed{\,P^{-1}=E^{-1}=\begin{bmatrix}1&amp;0\\ 2&amp;1\end{bmatrix},\qquad Q=F^{-1}=\begin{bmatrix}1&amp;2&amp;3\\ 0&amp;1&amp;0\\ 0&amp;0&amp;1\end{bmatrix},\qquad P=E,\qquad Q^{-1}=F\,}

最终验证:

P−1ΓQ=[1021][100000][123010001]=[100200][123010001]=[123246]=A ✓P^{-1}\Gamma Q=\begin{bmatrix}1&amp;0\\2&amp;1\end{bmatrix}\begin{bmatrix}1&amp;0&amp;0\\0&amp;0&amp;0\end{bmatrix}\begin{bmatrix}1&amp;2&amp;3\\0&amp;1&amp;0\\0&amp;0&amp;1\end{bmatrix} =\begin{bmatrix}1&amp;0&amp;0\\2&amp;0&amp;0\end{bmatrix}\begin{bmatrix}1&amp;2&amp;3\\0&amp;1&amp;0\\0&amp;0&amp;1\end{bmatrix}=\begin{bmatrix}1&amp;2&amp;3\\2&amp;4&amp;6\end{bmatrix}=A\ \checkmark

要点:读出的两块是「变换的总积」EE、FF,满足 EAF=ΓEAF=\Gamma;写分解式时要取逆(P−1=E−1P^{-1}=E^{-1}、Q=F−1Q=F^{-1})。用 EAF=ΓEAF=\Gamma 验算比用 P−1ΓQP^{-1}\Gamma Q 少两步矩阵乘法,也不容易把栏目记错。∎

C2 · 用相抵标准形解方程组

用 C1 的分解解 Ax=bA\boldsymbol{x}=\boldsymbol{b}(A=[123246]A=\begin{bmatrix}1&amp;2&amp;3\\2&amp;4&amp;6\end{bmatrix},b=(1,2)T\boldsymbol{b}=(1,2)^{\mathrm T}),写出通解;再判断 b′=(1,3)T\boldsymbol{b}&#x27;=(1,3)^{\mathrm T} 是否有解。

解答 · C2

取 E=[10−21]E=\begin{bmatrix}1&amp;0\\-2&amp;1\end{bmatrix}(EAF=ΓEAF=\Gamma)。方程 Γy=Eb\Gamma\boldsymbol{y}=E\boldsymbol{b},x=Fy\boldsymbol{x}=F\boldsymbol{y}(因为 A=E−1ΓF−1A=E^{-1}\Gamma F^{-1},故 Γ(F−1x)=Eb\Gamma(F^{-1}\boldsymbol{x})=E\boldsymbol{b},记 y=F−1x\boldsymbol{y}=F^{-1}\boldsymbol{x})。

(a) Eb=[10−21][12]=[10]E\boldsymbol{b}=\begin{bmatrix}1&amp;0\\-2&amp;1\end{bmatrix}\begin{bmatrix}1\\2\end{bmatrix}=\begin{bmatrix}1\\0\end{bmatrix},后 m−r=1m-r=1 个分量为 00 ✓ 有解。由 Γy=(y1,0)T\Gamma\boldsymbol{y}=(y_1,0)^{\mathrm T} 得 y1=1y_1=1、y2,y3y_2,y_3 任意(记 y2=s,y3=ty_2=s,y_3=t):

x=Fy=[1−2−3010001][1st]=[1−2s−3tst].\boldsymbol{x}=F\boldsymbol{y}=\begin{bmatrix}1&amp;-2&amp;-3\\0&amp;1&amp;0\\0&amp;0&amp;1\end{bmatrix}\begin{bmatrix}1\\s\\t\end{bmatrix}=\begin{bmatrix}1-2s-3t\\ s\\ t\end{bmatrix}.

验算:Ax=(1−2s−3t+2s+3t, 2−4s−6t+4s+6t)T=(1,2)T=bA\boldsymbol{x}=(1-2s-3t+2s+3t,\ 2-4s-6t+4s+6t)^{\mathrm T}=(1,2)^{\mathrm T}=\boldsymbol{b} ✓;自由参数 2=n−r2=n-r ✓。

(b) Eb′=[1−2+3]=[11]E\boldsymbol{b}&#x27;=\begin{bmatrix}1\\-2+3\end{bmatrix}=\begin{bmatrix}1\\1\end{bmatrix},第二个分量 ≠0\neq0,无解 ✓(与经典判据一致:rank⁡(A∣b′)=2>rank⁡A=1\operatorname{rank}(A|\boldsymbol{b}&#x27;)=2&gt;\operatorname{rank}A=1,因为 b′\boldsymbol{b}&#x27; 不满足「第二分量 =2×=2\times 第一分量」)。

C3 · 相抵分解重做(例 3.1.A)

求 A=[1232573710]A=\begin{bmatrix}1&amp;2&amp;3\\ 2&amp;5&amp;7\\ 3&amp;7&amp;10\end{bmatrix} 的相抵分解式(Γ\Gamma、EE、FF,并写出 A=P−1ΓQA=P^{-1}\Gamma Q)。

解答 · C3

过程见 例 3.1.A(分行、列两步;r=2r=2):

Γ=[100010000],E=[100−210−1−11],F=[1−2−101−1001],EAF=Γ ✓\Gamma=\begin{bmatrix}1&amp;0&amp;0\\ 0&amp;1&amp;0\\ 0&amp;0&amp;0\end{bmatrix},\qquad E=\begin{bmatrix}1&amp;0&amp;0\\ -2&amp;1&amp;0\\ -1&amp;-1&amp;1\end{bmatrix},\qquad F=\begin{bmatrix}1&amp;-2&amp;-1\\ 0&amp;1&amp;-1\\ 0&amp;0&amp;1\end{bmatrix},\qquad EAF=\Gamma\ \checkmark

A=E−1ΓF−1,E−1=[100210311],F−1=[123011001].A=E^{-1}\Gamma F^{-1},\qquad E^{-1}=\begin{bmatrix}1&amp;0&amp;0\\ 2&amp;1&amp;0\\ 3&amp;1&amp;1\end{bmatrix},\qquad F^{-1}=\begin{bmatrix}1&amp;2&amp;3\\ 0&amp;1&amp;1\\ 0&amp;0&amp;1\end{bmatrix}.

验算:E−1ΓF−1=[100210310][123011001]=[1232573710]=AE^{-1}\Gamma F^{-1}=\begin{bmatrix}1&amp;0&amp;0\\2&amp;1&amp;0\\3&amp;1&amp;0\end{bmatrix}\begin{bmatrix}1&amp;2&amp;3\\0&amp;1&amp;1\\0&amp;0&amp;1\end{bmatrix}=\begin{bmatrix}1&amp;2&amp;3\\2&amp;5&amp;7\\3&amp;7&amp;10\end{bmatrix}=A ✓。∎

C4 · 全部右逆

求 A=[120011]A=\begin{bmatrix}1&amp;2&amp;0\\ 0&amp;1&amp;1\end{bmatrix} 的全部右逆。

提示 · C4

先用公式 B0=AH(AAH)−1B_0=A^{\mathrm H}(AA^{\mathrm H})^{-1} 造一个,再考虑「差」:A(B−B0)=OA(B-B_0)=O,即 B−B0B-B_0 的每一列都属于 N(A)N(A)——把 N(A)N(A) 的基求出来后加上参数即可。

解答 · C4

基本解:AAH=[5222]AA^{\mathrm H}=\begin{bmatrix}5&amp;2\\2&amp;2\end{bmatrix},(AAH)−1=16[2−2−25](AA^{\mathrm H})^{-1}=\dfrac16\begin{bmatrix}2&amp;-2\\-2&amp;5\end{bmatrix},

B0=AH(AAH)−1=16[2−221−25],AB0=I2 ✓B_0=A^{\mathrm H}(AA^{\mathrm H})^{-1}=\dfrac16\begin{bmatrix}2&amp;-2\\ 2&amp;1\\ -2&amp;5\end{bmatrix},\qquad AB_0=I_2\ \checkmark

通式:设 B=B0+ZB=B_0+Z,则 AB=B0AB=B_0 与 A=IA=I 的关系为 AB=A(B0+Z)=I2+AZAB=A(B_0+Z)=I_2+AZ;要 AB=I2AB=I_2 只需 AZ=OAZ=O,即 ZZ 的每一列属于 N(A)N(A)。解 Ax=0A\boldsymbol{x}=\boldsymbol{0}:x1+2x2=0x_1+2x_2=0、x2+x3=0x_2+x_3=0,取 x3=1x_3=1 得 n=(2,−1,1)T\boldsymbol{n}=(2,-1,1)^{\mathrm T},N(A)=span⁡{n}N(A)=\operatorname{span}\{\boldsymbol{n}\}(dim⁡=n−r=1\dim=n-r=1 ✓)。

于是全部右逆为

B=16[2−221−25]+[2s2t−s−tst],s,t∈CB=\dfrac16\begin{bmatrix}2&amp;-2\\ 2&amp;1\\ -2&amp;5\end{bmatrix}+\begin{bmatrix}2s&amp;2t\\ -s&amp;-t\\ s&amp;t\end{bmatrix},\qquad s,t\in\mathbb C

(每个参数配一列)。验算:A[2s2t−s−tst]=[2s−2s2t−2t−s+s−t+t]=OA\begin{bmatrix}2s&amp;2t\\-s&amp;-t\\ s&amp;t\end{bmatrix}=\begin{bmatrix}2s-2s&amp;2t-2t\\ -s+s&amp;-t+t\end{bmatrix}=O ✓,故 AB=I2AB=I_2 ✓。两个自由参数 == n−rn-r(每列一个)×\times mm(m=2m=2 列)÷\div ⋯——不背计数公式,只记住「多余的自由度来自 N(A)N(A)」。

C5 · 判断并构造单侧逆

判断下列矩阵是否有右逆 / 左逆;对存在的,用公式构造一个:

(a) A1=[100010]A_1=\begin{bmatrix}1&amp;0&amp;0\\ 0&amp;1&amp;0\end{bmatrix};(b) A2=[100100]A_2=\begin{bmatrix}1&amp;0\\ 0&amp;1\\ 0&amp;0\end{bmatrix};(c) A3=[1224]A_3=\begin{bmatrix}1&amp;2\\ 2&amp;4\end{bmatrix}。

解答 · C5

(a) 行满秩(r=2=mr=2=m):有右逆。A1A1H=I2A_1A_1^{\mathrm H}=I_2(A1A_1 的行已标准正交),故 B=A1H(A1A1H)−1=A1H=[100100]B=A_1^{\mathrm H}(A_1A_1^{\mathrm H})^{-1}=A_1^{\mathrm H}=\begin{bmatrix}1&amp;0\\0&amp;1\\0&amp;0\end{bmatrix};验算 A1B=I2A_1B=I_2 ✓。

(b) 列满秩(r=2=nr=2=n):有左逆。A2HA2=I2A_2^{\mathrm H}A_2=I_2,故 B=(A2HA2)−1A2H=A2H=[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],即 [100010]\begin{bmatrix}1&amp;0&amp;0\\0&amp;1&amp;0\end{bmatrix};验算 BA2=I2BA_2=I_2 ✓。

(c) 秩 1<21&lt;2:两者都没有(若硬算会发现 A3HA3=[5101020]A_3^{\mathrm H}A_3=\begin{bmatrix}5&amp;10\\10&amp;20\end{bmatrix} 不可逆,公式失效——这与 定理 3.1.1 的判断一致)。

要点:(a)(b) 互为「转置镜像」,两个公式也互为镜像。

C6 · 右逆公式实操

设 A=[110011]A=\begin{bmatrix}1&amp;1&amp;0\\ 0&amp;1&amp;1\end{bmatrix}:判断 AA 是否有右逆;用公式构造一个右逆并验算;说明为什么不用相抵分解也能得到答案。

解答 · C6

判断:2×32\times3、两行不成比例 ⇒\Rightarrow rank⁡A=2=m\operatorname{rank}A=2=m,行满秩,有右逆。

构造:

AAH=[110011][101101]=[2112],(AAH)−1=13[2−1−12],AA^{\mathrm H}=\begin{bmatrix}1&amp;1&amp;0\\0&amp;1&amp;1\end{bmatrix}\begin{bmatrix}1&amp;0\\1&amp;1\\0&amp;1\end{bmatrix}=\begin{bmatrix}2&amp;1\\1&amp;2\end{bmatrix},\qquad (AA^{\mathrm H})^{-1}=\frac13\begin{bmatrix}2&amp;-1\\-1&amp;2\end{bmatrix},

B=AH(AAH)−1=13[101101][2−1−12]=13[2−111−12].B=A^{\mathrm H}(AA^{\mathrm H})^{-1}=\frac13\begin{bmatrix}1&amp;0\\1&amp;1\\0&amp;1\end{bmatrix}\begin{bmatrix}2&amp;-1\\-1&amp;2\end{bmatrix}=\frac13\begin{bmatrix}2&amp;-1\\1&amp;1\\-1&amp;2\end{bmatrix}.

验算:AB=13[110011][2−111−12]=13[3003]=I2AB=\dfrac13\begin{bmatrix}1&amp;1&amp;0\\0&amp;1&amp;1\end{bmatrix}\begin{bmatrix}2&amp;-1\\1&amp;1\\-1&amp;2\end{bmatrix}=\dfrac13\begin{bmatrix}3&amp;0\\0&amp;3\end{bmatrix}=I_2 ✓。

为什么不用相抵分解:法 2 的公式依赖 命题 3.1.1(rank⁡(AAH)=rank⁡A=m\operatorname{rank}(AA^{\mathrm H})=\operatorname{rank}A=m,故 AAHAA^{\mathrm H} 可逆),而这一命题的证明只用方程组同解 + 内积正定性,与相抵分解无关——所以公式法比相抵法短得多。

C7 · 满秩分解(取 AA 的列)

求 A=[123011]A=\begin{bmatrix}1&amp;2&amp;3\\ 0&amp;1&amp;1\end{bmatrix} 的满秩分解,并写出 R(A)R(A) 的一组基。

解答 · C7

定基:a1=(1,0)T\boldsymbol{a}_1=(1,0)^{\mathrm T}、a2=(2,1)T\boldsymbol{a}_2=(2,1)^{\mathrm T}、a3=(3,1)T\boldsymbol{a}_3=(3,1)^{\mathrm T}。a1,a2\boldsymbol{a}_1,\boldsymbol{a}_2 无关;a3=a1+a2\boldsymbol{a}_3=\boldsymbol{a}_1+\boldsymbol{a}_2 ✓。r=rank⁡A=2r=\operatorname{rank}A=2,故 a1,a2\boldsymbol{a}_1,\boldsymbol{a}_2 是 R(A)R(A) 的基。

满秩分解:

B=[a1,a2]=[1201],C=[101011](坐标:a3=a1+a2),B=[\boldsymbol{a}_1,\boldsymbol{a}_2]=\begin{bmatrix}1&amp;2\\ 0&amp;1\end{bmatrix},\qquad C=\begin{bmatrix}1&amp;0&amp;1\\ 0&amp;1&amp;1\end{bmatrix}\quad(\text{坐标:}\boldsymbol{a}_3=\boldsymbol{a}_1+\boldsymbol{a}_2),

BC=[1201][101011]=[123011]=A ✓BC=\begin{bmatrix}1&amp;2\\ 0&amp;1\end{bmatrix}\begin{bmatrix}1&amp;0&amp;1\\ 0&amp;1&amp;1\end{bmatrix}=\begin{bmatrix}1&amp;2&amp;3\\ 0&amp;1&amp;1\end{bmatrix}=A\ \checkmark

R(A)=span⁡{(1,0)T,(2,1)T}R(A)=\operatorname{span}\{(1,0)^{\mathrm T},(2,1)^{\mathrm T}\}。∎

C8 · 满秩分解(RREF 法)

用 用 RREF 批量求满秩分解 求 A=[112123521231]A=\begin{bmatrix}1&amp;1&amp;2&amp;1\\ 2&amp;3&amp;5&amp;2\\ 1&amp;2&amp;3&amp;1\end{bmatrix} 的满秩分解,并读出 AA 的列之间的线性关系。

解答 · C8

化 RREF:r2−2r1=(0,1,1,0)r_2-2r_1=(0,1,1,0);r3−r2=(0,0,0,0)r_3-r_2=(0,0,0,0);回代 r1−r2=(1,0,1,1)r_1-r_2=(1,0,1,1)。得

R=[101101100000],r=2, 主元列=1,2.R=\begin{bmatrix}1&amp;0&amp;1&amp;1\\ 0&amp;1&amp;1&amp;0\\ 0&amp;0&amp;0&amp;0\end{bmatrix},\qquad r=2,\ \text{主元列}=1,2 .

取 BB、CC:B=AB=A 的第 1、2 列 =[112312]=\begin{bmatrix}1&amp;1\\2&amp;3\\1&amp;2\end{bmatrix},C=RC=R 的前两行 =[10110110]=\begin{bmatrix}1&amp;0&amp;1&amp;1\\ 0&amp;1&amp;1&amp;0\end{bmatrix}。

验算:[112312][10110110]=[112123521231]=A\begin{bmatrix}1&amp;1\\2&amp;3\\1&amp;2\end{bmatrix}\begin{bmatrix}1&amp;0&amp;1&amp;1\\0&amp;1&amp;1&amp;0\end{bmatrix}=\begin{bmatrix}1&amp;1&amp;2&amp;1\\2&amp;3&amp;5&amp;2\\1&amp;2&amp;3&amp;1\end{bmatrix}=A ✓。

列关系:CC 的第三列 (1,1)T(1,1)^{\mathrm T} ⟹ a3=a1+a2\boldsymbol{a}_3=\boldsymbol{a}_1+\boldsymbol{a}_2;第四列 (1,0)T(1,0)^{\mathrm T} ⟹ a4=a1\boldsymbol{a}_4=\boldsymbol{a}_1。代入验算:a3=(2,5,3)\boldsymbol{a}_3=(2,5,3) 而 a1+a2=(1,1)+(1,2)=(2,3,1)\boldsymbol{a}_1+\boldsymbol{a}_2=(1,1)+(1,2)=(2,3,1)——不对!查 AA:a3=(2,5,3)T\boldsymbol{a}_3=(2,5,3)^{\mathrm T},a1+a2=(1,1,1)+(2,3,2)=(3,4,3)\boldsymbol{a}_1+\boldsymbol{a}_2=(1,1,1)+(2,3,2)=(3,4,3),都不等于。重新核对 RR 的行 r1=(1,0,1,1)r_1=(1,0,1,1)、r2=(0,1,1,0)r_2=(0,1,1,0),且 BB 取的是 AA 的第 1、2 列:a1=(1,2,1)\boldsymbol{a}_1=(1,2,1)、a2=(1,3,2)\boldsymbol{a}_2=(1,3,2)。则列关系应为 a3=1⋅a1+1⋅a2=(2,5,3)\boldsymbol{a}_3=1\cdot\boldsymbol{a}_1+1\cdot\boldsymbol{a}_2=(2,5,3) ✓(正是 AA 的第三列);a4=1⋅a1+0⋅a2=(1,2,1)\boldsymbol{a}_4=1\cdot\boldsymbol{a}_1+0\cdot\boldsymbol{a}_2=(1,2,1) ✓。上面「不对」的检查把 BB 的列抄着用了 AA 的行,属笔误——教训:BB 是 AA 的列拼成的矩阵,核对时要用列向量。∎

C9 · 两种满秩分解与过渡矩阵

对 A=[123011]A=\begin{bmatrix}1&amp;2&amp;3\\ 0&amp;1&amp;1\end{bmatrix}:写出两种不同的满秩分解 A=B1C1=B2C2A=B_1C_1=B_2C_2,求出 定理 3.2.2 中的 DD,并验证 B1=B2DB_1=B_2D、C1=D−1C2C_1=D^{-1}C_2。

解答 · C9

分解一(取第 1、2 列):B1=[1201]B_1=\begin{bmatrix}1&amp;2\\ 0&amp;1\end{bmatrix},C1=[101011]C_1=\begin{bmatrix}1&amp;0&amp;1\\ 0&amp;1&amp;1\end{bmatrix}(B1B_1 可逆,C1=B1−1AC_1=B_1^{-1}A)。

分解二(取第 1、3 列):B2=[1301]B_2=\begin{bmatrix}1&amp;3\\ 0&amp;1\end{bmatrix},C2=[1−10011]C_2=\begin{bmatrix}1&amp;-1&amp;0\\ 0&amp;1&amp;1\end{bmatrix}。验算:[1301][1−10011]=[123011]=A\begin{bmatrix}1&amp;3\\0&amp;1\end{bmatrix}\begin{bmatrix}1&amp;-1&amp;0\\0&amp;1&amp;1\end{bmatrix}=\begin{bmatrix}1&amp;2&amp;3\\0&amp;1&amp;1\end{bmatrix}=A ✓。

求 DD:由 B1=B2DB_1=B_2D 得

D=B2−1B1=[1−301][1201]=[1−101],det⁡D=1≠0 ✓D=B_2^{-1}B_1=\begin{bmatrix}1&amp;-3\\ 0&amp;1\end{bmatrix}\begin{bmatrix}1&amp;2\\ 0&amp;1\end{bmatrix}=\begin{bmatrix}1&amp;-1\\ 0&amp;1\end{bmatrix},\qquad \det D=1\neq0\ \checkmark

验算 C1=D−1C2C_1=D^{-1}C_2:D−1=[1101]D^{-1}=\begin{bmatrix}1&amp;1\\ 0&amp;1\end{bmatrix},

D−1C2=[1101][1−10011]=[101011]=C1 ✓D^{-1}C_2=\begin{bmatrix}1&amp;1\\ 0&amp;1\end{bmatrix}\begin{bmatrix}1&amp;-1&amp;0\\ 0&amp;1&amp;1\end{bmatrix}=\begin{bmatrix}1&amp;0&amp;1\\ 0&amp;1&amp;1\end{bmatrix}=C_1\ \checkmark

读 DD:DD 是「由 B2B_2 的列到 B1B_1 的列」的过渡矩阵——即 ([a1,a3])=([a1,a2])D([\boldsymbol{a}_1,\boldsymbol{a}_3])=([\boldsymbol{a}_1,\boldsymbol{a}_2])D,与 注 3 一致。∎

C10 · 相抵分解(3×23\times2、秩 11)

求 A=[1−1002−2]A=\begin{bmatrix}1&amp;-1\\ 0&amp;0\\ 2&amp;-2\end{bmatrix} 的相抵分解式。

解答 · C10

秩:r=1r=1(三行都与 (1,−1)(1,-1) 成比例;实际 r3=2r1\boldsymbol{r}_3=2\boldsymbol{r}_1、r2=0\boldsymbol{r}_2=\boldsymbol{0})。

消元:r3←r3−2r1r_3\leftarrow r_3-2r_1;c2←c2+c1c_2\leftarrow c_2+c_1。得 E=[100010−201]E=\begin{bmatrix}1&amp;0&amp;0\\ 0&amp;1&amp;0\\ -2&amp;0&amp;1\end{bmatrix}、F=[1101]F=\begin{bmatrix}1&amp;1\\ 0&amp;1\end{bmatrix},EAF=[100000]=ΓEAF=\begin{bmatrix}1&amp;0\\0&amp;0\\0&amp;0\end{bmatrix}=\Gamma ✓。

分解式:A=E−1ΓF−1A=E^{-1}\Gamma F^{-1},

E−1=[100010201],F−1=[1−101],E−1ΓF−1=[1−1002−2]=A ✓E^{-1}=\begin{bmatrix}1&amp;0&amp;0\\ 0&amp;1&amp;0\\ 2&amp;0&amp;1\end{bmatrix},\qquad F^{-1}=\begin{bmatrix}1&amp;-1\\ 0&amp;1\end{bmatrix},\qquad E^{-1}\Gamma F^{-1}=\begin{bmatrix}1&amp;-1\\ 0&amp;0\\ 2&amp;-2\end{bmatrix}=A\ \checkmark

∎

C11 · 四个秩的核对

对 A=[122436]A=\begin{bmatrix}1&amp;2\\ 2&amp;4\\ 3&amp;6\end{bmatrix}:计算 AHAA^{\mathrm H}A、AAHAA^{\mathrm H},并核对 命题 3.1.1 的四个秩。

解答 · C11

AHA=[123246][122436]=[14282856],AAH=[122436][123246]=[51015102030153045].A^{\mathrm H}A=\begin{bmatrix}1&amp;2&amp;3\\2&amp;4&amp;6\end{bmatrix}\begin{bmatrix}1&amp;2\\2&amp;4\\3&amp;6\end{bmatrix}=\begin{bmatrix}14&amp;28\\ 28&amp;56\end{bmatrix},\qquad AA^{\mathrm H}=\begin{bmatrix}1&amp;2\\2&amp;4\\3&amp;6\end{bmatrix}\begin{bmatrix}1&amp;2&amp;3\\2&amp;4&amp;6\end{bmatrix}=\begin{bmatrix}5&amp;10&amp;15\\ 10&amp;20&amp;30\\ 15&amp;30&amp;45\end{bmatrix}.

秩:rank⁡A=1\operatorname{rank}A=1(列成比例);AHAA^{\mathrm H}A 第二行 =2×=2\times 第一行 ⇒\Rightarrow 秩 11;AAHAA^{\mathrm H} 的各行成比例 ⇒\Rightarrow 秩 11;rank⁡(AH)=1\operatorname{rank}(A^{\mathrm H})=1。四项相等 ✓:

rank⁡(A)=rank⁡(AHA)=rank⁡(AAH)=rank⁡(AH)=1.\operatorname{rank}(A)=\operatorname{rank}(A^{\mathrm H}A)=\operatorname{rank}(AA^{\mathrm H})=\operatorname{rank}(A^{\mathrm H})=1 .

要点:AHAA^{\mathrm H}A、AAHAA^{\mathrm H} 与 AA 同秩,但尺寸不同(AHAA^{\mathrm H}A 是 2×22\times2、AAHAA^{\mathrm H} 是 3×33\times3)——秩相同、零空间不同(N(AHA)=N(A)N(A^{\mathrm H}A)=N(A),而 N(AAH)=N(AH)N(AA^{\mathrm H})=N(A^{\mathrm H}))。∎

C12 · 满秩分解与四个子空间

对 A=[101011112]A=\begin{bmatrix}1&amp;0&amp;1\\ 0&amp;1&amp;1\\ 1&amp;1&amp;2\end{bmatrix}:求满秩分解;由此写出 R(A)R(A) 的基、rank⁡A\operatorname{rank}A 与 N(A)N(A) 的基。

解答 · C12

满秩分解:a3=a1+a2\boldsymbol{a}_3=\boldsymbol{a}_1+\boldsymbol{a}_2,r=2r=2,取第 1、2 列:

B=[100111],C=[101011],BC=[101011112]=A ✓B=\begin{bmatrix}1&amp;0\\ 0&amp;1\\ 1&amp;1\end{bmatrix},\qquad C=\begin{bmatrix}1&amp;0&amp;1\\ 0&amp;1&amp;1\end{bmatrix},\qquad BC=\begin{bmatrix}1&amp;0&amp;1\\0&amp;1&amp;1\\1&amp;1&amp;2\end{bmatrix}=A\ \checkmark

R(A)R(A) 的基:BB 的两列,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}\},rank⁡A=2\operatorname{rank}A=2。

N(A)N(A) 的基:解 Ax=0A\boldsymbol{x}=\boldsymbol{0}:x1+x3=0x_1+x_3=0、x2+x3=0x_2+x_3=0,取 x3=1x_3=1:n=(−1,−1,1)T\boldsymbol{n}=(-1,-1,1)^{\mathrm T}。验算 An=(−1+1,−1+1,−1−1+2)T=0A\boldsymbol{n}=(-1+1,-1+1,-1-1+2)^{\mathrm T}=\boldsymbol{0} ✓。N(A)=span⁡{n}N(A)=\operatorname{span}\{\boldsymbol{n}\},dim⁡=1=n−r\dim=1=n-r ✓。

一句话:满秩分解给出「活的」部分(R(A)R(A) 的基在 BB 里),而「死的」部分(N(A)N(A))要另解方程——但两者的维数由 rr 一分为二:2+1=32+1=3 ✓。

C13 · 右逆与全部右逆(2×32\times3)

设 A=[123011]A=\begin{bmatrix}1&amp;2&amp;3\\ 0&amp;1&amp;1\end{bmatrix}:用公式构造一个右逆,并写出全部右逆。

提示 · C13

先算 AAHAA^{\mathrm H} 与其逆;「全部右逆」的差集是「每列属于 N(A)N(A)」的矩阵——N(A)N(A) 的基是 (−1,−1,1)T(-1,-1,1)^{\mathrm T}。

解答 · C13

基本解:AAH=[14552]AA^{\mathrm H}=\begin{bmatrix}14&amp;5\\5&amp;2\end{bmatrix},det⁡=3\det=3,(AAH)−1=13[2−5−514](AA^{\mathrm H})^{-1}=\dfrac13\begin{bmatrix}2&amp;-5\\-5&amp;14\end{bmatrix},

B0=AH(AAH)−1=13[102131][2−5−514]=13[2−5−141−1],AB0=I2 ✓B_0=A^{\mathrm H}(AA^{\mathrm H})^{-1}=\frac13\begin{bmatrix}1&amp;0\\2&amp;1\\3&amp;1\end{bmatrix}\begin{bmatrix}2&amp;-5\\-5&amp;14\end{bmatrix}=\frac13\begin{bmatrix}2&amp;-5\\ -1&amp;4\\ 1&amp;-1\end{bmatrix},\qquad AB_0=I_2\ \checkmark

N(A)N(A):x1+2x2+3x3=0x_1+2x_2+3x_3=0、x2+x3=0x_2+x_3=0 得 n=(−1,−1,1)T\boldsymbol{n}=(-1,-1,1)^{\mathrm T}(验算 An=(−1−2+3,−1+1)T=0A\boldsymbol{n}=(-1-2+3,-1+1)^{\mathrm T}=\boldsymbol{0} ✓)。

全部右逆:

B=B0+[−s−t−s−tst],s,t∈CB=B_0+\begin{bmatrix}-s&amp;-t\\ -s&amp;-t\\ s&amp;t\end{bmatrix},\qquad s,t\in\mathbb C

验算:A[−s−t−s−tst]=[−s−2s+3s−t−2t+3t−s+s−t+t]=OA\begin{bmatrix}-s&amp;-t\\-s&amp;-t\\s&amp;t\end{bmatrix}=\begin{bmatrix}-s-2s+3s&amp;-t-2t+3t\\ -s+s&amp;-t+t\end{bmatrix}=O ✓。∎

C14 · AXA=AAXA=A 的一个解

设 A=[1111]A=\begin{bmatrix}1&amp;1\\ 1&amp;1\end{bmatrix}:用满秩分解求出 AXA=AAXA=A 的一个解 XX 并验证。

提示 · C14

r=1r=1:取 B=[11]B=\begin{bmatrix}1\\1\end{bmatrix}、C=[1  1]C=[1\ \ 1],代 X=CH(CCH)−1(BHB)−1BHX=C^{\mathrm H}(CC^{\mathrm H})^{-1}(B^{\mathrm H}B)^{-1}B^{\mathrm H} 再验算 AXA=AAXA=A。

解答 · C14

满秩分解:B=[11]B=\begin{bmatrix}1\\1\end{bmatrix}、C=[1  1]C=[1\ \ 1],BC=ABC=A ✓。

构造 XX:

X=CH(CCH)−1(BHB)−1BH=[11]⋅12⋅12⋅[1  1]=14[1111].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&amp;1\\ 1&amp;1\end{bmatrix}.

验算:

AXA=14[1111][1111][1111]=14[2222][1111]=14[4444]=[1111]=A ✓AXA=\frac14\begin{bmatrix}1&amp;1\\1&amp;1\end{bmatrix}\begin{bmatrix}1&amp;1\\1&amp;1\end{bmatrix}\begin{bmatrix}1&amp;1\\1&amp;1\end{bmatrix} =\frac14\begin{bmatrix}2&amp;2\\2&amp;2\end{bmatrix}\begin{bmatrix}1&amp;1\\1&amp;1\end{bmatrix}=\frac14\begin{bmatrix}4&amp;4\\4&amp;4\end{bmatrix}=\begin{bmatrix}1&amp;1\\1&amp;1\end{bmatrix}=A\ \checkmark

补充核对(用 本编补注):CXB=[1  1]⋅14[1111]⋅[11]=[1  1]⋅14[22]=1=I1CXB=[1\ \ 1]\cdot\frac14\begin{bmatrix}1&amp;1\\1&amp;1\end{bmatrix}\cdot\begin{bmatrix}1\\1\end{bmatrix}=[1\ \ 1]\cdot\frac14\begin{bmatrix}2\\2\end{bmatrix}=1=I_1 ✓。

观察:解不唯一——所有元素之和为 11 的 2×22\times2 矩阵都满足 AXA=AAXA=A(因为 AXA=(∑ijxij) AAXA=(\sum_{ij}x_{ij})\,A)。∎

C15 · 判断是否为满秩分解

指出下列各式是否为满秩分解(说明理由):

(a) [1234]=[12][1  2]\begin{bmatrix}1&amp;2\\ 3&amp;4\end{bmatrix}=\begin{bmatrix}1\\2\end{bmatrix}[1\ \ 2];(b) [1234]=I2[1234]\begin{bmatrix}1&amp;2\\ 3&amp;4\end{bmatrix}=I_2\begin{bmatrix}1&amp;2\\ 3&amp;4\end{bmatrix};(c) [2412]=[21][1  2]\begin{bmatrix}2&amp;4\\ 1&amp;2\end{bmatrix}=\begin{bmatrix}2\\1\end{bmatrix}[1\ \ 2];(d) [2412]=[11][2  4]\begin{bmatrix}2&amp;4\\ 1&amp;2\end{bmatrix}=\begin{bmatrix}1\\1\end{bmatrix}[2\ \ 4]。

解答 · C15

(a) 不是:右边乘积 =[1224]≠A=\begin{bmatrix}1&amp;2\\2&amp;4\end{bmatrix}\neq A(且 AA 可逆、秩 22,而右边秩 ≤1\le1——秩都不等,必不是)。 (b) 是:I2I_2 列满秩(秩 22)、AA 行满秩(2×22\times2 秩 22),而 A=I2AA=I_2A——注意这正是 注 2(注 2)型的「平凡但合法」的分解。 (c) 是:B=[21]B=\begin{bmatrix}2\\1\end{bmatrix} 列满秩(11 列,非零)、C=[1  2]C=[1\ \ 2] 行满秩(11 行,非零),乘积 =[2412]=A=\begin{bmatrix}2&amp;4\\1&amp;2\end{bmatrix}=A ✓;r=1r=1 ✓。 (d) 不是:乘积 =[2424]≠[2412]=\begin{bmatrix}2&amp;4\\2&amp;4\end{bmatrix}\neq\begin{bmatrix}2&amp;4\\1&amp;2\end{bmatrix}(BB 的列不在 R(A)R(A) 的方向上:R(A)=span⁡{(2,1)T}R(A)=\operatorname{span}\{(2,1)^{\mathrm T}\},而 (1,1)T(1,1)^{\mathrm T} 不属于它)。

要点:判断满秩分解要过三关——乘积等于 AA、BB 列满秩、CC 行满秩;且 BB 的列必须落在 R(A)R(A) 里((d) 就是这一关不过)。

C16 · 从相抵分解读出满秩分解

设 AA 如 C3(3×33\times3、秩 22)。用相抵分解的「方法 1」(定理 3.2.1)从 E−1E^{-1}、F−1F^{-1} 直接读出 AA 的一个满秩分解,并验算。

解答 · C16

方法 1 的配方:A=P−1ΓQ=P−1[Ir0][Ir, 0]QA=P^{-1}\Gamma Q=P^{-1}\begin{bmatrix}I_r\\0\end{bmatrix}[I_r,\ 0]Q,故

B=P−1[Ir0]=(P−1 的前 r 列),C=[Ir, 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{ 行)}.

代入 C3 的结果 P−1=E−1=[100210311]P^{-1}=E^{-1}=\begin{bmatrix}1&amp;0&amp;0\\2&amp;1&amp;0\\3&amp;1&amp;1\end{bmatrix}、Q=F−1=[123011001]Q=F^{-1}=\begin{bmatrix}1&amp;2&amp;3\\0&amp;1&amp;1\\0&amp;0&amp;1\end{bmatrix},取前 r=2r=2 列 / 行:

B=[102131],C=[123011].B=\begin{bmatrix}1&amp;0\\ 2&amp;1\\ 3&amp;1\end{bmatrix},\qquad C=\begin{bmatrix}1&amp;2&amp;3\\ 0&amp;1&amp;1\end{bmatrix}.

验算:

BC=[102131][123011]=[1232573710]=A ✓BC=\begin{bmatrix}1&amp;0\\2&amp;1\\3&amp;1\end{bmatrix}\begin{bmatrix}1&amp;2&amp;3\\0&amp;1&amp;1\end{bmatrix}=\begin{bmatrix}1&amp;2&amp;3\\2&amp;5&amp;7\\3&amp;7&amp;10\end{bmatrix}=A\ \checkmark

且 BB 列满秩(两列无关)、CC 行满秩(两行无关)✓。

读信息: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}\}——注意这不是 AA 的前两列(AA 的前两列是 (1,2,3)T,(2,5,7)T(1,2,3)^{\mathrm T},(2,5,7)^{\mathrm T}):同一组基可以由不同的向量担任,正印证了注 2(注 2)「基不一定从 AA 的列中选取」。∎

C17 · GS 手算 QR 分解(2×22\times2)

求 A=[2112]A=\begin{bmatrix}2&amp;1\\ 1&amp;2\end{bmatrix} 的 QR 分解(QQ 正交、RR 正线上三角),并回代验证。

提示 · C17

照 定理 3.4.1 存在性证明的流程:先正交化得 yk\boldsymbol{y}_k,对角元取 ∥yk∥\|\boldsymbol{y}_k\|(残差的范数,不是 ∥ak∥\|\boldsymbol{a}_k\|),非对角元取 (ak,zi)(\boldsymbol{a}_k,\boldsymbol{z}_i)。

解答 · C17

det⁡A=3≠0\det A=3\neq0,满秩 ✓。a1=(2,1)T\boldsymbol{a}_1=(2,1)^{\mathrm T}、a2=(1,2)T\boldsymbol{a}_2=(1,2)^{\mathrm T}。

GS:∥y1∥=∥a1∥=5\|\boldsymbol{y}_1\|=\|\boldsymbol{a}_1\|=\sqrt5,z1=15(2,1)T\boldsymbol{z}_1=\tfrac{1}{\sqrt5}(2,1)^{\mathrm T};(a2,z1)=45(\boldsymbol{a}_2,\boldsymbol{z}_1)=\tfrac{4}{\sqrt5},

y2=a2−45z1=(1,2)T−45(2,1)T=(−35, 65)T,∥y2∥=35,z2=15(−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}.

拼装:

Q=15[2−112],R=[545035].Q=\frac{1}{\sqrt5}\begin{bmatrix}2&amp;-1\\ 1&amp;2\end{bmatrix},\qquad R=\begin{bmatrix}\sqrt5&amp;\tfrac{4}{\sqrt5}\\ 0&amp;\tfrac{3}{\sqrt5}\end{bmatrix}.

验算:QR=[285−35145+65]=[2112]=AQR=\begin{bmatrix}2&amp;\tfrac85-\tfrac35\\ 1&amp;\tfrac45+\tfrac65\end{bmatrix}=\begin{bmatrix}2&amp;1\\ 1&amp;2\end{bmatrix}=A ✓;QTQ=IQ^{\mathrm T}Q=I ✓、det⁡Q=4+15=1\det Q=\tfrac{4+1}{5}=1(这又是一个 Givens 旋转,Q=[cos⁡θ−sin⁡θsin⁡θcos⁡θ]Q=\begin{bmatrix}\cos\theta&amp;-\sin\theta\\\sin\theta&amp;\cos\theta\end{bmatrix},cos⁡θ=25\cos\theta=\tfrac2{\sqrt5}——对称矩阵的 QR 因子常常「长得好看」)。∎

C18 · 长方阵的完全 QR 分解(推论 3.4.1 实操)

设 B=[100111]B=\begin{bmatrix}1&amp;0\\ 0&amp;1\\ 1&amp;1\end{bmatrix}。(a) 求 B=Q1R1B=Q_1R_1(Q1Q_1 列正交规范、R1R_1 正线上三角);(b) 补齐 C3\mathbb C^3 的标准正交基,写出 B=QRB=QR(QQ 酉、R=[R1O]R=\begin{bmatrix}R_1\\O\end{bmatrix}),并验证 QHQ=I3Q^{\mathrm H}Q=I_3。

解答 · C18

(a) a1=(1,0,1)T\boldsymbol{a}_1=(1,0,1)^{\mathrm T}:∥a1∥=2\|\boldsymbol{a}_1\|=\sqrt2,z1=12(1,0,1)T\boldsymbol{z}_1=\tfrac{1}{\sqrt2}(1,0,1)^{\mathrm T}。(a2,z1)=12(\boldsymbol{a}_2,\boldsymbol{z}_1)=\tfrac{1}{\sqrt2}(a2=(0,1,1)T\boldsymbol{a}_2=(0,1,1)^{\mathrm T}):

y2=a2−12(1,0,1)T=(−12, 1, 12)T,∥y2∥=62,z2=16(−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}.

Q1=[12−160261216],R1=[212062],Q1HQ1=I2 ✓Q_1=\begin{bmatrix}\tfrac{1}{\sqrt2}&amp;-\tfrac{1}{\sqrt6}\\ 0&amp;\tfrac{2}{\sqrt6}\\ \tfrac{1}{\sqrt2}&amp;\tfrac{1}{\sqrt6}\end{bmatrix},\qquad R_1=\begin{bmatrix}\sqrt2&amp;\tfrac{1}{\sqrt2}\\ 0&amp;\tfrac{\sqrt6}{2}\end{bmatrix},\qquad Q_1^{\mathrm H}Q_1=I_2\ \checkmark

验算 Q1R1Q_1R_1:第一列 2z1=(1,0,1)T\sqrt2\boldsymbol{z}_1=(1,0,1)^{\mathrm T} ✓;第二列 12z1+62z2=(12,0,12)T+(−12,1,12)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} ✓。

(b) 取 z3\boldsymbol{z}_3 同时垂直于 z1∥(1,0,1)\boldsymbol{z}_1\parallel(1,0,1) 与 z2∥(−1,2,1)\boldsymbol{z}_2\parallel(-1,2,1):叉积 (1,0,1)×(−1,2,1)=(−2,−2,2)(1,0,1)\times(-1,2,1)=(-2,-2,2),单位化得 z3=13(1,1,−1)T\boldsymbol{z}_3=\tfrac{1}{\sqrt3}(1,1,-1)^{\mathrm T}(两两内积为零 ✓)。于是

Q=[12−1613026131216−13] (酉),B=Q[21206200] ✓Q=\begin{bmatrix}\tfrac{1}{\sqrt2}&amp;-\tfrac{1}{\sqrt6}&amp;\tfrac{1}{\sqrt3}\\ 0&amp;\tfrac{2}{\sqrt6}&amp;\tfrac{1}{\sqrt3}\\ \tfrac{1}{\sqrt2}&amp;\tfrac{1}{\sqrt6}&amp;-\tfrac{1}{\sqrt3}\end{bmatrix}\ (\text{酉}),\qquad B=Q\begin{bmatrix}\sqrt2&amp;\tfrac{1}{\sqrt2}\\ 0&amp;\tfrac{\sqrt6}{2}\\ 0&amp;0\end{bmatrix}\ \checkmark

QHQ=I3Q^{\mathrm H}Q=I_3 逐列验算 ✓。∎

C19 · 高次幂的两条路

设 A=[4−121]A=\begin{bmatrix}4&amp;-1\\ 2&amp;1\end{bmatrix}。(a) 求 fA(λ)f_A(\lambda);(b) 用 H–C 定理把 A10A^{10} 写成 AA 的一次多项式并算出显式结果;(c) 该矩阵可对角化吗?用对角化路复核 (b) 的系数。

解答 · C19

(a) tr⁡A=5\operatorname{tr}A=5、det⁡A=6\det A=6,fA(λ)=λ2−5λ+6=(λ−2)(λ−3)f_A(\lambda)=\lambda^2-5\lambda+6=(\lambda-2)(\lambda-3)。

(b) H–C:A2=5A−6IA^2=5A-6I,故 A10=αA+βIA^{10}=\alpha A+\beta I。在特征值处定系数(h(λ)=λ10−αλ−βh(\lambda)=\lambda^{10}-\alpha\lambda-\beta 在 2,32,3 处为零 ⇒fA∣h⇒h(A)=O\Rightarrow f_A\mid h\Rightarrow h(A)=O):

210=2α+β,310=3α+β ⟹ α=310−210=58025,β=211−310=−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 .

A10=58025 A−115026 I=[4×58025−115026−580252×5802558025−115026]=[117074−58025116050−57001].A^{10}=58025\,A-115026\,I=\begin{bmatrix}4\times58025-115026&amp;-58025\\ 2\times58025&amp;58025-115026\end{bmatrix} =\begin{bmatrix}117074&amp;-58025\\ 116050&amp;-57001\end{bmatrix}.

(本编数值复核:与双精度 A10A^{10} 逐位一致 ✓。)

(c) fAf_A 无重根 ⇒\Rightarrow 可对角化(不同特征值的特征向量线性无关 2.4.5):λ=2\lambda=2 时 (A−2I)x=0(A-2I)\boldsymbol{x}=\boldsymbol{0} 给 p1=(1,2)T\boldsymbol{p}_1=(1,2)^{\mathrm T};λ=3\lambda=3 给 p2=(1,1)T\boldsymbol{p}_2=(1,1)^{\mathrm T}。P=[p1,p2]P=[\boldsymbol{p}_1,\boldsymbol{p}_2],P−1AP=diag⁡(2,3)P^{-1}AP=\operatorname{diag}(2,3),A10=Pdiag⁡(210,310)P−1A^{10}=P\operatorname{diag}(2^{10},3^{10})P^{-1}。系数复核:α,β\alpha,\beta 满足的方程组正是「λm=αλ+β\lambda^m=\alpha\lambda+\beta 在 λ=2,3\lambda=2,3 处成立」,两条路的定系数步骤相同——对角化路把 PP 显式算出来,H–C 路根本不需要 PP(且对不可对角化的矩阵照用,见 例 3.5.2)。∎

C20 · 带余除法与最小多项式算 A4A^4

设 A=[11001000−1]A=\begin{bmatrix}1&amp;1&amp;0\\ 0&amp;1&amp;0\\ 0&amp;0&amp;-1\end{bmatrix}。(a) 求 fAf_A 与 mAm_A;(b) 用「除以 mAm_A」计算 A4A^4;(c) 直接计算 A4A^4 核对。

解答 · C20

(a) 分块三角:fA(λ)=(λ−1)2(λ+1)f_A(\lambda)=(\lambda-1)^2(\lambda+1)。候选 (λ−1)(λ+1)(\lambda-1)(\lambda+1)、(λ−1)2(λ+1)(\lambda-1)^2(\lambda+1)。试低次:(A−I)(A+I)=[020000000]≠O(A-I)(A+I)=\begin{bmatrix}0&amp;2&amp;0\\0&amp;0&amp;0\\0&amp;0&amp;0\end{bmatrix}\neq O,落选;故 mA=(λ−1)2(λ+1)=fAm_A=(\lambda-1)^2(\lambda+1)=f_A。

(b) λ4=mA⋅(λ+1)+(2λ2−1)\lambda^4=m_A\cdot(\lambda+1)+\big(2\lambda^2-1\big)(验算:(λ−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 ✓)。故

A4=mA(A)(A+I)+2A2−I=2A2−I.A^4=m_A(A)(A+I)+2A^2-I=2A^2-I .

(c) A2=[120010001]A^2=\begin{bmatrix}1&amp;2&amp;0\\0&amp;1&amp;0\\0&amp;0&amp;1\end{bmatrix}((−1)2=1(-1)^2=1),故 A4=2A2−I=[140010001]A^4=2A^2-I=\begin{bmatrix}1&amp;4&amp;0\\0&amp;1&amp;0\\0&amp;0&amp;1\end{bmatrix}。直接核对:A4=(A2)2A^4=(A^2)^2,上三角块 [1201]2=[1401]\begin{bmatrix}1&amp;2\\0&amp;1\end{bmatrix}^2=\begin{bmatrix}1&amp;4\\0&amp;1\end{bmatrix} ✓,右下 (−1)4=1(−1)^4=1 ✓。∎

C21 · 候选法求最小多项式(低次落选型)

求 A=[210021004]A=\begin{bmatrix}2&amp;1&amp;0\\ 0&amp;2&amp;1\\ 0&amp;0&amp;4\end{bmatrix} 的最小多项式。

解答 · C21

fA(λ)=(λ−2)2(λ−4)f_A(\lambda)=(\lambda-2)^2(\lambda-4)(上三角,对角 2,2,42,2,4)。候选:(λ−2)(λ−4)(\lambda-2)(\lambda-4) 与 (λ−2)2(λ−4)(\lambda-2)^2(\lambda-4)。

试低次:A−2I=[010001002]A-2I=\begin{bmatrix}0&amp;1&amp;0\\0&amp;0&amp;1\\0&amp;0&amp;2\end{bmatrix},A−4I=[−2100−21000]A-4I=\begin{bmatrix}-2&amp;1&amp;0\\0&amp;-2&amp;1\\0&amp;0&amp;0\end{bmatrix}。乘积第一行 =A−2I=A-2I 第一行 (0,1,0)(0,1,0) 乘开 =A−4I=A-4I 的第二行 =(0,−2,1)≠0=(0,-2,1)\neq\boldsymbol{0}——落选。

试高次:(A−2I)2=[001002004](A-2I)^2=\begin{bmatrix}0&amp;0&amp;1\\0&amp;0&amp;2\\0&amp;0&amp;4\end{bmatrix},(A−2I)2(A−4I)(A-2I)^2(A-4I) 的每一行都是 (0,0,∗)(0,0,*) 乘 A−4IA-4I(第三行为零行),故 =O=O ✓。

mA(λ)=(λ−2)2(λ−4).m_A(\lambda)=(\lambda-2)^2(\lambda-4).

读法:λ=2\lambda=2 的特征子空间 N(A−2I)N(A-2I) 由 e1\boldsymbol{e}_1 张成(解方程得 x2=x3=0x_2=x_3=0),几何重数 1<1&lt; 代数重数 22(几何重数的计算 2.4.6)——最小多项式里 (λ−2)(\lambda-2) 必须带平方,与 例 3.5.B 同型。∎

C22 · QR 迭代手算三轮(分数精确)

对 A=[2112]A=\begin{bmatrix}2&amp;1\\ 1&amp;2\end{bmatrix}(即 C17)执行 QR 算法三轮:写出 A2A_2、A3A_3(精确分数)与 A4A_4(四位小数),观察对角元向何值收敛、非对角元衰减多快。

解答 · C22

第一轮:A1=Q1R1A_1=Q_1R_1 用 C17 的结果,A2=R1Q1A_2=R_1Q_1:

A2=[545035]⋅15[2−112]=[145353565]=[2.80.60.61.2].A_2=\begin{bmatrix}\sqrt5&amp;\tfrac{4}{\sqrt5}\\0&amp;\tfrac{3}{\sqrt5}\end{bmatrix}\cdot\frac{1}{\sqrt5}\begin{bmatrix}2&amp;-1\\1&amp;2\end{bmatrix} =\begin{bmatrix}\tfrac{14}{5}&amp;\tfrac35\\ \tfrac35&amp;\tfrac65\end{bmatrix} =\begin{bmatrix}2.8&amp;0.6\\ 0.6&amp;1.2\end{bmatrix}.

第二轮:对 A2A_2 再 GS(z1=1205(14,3)T\boldsymbol{z}_1=\tfrac1{\sqrt{205}}(14,3)^{\mathrm T},(a2,z1)=12205(\boldsymbol{a}_2,\boldsymbol{z}_1)=\tfrac{12}{\sqrt{205}},y2∥(−3,14)T\boldsymbol{y}_2\parallel(-3,14)^{\mathrm T}、∥y2∥=320541\|\boldsymbol{y}_2\|=\tfrac{3\sqrt{205}}{41}):

Q2=1205[14−3314],R2=[2055122050320541],A3=R2Q2=[122419419414241]≈[2.97560.21950.21951.0244].Q_2=\frac{1}{\sqrt{205}}\begin{bmatrix}14&amp;-3\\ 3&amp;14\end{bmatrix},\qquad R_2=\begin{bmatrix}\tfrac{\sqrt{205}}{5}&amp;\tfrac{12}{\sqrt{205}}\\ 0&amp;\tfrac{3\sqrt{205}}{41}\end{bmatrix},\qquad A_3=R_2Q_2=\begin{bmatrix}\tfrac{122}{41}&amp;\tfrac9{41}\\ \tfrac9{41}&amp;\tfrac{42}{41}\end{bmatrix} \approx\begin{bmatrix}2.9756&amp;0.2195\\ 0.2195&amp;1.0244\end{bmatrix}.

第三轮(数值):A4≈[2.9973−0.0740−0.07401.0027]A_4\approx\begin{bmatrix}2.9973&amp;-0.0740\\ -0.0740&amp;1.0027\end{bmatrix}。

观察:特征值精确解 3,13,1(fA=(λ−3)(λ−1)f_A=(\lambda-3)(\lambda-1),模不同 ⇒\Rightarrow 收敛条件满足)。对角元 2.8→2.976→2.997→32.8\to2.976\to2.997\to3、1.2→1.024→1.003→11.2\to1.024\to1.003\to1(按模从大到小排定 ✓);非对角元 0.6→0.22→0.0740.6\to0.22\to0.074——每轮约缩小到三分之一(收敛率由相邻特征值模比 ∣λ2∣∣λ1∣=13\tfrac{|\lambda_2|}{|\lambda_1|}=\tfrac13 决定,本编注)。∎

C23 · 用推论 3.5.2 求逆

设 A=[120010002]A=\begin{bmatrix}1&amp;2&amp;0\\ 0&amp;1&amp;0\\ 0&amp;0&amp;2\end{bmatrix}。(a) 求 fA(λ)f_A(\lambda);(b) 用逆矩阵公式求 A−1A^{-1};(c) 回代验证。

解答 · C23

(a) 上三角:fA(λ)=(λ−1)2(λ−2)=λ3−4λ2+5λ−2f_A(\lambda)=(\lambda-1)^2(\lambda-2)=\lambda^3-4\lambda^2+5\lambda-2。

(b) n=3n=3,a2=−4a_2=-4、a1=5a_1=5、a0=−2≠0a_0=-2\neq0(det⁡A=2\det A=2,(−1)3det⁡A=−2(-1)^3\det A=-2 ✓):

A−1=−1a0(A2+a2A+a1I)=12(A2−4A+5I).A^{-1}=-\frac{1}{a_0}\big(A^2+a_2A+a_1I\big)=\frac12\big(A^2-4A+5I\big).

A2=[140010004]A^2=\begin{bmatrix}1&amp;4&amp;0\\0&amp;1&amp;0\\0&amp;0&amp;4\end{bmatrix},故

A2−4A+5I=[2−40020001],A−1=[1−200100012].A^2-4A+5I=\begin{bmatrix}2&amp;-4&amp;0\\0&amp;2&amp;0\\0&amp;0&amp;1\end{bmatrix},\qquad A^{-1}=\begin{bmatrix}1&amp;-2&amp;0\\ 0&amp;1&amp;0\\ 0&amp;0&amp;\tfrac12\end{bmatrix}.

(c) AA−1=[120010002][1−200100012]=[1−2+20010001]=IAA^{-1}=\begin{bmatrix}1&amp;2&amp;0\\0&amp;1&amp;0\\0&amp;0&amp;2\end{bmatrix}\begin{bmatrix}1&amp;-2&amp;0\\0&amp;1&amp;0\\0&amp;0&amp;\frac12\end{bmatrix}=\begin{bmatrix}1&amp;-2+2&amp;0\\0&amp;1&amp;0\\0&amp;0&amp;1\end{bmatrix}=I ✓。∎

C24 · 最小多项式与矩阵多项式综合

设 A=[410040002]A=\begin{bmatrix}4&amp;1&amp;0\\ 0&amp;4&amp;0\\ 0&amp;0&amp;2\end{bmatrix}。(a) 求 fAf_A 与 mAm_A;(b) 用 mAm_A 计算 A3A^3;(c) 直接乘 A3A^3 核对。

解答 · C24

(a) fA(λ)=(λ−4)2(λ−2)f_A(\lambda)=(\lambda-4)^2(\lambda-2)。试低次 (λ−4)(λ−2)(\lambda-4)(\lambda-2):(A−4I)(A−2I)(A-4I)(A-2I) 的第一行 =A−4I=A-4I 第一行 (0,1,0)(0,1,0) 乘开 =A−2I=A-2I 第二行 =(0,2,0)≠0=(0,2,0)\neq\boldsymbol{0},落选(λ=4\lambda=4 的几何重数 1<21&lt;2,同 例 3.5.B 型)。故 mA=(λ−4)2(λ−2)=fAm_A=(\lambda-4)^2(\lambda-2)=f_A。

(b) λ3=mA⋅1+(10λ2−32λ+32)\lambda^3=m_A\cdot1+\big(10\lambda^2-32\lambda+32\big)(验算:(λ−4)2(λ−2)=λ3−10λ2+32λ−32(\lambda-4)^2(\lambda-2)=\lambda^3-10\lambda^2+32\lambda-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 ✓)。故

A3=10A2−32A+32I.A^3=10A^2-32A+32I .

(c) A2=[16800160004]A^2=\begin{bmatrix}16&amp;8&amp;0\\0&amp;16&amp;0\\0&amp;0&amp;4\end{bmatrix},10A2−32A+32I=[160−128+3280−3200160−128+3200040−64+32]=[644800640008]10A^2-32A+32I=\begin{bmatrix}160-128+32&amp;80-32&amp;0\\0&amp;160-128+32&amp;0\\0&amp;0&amp;40-64+32\end{bmatrix}=\begin{bmatrix}64&amp;48&amp;0\\0&amp;64&amp;0\\0&amp;0&amp;8\end{bmatrix}。直接核对:[4104]3=[6448064]\begin{bmatrix}4&amp;1\\0&amp;4\end{bmatrix}^3=\begin{bmatrix}64&amp;48\\0&amp;64\end{bmatrix}(超对角元 3×42=483\times4^2=48 ✓)、23=82^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 及以后素材到位后本页续写。