蒙特卡罗法在理论力学中的应用.pdfVIP

第 3O卷 第 l期 高 师 理 科 学 刊 Vo1．30 No．1 2010年 1月 JournalofScienceofTeachers CollegeandUniversit? Jan． 20l0 文章编号：1007—9831(2010)0l一0064—04 蒙特卡罗法在理论力学中的应用 尹新国，公丕锋，朱孟正，张峰 (淮北煤炭师范学院 物理与电子信息学院，安徽 淮北 235000) 摘要：将蒙特卡罗方法应用到理论力学教学中，并借助于MATLAB软件编程解决简单的物理问题， 提高了学生学习蒙特卡罗方法的兴趣，在引导学生学习和理解蒙特卡罗方法9a-得了一定的成效． 关键词：蒙特卡罗方法；理论力学；MATLAB 中图分类号：0213．2 文献标识码 ：A 随着现代计算机技术的发展，蒙特卡罗方法已经在 自然科学研究中发挥了重要的作用，并正在 日益广 泛地应用于物理工程的各个方面．鉴于蒙特卡罗方法的重要性，将其应用到理论力学教学中，并借助于 MATLAB软件编程解决简单的物理问题，如求出物体的质心和转动惯量，提高了学生学习蒙特卡罗方法的 兴趣，在初步引导学生学习和掌握蒙特卡罗方法中取得了一定的成效，为进一步选修 《数值计算》课程和 掌握蒙特卡罗方法奠定了基础． 1 样本平均值法计算定积分 使用样本平均值法计算定积分“e )，7PI))~IN(a，6)上的一个概率密度函数g)，有 )ibPX()g( lj ) 对于均匀分布的密度函数g()=1／b(—a)，积分 (1)的估计值是 I2px()d~≈喜 i=1)b．-a- (2) 其中：Xi是区间a，b)内的，z个按照均匀分布的密度产生的随机数，可以通过MATLAB的mnd0m函数生成的，1 个 (0，1)区间的随机数计算得到，即 x=a+(6一口)rand (3) Xi通过式(3)可以一个一个产生计算得到，也可以一次性算出．如MATLAB的命令randm，，1)可产生 行 n列共mn个在0到1之间均匀分布的随机数的矩阵． 对于三重积分，可将式 (2)推广为 z ii袭 = ) (4) 式 (4)中的积分区域Q是长方体区域，若Q不是长方体区域，可以取包含区域Q的长方体区域 ，有 泐 喜条 = ) ㈩ 其中：JV是由均匀密度函数g=1／【一口)一cXf—P)】在长方体区域 内产生的随机点数； 是落在区域Q 内的点数，求和是对区域Q内的点进行的．本节的理论对于没有学过蒙特卡罗方法的学生来说，理解起来 收稿日期 ：2009—05—25 作者简介：尹新国 (1961一)，男，安徽霍山人，教授，主要从事理论物理研究．E—mail：msyyxg3081@126．2(0111 第 1驯 新 ，等：荣特 罗法在 论 ／J掌rl『n 川 65 一 — — — — — — — — — — — — — — — — — — — — — — — — — — — — — — — — — — — — — — — — — — — — — — — — — — — — — — — — — — — — — — — — — — — — — — — — — — — — — — — — — — — — — — — — — — — — 一 有一定的难度．为此，用物理的思维米诠释式(5)有必要．可以这样理解：n (，’，z)cL~-d3d’z表示的是 具有连续质量体密度分布 的，体积