數(shù)值分析課件第一章引論

1、數(shù)值分析數(shù)值分析任課教師：文有為電話： Email: 一、什么是數(shù)值分析一、什么是數(shù)值分析 數(shù)值分析是計(jì)算數(shù)學(xué)的一個(gè)主要部分?jǐn)?shù)值分析是計(jì)算數(shù)學(xué)的一個(gè)主要部分, ,計(jì)算數(shù)學(xué)是數(shù)計(jì)算數(shù)學(xué)是數(shù)學(xué)科學(xué)的一個(gè)分支學(xué)科學(xué)的一個(gè)分支, ,它研究用計(jì)算機(jī)求解各種數(shù)學(xué)問題它研究用計(jì)算機(jī)求解各種數(shù)學(xué)問題的數(shù)值計(jì)算方法及其理論與軟件實(shí)現(xiàn)的數(shù)值計(jì)算方法及其理論與軟件實(shí)現(xiàn). .實(shí)際問題實(shí)際問題數(shù)學(xué)模型數(shù)學(xué)模型數(shù)值計(jì)算方法數(shù)值計(jì)算方法 程序設(shè)計(jì)程序設(shè)計(jì)上機(jī)計(jì)算求出結(jié)果上機(jī)計(jì)算求出結(jié)果第第1 1章章 緒論緒論1 數(shù)值分析的研究對象與特點(diǎn)數(shù)值分析的研究對象與特點(diǎn)二、數(shù)值分析的基本內(nèi)容二、數(shù)值分析的基本內(nèi)容1 1、數(shù)值逼近、數(shù)

2、值逼近 插值法插值法 函數(shù)逼近與曲線擬和函數(shù)逼近與曲線擬和 數(shù)值積分與數(shù)值微分?jǐn)?shù)值積分與數(shù)值微分2 2、數(shù)值代數(shù)、數(shù)值代數(shù) 線性代數(shù)問題線性代數(shù)問題( (方程組和特征值方程組和特征值) ) 非線性方程非線性方程( (組組) )數(shù)值解法數(shù)值解法 3 3、常微方程數(shù)值解法和偏微方程數(shù)值解法、常微方程數(shù)值解法和偏微方程數(shù)值解法三、數(shù)值分析的特點(diǎn)三、數(shù)值分析的特點(diǎn),1) 1(312112ln1nn).31713151313131(22ln753vs,32)1ln(32xxxx1 1、面向計(jì)算機(jī)、面向計(jì)算機(jī) 2 2、可靠的理論分析、可靠的理論分析, ,保證收斂性、穩(wěn)定性保證收斂性、穩(wěn)定性3 3、良好的計(jì)

3、算復(fù)雜性、良好的計(jì)算復(fù)雜性4 4、數(shù)值實(shí)驗(yàn)、數(shù)值實(shí)驗(yàn)Cramer法則 vs Gauss消去法.31ln12()313xxxxx取四、如何學(xué)好數(shù)值分析四、如何學(xué)好數(shù)值分析1 1、注意掌握基本原理、處理技巧，誤差分析、注意掌握基本原理、處理技巧，誤差分析 3 3、積極動(dòng)手上機(jī)實(shí)踐、積極動(dòng)手上機(jī)實(shí)踐2 2、注重實(shí)際問題，練習(xí)、作業(yè)、注重實(shí)際問題，練習(xí)、作業(yè) 五、教學(xué)參考書五、教學(xué)參考書 數(shù)值計(jì)算引論 白峰杉 高等教育出版社 科學(xué)和工程計(jì)算基礎(chǔ) 施妙根等 清華大學(xué)出版社 數(shù)值分析，易大義等編，浙江科學(xué)技術(shù)出版社 數(shù)值方法教程，劉欽圣等編 , 冶金出版社,1998 計(jì)算方法，秦林祥等編, 兵器工業(yè)出版社

4、,1992 數(shù)值分析基礎(chǔ)，關(guān)治等編 , 高教出版社,1998一、誤差來源、分類一、誤差來源、分類 觀測誤差觀測誤差截?cái)嗾`差或方法誤差截?cái)嗾`差或方法誤差模型誤差模型誤差 2 數(shù)值計(jì)算的誤差數(shù)值計(jì)算的誤差截?cái)嗾`差：nnnxnfxfxffxPxf!)0(! 2)0(! 1)0()0()()()(2 1) 1()!1()()(nnnxnfxR舍入誤差舍入誤差.0000026. 014159. 3R數(shù)制轉(zhuǎn)換、機(jī)器數(shù). 在用數(shù)值方法解題過程中可能產(chǎn)生的誤差歸納起來有如下幾類： 1. 模型誤差 2. 觀測誤差 3. 截?cái)嗾`差 4. 舍入誤差誤差誤差 用數(shù)學(xué)方法解決一個(gè)具體的實(shí)際問題，首先要建立數(shù)學(xué)模型，這就

5、要對實(shí)際問題進(jìn)行抽象、簡化，因而數(shù)學(xué)模型本身總含有誤差，這種誤差叫做模型誤差 數(shù)學(xué)模型是指那些利用數(shù)學(xué)語言模擬現(xiàn)實(shí)而建立起來的有關(guān)量的描述 數(shù)學(xué)模型的準(zhǔn)確解與實(shí)際問題的真解不同實(shí)際問題的實(shí)際問題的真解真解數(shù)學(xué)模型的數(shù)學(xué)模型的真解真解為減化模型忽略次要為減化模型忽略次要因素因素定理在特定條件下建立與實(shí)定理在特定條件下建立與實(shí)際條件有別際條件有別模型誤差模型誤差 在數(shù)學(xué)模型中通常包含各種各樣的參變量，如溫度、長度、電壓等，這些參數(shù)往往是通過觀測得到的，因此也帶來了誤差，這種誤差叫觀測誤差 數(shù)學(xué)模型中的參數(shù)和原始數(shù)據(jù)，是由觀測和試驗(yàn)得到的 由于測量工具的精度、觀測方法或客觀條件的限制,使數(shù)據(jù)含有測量

6、誤差,這類誤差叫做觀測誤差或數(shù)據(jù)誤差 根據(jù)實(shí)際情況可以得到誤差上下界 數(shù)值方法中需要了解觀測誤差,以便選擇合理的數(shù)值方法與之適應(yīng)觀測誤差觀測誤差 精確公式用近似公式代替時(shí),所產(chǎn)生的誤差叫截?cái)嗾`差 例如, 函數(shù)f(x)用泰勒(Taylor)多項(xiàng)式 截?cái)嗾`差截?cái)嗾`差nnnxnfxfxffxp!)0(! 2)0(! 1)0()0()()(2 1) 1()!1()()()()(nnnnxnfxpxfxR（介于0與x之間）近似代替，則數(shù)值方法的截?cái)嗾`差是近似代替，則數(shù)值方法的截?cái)嗾`差是p 截?cái)嗾`差的大小直接影響計(jì)算結(jié)果的精度和計(jì)算截?cái)嗾`差的大小直接影響計(jì)算結(jié)果的精度和計(jì)算p 工作量，是數(shù)值計(jì)算中必須考慮

7、的一類誤差工作量，是數(shù)值計(jì)算中必須考慮的一類誤差 在數(shù)值計(jì)算中只能對有限位字長的數(shù)值進(jìn)行運(yùn)算 需要對參數(shù)、中間結(jié)果、最終結(jié)果作有限位字長的處理工作，這種處理工作稱作舍入處理 用有限位數(shù)字代替精確數(shù)，這種誤差叫做舍入誤差，是數(shù)值計(jì)算中必須考慮的一類誤差舍入誤差舍入誤差誤差誤差 例如在計(jì)算時(shí)用例如在計(jì)算時(shí)用3.141593.14159近似代替近似代替，產(chǎn)生的誤差產(chǎn)生的誤差R= R= -3.14159=0.0000026-3.14159=0.0000026就是舍入誤差。就是舍入誤差。 上述種種誤差都會影響計(jì)算結(jié)果的準(zhǔn)確上述種種誤差都會影響計(jì)算結(jié)果的準(zhǔn)確性，因此需要了解與研究誤差，在數(shù)值計(jì)算性，因此需

8、要了解與研究誤差，在數(shù)值計(jì)算中將著重研究截?cái)嗾`差、舍入誤差，并對它中將著重研究截?cái)嗾`差、舍入誤差，并對它們的傳播與積累作出分析們的傳播與積累作出分析二、誤差、有效數(shù)字二、誤差、有效數(shù)字定義定義1 1 絕對誤差，簡稱誤差：絕對誤差，簡稱誤差：.* ,*的近似值為準(zhǔn)確值其中xxxxe誤差限：誤差限：.|*|*的一個(gè)上界e相對誤差：相對誤差：,*xeer相對誤差限：相對誤差限：.|*的一個(gè)上界rre.*xeer或5 . 0765 x例如，毫米尺5.1000 1,10 yx例如，0.5%.| %,10|*yxyx.000008. 0 ,1416. 3 ,002. 0 ,14. 3 ,1415926.

9、3*5*5*3*3xxx取五位取三位定義定義2 2 .*,* 有效數(shù)字有效數(shù)字位有位，就說的第一位非零數(shù)字共有到該位的半個(gè)單位的誤差限是某一位數(shù)字若近似值nxnxx例例1 42.195, 0.0375551, 8.00033, 2.71828,1 42.195, 0.0375551, 8.00033, 2.71828,按四舍五按四舍五入寫出上述各數(shù)具有四位有效數(shù)字的近似數(shù)入寫出上述各數(shù)具有四位有效數(shù)字的近似數(shù). .(2.2) 1021* . 0(2.1) )1010(10* 11) 1(121nmnnmxxaaaax并且其中即例例2 2 考察三位有效數(shù)字重力加速度考察三位有效數(shù)字重力加速度g,

10、g,若以m/s2為單位, g9.80m/s2, 若以km/s2為單位, g0.00980km/s2,102180. 9g 23. 0, ) 1 . 2(nm，按,102100980. 0g 53. , 3 ) 1 . 2(nm，按.10212*1絕對誤差限.10215*2絕對誤差限.0.00980/0.0000050.005/9.80 *r而相對誤差限相同：11021*nm30 4 10.00009260.5 100.5 10 xx 解解1：若取近似值：若取近似值x*=3.1415，絕對誤差是，絕對誤差是0.0000926，有，有，即m=0,n4，故近似值x*=3.1415只有4位有效數(shù)字解2

11、：x*3.1415的絕對誤差限0.0005，它是x的小數(shù)后第3位的半個(gè)單位，故近似值x*=3.1415準(zhǔn)確到小數(shù)點(diǎn)后第3位故近似值x*=3.1415只有4位有效數(shù)字例例3 設(shè)設(shè)x= =3.1415926，求，求x*=3.1415的近似值及有效數(shù)字的近似值及有效數(shù)字1(1)121*(1)1 * *10(1010) (2.1)0 . *1 10 2mnnnrxxaaaaxna 設(shè)近似數(shù)表示為其中若具有 位有效數(shù)字，則其相對誤差限為；定理定理 *(1)11 * 10*2(1).nrxxan反之， 若的相對誤差限為，則至少具有 位有效數(shù)字 200.1%要使的相對誤差限小于，要取幾位有效數(shù)字？例例4 4

12、 111102na 1*111204.4,41100.125 100.1%42nnrana 只要取解：設(shè)取解：設(shè)取n位有效數(shù)字，相對誤差限位有效數(shù)字，相對誤差限*r=，1 31100.00252216110= 0.000 000 5629例例5 指出下列各數(shù)具有幾位有效數(shù)字，及其絕對誤差限和指出下列各數(shù)具有幾位有效數(shù)字，及其絕對誤差限和相對誤差限：相對誤差限：0.002 009 000.00解解 因?yàn)橐驗(yàn)閤1*=0.002 00， m=3絕對誤差限絕對誤差限0.000 005= 因?yàn)橐驗(yàn)閙=3，n=3， x1*= 0.002 00有有3位有效數(shù)字位有效數(shù)字. a1=2，相對誤差限相對誤差限r(nóng)=

13、x2*=9 000.00，絕對誤差限0.005，因?yàn)閙=3，n=6，x2*=9 000.00有6位有效數(shù)字，相對誤差限為r 如果認(rèn)為小數(shù)點(diǎn)后邊的0無用，將9 000.00隨便寫作90009103，那么它的絕對誤差就是=0.5=0.51034+1，即m=3，n=4，表明這個(gè)數(shù)有4位有效數(shù)字可見，小數(shù)點(diǎn)之后的0，不是可有可無的，它是有實(shí)際意義的.53 3 110.5 100.5 100.5 10m n 三、數(shù)值運(yùn)算的誤差估計(jì)三、數(shù)值運(yùn)算的誤差估計(jì)*1212,x xxx四則運(yùn)算，設(shè)為準(zhǔn)確值為近似值，則：誤差限.|)(|)(|)/( ),(|)(|)( ),()()( 2*2*1*2*2*1*2*1*

14、1*2*2*1*2*1*2*1*2*1xxxxxxxxxxxxxxxxx,*, ,*)(*)*)(*)()( ,*,)(22)(之間在公式由為近似值為準(zhǔn)確值，一元函數(shù)xxxxxxxfxfxfTaylorxxxff *).(|*)(|*)( *)(xxfxfxf的誤差限得).(*)( ),(,),(*1*11*11knkknnnnxxffxxfxxxxxxf的誤差限同理得的近似值為準(zhǔn)確值，多元函數(shù)* (8) ( *)( *)( *).sldsssldld場地面積：書上第 頁例例6 6ABC*(1000.10) , *(1200.10) ,*(600.02) ,ABCSobm cmA例 設(shè)觀測數(shù)據(jù)

15、為試估計(jì)面積的絕對誤差限和相對誤差限。21sin2()( )( )()11sin0.1sin0.1221cos0.0210.572180SbcASSSSbcAbcAcAbAbcAm解由則3( )10.57|( )| |2.035 101sin2rsssbcA誤差分析簡介誤差分析簡介 向后誤差分析法區(qū)間分析法概率分析法3 誤差定性分析、避免誤差危害誤差定性分析、避免誤差危害).,(),(111nnflnaagxaagx，xyyx ,一、病態(tài)問題與條件數(shù)一、病態(tài)問題與條件數(shù) ., ,)()()()(*)(條件數(shù)稱為計(jì)算函數(shù)值問題的考慮計(jì)算函數(shù)值問題ppxfxf xxxxfxfxfCC%.24%,2

16、,24. 1)02. 1 (, 1) 1 (,10,)(10函數(shù)值相對誤差為誤差為自變量相對例如ffCxxfp.10認(rèn)為是病態(tài)一般pC.,考慮是否病態(tài)條件數(shù)其他計(jì)算問題也要考慮二、算法的數(shù)值穩(wěn)定性二、算法的數(shù)值穩(wěn)定性考慮初始數(shù)據(jù)誤差在計(jì)算中的傳播問題. 1107 d , 0,1, .nxnIex ex n計(jì)算并估計(jì)誤差例例 ,.舍入義一個(gè)算法若輸入數(shù)據(jù)有誤差 而在計(jì)算過程中不增長 則稱此算法是數(shù)值穩(wěn)定的 否則是不誤差穩(wěn)定的定定3 3, 2 , 1 ,11nnIInn., 2 , 1 ,1,6321. 0)(10nInIIAnn.110eI. 1 , 8 , 9 ),1 ( ,0684. 0)(

17、*1*1*9nIIIBnnn)0684. 0)10101(21(19eI控制遞推公式中誤差的傳播控制遞推公式中誤差的傳播 對于一個(gè)數(shù)學(xué)問題的求解往往有多種數(shù)對于一個(gè)數(shù)學(xué)問題的求解往往有多種數(shù)值方法在選擇數(shù)值方法時(shí)，要注意所用的數(shù)值方法在選擇數(shù)值方法時(shí)，要注意所用的數(shù)值方法不應(yīng)將計(jì)算過程中難以避免的誤差放值方法不應(yīng)將計(jì)算過程中難以避免的誤差放大的較快，造成計(jì)算結(jié)果完全失真。大的較快，造成計(jì)算結(jié)果完全失真。例例 計(jì)算積分計(jì)算積分 并估計(jì)誤差并估計(jì)誤差解解 容易得到遞推公式容易得到遞推公式 10(0,1,2,10)10nnxIdxnx1 . 1ln)10ln(10110100 xdxxI 1.1ln

18、)10ln(10110100 xdxxI1011011011101010101010101010dxxxdxxxxdxxxxxdxxxInnnnnnnn11011011011010nnnIndxxxdxx)10, 2 , 1(n即即 為為 nI111 0(1, 2 ,1 0 )nnIInn則準(zhǔn)確的理論遞推式則準(zhǔn)確的理論遞推式 實(shí)際運(yùn)算的遞推式實(shí)際運(yùn)算的遞推式 兩式相減有兩式相減有 01101II*0*1101II)(10)(10*0*00*11IeIIII*2*112200( )10()( 10) ()( 1) 10 ()nnnnnnnnne IIIIIIIII 這就是說這就是說, ,若若 與

19、與 的誤差為的誤差為 = - , = - ,即即 ，則誤差的遞推規(guī)律為，則誤差的遞推規(guī)律為 0I*0I)(*0Ie0I*0I)(*0*00IeII于是于是 )(10)(10)(10)(*010*82*9*10IeIeIeIe計(jì)算計(jì)算 時(shí)的誤差被擴(kuò)大了時(shí)的誤差被擴(kuò)大了 倍倍, ,顯然算法是顯然算法是數(shù)值不穩(wěn)定的。數(shù)值不穩(wěn)定的。 如果將遞推公式如果將遞推公式 變換一種形式變換一種形式 *10I10101101nnInI101011nnInI準(zhǔn)確的理論遞推式準(zhǔn)確的理論遞推式實(shí)際運(yùn)算的遞推式實(shí)際運(yùn)算的遞推式從而有從而有 101011nnInI10101*1nnInI)(101*11nnnnIIII)(

20、10) 1()(101)(101*222*11*00nnnnIIIIIIII即即)(101)(101)(101)(*1010*22*1*0IeIeIeIe于是有于是有則這個(gè)算法的誤差傳遞規(guī)律為則這個(gè)算法的誤差傳遞規(guī)律為 *1()()10nne Ie I 即每計(jì)算一步的誤差的絕對值是上一步的十分即每計(jì)算一步的誤差的絕對值是上一步的十分之一，誤差的傳播逐步縮小，得到很好的控制，這之一，誤差的傳播逐步縮小，得到很好的控制，這個(gè)算法是數(shù)值穩(wěn)定的個(gè)算法是數(shù)值穩(wěn)定的 算法的數(shù)值穩(wěn)定性 算法優(yōu)劣的標(biāo)準(zhǔn) 從截?cái)嗾`差觀點(diǎn)看,算法必須是截?cái)嗾`差小,收斂斂速要快。即運(yùn)算量小,機(jī)器用時(shí)少. 從舍入誤差觀點(diǎn)看,舍入誤差

21、在計(jì)算過程中要能控制,即算法的數(shù)值要穩(wěn)定. 從實(shí)現(xiàn)算法的觀點(diǎn)看,算法的邏輯結(jié)構(gòu)不宜太復(fù)雜，便于程序編制和上機(jī)實(shí)現(xiàn). 設(shè)計(jì)算法時(shí)應(yīng)遵循的原則 要有數(shù)值要穩(wěn)定性,即能控制誤差的傳播. 避免大數(shù)吃小數(shù),即兩數(shù)相加時(shí),防止較小的數(shù)加不到較大的數(shù)上. 避免兩相近的數(shù)相減,以免有效數(shù)字的大量丟失. 避免分母很小(或乘法因子很大),以免產(chǎn)生溢出.三、避免誤差危害的若干原則三、避免誤差危害的若干原則除了分清問題是否病態(tài)和算法是否數(shù)值穩(wěn)定外,還要考慮避免誤差危害和防止有效數(shù)字損失的如下原則.1.避免大數(shù)除以小數(shù)例例8 仿計(jì)算機(jī)，采用仿計(jì)算機(jī)，采用3位十進(jìn)制，用消元法求解方程組位十進(jìn)制，用消元法求解方程組 51.

22、00 101.001.00 1.001.002.00 xyxy)1000. 100. 2()1000. 100. 1 ( 00. 100. 11000. 1555yyxx，得消00. 1 00. 100. 11000. 15yyx00. 1* ,00. 0* yx解：解：錯(cuò).為什么,怎么辦？9999899. 000001. 1555510110210110yx510) 1 () 2(減少運(yùn)算誤差原則減少運(yùn)算誤差原則2 2、兩個(gè)相近的數(shù)相減，會嚴(yán)重?fù)p失有效數(shù)字、兩個(gè)相近的數(shù)相減，會嚴(yán)重?fù)p失有效數(shù)字 例如例如x =1958.75x =1958.75，y =1958.32y =1958.32都具有五

23、位都具有五位 有效數(shù)字，但有效數(shù)字，但x-y=0.43x-y=0.43只有兩位有效數(shù)字只有兩位有效數(shù)字 通常采用的方法是改變計(jì)算公式通常采用的方法是改變計(jì)算公式, ,例如當(dāng)與例如當(dāng)與 很接近時(shí)很接近時(shí), ,由于由于2121lglglgxxxx用右端代替左端公式計(jì)算用右端代替左端公式計(jì)算, ,有效數(shù)字就不會損失有效數(shù)字就不會損失 減少運(yùn)算誤差原則減少運(yùn)算誤差原則當(dāng)當(dāng)x x很大時(shí)可作相應(yīng)的變換很大時(shí)可作相應(yīng)的變換 xxxx111) 1(11) 1(xxarctgarctgxxarctg則用右端來代替左端。則用右端來代替左端。 減少運(yùn)算誤差若干原則減少運(yùn)算誤差若干原則當(dāng)當(dāng)x x接近接近0 0時(shí)時(shí) x

24、xxxsin1sinsincos1一般情況，當(dāng)一般情況，當(dāng)f(x)f(xf(x)f(x* *) )時(shí)，可用泰勒展開時(shí)，可用泰勒展開 2*)(! 2)()()()(xxxfxxxfxfxf取右端的有限項(xiàng)近似左端。取右端的有限項(xiàng)近似左端。 如果計(jì)算公式不能改變，則可采用增加有效位如果計(jì)算公式不能改變，則可采用增加有效位數(shù)的方法保證精度數(shù)的方法保證精度 212 1610.x863 ,863xxx 求解例例9 972 A10 1 cos2 .1 cos2sin2xx。計(jì)算（）例例1010 xxxx11112301. 001. 001. 001. 001. 001. 0123100100項(xiàng)項(xiàng)例例11 仿

25、計(jì)算機(jī)在仿計(jì)算機(jī)在3位十進(jìn)制下，位十進(jìn)制下，、防止、防止大數(shù)大數(shù)吃吃小數(shù)小數(shù)例例 求二次方程求二次方程x2-105x+1=0 x2-105x+1=0的根的根 解：按二次方程求根公式解：按二次方程求根公式 x1=(105+(1010-4)1/2)/2 x1=(105+(1010-4)1/2)/2 x2=(105-(1010-4)1/2)/2 x2=(105-(1010-4)1/2)/2 在在8 8位浮點(diǎn)數(shù)計(jì)算得位浮點(diǎn)數(shù)計(jì)算得 x1=(105+105 )/2=105 ( x1=(105+105 )/2=105 (正確）正確）, , x2=(105-105 )/2=0 ( x2=(105-105 )

26、/2=0 (錯(cuò)誤）錯(cuò)誤）產(chǎn)生錯(cuò)誤的原因產(chǎn)生錯(cuò)誤的原因 出現(xiàn)大數(shù)出現(xiàn)大數(shù)10101010吃掉小數(shù)吃掉小數(shù)4 4的情況的情況 分子部分出現(xiàn)兩個(gè)相近數(shù)相減而喪失有分子部分出現(xiàn)兩個(gè)相近數(shù)相減而喪失有 效數(shù)位常稱為災(zāi)難性的抵消效數(shù)位常稱為災(zāi)難性的抵消4、絕對值太小的數(shù)不宜做除數(shù)當(dāng)分母為兩個(gè)相近數(shù)相減時(shí),會喪失有效數(shù)字4()()10 ()0.14560.14550.0001分子分子分子這里分子的誤差被擴(kuò)大這里分子的誤差被擴(kuò)大104104倍倍, ,再如再如若將分母變?yōu)槿魧⒎帜缸優(yōu)?.0011,0.0011,即分母只有即分母只有0.00010.0001的變化的變化時(shí)時(shí), ,計(jì)算結(jié)果卻有了很大變化計(jì)算結(jié)果卻有了

27、很大變化 減少運(yùn)算誤差若干原則減少運(yùn)算誤差若干原則3.14153141.50.0019 .28550011.01415.3例 計(jì)算0135. 00125. 00003. 00012. 00143. 00005. 0D 解: 分子分母分別計(jì)算后相除(取9位小數(shù))A=0.0005*0.0143*0.0012=0.00000715*0.0012 =0.000000009(有舍入)B=0.0003*0.0125*0.0135=0.00000375*0.0135 =0.000000051(有舍入)D=A/B=0.17647真值為0.16948148，所以D只準(zhǔn)確到小數(shù)后一位減少運(yùn)算誤差若干原則減少運(yùn)算誤

28、差若干原則 算法2。分成三組因子。每組只取六位小數(shù)計(jì)算 a=0.0005/0.0003=1.666667(有舍入) b=0.0143/0.0125=1.144000 c=0.0012/0.0135=0.088889 (有舍入) D=a*b*c=1. 666667* 1.144000* 0.088889 =0.169482,準(zhǔn)確到小數(shù)后5位。0135.00125.00003.00012.00143.00005.0Db bc ca a減少運(yùn)算誤差若干原則減少運(yùn)算誤差若干原則5、簡化計(jì)算步驟，減少運(yùn)算次數(shù)減少運(yùn)算次數(shù)可以不但節(jié)省時(shí)間，而且減少舍入誤差 例：x255=xx2x4x8x16x32x64x128 原先要做254次乘法現(xiàn)只需14次即可例 如計(jì)算多項(xiàng)式 p(x)=anxn an-1xn-1 a1x a0 的值 若直接計(jì)算akxk,再逐項(xiàng)相加，一共要做 n+(n-1)+2+1=n(n+1)/2次乘法和n次加法 減少運(yùn)算誤差若干原則減少運(yùn)算誤差若干原則如果