版權(quán)說(shuō)明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請(qǐng)進(jìn)行舉報(bào)或認(rèn)領(lǐng)
文檔簡(jiǎn)介
1、實(shí)驗(yàn)五一、實(shí)驗(yàn)名稱 復(fù)合梯形積分和復(fù)合Simpson積分計(jì)算數(shù)值積分二、實(shí)驗(yàn)?zāi)康呐c要求:實(shí)驗(yàn)?zāi)康? 掌握復(fù)合梯形積分和復(fù)合Simpson積分算法。實(shí)驗(yàn)要求:1.給出復(fù)合梯形積分和復(fù)合Simpson積分算法思路,2.用C語(yǔ)言實(shí)現(xiàn)算法,運(yùn)行環(huán)境為Microsoft Visual C+。三、算法思路:我們把整個(gè)積分區(qū)間a,b分成n個(gè)子區(qū)間xi,xi+1,i=0,1,2,n,其中x0=a,xn+1=b。這樣求定積分問(wèn)題就分解為求和問(wèn)題: 當(dāng)這n+1個(gè)結(jié)點(diǎn)為等距結(jié)點(diǎn)時(shí),即,i=0,1,2,n,復(fù)化梯形公式的形式是算法:input nfor i=1 to n do end do output S如果n還是
2、一個(gè)偶數(shù),則復(fù)合Simpson積分的形式是算法:input nfor i=1 to n/2 do end do output S 四、實(shí)驗(yàn)題目:五、問(wèn)題的解:編寫程序(程序見(jiàn)后面附錄),輸出結(jié)果如下:為了便于看清數(shù)值積分結(jié)果與原函數(shù)積分實(shí)際結(jié)果的差異。我在運(yùn)行程序時(shí)故意計(jì)算了一下原函數(shù)積分的實(shí)際結(jié)果。分析并比較得到的數(shù)據(jù)可以看出,當(dāng)k越來(lái)越大時(shí),數(shù)值積分的結(jié)果越來(lái)越靠近原函數(shù)積分實(shí)際結(jié)果,并且復(fù)合Simpson積分的結(jié)果更迅速地靠近原函數(shù)積分實(shí)際結(jié)果,這是有原因的,從兩種方法的誤差項(xiàng)即可看出。復(fù)合梯形積分的誤差項(xiàng)是,Simpson積分的誤差項(xiàng)是,當(dāng)h趨于零時(shí),顯然Simpson積分的誤差項(xiàng)更快
3、地趨于零,實(shí)驗(yàn)結(jié)果復(fù)符合這一結(jié)論。六、附錄:實(shí)驗(yàn)編程,運(yùn)行環(huán)境為Microsoft Visual C+#include <math.h>#include <stdio.h>#include <stdlib.h>double f(double x) /定義函數(shù)f(x)/double y;y=sin(x);return(y);double S1(int N,double a,double b) /建立復(fù)合梯形積分/double s,h;int i;h=(b-a)/N; s=0.0;for(i=1;i<=N;i+)s=s+h*(f(a+(i-1)*h)+f(
4、a+i*h)/2.0;return(s);double S2(int N,double a,double b) /建立符合Simpson積分/double s,h;int i;h=(b-a)/N; s=0.0;for(i=1;i<=N/2;i+)s=s+h*(f(a+(2*i-2)*h)+4*f(a+(2*i-1)*h)+f(a+2*i*h)/3.0;return(s);void main() /main函數(shù)進(jìn)行最終運(yùn)算并輸出結(jié)果/int k;printf("s1代表復(fù)合梯形積分結(jié)果,s2代表Simpson積分結(jié)果n");for(k=0;k<=12;k+) printf("k=%d時(shí),s1=%.12f , s2=%.12f n",k,S1(pow(2,k),0.0,4.0),
溫馨提示
- 1. 本站所有資源如無(wú)特殊說(shuō)明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請(qǐng)下載最新的WinRAR軟件解壓。
- 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請(qǐng)聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
- 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ì)用戶上傳內(nèi)容的表現(xiàn)方式做保護(hù)處理,對(duì)用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對(duì)任何下載內(nèi)容負(fù)責(zé)。
- 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容,請(qǐng)與我們聯(lián)系,我們立即糾正。
- 7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時(shí)也不承擔(dān)用戶因使用這些下載資源對(duì)自己和他人造成任何形式的傷害或損失。
最新文檔
- 含有防曬劑的美容貼片相關(guān)項(xiàng)目建議書
- 曲軸相關(guān)項(xiàng)目建議書
- 網(wǎng)絡(luò)虛擬化軟件
- 江蘇省無(wú)錫市輔仁中學(xué)2024-2025學(xué)年九上數(shù)學(xué)開學(xué)質(zhì)量跟蹤監(jiān)視試題【含答案】
- 2020-2021學(xué)年湖南師大附中博才實(shí)驗(yàn)中學(xué)七年級(jí)(下)期中數(shù)學(xué)試卷(附答案詳解)
- 衛(wèi)生毒理學(xué)練習(xí)試卷8(題后含答案及解析)
- 軟件產(chǎn)品分發(fā)運(yùn)輸合同
- 煤礦瓦斯氣體運(yùn)輸服務(wù)
- 花鳥市場(chǎng)裝修工程設(shè)計(jì)合約
- 環(huán)保工程居間協(xié)議模板
- 北京四中初一年級(jí)期中語(yǔ)文試題
- 運(yùn)營(yíng)主管崗位招聘筆試題與參考答案(某大型集團(tuán)公司)
- 2024年事業(yè)單位考試(綜合管理類A類)職業(yè)能力傾向測(cè)驗(yàn)試卷及答案指導(dǎo)
- 滬教版三年級(jí)上冊(cè)英語(yǔ)單詞表
- 【課件】跨學(xué)科實(shí)踐:制作隔音房間模型人教版物理八年級(jí)上冊(cè)
- 《外科學(xué)》教案:第四十二章 門靜脈高壓癥
- 2024二十屆三中全會(huì)知識(shí)競(jìng)賽題庫(kù)及答案
- 2024年江蘇省昆山市自然資源和規(guī)劃局招聘編外13人歷年(高頻重點(diǎn)復(fù)習(xí)提升訓(xùn)練)共500題附帶答案詳解
- -投標(biāo)技術(shù)標(biāo)書范文模板-人員配備與團(tuán)隊(duì)構(gòu)建
- 2024年輔警招聘考試試題庫(kù)及完整答案(全優(yōu))
- 統(tǒng)編版高一語(yǔ)文必修上冊(cè)期末復(fù)習(xí):文言文閱讀 練習(xí)題匯編(含答案解析)
評(píng)論
0/150
提交評(píng)論