淺議誤差理論在數值分析課程教學中的重要性

黃陳思，劉 蓉

（福州大學 數學與計算機科學學院，福建 福州 350002）

計算機與數學的有機結合形成了“科學計算”的研究方法，它的核心內容是以現代化的計算機以及數學軟件為工具，以數學模型為基礎進行模擬研究[1]。一個實際問題在計算機上的解決需要轉化成數學模型，當數學模型不能夠得到精確解時，通常要用數值方法求它的近似解，而后再編制成程序在計算機上運行。顯然，這些計算方法是數值分析課程研究的對象，只有構造出好的算法才能有效的在計算機上解決數學問題。然而，構造計算方法的手段是近似，研究并分析算法的核心問題則是近似對計算結果的影響。隨著授課經驗的積累以及對這門課程了解的深入，作者越來越覺得誤差的講解應該成為整門課程講授的重點和核心。實際問題通過建立數學模型在計算機上得到解決，每一步都會產生誤差。只有一個對誤差進行詳盡分析的算法才能稱之為一個好的算法。下面就誤差對算法的影響以及算法的誤差分析等幾個方面進行闡述，從中可以看出將誤差分析在課堂講授中應當特別作為重點是很有實際意義的。

一、算法的穩定性

用一個算法進行計算，由于初始數據誤差在計算中傳播使計算結果誤差增長很快，就是數值不穩定的[2]。在大部分數值計算教材中都能見到這樣一個關于描述算法穩定性的例子，式（1）與式（2）兩個公式是相互等價的，但是由于初始數據帶有誤差，計算結果就迥然不同。

當初始數據帶有誤差時，式（1）的計算過程使得這個誤差急劇擴散，速度為n！，假設初始數據I0的誤差為δ，則當n=7時I7的誤差為δ的5040倍。而式（2）卻能夠在計算的過程中逐步地縮小誤差，控制誤差的傳播。這就是算法的穩定性問題。在常微分方程初值問題的數值解法中也會發生這樣的問題。例如，分別用歐拉顯、隱式格式和改進的歐拉格式計算初值問題在區間[0，0.5]上的數值解。結果如圖1所示：

圖1 常微分方程初值問題

注意到在這個結果里，局部截斷誤差為2階的改進歐拉法在h=0.3，0.4，0.5的值相對于局部截斷誤差為1階的歐拉隱式公式在這些點的計算結果更加偏離了精確解。這是由于歐拉隱式公式的穩定性高于改進歐拉法。從這一方面也能夠看出算法的穩定性問題在實際計算中的重要性。理論上的截斷誤差并不能正確地反映一個算法計算結果的好壞。

二、病態問題

此時它的精確解為x=（1，1）T，從這個例子可以看到常數項的第2個分量只有萬分之一的微小變化，然而方程組的解卻變化很大。這樣的方程組就稱之為病態方程組。病態方程組的影響在《數值分析》的幾個章節中都是很有必要在課堂中加以解釋的。因此，必須要去尋求避免出現此類問題的解決方案。在這種情形下可以引導學生進行思考，既要保持原有的解法，又要避免求解線性方程組。解決的方法就是以正交的函數族代替φi（x），用正交多項式做最小二乘曲線擬合，系數矩陣變為對角矩陣，就可以避免利用線性方程組，而能夠直接計算出系數。同樣的，在插值法的推導以及數值微積分中，由于考慮到求解線性方程組可能會遇到病態方程組或者計算過程比較繁雜，都需要另外尋求解決方案，由此可見，清晰而明了地提醒學生時時需要注意到求解方程組的病態問題，才能夠盡可能地避免在實際求解中實際值偏離理論值的情況，這一點，無論是在算法的講解還是在實際問題的求解中都是不可避免的。

三、利用誤差創造新算法

誤差的分析既可以對算法的優劣進行判斷，還可以利用誤差來創造新的算法。例如，在數值微積分這一章節中的推導自適應的積分方法，復合梯形公式的余項公式中涉及到f（x）在η點的二階導數f"（η），在實際中由于f（x）復雜或未知，難以通過余項公式來估計誤差，而需要改變思路。也就是誤差的估計不能使用余項公式，但觀察余項公式后，可以引導學生考慮利用區間折半以后的復合梯形公式的余項中f"（η）與折半前復合梯形公式余項中的f"（ξ）近似相等，得到如下的式子

綜上所述，算法的誤差分析與課程每一個環節都息息相關，在講解算法的時候都很有必要提醒學生時刻注意到對誤差進行分析，如果脫離了誤差分析，而僅僅在于講授算法，就失去了算法應用于實際的意義。本文并沒有對算法的誤差分析方法進行討論，而是從誤差分析的另一個角度剖析了這個內容在課程講授中的重要性，引導學生在學習過程中要注意實際問題中的舍入誤差帶來的影響和理論誤差公式的應用。通過強調誤差分析，能夠在一定程度上提升學生對待問題的嚴謹性與探索性，使得數值分析課程的教與學更加豐富，更加貼近工程計算的實質。

[1]白峰杉.數值計算引論第2版[M].高等教育出版社，2010.

[2]李慶揚，等.數值分析第5版[M].清華大學出版社，2008.