线性方程组

相关概念

则线性方程组nn个未知数mm个方程的线性方程组的一般形式

{a11x1+a12x2++a1nxn=b1a21x1+a22x2++a2nxn=b2                                                am1x1+am2x2++amnxn=bn\begin{cases} a_{11}x_{1}+a_{12}x_{2}+\dots+a_{1n}x_{n}=b_{1} \\ a_{21}x_{1}+a_{22}x_{2}+\dots+a_{2n}x_{n}=b_{2} \\ \;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\vdots \\ a_{m1}x_{1}+a_{m_{2}}x_{2}+\dots+a_{mn}x_{n}=b_{n} \end{cases}

若常数项b1,b2,,bnb_{1},b_{2},\dots,b_{n}不全为零则称方程组为非齐次线性方程组
若常数项b1,b2,,bnb_{1},b_{2},\dots,b_{n}全为零则称方程组为齐次线性方程组
矩阵形式Ax=b\boldsymbol{Ax=b}
向量形式x1α1+x2α2++xnαn=bx_{1}\boldsymbol{\alpha_{1}}+x_{2}\boldsymbol{\alpha_{2}}+\cdots+x_{n}\boldsymbol{\alpha_{n}}=\boldsymbol{b}
其中A=(a11a12a1na21a22a2nam1am2amn)=(α1α2αn):系数矩阵\boldsymbol{A=}\begin{pmatrix}a_{11} & a_{12} & \dots & a_{1n} \\a_{21} & a_{22} &\dots & a_{2n} \\\vdots & \vdots & \ddots & \vdots \\a_{m1} & a_{m2} & \dots & a_{mn}\end{pmatrix}=\begin{pmatrix}\boldsymbol{\alpha_{1}} & \boldsymbol{\alpha_{2}} & \dots & \boldsymbol{\alpha_{n}}\end{pmatrix}:系数矩阵
x=(x1x2xn)b=(b1b2bn)y=(y1y2yn)\boldsymbol{x=}\begin{pmatrix}x_{1} \\x_{2}\\\vdots \\x_{n}\end{pmatrix},\boldsymbol{b=}\begin{pmatrix}b_{1} \\b_{2} \\\vdots \\b_{n}\end{pmatrix},\boldsymbol{y= }\begin{pmatrix}y_{1} \\y_{2} \\\vdots \\y_{n}\end{pmatrix}
(A,b)=(a11a12a1nb1a21a22a2nb2am1am2amnbm)=(α1,α2,,αn,b):增广矩阵(\boldsymbol{A,b})=\begin{pmatrix}a_{11} & a_{12} & \dots & a_{1n} & b_{1}\\a_{21} & a_{22} &\dots & a_{2n} & b_{2}\\\vdots & \vdots & \ddots & \vdots & \vdots \\a_{m1} & a_{m2} & \dots & a_{mn} & b_{m}\end{pmatrix}=\begin{pmatrix}\boldsymbol{\alpha_{1}},\boldsymbol{\alpha_{2}},\dots,\boldsymbol{\alpha_{n}} ,\boldsymbol{b}\end{pmatrix}:增广矩阵
如果 x1=c1,x2=c2,,xn=cnx_1 = c_1, x_2 = c_2, \cdots, x_n = c_n 满足 Ax=b\boldsymbol{Ax = b} 中的每一个方程,则称 x=(c1,c2,,cn)T\boldsymbol{x}= (c_1, c_2, \cdots, c_n)^T 为方程组的一个解(向量)
显然 x1=0,x2=0,,xn=0x_1 = 0, x_2 = 0, \cdots, x_n = 0Ax=0\boldsymbol{Ax = 0} 的一个解,称之为零解平凡解
解集:线性方程组的所有解构成的集合。
相容方程组:有解的线性方程组。
不相容(矛盾)方程组:无解的线性方程组。
同解方程组:有相同的解集的线性方程组。

高斯消元法

将方程组看成一个整体加以等价变形,对方程组做三种变换:方程组变形为“阶梯形方程组”,再回代求解。

{(1) 交换某两个方程的位置;(2) 某个方程除以非零常数;(3) 将某个方程的 k 倍加到另一个方程上。\left\{ \begin{aligned} &(1)\ \text{交换某两个方程的位置;} \\ &(2)\ \text{某个方程除以非零常数;} \\ &(3)\ \text{将某个方程的 } k \text{ 倍加到另一个方程上。} \end{aligned} \right.

这三种变换对应矩阵的三种初等行变换
上述过程可以在方程组的增广矩阵上实现:对方程组的增广矩阵进行初等行变换,使之化为行最简形,从而读出方程组的解

是否有解的判别方法

增广矩阵角度

设线性方程组Ax=b\boldsymbol{Ax = b}的增广矩阵的行最简形矩阵为

(A,b)(100b1,r+1b1nd1010b2,r+1b2nd2001br,r+1brndr00000dr+1000000000000)(\boldsymbol{A,b}) \sim \begin{pmatrix} 1 & 0 & \cdots & 0 & b_{1,r+1} & \cdots & b_{1n} & d_1 \\ 0 & 1 & \cdots & 0 & b_{2,r+1} & \cdots & b_{2n} & d_2 \\ \vdots & \vdots & \ddots & \vdots & \vdots & & \vdots & \vdots \\ 0 & 0 & \cdots & 1 & b_{r,r+1} & \cdots & b_{rn} & d_r \\ 0 & 0 & \cdots & 0 & 0 & \cdots & 0 & d_{r+1} \\ 0 & 0 & \cdots & 0 & 0 & \cdots & 0 & 0 \\ \vdots & \vdots & & \vdots & \vdots & & \vdots & \vdots \\ 0 & 0 & \cdots & 0 & 0 & \cdots & 0 & 0 \\ \end{pmatrix}

则有:

  1. dr+10d_{r+1} \neq 0 时,方程组无解
  2. dr+1=0d_{r+1} = 0r=nr = n 时,方程组有唯一解
  3. dr+1=0d_{r+1} = 0r<nr < n 时,方程组有无穷多解
  4. Ax=0\boldsymbol{Ax = 0} 总是有解,且当 r=nr = n 时,方程组只有零解,r<nr < n 时,有非零解

系数矩阵角度

A\boldsymbol{A}m×nm\times n矩阵,对齐次线性方程组Ax=0\boldsymbol{Ax=0}来说,下列命题等价:

  1. Ax=0\boldsymbol{Ax=0}有非零解
  2. A\boldsymbol{A}的列向量组线性相关
  3. r(A)<nr(\boldsymbol{A})<n
  4. 若有m=n,A=0m=n,则|\boldsymbol{A}|=0

A\boldsymbol{A}m×nm\times n矩阵,对非齐次线性方程组Ax=b\boldsymbol{Ax=b}来说,下列命题等价:

  1. Ax=b\boldsymbol{Ax=b}有解
  2. b\boldsymbol{b}可用A\boldsymbol{A}的列向量组线性表示
  3. r(A)=r(A,b)r(\boldsymbol{A})=r(\boldsymbol{A,b})
    r(A)<r(A,b)Ax=br(\boldsymbol{A})<r(\boldsymbol{A,b})\Leftrightarrow\boldsymbol{Ax=b}无解,且此时r(A,b)=r(A)+1r(\boldsymbol{A,b})=r(\boldsymbol{A})+1

齐次线性方程组的解

性质

ξ1,ξ2,,ξs\boldsymbol{\xi_{1},\xi_{2},\dots,\xi_{s}}Ax=0\boldsymbol{Ax=0}的解,则k1,k2,,ksR\forall k_{1},k_{2},\dots,k_{s}\in \mathbb{R}k1ξ1+k2ξ2++ksξsk_{1}\boldsymbol{\xi_{1}}+k_{2}\boldsymbol{\xi_{2}}+\dots+k_{s}\boldsymbol{\xi_{s}}也是Ax=0\boldsymbol{Ax=0}的解

基础解系

定义

α1,α2,,αs\boldsymbol{\alpha_{1}},\boldsymbol{\alpha_{2}},\dots,\boldsymbol{\alpha_{s}}是齐次线性方程组Ax=0\boldsymbol{Ax=0}ss个解,如果满足以下条件:

  1. α1,α2,,αs\boldsymbol{\alpha_{1}},\boldsymbol{\alpha_{2}},\dots,\boldsymbol{\alpha_{s}}线性无关
  2. Ax=0\boldsymbol{Ax=0}的任一解α\boldsymbol{\alpha}都可由α1,α2,,αs\boldsymbol{\alpha_{1}},\boldsymbol{\alpha_{2}},\dots,\boldsymbol{\alpha_{s}}线性表示
    则称α1,α2,,αs\boldsymbol{\alpha_{1}},\boldsymbol{\alpha_{2}},\dots,\boldsymbol{\alpha_{s}}Ax=0\boldsymbol{Ax=0}基础解系

Ax=0\boldsymbol{Ax=0}的基础解系就是它的解向量组的一个极大无关组
Ax=0\boldsymbol{Ax=0}的基础解系不唯一,但基础解系中向量的个数是唯一的

性质

nn元齐次线性方程组Am×nx=0\boldsymbol{A}_{m\times n}\boldsymbol{x=0}r(A)=r<nr(\boldsymbol{A})=r<n时存在基础解系,且基础解系中含有nrn-r个解向量

  1. r(A)=nr(\boldsymbol{A})=n时,方程组Am×nx=0\boldsymbol{A}_{m\times n}\boldsymbol{x=0}只有零解,此时没有基础解系
  2. r(A)=r<nr(\boldsymbol{A})=r<n时,方程组Am×nx=0\boldsymbol{A}_{m\times n}\boldsymbol{x=0}的基础解系为ξ1,ξ2,,ξnr\xi_{1},\xi_{2},\dots,\xi_{n-r}Am×nx=0\boldsymbol{A}_{m\times n}\boldsymbol{x=0}的解集SS可表示为S={x=k1ξ1+k2ξ2++knrξnr}S=\left\{ \boldsymbol{x}=k_{1}\boldsymbol{\xi_{1}}+k_{2}\boldsymbol{\xi_{2}}+\dots+k_{n-r}\boldsymbol{\xi_{n-r}} \right\}
    SSAm×nx=0\boldsymbol{A}_{m\times n}\boldsymbol{x=0}解空间(在线性变换Am×n\boldsymbol{A}_{m\times n}中被压缩至原点的空间,也叫零空间),它有基础解系张成,dimS=nr\dim S=n-r
  3. "nrn-r"的含义:
  4. nn:未知数的个数;系数矩阵A\boldsymbol{A}的列数
  5. rr:系数矩阵A\boldsymbol{A}的秩;同解方程组中有效方程组的个数;非自由未知数的个数
  6. nrn-r:自由未知数的个数;基础解系中解向量的个数;解空间的维数

求解方法

  1. 化行最简形(RREF)
    对系数矩阵A\boldsymbol{A}作初等行变换,化为 RREF。

  2. 确定主变量与自由变量

    • 主元列\to主变量
    • 非主元列\to自由变量
    • rank(A)=r\operatorname{rank}(A) = r,未知数个数为 nn,则自由变量个数为 nrn - r
    • r=nr = n,只有零解,无基础解系。
  3. 构造基础解系
    对每个自由变量:

    • 令该自由变量 = 1,其余自由变量 = 0;
    • 代入方程解出主变量,得到一个解向量。
  4. 结果
    得到 nrn - r 个线性无关解向量,即为基础解系。

eg:

求齐次线性方程组

{x1+x2x3x4=0,2x15x2+3x3+2x4=0,7x17x2+3x3+x4=0\begin{cases} x_1 + x_2 - x_3 - x_4 = 0, \\ 2x_1 - 5x_2 + 3x_3 + 2x_4 = 0, \\ 7x_1 - 7x_2 + 3x_3 + x_4 = 0 \end{cases}

的基础解系和通解。

:用初等行变换将方程组的系数矩阵化为行最简形:

A=(111125327731)(102/73/7015/74/70000)A = \begin{pmatrix} 1 & 1 & -1 & -1 \\ 2 & -5 & 3 & 2 \\ 7 & -7 & 3 & 1 \end{pmatrix} \sim \begin{pmatrix} 1 & 0 & -2/7 & -3/7 \\ 0 & 1 & -5/7 & -4/7 \\ 0 & 0 & 0 & 0 \end{pmatrix}

同解方程组为:

{x1=27x3+37x4,x2=57x3+47x4,\begin{cases} x_1 = \dfrac{2}{7}x_3 + \dfrac{3}{7}x_4, \\ x_2 = \dfrac{5}{7}x_3 + \dfrac{4}{7}x_4, \end{cases}

x3=k1, x4=k2x_3 = k_1,\ x_4 = k_2,代入得通解为:

(x1x2x3x4)=(27k1+37k257k1+47k2k1k2)=k1(2/75/710)+k2(3/74/701),(k1,k2R)\begin{pmatrix} x_1 \\ x_2 \\ x_3 \\ x_4 \end{pmatrix} = \begin{pmatrix} \dfrac{2}{7}k_1 + \dfrac{3}{7}k_2 \\ \dfrac{5}{7}k_1 + \dfrac{4}{7}k_2 \\ k_1 \\ k_2 \end{pmatrix} = k_1 \begin{pmatrix} 2/7 \\ 5/7 \\ 1 \\ 0 \end{pmatrix} + k_2 \begin{pmatrix} 3/7 \\ 4/7 \\ 0 \\ 1 \end{pmatrix}, \quad (k_1, k_2 \in \mathbb{R})

基础解系为

ξ1=(2/75/710),ξ2=(3/74/701).\xi_1 = \begin{pmatrix} 2/7 \\ 5/7 \\ 1 \\ 0 \end{pmatrix}, \quad \xi_2 = \begin{pmatrix} 3/7 \\ 4/7 \\ 0 \\ 1 \end{pmatrix}.

eg2:

Am×nBn×l=O\boldsymbol{A}_{m \times n} \boldsymbol{B}_{n \times l} = O,证明 r(A)+r(B)nr(\boldsymbol{A}) + r(\boldsymbol{B}) \leq n:
将矩阵 B\boldsymbol{B}0\boldsymbol{0} 按列分块,由 AB=O\boldsymbol{ AB = O } 可得:

A(b1,b2,,bl)=(0,0,,0)Abi=0(i=1,2,,l)\boldsymbol{A(b_1, b_2, \dots, b_l) = (0, 0, \dots, 0) \Rightarrow Ab_i = 0 \quad (i = 1,2,\dots,l)}

这说明 b1,b2,,bl\boldsymbol{b_1, b_2, \dots, b_l} 都是 $ Ax = 0 $ 的解向量。

  1. B=O\boldsymbol{B = O},则 r(A)+r(B)=r(A)nr(\boldsymbol{A}) + r(\boldsymbol{B}) = r(\boldsymbol{A}) \leq n
  2. BO\boldsymbol{B \neq O},则 Ax=0\boldsymbol{ Ax = 0 } 有非零解,从而有基础解系:

η1,η2,,ηnr,(r=r(A))\boldsymbol{\eta_1, \eta_2, \dots, \eta_{n - r}, \quad (r = r(A))}

b1,b2,,bl\boldsymbol{b_1, b_2, \dots, b_l} 可由 η1,η2,,ηnr\boldsymbol{\eta_1, \eta_2, \dots, \eta_{n - r}} 线性表示,从而:

r(B)=r(b1,b2,,bl)r(η1,η2,,ηnr)=nr(A)r(B) = r(b_1, b_2, \dots, b_l) \leq r(\eta_1, \eta_2, \dots, \eta_{n - r}) = n - r(A)

因此:

r(A)+r(B)nr(\boldsymbol{A}) + r(\boldsymbol{B}) \leq n

非齐次线性方程组的解

性质

  1. η1,η2\boldsymbol{\eta}_1, \boldsymbol{\eta}_2 都是 Ax=b\boldsymbol{Ax=b} 的解,则 η1η2\boldsymbol{\eta}_1 - \boldsymbol{\eta}_2 为其导出组 Ax=0\boldsymbol{Ax=0} 的解
  2. η\boldsymbol{\eta}Ax=b\boldsymbol{Ax=b} 的解,ξ\boldsymbol{\xi} 是对应的 Ax=0\boldsymbol{Ax=0} 的解,则 ξ+η\boldsymbol{\xi} + \boldsymbol{\eta} 仍为Ax=b\boldsymbol{Ax = b} 的解
  3. η1,η2,,ηn\boldsymbol{\eta}_1, \boldsymbol{\eta}_2, \cdots, \boldsymbol{\eta}_n 都是 Ax=b\boldsymbol{Ax=b} 的解,则k1η1+k2η2++knηn(k1+k2++kn=1)k_1\boldsymbol{\eta}_1 + k_2\boldsymbol{\eta}_2 + \cdots + k_n\boldsymbol{\eta}_n \quad (k_1 + k_2 + \cdots + k_n = 1)Ax=b\boldsymbol{Ax=b} 的解。
  4. Ax=b\boldsymbol{Ax=b} 的任一解都可表示成 ξ+η\boldsymbol{\xi} + \boldsymbol{\eta},其中 ξ\boldsymbol{\xi} 是其导出组Ax=0\boldsymbol{Ax=0} 的某个解,η\boldsymbol{\eta}Ax=bA\mathbf{x} = \mathbf{b} 的某个解。

结构

ξ1,ξ2,,ξnr\boldsymbol{\xi_{1},\xi_{2},\dots,\xi_{n-r}}Ax=0\boldsymbol{Ax=0}的基础解系,η\boldsymbol{\eta}^*Ax=b\boldsymbol{Ax=b}的一个特解,则Ax=b\boldsymbol{Ax=b}的通解是x=k1ξ1+k2ξ2++knrξnr+η\boldsymbol{x}=k_{1}\boldsymbol{\xi_{1}}+k_{2}\boldsymbol{\xi_{2}}+\dots+k_{n-r}\boldsymbol{\xi_{n-r}}+\boldsymbol{\eta}^*

即非齐次方程的通解=对应的齐次方程的通解+非齐次方程的特解

求解方法

  1. 将增广矩阵 (A,b)(\boldsymbol{A,b}) 做行变换化成行阶梯形,写出 r(A)r(\boldsymbol{A})r(A,b)r(\boldsymbol{A,b}),若 r(A)<r(A,b)r(\boldsymbol{A}) < r(\boldsymbol{A,b}),则方程组无解。
  2. r(A)=r(A,b)r(\boldsymbol{A}) = r(\boldsymbol{A,b}),继续将 (A,b)(\boldsymbol{A,b}) 化成行最简形。
  3. r(A)=r(A,b)=rr(\boldsymbol{A}) = r(\boldsymbol{A,b}) = r,将行最简形中 rr 个非零行的非零首元对应的未知数取作非自由未知数,其余 nrn - r 个未知数取作自由未知数,写出含有 nrn - r 个参数的通解。

eg1:

求解线性方程组

{x1x2x3+x4=0x1x2+x33x4=1x1x22x3+3x4=12\begin{cases} x_1 - x_2 - x_3 + x_4 = 0 \\ x_1 - x_2 + x_3 - 3x_4 = 1 \\ x_1 - x_2 - 2x_3 + 3x_4 = -\dfrac{1}{2} \end{cases}

对增广矩阵进行初等行变换:

(A,b)=(1111011131112312)(11011200121200000)(\boldsymbol{A,b}) = \begin{pmatrix} 1 & -1 & -1 & 1 & 0 \\ 1 & -1 & 1 & -3 & 1 \\ 1 & -1 & -2 & 3 & -\dfrac{1}{2} \end{pmatrix} \sim \begin{pmatrix} 1 & -1 & 0 & -1 & \dfrac{1}{2} \\ 0 & 0 & 1 & -2 & \dfrac{1}{2} \\ 0 & 0 & 0 & 0 & 0 \end{pmatrix}

因为 r(A)=r(A,b)=2<4r(\boldsymbol{A}) = r(\boldsymbol{A,b}) = 2 < 4,故方程组有无穷多解。
x2=x4=0x_2 = x_4 = 0 代入原方程组的同解方程组:

{x1=x2+x4+12x3=2x4+12\begin{cases} x_1 = x_2 + x_4 + \dfrac{1}{2} \\ x_3 = 2x_4 + \dfrac{1}{2} \end{cases}

得原方程组的一个特解为:

η=(12, 0, 12, 0)T\boldsymbol{\eta}^* = \left( \dfrac{1}{2},\ 0,\ \dfrac{1}{2},\ 0 \right)^T

再取 (x2x4)=(10),(01)\begin{pmatrix} x_2 \\ x_4 \end{pmatrix} = \begin{pmatrix} 1 \\ 0 \end{pmatrix}, \begin{pmatrix} 0 \\ 1 \end{pmatrix} 代入导出组的同解方程组:

{x1=x2+x4x3=2x4\begin{cases} x_1 = x_2 + x_4 \\ x_3 = 2x_4 \end{cases}

得导出组的基础解系为:

ξ1=(1,1,0,0)T,ξ2=(1,0,2,1)T\boldsymbol{\xi}_1 = (1, 1, 0, 0)^T,\quad \boldsymbol{\xi}_2 = (1, 0, 2, 1)^T

于是所求通解为:

x=k1(1,1,0,0)T+k2(1,0,2,1)T+(12, 0, 12, 0)T(k1,k2R)\boldsymbol{x} = k_1 (1, 1, 0, 0)^T + k_2 (1, 0, 2, 1)^T + \left( \dfrac{1}{2},\ 0,\ \dfrac{1}{2},\ 0 \right)^T \quad (k_1, k_2 \in \mathbb{R})

eg2:

试确定 λ\lambda 的值,使非齐次线性方程组

{(1+λ)x1+x2+x3=0x1+(1+λ)x2+x3=3x1+x2+(1+λ)x3=λ\begin{cases} (1 + \lambda)x_1 + x_2 + x_3 = 0 \\ x_1 + (1 + \lambda)x_2 + x_3 = 3 \\ x_1 + x_2 + (1 + \lambda)x_3 = \lambda \end{cases}

(1) 有唯一解,(2) 无解,(3) 有无穷多解,并求其通解。

方程组的系数行列式为:

A=1+λ1111+λ1111+λ=λ2(λ+3)|\boldsymbol{A}| = \begin{vmatrix} 1 + \lambda & 1 & 1 \\ 1 & 1 + \lambda & 1 \\ 1 & 1 & 1 + \lambda \end{vmatrix} = \lambda^2(\lambda + 3)

(1)A0|\boldsymbol{A}| \neq 0,即 λ0\lambda \neq 0λ3\lambda \neq -3 时,方程组有唯一解
(2)λ=0\lambda = 0 时,

(A,b)=(111011131110)(111000030000)(\boldsymbol{A,b}) = \begin{pmatrix} 1 & 1 & 1 & 0 \\ 1 & 1 & 1 & 3 \\ 1 & 1 & 1 & 0 \end{pmatrix} \sim \begin{pmatrix} 1 & 1 & 1 & 0 \\ 0 & 0 & 0 & 3 \\ 0 & 0 & 0 & 0 \end{pmatrix}

此时 r(A)=1r(\boldsymbol{A}) = 1r(A,b)=2r(\boldsymbol{A,b}) = 2,故方程组无解
(3)λ=3\lambda = -3 时,

(A,b)=(211012131123)(101101120000)(\boldsymbol{A,b}) = \begin{pmatrix} -2 & 1 & 1 & 0 \\ 1 & -2 & 1 & 3 \\ 1 & 1 & -2 & -3 \end{pmatrix} \sim \begin{pmatrix} 1 & 0 & -1 & -1 \\ 0 & 1 & -1 & -2 \\ 0 & 0 & 0 & 0 \end{pmatrix}

此时 r(A)=r(A,b)=2<3r(\boldsymbol{A}) = r(\boldsymbol{A,b}) = 2 < 3,方程组有无穷多解
同解方程组为:

{x1=x31x2=x32\begin{cases} x_1 = x_3 - 1 \\ x_2 = x_3 - 2 \end{cases}

x3=cx_3 = ccRc \in \mathbb{R}),得通解为:

{x1=c1x2=c2x3=c(cR)\begin{cases} x_1 = c - 1 \\ x_2 = c - 2 \\ x_3 = c \end{cases} \quad (c \in \mathbb{R})