找回密码
 注册
Simdroid-非首页
查看: 76|回复: 2

[符号计算] 请教有关Bessel积分

[复制链接]
发表于 2012-8-17 10:08:45 | 显示全部楼层 |阅读模式 来自 湖北武汉
在论坛上看到有一个求Bessel积分的例子:exp(-x^2)*sin(x)*x*besselj(1,x),当时版主说可以用Mathematica计算比较合适,并给出了若x积分区间为 (0,inf),则结果为0.144742。但是由于本人目前对该软件不熟,试了下出现NIntegrate::nlim: x = inf is not a valid limit of integration.请教各位牛人这个积分应该怎么做?
ps:我用matlab的fplot看了下该函数的图像,应该是可以积出来的

本帖子中包含更多资源

您需要 登录 才可以下载或查看,没有账号?注册

×
发表于 2012-8-17 11:16:14 | 显示全部楼层 来自 陕西西安
Simdroid开发平台
Infinity,不是inf。
回复 不支持

使用道具 举报

 楼主| 发表于 2012-8-17 16:49:30 | 显示全部楼层 来自 湖北武汉
谢谢版主,没想到和matlab还是有很大区别。最后在mathematica里面得到了结果0.144742。
f[x_] := Exp[-x^2]*Sin[x]*x*Besselj[1, x];
NIntegrate[f[x], {x, 0, Infinity}];
1/4 HypergeometricPFQ[{5/4, 7/4}, {3/2, 5/2}, -1]
N[1/4 HypergeometricPFQ[{5/4, 7/4}, {3/2, 5/2}, -1]]
0.144742
回复 不支持

使用道具 举报

您需要登录后才可以回帖 登录 | 注册

本版积分规则

Archiver|小黑屋|联系我们|仿真互动网 ( 京ICP备15048925号-7 )

GMT+8, 2024-4-20 00:03 , Processed in 0.033945 second(s), 11 queries , Gzip On, MemCache On.

Powered by Discuz! X3.5 Licensed

© 2001-2024 Discuz! Team.

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