1、線性代數應用實例求插值多項式右表給出函數 f (t)上4個點的值，試求三次插值多項式p(t )aoait a2t 2 a3t 3，并求f (1.5)的近似值。解：令三次多項式函數p(t)_ao a1t a2 t2 a3t 3過表中已知的 4點，可以得到四元線性方程組:tio123f(t i)31o-16aoaoa1a2a3 -oaoao_2a1 4a2 J3a1' 9a28a3 127 a3 -6應該用計算機求解了，鍵入:對于四元方程組，筆算就很費事了。>>A=1,0,0,0;1,1,1,1;1,2,4,8;1,3,9,27, b=3;0;-1;6, s=rref(A,b)

2、 得到x =1-2-2得到 _ao3,a12,a22,a31，三次多項函數為p(t)= = +3 2t 2t 2 t3，故 f (1.5)近似等于p(1.5)-3 2(1.5) 2(1.5)2(1.5)3 1.125。在一般情況下，當給出函數f (t)在 n+1 個點 ti (i 1,2,IB, 1)上的值f (ti)時，就可以用n次多項式p(t) ao a1ta2t二f 2 I " ant n對f (t )進行插值。在數字信號處理中的應用-數字濾波器系統函數數字濾波器的網絡結構圖實際上也是一種信號流圖。它的特點在于所有的相加節點都限 定為雙輸入相加器；另外，數字濾波器器件有一個遲延

3、一個節拍的運算，它也是一個線性算子， 它的標注符號為z1。根據這樣的結構圖，也可以用類似于例7.4的方法，求它的輸入輸出之間的傳遞函數，在數字信號處 理中稱為系統函數。圖1表示了某個數字濾波器的結構圖，4*1 *產11/4j.1/4I X3T1J3/8u2 x 1y圖1某數字濾波器結構圖現在要求出它的系統函數，即輸出y與輸入u之比。先在它的三個中間節點上標注信號的 名稱x1,x2,x3，以便對每個節點列寫方程。由于遲延算子 z 1不是數，要用符號代替，所以取q_z j按照圖示情況，可以寫出xi qx2 曲2u311X2 二 1X3 治u84丿 4X3 X1寫成矩陣形式為-1 ! oXiX =X

4、2| 二 0X3 J J- 101r n.七 1. 1X1p1 111 ；qX2J-貞84丿I04X301ux = Qx - Pu經過移項后，系統函數現在可以列寫計算系統函數的syms qQ(1,2)q; Q(2,3)=3/8*qQ(3,3)=0;P=2;1/4;0W二in v(eye (3)Q)*P程序運行的結果為W可以寫成：W = x/u = in v(I - Q)* PMATLAB 程序 ea705，%規定符號變量r/4; Q(3,1)=1;%給非零元素賦值%給右下角元素 Q( 3,3)賦值后，矩陣中未賦值元素都自動置零%給P賦值%用信號流圖求傳遞函數的公式W = 一16/(一8 3*q

5、A2 _2*q) _2*q/(飛3*qA2*q)一2*(3*q 一2)/( E 3*qA2 2*q) 一2/( _8 3*qA22*q)_ . ” bJk一16/(-8 3*qA2 一 2*q) _2*q/( _83*qA2 _2*q)我們關心的是以yx3作為輸出的系統函數，故再鍵入pretty(W(3)整理后得到y -1&-2qq 8z 18W (3)22 廿 +.- -f -令u 8 3q 2q 1.5q q 41.5z 2 z 14用線性代數方法的好處是適用于任何復雜系統，并能用計算機解決問題。信號與系統課程中的應用線性時不變系統的零輸入響應描述n階線性時不變,a1a 2dtnd

6、tLTI )連續系統的微分方程為汨護 dyd_u芽斗常duanan 1 y b1bmbm 1 u,n > mdtdtmdt已知y及其各階導數的初始值為y(0) , y(0) , ?, y(n-1)(0)，求系統的零輸入響應。解：當LTI系統的輸入為零時，其零輸入響應為微分方程的齊次解(即令微分方程等其中 p1， p2，?，y(t )C1ep1tC2ep2t扎 扎 pn是特征方程a1 n+a2 n-？e U|U ol!jp tCne n+ an + an+1 =0的根，它們可用roots(a)號右端為0)，其形式為(設特征根均為單根)C 1 + C2 + ? +Cn = y。P1C1+ P

7、2C2+? + pnCn=D y 0語句求得。各系數 C1，?， Cn由y及其各階導數的初始值來確定。對此有yo = y(0)(D y0表示y的導數的初始值y(0)pin 1Clp 2n 1 C2；3晉 pnn 1 C n D " 1 yo"1P11p2寫成矩陣形式為n -1n14p1p 2 事1S3-1 1C1 _ y o 1"亠p nC 2Dyoq ' :Bfa'ft>:n- 1pnCnJESTr n DyoaamC=VY 0即V C = Yo ,其解為式中C -C1 , C2 川 Cn T ; Yo - yo , Dyol , I ,

8、Dn 1yo T1 1 1p1p 2卩門V = t玉*嚴 * 小n小n小“P r p 11 p z12nV為范德蒙矩陣，在MATLAB的特殊矩陣庫中有vander函數可直接生成。MATLAB 程序 ea703.ma=input('輸入分母系數向量a=a1 ，a2，.=');n=len gth(a)-1;Y0=input('輸入初始條件向量Y0二y0，DyO， D2y0，二');p=roots(a);V=rot90(va nder(p);c= VY0'dt=i nput('dt二')；tf二in put('tf=')t=O:

9、dt:tf; y=zeros(1，len gth(t);for k=1: ny二 y+c(k)*exp(p(k)*t);e ndplot（t ，y），grid程序運行結果用這個通用程序來解一個三階系統， 并輸入a=3，5，7，1;dt=0.2; tf=8;而Y0取1，0，0;0 ，1，0;0 ，0，1運行此程序圖2三階系統的零輸入響應三種情況，用hold on語句使三次運行生成的圖形畫在一幅圖上，得到圖2減肥配方的實現營養每100g食物所含營養(g)減肥所要求的 每日營養量脫脂牛奶大豆面粉乳清蛋白質36511333碳水化合物52347445脂肪071.13設三種食物每 100克中蛋白質、碳水化

10、合物和脂肪的含量如下表，表中還給出了80年代美國流行的劍橋大學醫學院的簡捷營養處方。現在的問題是：如果用這三種食物作為每天的主要食物，那么它們的用量應各取多少？才能全面準確地實現這個營養要求。設脫脂牛奶的用量為X1個單位（100g ），大豆面粉的用量為X2個單位（100g ），乳清的用 量為X3個單位（100g ），表中的三個營養成分列向量為：3651 13 11rnIIa1 =：52 ,a2 二34Lu|a1 二 74,0 7 ；1.1- ria»則它們的組合所具有的營養為3651IE13 X1a1+x2a2 進X3a3 = X152： + X234+ X3740 71.1 1JI

11、一a1使這個合成的營養與劍橋配方的要求相等，就可以得到以下的矩陣方程:36L I520Ax r用MATLAB 解這個問題非常方便，列出程序ag763如下:A=36,51,13;52,34,74;0,7,1.1 b=33;45;3 x=Ab程序執行的結果為:0.2772x -0.39190.2332即脫脂牛奶的用量為27.7g，大豆面粉的用量為39.2g，乳清的用量為23.3g，就能保證所需的綜合 營養量。人口遷徙模型設在一個大城市中的總人口是固定的。人口的分布則因居民在市區和郊區之間遷徙而 變化。每年有 6%的市區居民搬到郊區去住，而有2%的郊區居民搬到市區。假如開始時有30%的居民住在市區，

12、70%的居民住在郊區，問十年后市區和郊區的居民人口比例是多少？30年、50年后又如何？這個問題可以用矩陣乘法來描述。把人口變量用市區和郊區兩個分量表示，即Xck廿Xk丨,其中Xc為市區人口所占比例，xs為郊區人口所占比例，k表示年份的次序。XskX Xcol0.31 在k=0的初始狀態：xo二J != J 。-Xs0 0.7一年以后，市區人口為 矩陣乘法來描述，可寫成：X1Xc10.94 0.020.3 一 AX0 -0.2960s10.06 0.980.70.7040Xc1 = (1-0.02) X c0+0.06x s0，郊區人口 Xs1= 0.02x c0 + (1-0.06)X s0，

13、用此關系可以從初始時間到用下列-k 年，擴展為 Xk|”_ Axk 1 A2 Xk 2 Ak xo ,MATLAB程序進行計算：A=0.94,0.02;0.06,0.98x0=0.3;0.7x1=A*x0,x10=AA10*x0x30=AA30*x0 x50=AA50*x0程序運行的結果為：嚴嚴厲嚴rirriri0.2960 J 0.27170.25411= J 0.2508X1, X10, X30, X500.70400.72830.74590.7492無限增加時間k，市區和郊區人口之比將趨向一組常數0.25/0.75。為了弄清為什么這個過程趨向于一個穩態值，我們改變一下坐標系統。在這個坐標

14、系統中可以更清楚地看到乘以矩陣 A的效果。選U1為穩態向量0.25,0.75 T的任意一個倍數， 令U1二1,3 T和U2二-1,1可以看到，用量的長度，不影響其相角（方向）0.940.02 勺11Au1:L. Ju10.060.98.33_r-1=0.940.02 f"I0.92Au2一 r0.060.9810.92初始向量x（J可以寫成這兩個基向量u1禾和=0.30 =11 一f 1X0J 0.25-1 10.05-0.7031Au2的線性組合;二0.92u 20.25u 10.05U2因此k>27，這第Xk Ak X00.25U10.05(0.82)k U2式中的第二項會

15、隨著k的增大趨向于零。如果只取小數點后兩位，則只要二項就可以忽略不計而得到- tI_沖 kI 0.25Xk k 27Ak X0 0.25u 10.75適當選擇基向量可以使矩陣乘法結果等價于一個簡單的實數乘子，避免相角項出現， 使得問題簡單化。這也是方陣求特征值的基本思想。這個應用問題實際上是所謂馬爾可夫過程的一個類型。所得到的向量序列X1,x2,.,Xk稱為馬爾可夫鏈。馬爾可夫過程的特點是k時刻的系統狀態xk完全可由其前一個時刻的狀態xk-1所決定，與k-1時刻之前的系統狀態無關。交通流的分析某城市有兩組單行道，構成了一個包含四個節點A,B,C,D的十字路口如圖 所示在交通繁忙時段的汽車從外部

16、進出此十字路口的流量(每小時的車流數)標于圖上。現要 求計算每兩個節點之間路段上的交通流量X1,X 2,x 3,x 4。解：在每個節點上，進入和離開的車數應該相等，這就決定了四個流通的方程：節點 A: x 1+450 = X2 +610節點 B: x 2+520 = X3 +480節點 C: x 3+390 = X4+600節點 D: x 4 +640 = X2 +310 將這組方程進行整理，寫成矩陣形式：xi - X2= 160X3 - X4=210FX4=-330其系數增廣矩陣為：1一 1:160 11 一 1:-40A,bf；1A "-1 2101廠11： -330X2 X 3

17、=-40或者直接調用 U0=rref(A,b),可以得出其精簡行階梯形式100100-1 ： 330 -1 :170 j為U0=001-1 :210.0000；0 i用消元法求其行階梯形式，注意這個系數矩陣所代表的意義，它的左邊四列從左至右依次為變量X1,X 2,x 3,x 4 的系數，第五列則是在等式右邊的常數項。把第四列移到等式右邊，可以按行列寫恢復為方程， 其結果為：X1=x 4+330,X2=X 4+170,X3=X 4+210由于最后一行變為全零，這個精簡行階梯形式只有三行有效，也就是說四個方程中有 一個是相依的，實際上只有三個有效方程。方程數比未知數的數目少，即沒有給出足夠的信息來

18、唯一地確定X1,X 2,X 3,和X4。其原因也不難從物理上想象，題目給出的只是進入和離開 這個十字路區的流量，如果有些車沿著這四方的單行道繞圈，那是不會影響總的輸入輸出流量的，但可以全面增加四條路上的流量。所以 不能完全自由，因為規定了這些路段都是單行道， 所以要準確了解這里的交通流情況，還應該在X4被稱為自由變量，實際上它的取值也X1,X 2,x 3,和X4。都不能取負值。X1,X 2,x 3 ,和X4中，再檢測一個變量。價格平衡模型在Leo ntiff 成為諾貝爾獎金獲得者的歷史中，線性代數曾起過重要的作用，我們來看看他的基本思路。假定一個國家或區域的經濟可以分解為n個部門，這些部門都有

19、生產產品或服務的獨立功能。設單列n元向量x是這些n個部門的產出，組成在R n空間的產出向量。先假定該社會是自給自足的經濟，這是一個最簡單的情況。因此各經濟部門生產出的產品，完全被自己部門和其它部門所消費。Leo ntiff提出的第一個問題是，各生產部門的實際產出的價格 p應該是多少，才能使各部門的收入和消耗相等，以維持持續的生產。Leon tiff的輸入輸出模型中的一個基本假定是：對于每個部門，存在著一個在Rn空間單位消耗列向量 vi，它表示第i個部門每產出一個單位(比如100萬美金)產品，由本部門和其他各個部門消耗的百分比。在自給自足的經濟中，這些列向量中所有元素的總和應該 為1。把這n個v

20、i，并列起來，它可以構成一個n x n的系數矩陣，可稱為內部需求矩陣V。舉一個最簡單的例子，假如一個自給自足的經濟體由三個部門組成，它們是煤炭業、電力業和鋼鐵業。它們的單位消耗列向量和銷售收入列向量p如下表：由下列部1每單位輸出的消耗分配銷售價格p門購買(收入)煤炭業電力業鋼鐵業煤炭業0.0.40.6Pc電力業0.60.10.2pe鋼鐵業0.40.50.2ps如果電力業產出了100個單位的產品，有40個單位會被煤炭業消耗， 10個單位被自己消耗，而被鋼鐵業消耗的是50個單位，各行業付出的費用為：0.4Pe r V2= Pe 0.1 I0.5 這就是內部消耗的計算方法，把幾個部門都算上，可以寫出

21、其中V 幣!, Ve ,Vs 丄0.0.40.60.60.10.2_0.40.50.2pcpc 1各部門消耗成本=Pc Vc * peVe * psVs二Vc, Ve,Vs -rpe銷售收入=pe |Ips-Ps于是總的價格平衡方程可以寫成為:p - Vp = 0(I - V ) p =0此等式右端常數項為零，是一個齊次方程。它有非零解的條件是系數行列式等于 零，或者用行階梯簡化來求解。用MATLAB語句寫出其解的表示式：V=0.,0.4,0.6;0.6,0.1,0.2;0.4,0.5,0.2,U0 = rref(eye (3)-V ,zeros(3,1)程序運行的結果為1.0000-0.84851.0000-0.9394這個結果是合理的，簡化行階梯形式只有兩行，說明I-V的秩是2，所以它的行列式必定為零。由于現在有三個變量，只有兩個方程，必定有一個變量可以作為自由變量。記住U0矩陣中各列的意義，它們分別是原方程中 表示的是下列方程：pc, Pe, Ps，的系數，所以簡化行階梯矩陣 U0pc - 0.9394 ppe -