復化梯形公式求積分.doc

《復化梯形公式求積分.doc》由會員上傳分享，免費在線閱讀，更多相關(guān)內(nèi)容在教育資源-天天文庫。

1、武漢工程大學計算機科學與工程學院《計算方法》實驗報告實驗項目復化梯形公式求積分實驗類別綜合實驗實驗目的及要求實驗目的：學會用復化梯形公式求函數(shù)的積分，并應用該算法于實際問題。實驗要求：要求能隨意輸入被積函數(shù)，進行算法設計，打印出誤差限例題：求被積函數(shù)在0

2、月日實驗內(nèi)容設計分析復化數(shù)值積分:將區(qū)間[a,b]n等分，取等距節(jié)點由定積分的區(qū)間可加性，有在每一個小區(qū)間上利用梯形積分公式有一般記稱做n+1點復化梯形積分公式。數(shù)學公式：算法描述：Step1:輸入a,b和正整數(shù)n；Step2:置h=(b-a)/n;Step3:F=f(a)+f(b);l=0;Step4:對j=1,2,…,n循環(huán)執(zhí)行5；Step5:置x=a+jh;l+=f(x);Step6:置T=h(F+2l)/2Step7:輸出T;程序源代碼:#include#includeusingna

3、mespacestd;doublef(doublex)//求函數(shù)的值；{return1/(1+pow(sin(x),2.0));}voidEchelonIntegral(intn)//梯形積分{doubley=0;doubleh=(1.0-0.0)/n;for(inti=0;i<=n;i++){doublea=0.0+i*h,b=0.0+(i+1)*h;y=y+h*(f(a)+f(b))/2.0;}cout<<"對應所求的梯形求積分為"<

4、;doublex=mid,p=2*sin(2*x)*sin(2*x)*(1+sin(x)*sin(x))-2*cos(2*x)*pow(1+sin(x)*sin(x),2);doublef2d=p/pow(1+sin(x)*sin(x),2);for(i=0;i

5、>m1;EchelonIntegral(m1);

6、cout<>answer1;while(answer1=='y'){cout<<"請輸入把0到1的范圍幾等分？"<<"t";cin>>m1;EchelonIntegral(m1);//3.直線求積分；cout<>answer1;}cout<

7、小區(qū)間上采取低階數(shù)值積分公式，然后把這些小區(qū)間上的數(shù)值積分結(jié)果加起來作為函數(shù)在整個區(qū)間上的近似，類似于分段差值。在此次實驗中，我學會了怎樣以專業(yè)的知識，從專業(yè)的角度來考慮所遇見的問題或者與專業(yè)相關(guān)的問題，并以自己的看法，給予實現(xiàn)。同時，在此次實驗中，也有一些不近人意的地方，也需要我不斷改進和練習。比如誤差怎樣取值，保證其最小。