声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 4082|回复: 6

[综合讨论] 为什么数值积分的结果为NaN

[复制链接]
发表于 2007-10-9 12:17 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
有这样一个程序:
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
为什么用数值积分算出来的结果为非数?是什么地方出错了?请高人指点
回复
分享到:

使用道具 举报

发表于 2007-10-9 12:59 | 显示全部楼层
quad(gg,0,10)
ans =
   0.499997178535039
数值积分上限取到inf当然不行了。事实上取到10被积函数就足够小了。

评分

1

查看全部评分

 楼主| 发表于 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进行积分。半周期最大,整个周期就是零了,我想道理应该是一样的
发表于 2007-10-10 12:01 | 显示全部楼层

回复 #4 花如月 的帖子

呵呵,积分项是个衰减的指数,不是周期函数,呵呵。
 楼主| 发表于 2007-10-10 12:19 | 显示全部楼层
被积函数非负,随着积分上限的增加,积分结果不可能减小呀。不知道为什么我那个积分会在20以后迅速减小
发表于 2007-10-10 12:28 | 显示全部楼层

回复 #5 w89986581 的帖子

谢谢你的意见,没细看积分公式:@L
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-24 19:31 , Processed in 0.100965 second(s), 27 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表