数值计算方法李维国课后答案（测绘学报王珂）

数值计算方法李维国课后答案（测绘学报王珂）WANG Ke LIU Guolin FU Zhengqing et al. A separable nonlinear least squares solution method based on Moore-Penrose generalized inverse and solid matrix[J]. Acta Geodaetica et Cartographica Sinica 2022 51(3): 340-350. DOI: 10.11947/j.AGCS.2022.20200597 引文格式：王珂 刘国林 付政庆 等. 基于Moore-Penrose广义逆及立体矩阵的可分离非线性最小二乘解算方法[J]. 测绘学报，2022，51(3)：340-350. DOI: 10.11947/j.AGCS.2022.20200597基金项目：国家自然科学基金(420740

本文内容来源于《测绘学报》2022年第3期（审图号GS（2022）1460号）

基于Moore-Penrose广义逆及立体矩阵的可分离非线性最小二乘解算方法

王珂^{1 2}

刘国林²

付政庆² 王路遥²

1. 山东理工大学建筑工程学院 山东 淄博 255000;2. 山东科技大学测绘与空间信息学院 山东 青岛 266590

基金项目：国家自然科学基金(42074009);山东省自然科学基金(ZR2020MD043)

摘要：针对测绘领域中函数模型为非线性函数的线性组合的特殊结构 本文提出了基于Moore-Penrose广义逆和立体矩阵的可分离非线性最小二乘解算方法。该方法首先利用变量投影算法消除可分离非线性模型中的线性参数 将包含两类参数的原非线性优化问题转化为仅含有非线性参数的最小二乘问题。然后 基于Moore-Penrose广义逆矩阵的微分和立体矩阵理论计算最小二乘目标函数的一阶导数 进而采用非线性优化的LM方法求解非线性参数的最优估值。最后 根据最小二乘方法求解线性参数的最优估值。通过指数函数模型拟合和机载LiDAR全波形参数求解试验与传统参数不分离优化方法进行对比 结果表明 基于Moore-Penrose广义逆和立体矩阵的可分离非线性最小二乘解算方法对待求参数初值依赖性低 同时避免了迭代过程中线性参数导致的病态问题 算法稳定性好 为测绘领域中可分离非线性最小二乘问题的解算提供了一种思路 也拓展了可分离非线性最小二乘方法的应用。

关键词：可分离非线性最小二乘方法 非线性模型 参数估计 变量投影算法

引文格式：王珂 刘国林 付政庆 等. 基于Moore-Penrose广义逆及立体矩阵的可分离非线性最小二乘解算方法[J]. 测绘学报，2022，51(3)：340-350. DOI: 10.11947/j.AGCS.2022.20200597

WANG Ke LIU Guolin FU Zhengqing et al. A separable nonlinear least squares solution method based on Moore-Penrose generalized inverse and solid matrix[J]. Acta Geodaetica et Cartographica Sinica 2022 51(3): 340-350. DOI: 10.11947/j.AGCS.2022.20200597

阅读全文：http://xb.sinomaps.com/article/2022/1001-1595/2022-3-340.htm

引 言

非线性最小二乘法作为基本的非线性优化方法，广泛应用于非线性函数模型的参数估计中^[1-6]。然而，根据实际问题的来源，许多非线性函数模型具有特殊结构，若函数模型中待求参数一部分变量是线性的，一部分变量是非线性的，并且函数模型为非线性函数的线性组合形式，数学上称为可分离非线性最小二乘(separable nonlinear least squares，SNLS)问题^[7-9]。测绘领域中存在许多这种模型结构的参数估计问题，如激光雷达(light detection and ranging，LiDAR)全波形回波信号高斯函数分解模型中，高斯特征参数的求解^[10-11]；空间直角坐标转换模型中，旋转参数、平移参数和尺度参数的求解^[12-13]；神经网络预测模型中，激活函数的中心位置、宽度和权值参数的求解^[14]等。

传统求解可分离非线性最小二乘问题的方法是将所有的待求参数(非线性和线性)均看作非线性参数，然后进行泰勒级数展开，利用非线性迭代法或直接搜索法寻找问题的最优解^[15-16]。文献[17]将LiDAR波形数据根据高斯分量的重要性进行筛选，然后利用非线性Levenberg-Marquardt(LM)优化算法求解各波形分量非线性参数(波峰位置、波形宽度)和线性参数(振幅)的最优值。文献[18]对机载激光测深波形分解中LM与Expectation Maximization(EM)参数优化方法进行了比较，就波形拟合效果而言，模拟数据分析和实测数据试验均说明LM方法通常可得到比EM方法更接近原始波形的拟合结果，有利于进一步的波形分析，但受波形数据质量的影响较大。文献[19]针对三维坐标转换模型参数求解中先验精度未知和观测数据质量不佳的问题，以坐标转换后的点位残差加权平方和最小为原则，提出了基于Nelder-Mead单纯形直接搜索的非线性参数求解方法，提高了三维坐标转换参数的求解质量。文献[20]针对非线性最小二乘求解病态问题的不稳定问题，提出了数值收敛解和收敛效率更优的Frozen-Barycentre迭代法。文献[21]针对参数迭代过程中雅可比矩阵计算复杂的问题，结合自动微分和隐式迭代预处理技术，采用最速下降法和信赖域法求解非线性最小二乘问题。文献[22]依赖对精度水平的控制，当检测到精度太低无法进行优化时，对函数逼近方法进行改进，提出了一种基于动态精度评价函数和梯度函数的大规模非线性最小二乘的LM迭代算法，并证明了该算法的全局和局部收敛性。文献[23]对基于变量投影算法参数分离的非线性最小二乘问题，提出采用LM优化方法对非线性参数进行求解以及采用线性最小二乘或基于L-曲线的谱修正迭代方法对线性参数进行估计的混合算法，丰富了可分离非线性最小二乘的参数估计方法，但对基于变量投影算法参数分离的目标函数中残差向量雅可比矩阵的计算方法并未进行深入研究。综上分析可知，非线性参数估计的研究主要针对模型的构建和参数优化过程进行改进，取得了很好的结果^[24-25]，但很少有考虑函数模型可分离的特殊结构，并且将所有待求参数均作为非线性参数进行解算时，对参数初值的要求较高^[26]。

本文从测绘领域中非线性函数模型的特殊结构出发，首先，将线性参数通过变量投影算法用非线性函数表示，基于Moore-Penrose广义逆矩阵微分和立体矩阵理论推导了残差向量雅可比矩阵的计算过程；然后，对参数分离后的非线性最小二乘问题利用LM优化方法进行求解，得到非线性参数的迭代方向和最优估值，进而采用最小二乘方法计算线性参数的最优解；最后，通过模拟试验和LiDAR波形参数求解试验验证了基于Moore-Penrose广义逆矩阵微分的可分离非线性最小二乘方法的有效性。

1 可分离非线性最小二乘法原理

设测量平差中的非线性函数模型和随机模型为

(1)

式中，f(x) 为非线性函数向量；x为待求参数向量；L为观测值向量；σ²为单位权方差；矩阵Q和P分别为观测向量协因数阵和权阵。在非线性最小二乘准则下建立目标函数为

(2)

为便于讨论，假设等精度观测，权阵为单位阵(P=I)，参数求解问题则等价于求解非线性最小二乘优化问题^[27-28]，即通过式(3)来确定最优参数

(3)

式中，‖·‖表示欧氏范数。

在式(1)的非线性函数模型中，若待求参数x由线性参数向量a和非线性参数向量b两部分组成，且fi(x)为非线性函数的线性组合，将其可以写成以下形式

(4)

式中，p是线性参数a的维数；q是非线性参数b的维数；令n=p q，n表示待求参数的总个数；m表示观测次数。

令Φ(b)的列对应的是与参数向量b相关的非线性函数ϕ_j(b) (j=1 2 … p)。那么，非线性函数模型写成矩阵形式为

(5)

则式(5)的非线性最小二乘优化问题为

(6)

若给定非线性参数b，则线性参数a的最小二乘最小范数解为

(7)

式中，Φ (b)是Φ(b)的Moore-Penrose广义逆^[29]。

令向量L在Φ(b)列向量张成的线性空间上的正交投影为PΦ(b)=Φ(b)Φ (b)，其正交补空间为PΦ(b)^⊥=I-PΦ(b)，将式(7)代入式(6)可得

(8)

则原问题转化为仅含有非线性参数的最小二乘问题。基于式(8)的算法也称为变量投影算法^[9]。

经参数分离后的平差函数模型为

(9)

2 基于Moore-Penrose广义逆矩阵微分的变量投影算法

根据最小二乘准则，待求参数的最优值通过求极值方法得到^[30]，其可分离非线性最小二乘的目标函数F(b)的一阶导数为

(10)

因此，只需计算变量投影算子PΦ(b)的一阶导数

。根据Moore-Penrose广义逆矩阵的性质：Φ(b)Φ (b)Φ(b)=Φ(b) Φ (b)Φ(b)Φ (b)=Φ (b)，有PΦ(b)=PΦ(b)PΦ(b)，由此可知PΦ(b)为幂等矩阵，那么

(11)

又由于PΦ(b)Φ(b)=Φ(b)，因此

(12)

将式(12)代入式(11)，得变量投影算子的一阶偏导为

(13)

式中，

是一个m×p×q的立体矩阵，m、p、q分别表示矩阵的行、列和面。设有m次观测，残差向量中非线性函数矩阵为

(14)

那么，的第i(i=1 2 … q)面由Φ(b)对b_i的微分构成，即

(15)

根据立体矩阵的计算规则^[31-32]，式(13)为

(16)

最终式(10)的目标函数一阶导数为

(17)

则残差向量的Jacobian矩阵J(b) 为

(18)

对目标函数F(b)用泰勒级数公式展开，并取至二阶项，由以下优化模型来得到迭代方向

(19)

由极值条件可得

(20)

式中

(21)

因此，式(20)迭代方向d_k为

(22)

由于矩阵J(b)对参数b求导计算复杂，对于小残差问题，常用Gauss-Newton算法通过忽略式(22)中的二次项信息，但对于大残差的问题，二次项信息的省略，导致该算法不能很好收敛。因此，为了兼顾计算量和算法的收敛性，LM算法通过增加μI项，其中μ>0，保证了J^T(b)J(b) μI始终为正定的，并且正参数μ的引入也避免了矩阵J^T(b)J(b) 接近奇异时，搜索方向的模‖d_k‖过大的问题。

设第k次迭代的非线性参数向量为b_k，从而基于变量投影的可分离非线性最小二乘LM算法的迭代方向d_k为

(23)

式中，μ_k>0，J(bk)是残差向量V(b_k)(式(9))关于b_k的一阶偏导构成的矩阵。

令V=[V₁V₂… V_m]^T，V_i的梯度为

(24)

残差向量的雅可比矩阵J(b) 的每一行由相应向量函数的梯度转置得到，即

(25)

直接由式(13)可得

(26)

然后计算新的迭代点

(27)

直到目标函数满足迭代终止条件，求解得到非线性参数向量b的最优解，结合式(7)计算线性参数的最优解。即

(28)

令系数矩阵

(29)

设观测值向量L的协因数阵为Q_LL，单位权中误差为

，根据协方差传播定律可知，非线性参数和线性参数的方差-协方差估计公式^[33-34]为

(30)

结合以上参数求解过程，基于变量投影算法参数分离的LM优化方法步骤如下。

(1) 给定非线性参数初值b₀以及阈值0 <ε< < 1，并置迭代次数k=0。

(2) 计算非线性函数矩阵Φ(b_k)和雅可比矩阵J(b_k)。

(3) 根据式(23)和式(27)计算可分离非线性参数的迭代方向d_k和第k 1次的参数值bk 1。

(4) 若目标函数梯度的范数‖ g_k‖₂²=‖J^T(b_k)V(b_k)‖₂²≤ε，停止迭代；否则k=k 1，转至步骤(2)。

(5) 根据式(7)计算线性参数的解，进而得到所有参数的最优解( )。

当非线性函数矩阵Φ(b_k)为列满秩且良态时，仅含有非线性参数的可分离非线性最小二乘问题(式(8))与原问题(式(6))参数估计结果一致，即，若非线性参数是可分离非线性最小二乘局部(全局)最优解，=Φ L是线性参数的局部(全局)最优解，那么，( )也是原非线性最小二乘的局部(全局)最优解^[22]。

当非线性函数矩阵Φ(b_k)秩亏或病态时，则需要通过解算不适定问题的方法进行求解^[35-36]，参数收敛结果也会有所不同。

3 试验分析

为了验证基于Moore-Penrose广义逆矩阵微分的可分离非线性最小二乘方法在参数估计中的有效性，采用指数函数拟合的模拟试验及LiDAR全波形分解试验与传统的参数不分离的非线性最小二乘LM优化方法进行了对比，试验环境为Matlab 2016b，2.80 GHz PC，Windows10系统。

3.1 指数函数模型的拟合

模拟试验源于放射性物质的衰变过程，其函数模型表达式为

(31)

式中，β=(β₁ β₂)^T，λ=(λ₁ λ₂)^T是待求衰减因子参数；t_i是观测时刻；ξ是观测误差。假设在没有观测误差的影响下，衰减因子的真值为β₁^*=3.1 β₂^*=2.8 λ₁^*=10.1 λ₁^*=1.4，令ξ是标准差为0.1的高斯噪声，试验随机生成时间t为[0 2]范围内均匀间隔的201个观测数据y_i i=1 2 … m。

为了获得待求参数的最优估计，首先根据实际观测值列出误差方程，并通过变量投影算法将线性参数进行消除，得到基于变量投影算法的误差方程，然后基于Moore-Penrose广义逆矩阵的微分和立体矩阵理论计算目标函数的雅可比矩阵，采用LM算法对非线性参数λ=[λ₁λ₂]^T进行估计，当目标函数(残差平方和)的梯度小于10^-6时，迭代终止。

结合衰变模型的结构特点，根据参数初值选择在真值附近或初值选择常用的0和1得到4组不同的参数初值为[3 2 10 1]^T、[0 1 10 1]^T、[3 2 0 1]^T和[0 1 0 1]^T，然后采用本文提出的参数分离的LM优化方法对非线性参数进行优化，与参数不分离的LM优化方法进行结果对比。给定4组参数初值，参数分离方法和参数不分离方法得到的衰减因子估值均为[3.30 2.68 9.54 1.32]^T。由该参数估值得到的拟合曲线如图 1所示。

图 1 参数不分离方法和参数分离方法的拟合曲线对比 Fig. 1 Comparison of fitted curves between parameter non-separation method and parameter separation method

图选项

由图 1可知，参数分离方法和参数不分离方法得到的指数函数模型曲线与观测值拟合较好。

从迭代次数、残差平方和、均方根误差、拟合优度、相关系数和最大差值等评价指标对参数分离的LM优化方法得到的结果进行分析，为了减少计算机性能等影响因素对时间的影响，将线性参数和非线性参数初值均在真值附近的试验中参数不分离方法的计算时间设为单位1，其余试验的时间对其进行比值，评价指标结果见表 1。

表 1 指数模型拟合试验参数不分离LM方法和参数分离LM方法结果的评价指标对比 Tab. 1 Comparison of evaluation indexes between LM method without parameter separation and LM method with parameter separation in exponential model fitting experiment

结合参数最终估值结果和表 2可知，4组试验的评价指标除迭代次数和迭代时间外，其余指标结果相同，并且从拟合优度和相关系数的数值可以看出，参数分离方法和参数不分离方法都能与观测值拟合较好；迭代次数方面，当目标函数(残差平方和)的梯度小于10^-6时，由于第1组和第2组的非线性参数初值相同，参数分离的LM优化方法不受线性参数初值的影响，因此迭代次数均为7，参数不分离方法由于线性参数的初值不同，迭代次数分别为23次和201次；第3组和第4组的非线性参数初值相同，参数分离的LM优化方法的迭代次数均为15，参数不分离方法的迭代次数分别为1808次和1794次；迭代时间方面，第1组试验中虽然参数分离的LM优化方法迭代次数少，但是由于迭代过程中目标函数的一阶导数计算复杂，导致计算时间比参数不分离的LM优化方法多了约23%，但第2组、第3组和第4组试验中，参数分离方法的计算时间均比参数不分离方法少，特别是第3组和第4组试验，由于参数分离的LM优化方法比参数不分离的LM优化方法的迭代次数分别减少了1793次和1779次，因此参数分离方法的计算时间明显减少，计算效率得到显著提高。另外，第1组和第2组试验参数分离方法的时间理论上应该完全相同，但是由于计算机实时状态的不同，所以计算时间有少许差别，第3组和第4组试验的计算时间也是如此。

表 2 不同参数初值情况下残差向量雅克比矩阵条件数对比 Tab. 2 Comparison of condition number of Jacobian matrix of the residual vector with different initial parameters

为了验证参数分离方法对原问题病态性的改进，对基于变量投影算法参数分离的残差向量雅可比矩阵的条件数进行分析，其4组初值对应的参数分离方法和参数不分离方法的条件数结果见表 2。

由表 2可知，同一组参数初值的参数分离方法雅克比矩阵的条件数比参数不分离方法更小，特别是第2组和第4组参数初值，原非线性最小二乘问题的雅可比矩阵的条件数均超过10¹⁸，存在明显病态问题，而基于变量投影参数分离方法的雅可比矩阵的条件数分别为16.66和5.23，原问题的病态性得到明显改善。

3.2 基于可分离非线性最小二乘的LiDAR波形分解

试验采用海南某部分海域2012年12月测量得到的机载LiDAR测深数据进行波形分解，其波形数据采用Optech Aquarius机载测深系统测量得到，激光脉冲频率为70 kHz，以1~2 GHz的高采样率记录各采样点的振幅，采样间隔为1 ns，采样数量为287。由于在反射率低的情况下，电压值会很小，因此，将观测值进行量化，得到数字信号值(digital number)，即每一条波形由(t_i DN_i) i=1 2 …，287组成。利用双高斯模型对海面回波和海底回波建立模型，即

(32)

式中，a=[a₀a₁a₂]^T，b=[μ₁μ₂σ₁σ₂]^T是待求线性参数和非线性参数；ξ是测量误差。

利用极大值检测法及原始观测值的图像特点假设线性参数和非线性参数的近似值为

=[200 700 150]^T，

=[30 130 3 3]^T，然后根据非线性参数和线性参数是否在近似值附近选取4组初值，分别采用参数分离的LM优化方法和参数不分离的LM优化方法进行参数估计，当目标函数(残差平方和)的梯度小于10^-6时，迭代终止，根据参数估值结果得到拟合曲线如图 2所示。

图 2 不同参数初值的参数不分离方法和参数分离方法的LiDAR波形拟合曲线 Fig. 2 Comparison of LiDAR waveform fitted curves between parameter non-separation method and parameter separation method with different initial values of parameters

图选项

由图 2可知，参数分离的LM优化方法根据第1、2和3组的参数初值得到LiDAR波形曲线与原始观测值拟合较好；参数不分离的LM优化方法第1组和第2组试验得到的LiDAR波形曲线与参数分离方法相比，在第1个波峰处相差较大，第2个波峰拟合较好，第3组试验得到的LiDAR波形在两个波峰处与参数分离方法得到的拟合曲线都存在差值；第4组试验无论是参数分离方法还是参数不分离方法都未能与观测值进行较好拟合，也说明了两种方法都对参数初值具有一定的依赖性。

从迭代次数、计算时间，残差平方和、均方根误差、拟合优度、相关系数和最大差值等评价指标与参数不分离的非线性最小二乘LM优化方法进行对比，与模拟试验一样，并将第1组参数初值在近似值附近的参数不分离方法的计算时间设为单位1，结果见表 3。

表 3 LiDAR波形拟合试验参数不分离的LM方法和参数分离的LM方法结果的评价指标对比 Tab. 3 Comparison of evaluation indexes between LM method without parameter separation and LM method with parameter separation in LiDAR waveform fitting experiment

表 3中加粗的值表示离近似值较远的初值。由表 3可知，当线性参数初值和非线性参数初值都选择在近似值附近时(第1组)，参数分离方法和参数不分离方法的迭代次数都是10，迭代终止时，参数不分离方法的残差平方和是6.31×10⁴，参数分离方法的残差平方和是6.13×10⁴，从均方根误差、拟合优度、相关系数和最大偏差4个评价指标也可以看出，参数分离方法得到的参数估计值精度更高，但是相同的迭代次数参数分离方法计算所需的时间更多。

当非线性参数初值选择在最优值附近，线性参数初值任意选择时(第2组)，参数分离方法由于只与非线性参数有关，所以与第1组试验结果相同；参数不分离方法与第1组试验相比，迭代次数增加1，均方根误差由14.82增大为16.67，结合其余评价指标可知，参数分离方法的拟合精度更高，但所需计算时间更多。

当线性参数初值选择在最优值附近，非线性的部分参数(波形宽度)初值任意选择时(第3组)，参数分离方法仍然可以得到与第1组试验相同的结果；而参数不分离方法均方根误差由第1组的14.82增加为21.52，拟合优度也由0.979 3降至0.956 4，由此可知，参数不分离方法更容易受参数初值的影响。

当线性参数初值选择在最优值附近，非线性的部分参数(波峰位置)初值任意选择时，参数不分离方法和参数分离方法都未得到与第1组精度一致的参数估值，结合图 2可知，参数分离方法和参数不分离方法都对波峰位置参数的初值比较敏感。

对上述4组试验中参数分离方法和参数不分离方法的残差平方和变化过程进行对比，结果如图 3所示。

图 3 LiDAR波形拟合试验参数不分离方法和参数分离方法残差平方和的变化 Fig. 3 The change of the residual sum of squares between the parameter non-separation method and the parameter separation method in LiDAR waveform fitting experiment

图选项

由图 3可知，第1、2和3组在迭代终止时，参数分离方法的残差平方和更小，也说明了参数分离方法在求解可分离非线性最小二乘问题时的参数估计精度更高，第4组试验在迭代终止时虽然参数不分离方法的残差平方和更小，但结合表 2可知，两种方法得到均方根误差(99.433 0和102.286 4) 都与参数初值在近似值附近(第1组)的结果(14.616 5)相差较大，说明两种方法均未得收敛到待估参数的最优值。

为了更全面地对比参数不分离方法和参数分离方法的收敛结果，根据LiDAR波形分解中常用的峰值探测法得到波峰位置参数和振幅参数的近似初值，波峰宽度参数的近似初值由峰值一半处对应的两点位置得到，然后以近似初值为中心，初值的一半为邻域设定线性参数和非线性参数的取值区间，在均匀分布区间内随机取200组线性参数和非线性参数的初值，进行迭代建模，进而分析200组试验中两种方法的均方根误差、拟合优度、相关系数和最大偏差等指标的分布情况，结果如图 4和表 4所示。

图 4 LiDAR波形拟合200次试验的效果评价指标分布情况 Fig. 4 The distribution of evaluation indexes in 200 experiments of LIDAR waveform fitting

图选项

表 4 效果评价指标的统计结果 Tab. 4 The statistical results of effect evaluation indexes

由图 4和表 4可知，参数分离方法的均方根误差均值比参数不分离方法小34.721 3，说明整体上参数分离方法的精度更高；200次试验中，拟合优度大于0.9的参数分离方法占74.5%，参数不分离方法占32.5%，相关系数大于0.9的参数分离方法占76.5%，参数不分离方法占35%；由两种方法相关系数的标准差数值可知，参数分离方法200次试验的相关系数标准差更小，说明其算法稳定性更好。

对200次试验中非线性参数初值相同的参数分离方法的残差向量雅可比矩阵的条件数与参数不分离方法进行对比，结果如图 5所示。

图 5 200次试验中残差向量雅可比矩阵的条件数对比 Fig. 5 The comparison of the condition number of the Jacobian matrix of residual vector in 200 experiments

图选项

综上分析可知：①参数分离的LM优化方法只对非线性参数初值具有一定的依赖性，比传统参数不分离的LM方法对待求参数初值的要求更低；②参数分离的LM优化方法在对参数优化过程中，所需的迭代次数更少，计算效率更高；③若两种方法的迭代次数相同，由于参数分离的LM优化方法在迭代过程中需要对Moore-Penrose广义逆矩阵的微分和立体矩阵进行计算，增加了计算复杂度，因此需要的时间更多。

4 结论

基于Moore-Penrose广义逆矩阵微分的可分离非线性最小二乘方法通过利用变量投影算法将原来含有两类参数的优化问题转化为仅含有非线性参数的最小二乘问题，降低了待求参数的维数，迭代过程只依赖非线性参数的初值，避免了线性参数初值导致的病态问题。机载LiDAR测深全波形数据拟合试验结果表明，参数分离后的非线性最小二乘方法在波形参数求解过程中，参数估计精度更高，算法稳定性也更强，为机载LiDAR波形分解参数的求解提供了新的方法，也拓展了可分离非线性最小二乘方法在测绘领域的应用，但基于Moore-Penrose广义逆矩阵的微分和立体矩阵的雅可比矩阵的显式表达增加了计算复杂度，导致得到相同参数估计精度时，计算时间更多，因此如何提高计算效率也是可分离非线性最小二乘解算需要进一步研究的内容。

作者简介

第一作者简介：王珂(1990—) 女 博士生 讲师 研究方向为测量数据处理理论及应用。E-mail: wke@sdut.edu.cn

通信作者：刘国林 E-mail: gliu@sdust.edu.cn

初审：张艳玲

复审：宋启凡

终审：金 君

往期推荐

资讯

○ 贵州省第二测绘院2022年公开招聘（地信、测绘遥感等相关专业）

○ 自然资源部国土卫星遥感应用中心正式对外发布卫星激光测高产品

○ 自然资源北京市卫星应用技术中心揭牌

○ 中国地质大学康志忠教授荣获国际摄影测量与遥感学会主席荣誉奖

○ 刚刚！自然资源部印发技术文件，扎实推进新型基础测绘与实景三维中国建设工作

○《测绘学报》青年科学家沙龙（第1辑）报告日程

会议

○ 关于举办《测绘学报》青年科学家沙龙(第1辑)的通知

○ “第二届新型基础测绘高峰论坛”会议日程来啦！

○ 第二届“地理资源青年创新论坛”一号通知

○ Geoinformatics 2022 暨CPGIS成立30周年年会（二号通知）

《测绘学报》

○《测绘学报》同济专刊 | 龚健雅：智能遥感深度学习框架与模型设计

○ 测绘学报 | 曹闻：顾及空间分布与注记相关性的点要素注记配置算法

○ 测绘学报 | 姚翔宇：数字环境下面状要素分级设色的适宜方案分析

○ 测绘学报 | 王昊：GNSS水汽层析的自适应非均匀指数分层方法

《测绘通报》

○ 《测绘通报》2022年第3期目录

○ 《测绘通报》2022年第2期目录

○ 地市级实景三维城市建设及应用

○ 基于GIS的东川区生态环境敏感性分析

《北京测绘》

○《北京测绘》2022年第2期摘要推荐

○《北京测绘》2022年第1期摘要推荐

○《北京测绘》2021年第12期摘要推荐

○《北京测绘》2021年第11期摘要推荐

《测绘科学技术学报》

○ 摘要 |《测绘科学技术学报》2021年第2期摘要推荐

○ 摘要 |《测绘科学技术学报》2021年第4期摘要推荐

○ 摘要 |《测绘科学技术学报》2021年第5期摘要推荐

○ 摘要 |《测绘科学技术学报》2021年第6期摘要推荐

《地球信息科学学报》

○ 《地球信息科学学报》2022年第3期佳文推介

○ 龚健雅院士：全球位置信息叠加协议与位置服务网技术研究进展与展望

○ 专题征稿：地球信息科学技术在旅游休闲领域的应用

○ 专刊征稿：社会感知与地理大数据挖掘（征稿中）

《测绘工程》

○ 摘要 |《测绘工程》2022年第2期摘要推荐

○ 摘要 |《测绘工程》2022年第1期摘要推荐

○ 摘要 |《测绘工程》2021年第6期摘要推荐

○ 佳文推介 | 单目视觉技术在室内定位中的应用研究

《中国空间科学技术》

○《中国空间科学技术》2022年第2期摘要

○《中国空间科学技术》2022年第1期摘要

○《中国空间科学技术》2021年第6期摘要

○《中国空间科学技术》2021年第5期摘要推荐

《卫星应用》

○《卫星应用》2022年第02期摘要○《卫星应用》2022年第01期摘要○《卫星应用》2021年第12期摘要○《卫星应用》2021年第11期摘要

《Journal of Geodesy and Geoinformation Science》

○《测绘学报（英文版）》专刊征稿 | 地图学与地球空间信息教育：理论与实践

○ 《测绘学报（英文版）》专刊征稿 | 用于三维地理信息的摄影测量和计算机视觉

○ InSAR专刊 | 《测绘学报（英文版）》（JGGS）2022年第1期发布

○ GFZ德国地球科学研究中心-贺玉芳 | 《测绘学报（英文版）》（JGGS）InSAR专刊

○ 长安大学-朱武教授 | 《测绘学报（英文版）》（JGGS）InSAR专刊

《Satellite Navigation》

○ 杨飞博士：GNSS天顶对流层精化模型的构建与分析| SANA佳文速递

○ 2022征文| SatNav“普适定位、室内导航与基于位置服务”专题

○ 2022征文| SatNav“GNSS地学应用”专题

○ 李昕博士：多频多星座PPP-RTK原理及性能| SANA佳文速递

《自然资源遥感》

○ 摘要 |《自然资源遥感》2022年第1期摘要推荐

○ 《自然资源遥感》征稿：“海岸带空间资源及生态健康遥感监测”专栏

○ 摘要 |《自然资源遥感》2021年第3期摘要推荐

○ 摘要 |《自然资源遥感》2021年第4期摘要推荐

《Journal of Geovisualization and Spatial Analysis》

○《Journal of Geovisualization and Spatial Analysis》入驻“智绘科服”融媒体平台！

○ JGSA国际期刊2021年第5卷第2期论文摘要

○ 高被引论文推荐 | Journal of Geovisualization and Spatial Analysis

○ JGSA论文推荐 | 地理信息科学研究在过去几十年中都在关注什么？

《全球定位系统》

○《全球定位系统》入驻“智绘科服”融媒体平台！

○《全球定位系统》论文推荐 | 刘光明：我国常用地心坐标系的现状与发展

○《全球定位系统》2022年第1期目次

○《全球定位系统》2022年“卫星导航与增强”专栏征文

《导航定位与授时》

○《导航定位与授时》入驻“智绘科服”融媒体平台！