《概率论与数理统计上机实习报告》由会员分享，可在线阅读，更多相关《概率论与数理统计上机实习报告（5页珍藏版）》请在金锄头文库上搜索。
1、概率论与数理统计上机实习报告姓名：学号：学部：管理与经济学部专业：班级：一、某人写了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(t
2、ime(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#inc
3、ludeusingnamespacestd;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;
4、部分截图如下：输入命令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;得到结果：
