主权项 |
一种实时北斗精密相对定位的方法,其特征在于,包括以下步骤:步骤一,整周模糊度参数初始化:利用基准站和移动站的同步北斗观测数据初始化载波相位整周模糊度参数;步骤二,周跳探测:对站星间的非同步双差载波相位进行周跳探测,对发生周跳的卫星进行标记;步骤三,伪距单点定位:利用基准站、移动站的北斗伪距观测信息,分别进行基准站和移动站伪距单点定位,获得基准站和移动站的绝对位置,获得基准站测量时的北斗卫星信号发送时刻的卫星坐标及钟差和移动站测量时的北斗卫星信号发送时刻的卫星坐标及钟差;步骤四,基准站与移动站的观测历元差判断:计算基准站的观测历元与移动站的观测历元之差,如果差值大于15s,则返回步骤一,否则进入步骤五;步骤五,周跳修复:采用附加模糊度参数方法修复发生周跳的卫星的整周模糊度,即利用未发生周跳的卫星的整周模糊度计算发生周跳的卫星的整周模糊度;步骤六,求解基准站与移动站的相对位置:利用基准站和移动站的非差非同步载波相位观测数据和卫星坐标以及钟差,建立非同步载波相位双差相对定位观测方程,计算基准站与移动站的相对位置,具体过程为:(S1)建立载波相位双差观测方程组:利用基准站和移动站的非差载波相位观测值及广播星历信息,建立双差载波相位相对定位观测方程;设基准站A在观测时刻t<sub>0</sub>,测得i号卫星的非差载波相位观测值为<img file="FDA0000743264550000021.GIF" wi="163" he="89" />和j号卫星的非差载波相位观测值为<img file="FDA0000743264550000022.GIF" wi="187" he="88" />移动站B在观测时刻t<sub>1</sub>,测得i号卫星的非差载波相位观测值为<img file="FDA0000743264550000023.GIF" wi="150" he="86" />和j号卫星的非差载波相位观测值为<img file="FDA0000743264550000024.GIF" wi="180" he="83" />i,j为整数;定义:<maths num="0001" id="cmaths0001"><math><![CDATA[<mfenced open='{' close=''><mtable><mtr><mtd><msubsup><mrow><mo>(</mo><mo>·</mo><mo>)</mo></mrow><mi>AB</mi><mi>ij</mi></msubsup><mrow><mo>(</mo><msub><mi>t</mi><mn>0</mn></msub><mo>,</mo><msub><mi>t</mi><mn>1</mn></msub><mo>)</mo></mrow><mo>=</mo><mrow><mo>(</mo><msubsup><mrow><mo>(</mo><mo>·</mo><mo>)</mo></mrow><mi>B</mi><mi>j</mi></msubsup><mrow><mo>(</mo><msub><mi>t</mi><mn>1</mn></msub><mo>)</mo></mrow><mo>-</mo><msubsup><mrow><mo>(</mo><mo>·</mo><mo>)</mo></mrow><mi>A</mi><mi>j</mi></msubsup><mrow><mo>(</mo><msub><mi>t</mi><mn>0</mn></msub><mo>)</mo></mrow><mo>)</mo></mrow><mo>-</mo><mrow><mo>(</mo><msubsup><mrow><mo>(</mo><mo>·</mo><mo>)</mo></mrow><mi>B</mi><mi>i</mi></msubsup><mrow><mo>(</mo><msub><mi>t</mi><mn>1</mn></msub><mo>)</mo></mrow><mo>-</mo><msubsup><mrow><mo>(</mo><mo>·</mo><mo>)</mo></mrow><mi>A</mi><mi>i</mi></msubsup><mrow><mo>(</mo><msub><mi>t</mi><mn>0</mn></msub><mo>)</mo></mrow><mo>)</mo></mrow></mtd></mtr><mtr><mtd><msubsup><mrow><mo>(</mo><mo>·</mo><mo>)</mo></mrow><mi>AB</mi><mi>ij</mi></msubsup><mrow><mo>(</mo><msub><mi>T</mi><mn>0</mn></msub><mo>,</mo><msub><mi>T</mi><mn>1</mn></msub><mo>)</mo></mrow><mo>=</mo><mrow><mo>(</mo><msubsup><mrow><mo>(</mo><mo>·</mo><mo>)</mo></mrow><mi>B</mi><mi>j</mi></msubsup><mrow><mo>(</mo><msubsup><mi>T</mi><mn>1</mn><mi>j</mi></msubsup><mo>)</mo></mrow><mo>-</mo><msubsup><mrow><mo>(</mo><mo>·</mo><mo>)</mo></mrow><mi>A</mi><mi>j</mi></msubsup><mrow><mo>(</mo><msubsup><mi>T</mi><mn>0</mn><mi>j</mi></msubsup><mo>)</mo></mrow><mo>)</mo></mrow><mo>-</mo><mrow><mo>(</mo><msubsup><mrow><mo>(</mo><mo>·</mo><mo>)</mo></mrow><mi>B</mi><mi>i</mi></msubsup><mrow><mo>(</mo><msubsup><mi>T</mi><mn>1</mn><mi>i</mi></msubsup><mo>)</mo></mrow><mo>-</mo><msubsup><mrow><mo>(</mo><mo>·</mo><mo>)</mo></mrow><mi>A</mi><mi>i</mi></msubsup><mrow><mo>(</mo><msubsup><mi>T</mi><mn>0</mn><mi>i</mi></msubsup><mo>)</mo></mrow><mo>)</mo></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000743264550000025.GIF" wi="1270" he="231" /></maths>其中,上标i和j代表卫星号,i号卫星为参考星,下标A和B代表基准站和移动站;T代表卫星信号发送时刻,如<img file="FDA0000743264550000026.GIF" wi="91" he="76" />表示在观测时刻t<sub>0</sub>时,对应i号卫星的信号发送时刻,单位为秒;则非同步载波相位双差观测方程为:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msubsup><mi>Φ</mi><mi>AB</mi><mi>ij</mi></msubsup><mrow><mo>(</mo><msub><mi>t</mi><mn>0</mn></msub><mo>,</mo><msub><mi>t</mi><mn>1</mn></msub><mo>)</mo></mrow><mo>=</mo><msubsup><mi>ρ</mi><mi>AB</mi><mi>ij</mi></msubsup><mrow><mo>(</mo><msub><mi>T</mi><mn>0</mn></msub><mo>,</mo><msub><mi>T</mi><mn>1</mn></msub><mo>)</mo></mrow><mo>+</mo><mi>λ</mi><msubsup><mi>N</mi><mi>AB</mi><mi>ij</mi></msubsup><mrow><mo>(</mo><msub><mi>t</mi><mn>0</mn></msub><mo>,</mo><msub><mi>t</mi><mn>1</mn></msub><mo>)</mo></mrow><mo>-</mo><mi>c</mi><msubsup><mi>δt</mi><mi>AB</mi><mi>ij</mi></msubsup><mrow><mo>(</mo><msub><mi>T</mi><mn>0</mn></msub><mo>,</mo><msub><mi>T</mi><mn>1</mn></msub><mo>)</mo></mrow><mo>+</mo><msubsup><mi>ϵ</mi><mi>AB</mi><mi>ij</mi></msubsup><mrow><mo>(</mo><msub><mi>t</mi><mn>0</mn></msub><mo>,</mo><msub><mi>t</mi><mn>1</mn></msub><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000743264550000027.GIF" wi="1544" he="85" /></maths>其中,<img file="FDA0000743264550000028.GIF" wi="230" he="85" />为不同时刻的双差载波相位观测值,单位为米;<img file="FDA0000743264550000029.GIF" wi="234" he="89" />为不同时刻的双差几何距离,单位为米;<img file="FDA00007432645500000210.GIF" wi="229" he="86" />为不同时刻的双差整周模糊度,单位为周;λ为某个频点的对应波长,单位为米;<img file="FDA00007432645500000211.GIF" wi="246" he="84" />为不同时刻的双差卫星钟差,单位为秒;c为光速,单位为米/秒;<img file="FDA00007432645500000212.GIF" wi="206" he="85" />为不同时刻的双差噪声,单位为米;式(1)中,基站与移动站的相对位置参数隐含于双差几何距离量<img file="FDA00007432645500000213.GIF" wi="234" he="84" />中,其具体为表达式为:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msubsup><mi>ρ</mi><mi>AB</mi><mi>ij</mi></msubsup><mrow><mo>(</mo><msub><mi>T</mi><mn>0</mn></msub><mo>,</mo><msub><mi>T</mi><mn>1</mn></msub><mo>)</mo></mrow><mo>=</mo><mrow><mo>(</mo><mo>|</mo><mo>|</mo><msub><mi>r</mi><mi>B</mi></msub><mrow><mo>(</mo><msub><mi>t</mi><mn>1</mn></msub><mo>)</mo></mrow><mo>-</mo><msup><mi>r</mi><mi>j</mi></msup><mrow><mo>(</mo><msubsup><mi>T</mi><mn>1</mn><mi>j</mi></msubsup><mo>)</mo></mrow><mo>|</mo><mo>|</mo><mo>-</mo><mo>|</mo><mo>|</mo><msub><mi>r</mi><mi>A</mi></msub><mrow><mo>(</mo><msub><mi>t</mi><mn>0</mn></msub><mo>)</mo></mrow><mo>-</mo><msup><mi>r</mi><mi>j</mi></msup><mrow><mo>(</mo><msubsup><mi>T</mi><mn>0</mn><mi>j</mi></msubsup><mo>)</mo></mrow><mo>|</mo><mo>|</mo><mo>)</mo></mrow><mo>-</mo><mrow><mo>(</mo><mo>|</mo><mo>|</mo><msub><mi>r</mi><mi>B</mi></msub><mrow><mo>(</mo><msub><mi>t</mi><mn>1</mn></msub><mo>)</mo></mrow><mo>-</mo><msup><mi>r</mi><mi>i</mi></msup><mrow><mo>(</mo><msubsup><mi>T</mi><mn>1</mn><mi>i</mi></msubsup><mo>)</mo></mrow><mo>|</mo><mo>|</mo><mo>-</mo><mo>|</mo><mo>|</mo><msub><mi>r</mi><mi>A</mi></msub><mrow><mo>(</mo><msub><mi>t</mi><mn>0</mn></msub><mo>)</mo></mrow><mo>-</mo><msup><mi>r</mi><mi>i</mi></msup><mrow><mo>(</mo><msubsup><mi>T</mi><mn>0</mn><mi>i</mi></msubsup><mo>)</mo></mrow><mo>|</mo><mo>|</mo><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA00007432645500000214.GIF" wi="1809" he="102" /></maths>式(2)中,<img file="FDA00007432645500000215.GIF" wi="159" he="100" />和<img file="FDA00007432645500000216.GIF" wi="152" he="97" />为j号卫星在信号发送时刻T<sub>1</sub>和T<sub>0</sub>的三维位置坐标,单位为米;<img file="FDA0000743264550000031.GIF" wi="148" he="95" />和<img file="FDA0000743264550000032.GIF" wi="144" he="97" />为i号卫星在信号发送时刻T<sub>1</sub>和T<sub>0</sub>的三维位置坐标,单位为米;||.||表示卫星到测站之间的几何距离,单位为米;r<sub>A</sub>(t<sub>0</sub>)为基准站在观测时刻t<sub>0</sub>的三维位置坐标;r<sub>B</sub>(t<sub>1</sub>)为移动站在观测时刻t<sub>1</sub>的三维位置坐标,单位为米;设基准站和移动站共同观测了m颗卫星,则有建立(m‑1)个双差载波相位观测方程:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mfenced open='{' close=''><mtable><mtr><mtd><msubsup><mi>Φ</mi><mi>AB</mi><mrow><mi>i</mi><mn>1</mn></mrow></msubsup><mrow><mo>(</mo><msub><mi>t</mi><mn>0</mn></msub><mo>,</mo><msub><mi>t</mi><mn>1</mn></msub><mo>)</mo></mrow><mo>=</mo><msubsup><mi>ρ</mi><mi>AB</mi><mrow><mi>i</mi><mn>1</mn></mrow></msubsup><mrow><mo>(</mo><msub><mi>T</mi><mn>0</mn></msub><mo>,</mo><msub><mi>T</mi><mn>1</mn></msub><mo>)</mo></mrow><mo>+</mo><mi>λ</mi><msubsup><mi>N</mi><mi>AB</mi><mrow><mi>i</mi><mn>1</mn></mrow></msubsup><mrow><mo>(</mo><msub><mi>t</mi><mn>0</mn></msub><mo>,</mo><msub><mi>t</mi><mn>1</mn></msub><mo>)</mo></mrow><mo>-</mo><mi>cδ</mi><msubsup><mi>t</mi><mi>AB</mi><mrow><mi>i</mi><mn>1</mn></mrow></msubsup><mrow><mo>(</mo><msub><mi>T</mi><mn>0</mn></msub><mo>,</mo><msub><mi>T</mi><mn>1</mn></msub><mo>)</mo></mrow><mo>+</mo><msubsup><mi>ϵ</mi><mi>AB</mi><mrow><mi>i</mi><mn>1</mn></mrow></msubsup><mrow><mo>(</mo><msub><mi>t</mi><mn>0</mn></msub><mo>,</mo><msub><mi>t</mi><mn>1</mn></msub><mo>)</mo></mrow></mtd></mtr><mtr><mtd><msubsup><mi>Φ</mi><mi>AB</mi><mrow><mi>i</mi><mn>2</mn></mrow></msubsup><mrow><mo>(</mo><msub><mi>t</mi><mn>0</mn></msub><mo>,</mo><msub><mi>t</mi><mn>1</mn></msub><mo>)</mo></mrow><mo>=</mo><msubsup><mi>ρ</mi><mi>AB</mi><mrow><mi>i</mi><mn>2</mn></mrow></msubsup><mrow><mo>(</mo><msub><mi>T</mi><mn>0</mn></msub><mo>,</mo><msub><mi>T</mi><mn>1</mn></msub><mo>)</mo></mrow><mo>+</mo><mi>λ</mi><msubsup><mi>N</mi><mi>AB</mi><mrow><mi>i</mi><mn>2</mn></mrow></msubsup><mrow><mo>(</mo><msub><mi>t</mi><mn>0</mn></msub><mo>,</mo><msub><mi>t</mi><mn>1</mn></msub><mo>)</mo></mrow><mo>-</mo><mi>cδ</mi><msubsup><mi>t</mi><mi>AB</mi><mrow><mi>i</mi><mn>2</mn></mrow></msubsup><mrow><mo>(</mo><msub><mi>T</mi><mn>0</mn></msub><mo>,</mo><msub><mi>T</mi><mn>1</mn></msub><mo>)</mo></mrow><mo>+</mo><msubsup><mi>ϵ</mi><mi>AB</mi><mrow><mi>i</mi><mn>2</mn></mrow></msubsup><mrow><mo>(</mo><msub><mi>t</mi><mn>0</mn></msub><mo>,</mo><msub><mi>t</mi><mn>1</mn></msub><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd></mtr><mtr><mtd><msubsup><mi>Φ</mi><mi>AB</mi><mi>ik</mi></msubsup><mrow><mo>(</mo><msub><mi>t</mi><mn>0</mn></msub><mo>,</mo><msub><mi>t</mi><mn>1</mn></msub><mo>)</mo></mrow><mo>=</mo><msubsup><mi>ρ</mi><mi>AB</mi><mi>ik</mi></msubsup><mrow><mo>(</mo><msub><mi>T</mi><mn>0</mn></msub><mo>,</mo><msub><mi>T</mi><mn>1</mn></msub><mo>)</mo></mrow><mo>+</mo><mi>λ</mi><msubsup><mi>N</mi><mi>AB</mi><mi>ik</mi></msubsup><mrow><mo>(</mo><msub><mi>t</mi><mn>0</mn></msub><mo>,</mo><msub><mi>t</mi><mn>1</mn></msub><mo>)</mo></mrow><mo>-</mo><mi>cδ</mi><msubsup><mi>t</mi><mi>AB</mi><mi>ik</mi></msubsup><mrow><mo>(</mo><msub><mi>T</mi><mn>0</mn></msub><mo>,</mo><msub><mi>T</mi><mn>1</mn></msub><mo>)</mo></mrow><mo>+</mo><msubsup><mi>ϵ</mi><mi>AB</mi><mi>ik</mi></msubsup><mrow><mo>(</mo><msub><mi>t</mi><mn>0</mn></msub><mo>,</mo><msub><mi>t</mi><mn>1</mn></msub><mo>)</mo></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000743264550000033.GIF" wi="1565" he="343" /></maths>式中,k表示卫星号,k=1,2,…,j,…m‑1;k≠i;k,m为整数;(S2)相对定位观测方程组线性化:根据泰勒展开原理,对非线性相对定位观测方程组进行线性化;在式(2)中,卫星的坐标以实时播发的广播星历参数计算;在基准站位置线性泰勒展开至一阶项,式(2)变为<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><mfenced open='{' close=''><mtable><mtr><mtd><msubsup><mi>ρ</mi><mi>AB</mi><mi>ij</mi></msubsup><mrow><mo>(</mo><msub><mi>T</mi><mn>0</mn></msub><mo>,</mo><msub><mi>T</mi><mn>1</mn></msub><mo>)</mo></mrow><mo>=</mo><msubsup><mi>ρ</mi><mn>0</mn><mi>j</mi></msubsup><mo>-</mo><mrow><mo>(</mo><msubsup><mi>u</mi><mi>B</mi><mi>j</mi></msubsup><mo>-</mo><msubsup><mi>u</mi><mi>B</mi><mi>i</mi></msubsup><mo>)</mo></mrow><mo>-</mo><msub><mi>r</mi><mi>AB</mi></msub></mtd></mtr><mtr><mtd><msubsup><mi>ρ</mi><mn>0</mn><mi>j</mi></msubsup><mo>=</mo><mrow><mo>(</mo><mo>|</mo><mo>|</mo><msub><mi>r</mi><mi>A</mi></msub><mrow><mo>(</mo><msub><mi>t</mi><mn>0</mn></msub><mo>)</mo></mrow><mo>-</mo><msup><mi>r</mi><mi>j</mi></msup><mrow><mo>(</mo><msubsup><mi>T</mi><mn>1</mn><mi>j</mi></msubsup><mo>)</mo></mrow><mo>|</mo><mo>|</mo><mo>-</mo><mo>|</mo><mo>|</mo><msub><mi>r</mi><mi>A</mi></msub><mrow><mo>(</mo><msub><mi>t</mi><mn>0</mn></msub><mo>)</mo></mrow><mo>-</mo><msup><mi>r</mi><mi>j</mi></msup><mrow><mo>(</mo><msubsup><mi>T</mi><mn>0</mn><mi>j</mi></msubsup><mo>)</mo></mrow><mo>|</mo><mo>|</mo><mo>)</mo></mrow><mo>-</mo><mrow><mo>(</mo><mo>|</mo><mo>|</mo><msub><mi>r</mi><mi>A</mi></msub><mrow><mo>(</mo><msub><mi>t</mi><mn>0</mn></msub><mo>)</mo></mrow><mo>-</mo><msup><mi>r</mi><mi>i</mi></msup><mrow><mo>(</mo><msubsup><mi>T</mi><mn>1</mn><mi>i</mi></msubsup><mo>)</mo></mrow><mo>|</mo><mo>|</mo><mo>-</mo><mo>|</mo><mo>|</mo><msub><mi>r</mi><mi>A</mi></msub><mrow><mo>(</mo><msub><mi>t</mi><mn>0</mn></msub><mo>)</mo></mrow><mo>-</mo><msup><mi>r</mi><mi>i</mi></msup><mrow><mo>(</mo><msubsup><mi>T</mi><mn>0</mn><mi>i</mi></msubsup><mo>)</mo></mrow><mo>|</mo><mo>|</mo><mo>)</mo></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000743264550000034.GIF" wi="1821" he="219" /></maths>其中,u代表单位矢量,<img file="FDA0000743264550000035.GIF" wi="56" he="80" />表示移动站B到j号卫星的单位矢量;<img file="FDA0000743264550000036.GIF" wi="64" he="77" />表示移动站B到i号卫星的单位矢量;<img file="FDA0000743264550000037.GIF" wi="70" he="76" />表示j号卫星的双差几何距离泰勒展开的零阶项;r<sub>AB</sub>为基准站到移动站的基线矢量;将式(4)代入式(3),则线性化的观测方程组为:<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><mfenced open='{' close=''><mtable><mtr><mtd><msubsup><mi>Φ</mi><mi>AB</mi><mrow><mi>i</mi><mn>1</mn></mrow></msubsup><mrow><mo>(</mo><msub><mi>t</mi><mn>0</mn></msub><mo>,</mo><msub><mi>t</mi><mn>1</mn></msub><mo>)</mo></mrow><mo>-</mo><msubsup><mi>ρ</mi><mn>0</mn><mn>1</mn></msubsup><mo>-</mo><mi>λ</mi><msubsup><mi>N</mi><mi>AB</mi><mrow><mi>i</mi><mn>1</mn></mrow></msubsup><mrow><mo>(</mo><msub><mi>t</mi><mn>0</mn></msub><mo>,</mo><msub><mi>t</mi><mn>1</mn></msub><mo>)</mo></mrow><mo>+</mo><mi>cδ</mi><msubsup><mi>t</mi><mi>AB</mi><mrow><mi>i</mi><mn>1</mn></mrow></msubsup><mrow><mo>(</mo><msub><mi>T</mi><mn>0</mn></msub><mo>,</mo><msub><mi>T</mi><mn>1</mn></msub><mo>)</mo></mrow><mo>=</mo><mo>-</mo><mrow><mo>(</mo><msubsup><mi>u</mi><mi>B</mi><mn>1</mn></msubsup><mo>-</mo><msubsup><mi>u</mi><mi>B</mi><mi>i</mi></msubsup><mo>)</mo></mrow><msub><mi>r</mi><mi>AB</mi></msub><mo>+</mo><msubsup><mi>ϵ</mi><mi>AB</mi><mrow><mi>i</mi><mn>1</mn></mrow></msubsup><mrow><mo>(</mo><msub><mi>t</mi><mn>0</mn></msub><mo>,</mo><msub><mi>t</mi><mn>1</mn></msub><mo>)</mo></mrow></mtd></mtr><mtr><mtd><msubsup><mi>Φ</mi><mi>AB</mi><mrow><mi>i</mi><mn>2</mn></mrow></msubsup><mrow><mo>(</mo><msub><mi>t</mi><mn>0</mn></msub><mo>,</mo><msub><mi>t</mi><mn>1</mn></msub><mo>)</mo></mrow><mo>-</mo><msubsup><mi>ρ</mi><mn>0</mn><mn>2</mn></msubsup><mo>-</mo><mi>λ</mi><msubsup><mi>N</mi><mi>AB</mi><mrow><mi>i</mi><mn>2</mn></mrow></msubsup><mrow><mo>(</mo><msub><mi>t</mi><mn>0</mn></msub><mo>,</mo><msub><mi>t</mi><mn>1</mn></msub><mo>)</mo></mrow><mo>+</mo><mi>cδ</mi><msubsup><mi>t</mi><mi>AB</mi><mrow><mi>i</mi><mn>2</mn></mrow></msubsup><mrow><mo>(</mo><msub><mi>T</mi><mn>0</mn></msub><mo>,</mo><msub><mi>T</mi><mn>1</mn></msub><mo>)</mo></mrow><mo>=</mo><mo>-</mo><mrow><mo>(</mo><msubsup><mi>u</mi><mi>B</mi><mn>2</mn></msubsup><mo>-</mo><msubsup><mi>u</mi><mi>B</mi><mi>i</mi></msubsup><mo>)</mo></mrow><msub><mi>r</mi><mi>AB</mi></msub><mo>+</mo><msubsup><mi>ϵ</mi><mi>AB</mi><mrow><mi>i</mi><mn>2</mn></mrow></msubsup><mrow><mo>(</mo><msub><mi>t</mi><mn>0</mn></msub><mo>,</mo><msub><mi>t</mi><mn>1</mn></msub><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd></mtr><mtr><mtd><msubsup><mi>Φ</mi><mi>AB</mi><mi>ik</mi></msubsup><mrow><mo>(</mo><msub><mi>t</mi><mn>0</mn></msub><mo>,</mo><msub><mi>t</mi><mn>1</mn></msub><mo>)</mo></mrow><mo>-</mo><msubsup><mi>ρ</mi><mn>0</mn><mi>k</mi></msubsup><mo>-</mo><mi>λ</mi><msubsup><mi>N</mi><mi>AB</mi><mi>ik</mi></msubsup><mrow><mo>(</mo><msub><mi>t</mi><mn>0</mn></msub><mo>,</mo><msub><mi>t</mi><mn>1</mn></msub><mo>)</mo></mrow><mo>+</mo><mi>cδ</mi><msubsup><mi>t</mi><mi>AB</mi><mi>ik</mi></msubsup><mrow><mo>(</mo><msub><mi>T</mi><mn>0</mn></msub><mo>,</mo><msub><mi>T</mi><mn>1</mn></msub><mo>)</mo></mrow><mo>=</mo><mo>-</mo><mrow><mo>(</mo><msubsup><mi>u</mi><mi>B</mi><mi>k</mi></msubsup><mo>-</mo><msubsup><mi>u</mi><mi>B</mi><mi>i</mi></msubsup><mo>)</mo></mrow><msub><mi>r</mi><mi>AB</mi></msub><mo>+</mo><msubsup><mi>ϵ</mi><mi>AB</mi><mi>ik</mi></msubsup><mrow><mo>(</mo><msub><mi>t</mi><mn>0</mn></msub><mo>,</mo><msub><mi>t</mi><mn>1</mn></msub><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="FDA0000743264550000041.GIF" wi="1656" he="372" /></maths>式中,k表示卫星号,k=1,2,…,j,…m‑1;k≠i;k,m为整数;式(5)中,卫星钟差双差<img file="FDA0000743264550000042.GIF" wi="249" he="85" />根据实时广播星历参数计算;(S3)线性观测方程组求解:根据最小二乘估计,即可求得基准站到移动站的基线矢量;记:<maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><mi>Y</mi><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msubsup><mi>Φ</mi><mi>AB</mi><mrow><mi>i</mi><mn>1</mn></mrow></msubsup><mrow><mo>(</mo><msub><mi>t</mi><mn>0</mn></msub><mo>,</mo><msub><mi>t</mi><mn>1</mn></msub><mo>)</mo></mrow><mo>-</mo><msubsup><mi>ρ</mi><mn>0</mn><mn>1</mn></msubsup><mo>-</mo><mi>λ</mi><msubsup><mi>N</mi><mi>AB</mi><mrow><mi>i</mi><mn>1</mn></mrow></msubsup><mrow><mo>(</mo><msub><mi>t</mi><mn>0</mn></msub><mo>,</mo><msub><mi>t</mi><mn>1</mn></msub><mo>)</mo></mrow><mo>+</mo><mi>cδ</mi><msubsup><mi>t</mi><mi>AB</mi><mrow><mi>i</mi><mn>1</mn></mrow></msubsup><mrow><mo>(</mo><msub><mi>T</mi><mn>0</mn></msub><mo>,</mo><msub><mi>T</mi><mn>1</mn></msub><mo>)</mo></mrow></mtd></mtr><mtr><mtd><msubsup><mi>Φ</mi><mi>AB</mi><mrow><mi>i</mi><mn>2</mn></mrow></msubsup><mrow><mo>(</mo><msub><mi>t</mi><mn>0</mn></msub><mo>,</mo><msub><mi>t</mi><mn>1</mn></msub><mo>)</mo></mrow><mo>-</mo><msubsup><mi>ρ</mi><mn>0</mn><mn>2</mn></msubsup><mo>-</mo><mi>λ</mi><msubsup><mi>N</mi><mi>AB</mi><mrow><mi>i</mi><mn>2</mn></mrow></msubsup><mrow><mo>(</mo><msub><mi>t</mi><mn>0</mn></msub><mo>,</mo><msub><mi>t</mi><mn>1</mn></msub><mo>)</mo></mrow><mo>+</mo><mi>cδ</mi><msubsup><mi>t</mi><mi>AB</mi><mrow><mi>i</mi><mn>2</mn></mrow></msubsup><mrow><mo>(</mo><msub><mi>T</mi><mn>0</mn></msub><mo>,</mo><msub><mi>T</mi><mn>1</mn></msub><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd></mtr><mtr><mtd><msubsup><mi>Φ</mi><mi>AB</mi><mi>ik</mi></msubsup><mrow><mo>(</mo><msub><mi>t</mi><mn>0</mn></msub><mo>,</mo><msub><mi>t</mi><mn>1</mn></msub><mo>)</mo></mrow><mo>-</mo><msubsup><mi>ρ</mi><mn>0</mn><mi>k</mi></msubsup><mo>-</mo><mi>λ</mi><msubsup><mi>N</mi><mi>AB</mi><mi>ik</mi></msubsup><mrow><mo>(</mo><msub><mi>t</mi><mn>0</mn></msub><mo>,</mo><msub><mi>t</mi><mn>1</mn></msub><mo>)</mo></mrow><mo>+</mo><mi>cδ</mi><msubsup><mi>t</mi><mi>AB</mi><mi>ik</mi></msubsup><mrow><mo>(</mo><msub><mi>T</mi><mn>0</mn></msub><mo>,</mo><msub><mi>T</mi><mn>1</mn></msub><mo>)</mo></mrow></mtd></mtr></mtable></mfenced><mo>,</mo><mi>Q</mi><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><mo>-</mo><mrow><mo>(</mo><msubsup><mi>u</mi><mi>B</mi><mn>1</mn></msubsup><mo>-</mo><msubsup><mi>u</mi><mi>B</mi><mi>i</mi></msubsup><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mo>-</mo><mrow><mo>(</mo><msubsup><mi>u</mi><mi>B</mi><mn>2</mn></msubsup><mo>-</mo><msubsup><mi>u</mi><mi>B</mi><mi>i</mi></msubsup><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd></mtr><mtr><mtd><mo>-</mo><mrow><mo>(</mo><msubsup><mi>u</mi><mi>B</mi><mi>k</mi></msubsup><mo>-</mo><msubsup><mi>u</mi><mi>B</mi><mi>i</mi></msubsup><mo>)</mo></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000743264550000043.GIF" wi="1598" he="372" /></maths>式中,k表示卫星号,k=1,2,…,j,…m‑1;k≠i;k,m为整数;根据最小二乘估计原理,求解式(5),得到:r<sub>AB</sub>=(Q<sup>T</sup>Q)<sup>‑1</sup>Q<sup>T</sup>Y (7)r<sub>AB</sub>即可确定基准站与移动站之间的相对位置。 |