




免費(fèi)預(yù)覽已結(jié)束,剩余1頁(yè)可下載查看
下載本文檔
版權(quán)說(shuō)明:本文檔由用戶(hù)提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請(qǐng)進(jìn)行舉報(bào)或認(rèn)領(lǐng)
文檔簡(jiǎn)介
計(jì)算機(jī)科學(xué)與工程學(xué)院計(jì)算方法實(shí)驗(yàn)報(bào)告學(xué)號(hào) 姓名班級(jí)實(shí)驗(yàn)項(xiàng)目名稱(chēng)實(shí)驗(yàn)二 數(shù)值積分一、實(shí)驗(yàn)名稱(chēng)實(shí)驗(yàn)二數(shù)值積分二、 實(shí)驗(yàn)?zāi)康模?1) 熟悉復(fù)化梯形方法、復(fù)化Simpson方法、梯形遞推算法、龍貝格算法;(2) 能編程實(shí)現(xiàn)復(fù)化梯形方法、復(fù)化Simpson方法、梯形遞推算法、龍貝格算法;(3) 理解并掌握自適應(yīng)算法和收斂加速算法的基本思想;(4) 分析實(shí)驗(yàn)結(jié)果體會(huì)各種方法的精確度,建立計(jì)算機(jī)求解定積分問(wèn)題的感性認(rèn)識(shí)三、 實(shí)驗(yàn)內(nèi)容及要求(1)設(shè)計(jì)復(fù)化梯形公式求積算法,編制并調(diào)試相應(yīng)的函數(shù)子程序(2)設(shè)計(jì)復(fù)化辛浦生求積算法,編制并調(diào)試相應(yīng)的函數(shù)子程序(3)用龍貝格算法計(jì)算輸入:積分區(qū)間,誤差限輸出:序列Tn,Sn,Cn,Rn及積分結(jié)果(參考書(shū)本P81的表2-5) 取n=2,4,8,16,精確解為0.9460831 四、 實(shí)驗(yàn)原理及算法描述在許多實(shí)際問(wèn)題中,常常需要計(jì)算定積分的值。根據(jù)微積分學(xué)基本定理,若被積函數(shù)f(x)在區(qū)間a,b上連續(xù),只要能找到f(x)的一個(gè)原函數(shù)F(x),便可利用牛頓-萊布尼茲公式求得積分值。但是在實(shí)際使用中,往往遇到如下困難,而不能使用牛頓-萊布尼茲公式。(1) 找不到用初等函數(shù)表示的原函數(shù)(2) 雖然找到了原函數(shù),但因表達(dá)式過(guò)于復(fù)雜而不便計(jì)算(3) f(x)是由測(cè)量或計(jì)算得到的表格函數(shù)由于以上種種困難,有必要研究積分的數(shù)值計(jì)算問(wèn)題。利用插值多項(xiàng)式 則積分轉(zhuǎn)化為,顯然易算。稱(chēng)為插值型求積公式。最簡(jiǎn)單的插值型求積公式是梯形公式和Simpson公式,。當(dāng)求積結(jié)點(diǎn)提供較多,可以分段使用少結(jié)點(diǎn)的梯形公式和Simpson公式,并稱(chēng)為復(fù)化梯形公式、復(fù)化Simpson公式。如步長(zhǎng)未知,可以通過(guò)誤差限的控制用區(qū)間逐次分半的策略自動(dòng)選取步長(zhǎng)的方法稱(chēng)自適應(yīng)算法。梯形遞推公式給出了區(qū)間分半前后的遞推關(guān)系。由梯形遞推公式求得梯形序列,相鄰序列值作線性組合得Simpson序列, Simpson序列作線性組合得柯特斯序列, 柯特斯序列作線性組合的龍貝格序列。若|R2-R1|e,則輸出R2;否則依此類(lèi)推。如此加工數(shù)據(jù)的過(guò)程叫龍貝格算法,如下圖所示:復(fù)化梯形公式復(fù)化Simpson公式梯形遞推公式加權(quán)平均公式: 龍貝格算法大大加快了誤差收斂的速度,由梯形序列O(h2) 提高到龍貝格序列的O(h8)五、 程序代碼及實(shí)驗(yàn)結(jié)果1 主程序int main()/cout longbeige(0, 1, 0.0000001);cout 請(qǐng)輸入 你的區(qū)間 和 誤差限度: x y z;cout 根據(jù)龍貝格算法 求出的精確值為: longbeige(x, y, z) endl;cout K T S C R endl;cout 0 tixing(0, 1, 1) endl;cout 1 tixing(0, 1, 2) xingbusheng(0, 1, 1) endl;cout 2 tixing(0, 1, 4) xingbusheng(0, 1, 2) (double (16/15)* xingbusheng(0, 1, 2) -(double (1/15)*xingbusheng(0, 1, 1) endl;cout 3 tixing(0,1,8) xingbusheng(0,1,4) (double(16 / 15)* xingbusheng(0, 1, 4) - (double(1 / 15)*xingbusheng(0, 1, 2) longbeige(0, 1, 0.000000001) endl;cout 4 tixing(0, 1, 16) xingbusheng(0, 1, 8) (double(16 / 15)* xingbusheng(0, 1, 4) - (double(1 / 15)*xingbusheng(0, 1, 2) longbeige(0, 1, 0.00000001) endl;cout the result 4 is longbeige(0, 1, 0.00000001);return 0;2 復(fù)化梯形公式子程序:double tixing(double a, double b, int n)double fa = f(a);double fb = f(b);double h = (b-a) / n;double fxk = 0.0;for (int k = 1; k = n - 1; k+)double xk = a + k*h;fxk = f(xk) + fxk;double res = (h / 2.0)*(fa + 2.0*fxk + fb);return res;3 復(fù)化辛浦生公式子程序:double xingbusheng(double a, double b, int n)double h = (b - a) /n;double fa = f(a);double fb = f(b);/double s = fb - fa;double x = a;double fxk12 = 0.0;double fxk = 0.0;/*for (int k = 1; k =n; k+)x = x + h / (2.0);double fx = f(x);s = s + 4.0*f(x);x = x + h / (2.0);s = s + 2.0*f(x);s = (h / 6.0)*s;return s;*/for (int k = 0; k = n - 1; k+)double xk1 = a + k*h;x = xk1 + h / 2.0;fxk12 = fxk12 + f(x);for (int k = 1; k = n - 1; k+)double xk = a + k*h;fxk = fxk + f(xk);double s = (h / 6)*(fa + 4 * fxk12 + 2 * fxk + fb);return s;4 龍貝格公式子程序:double longbeige(double a, double b, double wuchaxian)double h = b - a;double T1 = (h / 2.0)*(f(a) + f(b);int k = 1;double S,x,T2,S2,S1=0,C1=0,C2,R1=0,R2;double Tck100,Sck100,Cck100,Rck100;loop1: S = 0; x = a + h / 2.0;while (x= wuchaxian)R1 = R2;C1 = C2;k+;h = h / 2;T1 = T2;S1 = S2;goto loop1;elsereturn R2;實(shí)驗(yàn)結(jié)果: 圖1 六、 實(shí)驗(yàn)總結(jié)1. 梯形公式的收斂速度太慢,所以我們才會(huì)選
溫馨提示
- 1. 本站所有資源如無(wú)特殊說(shuō)明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請(qǐng)下載最新的WinRAR軟件解壓。
- 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請(qǐng)聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶(hù)所有。
- 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁(yè)內(nèi)容里面會(huì)有圖紙預(yù)覽,若沒(méi)有圖紙預(yù)覽就沒(méi)有圖紙。
- 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
- 5. 人人文庫(kù)網(wǎng)僅提供信息存儲(chǔ)空間,僅對(duì)用戶(hù)上傳內(nèi)容的表現(xiàn)方式做保護(hù)處理,對(duì)用戶(hù)上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對(duì)任何下載內(nèi)容負(fù)責(zé)。
- 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容,請(qǐng)與我們聯(lián)系,我們立即糾正。
- 7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時(shí)也不承擔(dān)用戶(hù)因使用這些下載資源對(duì)自己和他人造成任何形式的傷害或損失。
最新文檔
- 十堰活禽屠宰管理辦法
- 設(shè)計(jì)施工新技術(shù)
- 江蘇房屋維修管理辦法
- 發(fā)票違法檢舉管理辦法
- 職業(yè)規(guī)劃與就業(yè)指導(dǎo)教程
- 農(nóng)家栽培紅薯管理辦法
- 村級(jí)項(xiàng)目立項(xiàng)管理辦法
- 道路改造與混凝土管鋪設(shè)施工方案設(shè)計(jì)及舊路面拆除策略探討
- 醫(yī)用織物清洗管理辦法
- 杭州交警頭盔管理辦法
- 來(lái)料檢驗(yàn)規(guī)范
- 電鍍產(chǎn)品檢驗(yàn)記錄
- 2023-2024學(xué)年遼寧省大連市小學(xué)語(yǔ)文五年級(jí)期末評(píng)估試卷附參考答案和詳細(xì)解析
- 2023年小學(xué)數(shù)學(xué)必背定義和公式
- 2023年四川省宜賓市全科醫(yī)學(xué)專(zhuān)業(yè)實(shí)踐技能測(cè)試卷(含答案)
- 電梯井道腳手架施工方案
- 興平市生活垃圾焚燒發(fā)電項(xiàng)目環(huán)評(píng)報(bào)告
- 主令電器(課用)課件
- 湘少版英語(yǔ)六年級(jí)下冊(cè)全冊(cè)教案
- 湖南省長(zhǎng)郡中學(xué)“澄池”杯數(shù)學(xué)競(jìng)賽初賽試題(掃描版含答案)
- 消防系統(tǒng)施工總進(jìn)度計(jì)劃
評(píng)論
0/150
提交評(píng)論