24小时热门版块排行榜    

查看: 1034  |  回复: 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
回复此楼

» 猜你喜欢

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

zhangyq01

新虫 (小有名气)

【答案】应助回帖

好像函数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
2楼2020-03-14 16:46:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 黑色裂变 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 347求调剂 +4 啊欧欧欧 2026-03-03 4/200 2026-03-04 13:03 by 无际的草原
[考研] 欢迎采矿、地质、岩土、计算机、人工智能等专业的同学报考 +8 pin8023 2026-02-28 11/550 2026-03-04 00:12 by 不吃鱼肉
[考研] 289求调剂 +7 BrightLL 2026-03-02 9/450 2026-03-03 16:02 by tgxtgxtgx9
[考研] 成绩276,专业代码0856求调剂 +7 小陈朵 2026-03-03 7/350 2026-03-03 15:58 by tgxtgxtgx9
[考研] 一志愿中科大能动297求调剂,本科川大 +3 邵11 2026-03-03 3/150 2026-03-03 14:50 by 新能源达人
[考研] 环境工程专硕307 一志愿211 四六级已过 求调剂 +3 ccc! 2026-03-03 3/150 2026-03-03 12:55 by 热情沙漠
[考研] 299求调剂 +5 kkcoco25 2026-03-02 9/450 2026-03-03 12:55 by 公瑾逍遥
[考研] 调剂 +5 13853210211 2026-03-02 7/350 2026-03-03 11:42 by 13853210211
[考研] 清华大学 材料与化工 353分求调剂 +5 awaystay 2026-03-02 6/300 2026-03-03 09:03 by houyaoxu
[考研] 考研复试调剂,过国家线的同学都可报名 +7 黑!在干嘛 2026-02-28 8/400 2026-03-03 08:39 by ww哈
[考研] 0854复试调剂 276 +5 wmm9 2026-03-01 7/350 2026-03-03 02:49 by xiadaiyang
[考研] 0856调剂 +10 刘梦微 2026-02-28 10/500 2026-03-02 23:42 by hcy618
[考研] 295求调剂。一志愿报考郑州大学化学工艺学硕,总分295分 +8 yl1 2026-03-02 9/450 2026-03-02 21:21 by sunny81
[考研] 材料复试调剂 +5 学材料的点 2026-03-01 6/300 2026-03-02 20:01 by hypershenger
[考研] 一志愿东北大学材料专硕328,求调剂 +3 shs1083 2026-03-02 3/150 2026-03-02 17:27 by houyaoxu
[考研] 0856化工专硕求调剂 +15 董boxing 2026-03-01 15/750 2026-03-02 15:06 by 晃晃不许晃
[基金申请] 此成果不能导入原因:元数据必填信息不完整,可 进行补充。 +4 Kittylucky 2026-03-02 5/250 2026-03-02 11:07 by jurkat.1640
[考研] 化工专硕342,一志愿大连理工大学,求调剂 +6 kyf化工 2026-02-28 7/350 2026-03-02 10:56 by 无际的草原
[考研] 264求调剂 +4 巴拉巴拉根556 2026-02-28 4/200 2026-03-02 10:48 by yuchj
[考研] 调剂 +3 简木ChuFront 2026-02-28 3/150 2026-03-01 11:46 by 王伟要上岸啊
信息提示
请填处理意见