非齐次线性方程组
- 格式:doc
- 大小:94.00 KB
- 文档页数:3
非齐次线性方程组解法
非齐次线性方程组的通解=齐次线性方程组的通解+非齐次线性
方程组的一个特解(η=ζ+η*)。
非齐次线性方程组是常数项不全为零的线性方程组。
1、非齐次线性方程组解法
非齐次线性方程组Ax=b的求解步骤:
(1)对增广矩阵B施行初等行变换化为行阶梯形。
若R(A)<R(B),则方程组无解。
(2)若R(A)=R(B),则进一步将B化为行最简形。
(3)设R(A)=R(B)=r;把行最简形中r个非零行的非0首元所对应的未知数用其余n-r个未知数(自由未知数)表示,并令自由未知数分别等于C1,C2……,Cn-r,即可写出含n-r个参数的通解。
2、非齐次线性方程组解的判别
如果系数矩阵的秩小于增广矩阵的秩,方程组无解;如果系数矩阵的秩等于增广矩阵的秩,方程组有解。
在有解的情况下,如果系数矩阵的秩等于未知数的个数,非齐次线性方程组有唯一解。
如果系数矩阵的秩小于未知数的个数,非齐次线性方程组有无穷多解,如果有无穷多解,先求所对应齐次线性方程组的基础解系,再求出非齐次线性方程组的一个特解。
由此可知:如果非齐次线性方程组有无穷多解,则其对应的齐次线性方程组一定有非零解,且非齐次线性方程组的全部解(通解)可表示为:对应齐次线性方程组的通解+非齐次线性方程组的特解。
3.3 非齐次线性方程组3.3.1问λ 取何值时方程组1212122(4)70(2)2302560x x x x x x λλλ+-+=⎧⎪-++=⎨⎪++-=⎩有唯一解、无穷多个解、无解?并在有无穷多个解时求出其通解。
解:由于系数矩阵不是方阵,故只能使用初等行变换法。
22472562230112565686022A λλλλλλλλλ⎡⎤⎢⎥---⎡⎤⎢⎥⎢⎥=------⎢⎥⎢⎥⎢⎥⎢⎥---+⎣⎦⎢⎥⎣⎦① 当1λ=-时,2571115022000A -⎡⎤⎢⎥⎢⎥-⎢⎥⎢⎥⎣⎦,由()()2r A r A ==,知方程组有唯一解。
由 11011150111000A ⎡⎤-⎢⎥⎢⎥⎢⎥-⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦ 知唯一截为12111511x x ⎡⎤⎡⎤=-⎢⎥⎢⎥⎣⎦⎣⎦② 当1λ≠-时,256011(1)(12)002A λλλ⎡⎤⎢⎥-⎢⎥⎢⎥⎢⎥--⎢⎥⎣⎦,则若1λ=,则由()()2r A r A ==知有唯一解1251x x -⎡⎤⎡⎤=⎢⎥⎢⎥⎣⎦⎣⎦;若12λ=,则由()()2r A r A ==知也有唯一解121;21x x ⎡⎤⎡⎤⎢⎥=⎢⎥⎢⎥⎣⎦⎣⎦若1λ≠且12λ≠,则由()23()r A r A =≠=知方程组无解。
3.3.2 选择题(1)设A =1100011000111001⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦,1234a a b a a ⎡⎤⎢⎥⎢⎥=⎢⎥⎢⎥⎣⎦,Ax b =有解的充分必要条件为( D )。
(A )1234a a a a === (B )12341a a a a ==== (C )12340a a a a +++= (D )12340a a a a -+-=(2)非齐次线性方程组Ax b =,对应的导出组方程组0Ax =,则( D )正确。
(A ) 若0Ax =仅有零解,则Ax b =有唯一解 (B ) 若0Ax =有非零解,则Ax b =有无穷多组解 (C ) 若Ax b =有无穷多组解,则0Ax =仅有零解 (D ) 若Ax b =有无穷多组解,则0Ax =有非零解3.3.2设123,,a a a 是互不相同的常数,证明下面的方程组无解。
非齐次线性方程组Ax=b
一、基本理论
线性方程组Ax=b 有解条件: 系数矩阵A 的秩 = 增广矩阵(A,b )的秩.
非齐次线性方程组的解集结构:
若x 1是Ax=b 的一个特解, N (A )表示齐次线性方程组Ax=0的解空间, 则非齐次线性方程组Ax=b 的解集为x 1+N (A ).
解非齐次线性方程组的方法:
通过初等行变换将增广矩阵(A,b )化为最简行阶梯矩阵(A 1,b 1), 写出对应的方程组,根据方程组写出解.
二、Matlab 实现
调用rref(A )将A 化为最简行阶梯矩阵, 根据对应的方程组写出解.
若方程组有解, 且rank(A )=n ,即A 列满秩时, 方程组有唯一解. 此时可直接用A 左除b 求得唯一解:x=A\b .
三、例子
例1. 求解线性方程组
123452451234512351
2
3
4
5
343226333
434222026231
x x x x x x x x x x x x x x x x x x x x x x -++-=⎧⎪---=-⎪⎪-++-=⎨⎪++-=⎪-+-++=
⎪⎩
A=[3 -4 3 2 -1; 0 -6 0 -3 -3; 4 -3 4 2 -2; 1 1 1 0 -1; -2 6 -2 1 3]; b=[2; -3; 2; 0; 1]; A1=[A b]
A1 =
3 -
4 3 2 -1 2 0 -6 0 -3 -3 -3 4 -3 4 2 -2 2 1 1 1 0 -1 0 -2 6 -2 1 3 1
rref(A1)
ans =
1 0 1 0 -1 0 0 1 0 0 0 0 0 0 0 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0
化为方程组
324
1551
0x x x x x x ++=-⎧⎪=⎨⎪=
-⎩
所以解为
15233354555311000001100011010x x x x x x x x x x x x --⎛⎫⎛⎫⎛⎫⎛⎫⎛⎫ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪==++ ⎪ ⎪ ⎪ ⎪ ⎪-- ⎪ ⎪ ⎪ ⎪ ⎪
⎪ ⎪ ⎪ ⎪ ⎪
⎝⎭⎝⎭⎝⎭⎝⎭⎝⎭
++
例2. 设函数2
y ax
bx c =++经过点(1,1), (2,2), (3,0), 求系数a , b , c .
解
1422930a b c a b c a b c ++=⎧⎪
++=⎨⎪++=⎩
输入系数矩阵A 和右端项b
A=sym([1 1 1; 4 2 1; 9 3 1]); b=sym([1; 2; 0]);
增广矩阵1A A1=[A b]
A1 =
[ 1, 1, 1, 1] [ 4, 2, 1, 2] [ 9, 3, 1, 0]
利用rref 求解 R=rref(A1)
R =
[ 1, 0, 0, -3/2] [ 0, 1, 0, 11/2] [ 0, 0, 1, -3]
即解为
311
,,322
a b c =-==-
解二
判断方程组是否有解, 即系数矩阵A 的秩是否等于增广矩阵1A 的秩. rank(A)==rank(A1)
ans = 1 有解.
判断方程组是否有唯一解, 即系数矩阵 A 是否等于A 的列数n .
[m,n]=size(A); rank(A)==n
ans = 1
A 的秩等于列数n , 有唯一解.
直接用A 左除 b 求解 x=A\b
x = -3/2 11/2 -3
例 3. 设三种食物中每100g 中的蛋白质、碳水化合物、脂肪的含量如下表.
三种食物用量各为多少才能保证所需营养?
解. 设脱脂牛奶用量为1x , 大豆面粉用量为2x , 乳清用量为3x .
1231231
2
3
36 51 133352 34 74450 7 1.13
x x x x x x x x x ++=++=++=⎧⎪⎨⎪⎩
A=[36 51 13 33; 52 34 74 45; 0 7 1.1 3]
A =
36.0000 51.0000 13.0000 33.0000 52.0000 34.0000 74.0000 45.0000 0 7.0000 1.1000 3.0000 R=rref(A)
R =
1.0000 0 0 0.2772 0 1.0000 0 0.3919 0 0 1.0000 0.2332
所以脱脂牛奶的用量为27.72g ,大豆面粉的用量为39.19g ,乳清的用量为23.32g 。