【答案】应助回帖

» 收录本帖的淘帖专辑推荐

数学问题集锦

» 猜你喜欢

投稿返修后收到这样的回复，还有希望吗已经有7人回复
压汞仪和BET测气凝胶孔隙率已经有4人回复
博士申请都是内定的吗？已经有14人回复
博士申请已经有3人回复
谈谈两天一夜的“延安行” 已经有13人回复
氨基封端PDMS和HDI反应快速固化已经有11人回复
之前让一硕士生水了7个发明专利，现在这7个获批发明专利的维护费可从哪儿支出哈？已经有11人回复
论文投稿求助已经有4人回复
Applied Surface Science 这个期刊。有哪位虫友投过的能把word模板发给我参考一下嘛已经有3人回复
投稿精细化工已经有6人回复

» 本主题相关价值贴推荐，对您同样有帮助:

1楼2015-01-11 16:56:26

feixiaolin

荣誉版主 (文坛精英)

专家经验: +518
应助: 942 (博后)
贵宾: 1.275
金币: 2930
散金: 58785
红花: 532
沙发: 11
帖子: 24215
在线: 2601.8小时
虫号: 2139575
注册: 2012-11-21
专业: 光学信息获取与处理
管辖: 数学

【0， inf】积分区间映射到【0, 1】，然后再数值求解。

2楼2015-01-11 17:20:32

duyang2012

铜虫 (小有名气)

应助: 1 (幼儿园)
金币: 3.3
散金: 141
红花: 2
帖子: 131
在线: 119.5小时
虫号: 2257395
注册: 2013-01-23
专业: 通信理论与系统

引用回帖:

2楼: Originally posted by feixiaolin at 2015-01-11 17:20:32
【0， inf】积分区间映射到【0, 1】，然后再数值求解。

不懂，可以说详细一点吗，谢谢！！！

3楼2015-01-11 17:34:54

feixiaolin

荣誉版主 (文坛精英)

专家经验: +518
应助: 942 (博后)
贵宾: 1.275
金币: 2930
散金: 58785
红花: 532
沙发: 11
帖子: 24215
在线: 2601.8小时
虫号: 2139575
注册: 2012-11-21
专业: 光学信息获取与处理
管辖: 数学

引用回帖:

3楼: Originally posted by duyang2012 at 2015-01-11 17:34:54
不懂，可以说详细一点吗，谢谢！！！...

取 u=1/x, v=1/a

[ 发自小木虫客户端 ]

4楼2015-01-11 18:27:55

pippi6

铁杆木虫 (著名写手)

工程和科学数值计算咨询

数学EPI: 6
应助: 413 (硕士)
贵宾: 0.002
金币: 7116.5
散金: 15
红花: 63
帖子: 1639
在线: 798.9小时
虫号: 2469437
注册: 2013-05-14
专业: 计算数学与科学工程计算

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
duyang2012: 金币+20, ★★★★★最佳答案, 完美解决问题 2015-01-12 11:10:45

我用maple 7
> f:=BesselI(0, alpha*x)*alpha*exp(-alpha^2)*x*exp(-(x^2+alpha^2)/2)
> *(1-exp(-x^2/2))^1022;

> s:=int(int(f,alpha=0..infinity),x=0..infinity);
s := 3299122687584745137204367643230771106304997323142764671919892526213274482351735424595090152580805831487201081992151782359365957048266356808329891738576418898384550273099110184494910603827424479869838674181059409833013059735407792907696592569535276701736861158593056561375652527873922258647738643928389281331366095207532193114305843094745273866336452296051523189325911545359997748735545153365459326965824504140992990653748551563388621812295143879168563643028714331800616525107242938555105725864904443179487735448344512225826058337229629080065151365167343274055473303802181974301995992231288914893936583978868212030806302235382565661039319775858591452481276532369619717771944703319704609243884667963269698910969623461539798686254883/741993016396956735405613435356683782697512283069823991064277931426980443824577566664069913413257621498538182516652902989918523185791709453919677031772399696514284422802749913157553979575003207665736459377430450383919221006176370669384465219743309411041792052666184021088982010598590366702909953828519058138526150106114448075948267632468571674052634702160284009616244631715592753300091481834250539282210722819840060174111768646004068650228726292470427652007163617929736425185577083563672618801965072487519722221379722725217400914053294087849015460578362361120445241591716408147719927671948668690922637069198489929597905348227031903643587457869031130741184797205245991710558088926873306205286102419574747370548119522443416892882892800

数值
> evalf[15](s);
.00444629883931381
不知道是否是你期待的

5楼2015-01-11 19:09:55

duyang2012

铜虫 (小有名气)

应助: 1 (幼儿园)
金币: 3.3
散金: 141
红花: 2
帖子: 131
在线: 119.5小时
虫号: 2257395
注册: 2013-01-23
专业: 通信理论与系统

引用回帖:

5楼: Originally posted by pippi6 at 2015-01-11 19:09:55
我用maple 7
> f:=BesselI(0, alpha*x)*alpha*exp(-alpha^2)*x*exp(-(x^2+alpha^2)/2)
> *(1-exp(-x^2/2))^1022;

> s:=int(int(f,alpha=0..infinity),x=0..infinity);
s := 329912268758474513720 ...

非常感谢你的回复，你这个应该是解决了我的难题，万分感谢。。。我没用过maple，请问你最后那个“evalf[15](s); .00444629883931381”是什么含义？！谢谢

6楼2015-01-11 22:03:48