數(shù)值分析第一章 緒論

數(shù)學(xué)與信息科學(xué)學(xué)院數(shù)值分析教材

(TextBook)

數(shù)值分析(第5版)李慶揚(yáng)等，清華大學(xué)出版社，2008年數(shù)值計(jì)算方法，李維國(guó)等，石油大學(xué)出版社，2009年參考書目

(Reference)

數(shù)值計(jì)算方法(上、下冊(cè))

林成森編著（科學(xué)出版社1998年）

Principleof

NumericalAnalysis

數(shù)值分析原理封建湖、車剛明、聶玉峰編著（科學(xué)出版社，2001年）

MATLAB數(shù)值計(jì)算方法

,張德豐編(機(jī)械工業(yè)出版,2010年)

課程評(píng)分方法

(GradingPolicies)平時(shí)成績(jī)(30%左右)

期末考試成績(jī)

(70%左右)本課程稱為“數(shù)值分析”,也常被稱為“數(shù)值計(jì)算方法”、“計(jì)算方法”等。它是一門研究數(shù)學(xué)問(wèn)題的數(shù)值數(shù)值解法及其相關(guān)理論的課程

提問(wèn)：數(shù)值分析是做什么用的？研究對(duì)象：數(shù)值問(wèn)題——有限個(gè)輸入數(shù)據(jù)（問(wèn)題的自變量、原始數(shù)據(jù)）與有限個(gè)輸出數(shù)據(jù)（待求解數(shù)據(jù)）之間函數(shù)關(guān)系的一個(gè)明確無(wú)歧義的描述。如一階微分方程初值問(wèn)題求函數(shù)解析表達(dá)式數(shù)學(xué)問(wèn)題求函數(shù)在某些點(diǎn)的近似函數(shù)值數(shù)值問(wèn)題程序設(shè)計(jì)上機(jī)計(jì)算設(shè)計(jì)高效、可靠的數(shù)值方法數(shù)值問(wèn)題求解近似結(jié)果輸出重點(diǎn)討論數(shù)值問(wèn)題的來(lái)源：實(shí)際問(wèn)題建立數(shù)學(xué)模型數(shù)值問(wèn)題數(shù)值方法的設(shè)計(jì)原則收斂性：方法的可行性穩(wěn)定性：初始數(shù)據(jù)等產(chǎn)生的誤差對(duì)結(jié)果的影響便于編程實(shí)現(xiàn)：邏輯復(fù)雜度要小計(jì)算量要小：時(shí)間復(fù)雜度要小，運(yùn)行時(shí)間要短存貯量要盡量?。嚎臻g復(fù)雜度要小可靠性分析計(jì)算復(fù)雜性誤差估計(jì)：運(yùn)算結(jié)果不能產(chǎn)生太大的偏差且能夠控制誤差§1

誤差

/*Error*/一、誤差的來(lái)源與分類

/*Source&Classification*/

1、從實(shí)際問(wèn)題中抽象出數(shù)學(xué)模型

——模型誤差

/*ModelingError*/

2、通過(guò)觀測(cè)得到模型中某些參數(shù)（或物理量）的值

——觀測(cè)誤差

/*MeasurementError*/

3、數(shù)學(xué)模型與數(shù)值算法之間的誤差求近似解

——方法誤差

(截?cái)嗾`差

/*TruncationError*/)

4、由于機(jī)器字長(zhǎng)有限，原始數(shù)據(jù)和計(jì)算過(guò)程會(huì)產(chǎn)生新的誤差

——舍入誤差

/*RoundoffError*/二、誤差分析的基本概念

/*BasicConcepts*/設(shè)為真值（精確值），為的一個(gè)近似值稱為近似值的絕對(duì)誤差，簡(jiǎn)稱誤差。注：誤差可正可負(fù)，常常是無(wú)限位的絕對(duì)誤差限/*accuracy*/——絕對(duì)值的上界如：絕對(duì)誤差還不能完全表示近似值的好壞（絕對(duì)誤差/*absoluteerror*/）近似值的誤差與準(zhǔn)確值的比值：稱為近似值的相對(duì)誤差，記作注：實(shí)際計(jì)算時(shí)，相對(duì)誤差通常取因?yàn)椋ㄏ鄬?duì)誤差/*relativeerror*/）相對(duì)誤差也可正可負(fù)（有效數(shù)字/*SignificantDigits*/）相對(duì)誤差限——相對(duì)誤差的絕對(duì)值的上界/*relativeaccuracy*/如：3位6位若近似值與準(zhǔn)確值的誤差絕對(duì)值不超過(guò)某一位的半個(gè)單位，該位到的第一位非零數(shù)字共有位，則稱有位有效數(shù)字有效數(shù)字(另外一種定義形式)注：0.2300有4位有效數(shù)字，而00023只有2位有效數(shù)字。12300如果寫成0.123105，則表示只有3位有效數(shù)字。數(shù)字末尾的0不可隨意省去！用科學(xué)計(jì)數(shù)法，記其中，若（即的截取按四舍五入規(guī)則），則稱為有位有效數(shù)字，精確到。例1：?jiǎn)枺河袔孜挥行?shù)字？請(qǐng)證明你的結(jié)論。證明：有位有效數(shù)字，精確到小數(shù)點(diǎn)后第位。43規(guī)格化形式有效數(shù)字與相對(duì)誤差限的關(guān)系.

有效數(shù)字相對(duì)誤差限已知x*

有n位有效數(shù)字，則其相對(duì)誤差限為相對(duì)誤差限

有效數(shù)字已知x*

的相對(duì)誤差限可寫為則可見(jiàn)x*

至少有n

位有效數(shù)字。例2：為使

的相對(duì)誤差小于0.001%,至少應(yīng)取幾位有效數(shù)字？解：假設(shè)*取到n位有效數(shù)字，則其相對(duì)誤差上限為要保證其相對(duì)誤差小于0.001%，只要保證其上限滿足已知

a1=3，則從以上不等式可解得

n>6log6，即

n6，應(yīng)取

*=3.14159。例3計(jì)算下列多項(xiàng)式的值為已知數(shù)據(jù)分析：輸入數(shù)據(jù)為，輸出數(shù)據(jù)為，若直接由算出，再乘相應(yīng)的系數(shù)并相加，則要做次乘法和次加法，占用個(gè)存儲(chǔ)單元。秦九韶方法，也稱為Horner算法用遞推公式表示為只用次乘法和次加法，并占用個(gè)存儲(chǔ)單元

三、數(shù)值算法及穩(wěn)定性

/*NumericalAlgorithmandStability*/大家一起猜？11/e解法之一：將作Taylor展開(kāi)后再積分S4R4

/*Remainder*/取則稱為截?cái)嗾`差

/*TruncationError*/|

舍入誤差

/*RoundoffError*/|=0.747……由截去部分/*excludedterms*/引起由留下部分/*includedterms*/引起例4近似計(jì)算一個(gè)算法如果輸入數(shù)據(jù)有擾動(dòng)（即誤差），而計(jì)算過(guò)程中舍入誤差不增長(zhǎng),則稱此算法是數(shù)值穩(wěn)定的，否則此算法就稱為不穩(wěn)定的。（數(shù)值穩(wěn)定性/*NumericalStability*/）對(duì)數(shù)學(xué)問(wèn)題本身如果輸入數(shù)據(jù)有微小擾動(dòng)，引起輸出數(shù)據(jù)（即問(wèn)題真解）的很大擾動(dòng)，這就是病態(tài)問(wèn)題。（病態(tài)問(wèn)題/*ill-posedproblem*/）

它是數(shù)學(xué)問(wèn)題本身性質(zhì)所決定的，與算法無(wú)關(guān)，也就是說(shuō)對(duì)病態(tài)問(wèn)題，用任何算法（或方法）直接計(jì)算都將產(chǎn)生不穩(wěn)定性。

此公式精確成立記為則初始誤差????!!!Whathappened?!例5計(jì)算公式一：考察第n步的誤差我們有責(zé)任改變。造成這種情況的是不穩(wěn)定的算法

/*unstablealgorithm*/迅速積累，誤差呈遞增趨勢(shì)。初始的小擾動(dòng)公式二：注意此公式與公式一在理論上等價(jià)。方法：先估計(jì)一個(gè)IN

,再反推要求的In(n<<N)?？扇∪】疾旆赐埔徊降恼`差：以此類推，對(duì)n<N

有：誤差逐步遞減,這樣的算法稱為穩(wěn)定的算法

/*stablealgorithm*/

在我們今后的討論中，誤差將不可回避，算法的穩(wěn)定性將會(huì)是一個(gè)非常重要的話題。例6：蝴蝶效應(yīng)

——紐約的一只蝴蝶翅膀一拍，風(fēng)和日麗的北京就刮起臺(tái)風(fēng)來(lái)了？！紐約北京這是一個(gè)病態(tài)問(wèn)題關(guān)于蝴蝶效應(yīng)§2

誤差分析的方法和原則

/*ErrorAnalysis*/一、誤差分析的方法利用原始數(shù)據(jù)的誤差限，根據(jù)誤差的運(yùn)算法則，隨著計(jì)算過(guò)程逐步向前進(jìn)行分析，直至估計(jì)出最后的結(jié)果。注：兩個(gè)近似數(shù)，四則運(yùn)算得到的誤差限分別為（1）（2）對(duì)于函數(shù)y=f(x)，若用x*取代x，將對(duì)y產(chǎn)生什么影響？分析：e*(y)=f(x*)f(x)e*(x)=x*xMeanValueTheoremx*

與x非常接近時(shí)，可認(rèn)為，則有：即：產(chǎn)生的誤差經(jīng)過(guò)作用后被放大/縮小了

倍。故稱

為放大/縮小因子

/*amplificationfactor*/

或

絕對(duì)條件數(shù)

/*absoluteconditionnumber*/.相對(duì)誤差條件數(shù)

/*relativeconditionnumber*/

f的條件數(shù)在某一點(diǎn)是小\大，則稱

f在該點(diǎn)是好條件的

/*well-conditioned*/

\壞條件的

/*ill-conditioned*/。注：關(guān)于多元函數(shù)

可類似討論，理論工具：Taylor公式例7設(shè),試求函數(shù)的相對(duì)誤差限.解：由題設(shè)知:近似值為,絕對(duì)誤差限為

二、幾點(diǎn)注意事項(xiàng)

/*Remarks*/1、

避免相近二數(shù)相減例：a1=0.12345，a2=0.12346，各有5位有效數(shù)字。而a2

a1=0.00001，只剩下1位有效數(shù)字。

幾種經(jīng)驗(yàn)性避免方法：當(dāng)|x|<<1時(shí)：2、

避免小分母:分母小會(huì)造成浮點(diǎn)溢出

/*overflow*/3、避免大數(shù)吃小數(shù)例：用單精度計(jì)算的根。精確解為算法1：利用求根公式在計(jì)算機(jī)內(nèi)，109存為0.11010，1存為0.1101。做加法時(shí)，兩加數(shù)的指數(shù)先向大指數(shù)對(duì)齊，再將浮點(diǎn)部分相加。即1的指數(shù)部分須變?yōu)?010，則：1=0.00000000011010，取單精度時(shí)就成為：

109+1=0.100000001010+0.000000001010=0.100000001010大數(shù)吃小數(shù)算法2：先解出再利用注：求和時(shí)從小到大相加，可使和的誤差減小。例：按從小到大、以及從大到小的順序分別計(jì)算1+2+3+…+40+1094、

先化簡(jiǎn)再計(jì)算，減少步驟，避免誤差積累。一般來(lái)說(shuō)，計(jì)算機(jī)處理下列運(yùn)算的速度為5、選用穩(wěn)定的算法。稱為尾數(shù)，j為階§3計(jì)算機(jī)的數(shù)系結(jié)構(gòu)/*StructureofNumberSystem

*/計(jì)算機(jī)的數(shù)系是一個(gè)不完整的數(shù)系。計(jì)算機(jī)只能表示有限個(gè)數(shù)，即計(jì)算機(jī)的精度是有限的。每種計(jì)算機(jī)內(nèi)部運(yùn)算是按固定的有限位數(shù)進(jìn)行的，也就是按固定位數(shù)的有限位浮點(diǎn)數(shù)進(jìn)行運(yùn)算的。浮點(diǎn)數(shù)系統(tǒng)由四個(gè)整數(shù)表征：基，精度（尾數(shù)）位數(shù)，下溢界和上溢界計(jì)算機(jī)所表示的數(shù)的集合(規(guī)格化浮點(diǎn)數(shù)系)：

若