求助:关于泊松分布的蒙特卡罗模拟 excel

(window.slotbydup=window.slotbydup || []).push({
id: '2014386',
container: s,
size: '234,60',
display: 'inlay-fix'
&&|&&0次下载&&|&&总9页&&|
您的计算机尚未安装Flash,点击安装&
阅读已结束,如需下载到电脑,请使用积分()
下载:5积分
5人评价44页
1人评价15页
1人评价1页
0人评价13页
0人评价11页
所需积分:(友情提示:大部分文档均可免费预览!下载之前请务必先预览阅读,以免误下载造成积分浪费!)
(多个标签用逗号分隔)
文不对题,内容与标题介绍不符
广告内容或内容过于简单
文档乱码或无法正常显示
若此文档涉嫌侵害了您的权利,请参照说明。
评价文档:
下载:5积分什么是泊松分布?泊松分布是一种统计与概率学里常见到的离散概率分布,由法国数学家西莫恩·德尼·泊松(Siméon-Denis Poisson)在1838年时发表。泊松分布适合于描述单位时间(或空间)内随机事件发生的次数。如某一服务设施在一定时间内到达的人数,电话交换机接到呼叫的次数,汽车站台的候客人数,机器出现的故障数,自然灾害发生的次数,一块产品上的缺陷数,显微镜下单位分区内的细菌分布数等等。如果您对搜索结果不满意,可以在百度中搜索“
&nbsp&nbsp|
&nbsp&nbsp|
&nbsp&nbsp|
&nbsp&nbsp|
&nbsp&nbsp|
&nbsp&nbsp|
如有投资本站或合作意向,请联系(010-);
邮箱:service@pinggu.org
投诉或不良信息处理:(010-)
京ICP证090565号
论坛法律顾问:王进律师404 Not Found
404 Not Found
openresty/1.11.2.350中国农村水利水电#2002年第8期用MATLAB实现蒙特卡罗法计算结构可靠度冯晓波 杨 桦(武汉大学水利水电学院430072)摘要 针对应用蒙特卡罗直接抽样法解决结构可靠度所遇到的困难,提出利用MATLAB的强大数值计算功能解决此类问题。(]利用MATLAB进行蒙特卡罗抽样模拟,在一定程度上减少了对连续型分布采用直接抽样时的困难,大大提高了计算效率。关键词 蒙特卡罗法 MATLAB 可靠度 随机数系,例如常用的矩阵代数运算、数组运算等,使用户可以以多种形式快速地操作数据集。与Basic、Fortran、C等编程语言相比,MATLAB具有编程简单直观、用户界面友善、开放性强等特点,大大提高了编程效率。将MATLAB用于蒙特卡罗法的一个显著优点是它拥有功能强大的随机数发生器指令。1 蒙特卡罗法在结构可靠性分析中运用蒙特卡罗法[1]方法,首先考虑各基本变量相互独立的情况。设基本变量x1,x2,,,xn分别有分布函数Fx(x1),Fx(x2),,,Fx(xn),因为Fx(xi)为[0,1]区间12ni上的一个数,可以将其与由蒙特卡罗法产生的随机数rj对应。1这样,便可得到xi=F-对于每一个rj值,x(rj),i=1,2,,,n。i2.2 随机数的产生通常在用蒙特卡罗直接抽样法时,必须解x=F-1(r)以求得服从相应分布类型的随机变量x,再代入功能函数求解。而MATLAB(6.0版)提供了23种随机变量分布类型的随机数发生器,如正态分布、对数正态分布、泊松分布、威布尔分布等(这些分布基本包括了工程实际中出现的变量分布情况),可直接产生变量x以代入功能函数,省去了可能会带来很大麻烦的求分布函数反函数这一步,极大地提高了效率。常用指令举例如下:r=rand(m,n)产生m行n列的(0~1)间均布随机数组r。r=normrnd(Mu,Sigma,m,n)产生服从N(L,R2)分布的m行n列随机变量数组r。r=lgnrnd(Mu,Sigma,m,n)产生lnR服从N(L,R2)分布的m行n列随机变量数组r。 MATLAB目前版本尚无极值型分布的随机数发生器,故对极值型分布,仍须求解相应反函数。可以得到一组对应的基本变量x1,x2,,,xn。将这组值代入功能函数g(x1,x2,,,xn),便得到一个值,该值若小于等于0,则在程序中记录一次功能函数的实现,大于0则不记入。再对另一随机数重复进行这些计算,直到完成预定的循环次数。假定所进行的循环次数为K次,g(x1,x2,,,xn)[0的次数为m次,则可得失效概率为Pf=m/K(要求K足够大)。若基本变量相关,可利用条件概率密度,把多维问题化为一维问题来解决,具体见文献[2]。因此从理论上说,该方法的应用几乎没有什么限制。但是在实际问题中,连续型分布是很复杂的。有的只能给分布函数的解析表达式,但给不出其反函数的解析表达式,如著名的B分布;有的则连分布函数的解析表达式都给不出,如正态分布。所以通常情况下对连续型分布采用直接抽样是有一定困难的。本文利用MATLAB的强大数值计算功能,实现了在MATLAB中采用蒙特卡罗直接抽样计算结构可靠度,较好地解决了上述问题。2.3 编程求解由于MATLAB能直接产生服从各相应概率分布函数的随机变量数组,从而使编程过程大为简化。设功能函数为g(x1,x2,,,xn),xi(i=1,2,,,n)为服从各相应概率分布函数的随机变量。首先用相应随机数发生器指令产生mc@nc的随机变量数组xi(k,l),(i=1,2,,,n;k=1,2,,,l=1,2,,,nc)。然后,将各数组中各元素一一对应代入功能函数。由于MATLAB提供了数组运算指令,避免了使用循环语句,从而执行更快。最后得到功能函数结果数组,统计0j,:2 在MATLAB中实现蒙特卡罗法在结构可靠度中的计算2.1 MATLAB简介MATLAB是一种功能极其强大的科学和工程计算数学软件系统。它汇集了大量数学、统计、科学和工程所需的函数,除了具有类似于其他计算机编程语言的编程特性外,对计算数学领域的特定数学问题,MATLAB都给出了该问题的各种高效算法。范围的matlab 泊松分布 用MATLAB实现蒙特卡罗法计算结构可靠度用MATLAB实现蒙特卡罗法计算结构可靠度 冯晓波 杨 桦jmc@nc51值分布的随机变量。[]在MATLAB中只要2条指令即可完成上述过程,远比Basic、Fortran、C等语言简单高效。最后此例结果为Pf=4.27@10-5。Pf=如果其中某一随机变量为极值型分布,则先用均匀分布随机数发生器指令r=rand(mc,nc),产生一个mc行nc列的均布随机变量数组r,然后用式x=F-1(r)解出服从极值分布的随机变量数组,再与用前述方法产生的其余变量一起代入功能函数求解即可。程序框图如图1:4 结 语(1)蒙特卡罗法在结构可靠度的数值模拟中具有收敛速度与基本变量的维数无关,极限状态函数的复杂程度与模拟过程无关且无需将状态函数线性化和随机变量/当量正态化0的特点,具有直接解决问题的能力。它的普遍适用性使其成为当今可靠度理论研究的热点,必将在结构可靠性分析中发挥越来越大的作用。(2)在蒙特卡罗法计算中,当结构失效概率很小时,抽样模拟次数很大,会大大增加计算量,尤其当功能函数的计算需要进行结构整体有限元分析时更是如此。但随着各种改进方法[4~5]和高性能计算机的出现,这些问题将得到较大的解决。(3)MATLAB的强大功能为结构可靠度计算提供了便利,研究人员可迅速编出科学高效的计算程序,大大提高了效率。在西方,MATLAB的应用已遍及现代科学界和工程界,相比之下,我国工程界人员对其了解和应用就少很多,MATLAB的强大功能等待去发掘。参考文献1 徐钟济.蒙特卡罗法[M].上海:上海科学技术出版社,19852 何水清,王善.结构可靠性分析与设计[M].北京:国防工业出版社,19933 张志涌.精通MATLAB5[M]1第3版1北京:北京航空航天大学出版社,20004 李亚东.重要抽样模拟及其在结构可靠度计算中的应用[J].西南交通大学学报,1992(4)5 程耿东,蔡文学.结构可靠度计算的近似重要抽样模拟及其应用[J].工程力学,)t图1 用MATLAB实现蒙特卡罗法程序框图3 算 例(1)算例1。失效函数g=2.5-N[C/(1+E)]Hlg[(P+$P)/P]有6个相互独立的正态变量,而且包含对数运算,具体数值见表1。
免责声明:文章由投稿,不代表本站的观点和立场!如有问题,请与本站联系。
您可能还喜欢
 虽然很多人喝水的时候都是用自己常用的杯子,但是总有一些公开场合或者是外出的时候,我们需要购买矿泉水喝完之后就产生了很多塑料瓶子,为了减少这些难降解的塑料瓶子,有团队发明了Ooho,Ooho是什么?有了Ooho之后人就能把水连瓶吃进嘴...
 《进击的巨人》第二季的播出不出意料引来了很多人的关注,有不少喜欢这部动漫的人对于进击的巨人2很是期待,但是却不知道进击的巨人第二季在哪看?下面十万个为什么网小编将为你解答进击的巨人第二季在哪看(b站)?进击的巨人2有哪些疑点? ...
 北京时间12月15日凌晨3:00,美联储公布利率决议,将联邦基金利率提高25个基点,美联储利率决定上限为0.75%,下限为0.5%。这是时隔几乎整整一年之后,美联储又一次加息。按照美联储的预计,明年还可能进行三次加息。因为美联储加息的问题...
 充满戏剧效果的2016年美国大选终于落下帷幕,共和党候选人特朗普以黑马之姿杀入美国大...
 在印尼巴厘岛召开的国际刑警组织第85届全体大会上,中国公安部副部长孟宏伟高票当选新一任主席,任期4年!下面为什么网带你看看国际刑警组织有哪些权利 孟宏伟警衔是什么?副总警监 正部级。 孟宏伟当选国际刑警组织主席 比起吃瓜群众看着...
大家都在看}

我要回帖

更多关于 蒙特卡罗模拟方法 的文章

更多推荐

版权声明:文章内容来源于网络,版权归原作者所有,如有侵权请点击这里与我们联系,我们将及时删除。

点击添加站长微信