A Brief Introduction to Random and Its Applications In Some Physical Models

姓名:譚善
學(xué)號(hào):2014301020106
班級(jí):14級(jí)弘毅班

Abstract

A stochastic process (or random process) is a probability model used to describe phenomena that evolve over time or space. It plays quite important role in Statistical and Thermal Physics, especially Thermal Physics. It is so fundamental that it was put in the first chapter in most Thermal textbooks. In this paper, I mainly talked random walks in one, two and three dimension situa-tions. Besides, I disscussed another example about random. Here, I used different functions of python to realize our goals.


1. Backgroud


Stochastic processes were first studied rigorously in the late 19th century to aid in understanding financial markets and Brownian motion. The first person to describe the mathematics behind Brownian motion was Thorvald N. Thiele in a paper on the method of least squares published in 1880. This was followed independently by Louis Bachelier in 1900 in his PhD thesis "The theory of speculation", in which he presented a stochastic analysis of the stock and option markets. Albert Einstein (in one of his 1905 papers) and Marian Smoluchowski (1906) brought the solution of the problem to the attention of physicists, and presented it as a way to indirectly confirm the existence of atoms and molecules. Their equations describing Brownian motion were subsequently verified by the experimental work of Jean Baptiste Perrin in 1908. It boosted the development of physics, so, we are really interested in it.


2. Generation of random number

Random numbers are required in different techniques of statistics, such as when a representative sample is taken from the population as a whole, or when animals are assigned to a different test group or during a Monte Carlo simulation and so on. True random numbers are generated using physical phenomena. But, we can obtain Pseudo random number by using computer.

  • Just like following two groups date. We can get random number using "prob=random.random()" or "x=random.randrange(N)":



This is what we obtained, if you run the program again and again, you will find that there is no regularity about the number we get.

3. Random walks in one dimension

Assuming that a walker that is able to take steps of length unity along a line. And the walker can go left and right, we just want to study its position(or deviation from the origin).
In theory, setting each step as s_i, then, we have




Simplifying it, we have


Now, one can assume that the walker move one step during \Delta t=1, and the walker can go right and left with the same probability. Besides, in order to study the walker's mean position and mean value of x^2. We just looped over 10000 times.

  • Following pictures can represent two walkers and three walkers in one dimension, their tracks are random.

Click the Code


From here, we just know that we can predict the walker's track.
  • The following is the figure about the walker's mean displacement.

Click the Code

Obviously, the mean displacement of the walker will fluctuate near 0. In theory, the walker have same probability to go to right and left, so it is not surprising to get this result.
  • Instantly, we can plot the mean value of the square of the displacement.

Click the Code

From the figure, we know the mean value of the square of the displacement is in direct proportion to the time t. In theory, we have
So, here, D=1/2
  • Now, we change the probability about going to left and going to right. Going to right with the probability 1/3 and goint to left with the probability 2/3. We can study some values of this random process. Following two figures represent its track and mean displacement versus time.



  • We can change the probability again, and we study the average of the x squared. One situation is that going to right with the probability 1/3 and going to left with the probability 2/3, the other is that going to right with the probability 1/4 and going to left with the probability 3/4.


Although the probability changed, the average of the x squared and time t into a quadratic relationship still. The fits the theory very well.

4. Random walk in two dimension

Python has provided us all kinds of functions, so we can take advantage of them sufficiently. The turtle can help us simulate random walks of two dimension. Assuming that there is a particle in the orign, in the next time-\Delta t, it can go up and down, go left and right with the same probability, then, we can plot its track using random process theory.

  • When the length of each step is constant, we have the below picture(left). When the length of each step can be arbitrary range from 0 to 1, we can obtain the following figure(right).

Click the Code


Just like the situation of one dimension, its track is unpredicted. These two figures embodied randomness, it is helpful for us to understand the motion of particles that on the surface.

5. Random walk in three dimension

This situation is the most universal situation in our daily life, especially for us who are major in physics. That is because we live in 3-D world. Particles always move in this way. Using python to produce random number, we can simulation particles' motion.

  • First, we assume the particle can move in six directions, that is the particle can move along x-axis, y-axis and z-axis. And the particle move one step per time. The length of step is constant. Then, we can obtain these pictures.

Click the Code

Click the Code

Similarily, it given the same results as the situation in one dimension.
  • Now, if we consider that each step has different length, the result are given below, the difference is that the relation between the average of the x squared and time t.



  • In fact we can find a realer case, that is, the particle can move in all directions, and the length of one step has a range from 0 to 1.



In this case, the motion of particle is more ruleless and random. But in real environment, it is more accuracy. Besides, if we want to use python to get more exact result, we should consider the distribution of mean free path.

6. Random walks and diffusion

consider there is a cup of water, then, we put into a pack of coffee in the center of the water from statistic, we know the density of coffee is then proportional to the probability per unit volume per unit time. From the probability P, we can derive diffusion fuction from random walk.


which leads to


Then


where


If we consider one dimension, we have


  • Consider a particle is at the origin when t=0, then, we study its distribution with time. The following reflected time evolution calculated from the diffusion equation in one dimension.

Click the Code

A curious feature of the results for t>0 is that the density alternates between zero and nonzero values. This behavior is due to the initial density profile which we assumed.
  • We can average the results for adjacent spatial sites so as to "smooth over" the points where the density was zero. Doing this yields the results in the following.


  • I choosed some certain time to plot figures.


  • Consider two dimension case, yielding figures below.

Click the Code

In any case, the random-walk results exihibit the Gaussian spreading of the particle distribution.

7. Conclusion

1.For one, two or three dimension, if the particle has same probability move forward to all possible direction, the mean value of the square of the displacement is in direct proportion to the time t.
2. If there are different chances moving to different direction, the average of the x squared and time t into a quadratic relationship.
3. Python provides us all kinds of tools to accomplish our results, we should learn how to apply random number to them.
4. In any case, the random-walk results exihibit the Gaussian spreading of the particle distribution.

8. References

[1] Computational Physics(Second Edition), Nicholas J.Giordano, Hisao Nakanishi.

最后編輯于
?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請(qǐng)聯(lián)系作者
  • 序言:七十年代末煌茴,一起剝皮案震驚了整個(gè)濱河市宏悦,隨后出現(xiàn)的幾起案子挎春,更是在濱河造成了極大的恐慌仑嗅,老刑警劉巖,帶你破解...
    沈念sama閱讀 216,402評(píng)論 6 499
  • 序言:濱河連續(xù)發(fā)生了三起死亡事件痹束,死亡現(xiàn)場(chǎng)離奇詭異检疫,居然都是意外死亡,警方通過(guò)查閱死者的電腦和手機(jī)参袱,發(fā)現(xiàn)死者居然都...
    沈念sama閱讀 92,377評(píng)論 3 392
  • 文/潘曉璐 我一進(jìn)店門(mén)电谣,熙熙樓的掌柜王于貴愁眉苦臉地迎上來(lái),“玉大人抹蚀,你說(shuō)我怎么就攤上這事剿牺。” “怎么了环壤?”我有些...
    開(kāi)封第一講書(shū)人閱讀 162,483評(píng)論 0 353
  • 文/不壞的土叔 我叫張陵晒来,是天一觀的道長(zhǎng)。 經(jīng)常有香客問(wèn)我郑现,道長(zhǎng)湃崩,這世上最難降的妖魔是什么荧降? 我笑而不...
    開(kāi)封第一講書(shū)人閱讀 58,165評(píng)論 1 292
  • 正文 為了忘掉前任,我火速辦了婚禮攒读,結(jié)果婚禮上朵诫,老公的妹妹穿的比我還像新娘。我一直安慰自己薄扁,他們只是感情好剪返,可當(dāng)我...
    茶點(diǎn)故事閱讀 67,176評(píng)論 6 388
  • 文/花漫 我一把揭開(kāi)白布。 她就那樣靜靜地躺著邓梅,像睡著了一般脱盲。 火紅的嫁衣襯著肌膚如雪。 梳的紋絲不亂的頭發(fā)上日缨,一...
    開(kāi)封第一講書(shū)人閱讀 51,146評(píng)論 1 297
  • 那天钱反,我揣著相機(jī)與錄音,去河邊找鬼匣距。 笑死面哥,一個(gè)胖子當(dāng)著我的面吹牛,可吹牛的內(nèi)容都是我干的墨礁。 我是一名探鬼主播幢竹,決...
    沈念sama閱讀 40,032評(píng)論 3 417
  • 文/蒼蘭香墨 我猛地睜開(kāi)眼耳峦,長(zhǎng)吁一口氣:“原來(lái)是場(chǎng)噩夢(mèng)啊……” “哼恩静!你這毒婦竟也來(lái)了?” 一聲冷哼從身側(cè)響起蹲坷,我...
    開(kāi)封第一講書(shū)人閱讀 38,896評(píng)論 0 274
  • 序言:老撾萬(wàn)榮一對(duì)情侶失蹤驶乾,失蹤者是張志新(化名)和其女友劉穎,沒(méi)想到半個(gè)月后循签,有當(dāng)?shù)厝嗽跇?shù)林里發(fā)現(xiàn)了一具尸體级乐,經(jīng)...
    沈念sama閱讀 45,311評(píng)論 1 310
  • 正文 獨(dú)居荒郊野嶺守林人離奇死亡,尸身上長(zhǎng)有42處帶血的膿包…… 初始之章·張勛 以下內(nèi)容為張勛視角 年9月15日...
    茶點(diǎn)故事閱讀 37,536評(píng)論 2 332
  • 正文 我和宋清朗相戀三年县匠,在試婚紗的時(shí)候發(fā)現(xiàn)自己被綠了风科。 大學(xué)時(shí)的朋友給我發(fā)了我未婚夫和他白月光在一起吃飯的照片。...
    茶點(diǎn)故事閱讀 39,696評(píng)論 1 348
  • 序言:一個(gè)原本活蹦亂跳的男人離奇死亡乞旦,死狀恐怖贼穆,靈堂內(nèi)的尸體忽然破棺而出,到底是詐尸還是另有隱情兰粉,我是刑警寧澤故痊,帶...
    沈念sama閱讀 35,413評(píng)論 5 343
  • 正文 年R本政府宣布,位于F島的核電站玖姑,受9級(jí)特大地震影響愕秫,放射性物質(zhì)發(fā)生泄漏慨菱。R本人自食惡果不足惜,卻給世界環(huán)境...
    茶點(diǎn)故事閱讀 41,008評(píng)論 3 325
  • 文/蒙蒙 一戴甩、第九天 我趴在偏房一處隱蔽的房頂上張望符喝。 院中可真熱鬧,春花似錦甜孤、人聲如沸洲劣。這莊子的主人今日做“春日...
    開(kāi)封第一講書(shū)人閱讀 31,659評(píng)論 0 22
  • 文/蒼蘭香墨 我抬頭看了看天上的太陽(yáng)囱稽。三九已至,卻和暖如春二跋,著一層夾襖步出監(jiān)牢的瞬間战惊,已是汗流浹背。 一陣腳步聲響...
    開(kāi)封第一講書(shū)人閱讀 32,815評(píng)論 1 269
  • 我被黑心中介騙來(lái)泰國(guó)打工扎即, 沒(méi)想到剛下飛機(jī)就差點(diǎn)兒被人妖公主榨干…… 1. 我叫王不留吞获,地道東北人。 一個(gè)月前我還...
    沈念sama閱讀 47,698評(píng)論 2 368
  • 正文 我出身青樓谚鄙,卻偏偏與公主長(zhǎng)得像各拷,于是被迫代替她去往敵國(guó)和親。 傳聞我的和親對(duì)象是個(gè)殘疾皇子闷营,可洞房花燭夜當(dāng)晚...
    茶點(diǎn)故事閱讀 44,592評(píng)論 2 353

推薦閱讀更多精彩內(nèi)容

  • 太過(guò)突然傻盟,太過(guò)欣喜速蕊,人們卻往往害怕他的真實(shí),會(huì)是在做夢(mèng)嗎娘赴?伸出雙手规哲,竊竊的觸摸自己的夢(mèng),擰過(guò)耳垂诽表,扶過(guò)發(fā)絲唉锌,掠過(guò)鼻...
    我有一顆小鈴鐺閱讀 340評(píng)論 1 0
  • 加班中的我(?_?)
    Melvin寯鴻閱讀 143評(píng)論 0 0
  • 從前吃飯,講究食不言寢不語(yǔ)竿奏,立如松臥如龍袄简。貌似連吃相都是有講究的,是不能給外人看的议双,不然要挨橫眼痘番。 小時(shí)上學(xué)怕極了...
    思敏_e2ec閱讀 297評(píng)論 0 0