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