Ax=0 的解集 S0 构成向量空间(零空间),其结构为:
S0={k1ξ1+k2ξ2+⋯+ktξt∣ki∈R}
其中 ξ1,…,ξt 是基础解系,t=n−r(A)。
Ax=b 的解集 S 是仿射子空间:
S=η∗+S0={η∗+ξ∣ξ∈S0}
其中 η∗ 是 Ax=b 的特解,S0 是导出组 Ax=0 的解空间。
Ax = b 的通解 = 特解 + 导出组的通解
x = η* + k₁ξ₁ + k₂ξ₂ + ... + kₜξₜ
| 概念 | 作用 | 唯一性 |
|---|
| 特解 η∗ | 确定解集的”位置” | 不唯一,任意特解均可 |
| 基础解系 ξ1,…,ξt | 确定解集的”形状” | 不唯一,但等价 |
| 自由参数 k1,…,kt | 参数化解集中的每个解 | 依赖于基础解系的选取 |
设 η1∗ 和 η2∗ 是 Ax=b 的两个不同特解,则:
η2∗=η1∗+ξ0
其中 ξ0∈S0。两个通解表达式等价:
η2∗+∑kiξi=η1∗+ξ0+∑kiξi
由于 ξ0 可由基础解系表示,故两个通解描述的是同一个解集。
dim(S0)=n−r(A)
r(A)+dim(N(A))=n
这是线性代数中最基本的维数关系之一。
对于 m×n 矩阵 A:
dim(Row(A))+dim(N(A))=n
dim(Col(A))+dim(N(AT))=m
Ax=0 的解集是 Rn 中过原点的 t 维子空间(t=n−r(A))。
- t=0:解集为单点 {0}
- t=1:解集为过原点的直线
- t=2:解集为过原点的平面
- t=n:解集为整个 Rn(A=O 时)
Ax=b 的解集是 Rn 中的 t 维仿射子空间(平移后的子空间)。
- t=0:解集为单点(唯一解)
- t=1:解集为不过原点的直线
- t=2:解集为不过原点的平面
{x+y+z=1
r(A)=1,n=3,dim(S0)=2。
解集是 R3 中不过原点的平面 x+y+z=1。
特解:(1,0,0)T,基础解系:(−1,1,0)T,(−1,0,1)T。
设 A=UΣVT 为 A 的奇异值分解,则:
- N(A)=span{vr+1,…,vn}(V 的后 n−r 列)
- Col(A)=span{u1,…,ur}(U 的前 r 列)
若 A 可对角化为 A=PΛP−1,则:
- N(A) 由对应零特征值的特征向量张成
- Ax=b 的解可通过对角化后求解
矩阵 A 的条件数定义为:
κ(A)=∥A∥⋅∥A−1∥
条件数衡量了 Ax=b 的解对 b 的扰动的敏感程度。
- κ(A)≈1:良态问题,解稳定
- κ(A)≫1:病态问题,解不稳定
若 b 有扰动 δb,则解的扰动 δx 满足:
∥x∥∥δx∥≤κ(A)⋅∥b∥∥δb∥
设 A 为 3×4 矩阵,η1=(1,0,1,0)T,η2=(0,1,0,1)T 是 Ax=0 的基础解系,η∗=(1,1,1,1)T 是 Ax=b 的特解,求 Ax=b 的通解。
解:通解 x=η∗+k1η1+k2η2
x=1111+k11010+k20101=1+k11+k21+k11+k2
设 A 为 m×n 矩阵,r(A)=r,证明 Ax=b 的任意 n−r+1 个解线性相关。
证明:设 η0 为特解,ηi=η0+ξi(i=1,…,n−r),其中 ξi 是基础解系。
考虑 n−r+1 个解 η0,η1,…,ηn−r:
k0η0+k1η1+⋯+kn−rηn−r=0
(k0+k1+⋯+kn−r)η0+k1ξ1+⋯+kn−rξn−r=0
取 k0=1,k1=k2=⋯=kn−r=−n−r1,则 k0+∑ki=0,且 ∑kiξi 是 Ax=0 的解。
但需要更仔细的分析。实际上,n−r+1 个解向量 η0,η1,…,ηn−r 在 n 维空间中,它们位于一个 n−r 维仿射子空间上,故必线性相关。