发明名称 一种高精度的导航卫星TGD参数标定方法
摘要 本发明公开了一种高精度的导航卫星TGD参数标定方法,直接修正导航卫星频间相位中心的垂向误差,对于水平误差,采用多站加权的方式,充分利用地面站的几何布局进行消除;参数标定过程不需要利用卫星质心位置以及质心与相位中心的位置关系,求解相位中心的位置,而且多站加权处理对卫星星下点坐标、电离层穿刺点坐标的精度要求不高,因此TGD参数标定精度不受卫星质心位置的工厂标定误差、在轨卫星的质心在卫星机械坐标系内的变化、卫星姿态测量误差以及卫星轨道测量误差的影响。将导航卫星频间相位中心分离引起的TGD估计误差降低至原来的5%~15%。根据相位中心实际分离的大小,可将相位中心分离引起的TGD估计误差控制在4cm以下甚至9mm以下。
申请公布号 CN105044733A 申请公布日期 2015.11.11
申请号 CN201510523192.9 申请日期 2015.08.24
申请人 北京空间飞行器总体设计部 发明人 刘家兴;毕少筠;聂欣;崔小准;钟兆丰;武向军;郑晋军
分类号 G01S19/07(2010.01)I 主分类号 G01S19/07(2010.01)I
代理机构 北京理工大学专利中心 11120 代理人 李微微;仇蕾安
主权项 一种导航卫星TGD参数的标定方法,其特征在于,包括如下步骤:步骤1、测量待测导航卫星的频率为f<sub>1</sub>和f<sub>2</sub>的导航天线的相位中心在垂直方向的距离v,其中,垂直方向为导航卫星指向地心的方向;步骤2、选定n个地面站,其中n不小于3;各地面站所在的经度线和纬度线围成TGD参数标定观测区,其中,在选定的TGD参数标定时间段内,待测导航卫星的星下点轨迹落于所述观测区内;步骤3、确定导航信号频率为f<sub>1</sub>和f<sub>2</sub>的导航天线在各个地面站接收通道内的传输时延之差,标记为τ<sub>r1</sub>‑τ<sub>r2</sub>,即实现地面站接收通道IFB的标定;步骤4、在TGD参数标定时间段内,各个地面站天线分别持续跟踪待测导航卫星,分别通过地面站接收机获得多个观测时刻下导航信号频率为f<sub>1</sub>和f<sub>2</sub>的导航天线的伪码测距值ρ<sub>1</sub>和ρ<sub>2</sub>;步骤5、基于任意一个地面站的任意观测时刻建立观测方程:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mo>(</mo><msub><mi>&rho;</mi><mn>1</mn></msub><mo>-</mo><msub><mi>&rho;</mi><mn>2</mn></msub><mo>)</mo><mo>-</mo><mo>(</mo><msub><mi>&tau;</mi><mrow><mi>r</mi><mn>1</mn></mrow></msub><mo>-</mo><msub><mi>&tau;</mi><mrow><mi>r</mi><mn>2</mn></mrow></msub><mo>)</mo><mo>+</mo><mi>v</mi><mo>/</mo><mi>c</mi><mo>=</mo><mo>(</mo><msub><mi>&tau;</mi><mn>1</mn></msub><mo>-</mo><msub><mi>&tau;</mi><mn>2</mn></msub><mo>)</mo><mo>+</mo><mi>&kappa;</mi><mi>s</mi><mo>(</mo><msub><mi>c</mi><mn>1</mn></msub><mo>+</mo><msub><mi>c</mi><mn>2</mn></msub><mi>&phi;</mi><mo>+</mo><msub><mi>c</mi><mn>3</mn></msub><msub><mi>&lambda;</mi><mrow><mi>c</mi><mi>r</mi></mrow></msub><mo>+</mo><msub><mi>c</mi><mn>4</mn></msub><msup><mi>&phi;</mi><mn>2</mn></msup><mo>+</mo><msub><mi>c</mi><mn>5</mn></msub><msubsup><mi>&lambda;</mi><mrow><mi>c</mi><mi>r</mi></mrow><mn>2</mn></msubsup><mo>+</mo><msub><mi>c</mi><mn>6</mn></msub><msub><mi>&phi;&lambda;</mi><mrow><mi>c</mi><mi>r</mi></mrow></msub><mo>)</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000787198460000011.GIF" wi="1878" he="104" /></maths>其中,c为光速,τ<sub>1</sub>和τ<sub>2</sub>分别为导航信号频率为f<sub>1</sub>和f<sub>2</sub>的导航天线的发射通道时延,<img file="FDA0000787198460000012.GIF" wi="415" he="157" />s是以导航卫星对地面站天线仰角为变量的倾斜因子,c<sub>1</sub>,…,c<sub>6</sub>是电离层模型系数,φ是电离层穿刺点的地理纬度,λ<sub>cr</sub>是电离层穿刺点的太阳时角;先将同一个地面站在多个观测时刻对应的(1)式表示的观测方程形成方程组,求解τ<sub>1</sub>‑τ<sub>2</sub>;然后求得基于该地面站的TGD参数估计值:<img file="FDA0000787198460000013.GIF" wi="309" he="127" />其中<img file="FDA0000787198460000014.GIF" wi="191" he="157" />步骤6、根据星历计算在TGD参数标定时间段的中间时刻该导航卫星的星下点坐标,根据各个地面站地址以及地面站天线参数,求取各地面站天线的相位中心在地球海平面上投影点的坐标,分别计算所述导航卫星的星下点至各地面站天线投影点的距离;其中,用r<sub>i</sub>表示所述星下点至第i个地面站对应的投影点距离;步骤7、根据步骤5获得的各个地面站对应的TGD参数估计值以及步骤6获得的各地面站对应的星下点至投影点距离,求出最终的TGD参数估计结果:<img file="FDA0000787198460000021.GIF" wi="384" he="290" />其中,<img file="FDA0000787198460000022.GIF" wi="93" he="81" />表示第i个地面站对应的TGD参数估计值;然后用该最终估计结果T′<sub>GD</sub>对导航卫星进行TGD参数标定。
地址 100094 北京市海淀区友谊路104号