第一次用Turtle绘图，做了个简单的小玩意

额外减小

如题
演绎三体运动的程序
效果：

代码：

1. import turtle
   
2. import math as ma
   
3. 
   
4. #
   
5. tstep=0.002
   
6. #
   
7. def _angle(z):
   
8.     if z.real>0:
   
9.         return ma.atan(z.imag/z.real)
   
10.     elif z.real<0:
    
11.         return ma.atan(z.imag/z.real)+ma.pi
    
12.     else:
    
13.         if z.imag>=0:
    
14.             return ma.pi/2
    
15.         else:
    
16.             return -ma.pi/2
    
17. 
    
18. #参数列表
    
19. mlist=[5000,10000,5000]
    
20. rlist=[0-30j,-20+10j,10+20j]
    
21. vlist=[20+0j,-10+5j,10+0j]
    
22. t_0,t_1,t_2=turtle.Turtle(),turtle.Turtle(),turtle.Turtle()
    
23. G=0.1
    
24. #
    
25. 
    
26. t_0.up()
    
27. t_1.up()
    
28. t_2.up()
    
29. #color
    
30. t_0.color("red")
    
31. t_1.color("green")
    
32. t_2.color("blue")
    
33. #
    
34. t_0.goto(rlist[0].real,rlist[0].imag)
    
35. t_1.goto(rlist[1].real,rlist[1].imag)
    
36. t_2.goto(rlist[2].real,rlist[2].imag)
    
37. t_0.down()
    
38. t_1.down()
    
39. t_2.down()
    
40. 
    
41. def F(a,b):
    
42.     return G*mlist[a]*mlist[b]/abs(rlist[a]-rlist[b])*ma.cos(_angle(rlist[b]-rlist[a]))+G*mlist[a]*mlist[b]/abs(rlist[a]-rlist[b])*ma.sin(_angle(rlist[b]-rlist[a]))*1j
    
43. 
    
44. while True:
    
45.     for i in range(5):
    
46.         vlist[0]+=(F(0,1)+F(0,2))/mlist[0]*tstep
    
47.         vlist[1]+=(F(1,0)+F(1,2))/mlist[1]*tstep
    
48.         vlist[2]+=(F(2,0)+F(2,1))/mlist[2]*tstep
    
49.         rlist[0]+=vlist[0]*tstep
    
50.         rlist[1]+=vlist[1]*tstep
    
51.         rlist[2]+=vlist[2]*tstep
    
52.    
    
53.     t_0.goto(rlist[0].real,rlist[0].imag)
    
54.     t_1.goto(rlist[1].real,rlist[1].imag)
    
55.     t_2.goto(rlist[2].real,rlist[2].imag)

复制代码

一点沙

好厉害

额外减小

可更改参数：
line 5 ： 每次演绎间隔时间
line 19： 三个质点的质量
line 20： 三个质点的初始位矢
line 21： 三个质点的初始速度矢量
line 23： 万有引力常量
line 30-32： turtle 的 颜色

歌者文明清理员

                               
登录/注册后可看大图

哈哈

额外减小

歌者文明清理员 发表于 2023-2-18 15:58
哈哈

ha

歌者文明清理员

额外减小 发表于 2023-2-18 17:27
ha

像极了题目不会做的时候我抓头皮屑的路线

额外减小

歌者文明清理员 发表于 2023-2-18 20:26
像极了题目不会做的时候我抓头皮屑的路线

哈哈

额外减小

绿色星球想掺和到红色和蓝色星球的紧密联系中
奈何质量太小，无足轻重，只能紧随其后

歌者文明清理员

额外减小 发表于 2023-2-18 21:30
绿色星球想掺和到红色和蓝色星球的紧密联系中
奈何质量太小，无足轻重，只能紧随其后

这，万一球跑到右边了，可以设置按下某个键移动视角吗

歌者文明清理员

等质二体+行星

1. import turtle
   
2. import math as ma
   
3. import random as random
   
4. #
   
5. tstep=0.02
   
6. #
   
7. def _angle(z):
   
8.     if z.real>0:
   
9.         return ma.atan(z.imag/z.real)
   
10.     elif z.real<0:
    
11.         return ma.atan(z.imag/z.real)+ma.pi
    
12.     else:
    
13.         if z.imag>=0:
    
14.             return ma.pi/2
    
15.         else:
    
16.             return -ma.pi/2
    
17. 
    
18. def rand():return complex(random.uniform(-25, 25), random.uniform(-25, 25))
    
19. #参数列表
    
20. mlist=[5000,5000,0.00001]
    
21. rlist=[50+0j, -50+0j, 0+0j]
    
22. vlist=[0-16j]
    
23. vlist.insert(1, -vlist[0])
    
24. vlist.append(vlist[0]+(0+10j))
    
25. t_0,t_1,t_2=turtle.Turtle(),turtle.Turtle(),turtle.Turtle()
    
26. G=0.1
    
27. #
    
28. 
    
29. t_0.up()
    
30. t_1.up()
    
31. t_2.up()
    
32. #color
    
33. t_0.pencolor("red")
    
34. t_1.pencolor("green")
    
35. t_2.pencolor("blue")
    
36. t_0.fillcolor('yellow')
    
37. t_1.fillcolor('cyan')
    
38. t_2.fillcolor('orange')
    
39. #
    
40. t_0.goto(rlist[0].real,rlist[0].imag)
    
41. t_1.goto(rlist[1].real,rlist[1].imag)
    
42. t_2.goto(rlist[2].real,rlist[2].imag)
    
43. t_0.down()
    
44. t_1.down()
    
45. t_2.down()
    
46. 
    
47. def F(a,b):
    
48.     return G*mlist[a]*mlist[b]/abs(rlist[a]-rlist[b])*ma.cos(_angle(rlist[b]-rlist[a]))+G*mlist[a]*mlist[b]/abs(rlist[a]-rlist[b])*ma.sin(_angle(rlist[b]-rlist[a]))*1j
    
49. 
    
50. while True:
    
51.     for i in range(5):
    
52.         vlist[0]+=(F(0,1)+F(0,2))/mlist[0]*tstep
    
53.         vlist[1]+=(F(1,0)+F(1,2))/mlist[1]*tstep
    
54.         vlist[2]+=(F(2,0)+F(2,1))/mlist[2]*tstep
    
55.         rlist[0]+=vlist[0]*tstep
    
56.         rlist[1]+=vlist[1]*tstep
    
57.         rlist[2]+=vlist[2]*tstep
    
58.    
    
59.     t_0.goto(rlist[0].real,rlist[0].imag)
    
60.     t_1.goto(rlist[1].real,rlist[1].imag)
    
61.     t_2.goto(rlist[2].real,rlist[2].imag)

复制代码

额外减小

歌者文明清理员 发表于 2023-2-19 07:53
这，万一球跑到右边了，可以设置按下某个键移动视角吗

这...我试试

额外减小

歌者文明清理员 发表于 2023-2-19 08:20
等质二体+行星

这个更为绝佳：等质三体+行星（楽）

1. import turtle
   
2. import math as ma
   
3. 
   
4. tstep=0.0002
   
5. def _angle(z):
   
6.     if z.real>0:
   
7.         return ma.atan(z.imag/z.real)
   
8.     elif z.real<0:
   
9.         return ma.atan(z.imag/z.real)+ma.pi
   
10.     else:
    
11.         if z.imag>=0:
    
12.             return ma.pi/2
    
13.         else:
    
14.             return -ma.pi/2
    
15. 
    
16. mlist=[100000,100000,100000,1]
    
17. rlist=[0+60j,40+50j,-80-80j,0+0j]
    
18. vlist=[1+1j,0-5j,-5+2j,1-5j]
    
19. t_0,t_1,t_2,t_3=turtle.Turtle(),turtle.Turtle(),turtle.Turtle(),turtle.Turtle()
    
20. f=0
    
21. G=1
    
22. 
    
23. t_0.up()
    
24. t_1.up()
    
25. t_2.up()
    
26. t_3.up()
    
27. t_0.color("red")
    
28. t_1.color("green")
    
29. t_2.color("blue")
    
30. t_3.color("orange")
    
31. t_0.goto(rlist[0].real,rlist[0].imag)
    
32. t_1.goto(rlist[1].real,rlist[1].imag)
    
33. t_2.goto(rlist[2].real,rlist[2].imag)
    
34. t_3.goto(rlist[3].real,rlist[3].imag)
    
35. #t_0.down()
    
36. #t_1.down()
    
37. #t_2.down()
    
38. t_3.down()
    
39. 
    
40. def F(a,b):
    
41.     f=G*mlist[a]*mlist[b]/abs(rlist[a]-rlist[b])
    
42.     return f*ma.cos(_angle(rlist[b]-rlist[a]))+f*ma.sin(_angle(rlist[b]-rlist[a]))*1j
    
43. 
    
44. while True:
    
45.     for i in range(50):
    
46.         vlist[0]+=(F(0,1)+F(0,2)+F(0,3))/mlist[0]*tstep
    
47.         vlist[1]+=(F(1,0)+F(1,2)+F(1,3))/mlist[1]*tstep
    
48.         vlist[2]+=(F(2,0)+F(2,1)+F(2,3))/mlist[2]*tstep
    
49.         vlist[3]+=(F(3,0)+F(3,1)+F(3,2))/mlist[3]*tstep
    
50.         rlist[0]+=vlist[0]*tstep
    
51.         rlist[1]+=vlist[1]*tstep
    
52.         rlist[2]+=vlist[2]*tstep
    
53.         rlist[3]+=vlist[3]*tstep
    
54.     t_0.goto(rlist[0].real,rlist[0].imag)
    
55.     t_1.goto(rlist[1].real,rlist[1].imag)
    
56.     t_2.goto(rlist[2].real,rlist[2].imag)
    
57.     t_3.goto(rlist[3].real,rlist[3].imag)

复制代码

额外减小

目前正在研究关于tkinter的做法。

歌者文明清理员

额外减小 发表于 2023-2-20 20:00
这个更为绝佳：等质三体+行星（楽）

一般三体会分裂成二体和另外一颗恒星旋转或者弹出
这个是很稳的了

额外减小

歌者文明清理员 发表于 2023-2-20 20:26
一般三体会分裂成二体和另外一颗恒星旋转或者弹出
这个是很稳的了

你懂的真多，但我想问问为什么？

歌者文明清理员

额外减小 发表于 2023-2-20 20:31
你懂的真多，但我想问问为什么？

不知道，b 站很多三体模拟都是这样
除了三体的特解不会分裂

再说一句，turtle 的内核是 tkinter

额外减小

歌者文明清理员 发表于 2023-2-20 20:43
不知道，b 站很多三体模拟都是这样
除了三体的特解不会分裂

最后一句我知道。
不过你这个分裂的结论还需要理论依据呀。
比如用概率论证明啥的

歌者文明清理员

额外减小 发表于 2023-2-20 20:45
最后一句我知道。
不过你这个分裂的结论还需要理论依据呀。
比如用概率论证明啥的

啊哈视频里也有讲，自行百度试试？

额外减小

歌者文明清理员 发表于 2023-2-20 20:50
啊哈视频里也有讲，自行百度试试？

好
刚刚我试了一下，发现三颗日一起飞出屏幕

歌者文明清理员

额外减小 发表于 2023-2-20 20:58
好
刚刚我试了一下，发现三颗日一起飞出屏幕

《团结一致》