用Python编了一个轨道炮模拟器

实际上是解一个ODE

$[Q,L,R,Q',L',R']'=[Q',L',R',-RQ'/L-Q/CL,Lr^2Q'^2/2m,LrRrQ'^2/2m]$   

Lr是轨道电感梯度，Rr是轨道电阻梯度

这个方程由轨道炮弹丸受力公式

$F=Lr^2I^2/2$ 

推导得出（因为我可能会推错，所以希望有大佬能帮我验证一下。。。）

运行截图

轨道炮.py 1.19KB PY 82次下载