资源预览内容
第1页 / 共5页
第2页 / 共5页
第3页 / 共5页
第4页 / 共5页
第5页 / 共5页
亲,该文档总共5页全部预览完了,如果喜欢就下载吧!
资源描述
概率论与数理统计上机实习报告姓名:学号:学部:管理与经济学部专业:班级:一、某人写了n封信,又写了n个信封,然后将这n封信随机地装入这n个信封中,用Pn表示至少有一封信装对的概率。1 .编制程序,用随机数模拟至少20000次,求当n=10时,Pn的值。2 .重复第一步,画由n=2,3,,50时,Pn的散点图。实验过程:1.编制程序,用随机数模拟至少20000次,求当n=10时,Pn的值。C+程序源代码如下:#include#include#include#includeusingnamespacestd;inta2000111;intn=0;intmain()(inti,j,k;srand(time(0);for(k=1;k=20000;k+)(for(i=1;i=10;i+)(intflag=0;while(flag=0)(aki=rand()%10+1;if(i=1)break;for(j=1;ji;j+)(if(aki=akj)flag=0;break;flag=1;for(i=1;i=10;i+)if(aki=i)n+;i=11;cout至少有一封信装对的事件数为n次.实验次数为20000次.概率为n/20000.0endl;return0;2.重复第一步,画由 n=2, 3,50 时, Pn 的散点图。程序代码如下:#include#include#include#includeusingnamespacestd;inta2000151;intb52=0;intmain()inti,j,k,l;srand(time(0);for(l=2;l=51;l+)for(k=1;k=20000;k+)for(i=1;i=l;i+)intflag=0;while(flag=0)aki=rand()%l+1;if(i=1)break;for(j=1;ji;j+)if(aki=akj)flag=0;break;flag=1;for(i=1;i=l;i+)if(aki=i)bl+;i=l+1;for(l=2;l=51;l+)coutbl/20000.0endl;return0;部分截图如下:输入命令fo门=2:51x(1,i-1)=iend将带入scatter(x,y,*)得到(x,Pn)散点图如下二、设X1,X2,Xn相互独立且都服从区间0,1上的均匀分布,f(x)n不为区间0,1上的一个可积函数,由大数定律可知n-依概率收敛于=,编制程序,用随机数模拟至少40000次,近似地求下列两个积分的值:实验过程:1 .在matlab下输入以下命令:y=(x)exp(xA2);sum=0;x=rand(1,40000);fori=1:40000sum=sum+y(x(1,i);endss=sum/40000;得到结果:2 .在matlab下输入以下命令:y=(x)sin(x)/x;sum=0;x=rand(1,40000);fori=1:40000sum=sum+y(x(1,i);endss=sum/40000;得到结果:
收藏 下载该资源
网站客服QQ:2055934822
金锄头文库版权所有
经营许可证:蜀ICP备13022795号 | 川公网安备 51140202000112号