matlab詳細(xì)學(xué)習(xí)資料.ppt

第三章 行列式,3.1 行列式的定義 3.2 行列式的性質(zhì)及應(yīng)用 3.3 克萊姆（Cramer）法則 3.4 行列式的計(jì)算 3.5 應(yīng)用實(shí)例 3.6 習(xí)題,3.1 行列式的定義,3.1.1 二、三階行列式的定義 引入記號(hào)： ，稱(chēng)它為二階行列式 其值規(guī)定為：,把 的連線稱(chēng)為二階行列式的主對(duì)角線， 把 的連線稱(chēng)為二階行列式的副對(duì)角線， 那么二階行列式的值就等于主對(duì)角線上元素 的乘積減去副對(duì)角線上元素的乘積。 例3.2 在平面上有一個(gè)平行四邊形OACB， A、B兩點(diǎn)的坐標(biāo)分別為： 、 ，如圖 3.1所示，求平行四邊形OACB的面積。,圖3.1 二階行列式等價(jià)于平行四邊形面積 解：過(guò)點(diǎn)A做x軸垂線，交x軸于點(diǎn)E；過(guò)點(diǎn)B 做平行x軸直線與過(guò)點(diǎn)C做平行y軸直線相交于 點(diǎn)D。顯然可以得到三角形CDB和三角形AEO 全等，則有： （3-2）,根據(jù)二階行列式的定義，該平行四邊形的面 積剛好是以A、B兩點(diǎn)坐標(biāo)所構(gòu)成的二階行列 式： 例3.3 求下面三元線性方程組的解：,解：利用消元法可以得到： （3-3） 當(dāng) 之前的系數(shù)不為零時(shí)，可以解出 的值 但這個(gè)結(jié)果很難記憶，為此引進(jìn)三階行列式 的定義，我們稱(chēng)： 是一個(gè)三階行列式,其值規(guī)定為： （3-4） 圖3.2給出了它的圖示計(jì)算規(guī)則（稱(chēng)為沙路 法）。 圖3.2 三階行列式的計(jì)算規(guī)則,有了三階行列式的定義，我們可以把式（3- 3）寫(xiě)為： 當(dāng)方程組（3-2）的系數(shù)行列式 時(shí)，可以得到它的解。,3.1.2 n階行列式的定義 把三階行列式定義式（3-4）改寫(xiě)為如下形 式： 則有： （3-5）,定義3.1 在n階行列式中，劃去元素 所在的 第i行和第j列元素后，余下的元素按原來(lái)次序 構(gòu)成的n-1階行列式，稱(chēng)為元素 的余子式 記作 ，稱(chēng)為元素 的代數(shù)余子式。 根據(jù)定義3.1，可以把式（3-5）變?yōu)椋?定義3.2 由 個(gè)數(shù)組成的 階行列式 是一個(gè)算式，當(dāng) 時(shí)， ；當(dāng) 時(shí)， （3-6）,3.1.3 行列式定義的進(jìn)一步討論 定義3.3 由n個(gè)自然數(shù)1、2、3、n組成 的一個(gè)有序數(shù)組，稱(chēng)為一個(gè)n元排列。 例如，1 2 3、132、213、231、312、321都 是3元排列。 在n元排列的n！個(gè)排列中，123n是唯一一 個(gè)按從小到大排列的n元排列，稱(chēng)為標(biāo)準(zhǔn)排列 （或自然排列）,定義3.4 一個(gè)排列中任兩個(gè)數(shù)，如果排在前 面的數(shù)大于排在后面的數(shù)，則稱(chēng)這兩個(gè)數(shù)構(gòu) 成一個(gè)逆序。一個(gè)排列中逆序的總數(shù)，稱(chēng)為 這個(gè)排列的逆序數(shù)。 排列的逆序數(shù)記為 例如，五元排列25341，1和2、5、3、4有四 個(gè)逆序，4和5有一個(gè)逆序，3和5有一個(gè)逆 序，則排列25341的逆序數(shù)為4116。,定義3.5 （行列式的另一種定義方法）： 由 個(gè)數(shù)組成的 階行列式： (3-7) 其中 是一個(gè)n元排列， 表示對(duì) 所有n元排列（n！個(gè)）求和。,例3.4 寫(xiě)出四階行列式中含有 的項(xiàng)。 解：四階行列式共有24項(xiàng)，其中含有 的項(xiàng) 為： ，我們只要分析列標(biāo)排列 1x2y的各種情況，顯然有1324和1423兩種情況，1324逆序數(shù)為1，1423逆序數(shù)為2，則四階行列式中含有 的項(xiàng)為： 和 。,3.1.4 矩陣與行列式的關(guān)系 矩陣是一個(gè)數(shù)表，而行列式是一個(gè)算式，即 它是一個(gè)值。 在比較兩個(gè)矩陣是否相等時(shí)，不僅要求兩個(gè) 矩陣同型，而且要求兩個(gè)矩陣所有對(duì)應(yīng)元素 相等。 而兩個(gè)行列式是否相等，只需分析其值是否 相等。 矩陣是由一對(duì)方括號(hào)（或圓括號(hào)）括起，而 行列式是由一對(duì)豎線括起。,矩陣的行數(shù)和列數(shù)不做任何限制，而行列式 的行數(shù)與列數(shù)必須相等。 當(dāng)討論的矩陣A是方陣時(shí)，把A的一對(duì)方括號(hào) 去掉，加上一對(duì)豎線，就變成了行列式，我 們把這個(gè)行列式稱(chēng)為方陣A的行列式， 記作： 或 。,例3.5 證明n階下三角矩陣 的行列式 。 證明：用數(shù)學(xué)歸納法證明，當(dāng) 、 時(shí)，結(jié)論顯然成立。 假設(shè)結(jié)論對(duì) 階下三角行列式成立，則由定義3.2得：,右端行列式是 階下三角行列式，根據(jù)歸 納假設(shè)得： 同理可證，n階對(duì)角矩陣的行列式（也稱(chēng)n階 對(duì)角行列式）,3.1.5 行列式按行（列）展開(kāi) 定理3.1 n階行列式D等于它的任一行（列） 的各元素與其對(duì)應(yīng)的代數(shù)余子式乘積之和， 即： 或：,例3.6 計(jì)算行列式 解：此行列式第3列只有一個(gè)非零元素，故應(yīng) 把行列式按第三列展開(kāi)。得到的三階行列式 中的第3列又只有一個(gè)非零元素，再得：,3.2 行列式的性質(zhì)及應(yīng)用,3.2.1 行列式的性質(zhì) 性質(zhì)1 行列式 與它的轉(zhuǎn)置行列式 相等。 行列式的轉(zhuǎn)置和矩陣的轉(zhuǎn)置概念相同。 如：,性質(zhì)2 互換行列式的兩行（列），行列式變 號(hào)。 如： 推論1 如果行列式有兩行（列）完全相同， 則此行列式為零。 推論2 n階行列式D的任一行（列）的各元素 與另一行（列）對(duì)應(yīng)元素的代數(shù)余子式乘積 之和等于零，即：,或： 例3.7 已知四階行列式 ， 求 ，（其中 為行列式D 的代數(shù)余子式）。 解：可以先求出行列式D的第四行各元素的代 數(shù)余子式，然后再進(jìn)一步求得題目的答案。 也可以利用代數(shù)余子式的性質(zhì)來(lái)分析此題。,構(gòu)造行列式 ，行列式 按第四 行的展開(kāi)式，剛好就是,性質(zhì)3 用數(shù)k乘以行列式D，等于D中某一行 （列）的所有元素同乘以數(shù)k。 如下等式中，把數(shù)3乘到了行列式的第二列 中： 推論1 行列式某一行（列）中所有元素的公 因子可以提到行列式的外面。,推論2 如果行列式的任意兩行（列）對(duì)應(yīng)元 素成比例，則行列式為零。 下列行列式的第一行和第三行所有對(duì)應(yīng)元素 成比例，故知： 性質(zhì)4 行列式可以按某一行（列）拆分成兩 個(gè)行列式之和。,具體拆分方法用4階行列式說(shuō)明如下： 性質(zhì)5 把行列式的某一行（列）的各元素乘 以同一數(shù)后加到另一行（列）對(duì)應(yīng)的元素上 去，行列式值不變。如：,3.2.2 方陣行列式的性質(zhì) 設(shè)A、B為n階方陣，k是數(shù)，根據(jù)行列式的性 質(zhì)可以得到方陣的行列式有如下性質(zhì)： （1） （2） （3） （4） （5）,3.2.3 方陣可逆的充要條件 定義3.6 設(shè)矩陣 ,且 的代數(shù) 余子式為 ，則稱(chēng)矩 為 的伴 隨矩陣。記為 ，或 。,伴隨矩陣的重要性質(zhì)： 定理3.2 n階方陣 為可逆矩陣的充要條件 是 。當(dāng) 可逆時(shí)， 。 證：充分性，當(dāng) 時(shí)， 知 故結(jié)論成立； 必要性， 設(shè) 可逆，有 ，兩邊同取 行列式 ，故,推論 若 和 為同階方陣，且滿足 ， 則 ，即矩陣 和矩陣 互逆。 例3.8 判斷三階方陣 ，是否可逆； 若可逆，求 解：因?yàn)?，所以 可逆。 中各元素的代數(shù)余子式分別為,則： 例3.9 設(shè) 為n階可逆方陣， 證明（1） 也是可逆矩陣且 （2） 證：（1）因?yàn)榫仃?為可逆方陣，則 又根據(jù)伴隨矩陣的性質(zhì),知 ，故 是可逆矩 陣且 （2）對(duì)等式 兩邊取行列式， 有 又因?yàn)榫仃?為可逆方陣，則 故有,3.3 克萊姆（Cramer）法則,討論用行列式來(lái)求解含有n個(gè)方程n個(gè)變量的 線性方程組 （3-7） 方程組（3-7）也可以寫(xiě)成矩陣形式： （3-8）,其中： 行列式 ，稱(chēng)為方程組（3-7）的系數(shù)行 列式。 定理3.3（克萊姆法則） 若方程組（3-7）的 系數(shù)行列式 ，則該方程組有唯一解： ， ， ， （3-9）,其中 是把 中第 列的元素用 方程組右端的常數(shù)項(xiàng)代替后所得到的n階行列 式，即： 第 列 定理3.3的逆否定理為：如果線性方程組（3- 7）無(wú)解或有超過(guò)一個(gè)以上的解，則它的系數(shù) 行列式必為零。,把常數(shù)項(xiàng)全為零的線性方程組 稱(chēng)為齊次線性方程組； 把常數(shù)項(xiàng)不全為零的線性方程組 稱(chēng)為非齊次線性方程組。 推論1 對(duì)于nn齊次線性方程組 ，當(dāng) 系數(shù)行列式 時(shí)， 只有一個(gè)零解。 推論2 若nn齊次線性方程組 ，有 非零解，則必有 。,例3.10 已知齊次線性方程組 有非零解，問(wèn) 應(yīng)取何值？ 解：根據(jù)推論2，知方程組系數(shù)行列式必為 零，故有： 得： 或,3.4 行列式的計(jì)算,3.4.1 行列式的筆算技巧 主要的方法是把矩陣變換為行階梯形（三角 形），然后計(jì)算其主對(duì)角線元素的連乘積； 其次是充分利用含零元素較多的行或列進(jìn)行 展開(kāi)； 其他還有加邊法、公式法、遞推法、數(shù)學(xué)歸 納法等等。,例3.12 計(jì)算四階行列式 解：利用行列式性質(zhì)把行列式化為三角行列 式（性質(zhì)法、三角化法）,例3.13 證明： 證：利用行列式性質(zhì)及行列式按列展開(kāi)（性 質(zhì)法、展開(kāi)法）,此例中的四階行列式，稱(chēng)為四階范德蒙 （Vander Monde）行列式, n階范德蒙行列 式為：,例3.14 計(jì)算n階行列式（空白處都為零） ： 解：其中只有n個(gè)非0元素，這n個(gè)元素之積正 是行列式唯一的非零項(xiàng)，再由列下標(biāo)全排 列（n-1,n-2,2,1,n）的逆序數(shù)確定該項(xiàng)的 正負(fù)。,例3.15 計(jì)算5階行列式： 解：由分塊矩陣行列式公式： 則得,例3.16 計(jì)算5階行列式： 解：該行列式稱(chēng)為三對(duì)角行列式，通?？梢?用遞推法來(lái)求解,例3.17 設(shè) ， 均為n階方陣， 求： 解：由于 ， 則有：,例3.18 設(shè)矩陣 ， 矩陣 滿 足： ，其中 為單位矩陣， 是 的伴隨矩陣，求 。 解：由于 ，則存在 ，且有 即有： 兩邊取行列式，有： 而 則,例3.19 設(shè) ， 為三階方陣且 ， ，求 。 解：根據(jù)分塊矩陣的乘法概念，有： 則,3.4.2 用MATLAB計(jì)算行列式 考慮的問(wèn)題主要是計(jì)算速度和計(jì)算精度問(wèn)題 一初等矩陣的行列式 對(duì)于第一類(lèi)初等矩陣E1，即行交換變換，它 的行列式等于-1。 det(E1)=-1 (3-11) 對(duì)于第二類(lèi)初等矩陣E2，即行數(shù)乘變換，它 的行列式等于k。 det(E2)=k (3-12） 對(duì)于第三類(lèi)初等矩陣E3，即行的乘加變換， 它的行列式仍等于1。det(E3)=1 (3-13）,二把方陣變換為上三角矩陣LU分解 如果不考慮出現(xiàn)基準(zhǔn)元素為零或很小的情況 時(shí)，連第一類(lèi)初等變換都用不到。這樣，通 過(guò)N=(n-1)2/2次使用第三類(lèi)初等矩陣E3，就 可以把主對(duì)角線左下方的N個(gè)元素全變?yōu)榱恪?寫(xiě)成 (3-14） 其中U是一個(gè)上三角矩陣，所有的E3矩陣也 是上三角矩陣。,由于初等變換矩陣都是可逆的，其乘積也是 可逆的，設(shè)其逆矩陣為L(zhǎng)，它應(yīng)該是一個(gè)下三 角矩陣，于是此式可寫(xiě)成 (3-15） 這種把矩陣A通過(guò)第三種初等矩陣左乘分解為 一個(gè)下三角矩陣和一個(gè)上三角矩陣乘積的變 換稱(chēng)為L(zhǎng)U變換。 其中下三角矩陣L的行列式為1，因而上三角 矩陣的行列式就等于原矩陣的行列式，即 det(A)= det(L)*det(U)= det(U),在實(shí)際工程中為了保證計(jì)算結(jié)果的精度，計(jì) 算軟件在做行階梯變換時(shí)還是要把基準(zhǔn)元素 取為每列的絕對(duì)值最大項(xiàng)，所以還是使用了 行交換變換。因?yàn)槠湫辛惺降扔?1，每多一 次交換，就改變一次符號(hào)。它并不影響行列 式的絕對(duì)值，但影響其正負(fù)號(hào)。 另外在(3-14）式左端的連乘積中，多了若干 個(gè)交換矩陣。會(huì)使得最后的下三角矩陣L不那 末標(biāo)準(zhǔn)，各行有些顛倒，常常稱(chēng)之為準(zhǔn)下三 角矩陣。,MATLAB提供了矩陣的三角分解函數(shù)lu.m， 其調(diào)用格式為： L,U=lu(A) 它返回的結(jié)果就是一個(gè)準(zhǔn)下三角矩陣L和一個(gè) 上三角矩陣U。因?yàn)檫@個(gè)函數(shù)并不專(zhuān)為行列式 計(jì)算之用，有一定的普遍性，所以它不限定A 為方陣。 另一種調(diào)用格式能同時(shí)給出真正的下三角矩 陣L和交換矩陣P,其形式為： L,U,P = lu(A) 此時(shí)，它滿足 P*A = L*U （3-16）,三求出上三角方陣的行列式 由(3-15）式知道，det(U)決定了det(A)的絕 對(duì)值。因U是一個(gè)上三角矩陣，它的行列式為 其對(duì)角元素的連乘積。 在不計(jì)正負(fù)號(hào)的時(shí)候，可以寫(xiě)出： 用MATLAB語(yǔ)句表示為 D=prod(diag(U),在必須知道行列式的正負(fù)號(hào)時(shí)，必須知道lu分 解過(guò)程中進(jìn)行了多少次交換，每一次交換就 改變一次正負(fù)號(hào)。 lu子程序沒(méi)有給出這個(gè)數(shù)據(jù)，所以解決不了問(wèn) 題。其實(shí)MATLAB已經(jīng)把上述過(guò)程集成在一 起，給出了直接計(jì)算方陣行列式的函數(shù)det.m 其調(diào)用格式為： D=det(A) 這個(gè)函數(shù)要求輸入變?cè)仨毷欠疥?3.5 應(yīng)用實(shí)例,3.5.1 用LU分解計(jì)算行列式 例3.14 用化簡(jiǎn)為三角矩陣的方法求下列矩陣 的行列式 解：列出程序： A10,8,6,4,1;2,5,8,9,4;6,0,9,9,8;5,8,7,4,0;9,4,2,9,1; L,Ulu(A), % 分解為上三角矩陣U和準(zhǔn)下三角矩陣L dU diag(U); % 取上三角矩陣U主對(duì)角線上元素向量,Dprod(dU) % 求主對(duì)角元素的連乘積 程序運(yùn)行的結(jié)果如下： dU 10 4.8 10.625 9.4824 1.2349 D 5.9720e003 5972,3.5.2 行列式奇異性對(duì)計(jì)算精度的影響 例3.15 設(shè)線性方程組 中， 是一個(gè)6階 的hilbert矩陣，就是說(shuō)它的下標(biāo)為(i,j)的元素 值為1/(i+j-1)，系數(shù)A,b1及其增量b2=b1+b 如下： 計(jì)算解x1,x2，分析兩個(gè)解的差與系數(shù)差 之間的關(guān)系。,解：用MATLAB寫(xiě)出程序ea344如下： A=hilb(6), b1=1;2;1;1.732;1;2; b2=1;2;1;1.7321;1;2; x1=inv(A)*b1, x2= inv(A)*b2 dx=x2-x1, db=b2-b1 程序運(yùn)行的結(jié)果為：,由于系數(shù)b的萬(wàn)分之一的變化，引起解x的誤 差dx最大可達(dá)近400。 主要因?yàn)樾辛惺紻=det(A)很接近于零。本題 中的矩陣系數(shù)是hilbert矩陣，它的主要特點(diǎn)就 是行列式很接近于零。這樣的矩陣方程，我 們就稱(chēng)之為病態(tài)的，或很接近于奇異的，對(duì) 它的解是否正確，要保持一定的懷疑。 為了定量地分析解的誤差和可信度，應(yīng)該用 相對(duì)誤差做標(biāo)準(zhǔn)。b的相對(duì)誤差是 x的相對(duì)誤差是,兩者之間是以條件數(shù)(Condition Number)相關(guān) 聯(lián)的，條件數(shù)與行列式有關(guān)，它隨著行列式 的減小而減小： （3-17） 例3.16 設(shè) ，求其逆陣V 解：輸入A的數(shù)據(jù)后，鍵入Vinv(A)，程序 為： A=-16,-4,-6;15,-3,9;18,0,9, V=inv(A),運(yùn)行后得到警告信息： Warning: Matrix is close to singular or badly scaled Results may be inaccurate. RCOND 6.042030e018. det(A)=0，故它是一個(gè)奇異矩陣，其逆不存 在。在用數(shù)值方法求矩陣的逆時(shí)，由于計(jì)算 中存在方法和截?cái)嗾`差，故矩陣是否奇異， 結(jié)果不是絕對(duì)的。為了評(píng)價(jià)矩陣接近“奇異”的 程度MATLAB用了“逆條件數(shù)”作為衡量指標(biāo)。,3.5.3 用逆陣進(jìn)行保密編譯碼 把消息中的英文字母用一個(gè)整數(shù)來(lái)表示。然 后傳送這組整數(shù)。 例如“SEND MONEY”這九個(gè)字母就用下面九 個(gè)數(shù)來(lái)表示； 5,8,10,21,7,2,10,8,3。5代表S，8代表 E，等等。 這種方法是很容易被破譯的。在一個(gè)很長(zhǎng)的 消息中，根據(jù)數(shù)字出現(xiàn)的頻率，往