用蒙特卡洛方法估计积分方法及matlab编程实现.doc

用蒙特卡洛方法估计积分 方法及matlab编程实现 实验内容： 1用蒙特卡洛方法估计积分 ，和的值，并将估计值与真值进行比较。 2用蒙特卡洛方法估计积分 和的值，并对误差进行估计。 要求： （1）针对要估计的积分选择适当的概率分布设计蒙特卡洛方法； （2）利用计算机产生所选分布的随机数以估计积分值； （3）进行重复试验，通过计算样本均值以评价估计的无偏性；通过计算均方误差（针 对第1类题）或样本方差（针对第2类题）以评价估计结果的精度。 目的： （1）能通过 MATLAB 或其他数学软件了解随机变量的概率密度、分布函数及其期望、方差、协方差等2） 熟练使用 MATLAB 对样本进行基本统计，从而获取数据的基本信息3） 能用 MATLAB 熟练进行样本的一元回归分析。 ，积分真值为1 m=f11 输入一元被积函数如x.*sin(x):x.*sin(x) g1 = x.*sin(x) 输入积分下界a:0 输入积分上界b:pi/2 积分真值:1 输入样本容量 10^V1--10^V2: V1:1 V2:5 n = 10 n = 100 n = 1000 n = 10000 n = 100000 Real = 1 EY1 = 1.2635 1.0088 1.0066 1.0109 1.0018 D1 = 0.2635 0.0088 0.0066 0.0109 0.0018 RD1 = 0.2635 0.0088 0.0066 0.0109 0.0018 MSE1 = 0.6439 0.0205 0.0028 0.0006 0.0001 m= 1.2635 1.0088 1.0066 1.0109 1.0018 0.2635 0.0088 0.0066 0.0109 0.0018 0.2635 0.0088 0.0066 0.0109 0.0018 0.6439 0.0205 0.0028 0.0006 0.0001 %估计积分 真值为0.8862 M=f11 输入一元被积函数如x.*sin(x):exp(-x.^2) g1 = exp(-x.^2) 输入积分下界a:0 输入积分上界b:+inf 积分真值:pi^0.5/2%0.8862 输入样本容量 10^V1--10^V2: V1:1 V2:4 n = 10 n = 100 n = 1000 n = 10000 Real = 0.8862 EY1 = 0.9333 0.9077 0.8873 0.8871 D1 = 0.0470 0.0215 0.0010 0.0009 RD1 = 0.0531 0.0243 0.0012 0.0010 MSE1 = 0.1927 0.0112 0.0016 0.0000 M = 0.9333 0.9077 0.8873 0.8871 0.0470 0.0215 0.0010 0.0009 0.0531 0.0243 0.0012 0.0010 0.1927 0.0112 0.0016 0.0000 第一类二重积分程序代码： %%%构造示性函数，求不同区域上积分只需更改示性函数 function I=I2(x,y) if x^2+y^2=1 I=1; else I=0; end %保存为I2.m %第一类二重积分程序主体 %保存为f12.m function outf12=f12() g2=input(输入二元被积函数如exp(x.^2+y.^2):,s)%输入被积函数 g2=inline(g2,x,y); Real=input(积分真值:);%输入积分真值 fprintf(输入样本容量 10^V1*10^V1--10^V2*10^V2:\r) V=zeros(1,2); V(1)=input(V1:);%输入样本