MATLAB復(fù)化梯形法與龍貝格法計算定積分

1、精選優(yōu)質(zhì)文檔-傾情為你奉上姓名：樊元君 學(xué)號： 日期：2012.11.06 1.實驗?zāi)康模赫莆諒?fù)化梯形法與龍貝格法計算定積分。2.實驗內(nèi)容：分別寫出變步長梯形法與龍貝格法計算定積分的算法，編寫程序上機調(diào)試出結(jié)果，要求所編程序適用于任何類型的定積分，即能解決這一類問題，而不是某一個問題。實驗中以下列數(shù)據(jù)驗證程序的正確性。 求。3.程序流程圖： 變步長梯形法流程圖：龍貝格法流程圖：4.源程序： 變步長梯形法：function y =BTX(a,b,e)a=input(a=);b=input(b=);e=input(e=);h=(b-a);t1=(f(a)+f(b)*h)/2;t2=0; %準(zhǔn)備初值

2、% while abs(t2-t1)e %while語句控制計算精度以及控制計算時長 if t2=0 %按照變步長梯形法 h=h/2; %遞推公式求二分后積分值 t1=t2; end s=0; x=a+h/2; while x=4 s=0;x=a+h/2; while xb s=s+f(x);x=x+h; end t2=(t1/2)+(h/2)*s; s2=t2+(t2-t1)/3;if k=1 k=k+1;h=h/2;t1=t2;s1=s2;continueend c2=s2+(s2-s1)/15;if k=2c1=c2; k=k+1;h=h/2;t1=t2;s1=s2; continuee

3、ndr2=c2+(c2-c1)/63;if k=3 r1=r2;c1=c2; k=k+1;h=h/2;t1=t2;s1=s2;continueend if abs(r2-r1)e break else r1=r2;c1=c2; k=k+1;h=h/2;t1=t2;s1=s2; continue end end disp(插值結(jié)果=);disp(r2); endfunction y =f(x) %被積分函數(shù)% if x=0 y=sin(x)/x; else y=1; end end5.運行結(jié)果：6.實驗小結(jié)：（1）剛開始沒注意到積分中，被積函數(shù)在x=0點函數(shù)值為1，而編程中函數(shù)定義f(x)=sin(x)/x，當(dāng)x=0時，程序只能得到結(jié)果0/0=NaN（無窮大），導(dǎo)致后面積分計算算法完全失去意義，后來采取if和else以特殊情況對待x=0，才解決問題；（2）龍貝格法是對變步長梯形法的升華，在對復(fù)合梯形法理解透徹，并且編程合理，再去解決龍貝格法編程，容易很多；（3）充分理解兩個算法的內(nèi)涵和推理過程，對于編程正確性有很大幫助，剛開始由于對變步長梯形法理解不到位，一處循環(huán)采用if語句（本應(yīng)采用while）導(dǎo)致計算結(jié)果不受約束條件約束的錯誤；（4）由于龍貝格法至少要算4個T值后才有R值，但由于龍貝格法收斂