matlab如何对wt积分,Matlab数值积分程序集合

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数值积分程序集合的全部内容，希望文章能够帮你解决所遇到的问题。

如果觉得生活随笔网站内容还不错，欢迎将生活随笔推荐给好友。