统计计算实验报告.doc

文档介绍：-
. z.
- .
- - 考试资料.
实验步骤：
复****第四章有关知识，熟悉Buffon投针试验的原理及各种求定积分的方法
运用MATLAB，编写相关程序，输入代码
观察实验结果，进展讨论
撰写实验报告
实验内容〔包括实验具体内容、算法分析、源代码等等〕：
Buffon投针试验求圆周率：
function y=Buffon(a,l,N) //a为平行线间的距离；l为针长；N为投针次数
M=0;
i=0;
while i<=N
*i=rand()*a/2;
yi=rand()*pi;
if *i<=l/2*sin(yi)
M=M+1;
end
i=i+1;
end
phi=2*l*N/(a*M);
y=phi;
end
随机模拟计算定积分
随机投点法：
function y=R(N)
n=0;
m=0;
while n<=N
kesi=rand();
eta=rand();
if eta<=(e*p(-1+2*kesi)-e*p(-1))/(e*p(1)-e*p(-1))
m=m+1;
end
n=n+1;
end
y=2*(e*p(1)-e*p(-1))*m/N+2*e*p(-1)
end
平均值估计法:
function y=avge(N)
i=0;
s=0;
while i<=N
-
. z.
- - .
- - 考试资料.
r=rand();
f=(e*p(-1+2*r)-e*p(-1))/(e*p(1)-e*p(-1));
s=s+f;
i=i+1;
end
zeta=1/N*s;
y=2*(e*p(1)-e*p(-1))*zeta+2*e*p(-1);
End
重要抽样法：
function y=zhongyao(N) //第一步
i=0;
s=0;
while i<=N
r=rand();
*