版块导航: 正在加载中...

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜

返回列表

黑色裂变

铜虫 (小有名气)

应助: 0 (幼儿园)
金币: 165.8
散金: 23
帖子: 80
在线: 25.9小时
虫号: 4524196
注册: 2016-03-20
专业: 控制论中的数学方法

[求助] matlab计算组合数出错已有1人参与

最近写的一个程序要计算组合数，但用MATLAB自带的nchoosek当数值较大时计算非常慢，于是改用阶乘除阶乘的形式，但是当计算结果较大的时候MATLAB直接把结果当成发散的了，显示NAN，请问有什么好的处理方法吗？把nchoosek提速或者使结果不发散的都可以。贴上我的代码：
m=500;
w=10;
S1=1;
pk1=0.5;
pdw1=0.5;
Adp1=zeros(w+1,m+1);
for r=1:w+1
for c=1:r+m-w
      for i=0:c-1
         Adp1(r,c)=Adp1(r,c)+(prod(1:r+m-w-1)/(prod(1:i)*prod(1:r+m-w-1-i)))*(1-pdw1)^(i)*(pdw1)^(r-i-1+m-w)*(prod(1:r-i+m-w-1)/(prod(1:c-i-1)*prod(1:r-c+m-w)))*((1-pk1)^S1)^(c-i-1)*(1-(1-pk1)^S1)^(r-c+m-w);
      end
end
end

回复此楼

» 猜你喜欢

售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急已经有4人回复
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急已经有4人回复
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急已经有8人回复
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急已经有4人回复
售SCI文章，我:8O.5.5.1O.54,科目全,可十急已经有4人回复
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急已经有4人回复
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急已经有4人回复
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急已经有6人回复
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急已经有6人回复
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急已经有4人回复

1楼 2019-05-04 15:21:15

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

zhangyq01

新虫 (小有名气)

应助: 1 (幼儿园)
金币: 211.4
散金: 16
红花: 2
帖子: 118
在线: 22.7小时
虫号: 2731642
注册: 2013-10-17
性别: GG
专业: 金属非晶态、准晶和纳米晶

【答案】应助回帖

好像函数nchoosek可以计算。
c=nchoosek(500,10)
Warning: Result may not be exact. Coefficient is greater than 1.000000e+015 and is
only accurate to 15 digits
> In nchoosek at 66

c =

2.4581e+020

赞一下