數(shù)值分析第一章誤差

1、1數(shù)數(shù) 值值 分分 析析林甲富林甲富2教材教材丁麗娟丁麗娟, 程杞元程杞元,數(shù)值計(jì)算方法數(shù)值計(jì)算方法, 高等教育高等教育出版社出版社, 2011年年.3最后成績(jī)最后成績(jī)= =實(shí)驗(yàn)作業(yè)成績(jī)實(shí)驗(yàn)作業(yè)成績(jī)(20%)+(20%)+考試成績(jī)考試成績(jī)(80%)(80%)實(shí)驗(yàn)作業(yè)：下列實(shí)驗(yàn)作業(yè)：下列1 1和和2 2選擇一個(gè)選擇一個(gè), ,希望選希望選2 2 1. 1.課堂布置的課本上的數(shù)值實(shí)驗(yàn)題課堂布置的課本上的數(shù)值實(shí)驗(yàn)題 2.2.結(jié)合所學(xué)專業(yè)自選題結(jié)合所學(xué)專業(yè)自選題 (1)(1)敘述實(shí)際問題敘述實(shí)際問題(2)(2)建立數(shù)學(xué)模型建立數(shù)學(xué)模型( (解常微分方程組，數(shù)據(jù)擬合等解常微分方程組，數(shù)據(jù)擬合等) )(3)

2、(3)設(shè)計(jì)計(jì)算方法設(shè)計(jì)計(jì)算方法(4)(4)程序程序(matlab(matlab) )(5)(5)計(jì)算結(jié)果及分析計(jì)算結(jié)果及分析 3.3.交打印文件交打印文件( (截止至考試前一日交截止至考試前一日交) ) 4. 4.注意完全重復(fù)的實(shí)驗(yàn)作業(yè)沒有實(shí)驗(yàn)作業(yè)成績(jī)注意完全重復(fù)的實(shí)驗(yàn)作業(yè)沒有實(shí)驗(yàn)作業(yè)成績(jī)4 數(shù)值分析是做什么用的？數(shù)值分析是做什么用的？數(shù)值數(shù)值分析分析輸入復(fù)雜問題或運(yùn)算輸入復(fù)雜問題或運(yùn)算.),(,)(,ln,xfdxddxxfbxAxaxbax 計(jì)算機(jī)計(jì)算機(jī)近似解近似解5 研究對(duì)象研究對(duì)象 那些在理論上有解而又無法手工計(jì)算的那些在理論上有解而又無法手工計(jì)算的數(shù)學(xué)問題數(shù)學(xué)問題 例例 解解300階

3、的線性方程組階的線性方程組 求求6階矩陣的全部特征值階矩陣的全部特征值6主要內(nèi)容主要內(nèi)容 數(shù)值代數(shù)數(shù)值代數(shù)近似求解線性方程組近似求解線性方程組 (直接解法直接解法, 迭代解法迭代解法)矩陣特征值的計(jì)算矩陣特征值的計(jì)算 數(shù)值逼近：數(shù)值逼近： 插值法，函數(shù)逼近插值法，函數(shù)逼近 數(shù)值微分與數(shù)值積分?jǐn)?shù)值微分與數(shù)值積分 微分方程近似求解：微分方程近似求解：常微分方程數(shù)值解法常微分方程數(shù)值解法 非線性方程求解非線性方程求解 7第一章第一章 誤差誤差2 誤差的基本概念誤差的基本概念 3 數(shù)值計(jì)算中誤差的傳播數(shù)值計(jì)算中誤差的傳播4 數(shù)值計(jì)算中應(yīng)注意的問題數(shù)值計(jì)算中應(yīng)注意的問題 82 誤差的基本概念誤差的基本概

4、念 誤差按來源可分為：誤差按來源可分為： 模型誤差模型誤差 觀測(cè)誤差觀測(cè)誤差 截?cái)嗾`差截?cái)嗾`差 舍入誤差舍入誤差 誤差：精確解與近似解之間的差誤差：精確解與近似解之間的差9 模型誤差模型誤差 數(shù)學(xué)模型通常是由實(shí)際問題抽象得到數(shù)學(xué)模型通常是由實(shí)際問題抽象得到的，一般帶有誤差，這種誤差稱為的，一般帶有誤差，這種誤差稱為模型誤差模型誤差. 觀測(cè)誤差觀測(cè)誤差 數(shù)學(xué)模型中包含的一些參數(shù)通常是通數(shù)學(xué)模型中包含的一些參數(shù)通常是通過觀測(cè)和實(shí)驗(yàn)得到的，難免帶有誤差，這種誤差稱過觀測(cè)和實(shí)驗(yàn)得到的，難免帶有誤差，這種誤差稱為為觀測(cè)誤差觀測(cè)誤差. 截?cái)嗾`差截?cái)嗾`差 求解數(shù)學(xué)模型所用的數(shù)值方法通常求解數(shù)學(xué)模型所用的數(shù)值

5、方法通常是一種近似方法，這種因方法產(chǎn)生的誤差稱為是一種近似方法，這種因方法產(chǎn)生的誤差稱為截截?cái)嗾`差斷誤差或或方法誤差方法誤差.10 543251413121)1ln(xxxxxx實(shí)際計(jì)算時(shí)只能截取有限項(xiàng)代數(shù)和計(jì)算，如取前實(shí)際計(jì)算時(shí)只能截取有限項(xiàng)代數(shù)和計(jì)算，如取前5項(xiàng)有：項(xiàng)有：5141312112ln 這里產(chǎn)生誤差這里產(chǎn)生誤差 (記作記作R5 )截?cái)嗾`差截?cái)嗾`差 8171615R例如例如，利用，利用 ln(x+1) 的的Taylor公式計(jì)算公式計(jì)算 ln2,11 舍入誤差舍入誤差 由于計(jì)算機(jī)只能對(duì)有限位數(shù)進(jìn)行由于計(jì)算機(jī)只能對(duì)有限位數(shù)進(jìn)行, e原則保留有限位，這時(shí)產(chǎn)生的誤差稱為原則保留有限位，這時(shí)

6、產(chǎn)生的誤差稱為舍入誤舍入誤差差。, 231等都要按舍入等都要按舍入運(yùn)算，在運(yùn)算中像運(yùn)算，在運(yùn)算中像在數(shù)值分析中，均假定數(shù)學(xué)模型是準(zhǔn)確的，因而在數(shù)值分析中，均假定數(shù)學(xué)模型是準(zhǔn)確的，因而不考慮模型誤差和觀測(cè)誤差，只討論不考慮模型誤差和觀測(cè)誤差，只討論截?cái)嗾`差截?cái)嗾`差和和舍入誤差舍入誤差對(duì)計(jì)算結(jié)果的影響對(duì)計(jì)算結(jié)果的影響.12 設(shè)設(shè)x* 是準(zhǔn)確值是準(zhǔn)確值x 的一個(gè)近似值，記的一個(gè)近似值，記e=x x*稱稱 e為近似值為近似值 x* 的的絕對(duì)誤差絕對(duì)誤差，簡(jiǎn)稱誤差，簡(jiǎn)稱誤差.絕對(duì)誤差一般很難準(zhǔn)確計(jì)算絕對(duì)誤差一般很難準(zhǔn)確計(jì)算, 但可以估計(jì)上界但可以估計(jì)上界. 絕對(duì)誤差絕對(duì)誤差則稱則稱 為近似值為近似值 x

7、* 的的絕對(duì)誤差限絕對(duì)誤差限，簡(jiǎn)稱誤差限，簡(jiǎn)稱誤差限. 若若 滿足滿足 |e 絕對(duì)誤差、相對(duì)誤差和有效數(shù)字絕對(duì)誤差、相對(duì)誤差和有效數(shù)字13例例 用毫米刻度的米尺測(cè)量一長(zhǎng)度用毫米刻度的米尺測(cè)量一長(zhǎng)度 x, 如讀出的長(zhǎng)度如讀出的長(zhǎng)度是是 x*=765 mm, 由于誤差限是由于誤差限是 0.5 mm, 故準(zhǔn)確值故準(zhǔn)確值.mm5 .765,mm5 .764 x 精確值精確值x , 近似值近似值 x* 和誤差限和誤差限 之間滿足：之間滿足：通常記為通常記為 *xxx *xx 14 絕對(duì)誤差有時(shí)并不能完全地反映近似值的好壞，絕對(duì)誤差有時(shí)并不能完全地反映近似值的好壞，如測(cè)量如測(cè)量 100 m 和和 10 m

8、 兩個(gè)長(zhǎng)度，若它們的絕對(duì)誤兩個(gè)長(zhǎng)度，若它們的絕對(duì)誤差都是差都是 1 cm，顯然前者的測(cè)量結(jié)果比后者的準(zhǔn)確，顯然前者的測(cè)量結(jié)果比后者的準(zhǔn)確. 因此，決定一個(gè)量的近似值的精確度，除了因此，決定一個(gè)量的近似值的精確度，除了要看要看絕對(duì)誤差絕對(duì)誤差外，還必須考慮外，還必須考慮該量本身的大小該量本身的大小.15稱稱 er 為近似值為近似值 x* 的的相對(duì)誤差相對(duì)誤差. 記記,*xxxxeer 由于由于 x 未知，實(shí)際使用時(shí)總是將未知，實(shí)際使用時(shí)總是將 x* 的相對(duì)誤差取為的相對(duì)誤差取為*xxxxeer .|rre 相對(duì)誤差相對(duì)誤差 稱為近似值稱為近似值x*的的相對(duì)誤差限相對(duì)誤差限. |*| xr 16例

9、例 設(shè)設(shè) x*=1.24是由精確值是由精確值 x 經(jīng)過四舍五入得到的經(jīng)過四舍五入得到的近似值近似值, 求求x*的絕對(duì)誤差限和相對(duì)誤差限的絕對(duì)誤差限和相對(duì)誤差限.由已知可得由已知可得:所以所以 =0.005,245. 1235. 1 x%.4 . 024. 1005. 0 r 解解 一般地一般地, 凡是由準(zhǔn)確值經(jīng)過四舍五入得到的近似凡是由準(zhǔn)確值經(jīng)過四舍五入得到的近似值值, 其絕對(duì)誤差限等于該近似值末位的半個(gè)單位其絕對(duì)誤差限等于該近似值末位的半個(gè)單位.17有有 位有效數(shù)字，精確到小數(shù)點(diǎn)后第位有效數(shù)字，精確到小數(shù)點(diǎn)后第 位位* 有效數(shù)字有效數(shù)字 若近似值若近似值 x*滿足滿足 則稱則稱 x*準(zhǔn)準(zhǔn)確到

10、小數(shù)點(diǎn)后第確到小數(shù)點(diǎn)后第n位位. 并把從第一個(gè)非零數(shù)字到這并把從第一個(gè)非零數(shù)字到這一位的所有數(shù)字均稱為一位的所有數(shù)字均稱為有效數(shù)字有效數(shù)字.,1021|*|nxx 1415.3*.;8979321415926535.3 例例：?jiǎn)枺簡(jiǎn)枺?有幾位有效數(shù)字？有幾位有效數(shù)字？* 31050* .|解：解：4318例例 已知下列近似值的絕對(duì)誤差限都是已知下列近似值的絕對(duì)誤差限都是0.005, 問問它們它們具有幾位有效數(shù)字具有幾位有效數(shù)字? a=12.175, b=0.10, c=0.1, d=0.0032由于由于0.0050.5102,解解所以所以a 有有4位有效數(shù)字位有效數(shù)字1, 2, 1,7; b

11、有有2位有效數(shù)字位有效數(shù)字1, 0;c 有有1位有效數(shù)字位有效數(shù)字1;d 沒有有效數(shù)字沒有有效數(shù)字.19數(shù)數(shù)x*總可以寫成如下形式總可以寫成如下形式.10. 0*21mnaaax x* 作為作為x的近似值的近似值, 具有具有n位有效數(shù)字當(dāng)且僅當(dāng)位有效數(shù)字當(dāng)且僅當(dāng)nmxx 1021*其中其中m是整數(shù)是整數(shù), ai是是0到到9中的一個(gè)數(shù)字中的一個(gè)數(shù)字,. 01 a由此可見由此可見, 近似值的有效數(shù)字越多近似值的有效數(shù)字越多, 其絕對(duì)誤差越小其絕對(duì)誤差越小. 有效數(shù)字的另一等價(jià)定義有效數(shù)字的另一等價(jià)定義20故取故取 n=6，即取，即取 6 位有效數(shù)字位有效數(shù)字. 此時(shí)此時(shí) x*=1.41421.解解

12、則近似值則近似值x*可寫為可寫為由于由于 ,414. 12 ,10. 0*121 naaax. 011 a51101021*2 nx令令例例 為了使為了使 的近似值的絕對(duì)誤差不大于的近似值的絕對(duì)誤差不大于105，問應(yīng)取幾位有效數(shù)字？問應(yīng)取幾位有效數(shù)字？2 x21 相對(duì)誤差限與有效數(shù)字之間的關(guān)系相對(duì)誤差限與有效數(shù)字之間的關(guān)系.111211021.021010.01050 nnmnnmra.aaa.a.x* 有效數(shù)字有效數(shù)字 相對(duì)誤差限相對(duì)誤差限已知已知 x* = 0.a1a2an10m有有 n 位位有效數(shù)字有效數(shù)字，則其，則其相對(duì)誤差限相對(duì)誤差限為為22nmmnmnr.aaa.aaxxx 105

13、010)1()1(21010.0)1(210|*|*|11112111 相對(duì)誤差限相對(duì)誤差限 有效數(shù)字有效數(shù)字1110)1(21 nra已知已知 x* 的的相對(duì)誤差限相對(duì)誤差限可寫為可寫為則則可見可見 x* 至少有至少有 n 位有效數(shù)字位有效數(shù)字.23 基本運(yùn)算中基本運(yùn)算中( )的誤差估計(jì)的誤差估計(jì),105 . 0|414. 12|3 ,105 . 0|236. 25|3 問問?|414. 1236. 225| ?236. 2414. 152 3 數(shù)值計(jì)算中誤差的傳播數(shù)值計(jì)算中誤差的傳播如如24例例 計(jì)算計(jì)算 A=f (x1, x2). 如果如果x1, x2的近似值為的近似值為 x1*, x2

14、*, 則則A的近似值為的近似值為 A*=f (x1*, x2*), 用多元函數(shù)微分近用多元函數(shù)微分近似公式可以得到似公式可以得到*)(*)*,(*)(*)*,(*)(*)*,(*)(*)*,(*)*,(),(*)(2221112122221111212121xexxxfxexxxfxxxxxfxxxxxfxxfxxfAAAe 絕對(duì)誤差絕對(duì)誤差 e 運(yùn)算可近似看成微分運(yùn)算運(yùn)算可近似看成微分運(yùn)算.25由此可以得到基本運(yùn)算中由此可以得到基本運(yùn)算中( )的誤差估計(jì)的誤差估計(jì),),()()(2121xexexxe 和差的誤差限不超過各數(shù)的誤差限之和和差的誤差限不超過各數(shù)的誤差限之和.| )(| )(|

15、)(|2121xexexxe 26)()()()()(212121211221xexexxxexxxxexxxerrr ),()()(211221xexxexxxe | )(| )(| )(|2121xexexxerrr 乘法相對(duì)誤差限不超過各數(shù)相對(duì)誤差限之和乘法相對(duì)誤差限不超過各數(shù)相對(duì)誤差限之和.27,)()(22211221xxexxexxxe ).()()()(211222211221xexexxxxexxexxxerrr 乘除相對(duì)誤差限不超過各數(shù)相對(duì)誤差限之和乘除相對(duì)誤差限不超過各數(shù)相對(duì)誤差限之和. | )(| )(|2121xexexxerrr 28例例 設(shè)設(shè) y=xn, 求求 y

16、的相對(duì)誤差與的相對(duì)誤差與 x 的相對(duì)誤差之間的相對(duì)誤差之間的關(guān)系的關(guān)系.解解)()()(1xenxxeyenn )()()()()(1xnexxenxxenxyyeyernnr 所以所以xn 的相對(duì)誤差是的相對(duì)誤差是 x 的相對(duì)誤差的的相對(duì)誤差的n倍倍.x2的相對(duì)誤差是的相對(duì)誤差是 x 的相對(duì)誤差的的相對(duì)誤差的 2 倍倍,x的相對(duì)誤差是的相對(duì)誤差是 x 的相對(duì)誤差的的相對(duì)誤差的 1/2 倍倍.29 算法的數(shù)值穩(wěn)定性算法的數(shù)值穩(wěn)定性 一種數(shù)值算法一種數(shù)值算法, 如果其計(jì)算舍入誤差積累是可控如果其計(jì)算舍入誤差積累是可控制的制的, 則稱其為數(shù)值穩(wěn)定的則稱其為數(shù)值穩(wěn)定的, 反之稱為數(shù)值不穩(wěn)定的反之稱為

17、數(shù)值不穩(wěn)定的.30 101dxexIxnn利用分部積分法可得計(jì)算利用分部積分法可得計(jì)算In的遞推公式的遞推公式, 2 , 11, 1 nnIInn例例 計(jì)算積分計(jì)算積分算法算法1： 1010dxeIx, 2 , 11, 1 nnIInn6321. 0632120558. 011 e由此遞推計(jì)算由此遞推計(jì)算 I1, I2, , I9.解解311 , 2 , 8 , 9, )1(11 nInInn 10109919110110dxxIdxexe取近似值取近似值,0684. 0)10110(2119 eI由此計(jì)算由此計(jì)算 I8, I7, , I0.并將計(jì)算公式改寫為并將計(jì)算公式改寫為算法算法2：此時(shí)

18、此時(shí)10121|1*99 eII.0316. 0 32InI0I1I2I3I4I5I6I7I8I9算法算法10.63210.36790.26420.20740.17040.14800.11200.21600.72807.5520算法算法20.63210.36790.26420.20730.17090.14550.12680.11210.10350.0684真值真值0.63210.36790.26420.20730.17090.14550.12680.11240.10090.091633 對(duì)任何對(duì)任何 n都應(yīng)有都應(yīng)有In0, 但算法但算法1的計(jì)算結(jié)果顯示的計(jì)算結(jié)果顯示I8 (n+1)n!當(dāng)當(dāng)n=25時(shí)時(shí), 在每秒百億次乘除運(yùn)算計(jì)算機(jī)上求解時(shí)間為在每秒百億次乘除運(yùn)算計(jì)算機(jī)上求解時(shí)間為 首先首先, 若算法計(jì)算量太大若算法計(jì)算量太大, 實(shí)際計(jì)算無法完成實(shí)際計(jì)算無法完成(億年億年)13 42 其次，即使是可行算法，則計(jì)算量越大積累的誤其次，即使是可行算法，則計(jì)算量越大積累的誤差也越大差也越大. 因此，算法的計(jì)算量越小越好因此，算法的計(jì)算量越小越好.0111.)(axaxaxaxpnnnnn 若直接逐項(xiàng)計(jì)算，大約需要乘法運(yùn)算次數(shù)為若直接逐項(xiàng)計(jì)算，大約需要乘法運(yùn)算次數(shù)為2)1(12.)1( nnnn例例 計(jì)算計(jì)算n次多項(xiàng)式：次多項(xiàng)式：43一般地