數(shù)值積分的matlab實現(xiàn)-在線瀏覽
2024-10-03 01:55本頁面
【正文】 調(diào)試等）出現(xiàn)的問題及其處理等廣泛的問題）. 實驗四 數(shù)值積分的Matlab實現(xiàn)一、實驗?zāi)康呐c要求1．熟練梯形公式、Simpson公式、復(fù)化梯形公式、復(fù)化Simpson公式、Romberg求積公式；2．熟悉符號積分，能數(shù)值積分的值與精確值的比較。3．培養(yǎng)編程與上機(jī)調(diào)試能力。y=(ba)/6*(feval(fun,a)+4*feval(fun,c)+feval(fun,b))。h=(ba)/n。 x=a。 z2= z2+2*feval(fun,x1)。%復(fù)化simpson公式主程序%fun表示被積函數(shù)句柄%a,b分別為積分區(qū)間的左右端點%n為等分積分區(qū)間的個數(shù)%y為積分近似值Sn 復(fù)合Simpson公式程序代碼 程序代碼說明 function y=simpson(fun,a,b,n)； z1=feval(fun,a)+feval(fun,b)。h=(ba)/(2*m)。z2=0。 x2=0。for k=1:m1 x2=x+2*k*h。 endfor k=1:m x3=x+(2*k1)*h。 endy=(z1+z2+z3)*h/3。R(1,1)=h*(feval(f,a)+feval(f,b))/2。 J=0。while J=3 | erreps J=J+1。 S=0。 S=S+feval(f,x)。 M=2*M。 end err=abs(R(J+1,J)R(J+1,J+1))。%f表示被積函數(shù)句柄%a,b表示被積區(qū)間[a,b]的端點%eps表示精度%quad是用Romberg加速算法求得的積分值%R為Romberg表%err表示誤差的估計 具體函數(shù)編程如下：（1）； function y=f1(x) 然后編寫如下程序，%clear。b=1。S=chen_simpson(f1,a,b)。FS4=simpson(f1,a,b,8)。syms x f=x^2*sqrt(2*x^2+3)。 format long, Sj=double(S1)。積分準(zhǔn)確值為:39。積分近似值分別為:39。與精確值比較的誤差為:39。function y=f2(x)y=x*sin(x)。a=0。T=chen_trap(f2,a,b)。F
環(huán)評公示相關(guān)推薦
鄂ICP備17016276號-1