例3观察:用binornd模拟5000次投球过程,观察 球堆积的情况。 >clear;clf, 9n=5;p=0.5;m=5000;x=[0:1:n] rand('seed',3) =binornd(n,p,1,m);%模拟服从二项分布的随机数,相当于模拟投球m 次 for=1:n+1%开始计数 9 k=[]; k=find(R==(-1);%find是一个有用的指令,本语句的作用是找出R中等 于(仙-)元素下标,并赋予向量k中 9 h(=length(k)m;%计算落于编号(-)的格子中的小球频率 end bar(x,h),axis([-1601])%画频率图 9 title(八fontsize{18八fontname{华文新魏5000次投球小球堆积的频率图') 15/8715/87 例3 观察:用binornd模拟5000次投球过程,观察小 球堆积的情况。 >> clear;clf, n=5;p=0.5;m=5000;x=[0:1:n] rand('seed',3) R=binornd(n,p,1,m);%模拟服从二项分布的随机数,相当于模拟投球m 次 for I=1:n+1 %开始计数 k=[ ]; k=find(R==(I-1));%find是一个有用的指令,本语句的作用是找出R中等 于(I-1)元素下标,并赋予向量k中 h(I)=length(k)/m;%计算落于编号(I-1)的格子中的小球频率 end bar(x,h),axis([-1 6 0 1])%画频率图 title('\fontsize{18}\fontname{华文新魏}5000次投球小球堆积的频率图')