線性代數(shù)應用舉例16829

1、線性代數(shù)應用實例取自線性代數(shù)機算與應用指導（matlab）版2011.12例例 1 平板穩(wěn)態(tài)溫度的計算平板穩(wěn)態(tài)溫度的計算為了計算平板形導熱體的溫度分布，將平板劃分為許多方格，每一個節(jié)點上的穩(wěn)態(tài)溫度將等于其周圍四個節(jié)點溫度的平均值。由此可得出階數(shù)與節(jié)點數(shù)相同的線性方程組，方程的解將取決于平板的邊界條件。 這個方法可以用來計算飛行器的蒙皮溫度等。t1t2t3t44/50304/50104/30204/2010324413412321tttttttttttt8060503041101401104101144321tttt平板溫度計算的模型整理為a=4,-1,-1,0; -1,4,0,-1; -1,0

2、,4,-1; 0,-1,-1,4;b=30; 50; 60; 80;u=rref(a,b)matlab 程序運行結(jié)果為：u = 1.0000 0 0 0 21.2500 0 1.0000 0 0 26.2500 0 0 1.0000 0 28.7500 0 0 0 1.0000 33.7500向高階系統(tǒng)擴展則要解 25 階的線性方程組。運行書上的程序得溫度分布如下 將平板分割得愈細，求出的解就愈精確。如果把上述區(qū)域分成 25 個點如右12233414=100 =72 =37 =9xxxxxxxx例例 2 交通流的建模交通流的建模對于一個有雙向車流的十字路口，根據(jù)流出流入車數(shù)相等的規(guī)則，可以列出

3、下列方程組:節(jié)點a：x1360 x2260節(jié)點b：x2220 x3292節(jié)點c：x3320 x4357節(jié)點d：x4260 x1251相應的矩陣方程為：a=1,-1,0,0;0,1,-1,0;0,0,1,-1;-1,0,0,1;b=-100;72;37;-9;u=rref(a,b)142434910937xxxxxxmatlab 程序運行結(jié)果為：u = 1 0 0 -1 9 0 1 0 -1 109 0 0 1 -1 37 0 0 0 0 0 由于 u 的最后一行為全零，也就是說，四個方程中實際上只有三個獨立。x4 可以任設，因為如果有一些車沿此路口環(huán)行，對方程無影響，故方程組的解可如上表示.

4、把上述模型擴展到多個十字路，乃至整個城市，就構成高階的線性代數(shù)方程組。例如下面的 6 節(jié)點交通流圖，它就要由 6 個方程和 7 個變量來描述。用行最簡型方法可以知道，它的解將包括兩個自由變量。其物理意義類推。向高階系統(tǒng)擴展左圖描述了四個城市之間的航空航線圖，其中1、2、3、4 表示四個城市；帶箭頭線段表示兩個城市之間的航線。設行號表示起點城市，列號為到達城市，則定義鄰接矩陣 a 為：0111001100001100a例例 3 飛機航線問題飛機航線問題轉(zhuǎn)機航線的數(shù)學模型不難證明：矩陣 a2=a*a 表示一個人連續(xù)坐兩次航班可以到達的城市，矩陣 a3=a*a*a 表示連續(xù)坐三次航班可以到達的城市：

5、其中，第 i 行描述從城市 i 出發(fā)，可以到達各個城市的情況，若能到達第 j 個城市，記 a(i,j)=1，否則 a(i,j)=0，規(guī)定 a(i,i)=0 (其中 i=1,2,3,4)。如第 2 行表示：從城市 2出發(fā)可以到達城市 3 和城市 4 而不能到達城市 1 和 2。31222012200002211a2444123300003433at2多次轉(zhuǎn)機到達的城市分析矩陣 a3 的第二行，可以得出：某人從城市 2 出發(fā)，連續(xù)坐三次航班可以到達城市 2、3 和城市4，不能到達城市 1，而到達城市 3 和城市 4 的方法各有兩種。 不難看出，轉(zhuǎn)機兩次以下的航線的航路矩陣為at2= a+ a2 +

6、 a3程序為：a=0,1,1,1; 0,0,1,1; 0,0,0,0; 1,1,0,0;at2=a+a2+a3例例 4 行列式的幾何應用行列式的幾何應用二階行列式的幾何意義是兩個二維向量構成的平行四邊形的面積，三階行列式的幾何意義是三個 3 維向量構成的平行六面體的體積。如下圖所示，用 matlab 軟件來實現(xiàn)面積和體積的運算。 由向量 和 所構成的平行四邊形的面積為行列式 的絕對值。計算的matlab語句為： s=abs(a1*b2-a2*b1)如果給出的是三角形三個頂點坐標a1,b1， a2,b2，a3,b3，求該三角形面積，則有：matlab寫成s=abs(det(a2-a1,b2-b1

7、; a3-a1,b3-b1)11,a bu22,a bv1122abab212131310.5aa bbsabsaa bb平行四邊形面積計算 多邊形可以劃分為多個三角形來計算。 先對三角形面積計算構成一個函數(shù)程序； 這個子程序名為：cal_area3(a,b,c) a,b,c為三個頂點的二維坐標向量 凸多邊形面積只需多次調(diào)用這個函數(shù)程序； 例如五邊形abcde，可由s5= cal_area3(a,b,c)+ cal_area3(a,c,d)+ cal_area3(a,d,e)求得。也可由多邊形面積子程序cal_arean(a)計算。擴展至多邊形面積計算matlab 程序程序function s

8、=cal_area3(a,b,c)% a,b,c 應為同形的 2 維行向量或列向量，% 格式檢驗語句略去 ab=b-a; % 計算向量abac=c-a; % 計算向量acif size(ab)=1,2 % 判讀向量ab是否為行向量 a=ab;ac; % 構造矩陣aelse a=ab,ac;ends=abs(det(a)/2; % 根據(jù)公式計算三角形面積成藥1號成藥2號成藥3號成藥4號成藥5號成藥6號成藥7號成藥a10214122038100b1201225356055c531105140d79255154735e012255336f255355355550g94172523925h651610

9、103510i821200620例例 5 藥方配置問題藥方配置問題 （1）某醫(yī)院要購買這 7 種特效藥，但藥廠的第 3號和第 6 號特效藥已經(jīng)賣完，請問能否用其它特效藥配制出這兩種脫銷的藥品。分析：即 3, 6 向量與其他向量是否線性相關 （2）現(xiàn)在該醫(yī)院想用這 7 種中草藥配制三種新的特效藥，下表為新藥所需的成分質(zhì)量 (單位: 克) 。請問如何配制。分析：這是新藥向量與原來藥向量是否線性相關的問 題。問題及分析思路1號新藥2號新藥3號新藥a4016288b6214167c14278d4410251e53607f5015580g7111838h416821i145230新藥的成分要求新藥的成分

10、要求 u1=10;12;5;7;0;25;9;6;8; u2=2;0;3;9;1;5;4;5;2; u3=14;12;11;25;2;35;17;16;12; u4=12;25;0;5;25;5;25;10;0; u5=20;35;5;15;5;35;2;10;0; u6=38;60;14;47;33;55;39;35;6; u7=100;55;0;35;6;50;25;10;20; u1=u1,u2,u3,u4,u5,u6,u7 v1,r=rref(u1)問題 (1) 的 matlab 程序運行結(jié)果運行結(jié)果v1 = 1 0 1 0 0 0 0 0 1 2 0 0 3 0 0 0 0 1 0

11、1 0 0 0 0 0 1 1 0 0 0 0 0 0 0 1 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 0r = 1 2 4 5 7可見這七種特效藥是“相關的”, 3、6 兩種藥可用其它 5種藥線性配制出來, 但第1 、2 、 4、5 、7 種藥“無關”。因此，8，9 兩種藥可以配出，第 10 種藥則不能配出。v2 = 1 0 1 0 0 0 0 1 3 0 0 1 2 0 0 3 0 3 4 0 0 0 0 1 0 1 0 2 2 0 0 0 0 0 1 1 0 0 0 0 0 0 0 0 0 0 1 0 1 0 0 0 0

12、 0 0 0 0 0 0 1 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 0s = 1 2 4 5 7 10為求第二個問題, 把 3 種新藥與 7 種原藥組成矩陣 u2，求 rref，得：假設一個城市的總?cè)丝跀?shù)是固定不變，但人口的分布情況變化如下：每年都有 5 的市區(qū)居民搬到郊區(qū)；而有 15 的郊區(qū)居民搬到市區(qū)。若開始有 700000 人口居住在市區(qū)，300000 人口居住在郊區(qū)。請分析：（1）10 年后市區(qū)和郊區(qū)的人口各是多少？（2）30 年后、50 年后市區(qū)和郊區(qū)的人口各是多少？（3）分析（2）中數(shù)據(jù)相似的原因。例例

13、6 人口遷徙問題人口遷徙問題,nnnxyx110.950.150.050.85nnnnnnxxyyxy1110.950.150.050.85nnnnnnxxyyxax解 這個問題可以用矩陣乘法來描述。令人口變量其中 xn 為市區(qū)人口所占比例，yn 為郊區(qū)人口所占比例。在 n+1年的人口分布狀態(tài)為：用矩陣乘法可寫成：000700000,300000 xyx2nnn-1n-20x = ax= a x= a xa=0.95,0.15;0.05,0.85;x0=700000;300000;x10=a10*x0 開始市區(qū)和郊區(qū)的人口數(shù)為 可以得到 n 年后市區(qū)和郊區(qū)的人口分布：因此 (1) 10 年后的

14、人口可用程序計算如下：1010107.44637446301.0005,2.5537255370 xeyx運行結(jié)果為：故市區(qū)和郊區(qū)人口數(shù)約為：744630和255370。 kkk-1-1-1-1a = pp ppppp p 00kk-1kx= a x = p p x1122kkk k 無限增加時間 n，市區(qū)和郊區(qū)人口之比將趨向常數(shù)0.75/0.25。為了弄清為什么它趨向于一個穩(wěn)態(tài)值，需要可以求 ak, 為此可先將 a 對角化, 然后求其冪。對角矩陣的冪次可以化為元素的冪次 余下很容易計算。-1a = pp令其中為對角矩陣，則有% 分析 n 年后城市人口分布cleara=0.95,0.15; 0

15、.05,0.85;x0=700000; 300000;p,lambda=eig(a);syms n % 定義符號變量 nxn=p*lamda.n*inv(p)*x0 matlab 程序顯然, 隨 n 增大 (4/5)n 趨近于零, 而 xn 趨于75000050000*4/525000050000*4/5nnxn運行結(jié)果為：750000250000 x01234y-270210-75例例 7 多項式插值與擬合多項式插值與擬合234401234( )pxaa xa xa xa x求: (1) 過這五個點作一個四次多項式函數(shù)22012( )pxaa xa x(2) 請根據(jù)這五個點，擬合一個二次多項

16、式函數(shù)下表給出了平面坐標系中五個點的坐標。并求 x=5 時的函數(shù)值 p4(5)。用 matlab 繪制多項式函數(shù) p4(x) 的曲線、已知點及插值點 (5, p4(5)。并用 matlab 繪制 p2(x) 的曲線及已知的五個點。23401234234012342340123423401234234012340000271111022222133330444475aaaaaaaaaaaaaaaaaaaaaaaaa aa = y234023412342234323442710000011111,21122220133337514444aaaaaaay其中矩陣：解：(1) 根據(jù)已知條件，把五個點的

17、坐標值分別代入四次多項式函數(shù)，可以得到如下線性方程組：系數(shù)矩陣 a 的行列式為范德蒙 (vandermonde) 行列式，且五個坐標點的橫坐標各不相同，則該行列式不等于零，所以方程組有唯一解。 matlab 程序：x=0;1;2;3;4; % 輸入已知點坐標y=-27;0;21;0;-75;a=x.0, x.1, x.2, x.3, x.4; % 構造 vandermonde 矩陣a=ay; % 得到適定方程組的唯一解 a運行程序，得到a(1)=-27, a(2)=12, a(3)=26, a(4)=-12, a(5)= 1.2012201220122012201200271102221330

18、4475aaaaaaaaaaaaaaa aa = y 把五個點的坐標值分別代入二次多項式函數(shù)，可以得到如下線性方程組：22021222271000111,21122013375144aaaaya其中，(2) 多項式擬合要解一個超定方程該方程組有三個未知數(shù)，但有五個方程，進一步分析可以得到該方程組無解，即不存在一個二次多項式曲線剛好能過已知的五個點。matlab 軟件提供了一個利用最小二乘法解決超定方程組解的方法。求系數(shù)的公式也是 a = ay，以找到一條二次曲線來近似地描述已知 5 個點的變化情況。對比插值和擬合的曲線如下圖用平面坐標系中的一個閉合圖形來描述剛體，用一個矩陣 x 來表示它。x

19、的一列表示剛體一個頂點的坐標。為了使圖形閉合，x 的最后一列和第一列相同；為了實現(xiàn)剛體的平移運算，給矩陣 x 添加元素值都為 1 的一行，使 x 為 3n 矩陣。1210cossin001,sincos0 ,001001cttctt12mry = mxy = rx且，若有矩陣：則可以證明，矩陣 y1 是剛體 x 沿 x 軸正方向平移 c1，沿 y 軸正方向平移 c2 后的結(jié)果；矩陣 y2 是剛體 x 以坐標原點為中心逆時針轉(zhuǎn)動 t 弧度的結(jié)果。例例 8 剛體的平面運動剛體的平面運動x04610853.56.1 6.53.220y014140011664.54.500實 例用下列數(shù)據(jù)表示字母 a: 對 a 進行以下平面運動, 并繪制移動前后的圖形。(1) 向上移動 15, 向左移動 30;(2) 逆時針轉(zhuǎn)動 /3;(3) 先逆時針轉(zhuǎn)動3 /4, 然后向上平移 30, 向右平移 20。04610853.56.16.53.220014140011664.54.500111111111111x解 構造剛體矩陣 x，平移矩陣及轉(zhuǎn)動矩陣。m1=100010-30201m2=10001020301r1=3