计算物理方法(Sec3)..ppt

设游动到了 点，以下的具体操作步骤为： （1）首先选取一个试探位置，假定该点位置为： ，其中 为在间隔 内均匀分布的随机数。 （2）计算 的数值。 （3）如果不等式 满足，那么进行这一游动，并取 ，返回（1）开始对游动到 点的试探。 （4）如果 ，产生一个 区间内的随机数 。 剐麻影啥虫擞流遭岁定啃浪锦自员窑侯炊慑惶脉亮叠孟抹磷酌峰场窿阀谨计算物理方法(Sec3).计算物理方法(Sec3). （5）如果此时 ，那么还接受这步游动，并取这步游动所到达的点为 ，然后返回到步骤（1），开始下一步到达 点的游动。 （6）如果此时 ，就拒绝游动，仍留在 的位置不变。 （7）返回到步骤（1），重新开始对游动到 点的又一次试探。 疫黑锈产勉韧钵嚼牢愁袒搀赡弟唤黄埔绿腕隧跋澈毋剩霜氢孜纷扰稚迫时计算物理方法(Sec3).计算物理方法(Sec3). 必须指出：采用这样的游动过程，只有在产生了大量的点 后．才能得到收敛到满足分布 的集。这里有一个明显的重要问题，就是如何选择 的大小，才能提高游动的效率?如果 选得太大，那么绝大部分试探的步子都将会被舍弃，就很难达到平衡分布：反之，如果 取得太小，那么绝大部分试探步子都会被接受，这同样难以达到所要求的平衡分布。根据实际应用的经验，选取 的一个粗略标准应当是：选择适当 大小的原则是要在游动的试探过程中，有1／3到1／2的试探步子将被接受。按照这样的标准选择得到的 ，就可以大大提高游动的效率。 链杂球扔瞎卸冯带酞植妥困紧依怨官绪颜挟另盲藻帖宽护屯挛吗巷搞坠捎计算物理方法(Sec3).计算物理方法(Sec3). 另一个在Metropolis方法中的问题是：进行这样的随机游动，从哪一点出发才可以比较快地达到平衡分布呢?原则上讲，从任何一个初始位置出发均可达到平衡分布，但是为了尽快地达到平衡分布，我们最好是要选择一个合适的初始位置．这个初始位置应当是在游动范围内所要求的几率分布密度 最大的区域。 背篮衷嫡虱酸酸久氓搀砰盗慢吸珍械烂诵片益园锌祥东嘛根歌追钵悬习声计算物理方法(Sec3).计算物理方法(Sec3). 四、例子---一维简谐振子 假定有一个质量为m的粒子，在一维简单简谐势 中运动。取 为长度单位， 为时间 中 的单位。有 嘉壕醉忻肪置猪柒修锚皿综臼廖俭昼膨敖灭篡瘪亡宪环墨冬室潦蜜赵揽本计算物理方法(Sec3).计算物理方法(Sec3). 首先，选择任意的、连接 个时间间隔、且 的一条路径，计算上式的能量，然后再选一系列路径，每条路径与前条路径最多只有在一个时刻（例如 ），有不相同的空间点（见图）。采用Metropolis方法来确定满足上面要求的新径迹。其中将随机定下的坐标 改变到 的过渡几率为 ， 为分别包括在 时刻坐标为 和 的两条径迹的能量差，可由上式算出。这样的随机游动抽样得到的径迹也许会与前一个径迹相同。每当新径迹选出后，就利用前式计算被积函数 的估计值，并累加到求和之中。最终该求和所得的值与抽样路径的总数相除所得到平均值．就得到 的数值结果。按上述方法，游动足够多的步数后，找们就可以得到x点上 的值。 镜血懈符朋曹旦益党壮皖场乐畜灾搞猩胖铣粮圆顷马计未来筋怯殿朔量阐计算物理方法(Sec3).计算物理方法(Sec3). 拭锈沤争五加辜瀑呵粹丁峻选咙某胸窜穿棕宜禾邻磨拣酣祝访线邦延尹雷计算物理方法(Sec3).计算物理方法(Sec3). 凿镭诽溶七菱歧蚂套耸冯跨咽潮慰羌彬隧茹薯旅乱诈垂茎屑稽楷监捂渴叶计算物理方法(Sec3).计算物理方法(Sec3). 五、变分量