數(shù)學(xué)建模案例分析--線性代數(shù)建模案例(20例)

1、線性代數(shù)建模案例匯編目 錄案例一. 交通網(wǎng)絡(luò)流量分析問(wèn)題1案例二. 配方問(wèn)題4案例三. 投入產(chǎn)出問(wèn)題6案例四. 平板的穩(wěn)態(tài)溫度分布問(wèn)題8案例五. CT圖像的代數(shù)重建問(wèn)題10案例六. 平衡結(jié)構(gòu)的梁受力計(jì)算12案例七. 化學(xué)方程式配平問(wèn)題14案例八. 互付工資問(wèn)題16案例九. 平衡價(jià)格問(wèn)題18案例十. 電路設(shè)計(jì)問(wèn)題20案例十一. 平面圖形的幾何變換22案例十二. 太空探測(cè)器軌道數(shù)據(jù)問(wèn)題24案例十三. 應(yīng)用矩陣編制Hill密碼25案例十四. 顯示器色彩制式轉(zhuǎn)換問(wèn)題27案例十五. 人員流動(dòng)問(wèn)題29案例十六. 金融公司支付基金的流動(dòng)31案例十七. 選舉問(wèn)題33案例十八. 簡(jiǎn)單的種群增長(zhǎng)問(wèn)題34案例十九.

2、 一階常系數(shù)線性齊次微分方程組的求解36案例二十. 最值問(wèn)題38附錄 數(shù)學(xué)實(shí)驗(yàn)報(bào)告模板391案例一. 交通網(wǎng)絡(luò)流量分析問(wèn)題城市道路網(wǎng)中每條道路、每個(gè)交叉路口的車(chē)流量調(diào)查，是分析、評(píng)價(jià)及改善城市交通狀況的基礎(chǔ)。根據(jù)實(shí)際車(chē)流量信息可以設(shè)計(jì)流量控制方案，必要時(shí)設(shè)置單行線，以免大量車(chē)輛長(zhǎng)時(shí)間擁堵。【模型準(zhǔn)備】 某城市單行線如下圖所示, 其中的數(shù)字表示該路段每小時(shí)按箭頭方向行駛的車(chē)流量(單位: 輛). 5001234400300100200300x1x2x3x4圖3 某城市單行線車(chē)流量(1) 建立確定每條道路流量的線性方程組.(2) 為了唯一確定未知流量, 還需要增添哪幾條道路的流量統(tǒng)計(jì)? (3) 當(dāng)x

3、4 = 350時(shí), 確定x1, x2, x3的值.(4) 若x4 = 200, 則單行線應(yīng)該如何改動(dòng)才合理? 【模型假設(shè)】 (1) 每條道路都是單行線. (2) 每個(gè)交叉路口進(jìn)入和離開(kāi)的車(chē)輛數(shù)目相等. 【模型建立】 根據(jù)圖3和上述假設(shè), 在, , , 四個(gè)路口進(jìn)出車(chē)輛數(shù)目分別滿(mǎn)足500 = x1 + x2 400 + x1 = x4 + 300 x2 + x3 = 100 + 200 x4 = x3 + 300 【模型求解】根據(jù)上述等式可得如下線性方程組其增廣矩陣(A, b) =由此可得即. 為了唯一確定未知流量, 只要增添x4統(tǒng)計(jì)的值即可. 當(dāng)x4 = 350時(shí), 確定x1 = 250, x

4、2 = 250, x3 = 50.若x4 = 200, 則x1 = 100, x2 = 400, x3 = -100 < 0. 這表明單行線“¬”應(yīng)該改為“®”才合理. 【模型分析】(1) 由(A, b)的行最簡(jiǎn)形可見(jiàn), 上述方程組中的最后一個(gè)方程是多余的. 這意味著最后一個(gè)方程中的數(shù)據(jù)“300”可以不用統(tǒng)計(jì). (2) 由可得, , , 這就是說(shuō)x1, x2, x3, x4這四個(gè)未知量中, 任意一個(gè)未知量的值統(tǒng)計(jì)出來(lái)之后都可以確定出其他三個(gè)未知量的值. Matlab實(shí)驗(yàn)題某城市有下圖所示的交通圖, 每條道路都是單行線, 需要調(diào)查每條道路每小時(shí)的車(chē)流量. 圖中的數(shù)字表示

5、該條路段的車(chē)流數(shù). 如果每個(gè)交叉路口進(jìn)入和離開(kāi)的車(chē)數(shù)相等, 整個(gè)圖中進(jìn)入和離開(kāi)的車(chē)數(shù)相等. 300500150180350160220300100290400150x1x2x3x4x5x6x7x8x9x10x11x12圖4 某城市單行線車(chē)流量(1)建立確定每條道路流量的線性方程組. (2)分析哪些流量數(shù)據(jù)是多余的. (3)為了唯一確定未知流量, 需要增添哪幾條道路的流量統(tǒng)計(jì). 案例二. 配方問(wèn)題在化工、醫(yī)藥、日常膳食等方面都經(jīng)常涉及到配方問(wèn)題. 在不考慮各種成分之間可能發(fā)生某些化學(xué)反應(yīng)時(shí), 配方問(wèn)題可以用向量和線性方程組來(lái)建模. 【模型準(zhǔn)備】一種佐料由四種原料A、B、C、D混合而成. 這種佐料

6、現(xiàn)有兩種規(guī)格, 這兩種規(guī)格的佐料中, 四種原料的比例分別為2:3:1:1和1:2:1:2. 現(xiàn)在需要四種原料的比例為4:7:3:5的第三種規(guī)格的佐料. 問(wèn): 第三種規(guī)格的佐料能否由前兩種規(guī)格的佐料按一定比例配制而成? 【模型假設(shè)】 (1) 假設(shè)四種原料混合在一起時(shí)不發(fā)生化學(xué)變化. (2) 假設(shè)四種原料的比例是按重量計(jì)算的. (3) 假設(shè)前兩種規(guī)格的佐料分裝成袋, 比如說(shuō)第一種規(guī)格的佐料每袋凈重7克(其中A、B、C、D四種原料分別為2克, 3克, 1克, 1克), 第二種規(guī)格的佐料每袋凈重6克(其中A、B、C、D四種原料分別為1克, 2克, 1克, 2克).【模型建立】 根據(jù)已知數(shù)據(jù)和上述假設(shè),

7、 可以進(jìn)一步假設(shè)將x袋第一種規(guī)格的佐料與y袋第二種規(guī)格的佐料混合在一起, 得到的混合物中A、B、C、D四種原料分別為4克, 7克, 3克, 5克, 則有以下線性方程組 【模型求解】上述線性方程組的增廣矩陣(A, b) =, 可見(jiàn) 又因?yàn)榈谝环N規(guī)格的佐料每袋凈重7克, 第二種規(guī)格的佐料每袋凈重6克, 所以第三種規(guī)格的佐料能由前兩種規(guī)格的佐料按7:12的比例配制而成. 【模型分析】(1) 若令a1 = (2, 3, 1, 1)T, a2 = (1, 2, 1, 1)T, b = (4, 7, 5, 3)T, 則原問(wèn)題等價(jià)于“線性方程組Ax = b是否有解”, 也等價(jià)于“b能否由a1, a2線性表示

8、”. (2) 若四種原料的比例是按體積計(jì)算的, 則還要考慮混合前后體積的關(guān)系(未必是簡(jiǎn)單的疊加), 因而最好還是先根據(jù)具體情況將體積比轉(zhuǎn)換為重量比, 然后再按上述方法處理. (3) 上面的模型假設(shè)中的第三個(gè)假設(shè)只是起到簡(jiǎn)化運(yùn)算的作用. 如果直接設(shè)x克第一種規(guī)格的佐料與y克第二種規(guī)格的佐料混合得第三種規(guī)格的佐料, 則有下表表1 混合后四種原料的含量原料 佐料規(guī)格ABCD第一種xxxx第二種yyyy第三種(x + y)(x + y)(x + y)(x + y)因而有如下線性方程組 (*)【模型檢驗(yàn)】把x = 7, y = 12代入上述方程組(*), 則各等式都成立. 可見(jiàn)模型假設(shè)中的第三個(gè)假設(shè)不影

9、響解的正確性. Matlab實(shí)驗(yàn)題蛋白質(zhì)、碳水化合物和脂肪是人體每日必須的三種營(yíng)養(yǎng), 但過(guò)量的脂肪攝入不利于健康.人們可以通過(guò)適量的運(yùn)動(dòng)來(lái)消耗多余的脂肪. 設(shè)三種食物(脫脂牛奶、大豆面粉、乳清)每100克中蛋白質(zhì)、碳水化合物和脂肪的含量以及慢跑5分鐘消耗蛋白質(zhì)、碳水化合物和脂肪的量如下表.表2 三種食物的營(yíng)養(yǎng)成分和慢跑的消耗情況營(yíng)養(yǎng)每100克食物所含營(yíng)養(yǎng)(克)慢跑5分鐘消耗量(克)每日需要的營(yíng)養(yǎng)量(克)牛奶大豆面粉乳清蛋白質(zhì)3651131033碳水化合物5234742045脂肪1071153問(wèn)怎樣安排飲食和運(yùn)動(dòng)才能實(shí)現(xiàn)每日的營(yíng)養(yǎng)需求?案例三. 投入產(chǎn)出問(wèn)題在研究多個(gè)經(jīng)濟(jì)部門(mén)之間的投入產(chǎn)出關(guān)系時(shí)

10、, W. Leontief提出了投入產(chǎn)出模型. 這為經(jīng)濟(jì)學(xué)研究提供了強(qiáng)有力的手段. W. Leontief因此獲得了1973年的Nobel經(jīng)濟(jì)學(xué)獎(jiǎng). 【模型準(zhǔn)備】某地有一座煤礦, 一個(gè)發(fā)電廠和一條鐵路. 經(jīng)成本核算, 每生產(chǎn)價(jià)值1元錢(qián)的煤需消耗0.3元的電; 為了把這1元錢(qián)的煤運(yùn)出去需花費(fèi)0.2元的運(yùn)費(fèi); 每生產(chǎn)1元的電需0.6元的煤作燃料; 為了運(yùn)行電廠的輔助設(shè)備需消耗本身0.1元的電, 還需要花費(fèi)0.1元的運(yùn)費(fèi); 作為鐵路局, 每提供1元運(yùn)費(fèi)的運(yùn)輸需消耗0.5元的煤, 輔助設(shè)備要消耗0.1元的電. 現(xiàn)煤礦接到外地6萬(wàn)元煤的訂貨, 電廠有10萬(wàn)元電的外地需求, 問(wèn): 煤礦和電廠各生產(chǎn)多少才能

11、滿(mǎn)足需求? 【模型假設(shè)】假設(shè)不考慮價(jià)格變動(dòng)等其他因素. 【模型建立】設(shè)煤礦, 電廠, 鐵路分別產(chǎn)出x元, y元, z元?jiǎng)偤脻M(mǎn)足需求. 則有下表表3 消耗與產(chǎn)出情況產(chǎn)出(1元)產(chǎn)出消耗訂單煤電運(yùn)消耗煤00.60.5x0.6y + 0.5z60000電0.30.10.1y0.3x + 0.1y + 0.1z100000運(yùn)0.20.10z0.2x + 0.1y0根據(jù)需求, 應(yīng)該有, 即【模型求解】在Matlab命令窗口輸入以下命令>> A = 1,-0.6,-0.5;-0.3,0.9,-0.1;-0.2,-0.1,1; b = 60000;100000;0;>> x = Ab

12、 Matlab執(zhí)行后得x = 1.0e+005 * 1.9966 1.8415 0.5835可見(jiàn)煤礦要生產(chǎn)1.9966´105元的煤, 電廠要生產(chǎn)1.8415´105元的電恰好滿(mǎn)足需求. 【模型分析】令x =, A =, b =, 其中x稱(chēng)為總產(chǎn)值列向量, A稱(chēng)為消耗系數(shù)矩陣, b稱(chēng)為最終產(chǎn)品向量, 則Ax =根據(jù)需求, 應(yīng)該有x - Ax = b, 即(E - A)x = b. 故x = (E - A)-1b. Matlab實(shí)驗(yàn)題某鄉(xiāng)鎮(zhèn)有甲、乙、丙三個(gè)企業(yè). 甲企業(yè)每生產(chǎn)1元的產(chǎn)品要消耗0.25元乙企業(yè)的產(chǎn)品和0.25元丙企業(yè)的產(chǎn)品. 乙企業(yè)每生產(chǎn)1元的產(chǎn)品要消耗0.65

13、元甲企業(yè)的產(chǎn)品, 0.05元自產(chǎn)的產(chǎn)品和0.05元丙企業(yè)的產(chǎn)品. 丙企業(yè)每生產(chǎn)1元的產(chǎn)品要消耗0.5元甲企業(yè)的產(chǎn)品和0.1元乙企業(yè)的產(chǎn)品. 在一個(gè)生產(chǎn)周期內(nèi), 甲、乙、丙三個(gè)企業(yè)生產(chǎn)的產(chǎn)品價(jià)值分別為100萬(wàn)元, 120萬(wàn)元, 60萬(wàn)元, 同時(shí)各自的固定資產(chǎn)折舊分別為20萬(wàn)元, 5萬(wàn)元和5萬(wàn)元. (1) 求一個(gè)生產(chǎn)周期內(nèi)這三個(gè)企業(yè)扣除消耗和折舊后的新創(chuàng)價(jià)值. (2) 如果這三個(gè)企業(yè)接到外來(lái)訂單分別為50萬(wàn)元, 60萬(wàn)元, 40萬(wàn)元, 那么他們各生產(chǎn)多少才能滿(mǎn)足需求? 案例四. 平板的穩(wěn)態(tài)溫度分布問(wèn)題在熱傳導(dǎo)的研究中, 一個(gè)重要的問(wèn)題是確定一塊平板的穩(wěn)態(tài)溫度分布. 根據(jù)定律, 只要測(cè)定一塊矩形平板

14、四周的溫度就可以確定平板上各點(diǎn)的溫度. 圖8 一塊平板的溫度分布圖【模型準(zhǔn)備】如圖9所示的平板代表一條金屬梁的截面. 已知四周8個(gè)節(jié)點(diǎn)處的溫度(單位°C), 求中間4個(gè)點(diǎn)處的溫度T1, T2, T3, T4. T1T2T3T410080908060506050圖9 一塊平板的溫度分布圖【模型假設(shè)】假設(shè)忽略垂直于該截面方向上的熱傳導(dǎo), 并且每個(gè)節(jié)點(diǎn)的溫度等于與它相鄰的四個(gè)節(jié)點(diǎn)溫度的平均值. 【模型建立】根據(jù)已知條件和上述假設(shè), 有如下線性方程組【模型求解】將上述線性方程組整理得. 在Matlab命令窗口輸入以下命令>> A = 4,-1,-1,0;-1,4,0,-1;-1,

15、0,4,-1;0,-1,-1,4; b = 190;140;140;100;>> x = Ab; xMatlab執(zhí)行后得ans = 82.9167 70.8333 70.8333 60.4167可見(jiàn)T1 = 82.9167, T2 = 70.8333, T3 = 70.8333, T4 = 60.4167. 參考文獻(xiàn)陳懷琛, 高淑萍, 楊威, 工程線性代數(shù), 北京: 電子工業(yè)出版社, 2007. 頁(yè)碼: 15-16. Matlab實(shí)驗(yàn)題假定下圖中的平板代表一條金屬梁的截面, 并忽略垂直于該截面方向上的熱傳導(dǎo). 已知平板內(nèi)部有30個(gè)節(jié)點(diǎn), 每個(gè)節(jié)點(diǎn)的溫度近似等于與它相鄰的四個(gè)節(jié)點(diǎn)溫度

16、的平均值. 設(shè)4條邊界上的溫度分別等于每位同學(xué)學(xué)號(hào)的后四位的5倍, 例如學(xué)號(hào)為16308209的同學(xué)計(jì)算本題時(shí), 選擇Tl = 40, Tu = 10, Tr = 0, Td = 45. TuT1T5TlTlTdT2T6T7T10TrTrTuT26T30TdT27TuTrTdTl圖10 一塊平板的溫度分布圖(1) 建立可以確定平板內(nèi)節(jié)點(diǎn)溫度的線性方程組. (2) 用Matlab軟件求解該線性方程組. (3) 用Matlab中的函數(shù)mesh繪制三維平板溫度分布圖. 12案例五. CT圖像的代數(shù)重建問(wèn)題X射線透視可以得到3維對(duì)象在2維平面上的投影, CT則通過(guò)不同角度的X射線得到3維對(duì)象的多個(gè)2維

17、投影, 并以此重建對(duì)象內(nèi)部的3維圖像. 代數(shù)重建方法就是從這些2維投影出發(fā), 通過(guò)求解超定線性方程組, 獲得對(duì)象內(nèi)部3維圖像的方法. 圖11雙層螺旋CT 圖12 CT圖像這里我們考慮一個(gè)更簡(jiǎn)單的模型, 從2維圖像的1維投影重建原先的2維圖像. 一個(gè)長(zhǎng)方形圖像可以用一個(gè)橫豎均勻劃分的離散網(wǎng)格來(lái)覆蓋, 每個(gè)網(wǎng)格對(duì)應(yīng)一個(gè)像素, 它是該網(wǎng)格上各點(diǎn)像素的均值. 這樣一個(gè)圖像就可以用一個(gè)矩陣表示,其元素就是圖像在一點(diǎn)的灰度值(黑白圖像). 下面我們以3´3圖像為例來(lái)說(shuō)明. 表4 消耗與產(chǎn)出情況3´3圖像各點(diǎn)的灰度值水平方向上的疊加值x1 = 1x2 = 0x3 = 0x1 + x2 +

18、 x3 = 1x4 = 0x5 = 0.5x6 = 0.5x4 + x5 + x6 = 1x7 = 0.5x8 = 0x9 = 1x7 + x8 + x9 = 1.5豎直方向上的疊加值x1 + x4 + x7 = 1.5x2 + x5 + x8 = 0.5x3 + x6 + x9 = 1.5每個(gè)網(wǎng)格中的數(shù)字xi代表其灰度值, 范圍在0, 1內(nèi). 0表示白色, 1表示黑色, 0.5表示灰色. 如果我們不知道網(wǎng)格中的數(shù)值, 只知道沿豎直方向和水平方向的疊加值, 為了確定網(wǎng)格中的灰度值, 可以建立線性方程組(含有6個(gè)方程, 9個(gè)未知數(shù)) 顯然該方程組的解是不唯一的, 為了重建圖像, 必須增加疊加值.

19、 如我們?cè)黾訌挠疑戏降阶笙路降寞B加值, 則方程組將增加5個(gè)方程x1 = 1,x2 + x4 = 0, x3 + x5 + x7 = 1, x6 + x8 = 0.5, x9 = 1, 和上面的6個(gè)方程放在一起構(gòu)成一個(gè)含有11個(gè)方程, 9個(gè)未知數(shù)的線性方程組. 【模型準(zhǔn)備】設(shè)3´3圖像中第一行3個(gè)點(diǎn)的灰度值依次為x1, x2, x3, 第二行3個(gè)點(diǎn)的灰度值依次為x4, x5, x6, 第三行3個(gè)點(diǎn)的灰度值依次為x7, x8, x9. 沿豎直方向的疊加值依次為1.5, 0.5, 1.5, 沿水平方向的疊加值依次為1, 1, 1.5, 沿右上方到左下方的疊加值依次為1, 0, 1, 0.5

20、, 1. 確定x1, x2, , x9的值. 【模型建立】由已知條件可得(含有11個(gè)方程, 9個(gè)未知數(shù)的)線性方程組【模型求解】在Matlab命令窗口輸入以下命令>> A = 1,1,1,0,0,0,0,0,0;0,0,0,1,1,1,0,0,0;0,0,0,0,0,0,1,1,1;1,0,0,1,0,0,1,0,0;0,1,0,0,1,0,0,1,0;0,0,1,0,0,1,0,0,1;1,0,0,0,0,0,0,0,0;0,1,0,1,0,0,0,0,0;0,0,1,0,1,0,1,0,0;0,0,0,0,0,1,0,1,0;0,0,0,0,0,0,0,0,1;>>

21、 b = 1;1;1.5;1.5;0.5;1.5;1;0;1;0.5;1;>> x = Ab; xMatlab執(zhí)行后得Warning: Rank deficient, rank = 8 tol = 4.2305e-015.ans =1.0000 0.0000 0 -0.0000 0.5000 0.5000 0.5000 -0.0000 1.0000可見(jiàn)上述方程組的解不唯一. 其中的一個(gè)特解為x1 = 1, x2 = 0, x3 = 0, x4 = 0, x5 = 0.5, x6 = 0.5, x7 = 0.5, x8 = 0, x9 = 1.【模型分析】上述結(jié)果表明, 僅有三個(gè)方向

22、上的疊加值還不夠.可以再增加從左上方到右下方的疊加值. 在實(shí)際情況下, 由于測(cè)量誤差, 上述線性方程組可能是超定的. 這時(shí)可以將超定方程組的近似解作為重建的圖像數(shù)據(jù). Matlab實(shí)驗(yàn)題給定一個(gè)3´3圖像的2個(gè)方向上的灰度疊加值: 沿左上方到右下方的灰度疊加值依次為0.8, 1.2, 1.7, 0.2, 0.3; 沿右上方到左下方的灰度疊加值依次為0.6, 0.2, 1.6, 1.2, 0.6. (1) 建立可以確定網(wǎng)格數(shù)據(jù)的線性方程組, 并用Matlab求解. (2) 將網(wǎng)格數(shù)據(jù)乘以256, 再取整, 用Matlab繪制該灰度圖像. 案例六. 平衡結(jié)構(gòu)的梁受力計(jì)算在橋梁、房頂、鐵塔

23、等建筑結(jié)構(gòu)中, 涉及到各種各樣的梁. 對(duì)這些梁進(jìn)行受力分析是設(shè)計(jì)師、工程師經(jīng)常做的事情. 圖14 埃菲爾鐵塔局部下面以雙桿系統(tǒng)的受力分析為例, 說(shuō)明如何研究梁上各鉸接點(diǎn)處的受力情況. 【模型準(zhǔn)備】在圖15所示的雙桿系統(tǒng)中, 已知桿1重G1 = 200牛頓, 長(zhǎng)L1 = 2米, 與水平方向的夾角為q1 = p/6, 桿2重G2 = 100牛頓, 長(zhǎng)L2 = 米, 與水平方向的夾角為q2 = p/4. 三個(gè)鉸接點(diǎn)A, B, C所在平面垂直于水平面. 求桿1, 桿2在鉸接點(diǎn)處所受到的力.ABC桿1桿2p/6p/4圖15雙桿系統(tǒng)【模型假設(shè)】假設(shè)兩桿都是均勻的. 在鉸接點(diǎn)處的受力情況如圖16所示. 【模

24、型建立】對(duì)于桿1:水平方向受到的合力為零, 故N1 = N3, 豎直方向受到的合力為零, 故N2 + N4 = G1, 以點(diǎn)A為支點(diǎn)的合力矩為零, 故(L1sinq1)N3 + (L1cosq1)N4 = (L1cosq1)G1. ABC桿1桿2CN1N2N4N3N7N8N5N6G1G2圖16 兩桿受力情況對(duì)于桿2類(lèi)似地有N5 = N7, N6 = N8 + G2, (L2sinq2)N7 = (L2cosq2)N8 + (L2cosq2)G2.此外還有N3 = N7, N4 = N8. 于是將上述8個(gè)等式聯(lián)立起來(lái)得到關(guān)于N1, N2, , N8的線性方程組:【模型求解】在Matlab命令窗口

25、輸入以下命令>> G1=200; L1=2; theta1=pi/6; G2=100; L2=sqrt(2); theta2=pi/4;>> A = 1,0,-1,0,0,0,0,0;0,1,0,1,0,0,0,0;0,0,L1*sin(theta1),L1*cos(theta1),0,0,0,0;0,0,0,0,1,0,-1,0;0,0,0,0,0,1,0,-1;0,0,0,0,0,0,L2*sin(theta2),-L2*cos(theta2);0,0,1,0,0,0,-1,0;0,0,0,1,0,0,0,-1;>> b = 0;G1;0.5*L1*co

26、s(theta1)*G1;0;G2;0.5*L2*cos(theta2)*G2;0;0;>> x = Ab; xMatlab執(zhí)行后得ans =95.0962 154.9038 95.0962 45.0962 95.0962 145.0962 95.0962 45.0962 【模型分析】最后的結(jié)果沒(méi)有出現(xiàn)負(fù)值, 說(shuō)明圖16中假設(shè)的各個(gè)力的方向與事實(shí)一致. 如果結(jié)果中出現(xiàn)負(fù)值, 則說(shuō)明該力的方向與假設(shè)的方向相反. 參考文獻(xiàn)陳懷琛, 高淑萍, 楊威, 工程線性代數(shù), 北京: 電子工業(yè)出版社, 2007. 頁(yè)碼: 157- 158. Matlab實(shí)驗(yàn)題有一個(gè)平面結(jié)構(gòu)如下所示, 有13條梁(

27、圖中標(biāo)號(hào)的線段)和8個(gè)鉸接點(diǎn)(圖中標(biāo)號(hào)的圈)聯(lián)結(jié)在一起. 其中1號(hào)鉸接點(diǎn)完全固定, 8號(hào)鉸接點(diǎn)豎直方向固定, 并在2號(hào), 5號(hào)和6號(hào)鉸接點(diǎn)上, 分別有圖示的10噸, 15噸和20噸的負(fù)載. 在靜平衡的條件下,任何一個(gè)鉸接點(diǎn)上水平和豎直方向受力都是平衡的. 已知每條斜梁的角度都是45º.(1) 列出由各鉸接點(diǎn)處受力平衡方程構(gòu)成的線性方程組. (2) 用Matlab軟件求解該線性方程組, 確定每條梁受力情況. 圖17 一個(gè)平面結(jié)構(gòu)的梁案例七. 化學(xué)方程式配平問(wèn)題在用化學(xué)方法處理污水過(guò)程中, 有時(shí)會(huì)涉及到復(fù)雜的化學(xué)反應(yīng). 這些反應(yīng)的化學(xué)方程式是分析計(jì)算和工藝設(shè)計(jì)的重要依據(jù). 在定性地檢測(cè)出

28、反應(yīng)物和生成物之后,可以通過(guò)求解線性方程組配平化學(xué)方程式. 【模型準(zhǔn)備】某廠廢水中含KCN, 其濃度為650mg/L. 現(xiàn)用氯氧化法處理, 發(fā)生如下反應(yīng): KCN + 2KOH + Cl2 = KOCN + 2KCl + H2O.投入過(guò)量液氯, 可將氰酸鹽進(jìn)一步氧化為氮?dú)? 請(qǐng)配平下列化學(xué)方程式: KOCN + KOH + Cl2 = CO2 + N2 + KCl + H2O. (注: 題目摘自福建省廈門(mén)外國(guó)語(yǔ)學(xué)校2008-2009學(xué)年高三第三次月考化學(xué)試卷)【模型建立】設(shè)x1KOCN x2KOH x3Cl2 = x4CO2 x5N2 x6KCl x7H2O,則, 即【模型求解】在Matlab

29、命令窗口輸入以下命令>> A = 1,1,0,0,0,-1,0;1,1,0,-2,0,0,-1;1,0,0,-1,0,0,0;1,0,0,0,-2,0,0;0,1,0,0,0,0,-2;0,0,2,0,0,-1,0;>> x = null(A,r); format rat, xMatlab執(zhí)行后得ans = 1 2 3/2 1 1/2 3 1可見(jiàn)上述齊次線性方程組的通解為x = k(1, 2, 3/2, 1, 1/2, 3, 1)T. 取k = 2得x = (2, 4, 3, 2, 1, 6, 2)T. 可見(jiàn)配平后的化學(xué)方程式如下2KOCN + 4KOH + 3Cl2

30、= 2CO2 + N2 + 6KCl + 2H2O. 【模型分析】利用線性方程組配平化學(xué)方程式是一種待定系數(shù)法. 關(guān)鍵是根據(jù)化學(xué)方程式兩邊所涉及到的各種元素的量相等的原則列出方程. 所得到的齊次線性方程組Ax = q中所含方程的個(gè)數(shù)等于化學(xué)方程式中元素的種數(shù)s, 未知數(shù)的個(gè)數(shù)就是化學(xué)方程式中的項(xiàng)數(shù)n. 當(dāng)r(A) = n -1時(shí), Ax = q的基礎(chǔ)解系中含有1個(gè)(線性無(wú)關(guān)的)解向量. 這時(shí)在通解中取常數(shù)k為各分量分母的最小公倍數(shù)即可. 例如本例中1, 2, 3/2, 1, 1/2, 3, 1分母的最小公倍數(shù)為2, 故取k = 2. 當(dāng)r(A) £ n -2時(shí), Ax = q的基礎(chǔ)解

31、系中含有2個(gè)以上的線性無(wú)關(guān)的解向量. 這時(shí)可以根據(jù)化學(xué)方程式中元素的化合價(jià)的上升與下降的情況, 在原線性方程組中添加新的方程.Matlab實(shí)驗(yàn)題配平下列反應(yīng)式(1) FeS + KMnO4 + H2SO4 K2SO4 + MnSO4 + Fe2(SO4)3 + H2O + S(2) Al2(SO4)3 + Na2CO3 + H2O Al(OH)3+ CO2+ Na2SO4 案例八. 互付工資問(wèn)題互付工資問(wèn)題是多方合作相互提供勞動(dòng)過(guò)程中產(chǎn)生的. 比如農(nóng)忙季節(jié), 多戶(hù)農(nóng)民組成互助組, 共同完成各戶(hù)的耕、種、收等農(nóng)活. 又如木工, 電工, 油漆工等組成互助組, 共同完成各家的裝潢工作. 由于不同工種

32、的勞動(dòng)量有所不同, 為了均衡各方的利益, 就要計(jì)算互付工資的標(biāo)準(zhǔn). 【模型準(zhǔn)備】現(xiàn)有一個(gè)木工, 電工, 油漆工. 相互裝修他們的房子, 他們有如下協(xié)議:(1) 每人工作10天(包括在自己家的日子), (2) 每人的日工資一般的市價(jià)在6080元之間, (3) 日工資數(shù)應(yīng)使每人的總收入和總支出相等. 表5 工作天數(shù)在誰(shuí)家 工人木工電工油漆工木工家216電工家451油漆工家443求每人的日工資. 【模型假設(shè)】假設(shè)每人每天工作時(shí)間長(zhǎng)度相同. 無(wú)論誰(shuí)在誰(shuí)家干活都按正常情況工作, 既不偷懶, 也不加班. 【模型建立】設(shè)木工, 電工, 油漆工的日工資分別為x, y, z元, 則由下表表6 各家應(yīng)付工資和各人

33、應(yīng)得收入在誰(shuí)家 工人木工電工油漆工各家應(yīng)付工資木工家2x1y6z2x + y + 6z 電工家4x5y1z4x + 5y + z油漆工家4x4y3z4x + 4y + 3z各人應(yīng)得收入10x10y10z可得, 即【模型求解】在Matlab命令窗口輸入以下命令>> A = -8,1,6;4,-5,1;4,4,-7;>> x = null(A,r); format rat, xMatlab執(zhí)行后得ans = 31/36 8/9 1可見(jiàn)上述齊次線性方程組的通解為x = k(31/36, 8/9, 1)T. 因而根據(jù)“每人的日工資一般的市價(jià)在6080元之間”可知60 £

34、;k <k < k £ 80, 即 £ k £ 80. 也就是說(shuō), 木工, 電工, 油漆工的日工資分別為k元, k元, k元, 其中£ k £ 80. 為了簡(jiǎn)便起見(jiàn), 可取k = 72, 于是木工, 電工, 油漆工的日工資分別為62元, 64元, 72元. 【模型分析】事實(shí)上各人都不必付自己工資, 這時(shí)各家應(yīng)付工資和各人應(yīng)得收入如下表7 各家應(yīng)付工資和各人應(yīng)得收入在誰(shuí)家 工人木工電工油漆工各家應(yīng)付工資木工家01y6zy + 6z 電工家4x01z4x + z油漆工家4x4y04x + 4y 個(gè)人應(yīng)得收入8x5y7z由此可得, 即可見(jiàn)

35、這樣得到的方程組與前面得到的方程組是一樣的. Matlab實(shí)驗(yàn)題甲, 乙, 丙三個(gè)農(nóng)民組成互助組, 每人工作6天(包括為自己家干活的天數(shù)), 剛好完成他們?nèi)思业霓r(nóng)活, 其中甲在甲, 乙, 丙三家干活的天數(shù)依次為: 2, 2.5, 1.5; 乙在甲, 乙, 丙三家各干2天活, 丙在甲, 乙, 丙三家干活的天數(shù)依次為: 1.5, 2, 2.5. 根據(jù)三人干活的種類(lèi), 速度和時(shí)間, 他們確定三人不必相互支付工資剛好公平. 隨后三人又合作到鄰村幫忙干了2天(各人干活的種類(lèi)和強(qiáng)度不變), 共獲得工資500元. 問(wèn)他們應(yīng)該怎樣分配這500元工資才合理? 案例九. 平衡價(jià)格問(wèn)題為了協(xié)調(diào)多個(gè)相互依存的行業(yè)的

36、平衡發(fā)展, 有關(guān)部門(mén)需要根據(jù)每個(gè)行業(yè)的產(chǎn)出在各個(gè)行業(yè)中的分配情況確定每個(gè)行業(yè)產(chǎn)品的指導(dǎo)價(jià)格, 使得每個(gè)行業(yè)的投入與產(chǎn)出都大致相等. 【模型準(zhǔn)備】假設(shè)一個(gè)經(jīng)濟(jì)系統(tǒng)由煤炭、電力、鋼鐵行業(yè)組成, 每個(gè)行業(yè)的產(chǎn)出在各個(gè)行業(yè)中的分配如下表所示: 表7 行業(yè)產(chǎn)出分配表產(chǎn)出分配購(gòu)買(mǎi)者煤炭電力鋼鐵00.40.6煤炭0.60.10.2電力0.40.50.2鋼鐵每一列中的元素表示占該行業(yè)總產(chǎn)出的比例. 求使得每個(gè)行業(yè)的投入與產(chǎn)出都相等的平衡價(jià)格. 【模型假設(shè)】假設(shè)不考慮這個(gè)系統(tǒng)與外界的聯(lián)系. 【模型建立】把煤炭、電力、鋼鐵行業(yè)每年總產(chǎn)出的價(jià)格分別用x1, x2, x3表示, 則, 即. 【模型求解】在Matlab

37、命令窗口輸入以下命令>> A = 1,-0.4,-0.6;-0.6,0.9,-0.2;-0.4,-0.5,0.8;>> x = null(A,r); format short, xMatlab執(zhí)行后得ans = 0.9394 0.8485 1.0000可見(jiàn)上述齊次線性方程組的通解為x = k(0.9394, 0.8485, 1)T. 這就是說(shuō), 如果煤炭、電力、鋼鐵行業(yè)每年總產(chǎn)出的價(jià)格分別0.9394億元, 0.8485億元, 1億元, 那么每個(gè)行業(yè)的投入與產(chǎn)出都相等. 【模型分析】實(shí)際上, 一個(gè)比較完整的經(jīng)濟(jì)系統(tǒng)不可能只涉及三個(gè)行業(yè), 因此需要統(tǒng)計(jì)更多的行業(yè)間的分配數(shù)

38、據(jù). Matlab實(shí)驗(yàn)題假設(shè)一個(gè)經(jīng)濟(jì)系統(tǒng)由煤炭、石油、電力、鋼鐵、機(jī)械制造、運(yùn)輸行業(yè)組成, 每個(gè)行業(yè)的產(chǎn)出在各個(gè)行業(yè)中的分配如下表所示: 表8 行業(yè)產(chǎn)出分配表產(chǎn)出分配購(gòu)買(mǎi)者煤炭石油電力鋼鐵制造運(yùn)輸000.20.10.20.2煤炭000.10.10.20.1石油0.50.10.10.20.10.1電力0.40.10.200.10.4鋼鐵00.10.30.600.2制造0.10.70.100.40運(yùn)輸每一列中的元素表示占該行業(yè)總產(chǎn)出的比例. 求使得每個(gè)行業(yè)的投入與產(chǎn)出都相等的平衡價(jià)格. 案例十. 電路設(shè)計(jì)問(wèn)題電路是電子元件的神經(jīng)系統(tǒng). 參數(shù)的計(jì)算是電路設(shè)計(jì)的重要環(huán)節(jié). 其依據(jù)來(lái)自?xún)蓚€(gè)方面: 一是客

39、觀需要, 二是物理學(xué)定律. 圖22 USB擴(kuò)展板【模型準(zhǔn)備】假設(shè)圖23中的方框代表某類(lèi)具有輸入和輸出終端的電路. 用記錄輸入電壓和輸入電流(電壓v以伏特為單位, 電流i以安培為單位), 用記錄輸出電壓和輸入電流. 若= A, 則稱(chēng)矩陣A為轉(zhuǎn)移矩陣.輸入終端v1輸出終端v2i1i2電路圖23 具有輸入和輸出終端的電子電路圖圖24給出了一個(gè)梯形網(wǎng)絡(luò), 左邊的電路稱(chēng)為串聯(lián)電路, 電阻為R1(單位: 歐姆). 右邊的電路是并聯(lián)電路, 電路R2. 利用歐姆定理和楚列斯基定律, 我們可以得到串聯(lián)電路和并聯(lián)電路的轉(zhuǎn)移矩陣分別是和v1v2i1i2R1v3i2i3R2串聯(lián)電路 并聯(lián)電路圖24 梯形網(wǎng)絡(luò)設(shè)計(jì)一個(gè)梯

40、形網(wǎng)絡(luò), 其轉(zhuǎn)移矩陣是.【模型假設(shè)】假設(shè)導(dǎo)線的電阻為零. 【模型建立】設(shè)A1和A2分別是串聯(lián)電路和并聯(lián)電路的轉(zhuǎn)移矩陣, 則輸入向量x先變換成A1x, 再變換到A2(A1x). 其中A2A1 = 就是圖22中梯形網(wǎng)絡(luò)的轉(zhuǎn)移矩陣. 于是, 原問(wèn)題轉(zhuǎn)化為求R1, R2的值使得=. 【模型求解】由=可得. 根據(jù)其中的前兩個(gè)方程可得R1 = 8, R2 = 2. 把R1 = 8, R2 = 2代入上面的第三個(gè)方程確實(shí)能使等式成立. 這就是說(shuō)在圖22中梯形網(wǎng)絡(luò)中取R1 = 8, R2 = 2即為所求. 【模型分析】若要求的轉(zhuǎn)移矩陣改為, 則上面的梯形網(wǎng)絡(luò)無(wú)法實(shí)現(xiàn). 因?yàn)檫@時(shí)對(duì)應(yīng)的方程組是. 根據(jù)前兩個(gè)方程

41、依然得到R1 = 8, R2 = 2, 但把R1 = 8, R2 = 2代入上第三個(gè)方程卻不能使等式成立. 練習(xí)題根據(jù)基爾霍夫回路電路定律(各節(jié)點(diǎn)處流入和流出的電流強(qiáng)度的代數(shù)和為零, 各回路中各支路的電壓降之和為零), 列出下圖所示電路中電流i1, i2, i3所滿(mǎn)足的線性方程組, 并用矩陣形式表示: E1E2R1R2R3R4R5i1i2i3 圖25 簡(jiǎn)單的回路案例十一. 平面圖形的幾何變換隨著計(jì)算機(jī)科學(xué)技術(shù)的發(fā)展, 計(jì)算機(jī)圖形學(xué)的應(yīng)用領(lǐng)域越來(lái)越廣, 如仿真設(shè)計(jì)、效果圖制作、動(dòng)畫(huà)片制作、電子游戲開(kāi)發(fā)等. 圖形的幾何變換, 包括圖形的平移、旋轉(zhuǎn)、放縮等, 是計(jì)算機(jī)圖形學(xué)中經(jīng)常遇到的問(wèn)題. 這里暫

42、時(shí)只討論平面圖形的幾何變換. 【模型準(zhǔn)備】平面圖形的旋轉(zhuǎn)和放縮都很容易用矩陣乘法實(shí)現(xiàn), 但是圖形的平移并不是線性運(yùn)算, 不能直接用矩陣乘法表示. 現(xiàn)在要求用一種方法使平移、旋轉(zhuǎn)、放縮能統(tǒng)一用矩陣乘法來(lái)實(shí)現(xiàn). 【模型假設(shè)】設(shè)平移變換為(x, y) ® (x+a, y+b)旋轉(zhuǎn)變換(繞原點(diǎn)逆時(shí)針旋轉(zhuǎn)q角度)為(x, y) ® (xcosq - ysinq, xsinq + ycosq)放縮變換(沿x軸方向放大s倍, 沿y軸方向放大t倍)為 (x, y) ® (sx, ty)【模型求解】R2中的每個(gè)點(diǎn)(x, y)可以對(duì)應(yīng)于R3中的(x, y, 1). 它在xOy平面上方

43、1單位的平面上. 我們稱(chēng)(x, y, 1)是(x, y)的齊次坐標(biāo). 在齊次坐標(biāo)下, 平移變換(x, y) ® (x+a, y+b)可以用齊次坐標(biāo)寫(xiě)成(x, y, 1) ® (x+a, y+b, 1).于是可以用矩陣乘積=實(shí)現(xiàn). 旋轉(zhuǎn)變換(x, y) ® (xcosq - ysinq, xsinq + ycosq)可以用齊次坐標(biāo)寫(xiě)成(x, y, 1) ® (xcosq - ysinq, xsinq + ycosq, 1).于是可以用矩陣乘積=實(shí)現(xiàn). 放縮變換 (x, y) ® (sx, ty)可以用齊次坐標(biāo)寫(xiě)成(x, y, 1) ®

44、(sx, ty, 1). 于是可以用矩陣乘積=實(shí)現(xiàn). 【模型分析】由上述求解可以看出, R2中的任何線性變換都可以用分塊矩陣乘以齊次坐標(biāo)實(shí)現(xiàn), 其中A是2階方陣. 這樣, 只要把平面圖形上點(diǎn)的齊次坐標(biāo)寫(xiě)成列向量, 平面圖形的每一次幾何變換, 都可通過(guò)左乘一個(gè)3階變換矩陣來(lái)實(shí)現(xiàn). 參考文獻(xiàn)David C. Lay, 線性代數(shù)及其應(yīng)用, 沈復(fù)興, 傅鶯鶯等譯, 北京: 人民郵電出版社, 2009. 頁(yè)碼: 139-141. Matlab實(shí)驗(yàn)題在Matlab命令窗口輸入以下命令>>clear all, clc,>>t = 1,3,5,11,13,15*pi/8;>>

45、;x = sin(t); y=cos(t); >>fill(x,y,'r');>>grid on;>>axis(-2.4, 2.4, -2, 2)運(yùn)行后得圖25. 圖26 Matlab繪制的圖形(1) 寫(xiě)出該圖形每個(gè)頂點(diǎn)的齊次坐標(biāo); (2) 編寫(xiě)Matlab程序, 先將上面圖形放大0.9倍; 再逆時(shí)針旋轉(zhuǎn); 最后進(jìn)行橫坐標(biāo)加0.8, 縱坐標(biāo)減1的圖形平移. 分別繪制上述變換后的圖形. 案例十二. 太空探測(cè)器軌道數(shù)據(jù)問(wèn)題太空航天探測(cè)器發(fā)射以后, 可能需要調(diào)整以使探測(cè)器處在精確計(jì)算的軌道里. 雷達(dá)監(jiān)測(cè)到一組列向量x1, , xk, 它們給出了不同

46、時(shí)刻探測(cè)器的實(shí)際位置與預(yù)定軌道之間的偏差的信息.圖28 火星探測(cè)器【模型準(zhǔn)備】令Xk = x1, , xk. 在雷達(dá)進(jìn)行數(shù)據(jù)分析時(shí)需要計(jì)算出矩陣Gk = XkXkT. 一旦接收到數(shù)據(jù)向量xk+1, 必須計(jì)算出新矩陣Gk+1. 因?yàn)閿?shù)據(jù)向量到達(dá)的速度非常快, 隨著k的增加, 直接計(jì)算的負(fù)擔(dān)會(huì)越來(lái)越重. 現(xiàn)需要給出一個(gè)算法, 使得計(jì)算Gk的負(fù)擔(dān)不會(huì)因?yàn)閗的增加而加重. 【模型求解】因?yàn)镚k = XkXkT = x1, , xk=, Gk+1 = Xk+1= Xk, xk+1= XkXkT + xk+1= Gk + xk+1, 所以一旦接收到數(shù)據(jù)向量xk+1, 只要計(jì)算xk+1, 然后把它與上一步計(jì)

47、算得到的Gk 相加即可. 這樣計(jì)算Gk的負(fù)擔(dān)不會(huì)因?yàn)閗的增加而加重. 【模型分析】計(jì)算機(jī)計(jì)算加法的時(shí)間與計(jì)算乘法的時(shí)間相比可以忽略不計(jì). 因此在考慮計(jì)算矩陣乘積的負(fù)擔(dān)時(shí), 只要考察乘法的次數(shù)就可以了. 設(shè)xk的維數(shù)是n, 則Xk = x1, , xk是n´k的矩陣, Gk = XkXkT是n´n的矩陣. 直接計(jì)算Gk = XkXkT需要做n2k次乘法. 因而計(jì)算的負(fù)擔(dān)會(huì)隨著k的增加而增加. 但是對(duì)于每一個(gè)k, 計(jì)算xk始終只要做n2次乘法. Matlab實(shí)驗(yàn)題用Matlab編寫(xiě)一個(gè)程序用于處理這個(gè)問(wèn)題. 案例十三. 應(yīng)用矩陣編制Hill密碼密碼學(xué)在經(jīng)濟(jì)和軍事方面起著極其重要

48、的作用. 現(xiàn)代密碼學(xué)涉及很多高深的數(shù)學(xué)知識(shí). 這里無(wú)法展開(kāi)介紹. 信 源 噪聲信 宿 信 道 攻擊 解 密 請(qǐng)求重傳 加 冗 編 碼 加 密 識(shí) 錯(cuò) 糾 錯(cuò) 圖29 保密通信的基本模型密碼學(xué)中將信息代碼稱(chēng)為密碼, 尚未轉(zhuǎn)換成密碼的文字信息稱(chēng)為明文, 由密碼表示的信息稱(chēng)為密文. 從明文到密文的過(guò)程稱(chēng)為加密, 反之為解密. 1929年, 希爾(Hill)通過(guò)線性變換對(duì)待傳輸信息進(jìn)行加密處理, 提出了在密碼史上有重要地位的希爾加密算法. 下面我們略去一些實(shí)際應(yīng)用中的細(xì)節(jié), 只介紹最基本的思想.【模型準(zhǔn)備】若要發(fā)出信息action, 現(xiàn)需要利用矩陣乘法給出加密方法和加密后得到的密文, 并給出相應(yīng)的解密

49、方法. 【模型假設(shè)】(1) 假定每個(gè)字母都對(duì)應(yīng)一個(gè)非負(fù)整數(shù), 空格和26個(gè)英文字母依次對(duì)應(yīng)整數(shù)026(見(jiàn)下表). 表9 空格及字母的整數(shù)代碼表空格ABCDEFGHIJKLM012345678910111213NOPQRSTUVWXYZ14151617181920212223242526(2) 假設(shè)將單詞中從左到右, 每3個(gè)字母分為一組, 并將對(duì)應(yīng)的3個(gè)整數(shù)排成3維的行向量, 加密后仍為3維的行向量, 其分量仍為整數(shù). 【模型建立】設(shè)3維向量x為明文, 要選一個(gè)矩陣A使密文y = xA, 還要確保接收方能由y準(zhǔn)確地解出x. 因此A必須是一個(gè)3階可逆矩陣. 這樣就可以由y = xA得x = yA-

50、1. 為了避免小數(shù)引起誤差, 并且確保y也是整數(shù)向量, A和A-1的元素應(yīng)該都是整數(shù). 注意到, 當(dāng)整數(shù)矩陣A的行列式= ±1時(shí), A-1也是整數(shù)矩陣. 因此原問(wèn)題轉(zhuǎn)化為 (1) 把a(bǔ)ction翻譯成兩個(gè)行向量: x1, x2. (2) 構(gòu)造一個(gè)行列式= ±1的整數(shù)矩陣A(當(dāng)然不能取A = E).(3) 計(jì)算x1A和x2A. (4) 計(jì)算A-1. 【模型求解】(1) 由上述假設(shè)可見(jiàn)x1 = (1, 3, 20), x2 = (9, 15, 14). (2) 對(duì)3階單位矩陣E =進(jìn)行幾次適當(dāng)?shù)某醯茸儞Q(比如把某一行的整數(shù)被加到另一行, 或交換某兩行), 根據(jù)行列式的性質(zhì)可知,

51、 這樣得到的矩陣A的行列式為1或-1. 例如A =, |A| = -1. (3) y1 = x1A = (1, 3, 20)= (67, 44, 43), y2 = Ax2 = (9, 15, 14)= (81, 52, 43). (4) 由(A, E) =可得A-1 =. 這就是說(shuō), 接收方收到的密文是67, 44, 43, 81, 52, 43. 要還原成明文, 只要計(jì)算(67, 44, 43)A-1和(81, 52, 43)A-1, 再對(duì)照表9“翻譯”成單詞即可.【模型分析】如果要發(fā)送一個(gè)英文句子, 在不記標(biāo)點(diǎn)符號(hào)的情況下, 我們?nèi)匀豢梢园丫渥?含空格)從左到右每3個(gè)字符分為一組(最后不

52、足3個(gè)字母時(shí)用空格補(bǔ)上). 【模型檢驗(yàn)】(67, 44, 43) A-1 = (1, 3, 20), (81, 52, 43)A-1 = (9, 15, 14). 參考文獻(xiàn)楊威, 高淑萍, 線性代數(shù)機(jī)算與應(yīng)用指導(dǎo), 西安: 西安電子科技大學(xué)出版社, 2009. 頁(yè)碼: 98-102. Matlab實(shí)驗(yàn)題按照上面的加密方法, 設(shè)密文為: 112, 76, 57, 51, 38, 18, 84, 49, 49, 68, 41, 32, 83, 55, 37, 70, 45, 25, 問(wèn)恢復(fù)為原來(lái)的信息是什么? 案例十四. 顯示器色彩制式轉(zhuǎn)換問(wèn)題彩色顯示器使用紅(R)、綠(G)和藍(lán)(B)光的疊成效應(yīng)

53、生成顏色. 顯示器屏幕的內(nèi)表面由微粒象素組成, 每個(gè)微粒包括三個(gè)熒光點(diǎn): 紅、綠、藍(lán). 電子槍位于屏幕的后方, 向屏幕上每個(gè)點(diǎn)發(fā)射電子束. 計(jì)算機(jī)從圖形應(yīng)用程序或掃描儀發(fā)出數(shù)字信號(hào)到電子槍, 這些信號(hào)控制電子槍設(shè)置的電壓強(qiáng)度. 不同 RGB 的強(qiáng)度組合將產(chǎn)生不同的顏色. 電子槍由電磁石幫助瞄準(zhǔn)以確?？焖倬_地屏幕刷新. 圖30 彩色顯示器的工作原理顏色模型規(guī)定一些屬性或原色, 將顏色分解成不同屬性的數(shù)字化組合. 這就色彩制式的轉(zhuǎn)換問(wèn)題.【模型準(zhǔn)備】觀察者在屏幕上實(shí)際看到的色彩要受色彩制式和屏幕上熒光點(diǎn)數(shù)量的影響. 因此每家計(jì)算機(jī)屏幕制造商都必須在(R, G, B)數(shù)據(jù)和國(guó)際通行的CIE色彩標(biāo)準(zhǔn)

54、之間進(jìn)行轉(zhuǎn)換, CIE標(biāo)準(zhǔn)使用三原色, 分別稱(chēng)為X, Y和Z. 針對(duì)短余輝熒光點(diǎn)的一類(lèi)典型轉(zhuǎn)換是=.計(jì)算機(jī)程序把用CIE數(shù)據(jù)(X, Y, Z)表示的色彩信息流發(fā)送到屏幕. 求屏幕上的電子槍將這些數(shù)據(jù)轉(zhuǎn)換成(R, G, B)數(shù)據(jù)的方程.【模型建立】令A(yù) =, a =, b =, 則Aa = b. 現(xiàn)在要根據(jù)CIE數(shù)據(jù)(X, Y, Z)計(jì)算對(duì)應(yīng)的(R, G, B)數(shù)據(jù), 就是等式Aa = b中A和b已知, 求a. 如果A是可逆矩陣, 則由Aa = b可得a = A-1b. 【模型求解】在Matlab命令窗口輸入以下命令>> A = 0.61,0.29,0.15;0.35,0.59,0.063;0.04,0.12,0.787; >> if det(A)=0, 'A不可逆' elseif &