主权项 |
一种基于双向通信协同定位辅助的单星定位方法,其特征在于:包括以下步骤:步骤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>Σ</mo><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msubsup><mi>λ</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>·</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>Σ</mo><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msubsup><mi>λ</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>Σ</mo><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msubsup><mi>λ</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>·</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>Σ</mo><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msubsup><mi>λ</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>Σ</mo><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msubsup><mi>λ</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>·</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>Σ</mo><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msubsup><mi>λ</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>Δ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>Δ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>Δ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三个方向上的分量。 |