前置知识: 线性代数

解的结构

10 minIntermediate2026/6/14

线性方程组解的结构理论,通解、特解与基础解系的关系,解空间的维数与结构,解集的几何描述。

1. 解的结构总览

1.1 齐次方程组的解结构

Ax=0Ax = 0 的解集 S0S_0 构成向量空间(零空间),其结构为:

S0={k1ξ1+k2ξ2++ktξtkiR}S_0 = \{k_1\boldsymbol{\xi}_1 + k_2\boldsymbol{\xi}_2 + \cdots + k_t\boldsymbol{\xi}_t \mid k_i \in \mathbb{R}\}

其中 ξ1,,ξt\boldsymbol{\xi}_1, \ldots, \boldsymbol{\xi}_t 是基础解系,t=nr(A)t = n - r(A)

1.2 非齐次方程组的解结构

Ax=bAx = b 的解集 SS 是仿射子空间:

S=η+S0={η+ξξS0}S = \boldsymbol{\eta}^* + S_0 = \{\boldsymbol{\eta}^* + \boldsymbol{\xi} \mid \boldsymbol{\xi} \in S_0\}

其中 η\boldsymbol{\eta}^*Ax=bAx = b 的特解,S0S_0 是导出组 Ax=0Ax = 0 的解空间。

2. 通解、特解与基础解系的关系

2.1 核心关系

Ax = b 的通解 = 特解 + 导出组的通解
     x    =  η*  +  k₁ξ₁ + k₂ξ₂ + ... + kₜξₜ

2.2 各部分的作用

概念作用唯一性
特解 η\boldsymbol{\eta}^*确定解集的”位置”不唯一,任意特解均可
基础解系 ξ1,,ξt\boldsymbol{\xi}_1, \ldots, \boldsymbol{\xi}_t确定解集的”形状”不唯一,但等价
自由参数 k1,,ktk_1, \ldots, k_t参数化解集中的每个解依赖于基础解系的选取

2.3 特解的选取对通解的影响

η1\boldsymbol{\eta}_1^*η2\boldsymbol{\eta}_2^*Ax=bAx = b 的两个不同特解,则:

η2=η1+ξ0\boldsymbol{\eta}_2^* = \boldsymbol{\eta}_1^* + \boldsymbol{\xi}_0

其中 ξ0S0\boldsymbol{\xi}_0 \in S_0。两个通解表达式等价:

η2+kiξi=η1+ξ0+kiξi\boldsymbol{\eta}_2^* + \sum k_i\boldsymbol{\xi}_i = \boldsymbol{\eta}_1^* + \boldsymbol{\xi}_0 + \sum k_i\boldsymbol{\xi}_i

由于 ξ0\boldsymbol{\xi}_0 可由基础解系表示,故两个通解描述的是同一个解集。

3. 解空间的维数

3.1 维数公式

dim(S0)=nr(A)\dim(S_0) = n - r(A)

3.2 秩-零度定理

r(A)+dim(N(A))=nr(A) + \dim(N(A)) = n

这是线性代数中最基本的维数关系之一。

3.3 推广

对于 m×nm \times n 矩阵 AA

dim(Row(A))+dim(N(A))=n\dim(\text{Row}(A)) + \dim(N(A)) = n

dim(Col(A))+dim(N(AT))=m\dim(\text{Col}(A)) + \dim(N(A^T)) = m

4. 解集的几何描述

4.1 齐次方程组

Ax=0Ax = 0 的解集是 Rn\mathbb{R}^n 中过原点的 tt 维子空间(t=nr(A)t = n - r(A))。

  • t=0t = 0:解集为单点 {0}\{0\}
  • t=1t = 1:解集为过原点的直线
  • t=2t = 2:解集为过原点的平面
  • t=nt = n:解集为整个 Rn\mathbb{R}^nA=OA = O 时)

4.2 非齐次方程组

Ax=bAx = b 的解集是 Rn\mathbb{R}^n 中的 tt 维仿射子空间(平移后的子空间)。

  • t=0t = 0:解集为单点(唯一解)
  • t=1t = 1:解集为不过原点的直线
  • t=2t = 2:解集为不过原点的平面

4.3 示例

{x+y+z=1\begin{cases} x + y + z = 1 \end{cases}

r(A)=1r(A) = 1n=3n = 3dim(S0)=2\dim(S_0) = 2

解集是 R3\mathbb{R}^3 中不过原点的平面 x+y+z=1x + y + z = 1

特解:(1,0,0)T(1, 0, 0)^T,基础解系:(1,1,0)T(-1, 1, 0)^T(1,0,1)T(-1, 0, 1)^T

5. 解的结构与矩阵分解

5.1 与 SVD 的关系

A=UΣVTA = U\Sigma V^TAA 的奇异值分解,则:

  • N(A)=span{vr+1,,vn}N(A) = \text{span}\{v_{r+1}, \ldots, v_n\}VV 的后 nrn-r 列)
  • Col(A)=span{u1,,ur}\text{Col}(A) = \text{span}\{u_1, \ldots, u_r\}UU 的前 rr 列)

5.2 与特征值的关系

AA 可对角化为 A=PΛP1A = P\Lambda P^{-1},则:

  • N(A)N(A) 由对应零特征值的特征向量张成
  • Ax=bAx = b 的解可通过对角化后求解

6. 解的稳定性

6.1 条件数

矩阵 AA条件数定义为:

κ(A)=AA1\kappa(A) = \|A\| \cdot \|A^{-1}\|

条件数衡量了 Ax=bAx = b 的解对 bb 的扰动的敏感程度。

  • κ(A)1\kappa(A) \approx 1:良态问题,解稳定
  • κ(A)1\kappa(A) \gg 1:病态问题,解不稳定

6.2 扰动分析

bb 有扰动 δb\delta b,则解的扰动 δx\delta x 满足:

δxxκ(A)δbb\frac{\|\delta x\|}{\|x\|} \leq \kappa(A) \cdot \frac{\|\delta b\|}{\|b\|}

7. 典型例题

例1

AA3×43 \times 4 矩阵,η1=(1,0,1,0)T\boldsymbol{\eta}_1 = (1, 0, 1, 0)^Tη2=(0,1,0,1)T\boldsymbol{\eta}_2 = (0, 1, 0, 1)^TAx=0Ax = 0 的基础解系,η=(1,1,1,1)T\boldsymbol{\eta}^* = (1, 1, 1, 1)^TAx=bAx = b 的特解,求 Ax=bAx = b 的通解。

:通解 x=η+k1η1+k2η2x = \boldsymbol{\eta}^* + k_1\boldsymbol{\eta}_1 + k_2\boldsymbol{\eta}_2

x=(1111)+k1(1010)+k2(0101)=(1+k11+k21+k11+k2)x = \begin{pmatrix} 1 \\ 1 \\ 1 \\ 1 \end{pmatrix} + k_1\begin{pmatrix} 1 \\ 0 \\ 1 \\ 0 \end{pmatrix} + k_2\begin{pmatrix} 0 \\ 1 \\ 0 \\ 1 \end{pmatrix} = \begin{pmatrix} 1 + k_1 \\ 1 + k_2 \\ 1 + k_1 \\ 1 + k_2 \end{pmatrix}

例2

AAm×nm \times n 矩阵,r(A)=rr(A) = r,证明 Ax=bAx = b 的任意 nr+1n - r + 1 个解线性相关。

证明:设 η0\boldsymbol{\eta}_0 为特解,ηi=η0+ξi\boldsymbol{\eta}_i = \boldsymbol{\eta}_0 + \boldsymbol{\xi}_ii=1,,nri = 1, \ldots, n-r),其中 ξi\boldsymbol{\xi}_i 是基础解系。

考虑 nr+1n - r + 1 个解 η0,η1,,ηnr\boldsymbol{\eta}_0, \boldsymbol{\eta}_1, \ldots, \boldsymbol{\eta}_{n-r}

k0η0+k1η1++knrηnr=0k_0\boldsymbol{\eta}_0 + k_1\boldsymbol{\eta}_1 + \cdots + k_{n-r}\boldsymbol{\eta}_{n-r} = 0

(k0+k1++knr)η0+k1ξ1++knrξnr=0(k_0 + k_1 + \cdots + k_{n-r})\boldsymbol{\eta}_0 + k_1\boldsymbol{\xi}_1 + \cdots + k_{n-r}\boldsymbol{\xi}_{n-r} = 0

k0=1,k1=k2==knr=1nrk_0 = 1, k_1 = k_2 = \cdots = k_{n-r} = -\dfrac{1}{n-r},则 k0+ki=0k_0 + \sum k_i = 0,且 kiξi\sum k_i\boldsymbol{\xi}_iAx=0Ax = 0 的解。

但需要更仔细的分析。实际上,nr+1n - r + 1 个解向量 η0,η1,,ηnr\boldsymbol{\eta}_0, \boldsymbol{\eta}_1, \ldots, \boldsymbol{\eta}_{n-r}nn 维空间中,它们位于一个 nrn - r 维仿射子空间上,故必线性相关。