24小时热门版块排行榜    

查看: 1039  |  回复: 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的回帖
相关版块跳转 我要订阅楼主 黑色裂变 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 283求调剂,材料、化工皆可 +5 苏打水7777 2026-03-11 5/250 2026-03-12 10:44 by vgtyfty
[考研] 一志愿北京化工大学材料与化工296分求调剂 +15 稻妻小编 2026-03-09 17/850 2026-03-11 22:45 by 2020015
[考研] 材料与化工(0856)304求B区调剂 +6 邱gl 2026-03-10 9/450 2026-03-11 16:37 by @飒飒飒飒
[考研] 接受26届调剂生 +22 猪猪猪毛 2026-03-06 23/1150 2026-03-11 14:55 by Lly031128
[考研] 一志愿浙江大学0856材料与化工求调剂 +3 yansheng@211 2026-03-09 4/200 2026-03-11 10:52 by @飒飒飒飒
[考研] 08工科 +5 li李乐成 2026-03-06 5/250 2026-03-10 14:51 by 在风落中
[考研] 调剂的同学,走过路过,不要错过 +11 likeihood 2026-03-06 16/800 2026-03-10 11:38 by Djdjj12
[考研] 复试调剂 +6 呼呼?~+123456 2026-03-08 8/400 2026-03-10 11:36 by 呼呼?~+123456
[考研] 085600材料与化工 326 求调剂 +4 热爱生活ing 2026-03-09 4/200 2026-03-10 08:25 by Linda Hu
[考研] 294 英二数二物化 求调剂 +6 米饭团不好吃 2026-03-09 6/300 2026-03-09 23:55 by barlinike
[基金申请] 面上项目还需要AI说明吗? +3 liyundong 2026-03-08 3/150 2026-03-09 22:30 by kingkocxr
[考研] 数二英二309分请求调剂 +3 dtdxzxx 2026-03-09 4/200 2026-03-09 19:56 by yuningshan
[考研] 334求调剂 +5 Trying] 2026-03-07 5/250 2026-03-09 14:47 by ruiyingmiao
[考研] 313求调剂 +4 Yyt杨1 2026-03-07 5/250 2026-03-09 13:48 by macy2011
[考研] 求调剂,一志愿江南大学,食品科学与工程,总分,320 +3 yyyyyukino 2026-03-07 3/150 2026-03-08 23:07 by 清风月
[考研] 269求调剂 +3 朔朔话 2026-03-08 4/200 2026-03-08 20:39 by 热情沙漠
[考研] 288求调剂(0703)一志愿东北大学 +5 好好- 2026-03-07 5/250 2026-03-08 18:22 by zzpnuli111
[考研] 286求调剂 +12 Faune 2026-03-06 14/700 2026-03-08 14:54 by hmn_wj
[考研] 303求调剂 +8 forgman95 2026-03-05 10/500 2026-03-08 12:41 by 蓝莓都给你
[考研] 304求调剂 +4 52hz~~ 2026-03-05 5/250 2026-03-07 15:47 by lature00
信息提示
请填处理意见