24小时热门版块排行榜    

查看: 1042  |  回复: 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的回帖
相关版块跳转 我要订阅楼主 黑色裂变 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 332求调剂 +3 凤凰院丁真 2026-03-20 3/150 2026-03-21 10:27 by luoyongfeng
[考研] 306求调剂 +4 chuanzhu川烛 2026-03-18 4/200 2026-03-21 08:25 by laoshidan
[考研] 316求调剂 +6 梁茜雯 2026-03-19 6/300 2026-03-21 06:32 by Ecowxq666!
[考研] 310求调剂 +3 baibai1314 2026-03-16 3/150 2026-03-21 03:56 by JourneyLucky
[考研] 一志愿中国石油大学(华东) 本科齐鲁工业大学 +3 石能伟 2026-03-17 3/150 2026-03-21 02:22 by JourneyLucky
[考研] 华东师范大学-071000生物学-293分-求调剂 +3 研究生何瑶明 2026-03-18 3/150 2026-03-21 01:30 by JourneyLucky
[考研] 一志愿华中科技大学,080502,354分求调剂 +5 守候夕阳CF 2026-03-18 5/250 2026-03-21 01:06 by JourneyLucky
[考研] 一志愿南京理工大学085701资源与环境302分求调剂 +4 葵梓卫队 2026-03-18 6/300 2026-03-20 23:02 by JourneyLucky
[考研] 350求调剂 +5 weudhdk 2026-03-19 5/250 2026-03-20 22:04 by luoyongfeng
[考研] 求调剂 +3 @taotao 2026-03-20 3/150 2026-03-20 19:35 by JourneyLucky
[考研] 298-一志愿中国农业大学-求调剂 +9 手机用户 2026-03-17 9/450 2026-03-20 14:24 by 无懈可击111
[考研] 梁成伟老师课题组欢迎你的加入 +9 一鸭鸭哟 2026-03-14 11/550 2026-03-19 17:22 by !本暗一次!
[考研] 0703化学调剂 +4 18889395102 2026-03-18 4/200 2026-03-19 16:13 by 30660438
[考研] 085600材料与化工求调剂 +6 绪幸与子 2026-03-17 6/300 2026-03-19 13:27 by houyaoxu
[考研] 化学工程321分求调剂 +15 大米饭! 2026-03-15 18/900 2026-03-18 14:52 by haxia
[考研] 293求调剂 +11 zjl的号 2026-03-16 16/800 2026-03-18 08:10 by zhukairuo
[考博] 26博士申请 +3 1042136743 2026-03-17 3/150 2026-03-17 23:30 by 轻松不少随
[考研] 301求调剂 +4 A_JiXing 2026-03-16 4/200 2026-03-17 17:32 by ruiyingmiao
[考研] 26考研一志愿中国石油大学(华东)305分求调剂 +3 嘉年新程 2026-03-15 3/150 2026-03-15 13:58 by 哈哈哈哈嘿嘿嘿
[考研] 288求调剂 +4 奇点0314 2026-03-14 4/200 2026-03-14 23:04 by JourneyLucky
信息提示
请填处理意见