CN105549050B - 一种基于模糊置信度滤波的北斗变形监测定位方法 - Google Patents
一种基于模糊置信度滤波的北斗变形监测定位方法 Download PDFInfo
- Publication number
- CN105549050B CN105549050B CN201510896845.8A CN201510896845A CN105549050B CN 105549050 B CN105549050 B CN 105549050B CN 201510896845 A CN201510896845 A CN 201510896845A CN 105549050 B CN105549050 B CN 105549050B
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- mtd
- mtr
- msup
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/38—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
- G01S19/39—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/42—Determining position
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
Abstract
本发明公开了一种基于模糊置信度滤波的北斗变形监测定位方法，其特征是按如下步骤进行：首先将监测时段划分成多个子时段，分别计算出各子时段的定位结果；然后根据各个子时段的卫星分布、误差源、卫星历元数，采用模糊集合理论确定各个定位结果的模糊置信度；最后利用模糊置信度加权滤波出整个监测时段的最终定位结果。本发明细化了定位计算过程，充分考虑了影响定位精度的各项因素，使最终定位结果达到更高的定位精度，特别是在监测时段内天气变化的情况下定位精度优势显著，具有广阔的应用前景。
Description
技术领域
本发明属于卫星导航定位领域，具体地说是一种基于模糊置信度滤波的北斗变形监测定位方法。
背景技术
北斗变形监测是指利用北斗导航卫星“静态相对定位技术”获取监测点高精度坐标(毫米级)，从而分析判断大坝、桥梁、地基等的位移或沉降，具有重要的应用价值。
目前，采用导航卫星(GPS、GLONASS、北斗)技术进行大地变形监测时，均采用“静态相对定位方法”。该方法认为监测站静止不动，通过采集长时间的、大量的卫星观测数据，应用最小二乘平差原理求解基准站和监测站之间的基线矢量。因为拥有大量的卫星观测数据，所以可以通过迭代计算逐步求精，解算出高精度的基线矢量，从而得到“毫米级”的定位结果，进而可以对缓慢变形的对象(大坝、桥梁、地基等)进行长期监测。
如图1所示，传统的静态相对定位方法是利用一个较长时段的所有卫星观测数据，以后处理的方式计算出定位结果L。由于在这个时段内卫星空间分布会变化，天气也会变化，因此接收到卫星数据质量是变化的，有些子时段卫星数据质量好，有些子时段卫星数据质量一般或者差。传统方法并不区分这种数据质量的差异性，笼统地利用所有卫星数据计算定位结果，从而导致了定位结果误差增大，难以满足高精度应用需求。
发明内容
本发明是为克服上述现有技术的不足之处，提出了一种基于模糊置信度滤波的北斗变形监测定位方法，以期能提高定位结果的精度，特别是在监测时段内天气变化的情况下定位精度优势显著。
本发明为达到上述发明目的，采用如下技术方案：
本发明一种基于模糊置信度滤波的北斗变形监测定位方法，是应用于由北斗导航卫星、基准站和流动站组成的监测环境中；所述基准站接收所述北斗导航卫星发送的基准站卫星历元数据并发送给流动站；所述流动站接收所述北斗导航卫星发送的流动站卫星历元数据和所述基准站发送的基准站卫星历元数据并进行差分处理，获得差分定位数据，记为X；其特点是，所述北斗变形监测定位方法按如下步骤进行：
步骤1、对所述流动站的差分定位数据X按照时间段进行分割，获得N个子时间段的差分定位数据，记为X＝{X1,X2,…,Xn,…,XN}；Xn表示流动站的差分定位数据X中第n个子时间段的差分定位数据；1≤n≤N；
步骤2、对所述N个子时间段的差分定位数据X分别采用静态相对定位算法进行计算，获得N个定位结果，记为L＝{L1,L2,…,Ln,…,LN}；Ln表示流动站第n个子时间段的差分定位数据xn的定位结果；
步骤3、离线建立置信度的参数矩阵Sm×3；
流动站的定位结果的可信程度主要受三个参数的影响，它们可计算如下：
步骤3.1、利用式(1)获得第n个子时间段的定位结果Ln的几何精度因子平均值
式(1)中，表示流动站第n个子时间段的定位结果Ln的水平分量精度因子平均值；表示流动站第n个子时间段的定位结果Ln的垂直分量精度因子平均值；表示流动站第n个子时间段的定位结果Ln的钟差精度因子平均值；
步骤3.2、利用式(2)获得第n个子时间段的定位结果Ln的高空大气误差平均值
式(2)中，表示流动站第n个子时间段的定位结果Ln的电离层误差平均值；表示流动站第n个子时间段的定位结果Ln的对流层误差平均值；
步骤3.3、统计流动站第n个子时间段内的卫星历元数量Nn；
步骤3.4、定义置信度T∈{1,2,…,j,…,m}表示定位结果的可信程度，其中m是正整数；定义置信度的参数向量为S＝{S1,S2,…,Sj,…,Sm}；Sj表示置信度T＝j所对应的参数向量；并有表示几何精度因子参数值，通过式(3)计算得到；表示高空大气误差参数值，通过式(4)计算得到；表示卫星历元数量参数值，通过式(5)计算得到。
式(3)中，表示几何精度因子的最佳值；表示置信度T＝j的几何精度因子；式(4)中，表示高空大气误差的最佳值；表示置信度T＝j的高空大气误差；式(5)中，N*表示卫星历元数量的最佳值；Nj表示置信度T＝j的卫星历元数量；
步骤3.5、建立置信度的参数矩阵
步骤4、建立定位结果的评价矩阵EN×3；
步骤4.1、定义定位结果评价向量为E＝{E1,E2,…,En,…,EN}；En表示流动站第n个子时间段定位结果的评价向量；并有表示流动站第n个子时间段的几何精度因子评价值，通过式(6)计算得到；表示流动站第n个子时间段的高空大气误差评价值，通过式(7)计算得到；表示流动站第n个子时间段的卫星历元数量评价值，通过式(8)计算得到；
步骤4.2、建立定位结果的评价矩阵
步骤5、利用式(9)建立隶属度函数
式(9)中，
步骤6、建立关系矩阵
步骤7、利用式(10)获得阈值λ；
步骤8、利用式(11)对所述关系矩阵R进行处理，获得布尔矩阵
式(11)中，表示流动站第n个子时间段差分定位数据xn的定位结果Ln的置信度T＝j；
步骤9、利用式(12)对流动站的N个定位结果进行加权滤波处理，获得最终定位结果L*；
式(12)中，T(k)表示流动站第k个子时间段的定位结果Lk的置信度。
与现有技术相比，本发明的有益效果体现在：
1、本发明将一个较长监测时段划分成多个子时段，利用各个子时段的卫星观测数据分别计算出各自子时段的定位结果，细化了定位计算过程，避免了传统定位方法的“笼统”的计算理念，有助于利用各个子时段卫星数据质量的差异性，以提高定位精度。
2、本发明引入置信度的概念，依据各个子时段的卫星数据质量参数(卫星分布几何精度因子、高空大气误差、卫星历元数量)的变化划分不同的置信度等级，充分考虑了卫星空间分布变化，天气变化等各项因素对定位结果的影响，提高了定位结果的可靠性；同时，采用了模糊集合理论确定各个子时段定位结果的置信度，并以加权滤波的方式求解出最终定位结果，从而提高了定位结果的精度，具有广阔的应用前景。
附图说明
图1为传统静态相对定位方法；
图2为本发明的方法原理；
图3为本发明由数据质量确定定位结果置信度的内在关联图；
图4为本发明基于模糊置信度滤波的北斗变形监测定位方法流程图；
图5为本发明方法与传统静态相对定位方法的实验结果比较。
具体实施方式
一种基于模糊置信度滤波的北斗变形监测定位方法，是应用于由北斗导航卫星、基准站和流动站组成的监测环境中；基准站接收北斗导航卫星发送的基准站卫星历元数据并发送给流动站；流动站接收北斗导航卫星发送的流动站卫星历元数据和基准站发送的基准站卫星历元数据，二者进行差分处理获得差分定位数据，记为X；
如图4所示，北斗变形监测定位方法按如下步骤进行：
步骤1、对流动站的差分定位数据X按照时间段进行分割，获得N个子时间段的差分定位数据，例如将24小时观测时段划分成12个子时段，每个子时段2小时；记为X＝{X1,X2,…,Xn,…,XN}；Xn表示流动站的差分定位数据X中第n个子时间段的差分定位数据；1≤n≤N；
步骤2、对N个子时间段的差分定位数据X分别采用静态相对定位算法进行计算，获得N个定位结果，记为L＝{L1,L2,…,Ln,…,LN}；Ln表示流动站第n个子时间段的差分定位数据Xn的定位结果；静态相对定位算法的内容见可参考“王成.载波相位差分动态定位的方法研究[D].西安:长安大学,2010.”
步骤3、离线建立置信度的参数矩阵Sm×3；
流动站的定位结果的可信程度主要受三个参数的影响，它们可计算如下：
步骤3.1、利用式(1)获得第n个子时间段的定位结果Ln的几何精度因子平均值它反映了监测站和卫星之间的空间几何关系。值与监测站到卫星单位矢量端点所形成的多面体体积成反比，多面体体积越大的卫星组合的值越小。在观测误差一定时，值越小，定位精度越高。
式(1)中，表示流动站第n个子时间段的定位结果Ln的水平分量精度因子平均值；表示流动站第n个子时间段的定位结果Ln的垂直分量精度因子平均值；表示流动站第n个子时间段的定位结果Ln的钟差精度因子平均值；它们可通过指令从接收机获取。在流动站第n个子时间段内，作为该子时段定位结果可信程度的评价参数之一。
大量实验结果表明：当时，定位算法解算的结果精度高；否则，定位算法解算的结果误差较大；因此，本实施例中，令最佳值
步骤3.2、利用式(2)获得第n时间段的定位结果Ln的高空大气误差平均值它反映了监测站上空大气变化对定位结果的影响。
式(2)中，表示流动站第n个子时间段的定位结果Ln的电离层误差平均值；表示流动站第个子时间段的定位结果Ln的对流层误差平均值；它们可通过指令从接收机获取。在流动站第n个子时间段内，作为该子时段定位结果可信程度的评价参数之二。
大量实验结果表明：当时，定位算法解算的结果精度高；否则，定位算法解算的结果误差较大。因此，本实施例中，令最佳值
步骤3.3、统计流动站第n个子时间段内的卫星历元数量Nn；
卫星历元数量是否充足直接决定了定位结果的精度。对于静态相对定位，卫星历元数量越大，定位结果精度越高；反之，定位结果精度越低。在流动站第n个子时间段内，我们统计卫星历元数量Nn，作为该子时段定位结果可信程度的评价参数之三；如图3所示。
在正常情况下，接收机工作频率为1Hz，每2小时可以收到7200个卫星历元，可以保证定位结果的高精度。但是在卫星信号被遮挡、接收机供电不稳定甚至停电的情况下，接收机收到的卫星历元数量Nn<3600时，静态相对定位的结果误差较大。令最佳值N*＝7200。
步骤3.4、定义置信度T∈{1,2,…,j,…,m}表示定位结果的可信程度，其中m是正整数；T取值越大表示可信程度越高；定义置信度的参数向量为S＝{S1,S2,…,Sj,…,Sm}；Sj表示置信度T＝j所对应的参数向量；并有表示几何精度因子参数值，通过式(3)计算得到；表示高空大气误差参数值，通过式(4)计算得到；表示卫星历元数量参数值，通过式(5)计算得到：
式(3)中，表示几何精度因子的最佳值；表示置信度T＝j的几何精度因子；式(4)中，表示高空大气误差的最佳值；表示置信度T＝j的高空大气误差；式(5)中，N*表示卫星历元数量的最佳值；Nj表示置信度T＝j的卫星历元数量；
步骤3.5、建立置信度的参数矩阵
置信度的参数矩阵列出了不同置信度分别的3个参数值。体现了不同置信度的差异。
在实际应用中，可设置每个置信度的3个参数值如表1所示。采用最佳值对其进行数量级统一后的参数向量Sj如表2所示。
表1.置信度的参数值
表2.置信度的参数向量Sj
那么，一种有效的置信度的参数矩阵为：
步骤4、建立定位结果的评价矩阵EN×3；
步骤4.1、定义定位结果评价向量为E＝{E1,E2,…,En,…,EN}；En表示流动站第n个子时间段定位结果的评价向量；并有表示流动站第n个子时间段的几何精度因子评价值，通过式(6)计算得到；表示流动站第n个子时间段的高空大气误差评价值，通过式(7)计算得到；表示流动站第n个子时间段的卫星历元数量评价值，通过式(8)计算得到；
步骤4.2、建立定位结果的评价矩阵
定位结果的评价矩阵列出了N个子时段的N个定位结果分别的3个评价值；反映了各个定位结果在3个评价值上的异同。
步骤5、利用式(9)建立隶属度函数
式(9)中，把定位结果评价向量Ei对置信度参数向量Sj的隶属情况进行了综合。
步骤6、建立关系矩阵
矩阵R根据式(9)计算所得，其行代表N个定位结果评价向量，列代表m个置信度参数向量，这样就使得定位结果与置信度建立起关系，从而为定位结果置信度的确定提供了计算依据。
步骤7、利用式(10)获得阈值λ：
阈值λ的大小直接影响定位结果的置信度的确定。λ值越大，定位结果会被确定为越少的置信度，从而专注率越高；反之，定位结果会被确定为越多的置信度。考虑到一个定位结果至少要有一个置信度，则
步骤8、利用式(11)对所述关系矩阵R进行处理，获得布尔矩阵
式(11)中，表示流动站第n个子时间段差分定位数据xn的定位结果Ln的置信度T＝j；这里依据模糊聚类中的λ-截矩阵理论将模糊关系矩阵转化为布尔矩阵。
由于本专利方法是基于模糊集合理论确定的各子时段定位结果的置信度，因此该置信度称为“模糊置信度”。
步骤9、利用式(12)对流动站的N个定位结果进行加权滤波处理，获得最终定位结果L*如图2所示；
式(12)中，T(k)表示流动站第k个子时间段的定位结果Lk的模糊置信度。
实验验证：
2015年10月8日至11月8日，在合肥某处高填方变形监测项目中同时应用了传统静态相对定位方法和本专利方法，均进行了30天的变形监测，它们的沉降监测结果如图5所示。可见，传统静态相对定位方法的监测结果起伏较大，误差明显，它主要受到了天气变化和卫星历元数量等因素的影响；本专利方法的监测结果更准确，更稳定，更真实地反映了高填方的持续沉降过程。
Claims (1)
1.一种基于模糊置信度滤波的北斗变形监测定位方法，是应用于由北斗导航卫星、基准站和流动站组成的监测环境中；所述基准站接收所述北斗导航卫星发送的基准站卫星历元数据并发送给流动站；所述流动站接收所述北斗导航卫星发送的流动站卫星历元数据和所述基准站发送的基准站卫星历元数据并进行差分处理，获得差分定位数据，记为X；其特征是，所述北斗变形监测定位方法按如下步骤进行：
步骤1、对所述流动站的差分定位数据X按照时间段进行分割，获得N个子时间段的差分定位数据，记为X＝{X1,X2,…,Xn,…,XN}；Xn表示流动站的差分定位数据X中第n个子时间段的差分定位数据；1≤n≤N；
步骤2、对所述N个子时间段的差分定位数据X分别采用静态相对定位算法进行计算，获得N个定位结果，记为L＝{L1,L2,…,Ln,…,LN}；Ln表示流动站第n个子时间段的差分定位数据Xn的定位结果；
步骤3、离线建立置信度的参数矩阵Sm×3；
流动站的定位结果的可信程度主要受三个参数的影响，它们可计算如下：
步骤3.1、利用式(1)获得第n个子时间段的定位结果Ln的几何精度因子平均值
<mrow>
<mover>
<mrow>
<msub>
<mi>GDOP</mi>
<mi>n</mi>
</msub>
</mrow>
<mo>&OverBar;</mo>
</mover>
<mo>=</mo>
<msqrt>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<mover>
<mrow>
<msub>
<mi>HDOP</mi>
<mi>n</mi>
</msub>
</mrow>
<mo>&OverBar;</mo>
</mover>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<mover>
<mrow>
<msub>
<mi>VDOP</mi>
<mi>n</mi>
</msub>
</mrow>
<mo>&OverBar;</mo>
</mover>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<mover>
<mrow>
<msub>
<mi>TDOP</mi>
<mi>n</mi>
</msub>
</mrow>
<mo>&OverBar;</mo>
</mover>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
式(1)中，表示流动站第n个子时间段的定位结果Ln的水平分量精度因子平均值；表示流动站第n个子时间段的定位结果Ln的垂直分量精度因子平均值；表示流动站第n个子时间段的定位结果Ln的钟差精度因子平均值；
步骤3.2、利用式(2)获得第n个子时间段的定位结果Ln的高空大气误差平均值
<mrow>
<mover>
<mrow>
<msub>
<mi>&Delta;&tau;</mi>
<mi>n</mi>
</msub>
</mrow>
<mo>&OverBar;</mo>
</mover>
<mo>=</mo>
<msqrt>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<mover>
<mrow>
<msub>
<mi>&Delta;</mi>
<mi>I</mi>
</msub>
<msub>
<mi>&tau;</mi>
<mi>n</mi>
</msub>
</mrow>
<mo>&OverBar;</mo>
</mover>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<mover>
<mrow>
<msub>
<mi>&Delta;</mi>
<mi>T</mi>
</msub>
<msub>
<mi>&tau;</mi>
<mi>n</mi>
</msub>
</mrow>
<mo>&OverBar;</mo>
</mover>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
式(2)中，表示流动站第n个子时间段的定位结果Ln的电离层误差平均值；表示流动站第n个子时间段的定位结果Ln的对流层误差平均值；
步骤3.3、统计流动站第n个子时间段内的卫星历元数量Nn；
步骤3.4、定义置信度T∈{1,2,…,j,…,m}表示定位结果的可信程度，m是正整数；定义置信度的参数向量为S＝{S1,S2,…,Sj,…,Sm}；Sj表示置信度T＝j所对应的参数向量；并有 表示几何精度因子参数值，通过式(3)计算得到；表示高空大气误差参数值，通过式(4)计算得到；表示卫星历元数量参数值，通过式(5)计算得到：
<mrow>
<msubsup>
<mi>s</mi>
<mi>j</mi>
<mn>1</mn>
</msubsup>
<mo>=</mo>
<msup>
<mover>
<mrow>
<mi>G</mi>
<mi>D</mi>
<mi>O</mi>
<mi>P</mi>
</mrow>
<mo>&OverBar;</mo>
</mover>
<mo>*</mo>
</msup>
<mo>/</mo>
<msup>
<mover>
<mrow>
<mi>G</mi>
<mi>D</mi>
<mi>O</mi>
<mi>P</mi>
</mrow>
<mo>&OverBar;</mo>
</mover>
<mi>j</mi>
</msup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msubsup>
<mi>s</mi>
<mi>j</mi>
<mn>2</mn>
</msubsup>
<mo>=</mo>
<msup>
<mover>
<mrow>
<mi>&Delta;</mi>
<mi>&tau;</mi>
</mrow>
<mo>&OverBar;</mo>
</mover>
<mo>*</mo>
</msup>
<mo>/</mo>
<msup>
<mover>
<mrow>
<mi>&Delta;</mi>
<mi>&tau;</mi>
</mrow>
<mo>&OverBar;</mo>
</mover>
<mi>j</mi>
</msup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msubsup>
<mi>s</mi>
<mi>j</mi>
<mn>3</mn>
</msubsup>
<mo>=</mo>
<msup>
<mi>N</mi>
<mi>j</mi>
</msup>
<mo>/</mo>
<msup>
<mi>N</mi>
<mo>*</mo>
</msup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>5</mn>
<mo>)</mo>
</mrow>
</mrow>
式(3)中，表示几何精度因子的最佳值；表示置信度T＝j的几何精度因子；式(4)中，表示高空大气误差的最佳值；表示置信度T＝j的高空大气误差；式(5)中，N*表示卫星历元数量的最佳值；Nj表示置信度T＝j的卫星历元数量；
步骤3.5、建立置信度的参数矩阵
步骤4、建立定位结果的评价矩阵EN×3；
步骤4.1、定义定位结果评价向量为E＝{E1,E2,…,En,…,EN}；En表示流动站第n个子时间段定位结果的评价向量；并有 表示流动站第n个子时间段的几何精度因子评价值，通过式(6)计算得到；表示流动站第n个子时间段的高空大气误差评价值，通过式(7)计算得到；表示流动站第n个子时间段的卫星历元数量评价值，通过式(8)计算得到；
<mrow>
<msubsup>
<mi>e</mi>
<mi>n</mi>
<mn>1</mn>
</msubsup>
<mo>=</mo>
<mn>1</mn>
<mo>-</mo>
<mfrac>
<mrow>
<mo>|</mo>
<mover>
<mrow>
<msub>
<mi>GDOP</mi>
<mi>n</mi>
</msub>
</mrow>
<mo>&OverBar;</mo>
</mover>
<mo>-</mo>
<msup>
<mover>
<mrow>
<mi>G</mi>
<mi>D</mi>
<mi>O</mi>
<mi>P</mi>
</mrow>
<mo>&OverBar;</mo>
</mover>
<mo>*</mo>
</msup>
<mo>|</mo>
</mrow>
<msup>
<mover>
<mrow>
<mi>G</mi>
<mi>D</mi>
<mi>O</mi>
<mi>P</mi>
</mrow>
<mo>&OverBar;</mo>
</mover>
<mo>*</mo>
</msup>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>6</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msubsup>
<mi>e</mi>
<mi>n</mi>
<mn>2</mn>
</msubsup>
<mo>=</mo>
<mn>1</mn>
<mo>-</mo>
<mfrac>
<mrow>
<mo>|</mo>
<mover>
<mrow>
<msub>
<mi>&Delta;&tau;</mi>
<mi>n</mi>
</msub>
</mrow>
<mo>&OverBar;</mo>
</mover>
<mo>-</mo>
<msup>
<mover>
<mrow>
<mi>&Delta;</mi>
<mi>&tau;</mi>
</mrow>
<mo>&OverBar;</mo>
</mover>
<mo>*</mo>
</msup>
<mo>|</mo>
</mrow>
<msup>
<mover>
<mrow>
<mi>&Delta;</mi>
<mi>&tau;</mi>
</mrow>
<mo>&OverBar;</mo>
</mover>
<mo>*</mo>
</msup>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>7</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msubsup>
<mi>e</mi>
<mi>n</mi>
<mn>3</mn>
</msubsup>
<mo>=</mo>
<mn>1</mn>
<mo>-</mo>
<mfrac>
<mrow>
<mo>|</mo>
<msub>
<mi>N</mi>
<mi>n</mi>
</msub>
<mo>-</mo>
<msup>
<mi>N</mi>
<mo>*</mo>
</msup>
<mo>|</mo>
</mrow>
<msup>
<mi>N</mi>
<mo>*</mo>
</msup>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>8</mn>
<mo>)</mo>
</mrow>
</mrow>
步骤4.2、建立定位结果的评价矩阵
步骤5、利用式(9)建立隶属度函数
<mrow>
<msub>
<mi>&mu;</mi>
<msub>
<mi>S</mi>
<mi>j</mi>
</msub>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>E</mi>
<mi>n</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mn>3</mn>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>v</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mn>3</mn>
</munderover>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<msup>
<mrow>
<mo>(</mo>
<msubsup>
<mi>e</mi>
<mi>n</mi>
<mi>v</mi>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>s</mi>
<mi>j</mi>
<mi>v</mi>
</msubsup>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</msup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>9</mn>
<mo>)</mo>
</mrow>
</mrow>
式(9)中，
步骤6、建立关系矩阵
步骤7、利用式(10)获得阈值λ；
<mrow>
<mi>&lambda;</mi>
<mo>&le;</mo>
<munder>
<mrow>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
</mrow>
<mi>j</mi>
</munder>
<mo>{</mo>
<munder>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
<mi>n</mi>
</munder>
<mo>{</mo>
<msub>
<mi>&mu;</mi>
<msub>
<mi>S</mi>
<mi>j</mi>
</msub>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>E</mi>
<mi>n</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>}</mo>
<mo>}</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>10</mn>
<mo>)</mo>
</mrow>
</mrow>
步骤8、利用式(11)对所述关系矩阵R进行处理，获得布尔矩阵
<mrow>
<msup>
<mi>R</mi>
<mo>&prime;</mo>
</msup>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msup>
<msub>
<mi>&mu;</mi>
<msub>
<mi>S</mi>
<mn>1</mn>
</msub>
</msub>
<mo>&prime;</mo>
</msup>
<mrow>
<mo>(</mo>
<msub>
<mi>E</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<msup>
<msub>
<mi>&mu;</mi>
<msub>
<mi>S</mi>
<mn>2</mn>
</msub>
</msub>
<mo>&prime;</mo>
</msup>
<mrow>
<mo>(</mo>
<msub>
<mi>E</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<mrow>
<msup>
<msub>
<mi>&mu;</mi>
<msub>
<mi>S</mi>
<mi>j</mi>
</msub>
</msub>
<mo>&prime;</mo>
</msup>
<mrow>
<mo>(</mo>
<msub>
<mi>E</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<mrow>
<msup>
<msub>
<mi>&mu;</mi>
<msub>
<mi>S</mi>
<mi>m</mi>
</msub>
</msub>
<mo>&prime;</mo>
</msup>
<mrow>
<mo>(</mo>
<msub>
<mi>E</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msup>
<msub>
<mi>&mu;</mi>
<msub>
<mi>S</mi>
<mn>1</mn>
</msub>
</msub>
<mo>&prime;</mo>
</msup>
<mrow>
<mo>(</mo>
<msub>
<mi>E</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<msup>
<msub>
<mi>&mu;</mi>
<msub>
<mi>S</mi>
<mn>2</mn>
</msub>
</msub>
<mo>&prime;</mo>
</msup>
<mrow>
<mo>(</mo>
<msub>
<mi>E</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<mrow>
<msup>
<msub>
<mi>&mu;</mi>
<msub>
<mi>S</mi>
<mi>j</mi>
</msub>
</msub>
<mo>&prime;</mo>
</msup>
<mrow>
<mo>(</mo>
<msub>
<mi>E</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<mrow>
<msup>
<msub>
<mi>&mu;</mi>
<msub>
<mi>S</mi>
<mi>m</mi>
</msub>
</msub>
<mo>&prime;</mo>
</msup>
<mrow>
<mo>(</mo>
<msub>
<mi>E</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mtable>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
</mtable>
</mtd>
<mtd>
<mtable>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
</mtable>
</mtd>
<mtd>
<mtable>
<mtr>
<mtd>
<mrow></mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow></mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow></mrow>
</mtd>
</mtr>
</mtable>
</mtd>
<mtd>
<mtable>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
</mtable>
</mtd>
<mtd>
<mtable>
<mtr>
<mtd>
<mrow></mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow></mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow></mrow>
</mtd>
</mtr>
</mtable>
</mtd>
<mtd>
<mtable>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
</mtable>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msup>
<msub>
<mi>&mu;</mi>
<msub>
<mi>S</mi>
<mn>1</mn>
</msub>
</msub>
<mo>&prime;</mo>
</msup>
<mrow>
<mo>(</mo>
<msub>
<mi>E</mi>
<mi>n</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<msup>
<msub>
<mi>&mu;</mi>
<msub>
<mi>S</mi>
<mn>2</mn>
</msub>
</msub>
<mo>&prime;</mo>
</msup>
<mrow>
<mo>(</mo>
<msub>
<mi>E</mi>
<mi>n</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<mrow>
<msup>
<msub>
<mi>&mu;</mi>
<msub>
<mi>S</mi>
<mi>j</mi>
</msub>
</msub>
<mo>&prime;</mo>
</msup>
<mrow>
<mo>(</mo>
<msub>
<mi>E</mi>
<mi>n</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<mrow>
<msup>
<msub>
<mi>&mu;</mi>
<msub>
<mi>S</mi>
<mi>m</mi>
</msub>
</msub>
<mo>&prime;</mo>
</msup>
<mrow>
<mo>(</mo>
<msub>
<mi>E</mi>
<mi>n</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mtable>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
</mtable>
</mtd>
<mtd>
<mtable>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
</mtable>
</mtd>
<mtd>
<mtable>
<mtr>
<mtd>
<mrow></mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow></mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow></mrow>
</mtd>
</mtr>
</mtable>
</mtd>
<mtd>
<mtable>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
</mtable>
</mtd>
<mtd>
<mtable>
<mtr>
<mtd>
<mrow></mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow></mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow></mrow>
</mtd>
</mtr>
</mtable>
</mtd>
<mtd>
<mtable>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>.</mo>
</mtd>
</mtr>
</mtable>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msup>
<msub>
<mi>&mu;</mi>
<msub>
<mi>S</mi>
<mn>1</mn>
</msub>
</msub>
<mo>&prime;</mo>
</msup>
<mrow>
<mo>(</mo>
<msub>
<mi>E</mi>
<mi>N</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<msup>
<msub>
<mi>&mu;</mi>
<msub>
<mi>S</mi>
<mn>2</mn>
</msub>
</msub>
<mo>&prime;</mo>
</msup>
<mrow>
<mo>(</mo>
<msub>
<mi>E</mi>
<mi>N</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<mrow>
<msup>
<msub>
<mi>&mu;</mi>
<msub>
<mi>S</mi>
<mi>j</mi>
</msub>
</msub>
<mo>&prime;</mo>
</msup>
<mrow>
<mo>(</mo>
<msub>
<mi>E</mi>
<mi>N</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mo>...</mo>
</mtd>
<mtd>
<mrow>
<msup>
<msub>
<mi>&mu;</mi>
<msub>
<mi>S</mi>
<mi>m</mi>
</msub>
</msub>
<mo>&prime;</mo>
</msup>
<mrow>
<mo>(</mo>
<msub>
<mi>E</mi>
<mi>N</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>;</mo>
</mrow>
<mrow>
<msup>
<msub>
<mi>&mu;</mi>
<msub>
<mi>S</mi>
<mi>j</mi>
</msub>
</msub>
<mo>&prime;</mo>
</msup>
<mrow>
<mo>(</mo>
<msub>
<mi>E</mi>
<mi>n</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mn>1</mn>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>&mu;</mi>
<msub>
<mi>S</mi>
<mi>j</mi>
</msub>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>E</mi>
<mi>n</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>&GreaterEqual;</mo>
<mi>&lambda;</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mn>0</mn>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>&mu;</mi>
<msub>
<mi>S</mi>
<mi>j</mi>
</msub>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>E</mi>
<mi>n</mi>
</msub>
<mo>)</mo>
</mrow>
<mo><</mo>
<mi>&lambda;</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>11</mn>
<mo>)</mo>
</mrow>
</mrow>
式(11)中，表示流动站第n个子时间段差分定位数据Xn的定位结果Ln的置信度T＝j；
步骤9、利用式(12)对流动站的N个定位结果进行加权滤波处理，获得最终定位结果L*；
<mrow>
<msup>
<mi>L</mi>
<mo>*</mo>
</msup>
<mo>=</mo>
<mfrac>
<mrow>
<mi>T</mi>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<mi>T</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<msub>
<mi>L</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<mfrac>
<mrow>
<mi>T</mi>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<mi>T</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<msub>
<mi>L</mi>
<mn>2</mn>
</msub>
<mo>+</mo>
<mn>...</mn>
<mo>+</mo>
<mfrac>
<mrow>
<mi>T</mi>
<mrow>
<mo>(</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<mi>T</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<msub>
<mi>L</mi>
<mi>n</mi>
</msub>
<mo>+</mo>
<mn>...</mn>
<mo>+</mo>
<mfrac>
<mrow>
<mi>T</mi>
<mrow>
<mo>(</mo>
<mi>N</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<mi>T</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<msub>
<mi>L</mi>
<mi>N</mi>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>12</mn>
<mo>)</mo>
</mrow>
</mrow>
式(12)中，T(k)表示流动站第k个子时间段的定位结果Lk的置信度。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510896845.8A CN105549050B (zh) | 2015-12-04 | 2015-12-04 | 一种基于模糊置信度滤波的北斗变形监测定位方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510896845.8A CN105549050B (zh) | 2015-12-04 | 2015-12-04 | 一种基于模糊置信度滤波的北斗变形监测定位方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105549050A CN105549050A (zh) | 2016-05-04 |
| CN105549050B true CN105549050B (zh) | 2017-11-28 |
Family
ID=55828350
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510896845.8A Active CN105549050B (zh) | 2015-12-04 | 2015-12-04 | 一种基于模糊置信度滤波的北斗变形监测定位方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105549050B (zh) |
Families Citing this family (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106199662B (zh) * | 2016-06-28 | 2018-10-02 | 合肥工业大学 | 一种基于层次分析法的北斗高精度定位方法 |
| CN106443735B (zh) * | 2016-11-01 | 2019-02-26 | 国网新疆电力公司信息通信公司 | 基于综合分析法的北斗高精度测量方法 |
| CN106767664B (zh) * | 2016-11-29 | 2019-05-14 | 中南安全环境技术研究院股份有限公司 | 基于北斗系统的大跨径桥梁形变实时监测方法 |
| CN111352131B (zh) * | 2020-03-24 | 2022-04-12 | 广东星舆科技有限公司 | 一种高精度定位数据处理的方法和装置 |
| CN112884315B (zh) * | 2021-02-08 | 2022-08-05 | 合肥工业大学智能制造技术研究院 | 一种基于模糊评判的高大建筑倾斜预警方法 |
| CN113311460B (zh) * | 2021-07-28 | 2021-10-26 | 湖南联智科技股份有限公司 | 一种基于北斗的预警方法及其系统 |
| CN114758504B (zh) * | 2022-06-13 | 2022-10-21 | 之江实验室 | 一种基于滤波校正的网联车超速预警方法及系统 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103148813A (zh) * | 2013-01-31 | 2013-06-12 | 湖南致力地质资源环境工程科技有限公司 | 用于处理gps变形监测数据的方法 |
| CN103471544A (zh) * | 2013-09-18 | 2013-12-25 | 柳州天运寰通科技有限公司 | 一种基于北斗的高精度位移形变监测应用系统 |
| CN104180781A (zh) * | 2014-09-10 | 2014-12-03 | 中南大学 | 一种单双频gps混合网变形监测数据处理方法 |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| FR2948188B1 (fr) * | 2009-07-20 | 2011-09-09 | Soletanche Freyssinet | Procede de surveillance des mouvements d'un terrain |
| US9778360B2 (en) * | 2013-12-17 | 2017-10-03 | Fugro N.V. | Method and system for generating a geoid via three computation spaces and airborne-acquired gravity data |
-
2015
- 2015-12-04 CN CN201510896845.8A patent/CN105549050B/zh active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103148813A (zh) * | 2013-01-31 | 2013-06-12 | 湖南致力地质资源环境工程科技有限公司 | 用于处理gps变形监测数据的方法 |
| CN103471544A (zh) * | 2013-09-18 | 2013-12-25 | 柳州天运寰通科技有限公司 | 一种基于北斗的高精度位移形变监测应用系统 |
| CN104180781A (zh) * | 2014-09-10 | 2014-12-03 | 中南大学 | 一种单双频gps混合网变形监测数据处理方法 |
Non-Patent Citations (1)
| Title |
|---|
| "基于模糊算法的GPS定位误差估计";王美玲 等;《中国惯性技术学报》;20070430;第15卷(第2期);第185-187,192页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105549050A (zh) | 2016-05-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105549050B (zh) | 一种基于模糊置信度滤波的北斗变形监测定位方法 | |
| CN102508278B (zh) | 一种基于观测噪声方差阵估计的自适应滤波方法 | |
| CN102298151B (zh) | 一种gnss网络差分定位系统中的误差改正方法 | |
| CN104680008B (zh) | 一种基于多参考站的网络rtk区域大气误差建模方法 | |
| CN108919321B (zh) | 一种基于尝试法的gnss定位粗差探测方法 | |
| CN106959456B (zh) | 一种gnss测量控制网精度估计方法 | |
| CN105043342B (zh) | 单向精密测距三角高程测量方法 | |
| CN104050513B (zh) | 农作物种植面积监测空间抽样方案优化方法 | |
| CN107622304A (zh) | 一种基于bp神经网络技术的电离层球谐函数改进方法 | |
| CN100582661C (zh) | 精确确定区域高程异常的方法 | |
| CN108254773A (zh) | 一种多gnss的实时钟差解算方法 | |
| CN105005019A (zh) | 一种基于bds/gps rtk的地下管线测量方法 | |
| CN107976702A (zh) | 一种基于cors的位置修正方法、定位终端及定位系统 | |
| Waechter et al. | Modern glacier velocities across the Icefield Ranges, St Elias Mountains, and variability at selected glaciers from 1959 to 2012 | |
| CN112099069A (zh) | 实测气象参数修正对流层经验模型的rtk算法及应用 | |
| CN106199662B (zh) | 一种基于层次分析法的北斗高精度定位方法 | |
| Vestøl et al. | NKG2016LU, an improved postglacial land uplift model over the Nordic-Baltic region | |
| CN106934113A (zh) | 适用于区域电离层垂直总电子含量建模的改进多项式模型的建模方法 | |
| CN112987058A (zh) | 一种利用地表气象站增强短距离大高差rtk定位的方法 | |
| CN104392113B (zh) | 一种近海海面冷空气大风风速的估算方法 | |
| CN115421172B (zh) | 一种基于实时与准实时结合的北斗变形监测方法 | |
| CN104992054A (zh) | 基于时间序列二维化的电离层垂直总电子含量预报方法 | |
| CN101266153B (zh) | 测绘工程类陀螺全站仪精度评定方法 | |
| CN106443735A (zh) | 基于综合分析法的北斗高精度测量方法 | |
| CN103941266A (zh) | 一种基于互差的卫星完好性监测方法及系统 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |