【vpython模拟天体运动系列】公告

歌者文明清理员

成立了一个淘专辑去看看，希望鱼油们积极订阅
注意事项：
    1. 请先阅读淘专辑简介
因为整个程序在不断优化和改进中，所以每一期的代码非最终代码，所以现在不能放出代码

额外减小

额外减小 发表于 2023-3-5 15:08
写好了。这是模拟星体碎裂后绕中心天体绕转最后同时回到同一点的程序。
想知道怎么证明的扣6

1. from turtle import *
   
2. from math import *
   
3. 
   
4. M_center=100000
   
5. G=1
   
6. tstep=0.1/20
   
7. list_size=4
   
8. mlist=[10,10,10,10]
   
9. rlist=[100+100j,100+100j,100+100j,100+100j]
   
10. vlist=[10-20j,-10+20j,5+15j,-5-15j]
    
11. tlist=[Turtle(),Turtle(),Turtle(),Turtle()]
    
12. for i in range(list_size):
    
13.     tlist[i].speed(1000)
    
14.     tlist[i].up()
    
15.     tlist[i].goto(rlist[i].real,rlist[i].imag)
    
16.     tlist[i].down()
    
17. 
    
18. def _angle(z):
    
19.     if z.real>0:
    
20.         return atan(z.imag/z.real)
    
21.     elif z.real<0:
    
22.         return atan(z.imag/z.real)+pi
    
23.     else:
    
24.         if z.imag>=0:
    
25.             return pi/2
    
26.         else:
    
27.             return -pi/2
    
28. def _F(t_n):
    
29.     # |F|=G*M_center*mlist[n]/(abs(rlist[n])**2)
    
30.     # return : |F|*cos(_angle(rlist[n]))+|F|*sin(_angle(rlist[n]))*1j
    
31.     return -G*M_center*mlist[t_n]/(abs(rlist[t_n])**2)*cos(_angle(rlist[t_n]))-G*M_center*mlist[t_n]/(abs(rlist[t_n])**2)*sin(_angle(rlist[t_n]))*1j
    
32. # main
    
33. while True:
    
34.     for i in range(20):
    
35.         for j in range(list_size):
    
36.             rlist[j]+=vlist[j]*tstep
    
37.             tlist[j].goto(rlist[j].real,rlist[j].imag)
    
38.             vlist[j]+=_F(j)/mlist[j]*tstep
    
39. 
    

复制代码

sfqxx

我第一个耶

歌者文明清理员

sfqxx 发表于 2023-3-5 11:41
我第一个耶

算了我直说吧
请要订阅的鱼油把希望模拟的内容放在我的留言板里

歌者文明清理员

sfqxx 发表于 2023-3-5 11:41
我第一个耶

想当维护员是吧
成全

sfqxx

歌者文明清理员 发表于 2023-3-5 11:47
想当维护员是吧
成全

蟹蟹

歌者文明清理员

怎么一下就 20 查看量了？好奇怪

歌者文明清理员

注意，有时候电脑不在身边不能模拟，请众鱼油原谅
今天趁着可以出，大家赶紧想想有啥要模拟的，说不定还能写一起呢

额外减小

歌者文明清理员 发表于 2023-3-5 12:00
注意，有时候电脑不在身边不能模拟，请众鱼油原谅
今天趁着可以出，大家赶紧想想有啥要模拟的，说不定还能 ...

当然要模拟，史上最最经典的：
天体碎裂轨迹模拟！
超级的好看，超级的美妙

额外减小

额外减小 发表于 2023-3-5 13:36
当然要模拟，史上最最经典的：
天体碎裂轨迹模拟！
超级的好看，超级的美妙

用的动量守恒定律（）具体还没想好怎么做
待我去写写看

歌者文明清理员

额外减小 发表于 2023-3-5 13:36
当然要模拟，史上最最经典的：
天体碎裂轨迹模拟！
超级的好看，超级的美妙

我想想……不过请把请求发在留言板

额外减小

歌者文明清理员 发表于 2023-3-5 13:40
我想想……不过请把请求发在留言板

好的

额外减小

写好了。这是模拟星体碎裂后绕中心天体绕转最后同时回到同一点的程序。
想知道怎么证明的扣6

Mike_python小

可以模拟宇宙中所有的天梯运动轨迹吗

歌者文明清理员

Mike_python小 发表于 2023-3-5 19:38
可以模拟宇宙中所有的天梯运动轨迹吗

这个

一点沙

歌者文明清理员 发表于 2023-3-5 12:00
注意，有时候电脑不在身边不能模拟，请众鱼油原谅
今天趁着可以出，大家赶紧想想有啥要模拟的，说不定还能 ...

模拟一下黑洞呗

歌者文明清理员

一点沙 发表于 2023-3-12 14:32
模拟一下黑洞呗

ok

歌者文明清理员

一点沙 发表于 2023-3-12 14:32
模拟一下黑洞呗

由于怕被发现停更一周

一点沙

歌者文明清理员 发表于 2023-3-12 15:56
由于怕被发现停更一周

啥意思

liuhongrun2022

一点沙 发表于 2023-3-12 18:27
啥意思

怕被家长发现