24小时热门版块排行榜    

CyRhmU.jpeg
查看: 1514  |  回复: 10
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

握瑜怀瑾

木虫 (小有名气)

[求助] 一个简单的式子,却算不正确

各位高手大家好!

目前我遇到一个问题,就是附件文档中所示的式子,当K比较大时,MATLAB中的nchoosek函数就会报错;因此我采用了mathematica来计算

在mathematica中,当K=50时,基本能算正确;当K=100时就悲剧了,算出了一个很大的值,是显然错误的

求教:当K比较大时,比如说100,200,甚至更大时,有什么办法能算出这个式子的准确值?

[ Last edited by jjdg on 2011-9-29 at 12:44 ]
回复此楼

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

通信人生

» 猜你喜欢

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

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

inhaul

新虫 (正式写手)

【答案】应助回帖

握瑜怀瑾(金币+2): 呵呵,我也发现了,K=10时,是正确的 2011-10-10 08:18:47
引用回帖:
4楼: Originally posted by gsview at 2011-09-29 15:05:06:
1,      0.5963473623231725
10,    3.33204505399828
50,    5.383752988862083
100,  6.225487226306522
200,  7.048564754038644
直接做的数值积分,不知道结果对吗?

算了一下,你的结果好像不对,除了K=10是正确的
代码
K=50;
A=10;
f=NIntegrate[Log[1+A*x]*(1-Exp[-x])^(K-1)*Exp[-x],{x,0,Infinity}]*K

K   f
1,2.01464
10,3.33205
50,3.79314
100,3.94101
200,4.06953
10楼2011-10-10 00:16:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 11 个回答

握瑜怀瑾

木虫 (小有名气)

哦,貌似标题错了,应该是“却不能算正确”
3楼2011-09-29 09:06:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gsview

木虫 (小有名气)

【答案】应助回帖


握瑜怀瑾(金币+2): 谢谢您! 2011-10-05 10:48:13
xzhdty(金币+1): 欢迎常来程序语言看看 2011-10-05 18:48:05
1,      0.5963473623231725
10,    3.33204505399828
50,    5.383752988862083
100,  6.225487226306522
200,  7.048564754038644
直接做的数值积分,不知道结果对吗?
4楼2011-09-29 15:05:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

握瑜怀瑾

木虫 (小有名气)

引用回帖:
4楼: Originally posted by gsview at 2011-09-29 15:05:06:
1,      0.5963473623231725
10,    3.33204505399828
50,    5.383752988862083
100,  6.225487226306522
200,  7.048564754038644
直接做的数值积分,不知道结果对吗?

您所说的“直接做的数值积分”,是怎么个计算法啊?在MATLAB还是Mathematica里面算的呢?能说得详细点吗?是不是未采用推导出来的求和公式计算?
5楼2011-10-05 10:49:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
信息提示
请填处理意见