发明名称 基于双向通信协同定位辅助的单星定位方法
摘要 本发明提出一种基于双向通信协同定位辅助的单星定位方法,分为单星定位模块和协同定位模块两个部分,单星定位模块位于卫星和地面探测器,主要完成地面探测器的单星定位解算以及地面探测器之间距离的求解,协同定位模块位于地面探测器,主要完成地面探测器之间的协同定位。将地面探测器所得到的自身定位坐标发送到定位卫星,定位卫星根据所需定位的地面探测器自身定位坐标按照一定的距离限制条件确定协同定位的地面探测器,并通过卫星计算得到的地面探测器和协同地面探测器之间的距离信息建立权重系数,并通过该系数融合各个协同地面探测器的单星定位结果从而减小卫星定轨误差对地面探测器定位误差的影响,并提高定位精度。
申请公布号 CN106019346A 申请公布日期 2016.10.12
申请号 CN201610517247.X 申请日期 2016.07.04
申请人 西北工业大学 发明人 唐成凯;廉保旺;岳哲;张怡;张玲玲
分类号 G01S19/46(2010.01)I 主分类号 G01S19/46(2010.01)I
代理机构 西北工业大学专利中心 61204 代理人 陈星
主权项 一种基于双向通信协同定位辅助的单星定位方法,其特征在于:包括以下步骤:步骤1:采用以下过程获得每个地面探测器的初始坐标估计值:对于第i个地面探测器,在三个不同时刻接收单星发射信号,并获得三个不同时刻下单星的坐标,利用多普勒测频仪获得第i个地面探测器和单星之间,在第一时刻与第二时刻的伪距差,以及在第二时刻与第三时刻的伪距差;根据得到的两个伪距差以及第i个地面探测器的自身高度信息,结合三个不同时刻下单星的坐标,通过最小二乘法得到第i个地面探测器的初始坐标估计值<img file="FDA00010403477300000113.GIF" wi="203" he="64" />并将该初始坐标估计值发送给单星;步骤2:采用以下过程迭代更新每个地面探测器的坐标估计值:对于第i个地面探测器,其坐标估计值更新公式为:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msubsup><mover><mi>x</mi><mo>^</mo></mover><mi>i</mi><mi>k</mi></msubsup><mo>=</mo><mfrac><mn>1</mn><mrow><mn>1</mn><mo>+</mo><munderover><mo>&Sigma;</mo><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msubsup><mi>&lambda;</mi><mrow><mi>j</mi><mo>,</mo><mi>i</mi></mrow><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msubsup></mrow></mfrac><mo>&CenterDot;</mo><mrow><mo>(</mo><msubsup><mover><mi>x</mi><mo>^</mo></mover><mi>i</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msubsup><mo>+</mo><munderover><mo>&Sigma;</mo><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msubsup><mi>&lambda;</mi><mrow><mi>j</mi><mo>,</mo><mi>i</mi></mrow><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msubsup><msubsup><mover><mi>x</mi><mo>^</mo></mover><mrow><mi>j</mi><mo>,</mo><mi>i</mi></mrow><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msubsup><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001040347730000011.GIF" wi="694" he="191" /></maths><maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msubsup><mover><mi>y</mi><mo>^</mo></mover><mi>i</mi><mi>k</mi></msubsup><mo>=</mo><mfrac><mn>1</mn><mrow><mn>1</mn><mo>+</mo><munderover><mo>&Sigma;</mo><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msubsup><mi>&lambda;</mi><mrow><mi>j</mi><mo>,</mo><mi>i</mi></mrow><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msubsup></mrow></mfrac><mo>&CenterDot;</mo><mrow><mo>(</mo><msubsup><mover><mi>y</mi><mo>^</mo></mover><mi>i</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msubsup><mo>+</mo><munderover><mo>&Sigma;</mo><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msubsup><mi>&lambda;</mi><mrow><mi>j</mi><mo>,</mo><mi>i</mi></mrow><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msubsup><msubsup><mover><mi>y</mi><mo>^</mo></mover><mrow><mi>j</mi><mo>,</mo><mi>i</mi></mrow><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msubsup><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001040347730000012.GIF" wi="702" he="191" /></maths><maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msubsup><mover><mi>z</mi><mo>^</mo></mover><mi>i</mi><mi>k</mi></msubsup><mo>=</mo><mfrac><mn>1</mn><mrow><mn>1</mn><mo>+</mo><munderover><mo>&Sigma;</mo><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msubsup><mi>&lambda;</mi><mrow><mi>j</mi><mo>,</mo><mi>i</mi></mrow><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msubsup></mrow></mfrac><mo>&CenterDot;</mo><mrow><mo>(</mo><msubsup><mover><mi>z</mi><mo>^</mo></mover><mi>i</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msubsup><mo>+</mo><munderover><mo>&Sigma;</mo><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msubsup><mi>&lambda;</mi><mrow><mi>j</mi><mo>,</mo><mi>i</mi></mrow><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msubsup><msubsup><mover><mi>z</mi><mo>^</mo></mover><mrow><mi>j</mi><mo>,</mo><mi>i</mi></mrow><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msubsup><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001040347730000013.GIF" wi="686" he="196" /></maths>其中<img file="FDA0001040347730000014.GIF" wi="266" he="71" />为第i个地面探测器在第k‑1次迭代后的坐标估计值,<img file="FDA0001040347730000015.GIF" wi="180" he="70" />为第i个地面探测器在第k次迭代后的坐标估计值,n为与第i个地面探测器直线距离不大于设定阈值的地面探测器个数,<img file="FDA0001040347730000016.GIF" wi="78" he="70" />为第k‑1次迭代后n个地面探测器中的第j个地面探测器到所述第i个地面探测器的融合系数,融合系数与对应两个地面探测器间的直线距离负相关;<img file="FDA0001040347730000017.GIF" wi="267" he="71" />为n个地面探测器中的第j个地面探测器对所述第i个地面探测器的估计坐标:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msubsup><mover><mi>x</mi><mo>^</mo></mover><mrow><mi>j</mi><mo>,</mo><mi>i</mi></mrow><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msubsup><mo>=</mo><msubsup><mover><mi>x</mi><mo>^</mo></mover><mi>j</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msubsup><mo>+</mo><msubsup><mi>&Delta;x</mi><mrow><mi>j</mi><mo>,</mo><mi>i</mi></mrow><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msubsup></mrow>]]></math><img file="FDA0001040347730000018.GIF" wi="349" he="71" /></maths><maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msubsup><mover><mi>y</mi><mo>^</mo></mover><mrow><mi>j</mi><mo>,</mo><mi>i</mi></mrow><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msubsup><mo>=</mo><msubsup><mover><mi>y</mi><mo>^</mo></mover><mi>j</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msubsup><mo>+</mo><msubsup><mi>&Delta;y</mi><mrow><mi>j</mi><mo>,</mo><mi>i</mi></mrow><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msubsup></mrow>]]></math><img file="FDA0001040347730000019.GIF" wi="357" he="71" /></maths><maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><msubsup><mover><mi>z</mi><mo>^</mo></mover><mrow><mi>j</mi><mo>,</mo><mi>i</mi></mrow><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msubsup><mo>=</mo><msubsup><mover><mi>z</mi><mo>^</mo></mover><mi>j</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msubsup><mo>+</mo><msubsup><mi>&Delta;z</mi><mrow><mi>j</mi><mo>,</mo><mi>i</mi></mrow><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msubsup></mrow>]]></math><img file="FDA00010403477300000110.GIF" wi="341" he="71" /></maths><img file="FDA00010403477300000111.GIF" wi="262" he="79" />为n个地面探测器中的第j个地面探测器自身在第k‑1次迭代后的坐标估计值,<img file="FDA00010403477300000112.GIF" wi="342" he="77" />为第k‑1次迭代后n个地面探测器中的第j个地面探测器到所述第i个地面探测器的直线距离在x,y,z三个方向上的分量。
地址 710072 陕西省西安市友谊西路127号
您可能感兴趣的专利