《数据分析》实验报告六.docxVIP

课程名称 课程编号 实验地点 实验时间 校内指导教师 实验 6 数值模拟 评阅人签字 成绩 1、试用蒙特卡罗方法计算π 。构造一个定积分4?1 1? x dx ? ? ，其中积分S ? 1? x dx 2 2 0 （即 1/4 单位圆的面积）用随机模拟方法求得。在正方形区域内生成一些均匀分布的随机数，随着 由这个几何图形我们可得到计算π 的 Monte Carlo 方法：先生成 N 个正方形区域内均匀分布的 二维随机数，然后计算这 N 个随机数中有多少个落在半圆区域内，设其为 M 个，最后 M/N 就等于 半圆面积除以正方形面积——此比例是π /4。 （1）给定曲线 y =2 –x2 和曲线 y3 = x2，。曲线围成平面有限区域，用蒙特卡罗方法计算区域面 积； ? 2 2 所围体积，其中? ?{(x, y, z) | ?1? x ?1,?1? y ?1,0 ? z ? 2} 。 z ? 2 2 二、实验原理 参数说明 (0,1)上均匀 分布 Y=rand(m,n) 生成参数为 p 的 m 行 n 列的 m×n 个几 何分布随机数 生成区间(a,b)上的 m 行 n 列的 m×n 个均 匀分布随机数 生成参数为 μ,σ 的 m 行 n 列的 m×n 个正 态分布随机数 生成参数为 μ,σ 的 m×n 个对数正态分布 随机数 生成参数为 A,B, 的 m×n 个贝塔分布随 机数 3.1417 k=0; for j=1:100000 a=rand(1,2); if a(1)^2+a(2)^2=1 k=k+1; end end PI=4*k/j PI = 3.1429 2、计算面积、体积 P=rand(10000,2); x=2*P(:,1)-1; y=2*P(:,2); II=find(y=2-x.^2y.^3=x.^2); M=length(II); S=4*M/10000 S = 2.1804 plot(x(II),y(II),g.) 体积 建立函数 function [q,error] = MonteC(L) if nargin==0,L=7;end N=10000; x=2*P(:,1)-1; y=2*P(:,2); z=2*P(:,3); error=q-pi; 说明：此部分的内容和格式可根据实验课程的具体需要、要求自行设计和确定相关栏目。