基于MATLAB的BDS卫星轨道仿真

 一、实现功能

本代码实现了以下几个功能

• 根据广播星历信息提取轨道六根数和摄动参数。

• 计算并绘制BDS GEO卫星、IGSO卫星和MEO卫星的无摄轨道平面。

• 计算并绘制MEO卫星在无摄全轨道平面的位置分布。

• 添加摄动参数，推算两小时之内IGSO卫星位置和速度变化，并绘制位置和速度随时间变化曲线。

 二、了解RINEX格式的BDS广播星历格式和内容

以下是一段北斗星历文件

        注意图中END OF HEADER，此字符串以上部分我们不做赘述，可自行查找，我们只需要考虑此字符串以下的部分。

        要想通过卫星星历文件提取数剧，我们就要了解星历格式每个部分的含义，观察可以发现，文档中的数据是一块一块的，观察每一块数据，可得，它是一个8×5的表格，具体每一个数据对应的意义可以看下表：

卫星编号	时间
		轨道半径改正项Crs	平均角速度改正项delta_n	平近点角M0
	升交点角距改正项Cuc	轨道偏心率e	升交点角距改正项Cus	轨道长半轴平方根sqrta
	星历的参考时刻t_oe	轨道倾角改正项Cic	升交点经度omega	轨道倾角改正项Cis
	轨道倾角i0