第二次用turtle作图：相对论

额外减小

如题
该程序模拟水星近日点进动现象
原理：狭义相对论质增效应
具体见代码
c 表示光速（不要低于行星起始速度！）

1. # 水星近日点进动 | 相对论质增效应
   
2. # Jianbing_Juan
   
3. 
   
4. import turtle as t
   
5. import math as m
   
6. 
   
7. def _angle(z):
   
8.     if z.real>0:
   
9.         return m.atan(z.imag/z.real)
   
10.     elif z.real<0:
    
11.         return m.atan(z.imag/z.real)+m.pi
    
12.     else:
    
13.         if z.imag>0:
    
14.             return m.pi/2
    
15.         elif z.imag<0:
    
16.             return -m.pi/2
    
17. 
    
18. def F(z):
    
19.     return -G*Mass*mass/pow(abs(z),2)*m.cos(_angle(z))-G*Mass*mass/pow(abs(z),2)*m.sin(_angle(z))*1j
    
20. 
    
21. mass=100
    
22. mass_v=mass
    
23. Mass=100000
    
24. G=1
    
25. c=500
    
26. 
    
27. st_r=-100+100j
    
28. st_v=20-5j
    
29. tstep=0.00005
    
30. r=st_r
    
31. v=st_v
    
32. t.up()
    
33. t.goto(r.real,r.imag)
    
34. t.down()
    
35. while True:
    
36.     for i in range(200):
    
37.         r+=v*tstep
    
38.         v+=F(r)/mass_v*tstep
    
39.         mass_v=mass/m.sqrt(1-(abs(v)*abs(v))/(c*c))
    
40.     t.goto(r.real,r.imag)

复制代码

额外减小

公式如下

其中 m 为本征质量
m' 为非本征质量
如有谬误还请大佬指正，谢谢
@Mike_python小 @hveagle @sfqxx @高山 @liuhongrun2022 看看吗

歌者文明清理员

额外减小 发表于 2023-3-4 17:49
公式如下

其中 m 为本征质量

不叫我？

sfqxx

育碧

wyhpylc

草（三声）我物理本来就不好，而且我只有初二啊，没学这玩意

额外减小

歌者文明清理员 发表于 2023-3-4 17:55
不叫我？

az
最多我只能叫五个 所以我就胡乱选择了前五个
你放心，下次一定叫你

额外减小

wyhpylc 发表于 2023-3-4 20:46
草（三声）我物理本来就不好，而且我只有初二啊，没学这玩意

加油

歌者文明清理员

@元豪 @KeyError @学习编程中的Ben @zhangjinxuan @ExiaGN001

额外减小

歌者文明清理员 发表于 2023-3-5 07:18
@元豪 @KeyError @学习编程中的Ben @zhangjinxuan @ExiaGN001

az...
谢谢你

ExiaGN001

歌者文明清理员 发表于 2023-3-5 07:18
@元豪 @KeyError @学习编程中的Ben @zhangjinxuan @ExiaGN001

py的内容就别喊我过来了，我不会

liuhongrun2022

康康

学习编程中的Ben

力鼎！