概率论上机实验报告

西安交通大学

概率论与数理统计

实验报告

一． 实验问题

用蒙特卡方法估计积分值

1.估计.

的值并将估计值与真值进行比较

2.估计.

的值并对误差进行估计。

二．问题分析

蒙特卡罗方法求积分的一般规则如下：任何一个积分，都可看作某个随机变量的期望值，因此，可以用这个随机变量的平均值来近似它。

设欲求积分

其中，P＝P(x1，x2，…，xs) 表示 s 维空间的点，Vs表示积分区域。取Vs上任一联合概率密度函数 f (P)，令

则

即θ是随机变量 g(P) 的数学期望，P的分布密度函数为 f (P) 。现从 f (P) 中抽取随机向量 P 的 N 个样本：Pi，i＝1，2，…，N，则

就是θ的近似估计。

因此，我们可以利用matlab产生（0，1）区间均匀分布随机数以及标准正态分布随机数的指令得到理想的随机数。并对它进行一定的数学变换，即可得到需要的样本，进而得到积分值，以及对结果的误差估计

三． 程序设计

程序如下:

clc;

clear;

m=10;

n=10000;

d=0;

e=0;

for i=1:m

d=0;

a=rand(1,n);

for j=1:n

b=a(j)+2;

c=b^2;

e(i)=d+c/n;

d=e(i);

end

fprintf('e=%.8f\n',e(i))

end

p=sum(e)/m;

for j=1:m;

s(j)=(e(j)-p)^2;

end

q=sum(s);

fprintf('x=%.8f\ny=%.8f\n',p,q);

程序如下：

clc;

clear;

m=10;

n=10000;

d=0;

e=0;

for i=1:m

d=0;

a=rand(1,n);

for j=1:n

b=a(j)*3;

c=b*sin(b)*3

e(i)=d+c/n;

d=e(i);

end

fprintf('e=%.8f\n',e(i))

end

p=sum(e)/m;

for j=1:m;

s(j)=(e(j)-p)^2;

end

q=sum(s);

fprintf('x=%.8f\ny=%.8f\n',p,q);

程序如下：

clc;

clear;

m=10;

n=10000;

d=0;

e=0;

for i=1:m

d=0;

a=rand(1,n);

for j=1:n

b=exp(-a(j)^2/2);

c=b/n*(2*pi)^0.5;

e(i)=d+c/2;

d=e(i);

end

fprintf('e=%.8f\n',e(i))

end

p=sum(e)/m;

for j=1:m;

s(j)=(e(j)-p)^2;

end

q=sum(s);

fprintf('x=%.8f\ny=%.8f\n',p,q);

程序如下：

clc;

clear;

m=10;

n=10000;

d=0;

e=0;

for i=1:m

d=0;

a=rand(1,n);

for j=1:n

b=exp(a(j)^2);

c=b/n;

e(i)=d+c;

d=e(i);

end

fprintf('e=%.8f\n',e(i))

end

p=sum(e)/m;

for j=1:m;

s(j)=(e(j)-p)^2;

end

q=sum(s);

fprintf('x=%.8f\ny=%.8f\n',p,q);

程序如下：

clc;

clear;

m=10;

n=10000;

d=0;

e=0;

for i=1:m

d=0;

a=rand(1,n);

for j=1:n

b=a(i)*4;

c=1/((1+b^2)^0.5);

e(i)=d+c*4/n;

d=e(i);

end

fprintf('e=%.8f\n',e(i))

end

p=sum(e)/m;

for j=1:m;

s(j)=(e(j)-p)^2;

end

q=sum(s);

fprintf('x=%.8f\ny=%.8f\n',p,q);

四．实验结果

e=6.34879520

e=6.34068140

e=6.35081124

e=6.31353632

e=6.35586630

e=6.33058791

e=6.32419121

e=6.33707454

e=6.30357011

e=6.35063255

x=6.33557468

y=0.00273819

e=3.12211717

e=3.11373037

e=3.07484948

e=3.08660758

e=3.10052243

e=3.10475698

e=3.13762746

e=3.16481618

e=3.11552000

e=3.09615989

x=3.11167076

y=0.00602301

e=0.88617655

e=0.88538972

e=0.88635209

e=0.88575809

e=0.88653705

e=0.88606366

e=0.88634011

e=0.88613926

e=0.88573325

e=0.88644154

x=0.88609313

y=0.00000119

e=1.46211146

e=1.46154792

e=1.46327379

e=1.46256348

e=1.46318297

e=1.46235828

e=1.46241378

e=1.46316145

e=1.46203052

e=1.46280489

x=1.46254485

y=0.00000288

e=1.98511173

e=1.02167881

e=1.26713031

e=0.98770837

e=1.14216662

e=1.75642022

e=1.97055988

e=1.96227794

e=1.83229787

e=1.06190231

x=1.49872541

y=1.71205432

其中e为实验值，x为平均值，y为方差。

五实验总结与体会

通过这次实验，我掌握了利用计算机软件通过蒙特卡方法估计积分值的方法；体会到了概率论与数理统计这一学科在积分领域的独特作用以及估计的基本思想，加深了我对书本理论知识的理解；同时使我的matlab编程技巧有了一定的提高。

在实验过程中，不够注重程序语言的精准性，经常被一个小小的标点符号而影响了实验结果。在今后的实验中我要吸取教训。做到细心，谨慎。

第二篇：概率统计上机实验报告(电子版)

2.(1)BINOMDIST(2,15,0.05,FALSE)=0.13475

BINOMDIST(2,15,0.05,TRUE)=0.9638

(2)EXPONDIST(1,0.1,FALSE)=0.09048

EXPONDIST(4,0.1,TRUE)=0.32968

(3)NORMDIST(2,0,1, TRUE)=0.97725

NORMSDIST(2)-- NORMSDIST(--2)=0.9545

=NORMINV(0.98,0,1)=2.05

NORMSDIST(0.1)-- NORMSDIST(--1)=0.3812

=NORMINV(0.05,5,100)=--159.49

(4)POISSON(4,2,FALSE)=0.090

POISSON(4,2,TRUE)=0.9473

(5) BINOMDIST(2,15,0.05,FALSE)=0.13475

营业税金与社会商品总额关系

(1)打开EXCEL,建立数据文件如下图:

调用线性回归分析程序:单击工具/数据分析/回归/确定,填写对话框,确定后输出结果,分析结果知回归方程为:

Y=-2.258+0.0487X

(2)对数据调用相关分析程序:依次单击工具/数据分析/相关系数/确定,填写对话框后,单击确定得到下面表格:

所以,Y与X的皮尔逊相关系数为: 0.981069

(3)建立假设H0:b=0 ,H1:b=/0,统计检验量F=(SSR/k)/(SSE/n-k-1)

有数据分析结果知：F=179.6507

P(F(1,7)>179.6507)=3.02E-06<<0.05

所以认为回归方程是显著有效的。

（4）在（1）中表的B11中补充数据X=320

在A11中输入公式=-2.258+0.0487X320

运行课的到X=320的点预测值y=13.326

相关推荐

概率论上机实验报告

第二篇：概率统计上机实验报告(电子版)

专栏推荐