數(shù)值積分的matlab實現(xiàn)-文庫吧資料
2025-01-12 06:51本頁面
【正文】 _trap(f4,a,b)。 a=0。 首先編寫被積函數(shù)的 M 文件 function y=f4(x) y=4\(1+x^2)。與精確值比較的誤差為 :39。積分近似值分別為 :39。積分準(zhǔn)確值為 :39。 format long, Sj=double(S1)。 syms x f= x*sqrt(1+x^2)。 FS4=simpson(f3,a,b,8)。 S=chen_simpson(f3,a,b)。b=3。 然后編寫如下程序 ,并保存為 % clear。與精確值比較的誤差為 :39。積分近似值分別為 :39。積分準(zhǔn)確值為 :39。 format long, Sj=double(S1)。 syms x f=x*sin(x)。 FS4=simpson(f2,a,b,8)。 S=chen_simpson(f2,a,b)。b=2*pi。 然后編寫如下程序 ,并保存為 % clear。), wuchi=[double(abs(SjT)),double(abs(SjS)),double(abs(SjFT8)),double(abs(SjFS4)),double(abs(Sjquad))] 202 ??( ) sinx xdx。),[T,S,FT8,FS4,quad] disp(39。),Sj disp(39。 disp(39。 S1=int(f, x, 0,1)。 quad=Romberg(f1,a,b,10e6)。 FT8=trap(f1,a,b,8)。 T=chen_trap(f1,a,b)。 a=0。 end quad=R(J+1,J+1)。 for k=1:J R(J+1,k+1)=R(J+1,k)+(R(J+1,k) R(J,k))/(4^k1)。 end R(J+1,1)=R(J,1)/2+h*S。 for p=1:M x=a+h*(2*p1)。 h=h/2。 err=1。 M=1。 %復(fù)化 simpson 公式主程序 %fun 表示被積函數(shù)句柄 %a,b 分別為積分區(qū)間的左右端點 %n=2m,m 為等分積分區(qū)間的個數(shù) %y 為積分近似值 Sm Romberg求積公式程序代碼 程序代碼說明 function [quad,R]=Romberg(f,a,b,eps)； h=ba。 z3= z3+4*feval(fun,x3)。 z2=z2+2*feval(fun,x2)。 x3=0。 z3=0。 x=a。m=n/2。 end y=(z1+z2)*h/2。 for k=1:n1 x1=x+k*h。 z2=0。 %simpson 公式主程序 %fun 表示被積函數(shù)句柄 %a,b 分別為積分區(qū)間的左右端點 復(fù)合梯形公式程序代碼 程序代碼 說明
公司管理相關(guān)推薦
鄂ICP備17016276號-1