概率论实验报告

1、概率论与数理统计实验报告专业班级：XXX* ：xx* : XX日期ZXXXX、实验目的通过Matlab编程实验将抽象的理论转化为具体的图像，以便更 好的理解和记忆这些理论的内涵并将其应用于实践。二、实验内容及结果1.设 X N (呻2)；当 p = 1.5q = 0.5 时, 求 P1.8 X 2.9 ,P-2.5 1.6；当 p = 1.5q = 0.5 时，若 P X x = 0.95，求 ；分别绘制p = 1,2,3 , b = 0.5时的概率密度函数图形。解答：(1)源程序：clc;p1 = normcdf(2.9,1.5,0.5)-normcdf(1.8,1.5,0.5)p2=1-n

2、ormcdf(-2.5,1.5,0.5)p3=normcdf(0.1,1.5,0.5)+1-normcdf(3.3,1.5,0.5)运行结果：实验结论：尸1.8 x 2.9 =0.2717 ；P-2.5 1.6=0.0027。源程序：clc;*=0;p=normcdf(*,1.5,0.5);while(p0.95)*=*+0.001;p=normcdf(*,1.5,0.5);endP*运行结果：实验结论：此时*应为2.3230。(3)源程序：clc;clf;*=linspace(-1,5,1000); %(-1,5)等分为 1000 份p1=normpdf(*,1,0.5);p2=normpd

3、f(*,2,0.5);p3=normpdf(*,3,0.5);plot(*)p1),r,)*)p2),g,)*)p3,y,); %红色线表示 u=1，绿色线表示 u=2,黄色线表示u=3legend(,u=1,u=2,u=3,); %图线标记运行结果：已知每百份报纸全部卖出可获利14元，卖不出去将赔8元,设报纸的需求量x的分布律为X012345p0.05 0.10 0.25 0.35 0.150.10试确定报纸的最佳购进量。(要求使用计算机模拟)解答：源程序：clc;%假设报纸销售与购买均以百份为基本单位，不存在每百份中销售一部分、乘g余一部分的情况d=zeros(1,6);%用数组存储报纸销

4、售情况s=zeros(1,5);%s表示不同购进量下的盈利for(n=1:5)%至少应购进1的报纸(百份)，至多5,按照不同的购进量分别模拟规定次数的销售状况进行比较for(i = 1:365)%模拟一年的销售状况，也可以改变天数*=unifrnd(0,1);%模拟每日报纸销售量(百份)if(*0.05)%售出 0d(1)=d(1)+1;s(n)=s(n)-8*n;elseif(*0.15)%1d(2)=d(2)+1;s(n)=s(n)+14*1-8*(n-1);elseif(*0.4)%2d(3)=d(3)+1;if(n2)s(n)=s(n)+14;else s(n)=s(n)+14*2-8

5、*(n-2);endelseif(*0.75)%3d(4)=d(4)+1;if(n3)s(n)=s(n)+14*n;elses(n)=s(n)+14*3-8*(n-3); endelseif(*0.9)%4d(5)=d(5)+1;if(n4)s(n)=s(n)+14*n;elses(n)=s(n)+14*4-8*(n-4);endelse %5d(6)=d(6)+1;if(n5)s(n)=s(n)+14*n;elses(n)=s(n)+14*5;endendendendds运行结果：实验结论：由模拟结果可知，n=300时，收益最大为10666元, 故应取最佳购进量为300份。蒲丰投针实验取一大

6、白纸，在上面画出多条间距为d的平行直线，取一长度为 r (rd)的针，随机投到纸上n次，记针与直线相交的次数为 m.由此实验计算针与直线相交的概率。圆周率的近似值。解答：源程序：clc;d=2;%平行线间距r=1;%针长n = 10000;%试验次数m=0;%存储相交次数 for(i = 1:n)s=unifrnd(0,d/2);q = unifrnd(0,pi);if(s = r/2*sin(q)m = m + 1;endendp=m/n%针与直线相交的概率pai = 1/p%n(2rn)/(dm)运行结果：实验结论：针与直线相交的概率为0.3147 ；圆周率的近似值为3.1776.5. 5

7、.设* B(n, p),其中 np=2，对 n = 10,102111010，讨论用泊松 分布逼近二项分布的误差，画出逼近的图形。解答：源程序：clear;clc;k = 0:10;n = 10;p=0.2;lamda = n*p;B = binopdf(k,n,p);P = poisspdf(k,lamda);subplot(2,2,1)plot(k,B-o)title(二项分布)grid onsubplot(2,2,2)plot(k)P)-*)color,0 0.5 0)泊松分布)grid onsubplot(2,2,3)plot(k)abs(B-P),-r)绝对误差)grid onset

8、(gca1color,0.23110.44310.337)subplot(2,2,4)plotCk.B.-Q.k.P,-*)title(二项分布与泊松分布)legend(二项分布泊松分布)grid onset(gca1color,0.72910.83110.957)运行结果：实验结论：N越大，泊松分布的图像越逼近二项分布。6.对 n=102,10101 计算P(5 X 50) , P(20 x 90)1）用二项分布计算2）用泊松分布计算3）用正态分布计算比较用泊松分布逼近与正态分布逼近二项分布的优劣。解答：源程序：clc;clear;n = 10;y1=zeros(1,9);y2=zeros(

9、1,9);y3=zeros(1,9);y4=zeros(1,9);for i = 1:9n = n*10;y1(i)=binocdf(50)n)2/n)-binocdf(5)n,2/n);y2(i)=binocdf(90)n)2/n)-binocdf(20)n,2/n);y3(i)=normcdf(5012)sqrt(2*(1-2/n)-normcdf(512,sqrt(2*(1-2/n);y4(i)=normcdf(90)2)sqrt(2*(1-2/n)-normcdf(20)2,sqrt(2*(1-2/n);endy5=poisscdf(50,2)-poisscdf(512);y6=poisscdf(9012)-poisscdf(20,2);disp(5P=50 二项分布);disp(yl);disp(5P=50 正态分布);disp(y3);disp(5P=50 泊松分布);disp(y5);disp(20P=90 二项分布);disp(y2);disp(20P=90 正态分布);disp(y4);disp(20P=90 泊松分布);disp(y6);运行结果：实验结论：计算结果如上图所示。由图可得，p值较小时，泊松分布 与正态