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

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

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

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

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

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

6、復(fù)化求積公式2.3 Romberg求積公式2.4 高斯型求積公式3幾類數(shù)值積分方法的簡單比較評述4利用MATLAB編程應(yīng)用對幾類求積算法的分析比較研究方法：本論文主要通過對相關(guān)文獻(xiàn)和書籍的參考，合自己的見解，紀(jì)化求積公式，Newton-Cotes求積公式Romberg求積公式，高斯型求積公式進(jìn)行討論并進(jìn)行上機(jī)實(shí)驗(yàn)，從代數(shù)精度，求積公式誤差等角度對這些方法進(jìn)行分析比較完成期限和采取的主要措施：本論文計劃用6個月的時間完成，階段的任務(wù)如下：(1) 7月份查閱相關(guān)書籍和文獻(xiàn)；(2) 8月份完成開題報告并交老師批閱；(3) 9月份完成論文初稿并交老師批閱；(4) 10月份完成論文二搞并交老師批閱；(5

7、) 11月份完成論文三稿；(6) 12月份定稿.主要措施：考相關(guān)書籍和文獻(xiàn)，合自己的見解，老師的指導(dǎo)下和同學(xué)的幫助下完成主要參考文獻(xiàn)及資料名稱：1吳治陸金甫數(shù)學(xué)分析基礎(chǔ)(第二版)M.北京：等教育出版社.2010.72胡祖熾.林源渠.數(shù)值分析M北京：等教育出版社.1986.33薛毅.數(shù)學(xué)分析與實(shí)驗(yàn)M北京：業(yè)大學(xué)出版社2005.34徐士良.數(shù)值分析與算法M.北京：械工業(yè)出版社2007.15王開榮.楊大地.應(yīng)用數(shù)值分析M北京:等教育出版社2010.76楊一都.數(shù)值計算方法M.北京:等教育出版社.2008.47韓明.王家寶.李林.數(shù)學(xué)實(shí)驗(yàn)(MATLAB版M.上海：濟(jì)大學(xué)出版社2012.18圣寶建.關(guān)于

8、數(shù)值積分若干問題的研究J.南京信息工程大學(xué).2009.05.01:429劉緒軍.幾種求積公式計算精確度的比較J.南京職業(yè)技術(shù)學(xué)院.2009.10史萬明.吳裕樹.孫新.數(shù)值分析M.北京理工大學(xué)出版社.2010.4.指導(dǎo)教師苣見:簽名：開題報告會紀(jì)要時間2013年8月26日地點(diǎn)寧靜樓229教師辦公室與會人員姓名職務(wù)（職稱）姓名職務(wù)（職稱）姓名職務(wù)（職稱）蕤進(jìn)軍導(dǎo)師（講師）鄧喜才副教授李晟副教授組長會議記錄摘要：指導(dǎo)小組針對課題二次函數(shù)性質(zhì)的應(yīng)用提問了以下問題以及報告人的回答：雍老師問：擾擇此題目的目的？潘曉祥答：隨著計算機(jī)和計算方法的飛速發(fā)展，幾乎所有學(xué)科都走向定量化和精確化，計算數(shù)學(xué)中的數(shù)值計算

9、方法則是解決“計算”問題的橋梁和工具。鄧?yán)蠋熀危簩@個問題進(jìn)行研究有什么實(shí)際的意義？潘曉祥答：計算方法既有數(shù)學(xué)類課程中理論上的抽象性和嚴(yán)謹(jǐn)性，又有實(shí)用性和實(shí)驗(yàn)性的技術(shù)特征，計算方法是一門理論性和實(shí)踐性都很強(qiáng)的學(xué)科在科學(xué)研究和工程技術(shù)中都要用到各種計算方法.例如，在航天航空、地質(zhì)勘探、汽車制造、橋梁設(shè)計、天氣預(yù)報和漢字字樣設(shè)計中都有計算方法的蹤影。李老師何：對這個問題你有什么自己的看法？潘曉祥答：隨著計算機(jī)技術(shù)的迅速發(fā)展和普及，現(xiàn)在計算方法課程幾乎已成為所有理工科學(xué)生的必修課程我們知道，計算能力是計算工具和計算方法的效率的乘積，提高計算方法的效率與提高計算機(jī)硬件的效率向樣重要科學(xué)計算已用到科學(xué)技

10、術(shù)和社會生活的各個領(lǐng)域中所以，研究數(shù)值計算方法可以讓數(shù)學(xué)的應(yīng)用更大更廣。會議主持人簽名：記錄人簽名：在日口指導(dǎo)小組意見負(fù)責(zé)人簽名：年月日學(xué)院意見負(fù)責(zé)人簽名：年月日貴州師范學(xué)院數(shù)學(xué)與計算機(jī)科學(xué)學(xué)院指導(dǎo)教師指導(dǎo)本科畢業(yè)論文情況登記表論文（設(shè)計）題目幾種常用數(shù)值積分方法的比較用學(xué)生姓名潘曉祥學(xué)號1006020540200年級2010級所屬學(xué)院數(shù)計學(xué)院專業(yè)數(shù)學(xué)與應(yīng)用數(shù)學(xué)班級四班指導(dǎo)教師姓名雍進(jìn)軍職稱講師學(xué)歷碩士指導(dǎo)時間指導(dǎo)地點(diǎn)指導(dǎo)內(nèi)容指導(dǎo)教師簽名備注2013年06月10日致遠(yuǎn)樓416論文擾題，資料準(zhǔn)備面授2013年06月22日網(wǎng)上確定畢業(yè)論文擾題電子郵件2013年06月26日網(wǎng)上怎樣撰寫畢業(yè)論文開題報

11、告電子郵件2013年06月28日網(wǎng)上指導(dǎo)學(xué)生撰寫開題報告電子郵件2013年07月14日網(wǎng)上幫助學(xué)生查找有關(guān)參考文獻(xiàn)電子郵件2013年07月17日手機(jī)女M可構(gòu)思自己的畢業(yè)論文手機(jī)飛信2013年08月21日手機(jī)聽取學(xué)生畢業(yè)論文寫作進(jìn)展情況匯報手機(jī)飛信2013年08月28日網(wǎng)上解答學(xué)生在論文寫作中遇到的疑惑電子郵件2013年9月09日網(wǎng)上幫助學(xué)生查找有關(guān)參考文獻(xiàn)電子郵件2013年11月28日網(wǎng)上女M可規(guī)劃自己的論文電子郵件2013年12月04日手機(jī)怎樣寫好論文引言手機(jī)飛信2013年12月08日網(wǎng)上怎樣寫好論文引言電子郵件2013年12月12日網(wǎng)上怎樣寫論文摘要電子郵件2013年12月16日網(wǎng)上怎樣疑

12、取論文關(guān)鍵詞電子郵件2013年12月20日網(wǎng)上怎樣編輯論文中的公式電子郵件2014年01月05日手機(jī)督促學(xué)生在寒假中寫好論文的初稿電子郵件2014年02月27日寧靜樓219檢查學(xué)生論文完成情況面授2014年03月03日寧靜樓219對學(xué)生的論文初稿提出修改時意見面授2014年03月07日寧靜樓219解答學(xué)生在修改時的困惑面授2014年03月11日寧靜樓219指導(dǎo)學(xué)生修改論文面授貴州師范學(xué)院數(shù)學(xué)與計算機(jī)科學(xué)學(xué)院本科畢業(yè)論文（設(shè)計）交叉評閱表學(xué)院（蓋章）：學(xué)號1006020540200姓名潘曉祥專業(yè)數(shù)學(xué)與應(yīng)用數(shù)學(xué)論文（設(shè)計）題目幾種常用數(shù)值積分方法的比較班級四班指導(dǎo)教師意見評語：該向?qū)W在論文撰寫過程

13、中對相關(guān)文獻(xiàn)閱讀范圍廣泛，方法正確，內(nèi)容完整能綜合運(yùn)用所學(xué)知識分析和解決實(shí)際問題。畢業(yè)論文撰與過程中態(tài)度端正，勤奮刻苦。論文研究了Newton-Cotes求枳公式、復(fù)化求積公式、Romberg積分、高斯積分方法，通過算例分析，得出幾種常用數(shù)值積分方法是解決頭際計算問題的重要方法。論文結(jié)構(gòu)合理，符合邏輯，文章層次分明，語言準(zhǔn)確，文字通順達(dá)到規(guī)范性要求，建議作為學(xué)士論文答辯。成績：（滿分100分）指導(dǎo)教師（簽名）：年月日評閱教師意見評語：該向?qū)W具備較好的基礎(chǔ)理論與專業(yè)知識，學(xué)習(xí)態(tài)度認(rèn)真，閱讀教師指定的參考資料、文獻(xiàn)，較好的完成了任務(wù)書規(guī)定的工作量。論文研究了Newton-Cotes求積公式、復(fù)化求

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

15、的主要問題及回答的簡要情況記錄：1自己做的有哪些？答：第11頁至第12頁，總結(jié)進(jìn)行比較。2 .程序運(yùn)行過沒有？答：運(yùn)行過。3 .20頁程序代碼中，if后的是什么符號？答：連接作用的符號。4解釋一下什么時候用分號，什么時候不用？答：回答不清。5.摘要中央文拼與有錯答辯后修改答辯小組負(fù)責(zé)人簽名：左羽記錄人簽名：梅林林2014年5月11日答j_L辯小組意見評語：該生能在規(guī)定時間敘述論文的主要內(nèi)容，對提出的問題一般能回答，無原則錯誤。答辯小組經(jīng)過充分討論，根據(jù)該生論文質(zhì)量和答辯中的表現(xiàn)，向意評定論文成績?yōu)椤爸械取?。評定成績：77負(fù)責(zé)人（簽名）：左羽2014年5月11日Abstract21前言32數(shù)值積

16、分方法的基本思想33幾類常用數(shù)值積分方法的簡單分析51 .1NewtonCotes求積公式52 .2復(fù)化求積公式63 .3Romberg求積公式74 .4高斯型求積公式84幾類數(shù)值積分方法的簡單比較評述95利用MATLAB程應(yīng)用對幾類求積算法的分析比較10結(jié)束語錯誤！未定義書簽。致謝14附錄16貢州師范學(xué)院畢業(yè)論文（設(shè)計）我們在求函數(shù)的積分時，往往因?yàn)樵瘮?shù)非常復(fù)雜以至于難以求出或用初等函數(shù)表示，這讓我們計算起來非常困難，所以我們只能想辦法求它的近似值，因此直接借助牛頓一萊布尼茲公式計算定積分的情況是非常少見的。這時候數(shù)值積分就是解決這種問題的一種很好很有效的方法。本文從數(shù)值積分問題的產(chǎn)生出發(fā)

17、，詳細(xì)介紹了一些數(shù)值積分的常用方法（NewtonCotes求積公式，復(fù)化求積公式，Romberg求積公式，高斯型求積公式）并對其進(jìn)行了簡要的分析，在探討了這些數(shù)值積分算法的優(yōu)缺點(diǎn)的理論之外，我們還將這些數(shù)值積分算法在計算機(jī)上通過matlab軟件編程實(shí)現(xiàn)應(yīng)用，并分別用各自求積公式進(jìn)行運(yùn)算，以此來分析比較各種求積公式的代數(shù)精度和計算誤差。關(guān)鍵詞：數(shù)值積分；求積公式；代數(shù)精度貢州師范學(xué)院畢業(yè)論文（設(shè)計）Abstractfunctionisverycomplexthatitisdifficulttofindtheelementaryfunctions,whichmakesuWeinthefunctio

18、nfortheintegration,oftenbecausetheoriginalsverydifficulttocalculate,sowecanonlythinkofawaytofindtheapproximatevalue,thusdirectlywithNewton-Leibnizformulacalculatingdefiniteintegralsituationisveryrare.Whennumericalintegrationistosolvethisprobleminaveryeffectivemethod.Fromthenumericalintegrationproble

19、m,introducessomemethodsofnumericalintegration(Newton-Cotesquadratureformula,compositequadratureformulas,Longbeilatticequadratureformula,Gausstypequadratureformulas)andhascarriedonbriefanalysis,discussestheadvantagesanddisadvantagesofthesenumericalintegrationalgorithmtheory,wewillthesenumericalintegr

20、ationalgorithminthecomputerbyMATLABsoftwareprogrammingapplication,andseparatelywiththeirrespectivequadratureformulaforcomputing,inordertoanalyzethealgebraiccalculationprecisionanderrorcomparisonofvariousquadratureformulas.Keywords:Numericalintegration;Calculationmeth;numericalanalysis#員州加范學(xué)院畢業(yè)論文(設(shè)計)

21、、八1月I言微積分的發(fā)明是世界數(shù)學(xué)史上一項輝煌的成就。但在實(shí)際求積問題的時候，求解積分卻有著非常多局限性。比如對于定積分Ff(x)dx在求某函數(shù)的定積分時，在一定條件下，雖然有牛頓萊布里茨公式仁阡(x)dx二F(b)一F(a)可弋a(chǎn)以計算定積分的值，但在很多情況下f(x)的原函數(shù)不易求出或非常復(fù)雜。被積函數(shù)f(x)的原函數(shù)很難用初等函數(shù)表達(dá)出來，例如()二沁(，等；有的X函數(shù)f(x)的原函數(shù)F(x)存在，但其表達(dá)式太復(fù)雜，計算量太大，有的甚至無法有解析表達(dá)式。因此能夠借助牛頓萊布尼茲公式計算定積分的情形是不多的。另外，許多實(shí)際問題中的被積函數(shù)往往是列表函數(shù)或其他形式的非連續(xù)函數(shù)對這類函數(shù)的定積

22、分，也不能用不定積分方法求解，只能設(shè)法求其近似值。因此，探討近似計算的數(shù)值積分方法是有明顯的實(shí)際意義的即有必要研究定積分的數(shù)值計算方法，以解決定積分的近似計算。而數(shù)值積分就是解決此類問題的一種有效的方法，它的特點(diǎn)是利用被積函數(shù)f(x)在一些節(jié)點(diǎn)上的信息求出定積分的近似值。在很多實(shí)際應(yīng)用中，只能知道積分函數(shù)在某些特定點(diǎn)的取值比如天氣測量中的氣溫、濕度、氣壓等，醫(yī)學(xué)測量中的血壓、濃度等等。通過研究，我們將會更熟練掌握一些數(shù)值積分方法去計算一些特定條件的數(shù)值計算，以便我們得到自己想要的結(jié)果。2數(shù)值積分方法的基本思想在數(shù)學(xué)分析中，計算連續(xù)函數(shù)f(x)在區(qū)間a,b上的積分是通過f(x)的原函數(shù)F(x),

23、由下列定積分公式baf(x)dx二F(b)-F得到的。但由于大量被積函數(shù)的原函數(shù)不能用初等函數(shù)表示，因此，很難用求原函數(shù)的公式4(x)dx=F(b)-F(a)得到積分；有些被積函數(shù)f(x)不是明顯知道的，例如由數(shù)值表給出它的離散值，或者是它被定義為某個微分方程的解，而這個微分方程是不能顯示解出的。這說明按f(x)dx=F(b)一F(a)公式計算La定積分是有很大局限性的。因而常常采用在電子計算機(jī)上很有效的數(shù)值積分方法。我們從定積分的定義bnaf(x)dx:(b-a)C(n,k)f(xk)k=0bn出發(fā)。推導(dǎo)出兩個簡單的數(shù)值積分公式。Jaf(x)dx=n(kE)AXk式的幾何意義，-ka就是把整

24、塊曲線梯形的面積積分成若干個小曲邊梯形面積的和，當(dāng)無限細(xì)分時這個和取極限就是真正曲邊梯形面積。去掉取極限這一步，用有限個小曲邊梯形面積的和，代替整塊的曲邊梯形面積，從而求得一個近似值，這就是數(shù)值積分的基本思想。根據(jù)小區(qū)間的不同分割方法和各分點(diǎn)f(z)值的不同選擇，就得到不同的數(shù)值積分公式。數(shù)值求積公式是取lab1上若干個點(diǎn)Xk處的高度f(Xk)，通過加權(quán)人后，再n求和、Akf(Xk)k=0從而得到積分的近似值。數(shù)值求積公式寫成一般形式nf(x)dxAAJ(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的選取有關(guān)，而不依賴

25、被積函數(shù)f(x)的具體形式。記#貴州師范學(xué)院畢業(yè)的文（設(shè)計）Rff(x)dxA.f(xk)-k=0#員州師范學(xué)院畢業(yè)論文(設(shè)計)把RIf1稱為求積公式的截斷誤差或余項。數(shù)值求積方法的特點(diǎn)是直接利用積分區(qū)間l.a,b1上一些離散節(jié)點(diǎn)的函數(shù)值進(jìn)行線性組合來近似計算定積分的值，從而將定積分的計算歸結(jié)為函數(shù)值的計算，這就避開了牛頓萊布尼茲公式需要尋求原函數(shù)的困難，并為計算機(jī)求積分提供了可行性。3幾類常用數(shù)值積分方法的簡單分析3.1Newton-Cotes求積公式常用的梯形公式和Simpsor公式是低階的牛頓柯特斯公式，牛頓柯特斯公式是積分區(qū)間上等距節(jié)點(diǎn)的插值求積公式。插值求積公式在積分區(qū)間上，所取節(jié)點(diǎn)

26、是等距時稱為牛頓柯特斯公式，即nf(x)dx:(b-a)二C(n,k)f(xk)k=0其中c(n,k)為Cotes求積公式的系數(shù)，是n和k的函數(shù)。當(dāng)n=1時，為梯形公式：f(x)dx9ff(b)2梯形公式的代數(shù)精度為1，有兩個積分節(jié)點(diǎn)。當(dāng)n=2時，為Simpson公式：b(b-a)af(x)dx定1f(a)4f(；f(b)6Simpson公式的代數(shù)精度為3，有三個積分節(jié)點(diǎn)由于只增加一個節(jié)點(diǎn)，其代數(shù)精度增加2，由此可知，Simpson公式比梯形公式代數(shù)精度高。當(dāng)n=4時，NewtonCotes求積公式為Cotes公式：:f(x)dx烤7f(a)32f(專)12f(寧)32f(寧)7f(b)New

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

28、式有：復(fù)化梯形公式：b,a1nj1變步長梯形公式為：=+f(Xo)”)=豬12Mn二寧1f(X2jj)2nj#復(fù)化Simpson公式：b-a人盲“心)忖眸)變步長復(fù)化Simpson公式：nb_an心S2n八Ikf(X。)fgn)4、f(X2j2f(X2k)k：i6nkdkT3.3 Romberg求積公式Romberg積分方法也叫做逐次分半加速法，它是在復(fù)化梯形公式誤差估計的基礎(chǔ)上，應(yīng)用線性外推的方法構(gòu)造出的一種加速算法。將積分區(qū)間分成n等分和2n等分時，求得積分近似值和T.，并沒有誤差估計式1|.T2n3(T2n-Tn)3積分近似值T2n的誤差大致等于3仃2人)，當(dāng)用-)對T?.進(jìn)行修正時，3

29、3-(T2n-Tn)與T2n之和比T?n更接近于真值I,故陶-)是對T?n誤差的一種補(bǔ)償，33因此可以期望下式是一個更好的結(jié)果，即T二T2n.(J-)二今n2-；Tn333下面說明T即是分成n等分時Simpson公式的值5。將復(fù)化梯形公式f(a)2、f(xOf(b)梯形變步長求積公式T2ndin4f(Xf22心k2代人上式T表達(dá)式得-hn4n4T=;Jf(a)+4送f(x-)+2送f(xj+f(b)6-心遷心這就是說，用梯形法二分前后兩個梯形值和丁加作線性外推，結(jié)果得到Simpson法的積分值&。將誤差由。(忤)變?yōu)椤?幣)，從而提高了逼近精度。再考察Simpson法。其截斷誤差與h成正比，因

30、此，若將步長折半,則誤差減至，即有16I-En. -Sn16一一-7貴州師范學(xué)院畢業(yè)的文（設(shè)計）WWI：16S2-Sn1515不難驗(yàn)證，上式右端的值其實(shí)等于Cn，就是說，用Simpson法二分前后的兩個積分值&與s2n,按上式再作線性外推，結(jié)果得到柯特斯法的積分值G，即有161CnS2nSn1515這時將誤差由。（的變?yōu)?。（的，逼近精度又一次得以提高。同樣的方法，依?jù)柯特斯法的誤差公式，可進(jìn)一步導(dǎo)出下列龍貝格公式641“nC2nCn6363用逼近積分值的誤差為0（M，這樣Romberg公式將誤差由。（忖變?yōu)閛（h），逼近精度再次得以提高。Romberg公式有7次代數(shù)精度，這表明該公式不是牛頓柯

31、特斯公式。在步長二分的過程中運(yùn)用Sn、Cn、Rn表達(dá)式加工三次，就能將粗糙的積分值Tn逐步加工成精度較高的Romberg值R”，或者說，將收斂緩慢的梯形值序列入加工成收斂迅速的Romberg值序列R”，這種加速方法稱Romberg算法。3.4高斯型求積公式前面介紹的n-1個節(jié)點(diǎn)的Newton-Cotes求積公式，其特征是節(jié)點(diǎn)是等距的。這種特點(diǎn)使得求積公式便于構(gòu)造，復(fù)化求積公式易于形成。但同時也限制了公式的精度。n是偶數(shù)時，代數(shù)精度為n1，n是奇數(shù)時，代數(shù)精度為n；我們知道n1個節(jié)點(diǎn)的插值型求積公式的代數(shù)精確度不低于no能不能在區(qū)間a.b1上適當(dāng)選擇個節(jié)點(diǎn)Xo,X1,X2.x使插值求積公式的代數(shù)

32、精度高于n呢？答案是肯定的，適當(dāng)選擇節(jié)點(diǎn)，可使公式的精度最高達(dá)到2n1，這就是所學(xué)的高斯型求積公式。不失一般性，將求積公式/0f（x）dx八Axf（Xk）的求積區(qū)間La,b轉(zhuǎn)換成T,1ak=0的形式。對任意求積區(qū)間la,b作變換D-a,a夕ttb a _L a b *b a _L a b j 22 b 8 ?(t) dt1-1,1上，這時可以變換到區(qū)間af(X)dX=乜b-a2-baf出二2其中-仕）二f（batabt）。22高斯-勒讓德求積公式在這里簡稱高斯公式，它是在區(qū)間1-1,1上進(jìn)行討論的。4幾類數(shù)值積分方法的簡單比較評述由于我們在計算實(shí)際問題是往往要考慮到代數(shù)精度和計算量，所以不同類

33、型的求積公式有著不同的特點(diǎn)：Simpson積分方法和梯形積分方法雖然計算簡便，但是精度比較差，不理想。但對于光滑性較差的被積函數(shù)有時會比高精度的積分方法更為有效。特別是梯形積分方法對被積函數(shù)是周期函數(shù)的求積效果更為突出。n7時，Newton-Cotes公式是不穩(wěn)定的，然而復(fù)化梯形公式和復(fù)化Simpson公式不僅保留了低階公式的優(yōu)點(diǎn)還能夠獲得比較較高的精度，所以在實(shí)際計算中應(yīng)用得最為廣泛。Romberg積分方法的算法簡單，方便編程的實(shí)現(xiàn)。收斂速度快、計算精度較高，但是計算量較大。Gauss積分方法的精度較高，數(shù)值穩(wěn)定、收斂速度較快，但因?yàn)槠涔?jié)點(diǎn)不規(guī)則，計算比較麻煩5利用MATLAB編程應(yīng)用對幾類

34、求積算法的分析比較在簡單的認(rèn)識積分方法比且理論比較之后，則要進(jìn)行數(shù)學(xué)實(shí)驗(yàn)進(jìn)行驗(yàn)證，因此就要通9員州師范學(xué)院畢業(yè)論文(設(shè)計)過matlab軟件對各種積分方法進(jìn)行編程并運(yùn)算，然后對其各種方法的運(yùn)算結(jié)果進(jìn)行分析比較，掌握和理解各方法的優(yōu)缺點(diǎn)。規(guī)定各個程序都以I=0皿dx為例子進(jìn)行運(yùn)算。原積分的精確值為電XisinxoXdx”946083070367183例分別用不同的方法計算積分I二，沁dx，并作比較。x用以上介紹的幾類積分方法分別計算積分，得出誤差，并進(jìn)行比較：1、用Newton-Cotes公式當(dāng)n=1時，即用梯形公式，用程序一(程序見附錄)在MATLAB命令窗口中輸入NCotes(0,1,1,2

35、)得I.92703549240395Ir=0.01904757796323當(dāng)n=2時，即用Simpson公式，用程序一(程序見附錄)在MATLAB命令窗口中輸入NCotes(0,1,2,2)得I肝0.94614588227359Ir=0.000062811906407當(dāng)n=4時，即用科特斯公式，用程序一(程序見附錄)在MATLAB命令窗口中輸入NCotes(0,1,4,2)得I-0.94608300406367R=0.0000000663035132、用復(fù)化梯形公式令h=1/8=0.125,用程序二(程序見附錄)在MATLAB命令窗口中輸入trapri(0,1,8)，得J;叱dxf(0)+2f

36、(h)+.+f(7h)+f(l)=0.9456908635270x2R=0.0003922068401823、用復(fù)化Simpson公式令h=1/8=0.125,用程序三(程序見附錄)在MATLAB命令窗口中輸入simpr1(f,0,1,8),得1SnXdx：h(f(0)4f(h)f(7h)b：；,2f(2h)f(6h)IHf(10x3=0.94608308538495R=0.0000000150177674、用Romberg公式用程序四(程序見附錄)在MATLAB命令窗口中輸入romber(T,0,1,5,0.5*(10/(-8)得(沁dx4.94608307036718|R=0.000000

37、000000002X5、用高斯勒讓德求積公式令x=(t+1)/2，|=fR49Zdt=t+1(1)用2個節(jié)點(diǎn)的Gauss公式10.94604115827633(2)用3個節(jié)點(diǎn)的Gauss公式，用程序五(程序見附錄)在MATLAB命令窗口中輸入GuassLegendre(0,1,2,2)得1-0.946083134078473|R=0.000000063711290算法比較：1 .原積分的精確值為：1 sinx1 =dx=0.946083070367183x2 .由例題的各種求積算法可知：(1)對Newton-cotes公式，當(dāng)n=1時只有1位有效數(shù)字,當(dāng)n=2時有3位有效數(shù)字，當(dāng)n=4時有7位

38、有效數(shù)字。(4)用復(fù)化梯形公式有2位有效數(shù)字，對復(fù)化Simpson公式有7位有效數(shù)字。(5)用復(fù)化梯形公式，對積分區(qū)間0,1二分11次用了2049個函數(shù)值，才可以得7位有效數(shù)字。(6)用Romberg公式對區(qū)間0,1二分3次用了9個函數(shù)值，就可以得到7位有效數(shù)字；二分4次用了14個函數(shù)值，卻可以得到14位有效數(shù)字。(7)用高斯勒讓德求積公式僅僅用了3個函數(shù)值，就能得到比較精確的6位有效數(shù)字。11貢州師范學(xué)院畢業(yè)論文（設(shè)計）結(jié)束語本文主要研究了常用的幾類數(shù)值積分的求積算法并通過例題計算積分進(jìn)行分析比較。Newton-Cotes積分方法是一種非常普遍的積分方法，然而梯形積分方法的誤差最大，近似效果

39、最差，Simpson積分方法的精度比梯形積分方法高了一個數(shù)量級；Cotes積分方法精度比Simpson積分方法高兩個數(shù)量級。則Cotes代數(shù)精度比較高。由此可知一般情況下，積分公式代數(shù)精度越高，計算精度也越高。但是高階的Cotes積分方法收斂性沒有保證，因此實(shí)際應(yīng)用中很少用。復(fù)化梯形積分方法比梯形積分方法精度高，同樣的，復(fù)化Simpson積分方法比Simpson積分方法精度高，高了差不多7個數(shù)量級，所以復(fù)化積分方法比較優(yōu)越。Romberg積分方法收斂速度快、計算精度較高，但是計算量較大。Gauss積分方法精度高、數(shù)值穩(wěn)定、收斂速度較快，但是計算麻煩。經(jīng)研究可以知道Newton-Cotes方法的

40、代數(shù)精度越高，數(shù)值積分的效果越好、越精確。當(dāng)積分區(qū)間比較大的時候，積分?jǐn)?shù)值不穩(wěn)定，這個時候可以利用復(fù)化積分方法效果會更好；Romberg積分方法可以利用變步長復(fù)化積分公式得到更為精確的數(shù)值結(jié)果，是比較好的積分方法。高斯求積方法精確度高，收斂性快，比其他積分方法優(yōu)越。具有很廣泛15的運(yùn)用。參考文獻(xiàn)1關(guān)治.陸金甫.數(shù)學(xué)分析基礎(chǔ)（第二版）M.北京：等教育出版社.2010.72胡祖熾.林源渠.數(shù)值分析M北京:等教育出版社.1986.33薛毅.數(shù)學(xué)分析與實(shí)驗(yàn)M北京:業(yè)大學(xué)出版社2005.34徐士良.數(shù)值分析與算法M.北京版工業(yè)出版社2007.15王開榮.楊大地.應(yīng)用數(shù)值分析M北京:等教育出版社2010.

41、76楊一都.數(shù)值計算方法M.北京:等教育出版社.2008.47韓明.王家寶.李林.數(shù)學(xué)實(shí)驗(yàn)（MATLAB）版M.上海:濟(jì)大學(xué)出版社2012.18圣寶建.關(guān)于數(shù)值積分若干問題的研究J.南京信息工程大學(xué).2009.05.01.:429劉緒軍.幾種求積公式計算精確度的比較J.南京職業(yè)技術(shù)學(xué)院.2009.10史萬明.吳裕樹.孫新.數(shù)值分析M.北京理工大學(xué)出版社.2010.4.行文至此，我的這篇論文已接近尾聲；歲月如梭，我四年的大學(xué)時光也即將敲響結(jié)束的鐘聲。離別在即，站在人生的又一個轉(zhuǎn)折點(diǎn)上，心中難免思緒萬千，一種感恩之情油然而生。首先感謝貴州師范學(xué)院四年來對我的培養(yǎng)，是博學(xué)的老師們教會了我學(xué)習(xí)的方法、

42、鍛煉了我思考的能力、指明了我未來奮斗的方向，從而使我進(jìn)一步明確了人生的目標(biāo)。其次，我要感謝我的指導(dǎo)老師一雍進(jìn)軍老師，他的嚴(yán)謹(jǐn)細(xì)致、一絲不茍的作風(fēng)一直是我工作、學(xué)習(xí)中的榜樣；他的循循善誘的教導(dǎo)和不拘一格的思路給予我無盡的啟迪。在撰寫整個畢業(yè)論文的過程當(dāng)中，他為我們考慮到了每一個細(xì)節(jié)，從開題報告到畢業(yè)論文的擬定修改上，雍老師更是不厭其煩的為我們做好每一步的細(xì)心指導(dǎo)。對此，我表示衷心地感謝。沒有雍老師，我的論文也不可能這么順利的完成。同時，我也要感謝每一位給過我?guī)椭睦蠋熀屯瑢W(xué)，在我撰寫論文的過程當(dāng)中同樣給了我大量有益的建議，在此一并向他們表示真誠的感謝，感謝他們對我的支持和幫助。最后感謝這篇論文所

43、涉及到的各位學(xué)者，本文引用了數(shù)位學(xué)者的研究文獻(xiàn)，如果沒有各位學(xué)者的研究成果帶給我的的幫助和啟發(fā)，我將很難完成本篇論文的寫作。由于我的學(xué)術(shù)水平有限，所寫論文難免有不足之處，懇請各位老師和學(xué)友批評指正。最后，衷心感謝評閱論文及參加答辯的各位老師！員州師范學(xué)院畢業(yè)論文(設(shè)計)1 NewtonCotes求積公式的MATLAB實(shí)現(xiàn)先用M文件定義一個名為fl.m的函數(shù)：%i是要調(diào)用第幾個被積函數(shù)g(i)，x是自變量functionf=f1(i,x)g(1)=sqrt(x);ifx=0g(2)=1；elseg(2)=sin(x)/x;endg(3)=4/(1+xA2);仁g(i);程序一：functionC

44、,g=NCotes(a,b,n,m)%a，b分別為積分的上下限；%n是子區(qū)間的個數(shù)；%m是調(diào)用上面第幾個被積函數(shù)；%當(dāng)n=1時計算梯形公式；當(dāng)n=2時計算辛浦生公式，以此類推;i=n；h=(b-a)/i;z=0;forj=0:ix(j+1)=a+j*h;s=1;ifj=Os=s;elsefork=1:js=s*k;endendr=1;ifi-j=0r=r;elsefork=1:(i-j)r=r*k;endendifmod(i-j),2)=1q=-(i*s*r);elseq=i*s*r;endy=i;fork=0:iifk=jy=y*(sym(t)-k);endendl=int(y,0,i);C

45、(j+1尸l/q；z=z+C(j+1)*f1(m,x(j+1)；endg=(b-a)*z1)當(dāng)輸入a=0,b=1,n=1A=2時，即在MATLAB命令窗口輸入NCotes(0,1,1,2)即可得用梯形公式的積分值和相應(yīng)科特斯系數(shù)如圖1.12)當(dāng)輸入a=0,b=1，n=2,m=2時，即在MATLAB命令窗口輸入NCotes(0,1,2,2)即可得用辛浦生公式的積分值和相應(yīng)科特斯系數(shù)如圖1.23)當(dāng)輸入a=0,b=1，n=4,m=2時，即在MATLAB命令窗口輸入NCotes(0,1,4,2)即可得用科特斯公式的積分值和相應(yīng)科特斯系數(shù)如圖1.317貴州師范學(xué)院畢業(yè)論文（設(shè)計）CommandWind

46、owHCotesLL2)8293248040994423/9007199254740992ans=(1/2,1/2)8293248040994423/90071992547409920.92073549240395CommandWindowNCotes(01.2.2)8522124485690909/9007199254740992ans=(1/6,2/3,1/68522124485690909/9007199254740992ans=0.94614588227359圖1.219員州師范學(xué)院畢業(yè)論文(設(shè)計)CommandWindowNCotes(Qb42)g=38347011581064500

47、9/405323966463344640axis=7/90,16/45,2/15,16/45,7/90383470115810645009/405323966463344640an，=0.94608300406367圖1.32復(fù)化梯形求積公式的MATLAB實(shí)現(xiàn)通過f(x)的n1個等步長節(jié)點(diǎn)逼近積分:/(幼辦怎(/()+/(*)+方/Im)其中,xtkh,Xo=a.程序二:functions=trapr1(f,a,b,n)%f是被積函數(shù)；%a,b分別為積分的上下限；%n是子區(qū)間的個數(shù)；%S是梯形總面積；h=(b-a)/n;s=0;fork=1:(n-1)x=a+h*k;s=s+feval(f,x);endformatlongs=h*(feval(f,a)+feval(,f,b)/2+h*s;先用M文件定義一個名為f.m的函數(shù)：functiony=f(x)ifx=0y=i;elsey=sin(x)/x;end若取子區(qū)間的個數(shù)在MATLAB命令窗口中輸入trapri(