如何解决Sympy 与 astropy 结合使用单位进行数值求解?
有没有办法将 Sympy 的求解器与 Astropy 结合起来?例如,如果我想对所有其他变量的给定值求解开普勒第三定律:
from sympy import solveset,poly,Eq,Function,exp
from sympy import Symbol,pi
Ms = Symbol('Ms') # mass of sun
Mp = Symbol('Mp') # mass of planet
G = Symbol('G') # Gravitational consatant
a = Symbol('a') # semi jamor axes
P = Symbol('P') # period
solveset(G * (Ms + Mp) / 4 * pi**2 - a**3/P**2,a)
然后我如何使用 astropy 将具有相应单位的值添加到所有变量并获得数字输出?
import astropy.units as u
from astropy.constants import G
Ms = 1 * u.M_sun
Mp = 1 * u.M_jup
s = 5 * u.pc
P = 1 * u.yr
这是否可以与 Sympy 收到的结果相结合,以解决具有预期单位和数值结果的 a 中的值?
感谢任何建议。
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。