1、基 于 MATLAB 的 蒙 特 卡 洛 方 法 对 可 靠 度 的 计 算可靠性工程大作业2目录目录 .2摘要 .3绪论 .4一、编写 MONTE CARLO 模拟程序 5二、关于两个服从正态分布的可靠性验证 .8三、非正态分布的验证 .10四、总结 11参考文献 123摘要对于简单的概率计算，我们可以用离散或者连续的概率分布模型进行求解；但是对于复杂的模型的近似解的求解，蒙特卡洛方法是一种非常方便的方法。蒙特卡洛方法将最复杂的计算部分交给了电机计算机来完成，极大的方便了我们的求解过程。本文主要是用 MATLAB 编写蒙特卡洛的模拟程序，然后分别验证两个正态分布的模型和两个非正态分布的模型。
2、非正态分布的模型中的随机变量序列都是独立同分布的，这样我们可以方便的用列维-林德伯格中心极限定理进行处理。【关键字】：复杂模型、蒙特卡洛、MATLAB、正太分布、独立同分布的非正态模型、列维-林德伯格中心极限定理4绪论计算机技术的发展，促进了蒙特卡洛方法的推广、普及以及完善等。蒙特卡洛方法诞生之初是不被重视的，因为当时的计算机技术没有达到与之匹配的程度。蒙特卡洛模拟也称为随机模拟方法，或随机抽样技术。它是一种以概率论和数理统计为基础，通过对随机变量的统计实验、随机模拟来求解问题近似解的数值方法。它的主要思想是：为了求解数学、物理、化学及工程问题，建立一个概率模型或随机过程，使它的参数等于问解；
3、然后通过对模型或过程的观察或抽样来计算所求参数的统计特征（如均值、概率等） ，作为待解问题的数值解，最后给出所求解的近似值，而解的精度可用估计值的方差来表示。蒙卡洛模拟的步骤是：首先建立简单而又便于实现的概率分布模型，使分布模型的某些特征（如模型的概率分布或数学期望）恰好是所求问题的解；然后根据概率分布模型的特点和计算的需要改进模型，以便减少方差，降低费用，提高计算效率；再对分布模型进行随机模拟，其中包括建立产生伪随机数的方法和建立对所遇到的分布产生随机变量样本的随机抽样方法；最后建立各种统计量的估计，获得所求解的统计估计值及其方差。蒙特卡洛模拟方法可分为直接蒙特卡洛模拟、间接蒙特卡洛模拟和蒙
4、特卡洛积分。（1）直接蒙特卡洛模拟采用随机数来模拟本身具有复杂随机过程的效应。该方法是按照实际问题所遵循的概率统计规律，用计算机进行直接的抽样，然后计算其统计参数。直接蒙卡洛模拟法能充分体现蒙特卡洛方法的特殊性和优越性，因而在物理中得到了广泛的应用，该方法也就是通常所说的“计算机实验” 。（2）间接蒙特卡洛模拟是人为地构造出一个合适的概率模型，依照该模型进行大量的统计实验，使它的某些统计参数恰好是待求问题的解。Buffon 投针实验就是运用间接蒙特卡洛模拟来求解 。（3）蒙特卡洛积分是利用随机数系列计算积分的方法，积分维数越高，效率越高。定积分的计算是蒙特卡洛方法被引入计算数学的开端，这里以定
5、积分的计算说明其处理确定性问题的方法。如计算定积分： dxfks10)(1)(f此时，求定积分亦即求边长为 1 的正方形中一个曲边梯形的面积问题，如图 2 所示。可以随机地向正方形内投点，然后统计落在曲线下的点数 M，当总的投点N充分大时， NkM/就近似等于积分值 s。56一、编写 Monte Carlo 模拟程序1模型的建立 本章节根据抛掷骰子编制 Monte Carlo 模拟程序，验证各点出现的概率均为1/6。2模拟流程图绘制 初始化i=i+1K=？K=1 K=2 K=3 K=4 K=5 K=6K1+1 K2+1 K3+1 K4+1 K5+1 K6+1i0P(j)=P(j)+1;end
6、endR(j)=P(j)/N end 11三、非正态分布的验证对于非正态分布的强度-应力随机变量的可靠度计算，我们再 MATLAB 上用蒙的卡罗方法来验证。验证时我们取样本值 n=100000，分别验证强度服从期望为 10（及 = ）指数分布（x0p(j)=p(j)+1;end endR(j)=p(j)/n; end 12四、总结根据强度-应力干涉模型求解系统的可靠度，对于强度和应力都服从正态分布的干涉模型，查表计算法和蒙的卡罗方法都是正确有效的。但是当强度或者应力不再服从正态分布时，用蒙的卡罗方法可以干涉模型近似结果。13参考文献1朱陆陆. 蒙特卡洛方法及应用D. 华中师范大学,2014.2刘岚岚. 可靠性基础 D.中国质检出版社,2014.3张宇. 概率论与数理统计 9 讲J. 北京理工大学出版社.
