程序
Redball <- function(M, m, N) { ball = c(rep(1,m), rep(0,M-m)) ball = sample(ball, M) List = list() order = 1:M k = 0 for(i in 1:(M/N)) { B = sample(order,N) order = order[-B] List[[i]] = ball[B] if(max(ball[B])==1)k = k + 1 } list(组 = List, 含红球组数 = k) } pro.Redball <- function(M, m, N, n) #n为实验重复次数 { g = M/N + 1 b = rep(0,g) for(i in 1:n) { num = Redball(M,m,N)[[2]] + 1 b[num] = b[num] + 1 } p = b/n g = g - 1 barplot(p, names.arg = paste(0:g), col = rgb(0.3,0.5,0.7)) p } pro.Redball(24, 5, 6, 1000)
注:此处我和题主有一对符号概念互换了,但是我将错就错了,不想改了。声明: 为组数, 为一组球的个数。
设组数 , 表示第 组红球的个数, 表示第 组第 个球是否是红球,若是则 ;否则为 ,则满足不定方程:
则该方程解的全体 是全空间,显然有
而我们研究的 事件空间对应的不定方程:
设它的解集 ,
……
所求概率即为
其中