五矩陣運算與方程組求解

1、工程五 矩陣運算與方程組求解實驗3線性方程組實驗目的熟悉求解線性方程組的常用命令,能利用Mathematica命令各類求線性方程組的解.理解計算機求解的實用意義.根本命令1. 命令NullSpace A 1給出齊次方程組 AX =0的解空間的一個基.2. 命令LinearSolve A, b ,給出非齊次線性方程組AX =b的一個特解.3. 解一般方程或方程組的命令Solve見Mathematica入門.實驗舉例求齊次線性方程組的解空間設A為m n矩陣,X為n維列向量,那么齊次線性方程組 AX =0必定有解.假設矩陣A的 秩等于n,那么只有零解;假設矩陣A的秩小于n,那么有非零解,且所有解構(gòu)成

2、一向量空間.命令NullSpace給出齊次線性方程組 AX二0的解空間的一個基f X1 X2 -2X3 _X4 0,例3.1 教材 例3.1求解線性方程組1 3X1 -X2 -X3 “2X4 =0,5x2 - 7x3 3x4 =0,2X1 3x2 _5x3 - X4 0.輸入ClearA;A=1,1, -2,-1,3, 21,2,0,5,7,3,2, -3,-5,-1;NullSpaceA那么輸出 -2,1,-2,3說明該齊次線性方程組的解空間是一維向量空間,且向量,1,3是解空間的基注:如果輸出為空集 ,那么說明解空間的基是一個空集,該方程組只有零解.X! x2 ' 2x3 x4 =

3、0 例3.2求解線性方程組3X1 _2X2 _3X3 2X4 =0| 5X2 ' 7X3 3X4 =0J 2X1 -3X2 -5X3 -X4 =0輸入ClearA;A=1,1,2,-1,3,-2,-3,2,0,5,7,3,2,-3,-5,-1;NullspaceA輸岀為 因此解空間的基是一個空集,說明該線性方程組只有零解.例 3.3 (教材 例 3.2)向量組=(1,1,2,3), -.2 =(1,_1,1,1), ：.3 =(1,3,4,5), ：.4 =(3,1,5,7)是否線性相關(guān)？根據(jù)定義,如果向量組線性相關(guān),那么齊次線性方程組X1：1 *2用2=°有非零解.輸入Cl

4、earA,B;A=1,1,2,3,1,4,1,1,1,3,4,5,3,1,5,7;B=TransposeA;NullSpaceB輸岀為 2-1,0,1說明向量組線性相關(guān),且_：.2亠展4 =0非齊次線性方程組的特解捲 x2 2x3 % =4例3.4 教材 例3.3求線性方程組% 2X2 X3 +2X4 =2的特解.5x2 亠7x3 亠3% =-22X1 _3x2 -5x3 -X4 =4輸入ClearA,b;A=1,1, -2,-1,3, -2,-1,2,0,5,7,3,2,-3,-5,-1;b=4,2, -2,4LinearSolveA,b輸岀為1,1, -1,0注：命令LinearSolve

5、只給出線性方程組的一個特解| xi -X2 _2冷X44例3.5求線性方程組 嚴-2X2 -x3 +2滄=2的特解.5x2 十7x3 +3x4 =22x, 3x2 5x3 x4 =4輸入ClearA,b;A=1,1,2,-1,3,-2,-1,2,0,5,7,3,2,-3,-5,-1;b=4,2,2,4LinearsolveA,b輸岀為Linearsolve:nosol:Linear equation encountered which has no solution.說明該方程組無解.例3.6向量=(2, -1,3,4)是否可以由向量冷=(1,2,3,1), ：-2 =(5,5,12,11),

6、:-3 二 1, -3,6,3線性表示？根據(jù)定義,如果向量B可以由向量組 8,0(2,.3線性相關(guān),那么非齊次線性方程組x1 冷X2 ： 2 X3： 3 二：有解.輸入ClearA,B,b;A=1,2,-3,1,5,-5,12,11,0,5,7,3,1,-3,6,3;B=TransposeA;b=2,-1,3,4;LinearsolveB,b輸岀為1 1 1,1,03 311說明p可以由 耳笑,碼 線性表示,且卩=-8 +- a233例3.7 (教材 例3.4)求出通過平面上三點(0,7),(1,6)和(2,9)的二次多項式ax2 bx c,并畫岀其圖形.0 日 +0 b +c =7根據(jù)題設條

7、件有*1 a +1 b +c =6 ,輸入4 a 2 b c =9Clearx;A=0,0,1,1,1,1,4,2,1y=7,6,9p=LinearSolveA,yCleara,b,c,r,s,t;a,b,c.r,s,tfx_=p.xA2,x,1;Plotfx,x,0,2,GridLines>Automatic,PlotRange ->AII;那么輸出a,b,c的值為2, -3,7并畫出二次多項式2x2 _3x 7的圖形(略).非齊次線性方程組的通解用命令Solve求非齊次線性方程組的通解例3.8求出通過平面上三點(0,0),(1,1),(-1,3)以及滿足f ( V) =20,

8、f(1) =9的4次多項式 f (x).解 設 f (x) =ax4 - bx3 cx2 dx e,那么有|e =0a b cd e =1a -b - c -d e =3-4a 3b -2c d =204a亠3b亠2c亠d =9輸入Cleara,b,c,d,e;qx_=a*xA4+b*xA3+c*xA2+d*x+e; eqs=qO= =0,q1= =1,q-1= =3,q'1= =20,q '1= =9;A,y=LinearEquationsToMatriceseqs,a,b,c,d; p=LinearSolveA,y;fx_=p.xA4,xA3,xA2,x,1; Plotfx

9、,x,-1,1,GridLines->Automatic,PlotRange->AII;那么輸岀所求多項式f (x)19 4x435x,4非齊次線性方程組的通解用命令solve求非齊次線性方程組的通解fXix2 亠2x3 亠x4 =1例3.9解方程組2xi X2 X3 - 2x4 =3捲x3 +x4 =23xi X2 +3x4 =5輸入solvex-y+2z+w=1,2x-y+z+2w=3,x-z+w=2,3x-y+3w=5,x,y,z,w輸岀為x ； 2-w+z,yr 1+3z即X1 =2 -X4x3 , x2 =1 3x3 .于是,非齊次線性方程組的特解為2,1,0,0.對應的

10、齊次線性方程組的根底解系為1,3,1,0與-1,0,0,1.又一2x2 +3x3 4x4 =4例3.10解方程組1 X2 - X3 x 4 - -3X! +3x2 +% =1解法1用命令solve輸入solvex-2y+3z-4w=4, y-z+w=-3,x+3y+w=1,-7y+3z+3w=-3,x,y,z,w輸岀為x r -8,y r 3, z 6, w 0 即有唯一解 x1 = -8 , x23 , x3 =6 , x4 =0 .解法2這個線性方程組中方程的個數(shù)等于未知數(shù)的個數(shù),而且有唯一解,此解可以表示為x二Ab.其中A是線性方程組的系數(shù)矩陣,而b是右邊常數(shù)向量.于是,可以用逆陣計算

11、唯一解.輸入ClearA,b,x;A=1,-2,3,-4,0,1,-1,1,1,3,0,1,0,-7,3,1;b=4,-3,1,-3;x=lnverseA.b輸岀為-8,3,6,0解法3還可以用克拉默法計算這個線性方程組的唯一解.為計算各行列式,輸入未知數(shù)的系數(shù)向量,即系數(shù)矩陣的列向量.輸入Cleara,b,c,d,e;a=1,0,1,0;b=-2,1,3,-7;c=3,-1,0,3;d=-4,1,1,1;e=4,-3,1,-3;Dete,b,c,d/ Deta,b,c,dDeta,e,c,d/ Deta,b,c,dDeta,b,e,d/ Deta,b,c,dDeta,b,c,e/ Deta,

12、b,c,d輸岀為-8360bxt +x2 +x3 =1例3.10 教材例3.5當a為何值時,方程組峽+ax2 +x3 =1無解、有唯一解、有無窮Jxix2 ax3 =1多解？當方程組有解時,求通解.先計算系數(shù)行列式,并求a,使行列式等于0.輸入Cleara;Deta,1,1,1,a,1,1,1,a;Solve% =0,a那么輸出a ； 2,a ；1,a ；1當a =-2,a -1時,方程組有唯一解 輸入Solvea*x y z=1,x,a*y z =1,x'y,a*z =1,x,y,z那么輸出x- , y- , z- 2 +a2 +a2 +a當a亠2時,輸入Solve -2x+y+z=

13、1,x -2y+z=1,x+y -2z=1,x,y,z那么輸出 說明方程組無解.當a =1時,輸入Solvex+y+z=1,x+y+z=1,x+y+z=1,x,y,z那么輸出x ； 1 -y -z說明有無窮多個解.非齊次線性方程組的特解為1,0,0,對應的齊次線性方程組的根底解系為為二,1,0與,0,1.2X1 +X2 X3 +X4 =1例3.11 教材 例3.6求非齊次線性方程組3為_2x2 +x32x4 =4的通解Jx1 ._4x2 _3x3= -2解法1輸入A=2,1, -1,1,3, 21,用,1,4,5;b=1,4, /;particular=LinearSolveA,bnullsp

14、acebasis=NullSpaceA generalsolution=t*nullspacebasis1+k*nullspacebasis2+Flattenparticular generalsolution/MatrixForm解法2輸入B=2,1, -1,1,1,3, -2,1,書,4,1,4, £,5,-2RowReduceB/MatrixForm根據(jù)增廣矩陣的行最簡形,易知方程組有無窮多解.其通解為X2X31/71/76/7+k5/7-9/7-5/7k,t為任意常數(shù)實驗習題2x1 -X2 +3X3 =0,1. 解方程組2Xt +x2 +x3 =0,4X1 亠X2 亠2X3 =0.2xt 4x2 +5x3 +3x4 =0,2. 解方程組丿 3X1 -6X2 +4X3 +2X4 =0,4x1 一8x2 ' 17x3 11x4 =0.fX1 2x2 +3X3 4X4 =4,3. 解方程組"X2卞3 *X4 =£,X1 x3 -2x4 - -2Xi +2X2 +X3 X4 =2,4.