-
Notifications
You must be signed in to change notification settings - Fork 1
/
Brownian.py
74 lines (48 loc) · 1.5 KB
/
Brownian.py
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
#Python 实现布朗运动
#Brownian motion -- an example of a multi-threaded Tkinter program
from tkinter import *
import random
import threading
import time
import sys
#画布(窗口)大小
WIDTH = 400
HEIGHT = 300
#设置布朗粒子(小球)的相关参数
SIGMA = 50 #小球的分散程度
BUZZ = 10 #小球的扩散速度
RADIUS = 2 #小球半径(大小)
LAMBDA = 50 #小球的运动速度
FILL = 'green' #小球的颜色
stop = 0 # Set when main loop exits
def particle(canvas):
r = RADIUS
x = random.gauss(WIDTH/2.0, SIGMA) #随机生成符合高斯分布的随机数,WIDTH/2.0和SIGMA为高斯分布的两个参数
y = random.gauss(HEIGHT/2.0, SIGMA)
p = canvas.create_oval(x-r, y-r, x+r, y+r, fill = FILL)
while not stop:
dx = random.gauss(0, BUZZ)
dy = random.gauss(0, BUZZ)
dt = random.expovariate(LAMBDA)
try:
canvas.move(p, dx, dy)
except TclError:
break
time.sleep(dt)
def main(): #主函数
global stop
root = Tk()
canvas = Canvas(root, width = WIDTH, height = HEIGHT)
canvas.pack(fill = 'both', expand = 1)
#粒子数目
np = 30
if sys.argv[1:]:
np = int(sys.argv[1])
for i in range(np):
t = threading.Thread(target = particle, args = (canvas,))
t.start()
try:
root.mainloop()
finally:
stop = 1
main()