发明名称 数字卫星电视授时系统测轨方法
摘要 一种数字卫星电视授时系统卫星测轨方法,由地面站时间精确同步、上行地面南京站A时间戳生成与时间戳插播、上行地面南京站A时间戳生成与时间戳插播步骤组成。本发明采用现有数字电视系统达到卫星测轨的目的,通过标记电视信号中特殊标识位的收发时间,计算信号传播时延,采用一发多收的测轨模式,根据多个时延测量值,列出关于卫星坐标的方程组,计算得到卫星坐标。本发明与GNSS卫星导航授时系统相比,卫星电视信号具有信号强、抗干扰能力强等优点,避免了GNSS信号所面临的恶意干扰和攻击问题。
申请公布号 CN103983943A 申请公布日期 2014.08.13
申请号 CN201410242671.9 申请日期 2014.05.31
申请人 中国科学院国家授时中心 发明人 向渝;华宇;许林生;董道鹏;武晓亮;王善和;刘晓龙
分类号 G01S1/00(2006.01)I 主分类号 G01S1/00(2006.01)I
代理机构 西安永生专利代理有限责任公司 61201 代理人 申忠才
主权项 一种数字卫星电视授时系统测轨方法,其特征在于它是由下述步骤组成:(1)地面站时间精确同步以上行地面南京站为A、地面接收西安站记为B、地面接收乌鲁木齐站为C、地面接收昆明站为D、地面接收长春站为E、地面接收南京站为F,以每天凌晨00:00为起始时间,采用卫星双向时间频率传递方法,每隔1分钟测A与B、C、D、E、F之间的钟差一次,每天共测100次,每天B、C、D、E、F站分别与A站之间的钟差公式为:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>&Delta;T</mi><mi>Ai</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>a</mi><mi>i</mi></msub><mo>+</mo><msub><mi>b</mi><mi>i</mi></msub><mi>t</mi><mo>+</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><msub><mi>c</mi><mi>i</mi></msub><msup><mi>t</mi><mn>2</mn></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000514699990000011.GIF" wi="1284" he="128" /></maths>式中t为距00:00时的分钟数,a<sub>i</sub>为i站原子钟与A站原子钟的初始相位差,b<sub>i</sub>为i站原子钟的频偏,c<sub>i</sub>为i站原子钟的漂移率,i依次为B、C、D、E、F地面接收站;(2)上行地面南京站A时间戳生成与时间戳插播在上行地面南京站A电视上行系统中调制器的串行异步接口取出速率为270兆比特每秒的码流信号,对信号完成同步,生成纳秒级时间戳与时间戳插播;(3)电视信号接收与数据处理地面接收西安站B、地面接收乌鲁木齐站C、地面接收昆明站D、地面接收长春站E、地面接收南京站F接收卫星下行电视信号,将所接收的电视信号经解调通过网络输出到数据处理服务器,数据处理服务器对输入的数据进行数据处理,得到卫星的位置,具体如下:1)地面接收西安站B、地面接收乌鲁木齐站C、地面接收昆明站D、地面接收长春站E、地面接收南京站F通过卫星天线接收来自卫星的下行信号,经解调、输出270兆比特每秒的异步串行码流信号,进行信号同步过程,进行时间戳生成过程与上行地面南京站A发送过程完全相同;将传输流包中同步头的第七位到第八位的上升沿所对应的精确时刻记作t<sub>r</sub>,与对应的节目时钟参考值编码生成接收时刻的时间戳,存入先进先出的存储器1中;2)电视系统中的解码器通过包标识符为0x120以识别传输流包,解调0x120对应的传输流包数据得到发送时刻生成的时间戳。将发送时刻的时间戳中的节目时钟参考值作为关键字在存储器1中寻找相同的节目时钟参考,进行数据处理;数据处理步骤如下:将对应的时刻值相减,得到大环时延估计值;Δt′<sub>i</sub>=t<sub>r</sub>‑t<sub>s</sub>  (3)3)由公式(2)计算得到t<sub>r</sub>时刻钟差ΔT<sub>Ai</sub>(t<sub>r</sub>)修正,即可得到精确的大环时延值Δt<sub>i</sub>=Δt′<sub>i</sub>+ΔT<sub>Ai</sub>(t<sub>r</sub>)  (4)4)地面接收西安站B、地面接收乌鲁木齐站C、地面接收昆明站D、地面接收长春站E、地面接收南京站F将各自坐标和所测时延实时汇总至上行站,通过求解方程组(5)得到卫星的坐标:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mfenced open='{' close=''><mtable><mtr><mtd><mn>2</mn><mrow><mo>(</mo><msub><mi>X</mi><mi>B</mi></msub><mo>-</mo><msub><mi>X</mi><mi>A</mi></msub><mo>)</mo></mrow><mi>X</mi><mo>+</mo><mn>2</mn><mrow><mo>(</mo><msub><mi>Y</mi><mi>B</mi></msub><mo>-</mo><msub><mi>Y</mi><mi>A</mi></msub><mo>)</mo></mrow><mi>Y</mi><mo>+</mo><mn>2</mn><mrow><mo>(</mo><msub><mi>Z</mi><mi>B</mi></msub><mo>-</mo><msub><mi>Z</mi><mi>A</mi></msub><mo>)</mo></mrow><mi>Z</mi><mo>-</mo><mn>2</mn><mrow><mo>(</mo><msub><mi>&rho;</mi><mi>B</mi></msub><mo>-</mo><msub><mi>&rho;</mi><mi>A</mi></msub><mo>)</mo></mrow><mi>&lambda;</mi><mo>=</mo><msubsup><mi>L</mi><mi>B</mi><mn>2</mn></msubsup><mo>-</mo><msubsup><mi>L</mi><mi>A</mi><mn>2</mn></msubsup><mo>-</mo><mrow><mo>(</mo><msubsup><mi>&rho;</mi><mi>B</mi><mn>2</mn></msubsup><mo>-</mo><msubsup><mi>&rho;</mi><mi>A</mi><mn>2</mn></msubsup><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mn>2</mn><mrow><mo>(</mo><msub><mi>X</mi><mi>C</mi></msub><mo>-</mo><msub><mi>X</mi><mi>A</mi></msub><mo>)</mo></mrow><mi>X</mi><mo>+</mo><mn>2</mn><mrow><mo>(</mo><msub><mi>Y</mi><mi>C</mi></msub><mo>-</mo><msub><mi>Y</mi><mi>A</mi></msub><mo>)</mo></mrow><mi>Y</mi><mo>+</mo><mn>2</mn><mrow><mo>(</mo><msub><mi>Z</mi><mi>C</mi></msub><mo>-</mo><msub><mi>Z</mi><mi>A</mi></msub><mo>)</mo></mrow><mi>Z</mi><mo>-</mo><mn>2</mn><mrow><mo>(</mo><msub><mi>&rho;</mi><mi>C</mi></msub><mo>-</mo><msub><mi>&rho;</mi><mi>A</mi></msub><mo>)</mo></mrow><mi>&lambda;</mi><mo>=</mo><msubsup><mi>L</mi><mi>C</mi><mn>2</mn></msubsup><mo>-</mo><msubsup><mi>L</mi><mi>A</mi><mn>2</mn></msubsup><mo>-</mo><mrow><mo>(</mo><msubsup><mi>&rho;</mi><mi>C</mi><mn>2</mn></msubsup><mo>-</mo><msubsup><mi>&rho;</mi><mi>A</mi><mn>2</mn></msubsup><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mn>2</mn><mrow><mo>(</mo><msub><mi>X</mi><mi>D</mi></msub><mo>-</mo><msub><mi>X</mi><mi>A</mi></msub><mo>)</mo></mrow><mi>X</mi><mo>+</mo><mn>2</mn><mrow><mo>(</mo><msub><mi>Y</mi><mi>D</mi></msub><mo>-</mo><msub><mi>Y</mi><mi>A</mi></msub><mo>)</mo></mrow><mi>Y</mi><mo>+</mo><mn>2</mn><mrow><mo>(</mo><msub><mi>Z</mi><mi>D</mi></msub><mo>-</mo><msub><mi>Z</mi><mi>A</mi></msub><mo>)</mo></mrow><mi>Z</mi><mo>-</mo><mn>2</mn><mrow><mo>(</mo><msub><mi>&rho;</mi><mi>D</mi></msub><mo>-</mo><msub><mi>&rho;</mi><mi>A</mi></msub><mo>)</mo></mrow><mi>&lambda;</mi><mo>=</mo><msubsup><mi>L</mi><mi>D</mi><mn>2</mn></msubsup><mo>-</mo><msubsup><mi>L</mi><mi>A</mi><mn>2</mn></msubsup><mo>-</mo><mrow><mo>(</mo><msubsup><mi>&rho;</mi><mi>D</mi><mn>2</mn></msubsup><mo>-</mo><msubsup><mi>&rho;</mi><mi>A</mi><mn>2</mn></msubsup><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mn>2</mn><mrow><mo>(</mo><msub><mi>X</mi><mi>E</mi></msub><mo>-</mo><msub><mi>X</mi><mi>A</mi></msub><mo>)</mo></mrow><mi>X</mi><mo>+</mo><mn>2</mn><mrow><mo>(</mo><msub><mi>Y</mi><mi>E</mi></msub><mo>-</mo><msub><mi>Y</mi><mi>A</mi></msub><mo>)</mo></mrow><mi>Y</mi><mo>+</mo><mn>2</mn><mrow><mo>(</mo><msub><mi>Z</mi><mi>E</mi></msub><mo>-</mo><msub><mi>Z</mi><mi>A</mi></msub><mo>)</mo></mrow><mi>Z</mi><mo>-</mo><mn>2</mn><mrow><mo>(</mo><msub><mi>&rho;</mi><mi>E</mi></msub><mo>-</mo><msub><mi>&rho;</mi><mi>A</mi></msub><mo>)</mo></mrow><mi>&lambda;</mi><mo>=</mo><msubsup><mi>L</mi><mi>E</mi><mn>2</mn></msubsup><mo>-</mo><msubsup><mi>L</mi><mi>A</mi><mn>2</mn></msubsup><mo>-</mo><mrow><mo>(</mo><msubsup><mi>&rho;</mi><mi>E</mi><mn>2</mn></msubsup><mo>-</mo><msubsup><mi>&rho;</mi><mi>A</mi><mn>2</mn></msubsup><mo>)</mo></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000514699990000021.GIF" wi="1849" he="335" /></maths>式中,ρ<sub>i</sub>为c(Δt<sub>i</sub>‑τ′<sub>i</sub>),λ为cτ<sub>clk</sub>,<img file="FDA0000514699990000022.GIF" wi="60" he="79" />为X<sub>i</sub><sup>2</sup>+Y<sub>i</sub><sup>2</sup>+Z<sub>i</sub><sup>2</sup>,i为B、C、D、E、F,时延修正τ′<sub>i</sub>通过电离层和对流层数学模型得到,作为已知参数,i站坐标(X<sub>i</sub>,Y<sub>i</sub>,Z<sub>i</sub>)已知,τ<sub>clk</sub>是上行时延,c为光速。
地址 710600 陕西省西安市临潼区书院东路3号