- INPUT PROGRAM.
- LOOP #I=1 to 200000.
- COMPUTE x=UNIFORM(1).
- COMPUTE y=x.
- END CASE.
- END LOOP.
- END FILE.
- END INPUT PROGRAM.
- EXECUTE.
- DESCRIPTIVES VARIABLES=y
- /SORT SUM
- /STATISTICS=mean.
[td]
| 描述统计量 | ||
|
| N | 均值 |
y | 200000 | .50 |
有效的 N (列表状态) | 200000 |
|
第二个例子,稍微复杂点。在x在[0,1]上和y[0,1]上求x^2+y^2的二重积分。自己算结果下来为2/3,。用程序怎么实现?试试看。
- INPUT PROGRAM.
- LOOP #I=1 to 200000.
- COMPUTE x=UNIFORM(1).
- COMPUTE y=UNIFORM(1).
- COMPUTE z=x**2+y**2.
- END CASE.
- END LOOP.
- END FILE.
- END INPUT PROGRAM.
- EXECUTE.
- DESCRIPTIVES VARIABLES=z
- /SORT SUM
- /STATISTICS=mean.
描述统计量
N 均值
z 2000 .67
有效的 N (列表状态) 2000
第三个例子,更为复杂点;求在正负无穷上即整个实数范围内,求y=1/(1+x^2)的积分。这是个广义积分。是收敛的,有结果的。通过变量代换,转换为在[0,1]上的积分,即y=2/(t^2+(1-t)^2)。转换方法x=(1-t)/t。最后结果是pi。
- INPUT PROGRAM.
- LOOP #I=1 to 200000.
- COMPUTE t=UNIFORM(1).
- COMPUTE y=2/(t**2+(1-t)**2).
- END CASE.
- END LOOP.
- END FILE.
- END INPUT PROGRAM.
- EXECUTE.
- DESCRIPTIVES VARIABLES=y
- /SORT SUM
- /STATISTICS=mean.
描述统计量
N 均值
y 200000 3.14
有效的 N (列表状态) 200000
第四个例子,如何求体积。一个不规则的曲线构成的体积。如Y=Sinx+x。在[0,1]上求体积。
代码如下,
- INPUT PROGRAM.
- LOOP #I=1 to 200000.
- COMPUTE X=UNIFORM(1).
- COMPUTE Y=SIN(X)+X.
- COMPUTE Z=3.1415*Y**2.
- END CASE.
- END LOOP.
- END FILE.
- END INPUT PROGRAM.
- EXECUTE.
- DESCRIPTIVES VARIABLES=Z
- /STATISTICS=MEAN.



雷达卡



京公网安备 11010802022788号







