L8-求解 Ax=b 简化行阶梯形式

linear-algebra

L8-求解 Ax=b 简化行阶梯形式

参考

这一节主要内容是 Ax=bAx=b 的求解步骤:

  1. 判断方程组的可解性 solvability 对于方程组 Am×nx=bA_{m \times n}x=b
系数矩阵的简化行阶梯形式解数量解形式
列满秩 r=n<mr=n<mR=[I0]R=\begin{bmatrix} I \\ 0 \end{bmatrix}1 个或 0 个x=xpx=x_{p}
行满秩 r=m<nr=m<nR=[IF]R=\begin{bmatrix} I & F \end{bmatrix}无限x=xp+xnx=x_{p}+x_{n}
方阵满秩 r=m=nr=m=nR=IR=I有且只有 1 个x=xpx=x_{p}
不满秩 r<mr<mr<nr<nR=[IF00]R=\begin{bmatrix} I & F \\ 0 & 0 \end{bmatrix}0 个或无限(由 bb 决定)x=xp+xnx=x_{p}+x_{n}
  1. 求出特解 particular solutions xpx_{p}:课程中所取的特解是在所有自由变量设为 00 的情况下,再求出主变量
  2. 求出零空间 N(A)N(A),作为基础解系 special solutions xnx_{n}
  3. 这样就可以得到方程组(当有解时)所有解的形式 x=xp+xnx={\color{Green}x_{p}}+{\color{Blue}x_{n}}

消元

已知方程组 Ax=bAx=b 如下

{x1+2x2+2x3+2x4=b12x1+4x2+6x3+8x4=b23x1+6x2+8x3+10x4=b3\begin{aligned} \left\{\begin{matrix} x_{1} + 2x_{2} + 2x_{3} + 2x_{4} = b_{1} \\ 2x_{1} + 4x_{2} + 6x_{3} + 8x_{4} = b_{2} \\ 3x_{1} + 6x_{2} + 8x_{3} + 10x_{4} = b_{3} \end{matrix}\right. \end{aligned}

写成矩阵的形式

Ax=[1222246836810][x1x2x3x4]=[b1b2b3]=b\begin{aligned} Ax = \begin{bmatrix} 1 & 2 & 2 & 2 \\ 2 & 4 & 6 & 8 \\ 3 & 6 &8 & 10 \end{bmatrix} \begin{bmatrix} x_{1}\\ x_{2}\\ x_{3}\\ x_{4} \end{bmatrix} & = \begin{bmatrix} b_{1}\\ b_{2}\\ b_{3} \end{bmatrix} & = b \end{aligned}
提示

观察可知 row3 = row1 + row2,即系数矩阵第三行与第一、二行线性相关,因此当方程组有解时,右侧应该满足 b3=b1+b2b_{3}=b_{1}+b_{2}

使用消元法求解方程组 Ax=bAx=b 需要构建增广矩阵

[AB]=[1222b12468b236810b3]E[1222b10024b22b10024b33b1]E[1222b10024b22b10000b3b2b1]\begin{aligned} \begin{bmatrix} AB \end{bmatrix} & = \begin{bmatrix} 1 & 2 & 2 & 2 & b_{1} \\ 2 & 4 & 6 & 8 & b_{2} \\ 3 & 6 & 8 & 10 & b_{3} \end{bmatrix} \xrightarrow{E} \begin{bmatrix} {\color{Red}1 } & 2 & 2 & 2 & b_{1} \\ 0 & 0 & 2 & 4 & b_{2}-2b_{1} \\ 0 & 0 & 2 & 4 & b_{3}-3b_{1} \end{bmatrix} \xrightarrow{E} \begin{bmatrix} {\color{Red}1 } & 2 & 2 & 2 & b_{1} \\ 0 & 0 & {\color{Red}2 } & 4 & b_{2}-2b_{1} \\ 0 & 0 & 0 & 0 & b_{3}-b_{2}-b_{1} \end{bmatrix} \end{aligned}

当方程组有解时(三个等式都满足),由最后一行得到 0=b3b2b10=b_{3}-b_{2}-b_{1} 成立,即 b3=b1+b2b_{3}=b_{1}+b_{2}

若取

[b1b2b3]=[156]\begin{aligned} \begin{bmatrix} b_{1} \\ b_{2} \\ b_{3} \end{bmatrix} & = \begin{bmatrix} 1 \\ 5 \\ 6 \end{bmatrix} \end{aligned}

则增广矩阵为

[AB]E[122210024300000]\begin{aligned} \begin{bmatrix} AB \end{bmatrix} \xrightarrow{E} \begin{bmatrix} {\color{Red}1 } & 2 & 2 & 2 & 1 \\ 0 & 0 & {\color{Red}2 } & 4 & 3 \\ 0 & 0 & 0 & 0 & 0 \end{bmatrix} \end{aligned}
提示

其中第一列、第三列是主列 pivot column

可解性

使用消元法求解方程组 Am×nx=bA_{m \times n}x=b ,其可解性 solvability 可以从方程组右侧的 bb 向量考虑

  • 列角度:将 bb 看作是系数矩阵 AA 各列基于 xx 的线性组合,如果这个角度考虑,那么当方程组有解时,则 bb 必须在系数矩阵 AA 所构成的列空间 C(A)C(A)
  • 行角度:若系数矩阵 AA 在消元过程中存在某行各元素均为零,则相应地在增广矩阵中该行的 bb 的组合也要是零

为了便于讨论方程组的可解性,引入秩 rr 这一概念(它是指在阶梯形式的矩阵中,矩阵的主元数量

{rmrn\begin{aligned} \left\{\begin{matrix} r \le m \\ r \le n \end{matrix}\right. \end{aligned}

满秩 full rank 指各行或各列上都有主元,对于矩阵 Am×nA_{m \times n} 有两种情况:

  • r=nr=n 时,列满秩
  • r=mr=m 时,行满秩

当秩满足不同条件时,方程组有不同的解。

列满秩

当列满秩 r=n<mr=n<m,每一列都有一个主元,则没有自由变量,即自由变量的自由度为 00,即零空间的维度是零(几何形式是一个点)

例如对于方程组 Ax=bAx=b

[13216151][x1x2]=[b1b2b3b4]\begin{aligned} \begin{bmatrix} 1 & 3 \\ 2 & 1 \\ 6 & 1 \\ 5 & 1 \end{bmatrix} \begin{bmatrix} x_{1} \\ x_{2} \end{bmatrix} &= \begin{bmatrix} b_{1} \\ b_{2} \\ b_{3} \\ b_{4} \end{bmatrix} \end{aligned}

消元法求解方程组(为了演示简便,增广矩阵省略 bb 部分)

A=[13216151]E[1305017014]E[13050000]E[10010000]=R\begin{aligned} A= \begin{bmatrix} 1 & 3 \\ 2 & 1 \\ 6 & 1 \\ 5 & 1 \end{bmatrix} \xrightarrow{E} \begin{bmatrix} 1 & 3 \\ 0 & -5 \\ 0 & -17 \\ 0 & -14 \end{bmatrix} \xrightarrow{E} \begin{bmatrix} 1 & 3 \\ 0 & -5 \\ 0 & 0 \\ 0 & 0 \end{bmatrix} \xrightarrow{E} \begin{bmatrix} 1 & 0 \\ 0 & 1 \\ 0 & 0 \\ 0 & 0 \end{bmatrix} &=R \end{aligned}

通过系数矩阵的简化行阶梯形式可知,该方程组是否有解(唯一解),取决于变换后的增广矩阵第三、第四行是否成立,这与 bb 相关。

所以当系数矩阵是列满秩 r=n<mr=n<m 时,方程组 Ax=bAx=b 有 0 或 1 个解。如果存在唯一解 unique solution,其形式是 x=xpx=x_{p}

行满秩

当行满秩 r=m<nr=m<n,每一行都有一个主元,则自由变量的自由度为 nr=nmn-r=n-m

例如方程组 Ax=bAx=b

[12653111][x1x2x3x4]=[b1b2]\begin{aligned} \begin{bmatrix} 1 & 2 & 6 & 5 \\ 3 & 1 & 1 & 1 \end{bmatrix} \begin{bmatrix} x_{1} \\ x_{2} \\ x_{3} \\ x_{4} \end{bmatrix} &= \begin{bmatrix} b_{1} \\ b_{2} \end{bmatrix} \end{aligned}

消元法求解方程组(为了演示简便,增广矩阵省略 b 部分),在系数矩阵的简化行阶梯形式中 fnf_{n} 表示自由变量所对应的系数

A=[12653111]E[10f1f201f3f4]=R\begin{aligned} A = \begin{bmatrix} 1 & 2 & 6 & 5 \\ 3 & 1 & 1 & 1 \end{bmatrix} \xrightarrow{E} \begin{bmatrix} 1 & 0 & f_{1} & f_{2} \\ 0 & 1 & f_{3} & f_{4} \end{bmatrix} &=R \end{aligned}

根据 L7-求解 Ax=0 主变量和特解 这一节可知方程组 Ax=0Ax=0 的零空间是 N=[FI]N=\begin{bmatrix}-F \\ I\end{bmatrix},(通过确定自由变量的值后)可以方便地得到 Ax=0Ax=0 所有解的形式是 xnx_{n}

此时方程组 Ax=bAx=b 必有解,且是无限个,其形式是 x=xp+xnx=x_{p}+x_{n}

其中 xnx_{n} 表示方程组 Ax=0Ax=0 的解所构成零空间 N(A)N(A),它作为基础解系,只需要再结合方程组 Ax=bAx=b 其中任意一个特解 xpx_{p}(可以将该特解的作用看作是一个平移向量,将方程组 Ax=0Ax=0 的解所构成零空间 N(A)N(A) 平移到方程组 Ax=bAx=b 的解构成的向量空间)

方程组的解所构成的空间
方程组的解所构成的空间

说明

由于行满秩,在简化行阶梯形式中系数矩阵不存在一整行的元素全为 00 的情况,而由于自由变量是可变动的,因此对于方程组 Ax=bAx=b 无论 bb 取任意值都有解

方阵满秩

当方程组的系数矩阵是一个方阵时,且列满秩(同时也就行满秩),即 r=n=mr=n=m

例如方程组 Ax=bAx=b

[1231][x1x2]=[b1b2]\begin{aligned} \begin{bmatrix} 1 & 2 \\ 3 & 1 \end{bmatrix} \begin{bmatrix} x_{1} \\ x_{2} \end{bmatrix} & = \begin{bmatrix} b_{1} \\ b_{2} \end{bmatrix} \end{aligned}

消元法求解方程组

说明

为了演示简便,增广矩阵省略 bb 部分

A=[1231]E[1001]=R\begin{aligned} A = \begin{bmatrix} 1 & 2 \\ 3 & 1 \end{bmatrix} \xrightarrow{E} \begin{bmatrix} 1 & 0 \\ 0 & 1 \end{bmatrix} &=R \end{aligned}

此时方程组 Ax=bAx=b 有且只有 1 个解,其形式是 x=xpx=x_{p}

不满秩

当不满秩 r<mr<mr<nr<n 时,(类似于行满秩和列满秩的混合情况)此时方程组 Ax=bAx=b 有 0 个或无数个解,可解性由 bb 决定。当有无数个解时,其形式是 x=xp+xnx=x_{p}+x_{n}

求解方程组

在前面消元那一小节中,求解的方程组 Ax=bAx=b 得到了阶梯形式

[AB]E[122210024300000]\begin{aligned} \begin{bmatrix} AB \end{bmatrix} \xrightarrow{E} \begin{bmatrix} {\color{Red}1 } & 2 & 2 & 2 & 1 \\ 0 & 0 & {\color{Red}2 } & 4 & 3 \\ 0 & 0 & 0 & 0 & 0 \end{bmatrix} \end{aligned}

所以该矩阵并不满秩

根据前一小节对可解性的讨论可知,在不满秩情况下,方程组是否有解由 bb 决定

观察这个示例的增广矩阵,它的最后一行的 bb 可使得等式成立,因此该方程组有无数个解

  • 求出特解 xpx_{p}
    令所有自由变量为 00 {x2=0x4=0\begin{aligned} \left\{\begin{matrix} x_{2} = 0 \\ x_{4} = 0 \end{matrix}\right. \end{aligned}
    回代入增广矩阵可得 {x1+2x3=12x3=3\begin{aligned} \left\{\begin{matrix} x_{1}+2x_{3}=1 \\ 2x_{3}=3 \end{matrix}\right. \end{aligned}
    可得其中一个特解 xpx_{p} {x1=2x2=0x3=32x4=0\begin{aligned} \left\{\begin{matrix} x_{1} = -2 \\ x_{2} = 0 \\ x_{3} = \frac{3}{2} \\ x_{4} = 0 \end{matrix}\right. \end{aligned}
  • 以方程组 Ax=0Ax=0 的零空间 N(A)N(A) 作为基础解系
    求出系数矩阵的简化行阶梯形式 AE[122200240000]E[120200120000]=R\begin{aligned} A \xrightarrow{E} \begin{bmatrix} {\color{Red}1 } & 2 & 2 & 2 \\ 0 & 0 & {\color{Red}2 } & 4 \\ 0 & 0 & 0 & 0 \end{bmatrix} \xrightarrow{E} \begin{bmatrix} {\color{Red}1 } & 2 & 0 & -2 \\ 0 & 0 & {\color{Red}1 } & 2 \\ 0 & 0 & 0 & 0 \end{bmatrix} & = R \end{aligned}
    根据 Rx=[IF00][xpivotxfree]=0N=[FI]\begin{aligned} Rx = \begin{bmatrix} I & F \\ 0 & 0 \\ \end{bmatrix} \begin{bmatrix} x_{pivot} \\ x_{free} \\ \end{bmatrix} & = 0 \Rightarrow N= \begin{bmatrix} -F \\ I \end{bmatrix} \end{aligned}
    先将简化行阶梯形式 RR 进行改写 R=[102201020000]\begin{aligned} R^{'}= \begin{bmatrix} {\color{Red}1 } & 0 & 2 & -2 \\ 0 & {\color{Red}1 } & 0 & 2 \\ 0 & 0 & 0 & 0 \end{bmatrix} \end{aligned}
    Rx=[102201020000][x1x3x2x4]=0\begin{aligned} R^{'}x= \begin{bmatrix} {\color{Red}1 } & 0 & 2 & -2 \\ 0 & {\color{Red}1 } & 0 & 2 \\ 0 & 0 & 0 & 0 \end{bmatrix} \begin{bmatrix} x_{1} \\ x_{3} \\ x_{2} \\ x_{4} \end{bmatrix} & = 0 \end{aligned}
    可得 N=[FI]=[22021001]\begin{aligned} N = \begin{bmatrix} -F \\ I \end{bmatrix} & = \begin{bmatrix} -2 & 2 \\ 0 & -2 \\ 1 & 0 \\ 0 & 1 \end{bmatrix} \end{aligned}
    因此方程组 Ax=0Ax=0 的两个特解是(相应地需要调换位置以对应原向量中 x2x_{2}x3x_{3} 的顺序) xp1=[2100],xp2=[2020]\begin{aligned} x_{p1} = \begin{bmatrix} -2 \\ {\color{Green}1 } \\ {\color{Green}0 } \\ 0 \end{bmatrix} , x_{p2} = \begin{bmatrix} -2 \\ {\color{Green}0 } \\ {\color{Green}-2 } \\ 0 \end{bmatrix} \end{aligned}
    因此方程组 Ax=0Ax=0 的解所构成的零空间是 N(A)=α[2100]+β[2020]\begin{aligned} N(A) = \alpha \begin{bmatrix} -2 \\ {\color{Green}1 } \\ {\color{Green}0 } \\ 0 \end{bmatrix} & + \beta \begin{bmatrix} -2 \\ {\color{Green}0 } \\ {\color{Green}-2 } \\ 0 \end{bmatrix} \end{aligned}
    提示

    方程组 Ax=0Ax=0 的零空间维度为 2,需要两个特解张成

将方程组 Ax=bAx=b 一个特定的解 xpx_{p} 和方程组 Ax=0Ax=0 的所有解 xnx_{n} 结合,就是方程组 Ax=bAx=b 的通解,即满足以下等式

A(xp+xn)=b\begin{aligned} A(x_{p}+x_{n}) = b \end{aligned}

即对于方程组 Ax=bAx=b, 其中 b=[156]\begin{aligned} b = \begin{bmatrix} 1 \\ 5 \\ 6 \end{bmatrix}\end{aligned},它有无数个解,其任意解/通解的形式是

x=xparticular+xnull=[20320]+c1[2100]+c2[2020]\begin{aligned} x & = x_{particular} + x_{null} \\ & = \begin{bmatrix} -2 \\ 0 \\ \frac{3}{2} \\ 0 \end{bmatrix} + c_{1} \begin{bmatrix} -2 \\ 1 \\ 0 \\ 0 \end{bmatrix} + c_{2} \begin{bmatrix} -2 \\ 0 \\ -2 \\ 0 \end{bmatrix} \end{aligned}

Copyright © 2024 Ben

Theme BlogiNote

Icons from Icônes