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

实用文案 用蒙特卡洛方法估计积分 方法及 matlab 编程实现 专业班级： 材料 43 学生姓名： 王宏辉 学 号： 2140201060 指导教师：李耀武 完成时间： 2016 年 6 月 8 日 标准文档 实用文案 用蒙特卡洛方法估计积分 方法及 matlab 编程实现 实验内容： 2 2 2 2 1 用蒙特卡洛方法估计积分 ， - x 和 x y 的值， x sin xdx e dx e dxdy 0 0 x2 y2 1 并将估计值与真值进行比较。 1 2 用蒙特卡洛方法估计积分 ex2 dx 和 1 dxdy 的值，并对误 4 4 0 x2 y2 1 1 x y 差进行估计。 要求： （1）针对要估计的积分选择适当的概率分布设计蒙特卡洛方 法； （2）利用计算机产生所选分布的随机数以估计积分值； （3）进行重复试验，通过计算样本均值以评价估计的无偏性； 通过计算均方误差（针 对第 1 类题）或样本方差 （针对第2 类题） 以评价估计结果的精度。 目的： （1）能通过 MATLAB 或其他数学软件了解随机变量的概率密 度、分布函数 及其期望、方差、协方差等； （2） 熟练使用 MATLAB对样本进行基本统计，从而获取数据的 基本信息； 标准文档 实用文案 （3） 能用 MATLAB熟练进行样本的一元回归分析。 实验原理： 蒙特卡洛方法估计积分值， 总的思想是将积分改写为某个随机变 量的数学期望，借助相应的随机数，利用样本均值估计数学期望，从 而估计相应的积分值。 具体操作如下： b b b 一般地，积分 S g (x)dx 改写成 S g (x ) f (x )dx h(x) f (x )dx 的形 a a f(x) a 式， （其中为f(x ) 一随机变量 X 的概率密度函数，且 f(x ) 的支持域 g (x ) ｛x | f (x ) 0｝ （a, b） ) ，h(x ) ）; 令 Y=h(X) ，则积分 S=E （Y）；利用 f(x) matlab 软 件 ， 编 程 产 生 随 机 变 量 X 的 随 机 数 ， 在 由