Matlab數值積分程序集合[圖書館+網絡收集]

近來學習數值積分,手頭積累了不少程序,也拿來和各位朋友分享一下。。。主要是來自數值積分教材和網絡，基本的原理也就不打算多說了，隨便搜索一下就可以得到，那就開始上代碼了，呵呵，非原創，但是全部驗證過，有疑問可以給我e-mail：

1 梯形數值積分的MATLAB主程序

function T=rctrap(fun,a,b,m)

%fun 函數，a 積分上限 b積分下限 m 遞歸次數

n=1;h=b-a; T=zeros(1,m+1); x=a;

T(1)=h*(feval(fun,a)+feval(fun,b))/2;

for i=1:m

h=h/2; n=2*n; s=0;

for k=1:n/2

x=a+h*(2*k-1); s=s+feval(fun,x);

end

T(i+1)=T(i)/2+h*s;

end

T=T(1:m);

e.g

運行后屏幕顯示精確值F s，用rctrap計算的遞歸值T和T與精確值F s的絕對誤差w T

>> ) exp((-x^.2./2)./(sqrt(2*pi)))

T=rctrap(fun,0,pi/2,14), syms t

fi=int(exp((-t^2)/2)/(sqrt(2*pi)),t,0, pi/2);

Fs= double(fi), wT= double(abs(fi-T))

fun =

@(x)exp((-x^.2./2)./(sqrt(2*pi)))

T =

總結

以上是生活随笔為你收集整理的matlab如何对wt积分,Matlab数值积分程序集合的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。