很容易对二元或三元轨道进行蒙特卡罗模拟。
orbitutils的Python项目详细描述
该软件包允许对轨道总体进行简单快速的蒙特卡罗模拟，给出投影分离、组件间相对rvs等的即时分布，包括偏心轨道的计算。便于模拟成像或RV观察。这段代码是组成Morton (2012)验证过程的模块之一。
偶尔使用astropyQuantity对象跟踪一些数量，使用astropySkyCoord对象跟踪坐标。还使用pandas保存/恢复填充。
请参见下面的快速介绍，以及notebook demo了解更多信息。
安装
$ pip install [--user] orbitutils
或者克隆存储库并安装：
$ git clone https://github.com/timothydmorton/orbitutils.git $ cd orbitutils $ python setup.py install [--user]
基本用法
模拟给定的主质量、次质量和轨道周期的总体。偏心率默认为零，但可以设置。
fromorbitutilsimportOrbitPopulationpop=OrbitPopulation(1,1,1000,n=1e4)#Primary mass, secondary mass, orbital period (d)
还可以创建次质量分布和/或周期分布：
importnumpyasnpfromorbitutilsimportOrbitPopulationN=1e4M2s=np.linspace(0.1,1,N)Ps=np.logspace(1,3,N)#periods evenly log-spaced from 10 to 1000 dayspop=OrbitPopulation(1,M2s,Ps)
OrbitPopulation对象可以使用save_hdf方法保存：
pop.save_hdf('mypop.h5')
然后按如下方式读入/还原：
fromorbitutilsimportOrbitPopulation_FromH5pop_again=OrbitPopulation_FromH5('mypop.h5')
归因
如果你用这个软件做研究，请告诉我。最终会有一个更新的出版物可供引用；在此之前，请引用我的2012年星球验证方法论文，该论文使用了以下代码：
@ARTICLE{2012ApJ...761....6M,
author = {{Morton}, T.~D.},
title = "{An Efficient Automated Validation Procedure for Exoplanet Transit Candidates}",
journal = {\apj},
archivePrefix = "arXiv",
eprint = {1206.1568},
primaryClass = "astro-ph.EP",
keywords = {planetary systems, stars: statistics },
year = 2012,
month = dec,
volume = 761,
eid = {6},
pages = {6},
doi = {10.1088/0004-637X/761/1/6},
adsurl = {http://adsabs.harvard.edu/abs/2012ApJ...761....6M},
adsnote = {Provided by the SAO/NASA Astrophysics Data System}
}
欢迎加入QQ群-->： 979659372
