skysurfer 发表于 2007-10-9 12:17

为什么数值积分的结果为NaN

有这样一个程序：
syms x;
g=.39894*exp(-.50000*x^2); %标准正态分布的概率密度
int(g,0,inf)
>>
ans =
.49999714194064567021687270136125

如果用数值积分计算：
gg＝inline(g);
quad(gg,0,inf)
>>
Warning: Infinite or Not-a-Number function value encountered.
(Type "warning off MATLAB:quad:ImproperFcnValue" to suppress this warning.)
> In D:\MATLAB6p5\toolbox\matlab\funfun\quad.m at line 91
ans =
   NaN
为什么用数值积分算出来的结果为非数？是什么地方出错了？请高人指点

rocwoods 发表于 2007-10-9 12:59

quad(gg,0,10)
ans =
   0.499997178535039
数值积分上限取到inf当然不行了。事实上取到10被积函数就足够小了。

skysurfer 发表于 2007-10-10 10:49

多谢楼上的回答。
还有一个问题：我的被积函数为
f＝(49551./50000.*erf(22097./2500.*x-7071./1000)+49551./50000).*exp(-25./2.*(x-1).^2)
积分限为0～inf
用数值积分
quad(f,0,20)
>>0.8182
quad(f,0,30)
>>8.3732e-006
上限取得越大，结果反而越小。在20左右比较稳定。这是怎么回事呀？

[ 本帖最后由 eight 于 2007-10-10 11:40 编辑 ]

花如月 发表于 2007-10-10 10:55

回复 #3 skysurfer 的帖子

和被积分函数的曲线有关了，对sin进行积分。半周期最大，整个周期就是零了，我想道理应该是一样的

w89986581 发表于 2007-10-10 12:01

回复 #4 花如月 的帖子

呵呵，积分项是个衰减的指数，不是周期函数，呵呵。

skysurfer 发表于 2007-10-10 12:19

被积函数非负，随着积分上限的增加，积分结果不可能减小呀。不知道为什么我那个积分会在20以后迅速减小

花如月 发表于 2007-10-10 12:28

回复 #5 w89986581 的帖子

谢谢你的意见，没细看积分公式:@L
页: [1]
查看完整版本: 为什么数值积分的结果为NaN