| 查看: 1197 | 回复: 5 | |||
[交流]
蒙特卡洛法 求积分 被积函数含有gamma函数 已有3人参与
|
|
为了求解图片1中的tmax积分函数，由于被积函数中含有不完全gamma函数，故采用蒙特卡洛法。 由于积分区间是0-inf,所以重要性密度函数选择了gamma函数，以下是相关语句，采用不同的重要性密度函数和粒子数目，仿真结果就差好多，这个积分该如何求解呢？ clear clc N=100000; x=gamrnd(1,1,N,1);%生成gamma分布随机数 for i=1:N y(i)=(1-((1-gammainc(4,2*x(i)))/gamma(2*x(i))))/gampdf(x(i),1,1); %随机数带入被积函数 end D=mean(y)%求均值  1.jpg  2.jpg |
回复此楼» 猜你喜欢
焦虑
已经有8人回复
-
308求调剂
已经有4人回复
-
NSFC申报书里申请人简历中代表性论著还需要在申报书最后的附件里面再上传一遍吗
已经有14人回复
-
材料与化工一志愿南昌大学327求调剂推荐
已经有6人回复
-
化学调剂0703
已经有7人回复
-
327求调剂
已经有11人回复
-
调剂
已经有8人回复
-
梁成伟老师课题组欢迎你的加入
已经有7人回复
-
伙伴们，祝我生日快乐吧
已经有24人回复
-
中科院材料273求调剂
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求一个定积分的原函数。。。求说明过程。。不会算。。
已经有5人回复
- 特殊函数求解高手请进，这个类似于Gamma函数的定积分能否求解？
已经有12人回复
- 分布函数积分求助~
已经有4人回复
- 【求助】这个概率密度函数正确吗？
已经有5人回复
sciencejoy
新虫 (著名写手)
- 应助: 436 (硕士)
- 金币: 11207.1
- 红花: 90
- 帖子: 2974
- 在线: 498小时
- 虫号: 802149
- 注册: 2009-07-02
- 性别: GG
- 专业: 高分子物理与高分子物理化
2楼2013-11-12 18:39:33
3楼2013-11-13 21:13:23
锐利的碎片
木虫 (正式写手)
star watcher
- 应助: 136 (高中生)
- 金币: 3637.1
- 散金: 252
- 红花: 22
- 帖子: 988
- 在线: 1224.9小时
- 虫号: 961933
- 注册: 2010-03-05
- 专业: 凝聚态物性 II ：电子结构
4楼2013-11-13 22:21:26
yaozhq
金虫 (小有名气)
- 应助: 126 (高中生)
- 金币: 1137.8
- 红花: 7
- 帖子: 260
- 在线: 37.6小时
- 虫号: 1475454
- 注册: 2011-11-04
- 性别: GG
- 专业: 生物大分子结构与功能
5楼2013-11-14 11:21:09
6楼2013-11-24 10:40:29
