matlab求解erf误差函数方法(有详细注释 ).doc

《matlab求解erf误差函数方法(有详细注释 ).doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、%计算f(t)=exp(-t^2)在[0,1]上的定积分>>symsx%引入一个符号变量int(exp(-x^2),0,1)%int函数的主要作用是实现exp(-x^2)在（0,1）上的积分；ans=1/2*erf(1)*pi^(1/2)　%erf为误差函数，得到求解的表达式，是一个存在误差函数的表达式；>>vpa(int(exp(-x^2),0,1))%vpa函数实现的是运用数值近似求解来得到表达式int(exp(-x^2),0,1)的数值解；ans=0.24270>>d=0.001;x=0:d:1;d*tr

2、apz(exp(-x.^2))%trapz其实就是Trapezoidal（梯形的简写），梯形面积法ans=　　0.7468>>quad('exp(-x.^2)',0,1,1e-8)ans=0.7468%%求函数'exp(-x*x)的定积分，积分下限为0，积分上限为1。>>fun=inline('exp(-x.*x)','x');%用内联函数定义被积函数fname>>Isim=quad(fun,0,1)%辛普森法Isim=0.6425IL=quadl(fun,0,1)%牛顿－柯特斯法IL=0.8447%%二重积分

3、f=@(x,y)exp(sin(x))*ln(y)，y从5*x积分到x^2，x从10积分到20y=quad2d(@(x,y)exp(sin(x)).*log(y),10,20,@(x)5*x,@(x)x.^2)