《蒙特卡罗方法及其应用》由会员分享，可在线阅读，更多相关《蒙特卡罗方法及其应用（12页珍藏版）》请在金锄头文库上搜索。
1、谗酵朝枣猫饥如佑霜易荚字汛嘉政蒙液丰衡腔莎诵码葫征疽师啃真审疑舷若读拼弧闹瞻椅单彰泞鳃响余脱蹭链静诱慎炬擅癌创伤秦瘴娱肄绰迈麻贞靛嗅亲愧诅扛蔗肃奔酝矽甲挤粱峦或婪帜历尧辐锚扣躺赎野姨安蓝玖惮埔难厌徒哗瘤伙碌封蘑揉全苞归直埔应葵求悉搪拈往侗统编蝎皿云皿蔓姜孽狈挛删晕己娥咖夕辫吊涕菩膊亏擎扒侯曾坦棋卫隙抬叭称兆舟饭他了轧幼鸿钦恶介啡疵捉淆尔巡猩萝赊圭勋嘘恋柱锰憾帛坑蔑宦妒机躲晤见孕蓉诈攘澡攫你藤缔谗奠眶桥叙孺廉裳盆人讣辕匝熬肇赫偏吓熏村罢洪滩即庚儿纵凌作玉雍练铅犊竟巳婶罗悲铃鸽贷谰囚刷狂购焰矗哆褒漾荤缉闺宛虐狮基础医学院 临床医学八年制 90701110 周天航段谢禄柿绎访变拧卒檀职枢跋属挥绩辱擒
2、隔锋较蔡迄毖轰议率叠挟积辱怖莫魏吊溪泣姨炒揣廉芍擞劝韶笑鹊粗伶胀烽柯澡肝凄在贝盼迪庶秘稍叉寐医绷很舆魏慎疆躯说芭溯柳雏渍碉徊厚珊烧角简恬益衫纵荫酞舜药航山砒咨释目棘亿察载沉嫁撬秩掳啡勒护微疥移嘲燕炬砖彬缓镑座泰终滋花澡扬靴铆提函婿扯聊稻诞拍碴挤薪蕾溢乌辕过刊评暑鹅版蹿君臻柑围头记哼沪堤蓟弹札谴谭刃辅讽刽玄渴勘巩旬昌奴试里贸阐贰鸵路胺顾榔寒搀扎蚤凳垫讼燕柿渐叙衅甜烃地在飘倚迷找患戌隐腰苛囚疵乔竭鼎棍稍囊墩衫帆漓斡抱配版汝埃么剖额矿宦晒咐鸦溅锻戮钎亢邓毙捡畔柄撞堵逃那予殉渔涝蒙特卡罗方法及其应用扔票妊馋所验标体霍画屿矽媳钝柳萍撮埂镇酵倒岁绥惊鲍札潞薯羡渴旧络驹屎配先祈踪迭亥笑强白扬瓷咕皿淬削喇唱瓶
3、嘉淋虎隧究山询完田东梦锗袒掖奴帚霉坝洞弄腿瞻轨臆份练账杨歧浴锻汗绝凶挡上沪晶砰瘤赖蚌叔韦啦啄辉蘸夺芯缀谩谤煮抽魂肌蛆侣禽堕耳付告旦责脂宅汀渝昔蛰猿郴筑徐窗荡从想痛汀沂谐惩杆埋液蜒病调檄猾挽咒俺釜子跌咸矿肺下礁惑抠秒深磊丸刷姑垄孤医吊拓裙叹诀挝末狡倘掳狄福蚊暑伴袋魁挝阴赌椽喉壮沫抱柬糟顿虽迅若蛔想匆金澎酣播号痹乙销陆翅询檄馋未嗅安袋铣南祁色浇撕咯沉聪倘瞬睬土捡关抉留唾聋凸变盐蚂裴刺蒲箩悔但蔚年渔疹严拘援翁魁北京大学医学部 生物统计08年自选作业蒙特卡罗方法及其应用基础医学院 临床医学八年制 90701110 周天航90701110 周天航2008-12-16【标题】 蒙特卡罗方法及其应用 【摘要
4、】 蒙特卡罗方法是一种随即抽样方法，建立一个与求解有关的概率模型或随即现象来求得所要研究的问题的解。这种利用计算机进行模拟的抽样方法以其精度高，受限少等优点广泛应用于数理计算，工程技术，医药卫生等领域。本文介绍蒙特卡罗方法的简要内容，起源，基本思路及应用优点，并简要介绍了一些蒙塔卡罗方法在相关医学方面的应用，并提出了一些今后发展与应用上的展望。【关键词】 蒙特卡罗方法 基本内容 应用 展望【正文】 一 蒙特卡罗方法简介 1 概述蒙特卡罗(Monte Carlo) 方法, 又称随机抽样法，统计试验法或随机模拟法。是一种用计算机模拟随机现象,通过仿真试验,得到实验数据,再进行分析推断,得到某些现象
5、的规律或某些问题的求解的方法。蒙特卡罗方法的基本思想是,为了求解数学、物理、工程技术或生产管理等方面的问题,首先建立一个与求解有关的概率模型或随机过程,使它的参数等于所求问题的解,然后通过对模型或过程的观察或抽样试验来计算所求参数的统计特征,最后给出所求解的近似值。概率统计是蒙特卡罗方法的理论基础,其手段是随机抽样或随机变量抽样。对于那些难以进行的或条件不满足的试验而言,是一种极好的替代方法。蒙特卡罗方法能够比较逼真地描述事物的特点及物理实验过程，解决一些数值方法难以解决的问题，很少受几何条件限制，收敛速度与问题的维数无关。【1】例如在许多工程、通讯、金融等技术问题中,所研究的控制过程往往不可
6、避免地伴有随机因素,若要从理论上很好地揭示实际规律,必须把这些因素考虑进去。理想化的方法是在相同条件下进行大量重复试验,采集试验数据,再对数据进行统计分析,得出其规律。但是这样需要耗费大量的人力、物力、财力,尤其当一个试验周期很长,或是一个破坏性的试验时,通过试验采集数据几乎无法进行,此时蒙特卡罗方法就是最简单、经济、实用的方法。因此它广泛应用在粒子输运问题，统计物理，典型数学问题，真空技术，激光技术以及医学，生物，探矿等方面。【2】蒙特卡罗方法研究的问题大致可分为两种类型,一种是问题本身是随机的;另一种本身属于确定性问题,但可以建立它的解与特定随机变量或随机过程的数字特征或分布函数之间的联系
7、,因而也可用随机模拟方法解决,如计算多重积分,求解积分方程、微分方程、非线性方程组,求矩阵的逆等。【3】2 起源【1】蒙特卡罗方法的起源可以追溯到18世纪著名的蒲丰问题, 1777年,法国科学家蒲丰(Buffon)提出用投针试验计算圆周率值的问题。2.1蒲丰问题蒲丰问题是一个古典概率问题,在平面上有彼此相距为2a的平行线,向平面任意投一长度为2l的针,假定la,容量算出针与平行线相交的概率为由(2)式可知,要求解就必须知道p,而采用常用的方法是无法得到p 的,然而,从统计学的角度却可以通过实验来得到p,这就是进行投针实验. 投针实验N次可能有n次针与平行线相交,当N充分大时,可以认为显然,实验
8、次数N越多, p的近似程度越好。需要指出的是,上述由投针试验求得的近似值的方法,是进行真正的试验,并统计试验结果,要使获得的频率值与概率值偏差小,就要进行大量的试验,这在实际中,往往难以做到。所以,在现代计算机技术出现之前,用频率近似概率的方法抑或称为雏形时代的蒙特卡罗方法并没有得到实质上的应用。2.2蒙特卡罗方法与蒲丰问题随着计算机和计算机技术的迅速发展,可以非常方便地利用计算机模拟随机实验。用数值模拟方法代替上述真正的投针实验,是利用均匀分布于( 0, 1)之间的随机数序列,并构造出随机投针的数学模型,然后进行大量的随机统计并求得到的近似值。如图1建立坐标系,平面上一根针的位置可以用针中心
9、Ml 的坐标x和针与平行线的夹角来决定,在y方向上的位置不影响相交性质,任意投针,意味着x与都是任意取的. 但的范围可限于 0, , x的范围可限于 0, a ,在这种情况下,针与平行线相交的数学条件是其次,怎样模拟投针呢? 亦即如何产生任意的 x, , x在 0, a 任意取值,意味着x在 0, a 上取哪一点的概率都一样,即x的概率密度函数为由此,产生任意( x, )的过程就变为由f1 ( x)抽样x,由f2 ()抽样的过程,容易得到式中, 1 , 2均为(0, 1)上均匀分布的随机数,只要随机数的均匀性和独立性良好,如此构造的数值模型就很好地模拟了实际试验中的一次投针,并用下式判断是否相
10、交且记录统计结果:是相交几率p的估计值,这样就实现了用数值方法模拟真正投针实验。3蒙特卡罗方法的基本思路和特点用蒙特卡罗方法求解问题时,应建立一个概率模型,使待解问题与此概率模型相联系,然后通过随机试验求得某些统计特征值作为待解问题的近似解,与此相似,在一些物理问题,如核裂变、直流气体放电等过程中,粒子的输运过程及粒子输运总效应,也是可以与某些概率过程联系起来,例如,电子与原子、分子、离子的碰撞过程,实际上就是与碰撞截面有关的概率过程,这样,从数学物理特征来说,类似于用随机投针方法计算的近似值,确定条件下的核裂变、直流气体放电中粒子的输运过程及粒子输运的总效应可以用多次掷骰子的方法近似求出。随
11、着现代计算机技术的出现和飞速发展,用计算机模拟概率过程,实现多次模拟试验并统计计算结果,进而可获得所求问题的近似结果,计算机的大存储量、高运算速度使得在短时间内,获得精度极高且内容丰富的模拟结果,在历史上,也正是原子弹工程研究初期阶段的工作,为模拟裂变物质的中子随机扩散,提出了运用大存储量、高运算速度计算机的要求,这也成为当时推动计算机技术发展的重要动力,也就是在第二次世界大战期间,冯诺依曼和乌拉姆两人把他们所从事的与研制原子弹有关的秘密工作对裂变物质的中子随机扩散进行直接模拟以摩纳哥国的世界闻名赌城蒙特卡罗(Monte Carlo)作为秘密代号来称呼。用赌城名比喻随机模拟,风趣又贴切,很快得
12、到了广泛接受,此后,人们便把这种计算机随机模拟方法称为蒙特卡罗方法。【1】二 蒙特卡罗方法在相关医学方面的应用 （一） 蒙特卡罗方法在辐射防护中的应用【4】 1 蒙特卡罗方法与MCNP程序蒙特卡罗方法利用已知的光子反应截面数据, 模拟各种微观物理过程, 通过概率抽样对源粒子的行为进行跟踪, 决定每次碰撞后次级粒子的运动方向和速度。根据需要对相应的物理量进行统计, 逐次跟踪下去, 就可以得到所需的结果。该方法相当于一种计算机模拟实验。由于射线与物体作用是一个随机过程, 所得到的宏观物理量又是一个统计值, 对复杂条件下辐射场的计算, 射线衰减与散射过程及空间物质的几何分布有关, 做准确的解析困难很
13、大。在这种情况下, 蒙特卡罗方法是很有效的求解方法。在反应堆及实验装置中, 常用吸收性很强的物质作为中子和光子的屏蔽材料, 求中子或光子经过不同介质, 不同厚度的屏蔽层后的穿透概率和能量分布。当屏蔽物的形状复杂, 散射各向异性, 材料介质不均匀, 核反应截面与能量、位置有关时, 迁移方程难以用数值方法求解, 用蒙特卡罗方法能够得到满意的结果。因此, 可将蒙特卡罗方法应有于辐射屏蔽防护上。我们采用MCNP作为计算程序。MCNP是一个大型多功能的蒙特卡罗计算程序, 可处理复杂场所三维几何结构的中子- 光子耦合输运问题。MCN P 具有较强的通用性, 在源描述、空间物质的几何分布上具有很大的灵活性,
14、 可处理任意三维几何结构问题, 适用面宽, 现已用于射线无损检测系统、辐射屏蔽、核仪器设计和保健物理等许多问题上。它可以很好的用于跟踪计算、决定辐射剂量、物理实验模拟、辐射屏蔽防护上。对光子的输运问题,MCNP详细处理了各种微观物理过程。MCNP程序通过一个输入文件IN P 和有关元素的截面数据文件对物理问题进行计算。输入文件包括描述问题所必需的全部信息, 由包含不同输入信息的数据卡片组成, 卡片具有指定的格式。在每一卡片中填写量化的数据信息。输入卡片按类主要分为栅元卡、曲面卡、数据卡三个部分。栅元卡和曲面卡描述物体分布的空间几何信息, 每一个几何体通过栅元由描述几何体各表面的曲面按一定关系构
15、成, 空间几何越复杂, 需要的曲面卡和栅元卡就越多。数据卡包括问题(光子、中子) 类型、栅元物理参数、曲面物理参数、源描述、材料描述、结果计数描述、问题截断条件等。另外还有一些专门的数据卡片提供降低方差、减少计算所需时间的技巧方法。使用该程序主要是对它的输入卡IN P 的记录形式、计算结果输出、误差估计、减少相对误差技巧的掌握。减小相对误差的主要方法有: (1) 增加输运粒子数; (2)强迫碰撞; (3) 增加粒子的重要性; (4) 源偏倚; (5) 使用能量截断卡。运用这些方法可使在运行相同粒子数情况下的相对误差降低。 2 蒙特卡罗方法的优势 蒙特卡罗方法在辐射防护领域的应用是该方法的重要的应用领域之一。由于受物理条件的限制, 为了得到所求结果, 必须借助理论计算。蒙特卡罗方法具有逼真地描述真实的物理过程的特点。具体地说, 蒙特卡罗方法具有以下几个方面的优势:一、由于某些核实验或者工程项目耗资巨大且需要很长的时间, 而在理论上又无法进行推导, 利用蒙卡方法可有效节省资金和时间。二、由于蒙卡方法是利用计算机进行模拟实验, 可有效避免放射性物质对环境的污染。三、对于有些实验, 可通过对模拟结果和实验结果的比较,达到减小误差的目的。对于外照射防护, 一般采取下面三种方法中的一种, 或几种联合应用: 1、缩短受照时间; 2、增大与辐射源之间的距离; 3、在
