python 计算物理_计算物理期末报告

| 姓名 | 学号 | 班级 | 选题 | 论述 | 结论 | 总分 |

| ---- | ---- | ---- |

| 余康 | 2014301020117 | 弘毅班 | 随机系统 | | | | |

不同类型随机行走过程的python模拟

余康　2014301020117　14级物理学弘毅班

前言

这篇文章是作者(我)在阅读过《计算物理》教材第七章“随机系统”之后的作品，由于时间与知识水平的限制，所以无法涵盖第七章的全部内容。基于以上的考虑，作者选择了相对较为基础的内容作为自己的课题，即“随机行走的过程的python模拟”。尽管如此，但是作者还是尽力地完善了自己的作品，相当于对书中第七章第一二节的内容做了自己的补充。

摘要

从布朗运动我们开始接触随机行走，布朗运动是由于运动是由于颗粒受液体分子碰撞的不平衡力引起的。随机行走可以理解为布朗运动的理想状态，也可以称为无规则行走，或者是形象的解释为醉汉的行走。不止于布朗运动，随机行走过程在我们生活中的实例与应用是极其广泛的，比如实例有所有的扩散现象，应用于互联网的链接分析，金融的股票分析和高分子的构象(自回避行走)等领域。随机行走过程的模拟可以通过随机数生成器来实现，从而我们可以通过数值方法验证解析方法得出随机系统的概率描述(例如平均数，方差等)并且直观地认识随机系统的规律。对于一维与二维随机行走的模拟较为简单，但是对于三维的随机行走(全方向，不限于格点)的模拟，就涉及如何在三维球面上随机取点(均匀分布)，我采取了三维球面上的Marsaglia方法。

关键词：　随机行走　模拟　随机数　概率描述　Marsaglia方法

Ⅰ 介绍

随机行走是一个数学对象，它描述了由一系列随机步骤组成的路径。例如，跟踪液体或气体中分子的运动路径，搜索猎物的觅食动物的运动路径，超弦的行为，一支股票价格的波动和一个赌徒的财务情况都可以通过随机行走模型近似，尽管在现实中他们可能无法实现真正的随机。这些例子说明了随机行走在许多科学领域的应用， 包括生态学、心理学、计算机科学、物理学、化学和生物学，以及经济学。随机行走可以解释在这些领域的许多过程中所观察到的行为，从而作为随机活动的基本模型。谈到更数学的应用，π的值可以采用基于随机行走的编程方法数值逼近(Monte Carlo method，蒙特卡罗算法)。随机行走是由卡尔·皮尔逊在1905年第一次提出的。

人们对于随机行走有各种各样的不同方面的研究兴趣。这个词本身通常指的是一种特殊的马尔可夫链或马尔可夫过程。但是许多具有特定性质的随时间变化的过程也被称为随机行走。随机行走(无论是否马尔科夫)也能够发生在各种空间：大家通常研究的有图(图论)，整数(数论)，实线、平面、更高维的向量空间，曲面、更高维的黎曼流形，有限群、无限群、李群(群论)等。时间参数也可以被操控。在最简单的情形下，我们研究的是离散时间序列的随机行走。当然，我们也可以定义任意时间序列的随机行走，在这种情形下，所有时刻的位置都有定义。随机行走的实例包括莱维飞行，扩散模型(比如布朗运动)。

随机行走是马尔可夫过程讨论的一个基本主题.。它们的数学研究已十分广泛。随机行走的几个特性，包括分布(dispersal distributions)，首次超越(first-passage)，碰撞次数( hitting times)，会遇率(encounter rates)，循环(recurrence)，瞬态(transience)，已被引入来量化它们的行为。

随机行走所满足的扩散定律：

其中D为扩散常数。

以下是四幅图是维基百科上对二维与三维随机行走的模拟结果，也是我在编程过程中力求达到的目标。

二维：

二维随机行走2500步

二维随机行走25000步

二维随机行走2000000步

三维：

三维随机行走

Ⅱ 正文

一维随机行走的模拟

一维的随机行走，是最简单的随机行走模型，但是它的研究方法以及某些结论却并不局限于一维的情况，可以轻易地平行推广到更高维的情形，所以一维随机行走的讨论是本文的基础。

定步长，等概率

考虑最简单的情形，即一维定步长随机行走，也就是左右等可能随机行走，我们把模型简化为：一个醉汉在一条直线上行走，每次只能走一个单位长度，并且他向每一边走的概率都是完全相等的(1/2)。我采用random.random()函数生成一个从0到1的随机数作为醉汉每一步行走方向的标尺，如果这个随机数大于0.5，则醉汉下一步会向右走一个单位长度；否则，醉汉下一步会向左走一个单位长度。

借助模拟，我们可以观察醉汉行走一百步之后的位移，但是由于每一次醉汉的行走都是随机的，有意义的只是模拟次数足够多之后所表现出来的统计规律，所以我们可以假想有一个醉汉系综，实际上我们通过循环结构即可实现，他们都行走了一百步之后，对他们与原点的位移取平均值，才能验证我们对这个随机系统(醉汉系综)的性质的预测。

接下来我将展示模拟结果，首