LA@非齐次线性方程组解的结构
文章目录
非齐次线性方程组解的结构
非齐次线性方程组
- 本文讨论向量方程形式的非齐次线性方程组 A x = b \bold{Ax=b} Ax=b (1),并记方程 ( 1 ) (1) (1)的同系数矩阵的齐次线性方程组: ( 1 ) \bold{(1)} (1) (2)
导出组
- 将方程组 ( 1 ) (1) (1)常数项全置为0得到 ( 2 ) (2) (2),则 ( 2 ) (2) (2)是 ( 1 ) (1) (1)的导出组,也称为对应的齐次线性方程组
导出组的解和被导出组的解性质👺
- 设 ξ 1 , ξ 2 \xi_1,\xi_2 ξ1,ξ2是 ( 1 ) (1) (1)的两个解,则 ξ 1 − ξ 2 \xi_1-\xi_2 ξ1−ξ2是 ( 2 ) (2) (2)的解
- 若 ξ \xi ξ是 ( 2 ) (2) (2)的解,则 ξ 1 + ξ \xi_1+\xi ξ1+ξ是 ( 1 ) (1) (1)的解
证明
- A ξ 1 = b ; A ξ 2 = b \bold{A\xi_1=b};\bold{A\xi_2=b} Aξ1=b;Aξ2=b
- 显然, A ξ 1 − A ξ 2 = 0 A\xi_1-A\xi_2=0 Aξ1−Aξ2=0
- 所以 A ( ξ 1 − ξ 2 ) = 0 \bold{A(\xi_1-\xi_2)=0} A(ξ1−ξ2)=0, ξ 1 − ξ 2 \xi_1-\xi_2 ξ1−ξ2是 ( 2 ) (2) (2)的解
- A ( ξ 1 + ξ ) \bold{A(\xi_1+\xi)} A(ξ1+ξ)= b + 0 \bold{b+0} b+0,因此 ξ 1 + ξ \xi_1+\xi ξ1+ξ是 ( 1 ) (1) (1)的解
非齐次线性方程组的特解
- 设 ξ ∗ \xi^* ξ∗是方程 ( 1 ) (1) (1)的一个解,这个不能表示其他解的解称为特解
非齐次线性方程组的通解结构
-
D
:
ξ
1
,
ξ
2
,
⋯
,
ξ
s
D:\xi_1,\xi_2,\cdots,\xi_s
D:ξ1,ξ2,⋯,ξs是
(
2
)
(2)
(2)的基础解系,则
ξ
∗
+
∑
i
s
k
i
ξ
i
\xi^*+\sum\limits_{i}^{s}k_i\xi_i
ξ∗+i∑skiξi是方程
(
1
)
(1)
(1)的通解
- s = n − r s=n-r s=n−r, r = R ( A ) r=R(A) r=R(A), n n n是未知数的个数
- 其中 k 1 , ⋯ , k s k_1,\cdots,k_s k1,⋯,ks为任意常数
- 文字描述就是:非齐次线性方程的通解等于导出组的通解加上一个非齐次线性方程组本身的一个特解
证明
- 设
(
1
)
(1)
(1)的任意一个解为
x
x
x,则
x
,
ξ
∗
x,\xi^*
x,ξ∗都是
(
1
)
(1)
(1)的两个特解
- 由性质1: ξ = x − ξ ∗ \xi=x-\xi^* ξ=x−ξ∗是 ( 2 ) (2) (2)的一个解,则它可以用 D D D线性表示: ξ = ∑ i s k i ξ i \xi=\sum\limits_{i}^{s}k_i\xi_i ξ=i∑skiξi
- 所以 ( 1 ) (1) (1)的任意解 x x x可以表示为 x = ξ + ξ ∗ = ∑ i s k i ξ i + ξ ∗ x=\xi+\xi^*=\sum\limits_{i}^{s}k_i\xi_i+\xi^* x=ξ+ξ∗=i∑skiξi+ξ∗
- 也即是说 ξ ∗ + ∑ i s k i ξ i \xi^*+\sum\limits_{i}^{s}k_i\xi_i ξ∗+i∑skiξi,能够表示任何 ( 1 ) (1) (1)的解,因此是 ( 1 ) (1) (1)的通解
推论:导出组的解和原方程组解的关系
- 当 ( 1 ) (1) (1)有解时,“ ( 1 ) (1) (1)只有唯一解"的充要条件是” ( 2 ) (2) (2)只有唯一解".
- 该命题的逆否命题为:当 ( 1 ) (1) (1)有解时,“ ( 2 ) (2) (2)有多解"充要条件是” ( 1 ) (1) (1)具有多解"
- 证法1:
- 充分性
- 如果 ( 1 ) (1) (1)具有多解 ξ 1 , ξ 2 , ⋯ \xi_1,\xi_2,\cdots ξ1,ξ2,⋯, ξ i ≠ ξ j , i ≠ j \xi_{i}\neq{\xi_j},i\neq{j} ξi=ξj,i=j则 ξ = ξ 1 − ξ 2 \xi=\xi_1-\xi_2 ξ=ξ1−ξ2是 ( 2 ) (2) (2)的解,
- 因为 ξ 1 ≠ ξ 2 \xi_1\neq\xi_2 ξ1=ξ2,所以 ξ ≠ 0 \xi\neq{0} ξ=0,因此 ( 2 ) (2) (2)具有多解
- 必要性:
- 若 ( 2 ) (2) (2)有多解: ξ 1 , ξ 2 , ⋯ \xi_1,\xi_2,\cdots ξ1,ξ2,⋯,且 ξ ∗ \xi^* ξ∗是 ( 1 ) (1) (1)的一个特解
- 令 ξ 1 + ξ ∗ \xi_1+\xi^* ξ1+ξ∗以及 ξ 2 + ξ ∗ \xi_2+\xi^* ξ2+ξ∗是 ( 1 ) (1) (1)的两个不同的解
- 即 ( 1 ) (1) (1)有多解
- 充分性
- 证法2:
- 由于 ( 1 ) (1) (1)有解,所以 R ( A , b ) = R ( A ) = r R(\bold{A,b})=R(\bold{A})=r R(A,b)=R(A)=r
- 若 r < n r<n r<n则 ( 1 , 2 ) (1,2) (1,2)都有多解
- 若 r = n r=n r=n,则 ( 1 , 2 ) (1,2) (1,2)都只有唯一解
求非齐次线性方程组通解的步骤
将增广矩阵化为行最简形矩阵
- 将增广矩阵 B = ( A ∣ b ) \bold{B=(A|b)} B=(A∣b)在化简为行简化矩阵: B ~ = ( A ~ ∣ b ~ ) \bold{\widetilde{B}=(\widetilde{A}|\widetilde{b})} B =(A ∣b )
读出特解和通解
-
从 B ~ \widetilde{\bold{B}} B 中读出非自由未知数关于自由未知数的表达式(通解)
-
为 n − r n-r n−r个自由未知数取一组值(通常为全部取0最为方便),计算出一个(1)的一个特解
-
再由 B ~ \widetilde{\bold{B}} B 求出 ( 2 ) (2) (2)的通解
-
( 2 ) (2) (2)的通解加上 ( 1 ) (1) (1)的特解得到 ( 1 ) (1) (1)的通解
例
-
A x = b \bold{Ax=b} Ax=b
- A = ( 1 − 1 0 1 − 1 2 0 1 0 − 1 3 − 1 − 1 − 1 − 1 ) ; b = ( 1 2 0 ) \bold A=\begin{pmatrix} 1& -1& 0& 1& -1 \\ 2& 0& 1& 0& -1 \\ 3& -1& -1& -1& -1 \\ \end{pmatrix}; \bold b=\begin{pmatrix} 1 \\ 2 \\ 0 \\ \end{pmatrix} A= 123−10−101−110−1−1−1−1 ;b= 120
-
增广矩阵初等行变换为行最简形矩阵:
- A ‾ = ( 1 − 1 0 1 − 1 1 2 0 1 0 − 1 2 3 − 1 − 1 − 1 − 1 0 ) ∼ r ( 1 0 0 − 1 2 − 1 4 1 4 0 1 0 − 3 2 3 4 − 3 4 0 0 1 1 − 1 2 3 2 ) = U \overline{A}=\begin{pmatrix} 1& -1& 0& 1& -1& 1 \\ 2& 0& 1& 0& -1& 2 \\ 3& -1& -1& -1& -1& 0 \\ \end{pmatrix} \overset{r}{\sim} \begin{pmatrix} 1& 0& 0& -\frac{1}{2}& -\frac{1}{4}& \frac{1}{4} \\ 0& 1& 0& -\frac{3}{2}& \frac{3}{4}& -\frac{3}{4} \\ 0& 0& 1& 1& -\frac{1}{2}& \frac{3}{2} \\ \end{pmatrix} =U A= 123−10−101−110−1−1−1−1120 ∼r 100010001−21−231−4143−2141−4323 =U
-
容易看出,非自由未知数由 r = 3 r=3 r=3个,由 n − r = 5 − 3 = 2 n-r=5-3=2 n−r=5−3=2个自由未知数决定
-
为简单起见,我们将自由未知数 x 4 , x 5 x_4,x_5 x4,x5全部取0来求得一个特解
- x 1 = 1 2 x 4 + 1 4 x 5 + 1 4 x_1=\frac{1}{2}x_4+\frac{1}{4}x_5+\frac{1}{4} x1=21x4+41x5+41
- x 2 = 3 2 x 4 − 3 4 x 5 − 3 4 x_2=\frac{3}{2}x_4-\frac{3}{4}x_5-\frac{3}{4} x2=23x4−43x5−43
- x 3 = − x 4 + 1 2 x 5 + 3 2 x_3=-x_4+\frac{1}{2}x_5+\frac{3}{2} x3=−x4+21x5+23
- 走捷径就是直接构造用 U U U的最后一列再追加 r r r个0构造 ξ ∗ \xi^* ξ∗
-
从而该非齐次线性方程(1)的一个最简单特解就是
- ξ ∗ = ( 1 4 − 3 4 3 2 0 0 ) \xi^*=\begin{pmatrix} \frac{1}{4}\\ -\frac{3}{4}\\ \frac{3}{2}\\ 0\\ 0 \end{pmatrix} ξ∗= 41−432300
-
再求导出组,即方程 ( 2 ) (2) (2)的一个基础解系
-
A A A的行最简形矩阵 A ~ \widetilde{\bold{A}} A 包含在 U U U中
-
A ~ = ( 1 0 0 − 1 2 − 1 4 0 1 0 − 3 2 3 4 0 0 1 1 − 1 2 ) = ( E 3 , T 3 × 2 ) \widetilde{\bold A}=\begin{pmatrix} 1& 0& 0& -\frac{1}{2}& -\frac{1}{4} \\ 0& 1& 0& -\frac{3}{2}& \frac{3}{4} \\ 0& 0& 1& 1& -\frac{1}{2} \\ \end{pmatrix} =(E_3,T_{3\times{2}}) A = 100010001−21−231−4143−21 =(E3,T3×2)
-
因此可以立即得到 ( 2 ) (2) (2)的基础解系
- ξ 1 = ( 1 2 3 2 − 1 1 0 ) ; ξ 2 = ( 1 4 − 3 4 1 2 0 1 ) \xi_1=\begin{pmatrix} \frac{1}{2} \\ \frac{3}{2} \\ -1 \\ 1\\ 0 \end{pmatrix}; \xi_2=\begin{pmatrix} \frac{1}{4}\\ -\frac{3}{4}\\ \frac{1}{2}\\ 0\\ 1 \end{pmatrix} ξ1= 2123−110 ;ξ2= 41−432101
-
因此 ( 1 ) (1) (1)的通解为:
- ξ = ξ ∗ + ∑ i = 1 2 k i ξ i k i ( i = 1 , 2 ) 为任意常数 \xi=\xi^*+\sum\limits_{i=1}^{2}k_i\xi_i \\k_i(i=1,2)为任意常数 ξ=ξ∗+i=1∑2kiξiki(i=1,2)为任意常数
-
实际上为例书写紧凑,可以写成行向量转置的形式:
- ( ξ 1 , ξ 2 ) = ( − T T , E ) T = ( 1 2 3 2 − 1 1 0 1 4 − 3 4 1 2 0 1 ) T ξ 1 = ( 1 2 , 3 2 , − 1 , 1 , 0 ) T ξ 2 = ( 1 4 , − 3 4 , 1 2 , 0 , 1 ) T (\xi_1,\xi_2)=(-T^T,E)^T= \begin{pmatrix} \frac{1}{2}&\frac{3}{2}&-1 &1&0 \\ \frac{1}{4}&-\frac{3}{4}&\frac{1}{2} &0&1 \end{pmatrix}^T \\ \xi_1=(\frac{1}{2},\frac{3}{2},-1,1,0)^T \\ \xi_{2}=(\frac{1}{4},-\frac{3}{4},\frac{1}{2},0,1)^T (ξ1,ξ2)=(−TT,E)T=(214123−43−1211001)Tξ1=(21,23,−1,1,0)Tξ2=(41,−43,21,0,1)T
-
例👺
-
设 A x = b \bold{Ax=b} Ax=b的增广矩阵为 B B B
B = ( 1 − 1 − 1 1 0 1 − 1 1 − 3 1 1 − 1 − 2 3 − 1 2 ) ∼ r ( 1 − 1 0 − 1 1 2 0 0 1 − 2 1 2 0 0 0 0 0 ) = U B=\begin{pmatrix} 1&-1&-1&1&0\\ 1&-1&1&-3&1\\ 1&-1&-2&3&-\frac{1}{2} \end{pmatrix} \overset{r}{\sim} \begin{pmatrix} 1&-1&0&-1&\frac{1}{2}\\ 0&0&1&-2&\frac{1}{2}\\ 0&0&0&0&0 \end{pmatrix}=U B= 111−1−1−1−11−21−3301−21 ∼r 100−100010−1−2021210 =U -
可以看到行最简形矩阵 U U U不是强化行最简的
-
非自由未知数包含 r = 2 r=2 r=2个,自由未知数包含 n − r = 4 − 2 = 2 n-r=4-2=2 n−r=4−2=2个
-
比较推荐的做法是将第 1 , 3 1,3 1,3列(都是 n n n单位坐标列向量)对应的未知数作为非自由未知数,这样容易被其他自由未知数表达
- x 1 = x 2 + x 4 + 1 2 x_1=x_2+x_4+\frac{1}{2} x1=x2+x4+21
- x 3 = 2 x 4 + 1 2 x_3=2x_4+\frac{1}{2} x3=2x4+21
-
取 x 2 = x 4 = 0 x_2=x_4=0 x2=x4=0构造 ( 1 ) (1) (1)的特解
- ξ ∗ = ( 1 2 0 1 2 0 ) \xi^*=\begin{pmatrix} \frac{1}{2}\\ 0\\ \frac{1}{2}\\ 0 \end{pmatrix} ξ∗= 210210
-
导出组的基础解系:
-
x 1 = x 2 + x 4 x_1=x_2+x_4 x1=x2+x4
-
x 3 = 2 x 4 x_3=2x_4 x3=2x4
-
ξ 1 = ( 1 1 0 0 ) ξ 2 = ( 1 0 2 1 ) \xi_1=\begin{pmatrix} 1\\ 1\\ 0\\ 0 \end{pmatrix} \xi_2=\begin{pmatrix} 1\\ 0\\ 2\\ 1 \end{pmatrix} ξ1= 1100 ξ2= 1021
-
-
(2)的通解为 ξ ∗ + ∑ i = 1 2 k i ξ i \xi^{*}+\sum_{i=1}^{2}k_i\xi_i ξ∗+∑i=12kiξi, k i ∈ R k_i\in\mathbb{R} ki∈R
【推荐】国内首个AI IDE,深度理解中文开发场景,立即下载体验Trae
【推荐】编程新体验,更懂你的AI,立即体验豆包MarsCode编程助手
【推荐】抖音旗下AI助手豆包,你的智能百科全书,全免费不限次数
【推荐】轻量又高性能的 SSH 工具 IShell:AI 加持,快人一步
· 分享4款.NET开源、免费、实用的商城系统
· 全程不用写代码,我用AI程序员写了一个飞机大战
· MongoDB 8.0这个新功能碉堡了,比商业数据库还牛
· 记一次.NET内存居高不下排查解决与启示
· 白话解读 Dapr 1.15:你的「微服务管家」又秀新绝活了