sas课程论文(免费)(4)

武汉大学学习sas之后写的课程论文。。。。。

可以很大，譬如n=10000，甚至n=100000。

(2)对n对数据(xi,y)，i=1，2，…，n，记录满足如下不等式

i

y

的次数，这就是事件A发生的频数

i

f(xi)

n

。由此可得事件A发生的频率

n

n

，则J

n

n

。

以下我们用sas

计算

e

1

x

2

/2

，其精确值为0.341344。

在sas中运行如下程序：

data mtcl; N=0;

do i=1 to 100000;

x=uniform(1); /*产生随机数*/

y=uniform(1); /*产生随机数*/

if y<=exp(-x**2/2)/sqrt(6.28318530717959) then N=N+1; /*统计满足(2)中不等式的次数*/ J=N/100000; output; end; run;

程序运行之后得到一个数据集Work.Mtcl。从中我们可以看出随着i的增加，积分值J越来越接近精确值。由于数据集过长，只将最后一部分写在下边：