數(shù)值分析1-誤差及有效數(shù)字課件

1、數(shù)值分析： 研究各類數(shù)學(xué)問(wèn)題求解的數(shù)值計(jì)算及相關(guān)理論分析。 隨著計(jì)算機(jī)的產(chǎn)生和發(fā)展，數(shù)值分析越來(lái)越多地研究如何借助于計(jì)算機(jī)求解相關(guān)問(wèn)題。計(jì)算方法： 隨著計(jì)算機(jī)產(chǎn)生和發(fā)展而建立的一個(gè)重要數(shù)學(xué)分支，是研究建立計(jì)算機(jī)解決各種數(shù)學(xué)問(wèn)題的數(shù)值計(jì)算及相關(guān)理論分析。第一章 緒論1.1數(shù)值分析(計(jì)算方法)介紹：(Numerical Analysis)(Computational Method)主要內(nèi)容：（1）數(shù)值計(jì)算：非線性方程求根，（非）線性方程組求解，插值，逼近（最小二乘擬合），數(shù)值微分（積分），常微分方程，矩陣特征值求解，偏微分方程數(shù)值解,（2）理論分析：誤差分析，計(jì)算過(guò)程的收斂性、穩(wěn)定性（數(shù)學(xué)角度上）

2、，算法的計(jì)算時(shí)間復(fù)雜度，存儲(chǔ)容量大?。ㄓ?jì)算機(jī)角度上） 特點(diǎn) :具有數(shù)學(xué)的抽象性和邏輯嚴(yán)密性又具有廣泛的應(yīng)用性和高度的技術(shù)性（與計(jì)算機(jī)結(jié)合密切的一門(mén)課程）使用計(jì)算機(jī)進(jìn)行數(shù)值問(wèn)題求解是主要研究對(duì)象。如何學(xué)習(xí)這門(mén)課？這門(mén)課的學(xué)習(xí)意義，數(shù)值計(jì)算的重要性;如何上這門(mén)課（教材）, 學(xué)習(xí)方法;上課形式（授課、上機(jī)、大型實(shí)驗(yàn)）;成績(jī)?cè)u(píng)定（平時(shí)、實(shí)驗(yàn)、期中、期末）.1.2誤差基本概念1.2.1誤差定義及來(lái)源真實(shí)值與觀察、測(cè)量或計(jì)算的值之間存在差異，其差稱為誤差。結(jié)合實(shí)際問(wèn)題求解，誤差來(lái)源可分為：(1). 模型誤差（實(shí)際問(wèn)題數(shù)學(xué)問(wèn)題）， 如抽象化、忽略次要因素等.(2). 觀測(cè)誤差（數(shù)學(xué)問(wèn)題中的數(shù)據(jù)初始值觀察 測(cè)

3、量時(shí)產(chǎn)生）(Error)(3). 截?cái)嗾`差（計(jì)算過(guò)程中存在的一些無(wú)限計(jì)算），如無(wú)窮級(jí)數(shù)求和（無(wú)限次有限次： ，(4). 舍入誤差（計(jì)算結(jié)果中存在數(shù)據(jù)無(wú)限位，如Pi，無(wú)理數(shù)有理數(shù)，） 整個(gè)誤差來(lái)源可做圖表示: 總結(jié)：誤差是不可避免的，應(yīng)盡量減少誤差，提高精度（如選擇好的計(jì)算方法）1.2.2絕對(duì)誤差和絕對(duì)誤差限 定義：設(shè) 為準(zhǔn)確值， 是近似值 , 為絕對(duì)誤差分析：e可正可負(fù)（并不因?yàn)槭墙^對(duì)誤差，就以為是正值)e值實(shí)際上無(wú)法知道, 不知道， 但能知道誤差的某個(gè)范圍（即誤差限） 例：毫米刻度的尺子，正常情況下誤差不超過(guò) 0.5mm. 定義：若 ，則 稱為絕對(duì)誤差限， 為正數(shù)，有:1.2.3相對(duì)誤差和相

4、對(duì)誤差限為什么引入？因?yàn)橛美迕卓潭鹊某咦訙y(cè)量1米長(zhǎng)和10米長(zhǎng)的物體，其絕對(duì)誤差限都為0.5，但測(cè)量精度分別為1/100和1/1000，所以為了較好反應(yīng)測(cè)量精確度，引入相對(duì)誤差。定義： 為準(zhǔn)確值， 為近似值，則分析： (1). 可正可負(fù) (2). (3). 無(wú)法知道，因?yàn)?不知道， 也可表示為 和 之間關(guān)系為： （可作為習(xí)題）因?yàn)?無(wú)法求出，所以通常考慮相對(duì)誤差限若 或則稱 為相對(duì)誤差限。1.2.4 有效數(shù)字當(dāng) 有很多位數(shù)表示時(shí)，可按四舍五入取前幾位。定義：如果近似值 的誤差限是其末位上的半個(gè)單位，且該位直到 的第一個(gè)非零數(shù)字共有n位，則 有n位有效數(shù)字。具體計(jì)算：對(duì) ，從左往右數(shù)，從第一個(gè)非零

5、數(shù)字開(kāi)始，直到最右面的數(shù)共有n個(gè)，且其誤差限為末位的 個(gè)單位，則有效數(shù)字為n。 有效數(shù)字的位數(shù)確定.例：數(shù)0.00234711，取五位有效數(shù)字，例： =1.732050808若 =1.7321，但若 =1.7320，誤差限為則有5位有效數(shù)字，因?yàn)檎`差限則只有4位有效數(shù)字，因?yàn)檎`差限為0.0023471，1.2.5誤差傳播影響計(jì)算過(guò)程中（如四則運(yùn)算）的初始數(shù)據(jù)誤差會(huì)導(dǎo)致函數(shù)值誤差.泰勒級(jí)數(shù)展開(kāi)分析誤差傳播.設(shè) 為準(zhǔn)確值，準(zhǔn)確值為為近似值，近似值為先考慮絕對(duì)誤差：令利用二元函數(shù)一階泰勒展開(kāi)公式采用二元函數(shù)所以：再考慮相對(duì)誤差：根據(jù)以上兩公式，可得到兩數(shù)相加、減、乘、除的誤差傳播： （避免絕對(duì)值很大

6、的數(shù)為乘數(shù)） （避免 為很小的數(shù)為除數(shù)） （避免兩相近數(shù)相減運(yùn)算)1.3 機(jī)器數(shù)系. (略.主要防止計(jì)算機(jī)處理過(guò)程中的數(shù)字溢出和含入誤差) 這里，主要介紹計(jì)算機(jī)中浮點(diǎn)數(shù)的表示形式及表示范圍（4個(gè)參數(shù)）：其中， =0.a1a2a3at 稱為尾數(shù)-1,1, 中的正負(fù)號(hào)用一位數(shù)字區(qū)分；為基數(shù)，如取2、10、8、16；p為階數(shù)，有上限U和下限L, 由計(jì)算機(jī)存儲(chǔ)字節(jié)長(zhǎng)度決定。 1.4 誤差危害的防止（1）使用數(shù)值穩(wěn)定的計(jì)算公式數(shù)值穩(wěn)定是指計(jì)算過(guò)程中舍入誤差對(duì)計(jì)算影響不大的算法，若第n+1步的誤差en+1 與第n步的誤差en滿足，則稱該計(jì)算公式是絕對(duì)穩(wěn)定的例：建立積分In= （n=0,1.,20）遞推關(guān)系

7、式，并分析誤差傳播影響。解： In+5In-1= I0=ln6-ln5 遞推式： 在計(jì)算I0時(shí)，設(shè)近似值為I0為 可設(shè) e0=I0-In- =即初始誤差對(duì)第n步的影響是擴(kuò)大5n倍，誤差范圍變大，不穩(wěn)定. 對(duì)可改用另一種計(jì)算過(guò)程：（ 可通過(guò)積分第一中值定理算出）則 ，誤差范圍逐步減少。即若函數(shù)f(x)連續(xù), g(x)在區(qū)間a,b上不變號(hào)且可積, 則有設(shè)（2）避免兩相近數(shù)相減 例. 計(jì)算設(shè)和有六位有效數(shù)字,即x1=44.7325 x2=44.7102x1-x2=44.7325-44.7102（可以根據(jù)需要取任意位有效數(shù)字,這里取6位）方法1：直接相減：方法2:分子有理化：=0.0223 （事實(shí)上只

8、有2位有效數(shù)字）也可進(jìn)行理論分析，這里考慮絕對(duì)誤差：第一種方法只有2位有效數(shù)字理論上分析, 可以有6位有效數(shù)字(分子為常數(shù)2,分母為x1+x2兩變量之和)（3）避免絕對(duì)值大的數(shù)作乘數(shù)， 同樣，避免x2為很小的數(shù)作除數(shù)， （4）防止大數(shù)吃小數(shù)： （計(jì)算機(jī)硬件發(fā)展，浮點(diǎn)數(shù)表示位數(shù)增加，此問(wèn)題已很少出現(xiàn)） 主要原因是計(jì)算機(jī)運(yùn)算處理時(shí)，需對(duì)階處理（即取較大的階值運(yùn)算,較小數(shù)的尾數(shù)則會(huì)變的很小,計(jì)算機(jī)浮點(diǎn)數(shù)表示不出來(lái)）, 會(huì)出現(xiàn): 大數(shù)+小數(shù)=大數(shù)求和時(shí)，可先按絕對(duì)值從小到大排序，先對(duì)小數(shù)運(yùn)算，再對(duì)大數(shù)運(yùn)算。(5)簡(jiǎn)化計(jì)算步驟，減少計(jì)算次數(shù)例：計(jì)算方法1：直接計(jì)算30次乘法方法2：（這里4次乘法）（4次

9、乘法）共8次乘法空間上:需存儲(chǔ)x,x2,x4,x8,x16, 方法1只需要存儲(chǔ)x.例：計(jì)算常規(guī)方法：乘法： 加法：Horner方法（秦九韶方法）: 需n次乘法，n次加法空間上:除了an和x, 多存儲(chǔ)一個(gè)變量用來(lái)保存ai-1x+ai第2章 方程求根(Non-linear equation)2.1問(wèn)題提出對(duì)方程 ，若存在 ，使得 ， 則稱 為 的根，或稱為零點(diǎn)。當(dāng)為多項(xiàng)式形式時(shí)，即則 稱為代數(shù)方程。若可寫(xiě)成形式,為的m重根，或稱m重零點(diǎn)。則代數(shù)方程 的公式解（當(dāng)次數(shù) 時(shí)有）令 ，原方程又可寫(xiě)為：對(duì)三次方程（卡當(dāng)公式）：此類方程有公式解： 其中,（有可能出現(xiàn)復(fù)數(shù)根）對(duì)四次方程，可找相關(guān)文獻(xiàn)。對(duì)高次方程，使用數(shù)值方法求解，即在滿足一定精度的前提下，求根的近似值。具體步驟：找到根的隔離區(qū)間當(dāng) 在 內(nèi)連續(xù)，且則 內(nèi)有解； 當(dāng) 在 內(nèi)嚴(yán)格單調(diào)，則 內(nèi)