线性方程组的求解问题

本节内容可以联系之前的代数知识理解,因为本质上方程组问题就是向量组问题

该部分的完整解题流程是:

  1. 根据有解的条件判断是否有解
  2. 在有解的情况下确定通解的结构
  3. 解出来,且指明 k 1 . . . k s k_1...k_s k1...ks为任意常数

这部分内容涉及基础解系特解两方面,其中特解是非齐次线性方程组涉及的

齐次

对于齐次线性方程组 A x = 0 Ax=0 Ax=0,矩阵 A m × n A_{m×n} Am×n意味着有 m m m行个方程式子和 n n n个未知数。 x x x是一个具有 n n n个分量的列矩阵

其通解为 k 1 ξ 1 + k 2 ξ 2 + . . . + k s ξ s k_1\xi_1+k_2\xi_2+...+k_s\xi_s k1ξ1+k2ξ2+...+ksξs,其中 ξ s \xi_s ξs是一个 n n n维的列向量, s s s取决于自由度。

因为矩阵 A m × n A_{m×n} Am×n在进行初等行变换之后会化为一个秩 r ( A ) = r ≤ n r(A)=r\le n r(A)=rn的矩阵, r r r是独立方程(也称为约束)的个数,也就是该方程组能确定 n n n维中的 r r r个维度,但是仍然剩余 s = n − r s=n-r s=nr个维度不能被约束,所不能被约束的维度会张成一个 s s s维空间,即解的空间。


※对于非齐次方程组 A x = b Ax=b Ax=b的某两个解 η 1 \eta_1 η1 η 2 \eta_2 η2,作差 η 1 − η 2 \eta_1-\eta_2 η1η2就是齐次方程组 A x = 0 Ax=0 Ax=0的一个解

因为 A η 1 = 0 A\eta_1=0 Aη1=0 A η 2 = 0 A\eta_2=0 Aη2=0作差可以得到,类似地也有 A ( η 1 + η 2 ) = 2 b A(\eta_1+\eta_2)=2b A(η1+η2)=2b


【例】对于下面的方程组

{ x 1 + x 2 − 3 x 4 − x 5 = 0 x 1 − x 2 + 2 x 3 − x 4 = 0 4 x 1 − 2 x 2 + 6 x 3 + 3 x 4 − 4 x 5 = 0 2 x 1 + 4 x 2 − 2 x 3 + 4 x 4 − 7 x 5 = 0 \begin{cases} x_1+x_2-3x_4-x_5=0\\ x_1-x_2+2x_3-x_4=0\\ 4x_1-2x_2+6x_3+3x_4-4x_5=0\\ 2x_1+4x_2-2x_3+4x_4-7x_5=0\\ \end{cases} x1+x23x4x5=0x1x2+2x3x4=04x12x2+6x3+3x44x5=02x1+4x22x3+4x47x5=0
能够写成矩阵形式并化简

A = ∣ 1 1 0 − 1 − 1 1 − 1 2 − 1 0 4 − 2 6 3 − 4 2 4 − 2 4 − 7 ∣ → 初等行变换化简为 ∣ 1 1 0 − 3 − 1 0 − 2 2 2 1 0 0 0 3 − 1 0 0 0 0 0 ∣ A=\left| \begin{matrix} 1 & 1 & 0 & -1 & -1 \\ 1 & -1 & 2 & -1 & 0 \\ 4 & -2 & 6 & 3 & -4 \\ 2 & 4 & -2 & 4 & -7 \\ \end{matrix} \right| \xrightarrow {初等行变换化简为} \left| \begin{matrix} 1 & 1 & 0 & -3 & -1 \\ 0 & -2 & 2 & 2 & 1 \\ 0 & 0 & 0 & 3 & -1 \\ 0 & 0 & 0 & 0 & 0 \\ \end{matrix} \right| A= 11421124026211341047 初等行变换化简为 10001200020032301110

因为该阶梯矩阵秩为3,所以任意找出三个列组成秩为三的子矩阵即可,在这里选取一、二、四列。由这些列可以唯一地确定三维空间内的解(即确定五个未知量中的三个),但是由于空间是五维的(未知量的个数有五个),所以仍有两个维度不能确定,它们属于自由未知量。故取剩余第三、五列元素 x 1 x_1 x1 x 3 x_3 x3设为自由未知量,令 x 3 = k 1 x_3=k_1 x3=k1 x 5 = k 2 x_5=k_2 x5=k2,则 A x = 0 Ax=0 Ax=0即下式:

∣ 1 1 0 − 3 − 1 0 − 2 2 2 1 0 0 0 3 − 1 0 0 0 0 0 ∣ ∣ x 1 x 2 k 1 x 4 k 2 ∣ = 0 \left| \begin{matrix} 1 & 1 & 0 & -3 & -1 \\ 0 & -2 & 2 & 2 & 1 \\ 0 & 0 & 0 & 3 & -1 \\ 0 & 0 & 0 & 0 & 0 \\ \end{matrix} \right| \left| \begin{matrix} x_1 \\ x_2\\ k_1\\ x_4\\ k_2 \end{matrix} \right|=0 10001200020032301110 x1x2k1x4k2 =0
根据题意解的结构肯定是 k 1 ξ 1 + k 2 ξ 2 k_1\xi_1+k_2\xi_2 k1ξ1+k2ξ2则它大概长底下这样
ξ 1 = ( □ □  1 ⋅ k 1  □  0 ⋅ k 2   ) ξ 2 = ( □ □  0 ⋅ k 1  □  1 ⋅ k 2   ) \begin{matrix} \xi_1=(□\text{ }□\text{ }1·k_1\text{ }□\text{ }0·k_2\text{ })\\ \xi_2=(□\text{ }□\text{ }0·k_1\text{ }□\text{ }1·k_2\text{ }) \end{matrix} ξ1=(  1k1  0k2 )ξ2=(  0k1  1k2 )

其中 □ □ 是被约束的维度(即 x 1 x_1 x1 x 2 x_2 x2 x 4 x_4 x4),先不看他,主要是三和五这两个。因为基础解析的要求是线性无关,所以这俩位置一个填0另一个填1是最简单的形式(如上),但是未必是最优的,因为不一定方便计算,至少在本题用下面这种更方便:

ξ 1 = ( □ □  1 ⋅ k 1  □  0 ⋅ k 2   ) ξ 2 = ( □ □  0 ⋅ k 1  □  3 ⋅ k 2   ) \begin{matrix} \xi_1=(□\text{ }□\text{ }1·k_1\text{ }□\text{ }0·k_2\text{ })\\ \xi_2=(□\text{ }□\text{ }0·k_1\text{ }□\text{ }3·k_2\text{ }) \end{matrix} ξ1=(  1k1  0k2 )ξ2=(  0k1  3k2 )

随后解出 x 1 x_1 x1 x 2 x_2 x2 x 4 x_4 x4,即用 k 1 k_1 k1 k 2 k_2 k2表示 x 1 x_1 x1 x 2 x_2 x2 x 4 x_4 x4。按下述方法整理即可

∣ x 1 x 2 x 3 x 4 x 5 ∣ = ∣ − k 1 + 2 7 k 2 k 1 + 5 2 k 2 k 1 k 2 3 k 2 ∣ = k 1 ∣ − 1 1 1 0 0 ∣ + k 2 ∣ 2 7 5 2 0 1 3 ∣ \left| \begin{matrix} x_1 \\ x_2\\ x_3\\ x_4\\ x_5 \end{matrix} \right|=\left| \begin{matrix} -k_1+\frac27k_2 \\ k_1+\frac52k_2\\ k_1\\ k_2\\ 3k_2 \end{matrix} \right|=k_1\left| \begin{matrix} -1\\1\\1\\0\\0 \end{matrix} \right|+k_2\left| \begin{matrix} \frac27\\ \frac52\\0\\1\\3 \end{matrix} \right| x1x2x3x4x5 = k1+72k2k1+25k2k1k23k2 =k1 11100 +k2 7225013

其中 k 1 k_1 k1 k 2 k_2 k2是任意常数

非齐次

非齐次线性方程组的求解和上面的齐次形式类似,但是其解的结构是 k 1 ξ 1 + k 2 ξ 2 + . . . + k s ξ s + η k_1\xi_1+k2\xi_2+...+k_s\xi_s+\eta k1ξ1+k2ξ2+...+ksξs+η,其中 s = n − r ( A ) s=n-r(A) s=nr(A) η \eta η是一个特解

并且矩阵应当写为增广矩阵的形式,即对于方程 A x = b Ax=b Ax=b,要处理的矩阵是 [ A ∣ b ] [A|b] [Ab]

猜你喜欢

转载自blog.csdn.net/qq_39377889/article/details/129378255