readmeformatlabnotebook.docVIP

第0章 計算的誤差 本章旨在探討當我們執行算數計算時, 近似值取代真實值所產生的誤差之影響,  及如何避免誤差之擴展. 本章中沒有用到Matlab 的m-file, 往後如果有需要, 則將必要的m-file之檔案夾 加入搜尋路徑中, 例如: path(d:\numerical, path) 註: 如果你有安裝Matlab Notebook 要執行下列input cells (綠色指令敘述) 之前必須先執行上面的cell – [path (…) ] 藍色的內容是Matlab [output cells] 本章同時利用Extended Symbolic ToolBox( Maple)的指令或函數 例題1: Find the roots of equation x2 – x – 1/3 = 0 .  首先我們利用函數solve得到真實解 rt = solve( x^2 - x - 1/3 ) rt = [ 1/2+1/6*21^(1/2)] [ 1/2-1/6*21^(1/2)] 再利用digits函數得到10個有效數位之近似根 , digits(10) x1 = vpa(rt(1)) %vpa()說明於後 x2 = vpa(rt(2)) x1 = 1.263762616 x2 = -.2637626160 註: 函數vpa是variable-precision arithmetic的簡稱, 用來配合固定 有效數位-- digits()之計算. 如果將1/3取代為0.3333, 則兩近似根為 roots( [ 1 -1 -0.3333] ) ans = 1.2637 -0.2637 兩近似根維持4位有效數位. 註: Matlab 內定16位有效數位, 顯示時一般格式只有5個數字, Maple 內定32位有效數位. 例題2: Find the roots of equation 1.002x2 - 11.01x + 0.01265 =0 . 其近似根為 : format long ; roots( [1.002 -11.01 0.01265] ) format short ans = 10.98687487643590 0.00114907565991 假設我們所有的每個計算都維持4位有效數位, 如同一部4-digit的計算機, 下面套用二次方程式根的公式: 觀察每個步驟之數據: a =1.002 ; b = -11.01; c = 0.01265 ; digits(4) a2 = vpa(2 *a) bb = vpa(b^2) ac4 = vpa(4*a*c) dif = vpa(bb - ac4) sqr = vpa(sqrt(dif)) x1 = vpa((-b+sqr) /a2) x2 = vpa((-b-sqr) /a2) a2 = 2.004 bb = 121.2 ac4 = .5070e-1 dif = 121.1 sqr = 11.00 x1 = 10.98 x2 = .4990e-2 上面數據顯示近似根 x1蠻好, 但是 x2誤差蠻大的 .因為-b-sqrt(b2-4ac)  兩個相近的數值相減刪除重要的有效數字, 稱為 loss of significant digits . 下面兩個方法改善誤差的影響. ( I ) 將x2改採另一公式計算 : c2 = vpa(2*c