发明名称 一种基于升序排列优化的整周模糊度去相关方法
摘要 本发明公开了一种基于升序排列优化的整周模糊度去相关方法,包括以下几个步骤:步骤一,对整周模糊度浮点解的协方差矩阵进行Cholesky分解;步骤二,对分解后的矩阵进行行内积的升序排列;步骤三,对升序排列后的矩阵进行史密斯正交分解;步骤四,判定经史密斯正交分解后的矩阵G是否为单位矩阵,决定是否对矩阵进行反复迭代分解;步骤五,经过反复的迭代分解,最终得到去相关后的近似对角矩阵和反复迭代过程中产生的变换矩阵,完成对协方差矩阵的去相关。本发明采用升序方法对矩阵进行排序,使矩阵在去相关解算效果、解算成功率以及解算效率上都有显著提高。
申请公布号 CN103487819A 申请公布日期 2014.01.01
申请号 CN201310403119.9 申请日期 2013.09.06
申请人 哈尔滨工程大学 发明人 沈锋;陈潇;刘明凯;宋丽杰;李平敏;张桂贤;刘海峰
分类号 G01S19/44(2010.01)I 主分类号 G01S19/44(2010.01)I
代理机构 北京永创新实专利事务所 11121 代理人 赵文颖
主权项 1.一种基于升序排列优化的整周模糊度去相关方法,其特征在于,包括以下几个步骤:步骤一、根据GPS载波相位观测值,线性化后获取正定矩阵H,分解正定矩阵H;通过Cholesky分解,使正定矩阵H被分解为H=V<sup>T</sup>V   (6)其中,V为确定的满秩矩阵,v<sub>ij</sub>表示矩阵V的第i行第j列元素,则有:<img file="FDA0000378267520000011.GIF" wi="1205" he="261" />其中,V为下三角矩阵;步骤二、对矩阵V进行升序排列;采用排序矩阵S乘以Cholesky分解后的矩阵V得到:V<sub>k</sub>=SVS<sup>-1</sup>   (8)其中,S表示排序矩阵;排序方法为:若对矩阵进行下三角分解,则对下三角矩阵的行内积进行升序排列;若对矩阵进行上三角分解,则对上三角矩阵的行内积进行降序排列;则:B<sup>k</sup>=SB<sup>k+1</sup>S<sup>-1</sup>   (9)其中,B<sup>k</sup>表示排序后的三角矩阵,B<sup>k+1</sup>表示排序前的三角矩阵;当B<sup>k</sup>是下三角矩阵时,S是升序排列矩阵;当B<sup>k</sup>是上三角矩阵时,S是降序排列矩阵的逆,乘以排序矩阵后会使矩阵B<sup>k</sup>安升序或降序排列;步骤三、采用调整后的矩阵V<sub>k</sub>代替矩阵V进行史密斯正交化分解,V<sub>k</sub>=OG   (10)其中,G是单模矩阵,O矩阵为近似的正交化矩阵,整数的史密斯正交化分解过程如下:<maths num="0001"><![CDATA[<math><mrow><msub><mi>o</mi><mi>i</mi></msub><mo>=</mo><msub><mi>v</mi><mi>i</mi></msub><mo>-</mo><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mrow><mi>i</mi><mo>-</mo><mn>1</mn></mrow></munderover><msub><mi>g</mi><mi>ij</mi></msub><msub><mi>o</mi><mi>j</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>11</mn><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0002"><![CDATA[<math><mrow><msub><mi>g</mi><mi>ij</mi></msub><mo>=</mo><mfenced open='{' close=''><mtable><mtr><mtd><mo>[</mo><mo>&lt;</mo><msub><mi>v</mi><mi>i</mi></msub><mo>,</mo><msub><mi>o</mi><mi>j</mi></msub><mo>></mo><mo>/</mo><mo>&lt;</mo><msub><mi>o</mi><mi>j</mi></msub><mo>,</mo><msub><mi>o</mi><mi>j</mi></msub><mo>></mo><msub><mo>]</mo><mi>int</mi></msub></mtd><mtd><mi>j</mi><mo>&lt;</mo><mi>i</mi></mtd></mtr><mtr><mtd><mn>1</mn></mtd><mtd><mi>j</mi><mo>=</mo><mi>i</mi></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mi>j</mi><mo>></mo><mi>i</mi></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>12</mn><mo>)</mo></mrow></mrow></math>]]></maths>式中,o<sub>i</sub>、o<sub>j</sub>为矩阵O中的第i,j行元素,v<sub>i</sub>为矩阵V中的第i行元素,g<sub>ij</sub>为矩阵G中的第i行第j列元素;由式(4)和(6)可得:H=V<sup>T</sup>V=G<sup>T</sup>O<sup>T</sup>OG   (13)步骤四、判断,若矩阵G为单位阵则进入步骤五,若非单位阵则返回步骤一,对矩阵H<sub>1</sub>再次分解,直到分解后的矩阵G<sup>k</sup>为单位矩阵为止,其中k代表反复迭代的次数;步骤五、获取去相关矩阵H<sub>1k</sub>及变换矩阵Z<sub>k</sub>:H<sub>1k</sub>=Z<sub>k</sub>HZ<sub>k</sub><sup>T</sup>   (14)Z<sub>k</sub>=(G<sup>m</sup>)<sup>-1</sup>S<sup>m</sup>…(G<sup>k</sup>)<sup>-1</sup>S<sup>k</sup>…(G<sup>1</sup>)<sup>-1</sup>S<sup>1</sup>   (15)其中:S<sup>m</sup>表示经过m次反复迭代的排序矩阵,G<sup>m</sup>表示经过m次反复迭代的分解矩阵;根据去相关矩阵H<sub>1k</sub>以及变换矩阵Z<sub>k</sub>,进行模糊度搜索,实现GPS系统对高速运动的物体的定位。
地址 150001 黑龙江省哈尔滨市南岗区南通大街145号