幾種常用數(shù)值積分方法的比較匯總

1、貴州師范學(xué)院本科畢業(yè)論文（設(shè)計(jì)）誠(chéng)信聲明學(xué)科分類號(hào)110.3420GUIZHOU NORMAL COLLEGE本科畢業(yè)論文題 目一幾種常用數(shù)值積分方法的比較潘曉祥 學(xué)號(hào)1006020540200院（系）數(shù)學(xué)與計(jì)算機(jī)科學(xué)學(xué)院專業(yè)數(shù)學(xué)與應(yīng)用數(shù)學(xué)年級(jí)2010 級(jí)講師指導(dǎo)教師雍進(jìn)軍職稱二0四年五月本人鄭重聲明：所呈交的本科畢業(yè)論文（設(shè)計(jì)），是本人在指導(dǎo)老師 的指導(dǎo)下，獨(dú)立進(jìn)行研究工作所取得的成果，成果不存在知識(shí)產(chǎn)權(quán)爭(zhēng)議，除文中已經(jīng)注明引用的內(nèi)容外，本論文不含任何其他個(gè)人或集體已經(jīng)發(fā)表 或撰寫(xiě)過(guò)的作品成果。對(duì)本文的研究做出重要貢獻(xiàn)的個(gè)人和集體均已在文中以明確方式標(biāo)明。本人完全意識(shí)到本聲明的法律結(jié)果由本人

2、承擔(dān)。本科畢業(yè)論文作者簽名：年月曰貴州師范學(xué)院本科畢業(yè)論文（設(shè)計(jì)）任務(wù)書(shū)畢業(yè)設(shè)計(jì)題冃幾種常用數(shù)值積分方法的比較作者姓名潘曉祥學(xué)號(hào)1006020540200年級(jí)2010 級(jí)所屬學(xué)院數(shù)學(xué)與計(jì)算機(jī)科學(xué)學(xué)院專業(yè)數(shù)學(xué)與應(yīng)用數(shù)學(xué)班級(jí)四班指導(dǎo)教師簽名雍進(jìn)軍講師職稱講師開(kāi)題日期2013年7月10日主要目標(biāo)1了解什么數(shù)值積分基本思想和一些常用的數(shù)值積分方法；2對(duì)各種數(shù)值積分方法的誤差以及代數(shù)精度進(jìn)行分析；3對(duì)各積分方法進(jìn)行比較總結(jié)出優(yōu)缺點(diǎn)。通過(guò)對(duì)幾種常用的數(shù)值積分方法進(jìn)行了的分折*并用這幾種方法對(duì)被積函數(shù)是普通函數(shù)做了數(shù)值積分，并在計(jì)算機(jī)上進(jìn)行實(shí)驗(yàn)。數(shù)值積分是計(jì)算方法或數(shù)值分析理論中非常重要的內(nèi)容,數(shù)值積分方法

3、也 是解決實(shí)際計(jì)算問(wèn)題的重要方法'對(duì)幾種常用數(shù)值積分方法的分析很必要。貴州師范學(xué)院本科畢業(yè)論文（設(shè)計(jì)）開(kāi)題報(bào)告書(shū)本文通過(guò)對(duì)復(fù)化求積公式，Newton -Coles求積公式，Romberg求積公式，咼斯型求積公式 進(jìn)行分 析討論并在計(jì)算機(jī)上積分實(shí)驗(yàn)，從代數(shù)精度,求積公式i吳差等角度對(duì)這些方法進(jìn)行分析比較'并總結(jié)出毎種求積分法的優(yōu)缺點(diǎn)以及實(shí)用性。貴州師范學(xué)院本科畢業(yè)論文（設(shè)計(jì)）開(kāi)題報(bào)告書(shū)論文題冃幾種常用數(shù)值積分方法的比較作者姓名潘曉祥學(xué)號(hào)1006020540200年級(jí)2010 級(jí)所屬學(xué)院指導(dǎo)教師姓名數(shù)學(xué)與計(jì)算機(jī)科學(xué)學(xué)院雍進(jìn)軍專業(yè)職稱題印性質(zhì)應(yīng)用硏究數(shù)學(xué)與應(yīng)用數(shù)學(xué)班級(jí)數(shù)本（4 ）班講

4、師預(yù)計(jì)字?jǐn)?shù)5000.00 字日期2013年7月05日24選題的原由： 硏究意義：數(shù)值積分是數(shù)學(xué)上的重要課題之一,是數(shù)值分析申的重要內(nèi)容之一,也是數(shù)學(xué)的硏究重點(diǎn)并在實(shí)際問(wèn)題及應(yīng)用中有 著廣泛的應(yīng)用常用于科學(xué)與工程的計(jì)算中，如涉及到積分方程，工程計(jì)算，計(jì)算機(jī)圖形學(xué)，金融數(shù)學(xué)等應(yīng)用科學(xué)領(lǐng)域都有著相當(dāng)重要的應(yīng)用，所以研究數(shù)值積分問(wèn)題有很重要的意義數(shù)值積分是硏究如何求出一個(gè)積分的數(shù)值這一課題的超源可追溯到古代，其中一個(gè)突出的例子是希臘人用內(nèi)接與外接正多邊形推算出圖面積的方法也正是此法使阿基米德得以求出n值得上界與下界，若干世紀(jì)以來(lái)，尤其是十六世紀(jì)后，已提出了多彳中數(shù)值積分方法，其中有矩形求積法，內(nèi)橢求積

5、法， 牛頓科待斯公式，復(fù)化求積公式，龍貝格求積公式，高斯型求積公式但各種方法都有特點(diǎn)，在不冋的情況下試用程 度不冋，我們將著重從求積公式的代數(shù)精度和余項(xiàng)等角度對(duì)這些方法進(jìn)行分析比較 硏究動(dòng)態(tài)2這些年來(lái)，有矢數(shù)值積分的硏究已經(jīng)成為一個(gè)很活躍的硏究領(lǐng)域，歷史上，阿基米徳，牛頓，歐拉，高斯，切比雪夫等人都對(duì)此有過(guò)貢獻(xiàn)研究出各種各樣的數(shù)值求積公式，但一個(gè)好的數(shù)值求積公式應(yīng)該滿足：計(jì)算簡(jiǎn)單，誤差小，代數(shù)精度高我們將對(duì)矩形求積法，內(nèi)橢求積法，牛頓科特斯公式，化求積公式，貝格求積公式，斯型求積公式進(jìn)行比較對(duì)數(shù)值求積公式能有進(jìn)一步的了解和學(xué)習(xí)1數(shù)值積分方法的基本思想2幾類常用數(shù)值積分方法的基本分析2.1主要內(nèi)

6、容：Newt on - Cotes求積公式復(fù)化求積公式2.3Romberg求積公式高斯型求積公式3幾類數(shù)值積分方法的簡(jiǎn)單比較評(píng)述4利用MATLAB編程應(yīng)用對(duì)幾類求積算法的分析比較研究方法:本論文主雯通過(guò)對(duì)相矣文獻(xiàn)和書(shū)籍的參考，合自己的見(jiàn)解，復(fù)化求積公式，Newton - Cotes求積公式,Romberg求積公式，高斯型求積公式逬行討論并進(jìn)行上機(jī)實(shí)驗(yàn)- 角度對(duì)這些方法進(jìn)行分析比較完成期限和采取的主要措施：從代數(shù)精度，求積公式誤差等本論文計(jì)劃用6個(gè)月的時(shí)問(wèn)完成，階段的任務(wù)如下:7月份査閱相矢書(shū)籍和文獻(xiàn);8月份完成開(kāi)題報(bào)告并交老師批閱9月份完成論文初槁并交老師批閱10月份完成論文二搞并交老師批閱主

7、要措施：考相矣書(shū)藉和文獻(xiàn)，合自己的見(jiàn)解，老師的指導(dǎo)下和同學(xué)的幫助下完成11月份完成論文三搞;12月份定稿.主要參考文獻(xiàn)及資料名稱:1矢治陸金甫數(shù)學(xué)分析基礎(chǔ)（第二版）IM北京：等教育出版社-20107胡祖熾林源渠.數(shù)值分析M北京：等教育出版社.4 986.3薛毅數(shù)學(xué)分析與實(shí)驗(yàn)M北京：業(yè)大學(xué)出版社2005.34徐士良數(shù)值分析與算法M北京：械工業(yè)出版社2007.1王開(kāi)榮楊大地.應(yīng)用數(shù)值分析M北京:等教育出版社2010.7楊一都數(shù)值計(jì)算方法M北京:等教育出版社.2008478910韓明王家寶李林?jǐn)?shù)學(xué)實(shí)驗(yàn)（MATLAB版上海：濟(jì)大學(xué)出版社2012J 圣寶建矢于數(shù)值積分若干問(wèn)題的硏究J南京信息工程大學(xué).2

8、009.05.01. : 42 劉緒軍.幾種求積公式計(jì)算精確度的比較.南京職業(yè)技術(shù)學(xué)院.2009.史萬(wàn)明吳裕樹(shù)孫新數(shù)值分析ML北京理工大學(xué)出版社.2010 A指導(dǎo)教師意見(jiàn):簽名：開(kāi)題報(bào)告會(huì)紀(jì)要時(shí)間2013年8月26日地點(diǎn)寧?kù)o樓229教師辦公室與會(huì)人員姓名職務(wù)（職稱）姓名職務(wù)（職稱）姓名職務(wù)（職稱）雍進(jìn)軍導(dǎo)師（講師）鄧喜才副教授李晟副教授龍林林組長(zhǎng)會(huì)議記錄摘要：指導(dǎo)小組針對(duì)課題 二次函數(shù)性質(zhì)的應(yīng)用提問(wèn)了以下問(wèn)題以及報(bào)告人的回答-雍老師問(wèn)：選擇此題冃的冃的？潘曉祥答：隨著計(jì)算機(jī)和計(jì)算方法的飛速發(fā)展,幾乎所有學(xué)科都走向定量化和精確化中的數(shù)值計(jì)算方法則是解決“計(jì)算"問(wèn)題的橋梁和工具。鄧?yán)蠋焹?/p>

9、：對(duì)這個(gè)冋題進(jìn)行研究有什么實(shí)際的意義,又有實(shí)用性和實(shí)驗(yàn)性的技術(shù)待潘曉祥答：計(jì)算方法既有數(shù)學(xué)類課程中理論上的抽僉性和嚴(yán)謹(jǐn)性征，計(jì)算方法是一門(mén)理論性和實(shí)踐性都很強(qiáng)的學(xué)科在科學(xué)研究和工程技術(shù)中都要用到各種計(jì)算方法例如在航天航空、地質(zhì)勘探、汽車(chē)制造、 的蹤影。橋梁設(shè)計(jì)、天氣預(yù)報(bào)和漢字字樣設(shè)計(jì)中都有計(jì)算方法李老師冋：對(duì)這個(gè)冋題你有什么自己的看法？潘曉祥答：隨著計(jì)算機(jī)技術(shù)的迅速發(fā)展和普及,現(xiàn)在計(jì)算方法課程幾乎已成為所有理工科學(xué)生的必修課程我們知道，計(jì)算能力是計(jì)算工具和計(jì)算方法的效率的乘積,提高計(jì)算方法的效率與提高計(jì)算機(jī)硬件的效率冋樣重要科學(xué)計(jì)算已用到科學(xué)技術(shù)和社會(huì)生活的各個(gè)領(lǐng)域中所以，研究數(shù)值計(jì)算方法可以

10、讓數(shù)學(xué)的應(yīng)用更大更廣。會(huì)議主持人簽名：記錄人簽名：負(fù)責(zé)人簽名：負(fù)責(zé)人簽名：貴州師范學(xué)院數(shù)學(xué)與計(jì)算機(jī)科學(xué)學(xué)院指導(dǎo)教師指導(dǎo)本科畢業(yè)論文情況登記表論文（設(shè)計(jì)）題冃幾種常用數(shù)值積分方法的比較用學(xué)生姓名潘曉祥學(xué)號(hào)1006020540200年級(jí)2010 級(jí)所屬學(xué)院數(shù)計(jì)學(xué)院專業(yè)數(shù)學(xué)與應(yīng)用數(shù)學(xué)班級(jí)四班指導(dǎo)教師姓名雍進(jìn)軍職稱講師學(xué)歷碩士指導(dǎo)時(shí)間指導(dǎo)地點(diǎn)指導(dǎo)內(nèi)容導(dǎo)教師簽名備注2013 年 06 月 100致遠(yuǎn)樓416論文選題-資料準(zhǔn)備面授2013年06月22日網(wǎng)上確定畢業(yè)論文選題電子郵件2013年06月26日網(wǎng)上怎樣撰寫(xiě)畢業(yè)論文開(kāi)題報(bào)告電子郵件2013年06月28日網(wǎng)上指導(dǎo)學(xué)生撰寫(xiě)開(kāi)題報(bào)告電子郵件2013年07月

11、14日網(wǎng)上幫助學(xué)生查找有矢參考文獻(xiàn)電子郵件2013年07月仃日手機(jī)女M町構(gòu)思自己的畢業(yè)論文手機(jī)飛信2013年08月21日手機(jī)聽(tīng)取學(xué)生畢業(yè)論文寫(xiě)作進(jìn)展情況匯報(bào)手機(jī)飛信2013年08月28日網(wǎng)上解答學(xué)生在論文寫(xiě)作中遇到的疑惑電子郵件2013年9月09日網(wǎng)上幫助學(xué)生查找有矢參考文獻(xiàn)電子郵件2013年們?cè)?8日網(wǎng)上女M可規(guī)劃自己的論文電子郵件2013年12月04日手機(jī)怎樣寫(xiě)好論文引言手機(jī)飛信2013年12月08日網(wǎng)上怎樣寫(xiě)好論文引言電子郵件2013 年 12 月 120網(wǎng)上怎樣寫(xiě)論文摘要電子郵件2013 年 12 月 160網(wǎng)上怎樣選取論文矢鍵詞電子郵件2013年12月20日網(wǎng)上怎樣編輯論文中的公式電

12、子郵件2014年01月05日手機(jī)督促學(xué)生在寒假中寫(xiě)好論文的初稿電子郵件2014年02月27日寧辭樓2伯檢查學(xué)生論文完成悄況面授2014年03月03日寧辭樓2伯對(duì)學(xué)生的論文初犒提出修改時(shí)意見(jiàn)面授2014年03月07日寧辭樓2伯解答學(xué)生在修改時(shí)的閑惑面授2014年03月11日寧辭樓2伯指導(dǎo)學(xué)生修改論文面授貴州師范學(xué)院數(shù)學(xué)與計(jì)算機(jī)科學(xué)學(xué)院本科畢業(yè)論文（設(shè)計(jì)）學(xué)號(hào) 1006020540200姓名論文（設(shè)計(jì)）題冃交叉評(píng)閱表學(xué)院（蓋章）：潘曉祥專業(yè)幾種常用數(shù)值積分方法的比較數(shù)學(xué)與應(yīng)用數(shù)學(xué)班級(jí)四班評(píng)語(yǔ)：該冋學(xué)在論文撰寫(xiě)過(guò)程中對(duì)相矢文獻(xiàn)閱讀范圍廣泛，方法正確 >內(nèi)容完整-能綜合運(yùn)用所學(xué)知識(shí)分析和解決實(shí)際

13、冋題-畢業(yè)論文撰與過(guò)程中態(tài)度端正勒奮刻苦。論文硏究了Newt on- Cotes求枳公式 ' 復(fù)化求積公式、 Romberg積分、高斯積分方法-通過(guò)算例分析，得出幾種常用數(shù)值積分方 法是解決頭際計(jì)算冋題的重要方法°論文結(jié)構(gòu)合理' 符合邏輯文章層次分明*語(yǔ)言準(zhǔn)確文字通順達(dá)到規(guī)范性要求*建議作為學(xué)士論文答辯。（滿分100分）指導(dǎo)教師（簽名）：評(píng)語(yǔ)：該冋學(xué)具備較好的基礎(chǔ)理論與專業(yè)知識(shí) > 學(xué)習(xí)態(tài)度認(rèn)真 > 閱讀教師指定的參哮資料、文獻(xiàn)較好的完成了任務(wù)書(shū) 規(guī)定的工作量。論文硏究了 Newton- Cotes求積公式、復(fù)化求積公式、 Romberg積分、咼斯積分方法

14、-通過(guò)算 例分析'得出幾種常用數(shù)值積分方法是解決實(shí)際計(jì)算問(wèn)題的重要方法。論文結(jié)構(gòu)合理'符合邏輯文章層次分明*語(yǔ)言準(zhǔn)確文字通順達(dá)到本科畢業(yè)論文相矢要求。同意參加答辯。（滿分100分）評(píng)閱教師（簽名）：論文題目作者姓名所屬學(xué)院貴州師范學(xué)院本科畢業(yè)論文答辯記錄表幾種常用的多項(xiàng)式插值方法潘曉祥學(xué)號(hào)1006020540200年級(jí)2010 級(jí)數(shù)計(jì)專業(yè)數(shù)學(xué)與應(yīng)用數(shù)學(xué)班級(jí)本科（四）班雍進(jìn)軍講師指導(dǎo)教肺姓名' 幟稱答辯會(huì)紀(jì)要時(shí)間2014年5月11日地點(diǎn)致遠(yuǎn)樓406答辯小組成負(fù)姓名職務(wù)（職稱）姓名職務(wù)（職稱）姓名職務(wù)（職稱）左羽教授崔忠偉副教授麼玉梅講師答辯中提出的主要問(wèn)題及冋答的簡(jiǎn)要情況

15、記錄：1自己做的有哪些？答：第頁(yè)至第12頁(yè)，總結(jié)進(jìn)行比較。2程序運(yùn)行過(guò)沒(méi)有？ 答：運(yùn)行過(guò)。320頁(yè)程序代碼中，if后的是什么符號(hào)？ 答：連接作用的符號(hào)。4-解釋一下什么時(shí)候用分號(hào)-什么時(shí)候不用？ 答：回答不清。5 摘要中央文拼與有錯(cuò) 答辯后修改答辯小組負(fù)責(zé)人簽名：左羽記錄人簽名：梅林林2014年5月們?nèi)赵u(píng)語(yǔ)：該生能在規(guī)定時(shí)間敘述論文的主要內(nèi)容-對(duì)提出的問(wèn)題一般能回答-無(wú)原則錯(cuò)誤。答辯小組經(jīng)過(guò)充分討論- 根據(jù)該生論文質(zhì)S和答辯中的表現(xiàn)*冋意評(píng)定論文成績(jī)?yōu)橹械取?。答、丄辯評(píng)定成績(jī)：77負(fù)責(zé)人（簽名）：左羽2014年5月11日1Abstract1 前言2數(shù)值積分方法的基本思想3 幾類常用數(shù)值積分方法

16、的簡(jiǎn)單分析3.1 Newt on Cotes 求積公式3.2 復(fù)化求積公式3. 3 Romberg求積公式34高斯型求積公式4幾類數(shù)值積分方法的簡(jiǎn)單比較評(píng)述5 利用MATLAB®程應(yīng)用對(duì)幾類求積算法的分析比較10結(jié)束語(yǔ)錯(cuò)誤！未定義書(shū)簽。致謝14附錄16貴卅師范學(xué)院畢業(yè)論文（設(shè)計(jì)）我們?cè)谇蠛瘮?shù)的積分時(shí)'往往因?yàn)樵瘮?shù)非常復(fù)雜以至于難以求出或用初等函數(shù)表 示' 這讓我們計(jì)算起來(lái)非常困難' 所以我們只能想辦法求它的近 似值，因此直接借助牛頓 “菜布尼茲公式計(jì)算定積分的，睜況是非常少見(jiàn)的。這時(shí)候數(shù)值積分就是解決這種問(wèn)題的一 種很好很有效的方法。本文從數(shù)值積分問(wèn)題的產(chǎn)生

17、出發(fā)，詳細(xì)介紹了一些數(shù)值積分的常用 方法（Newton Cotes求積公式，復(fù)化求積公式，Romberg求積公式高斯型求積公式）并對(duì)其進(jìn)行了簡(jiǎn)要的分析，在探討了這些數(shù)值積分算法的優(yōu)缺點(diǎn)的理論之外，我們 還將這些數(shù)值積分算法在計(jì)算機(jī)上通過(guò)matbb軟件編程實(shí)現(xiàn)應(yīng)用，并分別用各自求積 公式進(jìn)行運(yùn)算'以此來(lái)分析比較各種求積公式的代數(shù)精度和計(jì)算誤差。矢鍵詞：數(shù)值積分；求積公式；代數(shù)精度貴卅師范學(xué)院畢業(yè)論文（設(shè)計(jì)）Abstractfunction is very complex that it is difficult to find the elementary functions, whic

18、h makes u We in the function for the integration, often because the original s very difficult to calculate, so we can only think of a way to find the approximate value, thus directly with Newt on - Leib niz formula calculati ng defi nite in tegral situation is very rare. When numerical integration i

19、s to solve this problem in a very effective method. From the numerical integration problem, introduces some methods of nu merical in tegrati on (Newt on - Cotes quadrature formula, composite quadrature formulas, Lon gbei lattice quadrature formula, Gauss type quadrature formulas) and has carried on

20、brief analysis, discusses the advantages and disadvantages of these numerical integration algorithm theory, we will these numerical integration algorithm in the computer by MATLAB software program ming applicati on, and separately with their res pective quadrature formula for com put in g, in order

21、to an alyze the algebraic calculati on precisi on and error comparis on of various quadrature formulas.Keywords: Numerical in tegrati on; Calculati on meth; nu merical an alysis#貴卅師范學(xué)院畢業(yè)論文(設(shè)計(jì)丄、八亠亠1冃I言微積分的發(fā)明是世界數(shù)學(xué)史上一項(xiàng)輝煌的成就。但在實(shí)際求積問(wèn)題的時(shí)候'求解積分卻有著非常多局限性。比如對(duì)于定積分Ff(x)dx衽求某函數(shù)的定3積分時(shí)，在一定條件下，雖然有牛頓萊布里茨公式1=叫X)

22、dx二F(b) 一 F(a)可弋a(chǎn)以計(jì)算定積分的值，但在很多情況下f(x)的原函數(shù)不易求出或非常復(fù)雜。被()=沁(，等；有的X積函數(shù)f(x)的原函數(shù)很難用初等函數(shù)表達(dá)出來(lái) ' 例如函數(shù)f(x)的原函數(shù)F(x)存在，但其表達(dá)式太復(fù)雜，計(jì)算量太大，有的甚至無(wú)法 有解析表達(dá)式。因此能夠借助牛頓萊布尼茲公式計(jì)算定積分的情形是不多的。另外，許多實(shí)際問(wèn)題中的被積函數(shù)往往是列表函數(shù)或其他形式的非連續(xù)函數(shù)對(duì)這類函 數(shù)的定積分，也不能用不定積分方法求解，只能設(shè)法求其近似值。因此，探討近似計(jì)算的 數(shù)值積分方法是有明顯的實(shí)際意義的即有必要研究定積分的數(shù)值計(jì)算方法，以解決定積分 的近似計(jì)算。而數(shù)值積分就是解決

23、此類問(wèn)題的一種有效的方法，它的特點(diǎn)是利用被積函數(shù)f M在一些節(jié)點(diǎn)上的信息求出定積分的近似值。在很多實(shí)際應(yīng)用中，只能知道積分函數(shù)在某些特定點(diǎn)的取值比如天氣測(cè)量中的氣溫'濕度'氣壓等，醫(yī)學(xué)測(cè)量中的血壓'濃度等等。通過(guò)研究，我們將會(huì)更熟練掌握一些 數(shù)值積分方法去計(jì)算一些特定條件的數(shù)值計(jì)算，以便 我們得到自己想要的結(jié)果。2數(shù)值積分方法的基本思想在數(shù)學(xué)分析中，計(jì)算連續(xù)函數(shù)f(x)在區(qū)間a,b上的積分是通過(guò)f(x)的原函數(shù)F(x),由下列定積分公式a f(x)dx 二 F(b)F (a)得到的。但由于大量被積函數(shù)的原函數(shù)不能用初等函數(shù)表示，因此，很難用求原函數(shù)的公式彳(x)dx =

24、 F(b)-F(a)得到積分；有些被積函數(shù)f(x)不是明顯知道的，例如由數(shù)值表給出它的離散值，或者是它被定義為某個(gè)微分方程的解，而這個(gè)微分方程是不能顯示解出的。這說(shuō)明，按f(x)dx = F(b) - F(a)公式計(jì)算定積分是有很大局限性的。因而常常采用在電子計(jì)算機(jī)上很有效的數(shù)值積分方法°我們從定積分的定義b n.(,f (x)dx :(b -a)' C(n.k) f(Xk)k=0出發(fā)。推導(dǎo)出兩個(gè)簡(jiǎn)單的數(shù)值積分公式。f(x)dx =沖(kE)AXk式的幾何意義，kn就是把整塊曲線梯形的面積積分成若干個(gè)小曲邊梯形面積的和，當(dāng)無(wú)限細(xì)分時(shí)這個(gè)和取極 限就是真正曲邊梯形面積。去掉取極

25、限這一步，用有限個(gè)小曲邊梯形面積的和，代替整塊 的曲邊梯形面積'從而求得一個(gè)近似值，這就是數(shù)值積分的基本思想。根據(jù)小區(qū)間的不同分割方法和各分點(diǎn)f( Z )值的不同選擇'就得到不同的數(shù)值積分公式。求和數(shù)值求積公式是取lab 1上若干個(gè)點(diǎn)Xk處的高度f(wàn)(XK)，通過(guò)加權(quán)&后，再、Akf(Xk)k=0從而得到積分的近似值。數(shù)值求積公式寫(xiě)成一般形式f(x)dx 八 AJ(Xk)k=0式中Xk稱求積節(jié)點(diǎn)，Ak稱求積系數(shù)，也稱伴隨節(jié)點(diǎn)Xk的權(quán)。當(dāng)積分區(qū)間a,bl確 定后 ' 求積系數(shù)人僅僅與節(jié)點(diǎn)Xk的選取有矢，而不依賴被積函數(shù)f(x)的具體貴卅師范學(xué)院畢業(yè)論文(設(shè)計(jì)Rjf

26、f (x)dxA f (Xk)k=0#貴州師范學(xué)院畢業(yè)論文（設(shè)計(jì)）把R If 1稱為求積公式的截?cái)嗾`差或余項(xiàng)。數(shù)值求積方法的特點(diǎn)是直接利用積分區(qū)間lab 1上一些離散節(jié)點(diǎn)的函數(shù)值 進(jìn)行線性組合來(lái)近似計(jì)算定積分的值 ' 從而將定積分的計(jì)算歸結(jié)為函數(shù)值的 計(jì)算，這就避開(kāi)了牛頓 菜布尼茲公式需要尋求原函數(shù)的困難，并為計(jì)算機(jī)求積分提供了可行性。3幾類常用數(shù)值積分方法的簡(jiǎn)單分析3.1 Newton- Cotes求積公式常用的梯形公式和Simpsor公式是低階的牛頓柯特斯公式，牛頓柯特斯 公式是積分區(qū)間上等距節(jié)點(diǎn)的插值求積公式。插值求積公式在積分區(qū)間上， 所取節(jié)點(diǎn)是等距時(shí)稱為牛頓柯特斯公式 

27、9; 即f (x)dx : (b-a)- C(n,k)f(Xk)k=0其中c(n,k)為Cotes求積公式的系數(shù)，是n和k的函數(shù)。當(dāng)時(shí)，為梯形公式：:f(x)dx © 9f(a)2f(b)7梯形公式的代數(shù)精度為1，有兩個(gè)積分節(jié)點(diǎn)。當(dāng)n=2時(shí) ' 為Simpson公式:b(b -a)/(b)af(x)cJx f(a) 4f(6Simpson公式的代數(shù)精度為3，有三個(gè)積分節(jié)點(diǎn)由于只增加一個(gè)節(jié)點(diǎn)，其代數(shù)精度增加2，由此可知，Simpson公式比梯形公式代當(dāng)n=4時(shí)' NewtonCotes求積公式為Cotes公式：:f(x)dx 烤7f(a) 32f (專)12f(寧)32

28、(寧)7f(b)Newton-Cotes公式的代數(shù)精度為5,有5個(gè)積分節(jié)點(diǎn)。所以對(duì)于Newt omCotes積分公式，n為偶數(shù)時(shí)的代數(shù)精度要比n為奇數(shù)時(shí)的積分公式效果比較優(yōu)越。但并不是n的值越大越好，當(dāng)n過(guò)大時(shí)（n=8），求積公式的數(shù)值穩(wěn)定性不好。3.2復(fù)化求積公式由于NewtomCotes的節(jié)點(diǎn)n越大對(duì)應(yīng)的精度就越高，但是n=8時(shí)公式的數(shù)值 是 不穩(wěn)定的'因此就不能用增加求積節(jié)點(diǎn)的方法來(lái)提高精度，因此，我們常常將求積區(qū)間a,b分成若干小區(qū)間，然后在每個(gè)小區(qū)間上采用數(shù)值穩(wěn)定的Cotes公式求小區(qū)間上的積分'然后把每個(gè)小區(qū)間上的結(jié)果加起來(lái)作為原定積 分的近視值，這種方法構(gòu)造的求積

29、公 式就叫做復(fù)化求積公式。常用的復(fù)化求積公式有: 復(fù)化梯形公式：n變步長(zhǎng)梯形公式為：T2計(jì) f（xo）：r）二豬Mn 二寧"f(X2ii)2n j#復(fù)化Simpson公式:b - a人盲“心)ygb_aS2n 八 Ik6n變步長(zhǎng)復(fù)化Simpson公式:"心f(X。)fgn) 4、f(X2 一 2' f(X2k)kdkT3.3 Romberg求積公式Romberg積分方法也叫做逐次分半加速法，它是在復(fù)化梯形公式誤差估計(jì)的基礎(chǔ)上'應(yīng)用線性外推的方法構(gòu)造出的一種加速算法。將積分區(qū)間分成n等分和2n等分時(shí)，求得積分近似值Tn和Tm，并沒(méi)有誤差估計(jì)式1I 'T

30、an. 3 （Tzn-Tn）3積分近似值Tzn的誤差大致等于3仃2人），當(dāng)用-）對(duì)丁?.進(jìn)行修正時(shí),33-<T2n-Tn）與n之和比T?n更接近于真值I，故Pn）是對(duì)Tth誤差的一種補(bǔ)償, 33因此可以期望下式是一個(gè)更好的結(jié)果，即T - T2n' (J -)二今 n2- ； Tn3 33F面說(shuō)明T即是分成n等分時(shí)Simpson公式的值&。將復(fù)化梯形公式nVf(a) 2 ' f(xO f(b)梯形變步長(zhǎng)求積公式T2ndin4f(X!2 2心n4n4代入上式T表達(dá)式得-hT= ; Jf(a)+4 送 f(x .)+2 送 f(xj+f(b)6 -心 遷心這就是說(shuō) &#

31、39; 用梯形法二分前后兩個(gè)梯形值Tn和丁加作線性外推，結(jié)果得到Simpson法的積分值&。將誤差由o（»）變?yōu)閛（hj，從而提高了逼近精度。再考察Simpson法。其截?cái)嗾`差與成正比，因此，若將步長(zhǎng)折半，則誤差減至】，即有 16:16 S2 j1515不難驗(yàn)證 ' 上式右端的值其實(shí)等于Cn，就是說(shuō)，用Simpson法二分前后的兩個(gè)積分值& 與S2n,按上式再作線性外推，結(jié)果得到柯特斯法的積分值即有16 1Cn S2n Sn1515這時(shí)將誤差由0（椚變?yōu)閛（h。），逼近精度又一次得以提高。同樣的方法'依據(jù)柯特斯法的誤差公式，可進(jìn)一步導(dǎo)出下列龍貝格公式6

32、41% C2nCn6363Rn逼近積分值的誤差為0（"，這樣Romberg公式將誤差由o（h。）變?yōu)閛（h），逼近精度再次得以提高。Romberg公式有7次代數(shù)精度，這表明該公式不是牛 頓柯特斯公 式。在步長(zhǎng)二分的過(guò)程中運(yùn)用SnG、Rn表達(dá)式加工三次，就能將粗糙的積 分值九逐步 加工成精度較高的Romberg值Rn，或者說(shuō)，將收斂緩慢的梯形值序列九加工成收斂迅 速的Romberg值序列Rn，這種加速方法稱Romberg算法。3.4高斯型求積公式前面介紹的n1個(gè)節(jié)點(diǎn)的Newton -Cotes求積公式，其特征是節(jié)點(diǎn)是等距 的。這 種特點(diǎn)使得求積公式便于構(gòu)造'復(fù)化求積公式易于形成

33、。但同時(shí)也限 制了公式的精度。n 是偶數(shù)時(shí) ' 代數(shù)精度為' n是奇數(shù)時(shí)，代數(shù)精度為n ；我們知道n個(gè)節(jié)點(diǎn)的插值型求n 0能不能在區(qū)間積公式的代數(shù)精確度不低于a,b 1上適當(dāng)選擇個(gè)節(jié)點(diǎn)Xo,Xi,X2 X.使插值求積公式的代數(shù)精度高于n呢?答案是肯定的'適當(dāng)選擇節(jié)點(diǎn)'可使公式的精度最高達(dá)到不失一般性，將求積公式Ja所學(xué)的高斯型求積公式°f（x）dx八Akf（Xk）的求積區(qū)間La,b轉(zhuǎn)換成TJk=0貴州師范學(xué)院畢業(yè)論文（設(shè)計(jì)）對(duì)任意求積區(qū)間ab 1作變換D-a, 322Xtt可以變換到區(qū)間1-1,1上，這時(shí)af(X) dX=乜"b-aT&quo

34、t;ba丄 ab' *b'a丄 ab*P(附'J (t) dt 二 2I 22b V(t) dt其中仕）二f （ b合tbt） o22高斯勒讓德求積公式在這里簡(jiǎn)稱高斯公式，它是在區(qū)間1-1,1上進(jìn)行討論 的。4幾類數(shù)值積分方法的簡(jiǎn)單比較評(píng)述所以不同類由于我們?cè)谟?jì)算實(shí)際問(wèn)題是往往要考慮到代數(shù)精度和計(jì)算最， 型的求積公式有著不同的特點(diǎn)：Simpson積分方法和梯形積分方法雖然計(jì)算簡(jiǎn)便，但是精度比較差，不理想。但對(duì)于光滑性較差的被積函數(shù)有時(shí)會(huì)比高精度的積分方法更為有效。特別是梯形積分方法對(duì)被 積函數(shù)是周期函數(shù)的求積效果更為突出。Newt on- Cotes公式是不穩(wěn)定的 &#

35、39; 然而復(fù)化梯形公式和復(fù)化Simpson公式不僅 保留 了低階公式的優(yōu)點(diǎn)還能夠獲得比較較高的精度，所以在實(shí)際計(jì)算中應(yīng)用 得最為廣泛。Romberg積分方法的算法簡(jiǎn)單.方便編程的實(shí)現(xiàn)。收斂速度快-計(jì)算精 度較高，但 是計(jì)算;較大。Gauss積分方法的精度較高，數(shù)值穩(wěn)定收斂速度較快，但因?yàn)槠涔?jié)點(diǎn)不規(guī)則'計(jì)算比較麻煩5利用MATLAB編程應(yīng)用對(duì)幾類求積算法的分析比較在簡(jiǎn)單的認(rèn)識(shí)積分方法比且理論比較之后，則要進(jìn)行數(shù)學(xué)實(shí)驗(yàn)進(jìn)行驗(yàn)證，因此就要通過(guò)matlab軟件對(duì)各種積分方法逬行編程并運(yùn)算，然后對(duì)其各種方法的運(yùn)算結(jié)果進(jìn)行分析比較'掌握和理解各方法的優(yōu)缺點(diǎn)。規(guī)定各個(gè)程序都以I = 0皿d

36、x為例子進(jìn)行運(yùn)算。原積分的精確值為isin XoX dx"9°®M"70367483例 分別用不同的方法計(jì)算積分I二，沁dx，并作比較。用以上介紹的幾類積分方法分別計(jì)算積分，得出誤差，并進(jìn)行比較:1 ' 用 Newt on-Cotes 公式當(dāng)時(shí)，即用梯形公式，用程序一（程序見(jiàn)附錄）11在MATLAB命令窗口中輸入>>NCotes （0,1,1,2）得92703549240395 Ir =0.01904757796323當(dāng)n=2時(shí)，即用Simpson公式，用程序一（程序見(jiàn)附錄）在MATLAB命令窗口中輸入>>NCotes （

37、0,1,2,2）得I 肝 0.94614588227359 Ir =0.000062811906407當(dāng)n=4時(shí)，即用科特斯公式，用程序一（程序見(jiàn)附錄）在MATLAB命令窗口中輸入>>NCotes （0,1,4,2）得I -0.94608300406367R =0.0000000663035132-用復(fù)化梯形公式令 h=1/8=0.125,ffl程序二（程序見(jiàn)附錄）在MATLAB命令窗口中輸入 » trapri Cf,0,1.8），得J ；叱0%(f (0) +2 f (h) +. .+ f(7h)+ f(l) =0.9456908635270 x2R =0.000392

38、2068401823 '用復(fù)化Simpson公式令h=1/8=0J25,用程序三（程序見(jiàn)附錄）怡nXdx： hf(0)4f(h) °x3=0.94608308538495R =0.000000015017767在MATLAB命令窗口中輸人>> simpr1（t0,8）,得f(7h)b： ； , 2f(2h)f(6h)IHf(1£4 '用Romberg公式用程序四（程序見(jiàn)附錄）在MATLAB命令窗口中輸人» romber(T,0,1,5,0.5*(10 八(8)，得(沁 dx 4.94608307036718X1R = 0.0000000

39、000000025 -用高斯勒讓德求積公式令 x=(t+1)/2，| =t+1（1）用2個(gè)節(jié)點(diǎn)的Gauss公式10.94604115827633（2）用3個(gè)節(jié)點(diǎn)的Gauss公式，用程序五（程序見(jiàn)附錄）在 MATLAB 命令窗口中輸人 » GuassLegendre (0,1,2,2)得1-0.946083134078473|R =0.000000063711290算法比較:1.原積分的精確值為：1 sin X/ = dx =0.946083070367183 oX2.由例題的各種求積算法可知：（1）對(duì)Newton-cotes公式，當(dāng)nh時(shí)只有1位有效數(shù)字，當(dāng)n=2時(shí)有3位有效數(shù)貴州師

40、范學(xué)院畢業(yè)論文(設(shè)計(jì))字，當(dāng)n=4時(shí)有7位有效數(shù)字。(4) 用復(fù)化梯形公式有2位有效數(shù)字 ' 對(duì)復(fù)化Simpson公式有7位有效數(shù)字。(5) 用復(fù)化梯形公式 ' 對(duì)積分區(qū)間0,1二分"次用了 2049個(gè)函數(shù)值，才可以得7位有效數(shù)字°(6) 用Romberg公式對(duì)區(qū)間0們二分3次用了 9個(gè)函數(shù)值，就可以得到7位有效數(shù)字；二分4次用了 14個(gè)函數(shù)值，卻可以得到14位有效數(shù)字。(7) 用高斯勒讓德求積公式僅僅用了 3個(gè)函數(shù)值，就能得到比較精確的6位有效數(shù)字。#貴州師范學(xué)院畢業(yè)論文（設(shè)計(jì)）結(jié)束語(yǔ)本文主要研究了常用的幾類數(shù)值積分的求積算法并通過(guò)例題計(jì)算積分進(jìn)行分析比

41、較。Newt on-Cotes積分方法是一種非常普遍的積分方法，然而梯形積分方法的誤差最大，近似效果最差，Simpson積分方法的精度比梯形積分方法高了一個(gè)數(shù)量級(jí)；Cotes 積分方法精度比Simpson積分方法高兩個(gè)數(shù)量級(jí)。則Cotes代數(shù)精度比較高。由此可知 一般情況下，積分公式代數(shù)精度越高，計(jì)算精度也越高。但是高階的Cotes積分方法收Simpson 積分?jǐn)啃詻](méi)有保證，因此實(shí)際應(yīng)用中很少用。復(fù)化梯形積分方法比梯形積分方法精度高，同樣的'復(fù)化 方法比Simpson積分方法精度高，高了差不多7個(gè)數(shù)最級(jí)，所以復(fù)化積分方法比較優(yōu) 越。Romberg積分方法收斂速度快-計(jì)算精度較高，但是計(jì)

42、算量較大。Gauss積分方法精度高 > 數(shù)值穩(wěn)定收斂速度較快，但是計(jì)算麻煩。經(jīng)研究可以知道Newt omCotes方法的代數(shù)精度越高，數(shù)值積分的效果越 好 ' 越精 確。當(dāng)積分區(qū)間比較大的時(shí)候，積分?jǐn)?shù)值不穩(wěn)定 ' 這個(gè)時(shí)候可以利用復(fù)化積分方法效果會(huì) 更好；Romberg積分方法可以利用變步長(zhǎng)復(fù)化積分公式得到更為精確的數(shù)值結(jié)果，是比 較好的積分方法。高斯求積方法精確度高，收斂性快，比其他積分方法優(yōu)越。具有很廣泛13的運(yùn)用。參考文獻(xiàn)1（第二版）32005.3薛毅.數(shù)學(xué)分析與實(shí)驗(yàn)M北京:業(yè)大學(xué)出版社4徐士良數(shù)值分析與算法M 北京:械工業(yè)出版社2007.15王開(kāi)榮.楊大地.應(yīng)用數(shù)

43、值分析M北京:等教育出版社2010.76楊一都數(shù)值計(jì)算方法M.北京:等教育出版社.2008.4韓明王家寶李林?jǐn)?shù)學(xué)實(shí)驗(yàn)（MATLAB）版M.上海:濟(jì)大學(xué)出版社2012.18圣寶建.矣于數(shù)值積分若干問(wèn)題的研究J 南京信息工程大學(xué).2009.05.01. :429劉緒軍.幾種求積公式計(jì)算精確度的比較J.南京職業(yè)技術(shù)學(xué)院.2009.10史萬(wàn)明.吳裕樹(shù).孫新.數(shù)值分析M.北京理工大學(xué)出版tt.2010.4.矢治陸金甫數(shù)學(xué)分析基礎(chǔ)M.北京：等教育出版社.2010.7胡祖熾林源渠數(shù)值分析M北京:等教育出版社.1986.3貴卅師范學(xué)院畢業(yè)論文（設(shè)計(jì)行文至此，我的這篇論文已接近尾聲：歲月如梭，我四年的大學(xué)時(shí)光也

44、 即將敲 響結(jié)束的鐘聲。離別在即，站在人生的又一個(gè)轉(zhuǎn)折點(diǎn)上，心中難免思緒萬(wàn)千' 一種感 恩之情油然而生。首先感謝貴州師范學(xué)院四年來(lái)對(duì)我的培養(yǎng)，是博學(xué)的老師們教會(huì)了我學(xué)習(xí)的方 法、鍛煉了我思考的能力指明了我未來(lái)奮斗的方向，從而使我進(jìn)一步明確了人生的 目標(biāo)。其次，我要感謝我的指導(dǎo)老師一雍進(jìn)軍老師，他的嚴(yán)謹(jǐn)細(xì)致' 一絲不茍的作風(fēng)一直是 我工作、學(xué)習(xí)中的榜樣；他的循循善誘的教導(dǎo)和不拘一格的思 路給予我無(wú)盡的啟迪。在撰 寫(xiě)整個(gè)畢業(yè)論文的過(guò)程當(dāng)中 ' 他為我們考慮到了每一個(gè)細(xì)節(jié)，從開(kāi)題報(bào)告到畢業(yè)論文的擬 定修改上，雍老師更是不厭其煩的 為我們做好每一步的細(xì)心指導(dǎo)。對(duì)此，我表示衷心地

45、感 謝。沒(méi)有雍老師' 我的論文也不可能這么順利的完成。同時(shí)，我也要感謝每一位給過(guò)我?guī)椭?的老師和同學(xué)，在我撰寫(xiě)論文的過(guò)程當(dāng)中同樣給了我大量有益的建議，在此一并向他們表 示真誠(chéng)的感謝一感謝他們對(duì)我的支持和幫助。最后感謝這篇論文所 涉及到的各位學(xué)者' 本 文引用了數(shù)位學(xué)者的研究文獻(xiàn)' 如果沒(méi)有各位學(xué)者的 研究成果帶給我的的幫助和啟發(fā)，我 將很難完成本篇論文的寫(xiě)作。由于我的學(xué)術(shù)水平有限，所寫(xiě)論文難免有不足之處，懇請(qǐng)各位老師和學(xué) 友批評(píng)指 正。最后，衷心感謝評(píng)閱論文及參加答辯的各位老師！15貴州師范學(xué)院畢業(yè)論文(設(shè)計(jì))1 NewtonCotes求積公式的MATLAB實(shí)現(xiàn)先用M文

46、件定義一個(gè)名為fl.m的函數(shù):% i是要調(diào)用第幾個(gè)被積函數(shù)g(i)，X是自變最 g(1)=sqrt(x);if x=0g(2)=1： else g(2)=si n(x)/x;end g(3)=4/(1+x2);匸 g(i);程序一：fun cti on C,g=NCotes(a,b, n,m)%a， b分別為積分的上下限；%n是子區(qū)間的個(gè)數(shù);% m是調(diào)用上面第幾個(gè)被積函數(shù)；%當(dāng)n=1時(shí)計(jì)算梯形公式；當(dāng)n=2時(shí)計(jì)算辛浦生公式，以此類推;i=n：h=(b-a)/i;貴卅師范學(xué)院畢業(yè)論文（設(shè)計(jì)19z=0;forj=0:ixG+1)=a+j*h; s=1;ifj=Os=s;elsefor k=1:js

47、=s*k;endendif i-j=0r=r：elser=r*k;endendif mod(i-j),2)=1elseq=i*s*r;end貴州師范學(xué)院畢業(yè)論文(設(shè)計(jì))y=i： for k=0:iif k =jy=y*(sym(T)-k);end end l=int(y,0J);CG+1)=l/q；2=z+CG+1)*f1(m,xG+1): endg=(b-a)*21.11) 當(dāng)輸入, n=1=2時(shí)，即在MATLAB命令窗口輸入 » NCotes(0,1,1,2)即可得用梯形公式的積分值和相應(yīng)科特斯系數(shù)如圖2) 當(dāng)輸入a=0,b=1 / n =2,m=2時(shí)即在MATLAB命令窗口輸入

48、 » NCotes(0,1,2,2)即可得用辛浦生公式的積分值和相應(yīng)科特斯系數(shù)如圖1.23) 當(dāng)輸入a=0,b=1 / n =4,m=2時(shí)即在MATLAB命令窗口輸入 » NCotes(0,1,4,2)即可得用科特斯公式的積分值和相應(yīng)科特斯系數(shù)如圖1.3#貴州師范學(xué)院畢業(yè)論文（設(shè)計(jì)）23圖"圖1.2圖1.32復(fù)化梯形求積公式的MATLAB實(shí)現(xiàn)通過(guò)f（X）的n 1個(gè)等步長(zhǎng)節(jié)點(diǎn)逼近積分：/ (幼辦怎 £(/()+/(*) )+方 £ yl匕)其中,xP kh,Xo=a程序二:fun cti on s=trapM (f,a,b, n)% f是被積函數(shù)

49、;%a,b分別為積分的上下限;% n是子區(qū)間的個(gè)數(shù)；%S是梯形總面積；h=(b-a)/n;s=0;for k=1:( n-1)x=a+h*k;s=s+feval(I,x);endformat longs=h*(feval(T,a)+feval(T,b)/2+h*s;先用M文件定義一個(gè)名為f.m的函數(shù):fun cti on y=f(x)if x=0y=1：elsey=s in( x)/x;end若取子區(qū)間的個(gè)數(shù) 在MATLAB命令窗口中輸入»trapr1(t0,1.8)回車(chē)得到如圖2.1貴卅師范學(xué)院畢業(yè)論文（設(shè)計(jì)）圖2.13復(fù)化Simpson求積公式的MATLAB實(shí)現(xiàn)程序三：fun c

50、ti on s=simpr1(f,a,b ,n)% f是被積函數(shù);% a,b分別為積分的上下限；%n是子區(qū)間的個(gè)數(shù)；%s是梯形總面積，即所求積分?jǐn)?shù)值；h=(b-a)/(2* n);s仁0;s2=0;for k=1: nx=a+h*(2*k-1);si 二 s1+feval(T,x);end25貴州岀碑堂聲f禪耳丫設(shè)計(jì)）s=0;27endx=a+h*2*k;s2=s2+feval(T,x);s 二 hWal.a)+feval(T,b)+4*s1+2*s2)/3：先用M文件定義一個(gè)名為f.m的函數(shù):fun cti on y=f(x)if x=0y=1：elsey=s in( x)/x;end若取子