Zongbo Huyan, Dapeng Zhang, Pengbin Ma, Hengnian Li, Zhai Min
{"title":"天文一号轨道确定的并行观测过程","authors":"Zongbo Huyan, Dapeng Zhang, Pengbin Ma, Hengnian Li, Zhai Min","doi":"10.1515/astro-2022-0202","DOIUrl":null,"url":null,"abstract":"To meet various requirements of the future deep space missions of China, State Key Laboratory of Astronautic Dynamics constructs a new orbit determinatin software with parallel observations process. By using 32 threads, the computational efficiency per iteration (including spacecraft’s integration) could be promoted to 10 times as that of single-threaded orbit determination. Suppose the number of observations is <jats:inline-formula> <jats:alternatives> <jats:inline-graphic xmlns:xlink=\"http://www.w3.org/1999/xlink\" xlink:href=\"graphic/j_astro-2022-0202_eq_001.png\"/> <m:math xmlns:m=\"http://www.w3.org/1998/Math/MathML\"> <m:mi>p</m:mi> </m:math> <jats:tex-math>p</jats:tex-math> </jats:alternatives> </jats:inline-formula>, the number of estimated parameters (including spacecraft’s state) is <jats:inline-formula> <jats:alternatives> <jats:inline-graphic xmlns:xlink=\"http://www.w3.org/1999/xlink\" xlink:href=\"graphic/j_astro-2022-0202_eq_002.png\"/> <m:math xmlns:m=\"http://www.w3.org/1998/Math/MathML\"> <m:mi>q</m:mi> </m:math> <jats:tex-math>q</jats:tex-math> </jats:alternatives> </jats:inline-formula>, the amount of computation of one observation is <jats:inline-formula> <jats:alternatives> <jats:inline-graphic xmlns:xlink=\"http://www.w3.org/1999/xlink\" xlink:href=\"graphic/j_astro-2022-0202_eq_003.png\"/> <m:math xmlns:m=\"http://www.w3.org/1998/Math/MathML\"> <m:mi>x</m:mi> </m:math> <jats:tex-math>x</jats:tex-math> </jats:alternatives> </jats:inline-formula>, the amount of computation of one Givens transformation is <jats:inline-formula> <jats:alternatives> <jats:inline-graphic xmlns:xlink=\"http://www.w3.org/1999/xlink\" xlink:href=\"graphic/j_astro-2022-0202_eq_004.png\"/> <m:math xmlns:m=\"http://www.w3.org/1998/Math/MathML\"> <m:mi>y</m:mi> </m:math> <jats:tex-math>y</jats:tex-math> </jats:alternatives> </jats:inline-formula> and the best number of threads is proved to be <jats:inline-formula> <jats:alternatives> <jats:inline-graphic xmlns:xlink=\"http://www.w3.org/1999/xlink\" xlink:href=\"graphic/j_astro-2022-0202_eq_005.png\"/> <m:math xmlns:m=\"http://www.w3.org/1998/Math/MathML\"> <m:msup> <m:mrow> <m:mrow> <m:mo>[</m:mo> <m:mrow> <m:mi>p</m:mi> <m:mo>⋅</m:mo> <m:mrow> <m:mo>(</m:mo> <m:mrow> <m:mi>x</m:mi> <m:mo>+</m:mo> <m:mi>y</m:mi> </m:mrow> <m:mo>)</m:mo> </m:mrow> <m:mspace width=\"0.1em\"/> <m:mtext>/</m:mtext> <m:mspace width=\"0.1em\"/> <m:mrow> <m:mo>(</m:mo> <m:mrow> <m:mi>q</m:mi> <m:mo>⋅</m:mo> <m:mi>y</m:mi> </m:mrow> <m:mo>)</m:mo> </m:mrow> </m:mrow> <m:mo>]</m:mo> </m:mrow> </m:mrow> <m:mrow> <m:mn>1</m:mn> <m:mspace width=\"0.1em\"/> <m:mtext>/</m:mtext> <m:mspace width=\"0.1em\"/> <m:mn>2</m:mn> </m:mrow> </m:msup> </m:math> <jats:tex-math>{\\left[p\\cdot \\left(x+y)\\text{/}\\left(q\\cdot y)]}^{1\\text{/}2}</jats:tex-math> </jats:alternatives> </jats:inline-formula> for one-step threads combination. The root mean square of the postfit residuals of China’s deep space monitoring network (CDSMN) and China’s Very Long Baseline Interferometry (VLBI) network (CVN) observations in the Earth-Mars transfer phase and the Mars-orbiting phase are almost the same: about 0.3 m for Ranging, about 0.3 mm/s for Doppler, about 3 cm for VLBI delay and about 0.5 mm/s for VLBI delay rate. It is also found that all the four types of observations of CDSMN and CVN are needed in orbit determination for deep space maneuver and braking at periareion calculation. In the Mars-orbiting phase, the position accuracy after orbit determination under CDSMN-only tracking mode can reach about 1 km.","PeriodicalId":19514,"journal":{"name":"Open Astronomy","volume":"33 1","pages":""},"PeriodicalIF":0.5000,"publicationDate":"2024-04-26","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":"{\"title\":\"Parallel observations process of Tianwen-1 orbit determination\",\"authors\":\"Zongbo Huyan, Dapeng Zhang, Pengbin Ma, Hengnian Li, Zhai Min\",\"doi\":\"10.1515/astro-2022-0202\",\"DOIUrl\":null,\"url\":null,\"abstract\":\"To meet various requirements of the future deep space missions of China, State Key Laboratory of Astronautic Dynamics constructs a new orbit determinatin software with parallel observations process. By using 32 threads, the computational efficiency per iteration (including spacecraft’s integration) could be promoted to 10 times as that of single-threaded orbit determination. Suppose the number of observations is <jats:inline-formula> <jats:alternatives> <jats:inline-graphic xmlns:xlink=\\\"http://www.w3.org/1999/xlink\\\" xlink:href=\\\"graphic/j_astro-2022-0202_eq_001.png\\\"/> <m:math xmlns:m=\\\"http://www.w3.org/1998/Math/MathML\\\"> <m:mi>p</m:mi> </m:math> <jats:tex-math>p</jats:tex-math> </jats:alternatives> </jats:inline-formula>, the number of estimated parameters (including spacecraft’s state) is <jats:inline-formula> <jats:alternatives> <jats:inline-graphic xmlns:xlink=\\\"http://www.w3.org/1999/xlink\\\" xlink:href=\\\"graphic/j_astro-2022-0202_eq_002.png\\\"/> <m:math xmlns:m=\\\"http://www.w3.org/1998/Math/MathML\\\"> <m:mi>q</m:mi> </m:math> <jats:tex-math>q</jats:tex-math> </jats:alternatives> </jats:inline-formula>, the amount of computation of one observation is <jats:inline-formula> <jats:alternatives> <jats:inline-graphic xmlns:xlink=\\\"http://www.w3.org/1999/xlink\\\" xlink:href=\\\"graphic/j_astro-2022-0202_eq_003.png\\\"/> <m:math xmlns:m=\\\"http://www.w3.org/1998/Math/MathML\\\"> <m:mi>x</m:mi> </m:math> <jats:tex-math>x</jats:tex-math> </jats:alternatives> </jats:inline-formula>, the amount of computation of one Givens transformation is <jats:inline-formula> <jats:alternatives> <jats:inline-graphic xmlns:xlink=\\\"http://www.w3.org/1999/xlink\\\" xlink:href=\\\"graphic/j_astro-2022-0202_eq_004.png\\\"/> <m:math xmlns:m=\\\"http://www.w3.org/1998/Math/MathML\\\"> <m:mi>y</m:mi> </m:math> <jats:tex-math>y</jats:tex-math> </jats:alternatives> </jats:inline-formula> and the best number of threads is proved to be <jats:inline-formula> <jats:alternatives> <jats:inline-graphic xmlns:xlink=\\\"http://www.w3.org/1999/xlink\\\" xlink:href=\\\"graphic/j_astro-2022-0202_eq_005.png\\\"/> <m:math xmlns:m=\\\"http://www.w3.org/1998/Math/MathML\\\"> <m:msup> <m:mrow> <m:mrow> <m:mo>[</m:mo> <m:mrow> <m:mi>p</m:mi> <m:mo>⋅</m:mo> <m:mrow> <m:mo>(</m:mo> <m:mrow> <m:mi>x</m:mi> <m:mo>+</m:mo> <m:mi>y</m:mi> </m:mrow> <m:mo>)</m:mo> </m:mrow> <m:mspace width=\\\"0.1em\\\"/> <m:mtext>/</m:mtext> <m:mspace width=\\\"0.1em\\\"/> <m:mrow> <m:mo>(</m:mo> <m:mrow> <m:mi>q</m:mi> <m:mo>⋅</m:mo> <m:mi>y</m:mi> </m:mrow> <m:mo>)</m:mo> </m:mrow> </m:mrow> <m:mo>]</m:mo> </m:mrow> </m:mrow> <m:mrow> <m:mn>1</m:mn> <m:mspace width=\\\"0.1em\\\"/> <m:mtext>/</m:mtext> <m:mspace width=\\\"0.1em\\\"/> <m:mn>2</m:mn> </m:mrow> </m:msup> </m:math> <jats:tex-math>{\\\\left[p\\\\cdot \\\\left(x+y)\\\\text{/}\\\\left(q\\\\cdot y)]}^{1\\\\text{/}2}</jats:tex-math> </jats:alternatives> </jats:inline-formula> for one-step threads combination. The root mean square of the postfit residuals of China’s deep space monitoring network (CDSMN) and China’s Very Long Baseline Interferometry (VLBI) network (CVN) observations in the Earth-Mars transfer phase and the Mars-orbiting phase are almost the same: about 0.3 m for Ranging, about 0.3 mm/s for Doppler, about 3 cm for VLBI delay and about 0.5 mm/s for VLBI delay rate. It is also found that all the four types of observations of CDSMN and CVN are needed in orbit determination for deep space maneuver and braking at periareion calculation. In the Mars-orbiting phase, the position accuracy after orbit determination under CDSMN-only tracking mode can reach about 1 km.\",\"PeriodicalId\":19514,\"journal\":{\"name\":\"Open Astronomy\",\"volume\":\"33 1\",\"pages\":\"\"},\"PeriodicalIF\":0.5000,\"publicationDate\":\"2024-04-26\",\"publicationTypes\":\"Journal Article\",\"fieldsOfStudy\":null,\"isOpenAccess\":false,\"openAccessPdf\":\"\",\"citationCount\":\"0\",\"resultStr\":null,\"platform\":\"Semanticscholar\",\"paperid\":null,\"PeriodicalName\":\"Open Astronomy\",\"FirstCategoryId\":\"101\",\"ListUrlMain\":\"https://doi.org/10.1515/astro-2022-0202\",\"RegionNum\":4,\"RegionCategory\":\"物理与天体物理\",\"ArticlePicture\":[],\"TitleCN\":null,\"AbstractTextCN\":null,\"PMCID\":null,\"EPubDate\":\"\",\"PubModel\":\"\",\"JCR\":\"Q4\",\"JCRName\":\"ASTRONOMY & ASTROPHYSICS\",\"Score\":null,\"Total\":0}","platform":"Semanticscholar","paperid":null,"PeriodicalName":"Open Astronomy","FirstCategoryId":"101","ListUrlMain":"https://doi.org/10.1515/astro-2022-0202","RegionNum":4,"RegionCategory":"物理与天体物理","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"Q4","JCRName":"ASTRONOMY & ASTROPHYSICS","Score":null,"Total":0}
引用次数: 0
摘要
为满足中国未来深空任务的各种要求,航天动力学国家重点实验室构建了一种新型并行观测过程的定轨软件。通过使用 32 个线程,每次迭代的计算效率(包括航天器的积分)可提高到单线程定轨的 10 倍。假设观测数据的数量为 p p,估计参数(包括航天器状态)的数量为 q q,一个观测数据的计算量为 x x 、一步线程组合的最佳线程数为 [ p ⋅ ( x + y ) / ( q ⋅ y ) ] 1 / 2 {\left[p\cdot \left(x+y)\text{/}left(q\cdot y)]}^{1/text{/}2}。中国深空监测网(CDSMN)和中国甚长基线干涉测量(VLBI)网(CVN)在地-火星转移阶段和火星入轨阶段的观测后拟合残差的均方根基本相同:Ranging约为0.3 m,Doppler约为0.3 mm/s,VLBI延迟约为3 cm,VLBI延迟率约为0.5 mm/s。研究还发现,在进行深空机动和近地轨道制动计算的轨道确定时,需要对 CDSMN 和 CVN 进行所有四种类型的观测。在火星轨道阶段,仅在 CDSMN 跟踪模式下确定轨道后的位置精度可达 1 公里左右。
Parallel observations process of Tianwen-1 orbit determination
To meet various requirements of the future deep space missions of China, State Key Laboratory of Astronautic Dynamics constructs a new orbit determinatin software with parallel observations process. By using 32 threads, the computational efficiency per iteration (including spacecraft’s integration) could be promoted to 10 times as that of single-threaded orbit determination. Suppose the number of observations is pp, the number of estimated parameters (including spacecraft’s state) is qq, the amount of computation of one observation is xx, the amount of computation of one Givens transformation is yy and the best number of threads is proved to be [p⋅(x+y)/(q⋅y)]1/2{\left[p\cdot \left(x+y)\text{/}\left(q\cdot y)]}^{1\text{/}2} for one-step threads combination. The root mean square of the postfit residuals of China’s deep space monitoring network (CDSMN) and China’s Very Long Baseline Interferometry (VLBI) network (CVN) observations in the Earth-Mars transfer phase and the Mars-orbiting phase are almost the same: about 0.3 m for Ranging, about 0.3 mm/s for Doppler, about 3 cm for VLBI delay and about 0.5 mm/s for VLBI delay rate. It is also found that all the four types of observations of CDSMN and CVN are needed in orbit determination for deep space maneuver and braking at periareion calculation. In the Mars-orbiting phase, the position accuracy after orbit determination under CDSMN-only tracking mode can reach about 1 km.
Open AstronomyPhysics and Astronomy-Astronomy and Astrophysics
CiteScore
1.30
自引率
14.30%
发文量
37
审稿时长
16 weeks
期刊介绍:
The journal disseminates research in both observational and theoretical astronomy, astrophysics, solar physics, cosmology, galactic and extragalactic astronomy, high energy particles physics, planetary science, space science and astronomy-related astrobiology, presenting as well the surveys dedicated to astronomical history and education.