非齐次线性方程组
- 格式:doc
- 大小:459.50 KB
- 文档页数:6
非齐次线性方程组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. 求解线性方程组1234524512345123512345343226333434222026231x 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 1rref(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化为方程组32415510x x x x x x ++=-⎧⎪=⎨⎪=-⎩所以解为15233354555311000001100011010x x x x x x x x x x x x --⎛⎫⎛⎫⎛⎫⎛⎫⎛⎫ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪==++ ⎪ ⎪ ⎪ ⎪ ⎪-- ⎪ ⎪ ⎪ ⎪ ⎪⎪ ⎪ ⎪ ⎪ ⎪⎝⎭⎝⎭⎝⎭⎝⎭⎝⎭++例2. 设函数2y axbx c =++经过点(1,1), (2,2), (3,0), 求系数a , b , c .解1422930a b c a b c a b c ++=⎧⎪++=⎨⎪++=⎩输入系数矩阵A 和右端项bA=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,,322a b c =-==-解二判断方程组是否有解, 即系数矩阵A 的秩是否等于增广矩阵1A 的秩. rank(A)==rank(A1)ans = 1 有解.判断方程组是否有唯一解, 即系数矩阵 A 是否等于A 的列数n .[m,n]=size(A); rank(A)==nans = 1A 的秩等于列数n , 有唯一解.直接用A 左除 b 求解 x=A\bx = -3/2 11/2 -3例 3. 设三种食物中每100g 中的蛋白质、碳水化合物、脂肪的含量如下表.三种食物用量各为多少才能保证所需营养?解. 设脱脂牛奶用量为1x , 大豆面粉用量为2x , 乳清用量为3x .12312312336 51 133352 34 74450 7 1.13x 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 。
非齐次线性方程组解法
非齐次线性方程组的通解=齐次线性方程组的通解+非齐次线性
方程组的一个特解(η=ζ+η*)。
非齐次线性方程组是常数项不全为零的线性方程组。
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 是互不相同的常数,证明下面的方程组无解。
非齐次线性方程组解的结构的进一步讨论摘要:本文通过矩阵的初等变换及非齐次线性方程组的解的有关性质进一步讨论了非齐次线性方程组的解的结构问题,虽然非齐次线性方程组的解向量的全体不能构成向量空间,也没有基础解系,但我们找到了类似齐次线性方程组的基础解系的解向量组,这个解向量组线性无关。
并且的任意一个解都可以由这个解向量组线性表示。
最后,给出了非齐次线性方程组有全非零解的充要条件,并给出了相应例题。
关键字:非零解,基础解系,线性无关,初等变换引言非其次线性方程组⎪⎪⎩⎪⎪⎨⎧=+++=+++=+++nn mn m m n n n n b x a x a x a b x a x a x a b x a x a x a 22112222222********* (Ⅰ)的矩阵形式为B AX =.取0=B ,得到其次线性方程组0=AX 称为非其次线性方程组B AX =的导出组。
我们知道非其次线性方程组B AX =的解有以下的一些性质:(1) 若1u 是非其次线性方程组B AX =的一个解,1v 是其导出组0=AX 的一个解,则11v u +也是0=AX 的一个解。
证明:因为1u 是非其次线性方程组B AX =的一个解,所以有B Au =1,同理有01=Av ,则由()B B Av Au v u A =+=+=+01111.所以11v u +是非其次线性方程组B AX =的解。
(2) 若21,v v 是非其次线性方程组的两个解,则21v v -是其导出组的解证明:由B Av =1,B Av =2,所以有()02121=-=-=-B B Av Av v v A ,故21v v -为其导出组的解。
2.定理(非其次线性方程组解的结构定理)若1v 是非其次线性方程组B AX =的一个解,v 是其导出组的通解,则11v v u +=是非其次线性方程组的通解。
证明:由性质(1)可知1u 加上其导出组的一个解仍是非其次线性方程组的一个解,所以只需证明,非其次线性方程组的任意一个解*v ,一定是1u 与其导出组某一个解1v 的和,取1*1u v v -=由性质(2)可知,1v 是导出组的一个解,于是得到11*v u v +=,即非其次线性方程组的任意一个解与其导出组的某一个解的和。
由上面这个定理我们可以知道,一个其次线性方程组的解的全体可以用基础解系来表示。
因此,根据定理我们可以用导出组的基础解系来表示出一般方程组的一般解,如果0v 是方程组(Ⅰ)的一个特解,r n -ηηη,,21 是其导出组的一个基础解系,那么(Ⅰ)的任一个解都可以表示成:r n r n k k k u u --++++=ηηη 22110*3.由上面2的证明过程,我们可以知道其次线性方程组0=AX 的全部解可由基础解系r n -ηηη,,21 线性表示出(其基础解系含有r n -个解向量),即()r n r n r n k k k k k k ---+++ ,,212211ηηη为任意实数。
那么,当非其次线性方程组B AX =有解时,则B AX =至多有多少个线性无关的解向量?B AX =的全部解又如何表示?定理若其次线性方程组0=AX 的基础解系为r n -ηηη,,21 ,当非其次线性方程组0≠=B AX 有解时,则它至多且一定有1+-r n 个线性无关的解向量121,,,+--r n r n ξξξξ ,B AX =的通解可以表示为112211+-+---++++r n r n r n r n k k k k ξξξξ 为满足关系式1121=+++++--r n r n k k k k ,的任意实数。
证明:(ⅰ)若ξ是非其次线性方程组B AX =的解,则ξ为非零解向量,那么向量组ξ,r n -ηηη,,21 线性无关(否则ξ可由r n -ηηη,,21 线性表示,与ξ是B AX =的解矛盾)。
那么,易证r n r n -+-+=+=+==ηξξηξξηξξξξ123121,,, 都是B AX =的解,并且121,,,+--r n r n ξξξξ 线性无关。
这说明B AX =至少有1+-r n 个线性无关的解向量。
下面再证B AX =至多有1+-r n 个线性无关的解向量。
反证:若B AX =有2+-r n 个线性无关的解向量2121,,,+-+-r n r n ξξξξ ,那么易证212221,,+-+-+-+----r n r n r n r n ξξξξξξ 均为0=AX 的解,并且线性无关。
这样0=AX 具有1+-r n 线性无关的解向量矛盾,所以,B AX =至多且一定有1+-r n 个线性无关的解向量B AX =。
(ⅱ)对于B AX =的任意一个解,一定可以表示成它的一个特解ξ与其导出组0=AX 的基础解系的线性组合,即()r n r n r n k k k k k k ---++++ ,,212211ηηηξ为任意常数 那么()()()()()13221121221121221111+--------++++----=+++++++----=++++r n r n r n r n r n r n rn r n k k k k k k k k k k k k k k k ηξξξηξηξηξαξηηηξ(r n k k k - ,,21为任意实数,且组合系数()r n r n k k k k k k ------ ,,,12121之和等于1.这说明,B AX =的任意解都可以表示成这样的形式。
另一方面,由于121,,,+--r n r n ξξξξ 都是B AX =的解,对于112211+-+---++++r n r n r n r n k k k k ξξξξ ,只要满足1121=+++++--r n r n k k k k 仍然是B AX =的解,所以,B AX =的通解可以表示成112211+-+---++++r n r n r n r n k k k k ξξξξ ,且121,,,+--r n r n k k k k 为满足关系式1121=+++++--r n r n k k k k ,的任意实数。
例2设0η是线性方程组的一个解,t ηηηη ,,,321是它导出组的一个基础解系,令010201,,ηηηη+===+t t r r r 。
证明:线性方程组的任一一个解112211+++++=t t r u r u r u ν,其中1121=++++t u u u 。
证明:由题可设方程组的任一解ν可以表示成t t u u ηηην1120++++= (12,+t u u 为常数) 令1211+---=t u u u ,则()0101201112011)()(ηηηηηηηην+++++=++++=+++t t tt t u u u u u u u 112211+++++=t t r u r u r u(1) 引理:设A 为m n ⨯矩阵,用初等行变换,把A 化为阶梯形矩阵,并使该梯形矩阵的每一个非零行的第一个非零元素(从左算起)为1,且该元素所在列的其他元素为零,这样的阶梯形矩阵的为A 的行简化阶梯形矩阵。
定理:非齐次线性方程组存在全非零解的充要条件是,它的增广矩阵A 的秩()r A 与系数矩阵A 的秩()r A 相等,且A 的行简化阶梯型矩阵中每个非零行的非零元素个数大于或等于2.证明:必要性方程组有全非零解,则必须满足方程组的条件,因而,()()r A r A =.不妨设其秩为r 且A 的简化阶梯矩阵为:121,11,21,12,12,12,2,1,2100001000001000000000000n r r n r r n r r r r r r l l l d l l l d l l l d ++++++⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦(2) 且其对应的方程组为11,111,221122,112,2222,11,22r r r r n n r r r r n n r r r r r r r rn n rx l x l x l x d x l x l x l x d x l x l x l x d ++++++++++++=----+=----+=----+若对某个r ()i i r ≤≤ 有,1,2,0i r i r i n i l l l d ++=====则0i x =,这和方程组(2)有全非零全部解矛盾,故对每个i (1,2,r n =),至少存在一个j (1r j n +≤≤)使0ij l ≠或0i d ≠,即(2)中第i (1,2,r r =)行至少有两个非零元素。
充分性:设N 是充分大的正数,令1n r r x N -+=,1n rr x N -+=,n x N =将其带入(2)得:1.1,2,n rn r i i r i r i N i x l Nl N l d --+++=----+(1,2,r r =),当0ij l =(1r j n +≤≤),0i d ≠时,显然成立;当上式右端至少存在一个非零系数,设第一个非零系数为(),1i r k l k n r +≤≤-,则1,,1,1,1,1,1n r k n r k i i r k i r k in in r k i r k i n i n r k i r k n r k i r k x l N l N l N d l N l N d l Nl N --+--+++--++--++--++=----+⎡⎤---+=-+⎢⎥-⎢⎥⎣⎦因为,1,1,lim 0i r k n k r i n in r k N i r k l N l N d l N++----+→+∞+---+=-所以lim i N x →+∞=∞,1,2,i r =,故存在充分大的正数i N ,使1,,10n r k n r k i i r k i r k in i x l N l N l N d --+--+++=----+≠(1,2,i r =);取{}12max ,,r N N N N =,可使1,,10n r k n r k i i r k i r k in i x l N l N l N d --+--+++=----+≠(1,2,i r =)这样,就得到方程组的一个全非零解()112,,,,,,Tn r n r n r r x x x x N N N N ----=例1 方程组()()134512345123451245123452412626(1)03413251035x x x x x x x x x x x x x c x x x x c x x x x x c x ⎧+++=⎪-++-=⎪⎪-++--=⎨⎪+++-=⎪⎪-+++-=⎩ 有全非零解的充要条件? 解:其增广矩阵A 的简化阶梯形矩阵为102411010221001032000000000000c ⎡⎤⎢⎥--⎢⎥⎢⎥-⎢⎥⎢⎥⎢⎥⎣⎦故由上述定理可知,该方程组有全非零解的充要条件是c 为任意实数。