实例解析,【例12-3】设袋中有10个球，其中3个是红球，7个是白球，从中取5个球，设X表示取到红球的个数，求X的概率分布。,解：易知X服从参数为10,3,5的超几何分布，故编写如下简单程序： for k=1:4 p(k)=hygepdf(k-1,10,3,5); % 生成服从参数为10，3，5的超几何分布随机数 end p,运行结果： p =0.0833 0.4167 0.4167 0.0833,下面利用蒙特卡罗方法求解该例，编写如下语句： rand(state,2009) % 设置随机数状态 s1=0;s2=0;s3=0;s4=0; % 设置计数器 N=50000; % 模拟次数 for i=1:N x=randperm(10); % 产生1：10的一个排列 if sum(x(1:5)=3)=0 % 取到0个红球的情形 s1=s1+1; elseif sum(x(1:5)=3)=1 % 取到1个红球的情形 s2=s2+1; elseif sum(x(1:5)=3)=2 % 取到两个红球的情形 s3=s3+1; else % 取到3个红球的情形
