本届研赛助攻题目 C D F三题论文均已经全部完成。后更新计划 如图所示。
免费给大家分享 三个问题的论文+部分代码
2024年华为杯-研赛分享资料（论文+部分代码）（已更新部分代码）：
链接：https://pan.baidu.com/s/1HGIYjV3lqzUc_3H0vg5H8w
提取码：sxjm
X射线脉冲星光子到达时间建模
摘要
脉冲星是一类高速自转的中子星，其自转形成规律性脉冲信号，类似于“宇宙中的灯塔”，因此被认为是极为精确的时钟。X射线脉冲星导航利用脉冲星信号为航天器提供时间和空间参考。通过比较脉冲信号到达航天器和太阳系质心的时间差，能够实现航天器的精确定位。为了准确计算该时间差，需要考虑卫星的轨道运动、脉冲星光子传播过程中的几何时延、Shapiro时延、引力红移效应及相对论效应等多种复杂因素。
2016年，我国发射的XPNAV-1卫星观测了蟹状星云脉冲星（PSR B0531+21）的X射线信号。该任务要求通过光子到达时间模拟和折叠，得到脉冲星的精确计时轮廓，为脉冲星导航系统的发展提供重要数据支持。
针对问题1，首先，利用给定的轨道根数 (偏心率 、角动量 、轨道倾角 、升交点赤经 、近地点幅角 、和真近点角 )，通过公式计算卫星在轨道平面内的径向距离 以及径向速度 和切向速度 。接着，通过旋转矩阵将卫星的二维轨道位置和速度从轨道平面转换到地心天球参考系 (GCRS) 中的三维位置和速度。最终，计算出卫星在特定时刻的三维位置 和速度
一、模型假设
为了方便模型的建立与模型的可行性，我们这里首先对模型提出一些假设，使得模型更加完备，预测的结果更加合理。
1.假设给出的数据均为真实数据，真实有效。
2.假设卫星的运动符合经典轨道力学，遵循开普勒定律，轨道为椭圆形。不考虑摄动效应，如地球的非球形引力场（J2项）、大气阻力、月球和太阳的引力影响等。
3.假设卫星的质量远小于地球质量，因此地球的引力场完全主导卫星的轨道，卫星的质量对轨道运动无影响。
4.轨道的椭圆性由偏心率决定，轨道参数是理想化的，不受任何其他外部因素的影响（如大气摩擦、引力摄动等）。
5.假设脉冲星发出的光子在空间中是平行传播的（即远距离的光子传播路径被近似为平行线）。
6.假设光子从脉冲星到卫星及太阳系质心的传播速度为常数 ccc，并且没有受到介质的阻挡和折射，光子传播路径不受天体间介质的影响。
7.假设光子传播过程中不受天体（如地球、太阳或其他天体）的引力作用影响，传播路径为直线。
一、符号说明
为了方便我们模型的建立与求解过程 ，我们这里对使用到的关键符号进行以下说明：
| 符号 | 符号说明 |
| 地球引力常数，单位 ，通常取 | |
| 轨道偏心率 (eccentricity)，描述椭圆轨道的扁平程度，取值范围为 | |
| 角动量（specific angular momentum），单位为 ，由公式 定义 | |
| 轨道倾角 (inclination)，单位为弧度或度，描述卫星轨道平面与地球赤道平面的夹角，取值范围为 | |
| 升交点赤经（right ascension of ascending node, RAAN），单位为弧度或度，描述升交点在天球赤道上的位置，取值范围为 | |
| 近地点幅角 (argument of perigee)，单位为弧度或度，描述轨道在轨道平面中的方向，取值范围为 |
| 结果项 | 计算值 | 单位 | 给定值 |
| 三维位置 (X, Y, Z) | (4199.74, 3428.29, 4242.30) | km | - |
| 三维速度 (Vx, Vy, Vz) | (4.36, 1.96, -5.91) | km/s | - |
| 角动量 h | 52330.85 | 52330.85 |
