復(fù)合梯形公式_第1頁
復(fù)合梯形公式_第2頁
復(fù)合梯形公式_第3頁
復(fù)合梯形公式_第4頁
復(fù)合梯形公式_第5頁
已閱讀5頁,還剩12頁未讀, 繼續(xù)免費(fèi)閱讀

下載本文檔

版權(quán)說明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請進(jìn)行舉報(bào)或認(rèn)領(lǐng)

文檔簡介

復(fù)合梯形公式0《復(fù)合梯形公式》實(shí)驗(yàn)報(bào)告_________?一、實(shí)驗(yàn)?zāi)康?1)掌握數(shù)值積分函數(shù)的調(diào)用格式(2)掌握復(fù)合梯形公式的思想和構(gòu)造過程(3)能夠應(yīng)用matlab軟件編寫復(fù)合梯形公式的程序并能熟悉應(yīng)用,以此來解決相關(guān)例題(4)利用復(fù)合梯形公式求數(shù)值積分的近似值,以解決其它科學(xué)實(shí)驗(yàn)的計(jì)算問題二、實(shí)驗(yàn)內(nèi)容驗(yàn)證積分函數(shù)的調(diào)用格式并求函數(shù)積分值編寫復(fù)合梯形公式程序并驗(yàn)證0能得到什么結(jié)論三、實(shí)驗(yàn)環(huán)境該實(shí)驗(yàn)應(yīng)用matlab2014來進(jìn)行實(shí)驗(yàn)的驗(yàn)證和設(shè)計(jì).四、實(shí)驗(yàn)步驟和實(shí)驗(yàn)結(jié)果驗(yàn)證積分函數(shù)的調(diào)用格式并求函數(shù)積分值例3:(數(shù)值積分)用兩種不同的方法求:I=j1ex2dx.0法一:建立被積函數(shù)文件法:functionex=ex(x)ex=exp(-x.^2);formatlongI=quad('ex',0,1)I=quadl('ex',0,1)法二:不建立關(guān)于被積函數(shù)的函數(shù)文件時(shí):g=inline('exp(-x.^2)');I=quadl(g,0,1)=例4:(數(shù)值積分)計(jì)算二重積分I=j1j2ex22sin(x2+y)dxdy(1)建立一個(gè)m文件functionf=fxy(x,y)globalki;ki=ki+1;f=exp(-x.^2/2).*sin(x.^2+y);globalki;ki=0;I=dblquad('fxy',-2,2,-1,1)iki=f=inline('exp(-x.^2/2).*sin(x.^2+y)','x','y');I=dblquad(f,-2,2,-1,1)=編寫復(fù)合梯形公式程序并驗(yàn)證0能得到什么結(jié)論建立復(fù)合梯形公式的文件:function[H]=fhTX(a,b,n)%ab分別是積分函數(shù)的上下限%%n為區(qū)間等分份數(shù)%%I返回的是得到的積分近似值%h=(b-a)/n;%h為步長%x=a:h:b;%等分點(diǎn)出的節(jié)點(diǎn)Xi%f=exp(x.^2).*sin(x);%被積函數(shù)表達(dá)式%fork=1:n;%節(jié)點(diǎn)xi開始循環(huán)%J=f(k)+f(k+1);%復(fù)合梯形公式的遞推公式%I=I+J;%遞推公式求和%H=I*h/2;EndH=fhTX(0,1,3)%函數(shù)調(diào)用%H=H%得到近似解%I1=quad('exp(x.^2).*sin(x)',0,1)%quad計(jì)算真值%R=abs(I1-H)%近似解和真值之間的誤差%ans=所以,將區(qū)間3等分后得到的結(jié)果為,和真值之間存在的誤差為.(2)將區(qū)間4等分得到的結(jié)果H=fhTX(0,1,4)%函數(shù)調(diào)用%H=H%得到近似解%abs(I1-H)ans=%近似解和真值之間的誤差%通過結(jié)果可以看出,將區(qū)間4等分后的結(jié)果和真值之間的誤差縮小為,較三等分得到的結(jié)果略微良好.(3)將區(qū)間6等分得到的結(jié)果H=fhTX(0,1,6)%函數(shù)調(diào)用%H=H%得到近似解%abs(I1-H)ans=%近似解和真值之間的誤差%精確(4)將區(qū)間9等分得到的結(jié)果五、實(shí)驗(yàn)討論、結(jié)論H=fhTX(0,1,9)%函數(shù)調(diào)用%H=H%得到近似解%abs(I1-H)ans=%近似解和真值之間的誤差%這時(shí)的誤差縮小為,更加靠近真值了.通過四種不同的等分區(qū)間來對積分值進(jìn)行近似估計(jì),我們可以發(fā)現(xiàn),區(qū)間等分份數(shù)越多,得到的結(jié)果越靠近真值,得到的結(jié)果越精確我們用相同的方法得到更多數(shù)據(jù)的結(jié)果入下:0nn從表格中我們很清楚的可以看到在某個(gè)較大的范圍內(nèi),當(dāng)n越大時(shí),得到的結(jié)果和真值越接近,當(dāng)達(dá)到一定程度時(shí),近似值將不會再變化.化通過四種不同的等分區(qū)間來對積分值進(jìn)行近似估計(jì),我們可以發(fā)現(xiàn),區(qū)間等分份數(shù)越多,得到的結(jié)果越靠近真值,得到的結(jié)果越精確,但符合梯形公式在對積分值進(jìn)行估計(jì)的過程中存在的問題是:在某個(gè)較大的范圍內(nèi),無論區(qū)間等分多少份,或者說無論步長取得多么小,總會存在一個(gè)更小的步長h,使得結(jié)果過更加精確,但當(dāng)超過這個(gè)規(guī)定的范圍,當(dāng)區(qū)間等分份數(shù)過大時(shí),得到的近似值將不會再變價(jià).所以說,采用梯形公式得到的結(jié)果,隨著區(qū)間等分份數(shù)的增多,會越來越接近真值,但永遠(yuǎn)不能得到最精確的結(jié)果.所以,在對某些函數(shù)進(jìn)行積分計(jì)算的時(shí)候,一些簡單的函數(shù)我們可以直接進(jìn)行積分到精確的結(jié)果,但是對于一些

溫馨提示

  • 1. 本站所有資源如無特殊說明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
  • 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
  • 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁內(nèi)容里面會有圖紙預(yù)覽,若沒有圖紙預(yù)覽就沒有圖紙。
  • 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
  • 5. 人人文庫網(wǎng)僅提供信息存儲空間,僅對用戶上傳內(nèi)容的表現(xiàn)方式做保護(hù)處理,對用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對任何下載內(nèi)容負(fù)責(zé)。
  • 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容,請與我們聯(lián)系,我們立即糾正。
  • 7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時(shí)也不承擔(dān)用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。

最新文檔

評論

0/150

提交評論