概率论上机

发布 2022-10-11 15:21:28 阅读 4476

概率论实验报告。

班级:医电21

学号:2121201019

姓名:闫创。

实验目的:通过使用matlab,实现数学概论问题的求解,达到精确计算,减小计算量。

实验问题:1. 设~;

1) 当时,求,,;

2) 当时,若,求;

3) 分别绘制, 时的概率密度函数图形。

p1=normcdf(2.9,1.5,0.5)-normcdf(1.8,1.5,0.5)

p2=1-normcdf(-2.5,1.5,0.5)

p3=normcdf(0.1,1.5,0.5)+1-normcdf(3.3,1.5,0.5)

x=0;p=normcdf(x,1.5,0.5);

while(p<0.95)

x=x+0.001;

p=normcdf(x,1.5,0.5);endp

xx=linspace(-1,5,1000); 1,5)等分为1000份。

p1=normpdf(x,1,0.5);

p2=normpdf(x,2,0.5);

p3=normpdf(x,3,0.5);

plot(x,p1,'y',x,p2,'g',x,p3,'r');红色线表示u=1,绿色线表示u=2,黄色线表示u=3

legend('u=1','u=2','u=3');图线标记。

2. 已知每百份报纸全部卖出可获利14元,卖不出去将赔8元,设报纸的需求量的分布律为

试确定报纸的最佳购进量。(要求使用计算机模拟)

2、clc; %假设报纸销售与购买均以百份为基本单位,不存在每百份中销售一部分、剩余一部分的情况。

d=zeros(1,6); 用数组存储报纸销售情况。

s=zeros(1,5); s表示不同购进量下的盈利。

for(n=1:5) %至少应购进1的报纸(百份),至多5,按照不同的购进量分别模拟规定次数的销售状况进行比较。

for(i=1:365) %模拟一年的销售状况,也可以改变天数。

x=unifrnd(0,1); 模拟每**纸销售量(百份)

if(x<0.05) %售出0

d(1)=d(1)+1;

s(n)=s(n)-8*n;

elseif(x<0.15) %1

d(2)=d(2)+1;

s(n)=s(n)+14*1-8*(n-1);

elseif(x<0.4) %2

d(3)=d(3)+1;

if(n<2)

s(n)=s(n)+14;

elses(n)=s(n)+14*2-8*(n-2);

endelseif(x<0.75) %3

d(4)=d(4)+1;

if(n<3)

s(n)=s(n)+14*n;

elses(n)=s(n)+14*3-8*(n-3);

endelseif(x<0.9) %4

d(5)=d(5)+1;

if(n<4)

s(n)=s(n)+14*n;

elses(n)=s(n)+14*4-8*(n-4);

endelse %5

d(6)=d(6)+1;

if(n<5)

s(n)=s(n)+14*n;

elses(n)=s(n)+14*5;

endend

endendds

3. 蒲丰投针实验。

取一张白纸,在上面画出多条间距为d的平行直线,取一长度为r(r1) 针与直线相交的概率。

2) 圆周率的近似值。

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;

endend

p=m/n %针与直线相交的概率。

pai=1/p %π2rn)/(dm)

证明。下面就是一个简单而巧妙的证明。

找一根铁丝弯成一个圆圈,使其直径恰恰等于平行线间的距离d。

可以想象得到,对于这样的圆圈来说,不管怎么扔下,都将和平行线有两个交点。

因此,如果圆圈扔下的次数为n次,那么相交的交点总数必为2n。

现在设想把圆圈拉直,变成一条长为πd的铁丝。

显然,这样的铁丝扔下时与平行线相交的情形要比圆圈复杂些,可能有4个交点,3个交点,2个交点,1个交点,甚至于都不相交。

由于圆圈和直线的长度同为πd,根据机会均等的原理,当它们投掷次数较多,且相等时,两者与平行线组交点的总数期望也是一样的。

这就是说,当长为πd的铁丝扔下n次时,与平行线相交的交点总数应大致为2n。

现在转而讨论铁丝长为l的情形。当投掷次数n增大的时候,这种铁丝跟平行线相交的交点总数m应当与长度l成正比,因而有:m=kl,式中k是比例系数。

为了求出k来,只需注意到,对于l=πd的特殊情形,有m=2n。于是求得k=(2n)/(d)。代入前式就有:m≈(2ln)/(d)从而π≈(2ln)/(dm)

概率论上机作业

西安交通大学实验报告。课程名称 概率论与数理统计实验名称 蒙特卡洛方法估计积分。学院理学院实验日期 2014 年 12 月 14日。班级 应化31 姓名 张迟学号 2130904023 一 实验内容。1.用蒙特卡洛方法估计积分 和的值,并将估计值与真值进行比较。2.用蒙特卡洛方法估计积分和的值,并对...

概率论上机作业

学部。班级。学号。姓名。概率论与数理统计 上机题目。1.假设车辆牌号的后两位数都是00 99中的某一个,任取m辆车 编制程序求出这m辆车中至少有两辆车后两位数相同的概率,画出m 10,11,60时,的散点图。2.某人写了封信,又写了个信封,然后将这封信随机地装入这个信封中,用表示至少有一封信装对的概...

概率论上机实验报告

假如某保险公司有10000个同阶层的人参加人寿保险,每人每年付1200元保险费,在一年内投保人的死亡概率为0.006,死亡时保险公司需赔付10万元。试问 1 保险公司一年的总赔偿金在590万元至610万元间的概率是多少?2 保险公司亏本的概率有多大?3 保险公司每年利润大于400万元的概率是多少?设...