用Matlab學(xué)習(xí)線性代數(shù)_線性方程組與矩陣代數(shù)

1、 用Matlab學(xué)習(xí)線性代數(shù)線性方程組與矩陣代數(shù)實(shí)驗(yàn)?zāi)康模菏煜ぞ€性方程組的解法和矩陣的基本運(yùn)算及性質(zhì)驗(yàn)證。Matlab命令：本練習(xí)中用到的Matlab命令有：inv，floor，rand，tic，toc，rref，abs，max，round，sum，eye，triu，ones，zeros。本練習(xí)引入的運(yùn)算有：+，-，*，。其中+和-表示通常標(biāo)量及矩陣的加法和減法運(yùn)算；*表示標(biāo)量或矩陣的乘法；對(duì)所有元素為實(shí)數(shù)的矩陣，運(yùn)算對(duì)應(yīng)于轉(zhuǎn)置運(yùn)算。若為一個(gè)非奇異矩陣(det!=0)且為一個(gè)矩陣，則運(yùn)算等價(jià)于。實(shí)驗(yàn)內(nèi)容：1 用Matlab隨機(jī)生成的矩陣和。求下列指定的，并確定那些矩陣是相等的。你可以利用Mat

2、lab計(jì)算兩個(gè)矩陣的差來(lái)測(cè)試兩個(gè)矩陣是否相等。（1）C=A*B,D=B*A,G=(A*B),H=(B*A) C=H;D=G；（2）C=A*B,D=(A*B),G=B*A,H=(B*A) C=H;D=G；（3）C=inv(A*B),D=inv(A)*inv(B),G=inv(B*A),H=inv(B)*inv(A)（4）C=inv(A*B),D=inv(A*B),G=inv(A)*inv(B),H=(inv(A)*inv(B)（3）（4）中無(wú)相等的2 令 n=200，并使用命令A(yù)=floor(10*rand(n); b=sum(A)z=ones(n,1); 注釋：（n行一列全為1的矩陣）生成一個(gè)

3、矩陣和兩個(gè)中的向量，它們的元素均為整數(shù)。（因?yàn)榫仃嚭拖蛄慷己艽?，我們添加分?hào)來(lái)控制輸出。（1） 方程組 的真解應(yīng)為。為什么？ 【A中的每一行的元素之和正好等于對(duì)應(yīng)b的每一列，故z為其一解，又det不等于0，RA=RAb=n，故z為其解】試說明，可在Matlab中利用”運(yùn)算或計(jì)算，然后用計(jì)算來(lái)求解。比較這兩種計(jì)算方法的速度和精度。我們將使用Matlab命令tic和toc來(lái)測(cè)量每一個(gè)計(jì)算過程消耗的時(shí)間。只需要用下面的命令：tic,x=Ab；toctic,y=inv(A)*b; toc 哪一種方法更快？ tic,x=Ab；更快！為了比較這兩種方法的精度，可以測(cè)量求得的解x和y與真解z接近的程度。利用

4、下面的命令：max(abs(x-z)max(abs(y-z)哪種方法的到的解更精確？>> max(abs(x-z)= 4.0168e-013 更精確！>> max(abs(y-z) = 6.1107e-013（2） 用n=500和n=1000替換（1）中的。 如（1）結(jié)果一樣！3 令A(yù)=floor(10*rand(6)。根據(jù)構(gòu)造，矩陣A將有整數(shù)元。將矩陣A的第六列更改，使得矩陣A為奇異的。令B=A，A(:,6)=-sum(B(1:5,:)（1） 設(shè)x=ones(6,1),并利用Matlab計(jì)算Ax。為什么我們知道A必為奇異的？【因化簡(jiǎn)列，>列成比例】試說明。通過化

5、為行最簡(jiǎn)形來(lái)判斷A是奇異的。（2） 令B=x*1:6，乘積AB應(yīng)為零矩陣。為什么？【因A的每一行的前五個(gè)元素之和等于第六個(gè)元素的相反數(shù)，且在A上的每一行的元素同乘以相同的數(shù)，則仍等于0】試說明。用Matlab的*運(yùn)算計(jì)算AB進(jìn)行驗(yàn)證。（3） 令C=floor(10*rand(6)和D=B+C，盡管,但乘積AC和AD是相等的。為什么？試說明。計(jì)算A*C和A*D，并驗(yàn)證它們確實(shí)相等?！敬颂嶣為令B=x*1:6；A為A(:,6)=-sum(B(1:5,:)】由于A*B=0；故AC=AD；A(B+C)=AB+AC;4 采用如下方式構(gòu)造一個(gè)矩陣。令B=eye(10)-triu(ones(10),1)，

6、參見最后附表二：為什么我們知道B必為非奇異的？【上三角矩陣的行列式的值等于對(duì)角線上的元素相乘】令C=inv(B)且x=C(:,10)，現(xiàn)在用B(10,1)=-1/256將B進(jìn)行微小改變。利用Matlab計(jì)算乘積Bx。由這個(gè)計(jì)算結(jié)果，你可以得出關(guān)于新矩陣B的什么結(jié)論？【化簡(jiǎn)此時(shí)B，得行最簡(jiǎn)式，RB=9<10，可以得出B的第10列（從19行）與x互為相反數(shù)，且都是2的指數(shù)冪數(shù)，且第十行為0，】 它是否為奇異的？【是】 試說明。用Matlab計(jì)算它的行最簡(jiǎn)形。5 生成一個(gè)矩陣A：A=floor(20*rand(6)并生成一個(gè)向量b：B=floor(20*rand(6,1)-10（1） 因?yàn)锳是

7、隨機(jī)生成的，我們可以認(rèn)為它是非奇異的。那么方程組應(yīng)有唯一解。用運(yùn)算“”求解。用Matlab計(jì)算A b的行最簡(jiǎn)形U。比較U的最后一列和解x，結(jié)果是什么？【相等】在精確算術(shù)運(yùn)算時(shí)，它們應(yīng)當(dāng)是相等的。為什么？【行最簡(jiǎn)式中可寫出對(duì)應(yīng)元素的實(shí)際含義，對(duì)應(yīng)處的未知元就等于最后的數(shù)】試說明。為比較他們兩個(gè)，計(jì)算差U(:,7)-x或用format long考慮它們。（2） 現(xiàn)在改變A，試它成為奇異的。令 A(:,3)=A(:,1:2)*4 3【第一列乘以4加上第二列乘以3替換到第三列上】，利用Matlab計(jì)算rref(A b)。方程組有多少組解？【無(wú)解】試說明。【RA<RAB】（3） 令 y=floor

8、(20*rand(6,1)-10 且 c=A*y，為什么我們知道方程組Ax=c必為 相容？的？【x此時(shí)必有一解y，故為相容的】試說明。計(jì)算A c的行最簡(jiǎn)形U。方程組有多少組解？【無(wú)窮多解】試說明?！綬A=RA c<6】（4） 由行最簡(jiǎn)形確定的自由變量應(yīng)為。通過考察矩陣U對(duì)應(yīng)的方程組，可以求得時(shí)所對(duì)應(yīng)的解。將這個(gè)解作為列向量輸入Matlab中。為檢驗(yàn)，計(jì)算剩余向量。（5） 令。矩陣U應(yīng)對(duì)應(yīng)于的行最簡(jiǎn)形。用求自變量時(shí)齊次線性方程組的解（手工計(jì)算），并將你的結(jié)果輸入為向量Z。用A*Z檢驗(yàn)?zāi)愕慕Y(jié)論。（6） 令。向量v應(yīng)為方程組的解。為什么？試說明。用Matlab計(jì)算剩余向量來(lái)驗(yàn)證v為方程組的解。

9、在這個(gè)解中，自由變量的取值是什么？ 【=3】 如何使用向量w和z來(lái)求所有可能的方程組的解？【v=w+n*z,其中n為任意實(shí)數(shù)】試說明。6 考慮下圖： （1） 確定圖的鄰接矩陣A，將其輸入Matlab；（2） 計(jì)算A2并確定長(zhǎng)度為2的路的條數(shù)【72】，其起止點(diǎn)分別為：【A2+A中的數(shù)值之和，數(shù)字表示有幾種路徑，具體看程序】（3） 計(jì)算A4、A6、A8并回答(2)中各種情況長(zhǎng)度為4、【368】6、【2362】8、【15800】的路的條數(shù)。試推測(cè)什么時(shí)候從頂點(diǎn)Vi到Vj沒有長(zhǎng)度為偶數(shù) 【即為0】 的路。 【i=1，j=6; i=2,j=5; i=3,j=6或8； i=4,j=7; i=5,j=8;i

10、=6,j=1或3; i=7,j=4; i=8,j=3或6;】（4） 計(jì)算A3、A5、A7并回答（2）中各情況長(zhǎng)度為3、【154】5、【922】7【6098】的路的條數(shù)。你由（3）得到的推測(cè)對(duì)長(zhǎng)度為奇數(shù)的路是否成立？【不成立】，試說明【見程序】。推測(cè)根據(jù)i+j+k的奇偶性，是否存在長(zhǎng)度為k的路。【若i+j+k為偶數(shù)，不存在；相反，則存在】 【路徑見程序】（5） 如果我們?cè)趫D中增加邊V3,V6，V5,V8，新圖的鄰接矩陣B可首先令B=A，然后令B(3,6)=144, B(6,3)=1, B(5,8)=1, B(8,5)=1，對(duì)k=2,3,4,5計(jì)算Bk。（4）中的推測(cè)在新的圖形中是否還是成立的？【

11、不成立】見程序】 （6） 在圖中增加V6,V8，并構(gòu)造得到的圖的鄰接矩陣C，計(jì)算C的冪次，并驗(yàn)證你在(4)中的推測(cè)對(duì)這個(gè)新圖是否仍然成立?！静怀闪ⅰ俊疽姵绦颉?令A(yù)=magic(8)，然后計(jì)算其行最簡(jiǎn)形。使得首1對(duì)應(yīng)于前三個(gè)變量，且其余的五個(gè)變量均為自由的。（1）令c=1:8，通過計(jì)算矩陣A c的行最簡(jiǎn)形確定方程組Ax=c是否相容。方程組是相容的嗎？ 【不相容】 試說明。 【RA<RAc】（2）令 b=8 -8 -8 8 8 -8 -8 8;并考慮方程組Ax=b。該方程組應(yīng)為相容的。通過U=rref(A b)驗(yàn)證。對(duì)五個(gè)自由變量的任一組取值，我們都應(yīng)可以得到一組解。事實(shí)上，令x2=flo

12、or(10*rand(5,1)，若x2表示方程組解的最后5個(gè)坐標(biāo)，則我們由x2求得x1=(x1,x2,x3)。要這樣做，只需要令U=rref(A b)。U的非零行對(duì)應(yīng)于分塊形式的線性方程組為解此方程組，令V=U(1:3,4:8),c=U(1:3,9)并利用Matlab，根據(jù)x2，c和V計(jì)算x1。令x=x1;x2，驗(yàn)證x是方程組的解。8令 B=-1,-1;1,1和A=zeros(2),eye(2);eye(2),B 驗(yàn)證B2=0。（1）用Matlab計(jì)算A2，A4，A6，A8。猜想用子矩陣E，O和B如何表示分塊形式的A2k。用數(shù)學(xué)歸納法證明你的猜想對(duì)任何正整數(shù)k都是成立的。（2）用Matlab計(jì)

13、算A3，A5，A7和A9。猜想用子矩陣E，O和B如何表示分塊形式的A2k-1。用數(shù)學(xué)歸納法證明你的猜想對(duì)任何正整數(shù)k都是成立的。9（1） Matlab命令A(yù)=floor(10*rand(6),B=A*A將得到元素為整數(shù)的對(duì)稱矩陣。為什么？試說明?！镜趇行第j列的數(shù)等于第i列的數(shù)分別乘以第j列的數(shù)之和；第j行第i列的數(shù)等于第j列的數(shù)分別乘以第i列的數(shù)之和，故為對(duì)稱矩陣】用這種方法計(jì)算B來(lái)驗(yàn)證結(jié)論，然后將B劃分成四個(gè)3x3的子矩陣。在Matlab中求子矩陣，令 B11=B(1:3,1:3),B12=B(1:3,4:6)并用B的第四行到第6行類似定義B21和B22。（2）令 C=inv(B11)。應(yīng)

14、有CT=C和B21T=B12。為什么？【對(duì)稱陣的逆矩陣與該逆矩陣的轉(zhuǎn)置是相等的，B12的第i行的數(shù)等于B21的第i列的數(shù)】試說明。用Matlab運(yùn)算符計(jì)算轉(zhuǎn)置，并驗(yàn)證結(jié)論。然后，令G=B21*C 和 H=B22-B21*C*B21利用Matlab函數(shù)eye和zeros構(gòu)造計(jì)算W=L*D*L，并通過計(jì)算W-B與B進(jìn)行比較。證明：若用算術(shù)運(yùn)算精確計(jì)算LDLT，它應(yīng)準(zhǔn)確等于B。附表：第一題：（1）>> A=rand(4);>> B=rand(4);>> C=A*B;>> D=B*A;>> G=(A'*B')'>

15、> H=(B'*A')'>> C-Dans = 2.2376e-001 4.7289e-001 1.3979e+000 1.3204e+000 -6.3633e-001 -3.0354e-001 2.2485e-002 -1.5056e-001 -1.7227e-001 -1.1938e-001 2.9484e-001 2.3624e-001 -8.7955e-001 -6.5016e-001 8.0370e-002 -2.1506e-001>> C-Gans = 2.2376e-001 4.7289e-001 1.3979e+000 1.

16、3204e+000 -6.3633e-001 -3.0354e-001 2.2485e-002 -1.5056e-001 -1.7227e-001 -1.1938e-001 2.9484e-001 2.3624e-001 -8.7955e-001 -6.5016e-001 8.0370e-002 -2.1506e-001>> C-Hans = 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0>> D-Gans = 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0>> D-Hans = -2.2376e-001 -4.7289e-0

17、01 -1.3979e+000 -1.3204e+000 6.3633e-001 3.0354e-001 -2.2485e-002 1.5056e-001 1.7227e-001 1.1938e-001 -2.9484e-001 -2.3624e-001 8.7955e-001 6.5016e-001 -8.0370e-002 2.1506e-001>> G-Hans = -2.2376e-001 -4.7289e-001 -1.3979e+000 -1.3204e+000 6.3633e-001 3.0354e-001 -2.2485e-002 1.5056e-001 1.722

18、7e-001 1.1938e-001 -2.9484e-001 -2.3624e-001 8.7955e-001 6.5016e-001 -8.0370e-002 2.1506e-001>>（2）>> C=A'*B'>> D=(A*B)'>> G=B'*A'>> H=(B*A)'>> C-Dans = -2.2376e-001 6.3633e-001 1.7227e-001 8.7955e-001 -4.7289e-001 3.0354e-001 1.1938e-001 6.

19、5016e-001 -1.3979e+000 -2.2485e-002 -2.9484e-001 -8.0370e-002 -1.3204e+000 1.5056e-001 -2.3624e-001 2.1506e-001>> C-Gans = -2.2376e-001 6.3633e-001 1.7227e-001 8.7955e-001 -4.7289e-001 3.0354e-001 1.1938e-001 6.5016e-001 -1.3979e+000 -2.2485e-002 -2.9484e-001 -8.0370e-002 -1.3204e+000 1.5056e-

20、001 -2.3624e-001 2.1506e-001>> C-Hans = 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0>> D-Gans = 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0>> D-Hans = 2.2376e-001 -6.3633e-001 -1.7227e-001 -8.7955e-001 4.7289e-001 -3.0354e-001 -1.1938e-001 -6.5016e-001 1.3979e+000 2.2485e-002 2.9484e-001 8.0370e-002 1.3204

21、e+000 -1.5056e-001 2.3624e-001 -2.1506e-001>> G-Hans = 2.2376e-001 -6.3633e-001 -1.7227e-001 -8.7955e-001 4.7289e-001 -3.0354e-001 -1.1938e-001 -6.5016e-001 1.3979e+000 2.2485e-002 2.9484e-001 8.0370e-002 1.3204e+000 -1.5056e-001 2.3624e-001 -2.1506e-001>>(3)>> C=inv(A*B);>>

22、D=inv(A)*inv(B);>> G=inv(B*A);>> H=inv(B)*inv(A);>> C-Dans = -3.9602e+001 -1.4016e+001 1.4537e+001 2.2261e+001 1.5266e+001 1.5778e+001 -1.9398e+001 -3.9304e+001 1.0821e+001 1.4313e+000 -2.7296e+001 -4.8956e+001 1.3845e+001 -5.5182e-001 2.6289e+001 5.1120e+001>> C-Gans = -3.96

23、02e+001 -1.4016e+001 1.4537e+001 2.2261e+001 1.5266e+001 1.5778e+001 -1.9398e+001 -3.9304e+001 1.0821e+001 1.4313e+000 -2.7296e+001 -4.8956e+001 1.3845e+001 -5.5182e-001 2.6289e+001 5.1120e+001>> C-Hans = -5.6843e-014 -1.2879e-014 3.0198e-014 7.1054e-014 -6.5370e-013 -1.4744e-013 3.3396e-013 8

24、.2423e-013 -1.5774e-012 -3.5527e-013 7.8870e-013 1.9895e-012 1.8758e-012 4.2988e-013 -9.4502e-013 -2.4016e-012>> D-Gans = 4.9738e-013 1.1013e-013 -8.3489e-014 -3.1264e-013 1.7053e-013 3.7303e-014 -2.4869e-014 -1.0747e-013 5.8265e-013 1.3145e-013 -9.4147e-014 -3.8369e-013 -1.0516e-012 -2.3448e-

25、013 1.7053e-013 6.6791e-013>> D-Hans = 3.9602e+001 1.4016e+001 -1.4537e+001 -2.2261e+001 -1.5266e+001 -1.5778e+001 1.9398e+001 3.9304e+001 -1.0821e+001 -1.4313e+000 2.7296e+001 4.8956e+001 -1.3845e+001 5.5182e-001 -2.6289e+001 -5.1120e+001>> G-Hans = 3.9602e+001 1.4016e+001 -1.4537e+001

26、-2.2261e+001 -1.5266e+001 -1.5778e+001 1.9398e+001 3.9304e+001 -1.0821e+001 -1.4313e+000 2.7296e+001 4.8956e+001 -1.3845e+001 5.5182e-001 -2.6289e+001 -5.1120e+001>>（4）>> c=inv(A*B)');>> d=inv(A'*B');>> g=inv(A')*inv(B');>> h=(inv(A)*inv(B)'>&

27、gt; c-dans = -3.9602e+001 1.5266e+001 1.0821e+001 1.3845e+001 -1.4016e+001 1.5778e+001 1.4313e+000 -5.5182e-001 1.4537e+001 -1.9398e+001 -2.7296e+001 2.6289e+001 2.2261e+001 -3.9304e+001 -4.8956e+001 5.1120e+001>> c-gans = -1.6875e-014 -5.4712e-013 -1.3216e-012 1.5774e-012 -2.8866e-015 -1.3145

28、e-013 -3.1264e-013 3.7659e-013 8.8818e-015 2.6290e-013 6.3949e-013 -7.6028e-013 2.5757e-014 7.1765e-013 1.7195e-012 -2.0606e-012>> c-hans = -3.9602e+001 1.5266e+001 1.0821e+001 1.3845e+001 -1.4016e+001 1.5778e+001 1.4313e+000 -5.5182e-001 1.4537e+001 -1.9398e+001 -2.7296e+001 2.6289e+001 2.226

29、1e+001 -3.9304e+001 -4.8956e+001 5.1120e+001>> d-gans = 3.9602e+001 -1.5266e+001 -1.0821e+001 -1.3845e+001 1.4016e+001 -1.5778e+001 -1.4313e+000 5.5182e-001 -1.4537e+001 1.9398e+001 2.7296e+001 -2.6289e+001 -2.2261e+001 3.9304e+001 4.8956e+001 -5.1120e+001>> d-hans = -2.4158e-013 -1.1724

30、e-013 -2.7711e-013 5.2580e-013 -5.6843e-014 -1.8652e-014 -5.3291e-014 1.0658e-013 4.2633e-014 1.7764e-014 4.7962e-014 -8.8818e-014 1.5987e-013 6.7502e-014 1.8474e-013 -3.3396e-013>> g-hans = -3.9602e+001 1.5266e+001 1.0821e+001 1.3845e+001 -1.4016e+001 1.5778e+001 1.4313e+000 -5.5182e-001 1.45

31、37e+001 -1.9398e+001 -2.7296e+001 2.6289e+001 2.2261e+001 -3.9304e+001 -4.8956e+001 5.1120e+001>>第二題：（1）>> n=200;>> A=floor(10*rand(n);>> b=sum(A')'>> z=ones(n,1);>> c=linsolve(A,b);>> d=c-z；精度為 1e-14-1e-13;tic,x=Ab,toc = Elapsed time is 0.016000 sec

32、onds.tic,inv(A)*b,toc = Elapsed time is 0.031000 seconds. （2）n=500;>> tic,x=Ab;tocElapsed time is 0.187000 seconds. 更快！>> tic,y=inv(A)*b;tocElapsed time is 0.343000 seconds.>> max(abs(x-z)=4.3987e-013 更精確！>> max(abs(y-z) = 2.2524e-012>>>> n=1000;>> tic,x=Ab;

33、tocElapsed time is 0.920000 seconds. 更快！>> tic,y=inv(A)*b;tocElapsed time is 1.404000 seconds.>> max(abs(x-z) =1.8221e-012 更精確！>> max(abs(y-z) =2.0862e-011>>（3）>> A=floor(10*rand(6);>>B=A;>> A(:,6)=-sum(B(1:5,:)'A = 0 6 7 7 0 -20 5 8 4 7 0 -24 6 7 4 3 3

34、-23 8 5 8 3 3 -27 1 8 3 9 4 -25 7 3 2 8 8 -28>> x=ones(6,1);>> b=A*xb = 0 0 0 0 0 0>> det(A)= 0>> rref(A)ans = 1 0 0 0 0 -1 0 1 0 0 0 -1 0 0 1 0 0 -1 0 0 0 1 0 -1 0 0 0 0 1 -1 0 0 0 0 0 0>>3.2>> A=floor(10*rand(6)；>> B=A'>> A(:,6)=-sum(B(1:5,:)'

35、;；>> B=x*1:6B = 1 2 3 4 5 6 1 2 3 4 5 6 1 2 3 4 5 6 1 2 3 4 5 6 1 2 3 4 5 6 1 2 3 4 5 6>> A*Bans = 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 03.3>> A=floor(10*rand(6);>> B=A'>> A(:,6)=-sum(B(1:5,:)'>> C=floor(10*rand(6);>>

36、 B=x*1:6;>> D=B+C;>> A*C-A*Dans = 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0>>4題：>> B=eye(10)-triu(ones(10),1);>> C=inv(B);>> B(10,1)=-1/256;>> B=eye(10)-triu(ones(10),1);>> B(10,1)=-1/256;行最簡(jiǎn)形：>> rref(B)ans = 1 0 0 0

37、 0 0 0 0 0 -256 0 1 0 0 0 0 0 0 0 -128 0 0 1 0 0 0 0 0 0 -64 0 0 0 1 0 0 0 0 0 -32 0 0 0 0 1 0 0 0 0 -16 0 0 0 0 0 1 0 0 0 -8 0 0 0 0 0 0 1 0 0 -4 0 0 0 0 0 0 0 1 0 -2 0 0 0 0 0 0 0 0 1 -1 0 0 0 0 0 0 0 0 0 0>> d=B*C; >>C = 1 1 2 4 8 16 32 64 128 256 0 1 1 2 4 8 16 32 64 128 0 0 1 1 2 4

38、8 16 32 64 0 0 0 1 1 2 4 8 16 32 0 0 0 0 1 1 2 4 8 16 0 0 0 0 0 1 1 2 4 8 0 0 0 0 0 0 1 1 2 4 0 0 0 0 0 0 0 1 1 2 0 0 0 0 0 0 0 0 1 1 0 0 0 0 0 0 0 0 0 1行最簡(jiǎn)形：>> rref(d)ans = 1 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0

39、0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0>> det(B)ans = 0>> rref(B c)ans = 1 0 0 0 0 0 0 0 0 -256 0 0 1 0 0 0 0 0 0 0 -128 0 0 0 1 0 0 0 0 0 0 -64 0 0 0 0 1 0 0 0 0 0 -32 0 0 0 0 0 1 0 0 0 0 -16 0 0 0 0 0 0 1 0 0 0 -8 0 0 0 0 0 0 0 1 0 0 -4 0 0 0 0 0

40、 0 0 0 1 0 -2 0 0 0 0 0 0 0 0 0 1 -1 0 0 0 0 0 0 0 0 0 0 0 1>>第五題：（1）>> A=floor(20*rand(6);>> B=floor(20*rand(6,1)-10;>> x=ABx = 0.32683 -2.7613 1.6764 0.76772 -0.33957 -1.5678>> C=A B;>> a=rref(C)a = 1 0 0 0 0 0 0.32683 0 1 0 0 0 0 -2.7613 0 0 1 0 0 0 1.6763 0 0

41、0 1 0 0 0.76772 0 0 0 0 1 0 -0.33957 0 0 0 0 0 1 -1.5678>> a(:,7)-xans = 1.3288e-007 9.5485e-006 -4.8037e-006 -2.7883e-006 5.0983e-007 -1.9237e-006>>（2）>> A=floor(20*rand(6);>> B=floor(20*rand(6,1)-10;>> A(:,3)=A(:,1:2)*4 3'>> rref(A B)ans = 1 0 4 0 0 0 0 0 1

42、3 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 1(3)>> y=floor(20*rand(6,1)-10;>> A=floor(20*rand(6);>> A(:,3)=A(:,1:2)*4 3'>> c=A*yc = 184 289 139 197 147 292>> rref(A c)ans = 1 0 4 0 0 0 1 0 1 3 0 0 0 6 0 0 0 1 0 0 8 0 0 0 0 1 0 -1 0 0 0 0 0 1 6 0 0 0 0 0 0 0>>(4)>> A=floor(20*rand(6);>> A(:,3)=A(:,1:2)*4 3'>> y=floor(20*rand(6,1)-10;>> c=A*y;>> U=rref(A c)U = 1 0 4 0 0 0 23 0 1 3 0 0 0 17 0 0 0 1 0 0 9 0 0 0