版塊導(dǎo)航: 正在加載中...

登錄注冊

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號能夠正常使用，請盡快對帳號進(jìn)行手機(jī)號驗(yàn)證，感謝您的理解與支持！

24小時熱門版塊排行榜

返回列表

黑色裂變

銅蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 165.8
散金: 23
帖子: 80
在線: 25.9小時
蟲號: 4524196
注冊: 2016-03-20
專業(yè): 控制論中的數(shù)學(xué)方法

[求助] matlab計(jì)算組合數(shù)出錯已有1人參與

最近寫的一個程序要計(jì)算組合數(shù)，但用MATLAB自帶的nchoosek當(dāng)數(shù)值較大時計(jì)算非常慢，于是改用階乘除階乘的形式，但是當(dāng)計(jì)算結(jié)果較大的時候MATLAB直接把結(jié)果當(dāng)成發(fā)散的了，顯示NAN，請問有什么好的處理方法嗎？把nchoosek提速或者使結(jié)果不發(fā)散的都可以。貼上我的代碼：
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

回復(fù)此樓

» 猜你喜歡

物理學(xué)求調(diào)劑已經(jīng)有4人回復(fù)
找博導(dǎo) 已經(jīng)有4人回復(fù)
282分材料專業(yè)求調(diào)劑院校已經(jīng)有18人回復(fù)
一志愿清華深研院材料專碩294分，專業(yè)課111分，本科中南大學(xué)材料，有六級，有工作經(jīng)驗(yàn) 已經(jīng)有4人回復(fù)
288求調(diào)劑已經(jīng)有15人回復(fù)
材料與化工一志愿山大 321分求調(diào)劑已經(jīng)有7人回復(fù)
哈工大材料324求調(diào)劑已經(jīng)有3人回復(fù)
304求調(diào)劑（085602一志愿985）已經(jīng)有5人回復(fù)
環(huán)境調(diào)劑已經(jīng)有5人回復(fù)
材料與化工，291，求調(diào)劑已經(jīng)有9人回復(fù)

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

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

zhangyq01

新蟲 (小有名氣)

應(yīng)助: 1 (幼兒園)
金幣: 211.4
散金: 16
紅花: 2
帖子: 118
在線: 22.7小時
蟲號: 2731642
注冊: 2013-10-17
性別: GG
專業(yè): 金屬非晶態(tài)、準(zhǔn)晶和納米晶

【答案】應(yīng)助回帖

好像函數(shù)nchoosek可以計(jì)算。
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

贊一下

回復(fù)此樓

2樓2020-03-14 16:46:02

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主黑色裂變的主題更新