版權(quán)說明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請進(jìn)行舉報或認(rèn)領(lǐng)
文檔簡介
1、數(shù)值分析實驗分析報告姓名:學(xué)號:S20160258日期:2016.10.15班級:1602、實驗名稱誤差傳播與算法穩(wěn)定性、實驗?zāi)康捏w會穩(wěn)定性在選擇算法中的地位。誤差擴(kuò)張的算法是不穩(wěn)定的,是我們 所不期望的;誤差衰竭的算法是穩(wěn)定的,是我們努力尋求的,這是貫穿本課 程的目標(biāo)。、實驗內(nèi)容計算EnJ n x 1 |,0x e dx, n = 1,2,四、算法描述1n x_1En = 0x e dx, n =1,2,利用分部積分可得:En 丄xnexdxJldee 0e 01 、f x nA.An e x dx = 10Jn x nJe x dx = 1 - nEnJ,n e 02,3,可得遞推公式為:
2、1、En 1 = 1 - (n 1)En, n =1,2,3,n =N -1,N -2,3,2F面分別以1 , 2遞推關(guān)系求解方案一:En 1 = 1 -( n 1)En, n = 1,2,3,1當(dāng) n = 1 時 E1 = 1 _ _e1x 0e x dx0二-0.367879 e(保留六位有效數(shù)字)En =1 - nEn J巳=0.367879n =2,3,萬案一;n = N -1, N -2,3,2當(dāng) x (0,1)時,nXn x4. nv x e v xe1 n1x , dxv x0 e 01x'dxv xndx二0亠V EnV丄e(n 1)n 1當(dāng)n時,En > 0這里
3、取n = 20111e+1=E20、0.0325685 (保留六位有效數(shù)字)2 21e 2142eEnE20 -n 二 20,19,3,21 £n0.0325685五、程序流程圖由于實驗方案明顯、簡單,實現(xiàn)步驟及流程圖省略六、實驗結(jié)果計算結(jié)果如表1-1 :1-1計算結(jié)果表nEn*En10.3678790.36787920.2642410.26424130.2072770.20727740.1708930.17089350.1455330.14553360.1268020.12680270.1123840.11238480.1009320.10093290.09161200.09161
4、23100.08387700.0838771110.07735200.0773522120.07177500.0717733130.06692800.0669477140.06300200.0627322150.05496870.0590175160.1205000.055719217-1.0485070.0527731819.8731220.050085719-376.5893160.0483716207532.790.0325685七、實驗結(jié)果分析1、通過表1-1可以看出,算法一在前15項中迭代值基本保持一致,但是從16項開始就有了較大的差別。同時:f(x) _0二 f(x)dx_O1 又
5、當(dāng)x (0,1)時,xnexdx _ 0恒成立。0但是算法一中當(dāng)n=17、19時明顯不滿足。算法二相比較于算法一在目前的迭代次數(shù)下比較精確。2、設(shè)算法一中E1的計算誤差為e1,由E1遞推計算到En的誤差為en ;算法二中En的計算誤差為飛,由En向前遞推計算到E.( n : N )的誤差為。設(shè)算法一中的精確值為E-n,計算值為En ;算法二中的精確值為En* :,計 算值為En*根據(jù)題設(shè)可得:算法一:e?二 E2: E2 =2E/ 一(1 2EJ 二-2q遞推可得:en - -neni, n _ 2二 en =(-1 f'nen!NTN算法二:同理可得:1;nN MS3、根據(jù)算法一的公
6、式可得:當(dāng)n-;門時,en-;門,由此可見最后的結(jié)果誤差不僅僅取決于初始值的誤 差,還會隨著遞推公式的不斷運行,誤差在不斷的增大。當(dāng)N 時,;n 0,由此可見最后的結(jié)果誤差不僅僅取決于初始值的誤差,還會隨著遞推公式的不斷運行,誤差在不斷的減小。4、通過前三問,可以發(fā)現(xiàn)算法二比算法一更加穩(wěn)定八、附錄(程序)算法一:fun ctio n ek= jife ng(e0 ,n )e(1)=vpa(e0,6);B=zeros( n, 2);B(:,1)=1: n;for i=1: ne(i+1)=1-(i+1)*e(i);ek=e(i+1);B(i,2)=vpa(ek,6);endxlswrite(
7、39;算法一 ',B,'sheet1')算法二:function i,ek = jife ng1( eO,n )e(n+1)=e0;B=zeros( n, 2);B(:,1)=1: n;for i=n:-1:1e(i)=(1-e(i+1)/(i+1);ek=e(i);B(i,2)=vpa(ek,6);endxlswrite('算法二',B,'sheet1')輸出結(jié)果:算法一計算截圖:A1BC1' 10. 264241220. 207277330* 170893446 1455335112384770 1009328 1L8二門 9129 30. 083S7710100. 07735211 110. 07177512120. 06692813130- 0630021414CL 05496915150. 12051616一104匪1171719 8731218匚1呂-376, 58919匚 197532. 786算法二計算截圖:1O. 36787920. 2642413O. 2072774O.17089350.1455336O. 12S8O27O. 11238480 100
溫馨提示
- 1. 本站所有資源如無特殊說明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
- 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
- 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁內(nèi)容里面會有圖紙預(yù)覽,若沒有圖紙預(yù)覽就沒有圖紙。
- 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
- 5. 人人文庫網(wǎng)僅提供信息存儲空間,僅對用戶上傳內(nèi)容的表現(xiàn)方式做保護(hù)處理,對用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對任何下載內(nèi)容負(fù)責(zé)。
- 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容,請與我們聯(lián)系,我們立即糾正。
- 7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時也不承擔(dān)用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。
最新文檔
- 《模具制造工藝學(xué)》教學(xué)大綱
- 教案裝訂順序
- 四個自信課件
- 玉溪師范學(xué)院《現(xiàn)代教育技術(shù)》2022-2023學(xué)年第一學(xué)期期末試卷
- 玉溪師范學(xué)院《田徑》2021-2022學(xué)年第一學(xué)期期末試卷
- 教練員繼續(xù)教育考試題目及答案-知識題庫
- 湖南師大附中2024-25屆高三年級月考試卷(二)(英語)
- 電商公司整體薪酬設(shè)計(早期)
- 《信號基礎(chǔ)設(shè)備》全套教學(xué)課件
- 2023年雙頻、雙模移動通信手機項目綜合評估報告
- 《老年人生活照護(hù)》試卷A卷及答案
- 消防安全知識培訓(xùn)課件
- 高中歷史選擇性必修2知識點總結(jié)歸納
- 16J914-1 公用建筑衛(wèi)生間
- 物聯(lián)網(wǎng)應(yīng)用技術(shù)職業(yè)生涯規(guī)劃
- 2024年廣東恒健投資控股有限公司招聘筆試參考題庫含答案解析
- 2型糖尿病人生活質(zhì)量評定量表(中文版)
- GR326CORE規(guī)范講解
- 新訓(xùn)工作總結(jié)(共5篇)
- 五年級下冊牛津英語期中試卷【精】
- 《垃圾分類》ppt課件
評論
0/150
提交評論