以反矩阵解联立方裎组.DOC

PAGE PAGE 1 計算機概論作業─聯立方程式與重算 計概作業大多是老師在研究與工作過程中所碰到的例子或問題，也有些是純屬好玩弄出來的。為了準備計概作業，去年暑假期間老師花了約三千大洋從亞馬遜（amazon）網路書店訂購了一批 Excel 工程與數學應用書籍（因為國內沒有相關書籍），看看有沒有新的問題可出。結果發現幾乎所有書上的問題都已想到過（適合各位程度已出現或以後會出現），唯有本次作業中的「重算」問題是以前沒想到過的。這是解決工程科學問題常用的數值方法（作業中的問題當然是最簡單的例子），但必須解超大型聯立方程組，一般需要撰寫程式才能解決。老師總算知道 Excel 也可以用來解這種問題。 步驟 1. 這個作業先要完成如下數學遊戲： 如下圖，有 16 個格子，其中八個已填上數字，其餘空格中的數字必須是相鄰空格（上下左右）的平均，角落四個空格中的數字則是相鄰兩個已有數字空格的平均。 100 100 100 100 100 100 200 200 2. 角落四個數字可以直接求得，中央四個格子中的數字則分別假設為 x1，x2，x3 與 x4（如下圖）。 100 100 100 100 100 x1 x2 100 100 x3 x4 100 150 200 200 150 3. 依題意，我們可建立以下的聯立方程組： 經整理後得： 保留上面帶 0 的各項是為方便後面的工作。 4. 一般用三個方程式所構成的聯立方程式還勉強可用手算，要用手算四個方程式所構成的聯立方程組就有點慘忍了，必須想想別的方法。好在 Excel 提供了方便的方法，這需要矩陣運算觀念，請參考後面「以矩陣解聯立方裎組」一節。 5 如下圖，在適當的位置鍵入矩陣 [A] 與 [B] 的數字（分別為 B6 到 E9 以及 I6 到 I9），並加上文字。 6 於 B12 到 E15 建立 [A]-1： (1) 選取 B12 到 E15（四欄與四列）。 (2) 直接鍵入 = MINVERSE(B6:E9)。（B6:E9 就是 [A] 的範圍；MINVERSE 中的 M 代表 matrix，INVESE 是相反） (3) 按 [Ctrl]+[Shift]+[Enter] 鍵（不要只按 [Enter] 鍵）後便完成了。 (4) 加上文字說明後如下圖。 假如你不小心只按了 [Enter] 鍵，可依以下步驟解決： (1) 繼續選取 B12 到 E15。 (2) 按 [F2] 鍵（編輯鍵）。 (3) 按 [Ctrl]+[Shift]+[Enter] 鍵。 7. 於 C18 到 C21 計算 [A]-1[B]： (1) 選取 C18 到 C21（一欄四列）。 (2) 直接鍵入 = MMULT(B12:E15,I6:I9)。（MULT 是 multiplication代表乘法） (3) 按 [Ctrl]+[Shift]+[Enter] 鍵。 (4) 加上文字說明後如下圖。 8. 於是你得到的答案是 x1 = x2 = 112.5，x3 = x4 = 137.5。 注意 (1) 矩陣的大小是以 (m ?n) 的形式表示，m 是（橫）列數，n是（縱）欄數。前面 [A] 與 [B] 的大小分別是 (4 ? 4) 與 (4 ? 1)。 (2) 若 [A] 與 [B] 兩矩陣能進行 [A][B] = [C] 相乘運算，[A] 的欄數必須等於 [B] 的列數。相乘所得的矩陣 [C] 的列數與 [A] 一樣，欄數與 [B] 一樣。 (3) 不好記，對不對？用這個辦法：[A] 的大小是 (m ? k)，[B] 的大小是 (k ? n)，相乘所得結果 [C] 的大小是 (m ? n)。因此 (m ? k) (k ? n) = (m ? n)。用這個關係驗證作業中所有的運算。 (4) 不要忘記，只有方陣（列數等於欄數的矩陣）才有反矩陣，而且反矩陣是與原矩陣大小相同的方陣。 9. 下面我們要試試看 [A]-1[A] 是否真能得到單位矩陣。以相同的做法在 B24 到 E27 執行 = MMULT(B6:E9,B12:E15)，加上文字說明如下圖。其中有些本來應該是 0 的儲存格出現如 2.01817E-17 （就是 ）的數字，這是因為電腦本身對數字精確度的限制。 10. 將工作表更名為「聯立方程式」後，切換到「Sheet2」，將 「Sheet2」 更名為「重算」。 11. 「循環參照」是在 Excel 中相當致命的問題。為示範這個問題，在工作表「重算」上的儲存格 A1 填入 =A2+1（A1 此時顯示 1）；再於 A2 填入 =A1+1，Excel 會給你以下的訊息： 因為計算 A1 要用到 A2 的值，計算 A2 要用到 A1 的值，形成了循環參照。通常若非有特殊目的，這個警告通常表示公式可能有問題。不過