用vs做音乐网站,医药网站前置审批,家庭网络如何做网站服务器,网站模版防被偷《测绘学报》构建与学术的桥梁 拉近与权威的距离等价条件平差模型的方差-协方差分量最小二乘估计方法刘志平1, 朱丹彤1, 余航1, 张克非1,21. 中国矿业大学环境与测绘学院, 江苏 徐州 221116; 2. 皇家墨尔本理工大学空间科学研究中心, 澳大利亚 维多利亚州 墨尔本 3001收稿日期 2. 皇家墨尔本理工大学空间科学研究中心, 澳大利亚 维多利亚州 墨尔本 3001收稿日期2018-01-22修回日期2018-11-28基金项目国家自然科学基金重点项目(41730109);国家自然科学基金(41771416;41204011);江苏省双创团队项目(CUMT07180005);江苏省双创个人项目(CUMT07180003);自然资源部精密工程与工业测量重点实验室开放基金(PF2017-12)第一作者简介刘志平(1982—)男博士副教授研究方向为混合误差理论与反演、全源定位导航与应用。E-mail:zhpliucumt.edu.cn; zhpnliu163.com摘要提出等价条件闭合差的方差-协方差分量最小二乘估计方法简称LSV-ECM法。首先利用等价条件平差模型建立了基于等价条件闭合差二次型的方差-协方差分量估计方程由矩阵半拉直算子将其变换为线性Gauss-Markov形式进而通过最小二乘准则导出了具有模型通用性、形式简洁性且满足无偏性和最优性的方差-协方差分量估计公式。其次证明了LSV-ECM方法与残差型VCE方法的等价性并在此基础上通过计算复杂度定量分析了所提方法的计算高效性。最后通过边角网平差和中国区域GNSS站坐标时序建模及其结果分析验证了所提新方法的正确性和计算高效性。关键词等价条件平差模型 方差-协方差分量估计 LSV-ECM法 边角网 GNSS站坐标时间序列 Least-square variance-covariance component estimation method based on the equivalent conditional adjustment modelLIU Zhiping1, ZHU Dantong1, YU Hang1, ZHANG Kefei1,21. School of Environment Science and Spatial Informatics, China University of Mining and Technology, Xuzhou 221116, China;2. Space Research Centre, RMIT University, Australia VIC 3001Foundation support: The State Key Program of National Natural Science Foundation of China (No. 41730109); The National Natural Science Foundation of China (Nos. 41771416; 41204011); The Jiangsu Dual Creative Teams Programme Projects Awarded In 2017 (No. CUMT07180005); The Jiangsu Dual Creative Tanlents Programme Projects Awarded In 2017 (No. CUMT07180003); The Open Foundation of Precise Engineering Industry Surveying Key Laboratory of Natural Resources Ministry (No. PF2017-12)First author: LIU Zhiping(1982—), male, PhD, associate professor, majors in mixed error theory and geodetic inversion, all source navigation positioning and application. E-mail:zhpliucumt.edu.cn; zhpnliu163.com.Abstract: A VCE method termed the least-square variance-covariance component estimation method based on the equivalent conditional misclosure (LSV-ECM) is developed. Three steps are involved. First, the equivalent conditional misclosure is extracted using the projection matrix in the equivalent conditional adjustment model, of which the quadratic equations are established for variance-covariance component estimation. The quadratic equations in the form of matrix are then transformed to the linearized Gauss-Markov form using the half-vectorization operator. A simplified and generalized LSV-ECM method is derived using the least-square principle with an unbiased and optimal estimation.Furthermore, the equivalence between the LSV-ECM and the existing VCE methods is proven mathematically, and computational complexities of the LSV-ECM and the existing VCE methods are quantitatively analyzed and investigated in the indirect adjustment model. It is shown that the new method gives the highest computational efficiency. Finally, the performance and superiority of the new method is evaluated through an adjustment of a triangulateration network and an analysis of a coordinate time series of GNSS stations.Key words: equivalent conditional adjustment model variance-covariance component estimation LSV-ECM method triangulateration network GNSS station coordinate time series 随着现代测量技术的发展测量数据处理的对象已经由单一同类观测数据转变为多源、多类或同类多因素的异方差结构观测数据在这种情况下随机模型的准确性对参数解估计的统计特性具有重要影响[1-2]。方差-协方差分量估计能够极大地改善随机模型不精确的问题从而成为国内外学者关注的焦点并形成了一整套的理论框架[3-8]。近年来主要理论研究包括①以平差结果之一的残差向量为基本输入量的残差型VCE方法如Helmert法[1-2, 7]、最小二乘方差分量估计法(简称LS-VCE)[9-11]、MINQUE法[12]、BIQUE法[13]、MLE-VCE法[14]、等效残差法[15-16]及EIV模型的VCE方法[4, 8]等②以等价条件闭合差为基本输入量的解析型VCE方法如VCE-ECM法[17-18]。近年来的主要应用研究有文献[19— 20]研究了常规导线控制网中角度观测值与距离观测值的定权问题文献[21—22]探讨了GNSS观测值随机模型的精化问题文献[23—24]研究了坐标转换中不同公共点坐标的精度问题文献[11, 25—26]开展了GNSS站坐标时间序列的噪声估计及评价研究。大量研究表明LS-VCE方法[9-10, 15]可导出Helmert、MINQUE、MLE-VCE等方法的方差-协方差分量估计公式具有最优无偏特性。与此同时基于等价条件平差模型的VCE-ECM法[17-18]以等价条件闭合差作为基本输入量摒弃了平差值求解或残差估计过程。鉴于此本文从方差-协方差分量最小二乘估计准则和等价条件平差模型出发利用等价条件闭合差的二次型构造方差分量估计公式进而将方差分量估计方程变换为线性Gauss-Markov形式并建立了基本估计方程然后顾及矩阵拉直、半拉直算子和Kronecker积的性质进行化简导出了基于等价条件闭合差的方差-协方差分量最小二乘估计公式(least-square variance-covariance estimation based on the equivalent condition misclosure)简称为LSV-ECM法。在此基础上证明了本文LSV-ECM法与残差型VCE方法的等价性定量分析了LSV-ECM法、Helmert法和LS-VCE法的计算复杂度以边角网平差和中国区域GNSS站坐标时序建模实例验证本文方法的正确性和高效性。1 等价条件闭合差的方差分量最小二乘估计1.1 等价条件平差模型设概括平差模型为(1)式中A、[BTCT]均为行满秩系数矩阵W为具有参数的条件方程闭合差Z为限制条件方程闭合差DL为观测值L的方差阵V、x分别为待求的残差与参数向量下标c、s分别为条件方程个数和限制条件方程个数n、u分别为观测数和参数个数平差模型自由度r(cs)u。令矩阵BT[BTCT]利用正交投影矩阵H消去参数向量即(2)式中H称为零空间算子Bu、Br为B的上u行与下r行分块矩阵即BT [BuTBrT]式(1)等式两边均左乘H根据式(1)可将式(2)概括平差模型化为等价的条件平差模型(3)式中AHcA、wHcWHsZ称为等价条件闭合差。其中Hc、Hs为H的左c列与右s列分块矩阵即H [HcHs]。1.2 方差-协方差分量最小二乘估计方法假设概括平差模型式(1)中随机模型DL可表示为(4)式中∑(·)表示累加运算k为方差-协方差分量的个数Qi为第i个协因数分量矩阵。由式(3)可将等价条件闭合差W的随机模型表示为(5)令FDW根据式(4)、式(5)可得(6)式中QiAQiAT。进一步令θ[σ0, 12σ0, 22… σ0,k2]T并顾及F、Qi的矩阵对称性质利用下三角矩阵拉直算子(简称半拉直算子)vh(·)对式(6)等式两端F,Qi进行变换可得(7)式中mr(r1)/2Fvhvh(F)Qvh[vh(Q) vh(1 Q) … vh(2 Qk)]。记向量Fvh的方差阵为ΣvhN令2QvhTΣvh1QvhU2QvhTΣvh1Fvh则由最小二乘准则可得式(7)的方差-协方差分量估计式(8)(9)式中D为复制矩阵(duplication matrix)[9]也即半拉直算子与拉直算子的映射矩阵(10)分析可知式(8)—式(10)含有大量的高维矩阵迭代运算(r2阶)因而需进一步简化矩阵运算。考虑到矩阵拉直算子vec(·)、Kronecker积⊗和矩阵迹算子tr(·)具有如下性质[9](11)利用式(9)—式(11)对方差-协方差分量估计式N、U中的元素进行简化整理可得(12)(13)式(8)和式(12)—式(13)构成了本文方法方差-协方差分量估计式的最终形式。由上述推导过程可知本文方法仅含等价条件平差模型的DW、Qi和W三类输入量降低了矩阵运算维数(由r2阶降为r阶)且在最小二乘迭代估计过程中仅需更新DW(Qi和W为已知量)。对于特定的平差模型只需变换DW、Qi和W计算式便可得到相应的方差-协方差分量估计公式(表 1)。因此本文称为等价条件闭合差的方差-协方差最小二乘估计方法简称LSV-ECM方法。该方法采用概括平差模型导出的等价条件闭合差进行表达不同于Helmert型通用VCE法[7]等采用基于平差因子的残差向量表示也区别于等效残差型VCE法[15-16]采用基于正交矩阵分解的等效残差表示所得估计公式具有更好的模型通用性和形式简洁性。其次本文LSV-ECM方法利用矩阵半拉直算子将方差-协方差分量估计方程变换为线性Gauss-Markov形式并利用最小二乘准则进行方差-协方差分量估计从而保证了估计结果具有无偏性和最优性。表 1 不同平差模型的Qi、W、DW计算式Tab. 1Qi、W、DWfor mulas for different adjustment models参量模型条件平差具有参数的条件平差间接平差附有限制的间接平差B0,C0C0A-I,C0A-IQiAQiATHAQiATHTHQiHTHcQiHcTWWHWHWHcWHsZDWADLATHADLATHTHDLHTHcDLHcT表选项现有文献已证明残差型VCE方法之间的等价性并形成了一套较为完整的理论框架。通过概括平差因子矩阵R也可以证明本文LSV-ECM方法与残差型VCE方法的等价性。对于概括平差模型LSV-ECM法与Helmert型通用VCE方法[7](简称通用Helmert法)等价对于间接平差模型LSV-ECM法与LS-VCE法[9, 11]、基于等效残差的LS-VCE法[15]等价对于具有参数的条件平差模型LSV-ECM法与文献[10]相应算法等价。为节省篇幅这里仅给出LSV-ECM法与通用Helmert法[7]等价关系的证明过程。根据概括平差模型[18]DW与概括平差因子矩阵R存在如下等式关系(14)顾及式(14)及QiHcAQiATHcT和WHcAV可将式(12)—式(13)表达式Nij、Ui化简整理为关于概括平差因子矩阵R的表达式(15)(16)分析对比式(15)—式(16)和式(12)—式(13)可知LSV-ECM方法可以导出通用Helmert法[7]等价性得证。其次残差型VCE方法(如Helmert法、LS-VCE法等)均以残差向量为输入量而本文LSV-ECM法以等价条件闭合差为输入量后者实现了平差值求解(残差为平差结果之一)与随机模型估计的分离。再次等价条件闭合差维数低于残差向量维数且前者为不需求解的已知量因此计算效率方面明显优于前者。表 2以间接平差模型为例分析了3种方法在构建N、U时所需的计算复杂度(迭代一次所需的加法乘法和求逆复杂度[27])。由表 2可看出通用Helmert法与LS-VCE法以残差向量作为基本输入量故二者的矩阵求逆复杂度完全相同仅在加法乘法复杂度方面略有差异故通用Helmert法和LS-VCE法在计算效率方面基本相等而本文LSV-ECM法以等价条件闭合差为基本输入量矩阵求逆复杂度和加法乘法复杂度均小于通用Helmert法和LS-VCE法因此LSV-ECM法的计算效率高于通用Helmert法和LS-VCE法。表 2 不同VCE方法的计算复杂度Tab. 2 A comparison of computational complexities of different VCE methodsVCE方法矩阵求逆复杂度加法乘法复杂度通用Helmert[7]4O(n3)k2(6n33n2n1)k(6n2n1)LS-VCE[9, 11]4O(n3)k2(6n35n2n1)k(6n2n1)LSV-ECM4O(r3)k2(6r33r2r1)k(6r2r1)表选项2 应用结果及分析2.1 边角网平差为验证本文方法的正确性和有效性利用边角网实例进行方差分量估计和参数平差计算。设有边角网[17]其中A、B、C是已知点P1、P2是待定点网中观测了12个角度编号为1~12观测了6条边长编号为13~18。采用间接平差模型先验测角中误差为1.5″测边中误差为2 cm记后验测角方差为方案1分别采用通用Helmert法[7]、LS-VCE法[9, 11]、本文LSV-ECM法进行方差分量估计和参数平差计算。其中测角和测边单位权中误差的迭代初值均取先验中误差迭代终止条件为测角和测边的单位权中误差相等。方案2在方案1的基础上剔除含有粗差的12号角度观测值和16号边长观测值[18]基于余下的16个观测值进行方差分量估计和参数平差计算。按上述方案分别进行计算统计不同方法迭代收敛时的方差分量估值、参数估值及其中误差。另外为避免单次计算的随机误差将上述方案重复计算100次并统计LS-VCE法和本文LSV-ECM法相对于通用Helmert法的耗时比值Tratio所得结果如表 3所示。表 3 不同VCE方法的估计结果Tab. 3 A comparison of the results by different VCE methods结果方案1方案2HelmertLS-VCELSV-ECMHelmertLS-VCELSV-ECM3.64, 5.923.64, 5.923.64, 5.920.79, 3.070.79, 3.070.79, 3.07-1.56, 0.98-1.56, 0.98-1.56, 0.98-2.90, 0.53-2.90, 0.53-2.90, 0.530.89.1.030.89.1.030.89.1.030.09, 0.550.09, 0.550.09, 0.555.64, 1.645.64, 1.645.64, 1.643.48, 1.203.48, 1.203.48, 1.20-12.38, 1.86-12.38, 1.86-12.38, 1.86-17.80, 1.28-17.80, 1.28-17.80, 1.28Tratio196%65%197%71%表选项从方案1方差分量和待估参数的估计结果看3种结果完全相同表明3种方法均能获得一致的估计结果。同时比较3种方法的计算效率通用Helmert法和LS-VCE的计算效率处于同一水平本文所提出的LSV-ECM法的计算效率最高其计算时间约为通用Helmert法的65%较通用Helmert法和LS-VCE法均有明显提升从而验证本文方法的有效性。方案2中同样可以验证本文所提方法的正确性和有效性不再赘述。同时方案2在剔除两个粗差观测值以后两类观测值的单位权方差、参数的估计精度均有较大提升表明3种方差分量估计方法都无法抑制粗差影响必须引入质量控制方法进行粗差识别与剔除。2.2 GNSS站坐标时序建模为进一步检验本文LSV-ECM法的正确性和高效性选用中国大陆构造环境监测网络(crustal movement observation network of ChinaCMONOC)中16个长期线性趋势稳定的GNSS基准站2005.0014(DOY)~2015.0014(DOY)共10年的原始坐标时间序列进行噪声特性分析(方差-协方差分量估计)和三维运动速度估计(平差参数计算)测站位置分布如图 1所示。图 1 CMONOC所选测站分布Fig. 1 Geographical distribution of selected stations in the CMONOC network图选项GNSS站坐标时序的观测方程和随机模型分别为[11](17)(18)式中ti是以年为单位的时间序列历元点a为常数项b为线性速度项c、d组合表示全年性周期运动e、f组合表示半年性周期运动vi为残差σWN、σFN为所求的噪声分量大小QWN、QFN分别为白噪声和闪烁噪声的协因数阵具体形式见文献[11]。基于上述间接平差模型分别利用通用Helmert法[7]、LS-VCE法[9, 11]和本文LSV-ECM法进行方差分量估计和参数平差计算并统计北(N)、东(E)、竖直(U)3个方向的噪声分量Tratio结果如表 4所示。表 4 噪声分量(WN, FN)估计结果Tab. 4 Noise analysis results of WN and FN using different VCE methodsmm站点NEUHelmertLS-VCELSV-ECMHelmertLS-VCELSV-ECMHelmertLS-VCELSV-ECMWN, FNWN, FNWN, FNWN, FNWN, FNWN, FNWN, FNWN, FNWN, FN注高斯白噪声(WN)的单位为mm闪烁噪声(FN)的单位为mm/a0.25ZHZC0.80, 2.110.80, 2.110.80, 2.110.78, 2.980.78, 2.980.78, 2.983.43, 9.733.43, 9.733.43, 9.73YANC0.52, 1.670.52, 1.670.52, 1.670.45, 1.870.45, 1.870.45, 1.872.10, 6.732.10, 6.732.10, 6.73XIAM0.79, 3.570.79, 3.570.79, 3.570.93, 3.190.93, 3.190.93, 3.193.78, 13.133.78, 13.133.78, 13.13WUHN0.76, 3.230.76, 3.230.76, 3.230.77, 3.030.77, 3.030.77, 3.033.01, 13.033.01, 13.033.01, 13.03TAIN0.96, 2.310.96, 2.310.96, 2.310.92, 3.100.92, 3.100.92, 3.103.18, 8.773.18, 8.773.18, 8.77QION0.98, 5.200.98, 5.200.98, 5.201.27, 4.631.27, 4.631.27, 4.634.7, 15.014.7, 15.014.7, 15.01LUZH0.68, 2.190.68, 2.190.68, 2.190.62, 2.630.62, 2.630.62, 2.632.73, 12.192.73, 12.192.73, 12.19KMIN0.67, 3.730.67, 3.730.67, 3.730.65, 4.940.65, 4.940.65, 4.943.29, 14.573.29, 14.573.29, 14.57JIXN0.58, 2.520.58, 2.520.58, 2.520.61, 1.980.61, 1.980.61, 1.982.37, 6.752.37, 6.752.37, 6.75HRBN0.61, 3.150.61, 3.150.61, 3.150.33, 3.450.33, 3.450.33, 3.451.50, 11.831.50, 11.831.50, 11.83HLAR0.75, 2.970.75, 2.970.75, 2.970.70, 2.650.70, 2.650.70, 2.652.32, 11.502.32, 11.502.32, 11.50GUAN1.07, 3.411.07, 3.411.07, 3.411.21, 5.081.21, 5.081.21, 5.085.23, 16.225.23, 16.225.23, 16.22DLHA0.41, 2.500.41, 2.500.41, 2.500.51, 2.310.51, 2.310.51, 2.311.51, 10.491.51, 10.491.51, 10.49CHUN0.55, 3.090.55, 3.090.55, 3.090.44, 3.050.44, 3.050.44, 3.051.92, 13.271.92, 13.271.92, 13.27BJSH0.76, 1.990.76, 1.990.76, 1.990.76, 1.350.76, 1.350.76, 1.352.81, 7.452.81, 7.452.81, 7.45BJFS0.61, 3.490.61, 3.490.61, 3.490.62, 2.190.62, 2.190.62, 2.192.73, 7.312.73, 7.312.73, 7.31Tratio(WN)199.5%74.0%199.8%73.8%199.9%72.1%Tratio(FN)199.5%74.0%199.8%73.8%199.9%72.1%表选项对比表 4中通用Helmert法、LS-VCE法和LSV-ECM法的计算结果可知3种方法所得的噪声分量结果完全相同验证了本文方法的正确性。同时对比3种方法的计算效率可知通用Helmert法与LS-VCE法基本一致而本文LSV-ECM法有较大提高三维方向计算时间比分别为通用Helmert法的74.0%、73.8%、72.1%表明本文方法的高效性。由表 4噪声分量的估计结果可知白噪声分量均远小于闪烁噪声分量表明有色噪声为中国区域GNSS站坐标时间序列的主要噪声在参数估计时若直接采用白噪声模型会导致估计结果有偏并会产生较高的虚假估计精度。另外对比不同方向的噪声分量估计结果北方向和东方向的噪声分量相差不大且量级较小94%和88%的测站白噪声分量在1 mm以内50%的测站闪烁噪声分量在3 mm/a0.25以内而竖直方向的噪声分量远高于水平方向这与现有结论相一致约有56%的测站白噪声分量在3 mm以内31%的测站闪烁噪声分量在9 mm/a0.25以内。为进一步分析噪声大小与经纬度之间的关系绘制了白噪声分量和闪烁噪声分量随经纬度的变化情况。由于噪声大小与纬度相关性较强、与经度相关性较弱为节省篇幅此处仅讨论噪声大小与纬度变化的关系见图 2。从图 2可以看出水平方向的白噪声和闪烁噪声分量大小较为平稳竖直方向的噪声大小波动较大。从整体趋势上看白噪声和闪烁噪声大小随纬度变化的趋势较为明显。具体表现为噪声大小随纬度增加而逐渐减小减少过程中有轻微波动且在中纬度地区有翘尾现象。该现象原因分析可能是GNSS的GDOP值随纬度增大而降低[28]从而使噪声分量表现出随纬度增加而减小的现象而且噪声分量在中纬度地区呈现一定的翘尾现象。图 2 噪声分量大小与纬度的关系Fig. 2 The relationship between noise components and latitude图选项此外为检验方差-协方差分量估计结果的正确性和有效性以CMONOC的公布数据(http://www.cgps.ac.cn/cgs/index.action)作为对比参考值表 5统计了中国区域GNSS站的三维运动速度和不确定度。由该表可知中国区域不同站点水平方向均呈现东南方向运动趋势这与亚欧板块的整体运动趋势相一致但竖直方向运动变化差异性较大其中约40%的测站呈现下沉趋势。对比本文结果与CMONOC的公布结果除JIXN、HLAR站N方向KMIN、GUAN站E方向CHUN、BJFS站U方向以外其余计算结果均在2倍中误差范围内。此外参照文献[29]的计算结果对比发现除CHUN站U方向外其余测站的计算结果与本文结果具有较好的相符性进一步验证了本文方法的正确性。需要指出的是CHUN站U方向的差异性可能是数据预处理和共模误差提取策略不同所致。表 5 站点运动速度和不确定度Tab. 5 Velocity and uncertainty of different stationsmm/a站点NEU本文方法CMONOC本文方法CMONOC本文方法CMONOCZHNZ-11.35±0.07-11.18±0.1532.90±0.1133.17±0.231.21±0.311.24±0.38YANC-9.30±0.05-8.69±0.1432.57±0.1032.46±0.051.06±0.221.02±0.13XIAM-12.16±0.11-12.48±0.1732.47±0.0632.82±0.171.4±0.420.78±0.36WUHN-10.93±0.10-11.08±-0.1232.50±0.1133.60±0.74-0.79±0.420.47±0.41TAIN-11.54±0.07-11.58±0.3730.98±0.1031.38±0.160.92±0.281.21±0.30QION-11.98±0.17-10.24±0.7831.57±0.1031.94±0.15-0.6±0.48-0.46±0.32LUZH-9.75±0.070-9.61±0.1334.96±0.1535.76±0.280.46±0.390.35±0.30KMIN-17.23±0.12-16.18±0.8733.09±0.0831.13±0.53-1.23±0.47-0.48±0.31JIXN-9.72±0.08-10.35±0.0629.13±0.1628.66±0.111.79±0.221.53±0.17HRBN-12.56±0.10-12.38±-0.2125.79±0.0625.95±0.51-0.49±0.38-0.09±0.21HLAR-10.51±0.09-11.35±0.0425.76±0.1125.88±0.072.05±0.371.44±0.16GUAN-11.11±0.11-11.23±0.1031.30±0.1033.11±0.19-0.33±0.52-1.97±0.46CHUN-11.58±0.10-12.21±0.1427.37±0.0726.53±0.52-1.81±0.42-0.15±0.29BJSH-11.45±0.06-11.28±0.1830.10±0.1029.94±0.131.36±0.241.05±0.29BJFS-9.94±0.11-10.18±0.1430.58±0.0430.10±0.192.63±0.24-0.11±0.63表选项3 结论本文基于等价条件平差模型和最小二乘准则利用等价条件闭合差的二次型构建方差-协方差分量估计方程并通过矩阵半拉直算子将其变换为线性Gauss-Markov形式进而顾及矩阵拉直算子、半拉直算子和Kronecker积运算性质导出了基于等价条件平差模型的方差-协方差分量最小二乘估计公式简称LSV-ECM法。该方法实现了平差值求解(残差为平差结果之一)与随机模型估计的分离有效兼顾了等价条件闭合差(已知)和最小二乘的特性。边角网平差实例的结果表明本文的LSV-ECM法与通用Helmert法、LS-VCE法的结果完全相同但计算效率更高验证了该方法与残差型VCE方法的等价性和计算高效性。分析指出了方差-协方差分量估计方法相较于常规GNSS-MLE法进行GNSS站坐标时间序列噪声分析的优势并利用LSV-ECM法、通用Helmert法、LS-VCE法计算分析了中国区域16个GNSS站坐标时间序列在WNFN模型下的噪声估计和站点速度。估计结果表明有色噪声是中国GNSS站坐标时序的主要噪声且白噪声和有色噪声分量随纬度增大而减小速度估计结果与陆态网络的公布结果基本一致水平方向的整体运动趋势呈东南方向而竖直方向的运动趋势差异性较大。因此计算结果也进一步验证了本文LSV-ECM法的正确性和高效性。【引文格式】刘志平, 朱丹彤, 余航, 等. 等价条件平差模型的方差-协方差分量最小二乘估计方法. 测绘学报201948(9)1088-1095. DOI: 10.11947/j.AGCS.2019.20180227权威 | 专业 | 学术 | 前沿微信投稿邮箱 | song_qi_fan163.com欢迎加入《测绘学报》作者QQ群 297834524进群请备注姓名单位稿件编号
