| 查看: 1223 | 回复: 1 | ||
[求助]
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一区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人回复
-
某top大学教授说“能够在市场中兑现的能力才是真能力”无比同意！
已经有13人回复
-
同一个分子，一条来自文献，一条来自AI——不告诉你答案，你会选哪条？
已经有7人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有4人回复
2楼2020-03-14 16:46:02
