一种基于一致性的分布式天基雷达组网空中目标高度估计与定位方法

王增福 邵 毅 祁登亮 金术玲

①(西北工业大学自动化学院 西安 710072)

②(西安航天动力试验技术研究所 西安 710100)

③(中国电子科技集团公司第三十八研究所 合肥 230088)

1 引言

天基探测雷达具有全天候、全天时的战略、战术预警能力，且具有不受地球曲率限制、不易受攻击等优势，在预警防御系统中具有巨大潜力。本文考虑面向大范围自主搜索发现任务的天基雷达。设计过程中，为了实现较好的最小可检测速度，通常方位向孔径尺寸较大；而综合考虑重量、收拢等因素，俯仰向口径尺寸较小，俯仰向测角能力较弱。因此，单个天基探测雷达对空中目标跟踪定位时，相比于目标径向距量测及方位角量测误差，目标俯仰角量测误差较大；仅利用单雷达径向距量测、方位角量测、俯仰角量测难以有效估计空中目标高度。空中目标高度信息的缺失严重放大了目标定位误差，同时不利于目标属性判别与威胁估计。天基探测雷达空中目标高度估计研究较为少见。Graves[1]提出了一种基于多路径干扰模式的天基雷达空中目标高度估计方法，理论上难以适用于低空目标。而通过多个天基探测雷达组网，融合利用多个天基探测雷达径向距、方位角量测，可提高空中目标高度估计准确度，对改善目标定位精度及威胁估计能力等具有重要意义。

天基探测雷达目标定位涉及多次坐标转换，具有较强的非线性，需要采用非线性滤波算法估计目标运动状态。常见的非线性滤波方法有扩展卡尔曼滤波(Extended Kalman Filter,EKF)、无迹卡尔曼滤波(Unscented Kalman Filter,UKF)、转换量测的卡尔曼滤波(Converted Measurement Kalman Filter,CMKF)和粒子滤波(Particle Filter,PF)等。Kulikov等人[2]基于全局误差控制常微分方程(Ordinary Differential Equations,ODE)求解器提出了一种精确的连续离散EKF方法，在高维雷达目标跟踪场景下对所提算法进行了验证，数值结果表明该方法可用于实际的目标跟踪，灵活度和鲁棒性更高。Kulikov等人[3]和Zhang等人[4]采用无迹变换(Unscented Transform,UT)来近似非线性系统的状态转移过程和量测方程，结果表明其滤波精度要优于EKF方法。然而，EKF和UKF滤波器忽略了状态协方差矩阵的欠估计问题，导致滤波精度下降。张连仲等人[5]采用一种基于期望最大化的转换量测卡尔曼滤波方法，将球坐标系下的雷达量测值转换到直角坐标系下，利用期望最大化方法修正了量测噪声协方差，最后利用标准卡尔曼滤波器进行滤波，提高了野值干扰下的目标跟踪精度，但该方法需要利用真值计算量测转换误差矩阵，进而导致兼容性问题。Wu等人[6]提出了一种重新加权的鲁棒粒子滤波方法来降低量测野值对目标跟踪精度的影响，结果表明该方法具有较高的跟踪精度，但是，由于使用蒙特卡罗采样，粒子滤波所需计算时间成本过高。Ait-El-Fquih等人[7]提出了一种基于变分贝叶斯(Variational Bayes,VB)的粒子滤波方法，该方法将VB方法与随机抽样技术相结合以分割状态空间，从而推导出一种新的计算效率高的粒子滤波算法。闫文旭等人[8]在VB框架下，将系统状态的后验分布近似于一个可解的变分分布，通过迭代优化证据下界，并考虑惩罚因子调整迭代的步长，在非线性基础上推导出闭合形式的迭代EKF。

近年来，概率图模型及运行其上的消息传递方法得到了极大发展，并被应用于目标跟踪定位问题中[9–11]。概率图模型结合概率论和图论，利用图模型表示变量的联合概率分布，利用条件独立性将高维推断问题分解，从而构造高效的学习和推断算法。常见图模型包括贝叶斯网络、马尔可夫随机场以及因子图，且它们之间可以相互转化。同时，对于复杂推断问题，具体的消息传递方法近年来也备受关注[12–15]，常见算法包括信度传播(Belief Propagation,BP)、环信度传播(Loopy BP,LBP)、变分消息传递(Variational Message Passing,VMP)、期望传播(Expectation Propagation,EP)、非参数信度传播(Nonparametric BP,NBP)、粒子信度传播(Particle BP,PBP)、广义信度传播(Generalized BP,GBP)、平均场近似与信度传播组合的消息传递等。

多天基雷达组网融合跟踪本质上是一个复杂的推断问题。考虑到天基平台计算资源的稀缺性等实际需求，需要准确、低计算复杂度的融合跟踪方法。鉴于因子图模型对复杂推断问题强大的表征能力，本文采用因子图对多部雷达目标跟踪问题进行建模，首先基于因子图描述传统的一致卡尔曼滤波，随后引入耦合参数，并考虑到目标跟踪模型的非线性特性，推导了分布式一致扩展卡尔曼滤波(Distributed Consistent Extended Kalman Filter,DCEKF)算法。随后，针对EKF线性化非线性模型可能导致误差变大的问题，在建立的因子图的基础上，利用分布式NBP (DCNBP)算法，有效解决了多雷达目标跟踪的非线性问题，相较于DCEKF算法跟踪精度明显提升。

2 系统建模与问题描述

2.1 天基雷达目标跟踪坐标系转换

天基雷达提供的目标测量(径向距、方位角和俯仰角)位于雷达测量坐标系，而常见的空中目标经度、纬度、高度信息在大地坐标系中描述。为了量测信息的统一，将大地坐标系向雷达测量坐标系转换。首先是大地坐标系到地心(第4)赤道坐标系(简称ECEF坐标系)的转换，本文采用近似迭代方法[16]，该方法具有更好的转换效率和精度，并且在地球半径50 km以内，不存在奇点和不收敛情况。下面介绍本文涉及的几种坐标系及坐标转换[8]。

如图1所示，ECEF坐标系以地球球心O为原点，OXE轴指向0经度方向，OZE轴垂直赤道平面并指向地球自转角速度方向，OXE,OYE与OZE轴构成右手法则坐标系。为了坐标转换需要，同样在图1中建立地心(第1)赤道坐标系(简称ECI，又称J2000.0坐标系)，其以地球球心O为原点，OXJ轴指向春分点，OZJ轴垂直赤道平面并指向地球自转角速度方向，OXJ,OYJ与OZJ轴构成右手法则坐标系。天基雷达通常安装在卫星上，故建立卫星轨道坐标系(简称LVLH坐标系)，其以卫星质心o为原点，oxL轴为地心指向卫星质心的卫星矢径方向，ozL轴为卫星运动速度方向，oxL,oyL和ozL轴构成右手法则坐标系，下面介绍从ECEF坐标系到雷达量测坐标系的转换步骤，坐标转换的流程如图2所示。

图1 坐标系示意图Fig.1 Coordinate systems

图2 坐标系转换流程图Fig.2 Procedure of coordinate systems transformation

为了得到较为精确的坐标转换，考虑地球的极移、自转、章动、岁差的影响，从ECEF到ECI的坐标变换为

其中，[A],[B],[C] 和 [D] 分别为t0时刻的极移、自转、章动和岁差矩阵，且

其中，xp和yp为地极坐标，θG为格林尼治恒星时角，εM,Δε和 Δϕ分别为平黄赤交角、交角章动和黄经章动，Zp,θp和ξp为岁差参数，t0为历元起始时刻，Rx(·),Ry(·)和Rz(·) 分别为绕x轴、y轴和z轴旋转的方向余弦矩阵。

从ECI坐标系到卫星LVLH坐标系的旋转矩阵为

其中，ν为纬度幅角，Ω和i为轨道倾角与升交点赤经。

卫星本体坐标系oxbybzb以卫星质心为坐标原点o，3轴为卫星的3个惯性主轴，oyb轴垂直于星箭分离面，ozb轴垂直指向星体对地面，oxb,oyb与ozb构成右手直角坐标系。理想情况下，卫星入轨姿态角全为 0°，卫星LVLH坐标系与卫星本体坐标系重合。卫星LVLH坐标系到卫星本体坐标系的转换为

其中，γs,βs和αs分别为卫星在本体坐标系中姿态的滚动角、俯仰角和偏航角，M1,M2,M3分别为按照右手法则绕x,y,z轴逆时针旋转一定角度的旋转矩阵，xc,yc和zc为目标在卫星LVLH坐标系下的位置，xb,yb和zb为目标在卫星本体坐标系中的位置。

如图3所示，天基雷达天线阵面坐标系oxayaza坐标原点o在天线中心，oxa轴指向天线横轴，oya轴指向天线纵轴，oza轴指向天线阵面法向，oxa,oya与oza构成右手直角坐标系。一般天线安装在卫星的对地面(通常天线阵面坐标系oxa的轴与卫星本体坐标系的oxb轴重合)，主要用于对地任务，入轨后转动角度满足对地工作要求，当所有阵面姿态角均为 0°时，天线阵面坐标系与卫星本体坐标系重合。卫星本体坐标系到天线阵面坐标系的转换为

图3 天基雷达天线阵面坐标系Fig.3 Space-based radar antenna array coordinate system

其中，γc,βc和αc分别为雷达天线阵面的滚动角、俯仰角和偏航角，xa,ya和za为目标在天线阵面坐标系下的位置，T1,T2,T3分别为按照右手法则绕x,y和z轴逆时针旋转一定角度的旋转矩阵。

雷达阵面测量坐标系为以天线中心为原点的球坐标系，从天线阵面坐标系到雷达量测坐标系的转换为

其中，rz,θz和ϕz分别为雷达阵面测量坐标系下的空中目标的径向距、方位角、俯仰角。

2.2 目标运动模型

不失一般性，考虑空中弱机动目标，此类目标的运动可以近似为匀速运动，其运动模型可以描述为

其中，xk为k时刻目标在ECEF坐标系下的位置速度矢量，，F为目标状态转移矩阵

其中，T为采样间隔，ωk-1为过程噪声，且假设ωk-1～N(0,Q)。

2.3 雷达量测模型

假设在k时刻共有N部天基雷达对目标进行测量，每部雷达得到的目标量测信息为径向距、方位角。为了方便描述多雷达组网分布式融合中每个雷达的局部目标跟踪，重写目标运动模型式(7)如下。在第n个雷达节点，目标的运动模型可以近似为

由于俯仰角误差较大，仅利用雷达径向距和方位角量测估计目标高度及位置，第n个雷达在k时刻的量测方程表示为

2.4 基于因子图的天基雷达组网目标跟踪建模

因子图是概率图模型的一种。一般情况下，高维复杂随机系统的推断问题可以通过带有局部交互函数的概率图模型来建模，即利用因子分解方式和条件独立性假设紧凑表示多个变量的联合分布。因子图模型中包含变量节点和因子节点，其中每个变量节点表示一个随机变量，每个因子节点表示一个局部函数，连接变量节点和因子节点的边表示该因子是该变量的函数。因子图模型结合消息传递方法为解决复杂环境下的目标跟踪与多源信息融合问题提供了统一框架，将复杂高维推断问题转换为优化问题，在估计精度、计算复杂度和实现的灵活性方面具有显著优势，能够为目标跟踪问题提供效果好、效率高和可扩展的统一解决方案。为此，在上述空中目标与天基雷达量测模型的基础上，本文采用因子图方法对天基雷达组网目标跟踪进行建模。

其中，p(xk |xk-1)为状态转移概率分布函数，p(yk |xk)为似然函数。式(12)可以用因子图4表示，其中，fk和hk分别表示k时刻状态转移函数和似然函数对应的因子节点。

图4 目标状态与量测联合概率分布函数的因子图分解Fig.4 Factor graph decomposition of the joint probability distribution function of target states and measurements

图5 k 时刻局部因子图模型Fig.5 A local factor graph model at time k

3 基于因子图的分布式滤波算法

3.1 分布式一致扩展卡尔曼滤波

图4所示因子图模型中，k时刻，通过因子图上的递归计算得到状态xk的后验边缘分布和状态xk+1的预测分布，即

式(13)与式(14)为目标全局状态变量的更新与预测。根据2.3节的雷达观测模型，在局部状态变量中，引入一个耦合因子节点gji来表示全局变量xk的两个复制状态变量和在相邻节点i和j上的关系，其中，i ∈Dj，Dj表示j的邻居集合。将gji定义为一个指数函数[17]

因此，k时刻局部变量的后验边缘分布为

因子节点gji到变量节点的消息为

将式(15)和式(21)代入到式(20)中，可以得到迭代方程。初始化后，第l次迭代时的方程为

至此，完成了DCEKF算法推导，具体流程总结如算法1所示。

3.2 分布式非参数信度传播算法

在DCEKF算法中，将量测方程式(10)中的非线性进行准线性化近似，以达到采用均值和方差表示消息的目的。但线性化近似可能产生较大损失，因此针对因子图模型中的非线性、非高斯概率分布，结合信度传播与粒子滤波，采用分布式非参数信度传播算法(DCNBP)计算因子图中的消息传递。在2.4节基础上加入耦合函数后，k时刻的因子图如图6所示。

算法 1 天基雷达组网目标跟踪DCEKF算法Alg.1 DCEKF algorithm for target tracking in space-based radar networks

图6 k时刻因子图模型Fig.6 The factor graph at time k

其中，g-1(·)表示从雷达阵面坐标系到ECEF坐标系的转换函数，且ϕz～U(0,2π)。因此，从因子节点传递至变量节点的消息为，其粒子化形式为

其中，N(x;xi,s,Pi)表示归一化的高斯密度函数，均值和协方差分别为xi,s和Pi。式(41)中，S个高斯密度函数选取相同的协方差

因此，D个输入消息的乘积可以表示为

易知，D个高斯密度函数的乘积仍然是一个高斯密度函数

因此，式(45)所示的高斯密度函数N(x;)的权重可以表示为

由式(40)可知，p(x)表示D个高斯混合式的乘积，每个高斯混合又含有S个高斯分量，则p(x)一共含有SD个高斯分量，直接采样的计算复杂度为O(SD) 。因此，对D个高斯混合式的乘积采用重要性采样[18]，引入辅助变量用来表示抽取样本的高斯分量

且ζ1:D={ζ1,ζ2,...,ζD}。离散随机变量ζi=li表示从第i个高斯混合式中的第li个高斯分量中抽取一个样本，则式(40)可写成

令重要性采样的建议性分布为p(ζ1:D)，且

综上，以线图拓扑下的天基雷达网络为例，基于DCNBP的天基雷达组网分布式目标融合跟踪结构图如图7所示。每个雷达节点包含利用本地量测的目标跟踪模块以及利用其他节点估计值的融合模块。初始化后，在k时刻，每部雷达在目标跟踪模块基于本地量测进行本地跟踪。本地跟踪后，每部雷达再根据相邻雷达传递的本地跟踪值与本地量测进行消息融合更新，经过多次迭代更新，实现多部雷达目标跟踪的一致性融合。

算法 2 基于重要性采样的消息乘积粒子化算法Alg.2 Message product particleization algorithm based on importance sampling

3.3 算法复杂度分析

本节通过计算执行的基本运算，比较分析DCEKF算法与DCNBP算法的计算复杂度。记目标的局部状态向量大小为s，局部量测向量大小为M，DCNBP算法中粒子数为S，消息迭代次数为L。

对于DCNBP算法，每个节点的消息均为单独计算，因此其状态向量和量测向量维度为s和L。每个时刻，首先采用式(28)和式(29)生成带权值粒子进行粒子初始化，生成粒子和权值的计算复杂度分别为s2S和s2+sS+s，从而该步每个节点的计算复杂度为O(s2S)。然后，采用算法2计算算法3中的所有消息乘积的粒子化形式。对于算法2，记节点的输入消息个数为D，计算所有高斯混合消息的均值与协方差的复杂度为sS和s2DS2。根据设置辅助变量采样粒子，式(46)和式(47)的计算复杂度分别为s3D+s2D+s2和sD+s2+2s。因此，D个输入消息乘积的粒子化形式的计算复杂度为DS(s3D+s2D+sD+2s2+2s)。从而，算法2的计算复杂度为O(s2DS2+s3D2S)。因此，经过多次迭代计算每个节点消息的计算复杂度为O(s2LDS2+s3LD2S)。最后通过算法2计算每个节点的后验分布。因此DCNBP算法的计算复杂度为O(s2LDS2+s3LD2S)。

算法 3 天基雷达组网目标跟踪DCNBP算法Alg.3 DCNBP algorithm for target tracking in space-based radar networks

本文选择3颗天基雷达测量空中目标的径向距和方位角，单颗雷达的局部量测向量维度M与消息输入个数大小D均较小。因此，DCEKF和DCNBP算法的计算复杂度分别为O(s3L)和O(s2S2L+s3LS)。可以看出，DCNBP算法的粒子数越多，计算复杂度相较DCEKF算法越大。

3.4 滤波初始状态的选取与俯仰角先验信息

在实际跟踪中，滤波器初值的选取会对跟踪性能造成极大的影响，空中目标俯仰角误差过大是天基雷达目标定位的一个关键问题。

雷达得到的俯仰角误差越大，得到的滤波器初值与真实值偏差就越大。本文利用几何法[19]先求解目标高度h，再结合地球半径re、径向距rz和方位角θz估计目标俯仰角ϕz，由雷达与目标空间位置关系可知

其中，rs为雷达到地心的距离，(xc,yc,zc)为目标在LVLH坐标系下的坐标。实际天基雷达天线较多安装在对地面，即阵面姿态角中仅滚动角变化较多，为简化推导，仅假设滚动角γc非零，则

其中，(xa,ya,za)为目标在雷达阵面坐标系下的坐标，将式(52)代入到式(51)中，得

整理得

再将式(6)代入到式(54)，整理得

通过上述方法可估计出雷达量测的俯仰角信息，并估计出目标高度作为滤波算法初值。一方面提高了滤波算法初值精度，另一方面对于单部雷达将俯仰角估计值初值加入到量测方程中，可以部分解决单部雷达因为俯仰角缺失或误差太大无法跟踪目标的问题，但目标高度估计误差仍然较大。

4 仿真与分析

本节通过仿真验证DCNBP算法的性能，与几何法、传统集中式EKF (CEKF)、DCEKF算法、集中式粒子滤波(CPF)算法进行对比，基于高度求解结果对俯仰角进行估计并将其加入到量测，仿真验证高度估计误差对目标定位精度的影响。

4.1 仿真场景设置

考虑一个空中匀速直线飞行的目标，高度为9 km，采用3颗天基探测雷达对此空中目标进行跟踪。3颗卫星轨道六根数，即轨道半长轴a，偏心率e，轨道倾角η，近地点幅角ω，升交点赤经Ω,平近点角M参数设置如下。卫星1：6978.14 km,1.369e–15,20°,0°,0°,50°；卫星2：7078.14 km,1.369e–15,20°,0°,0°,60°；卫星3：7178.14 km,1.369e–15,20°,0°,0°,70°；3颗卫星姿态角均为0°。3部天基雷达阵面姿态角均为滚动角 30°、俯仰角 0°、偏航角 0°。雷达量测噪声标准差相同，根据典型的雷达与目标参数(脉冲宽度、方位向波束宽度等)取值，设置为径向距0.09 km，方位角 0.03°。

设置迭代次数L=4 ；系统噪声协方差均为Q=diag{10-3,10-3,10-3,10-3/T,10-3/T,10-3/T}，T=15 s；仿真得到目标在ECEF坐标系下的位置，并将其转换到大地坐标系，得到目标高度估计。算法性能评估指标为空中目标高度的均方根误差，即

其中，N表示蒙特卡罗仿真次数，取N=100，和hk分别表示k时刻目标估计高度和真实高度。在CPF算法中，令粒子数S=800。

4.2 结果分析

首先令耦合参数κ=500，粒子数S=800。CEKF、DCEKF、DCNBP、几何法估计以及CPF的目标高度RMSE结果在图8中给出，从图中可以看出，几何法的目标高度估计误差最大，在1.5～3.0 km，这是因为仅采用了量测信息。集中式算法由于能够利用全局信息，一般情况下，相比于分布式估计方法，具有更高的估计精度。因此，可以看到，传统CEKF算法和CPF算法的精度最高，目标高度估计的RMSE为0.50 km和0.45 km。但实际工作时，集中式算法需要将所有量测信息与状态信息发送到一颗卫星上进行数据处理，对该单星的通信负载与计算速率要求高，并且如果该星发生故障后，可能会导致整个系统失效，从而，对于天基雷达应用，包括CEKF,CPF算法在内的传统集中式滤波器鲁棒性与灵活性较差。而DCEKF,DCNBP等分布式滤波算法能更好地克服此类问题。图8中，DCEKF算法、DCNBP算法高度估计RMSE分别为0.85 km和0.55 km。相比于DCEKF算法，DCNBP算法估计精度提升了35.3%，同时曲线波动更小，可见DCNBP算法不仅具有较高的估计精度，也具有较好的灵活性和鲁棒性。本文未实现分布式PF算法，理论上，与DCNBP算法相比，分布式PF算法具有更高的估计精度，但计算量也更大；天基平台计算资源较为稀缺，有必要在目标高度估计精度与计算量之间做出权衡。

图8 空中目标高度估计RMSEFig.8 RMSE of aerial target altitude estimation

DCNBP算法中，耦合参数κ和采样粒子数S为关键参数，能够影响算法估计性能。图9和图10分别为不同κ和S下DCNBP算法的目标高度估计结果，表1中列出了不同参数下算法估计的RMSE平均值。可以看出，随着κ和S增大，DCNBP算法高度估计的RMSE变小，精度也越高。这是由于随着κ的增加，相邻变量的相关性越强，邻居节点量测信息利用越充分，估计精度越好。而根据3.1节，当κ →∞时，局部变量与全局变量分布相同，高度估计精度并不会随着κ的增加一直提高，其最终会趋于某一定值。同时当粒子数S越多，对非线性分布的近似误差越小，最终的估计精度也更高，但S的增加也会提高算法的计算复杂度与通信负载，在应用当中需要根据实际情况进行权衡。

表1 不同参数下DCNBP算法估计的RMSE平均值(km)Tab.1 Average RMSE obtained by DCNBP with different parameters (km)

图9 空中目标高度估计RMSE(粒子数 S=400)Fig.9 RMSE of aerial target altitude estimation(particle number S=400)

图10 空中目标高度估计RMSE(粒子数 S=800)Fig.10 RMSE of aerial target altitude estimation(particle number S=800)

在天基雷达跟踪目标过程中，雷达测量误差与雷达到目标距离相关，因此分别取雷达对空中目标量测标准差为：径向距0.07 km，方位角 0.02°；径向距0.09 km，方位角 0.03°；径向距0.11 km，方位角 0.04°。图11所示为不同雷达测量误差下DCNBP算法空中目标高度估计结果。由图可知，当目标与雷达距离越近时，雷达量测误差越小，目标最终的高度估计结果越精确。

前文中推导了基于高度估计的雷达俯仰角计算方法，该计算结果可以被认为是一个雷达伪测量值，高度估计越准确，该俯仰角伪量测值误差就越小。为了进一步提高目标定位精度，将该伪测量值当作量测信息加入到滤波中，得到雷达径向距、方位角与俯仰角信息，对空中目标进行定位估计。图12为不同高度估计误差下空中目标定位结果的RMSE。可以看出，高度估计误差越大，俯仰角伪量测误差越大，目标定位结果的RMSE也就越大；当高度估计误差很大时，可等效为俯仰角信息缺失，此时，目标定位RMSE最大；随着高度估计误差的减小，俯仰角伪量测误差逐渐减小，目标定位RMSE也逐渐减小，可知对目标高度估计越准确，加入俯仰角估计后的目标位置估计精度越高。由此可得，优先对目标高度进行估计进而得到雷达俯仰角伪测量值能够弥补雷达在实际情况中俯仰角量测误差大的缺点，进一步提高目标定位精度。

图12 不同高度误差空中目标定位估计RMSEFig.12 RMSE of aerial target position estimation with different altitude errors

5 结语

本文针对天基探测雷达空中目标定位时面临的俯仰角信息缺失、量测非线性等问题，提出了一种基于一致性的分布式天基雷达组网空中目标高度估计与定位方法。首先建立了空中目标运动模型与天基雷达量测模型，并基于因子图对传感器目标跟踪进行建模。考虑到雷达量测非线性，推导了基于因子图消息传递的DCEKF算法；针对DCEKF对高非线性问题近似误差大的问题，提出了基于采样的DCNBP算法。仿真验证表明，相较于DCEKFF、几何法，DCNBP算法对目标高度估计精度更高。同时，分析了DCNBP算法参数对估计性能的影响。最后，根据得到的目标高度值估计出目标的俯仰角，并将该估计值当作伪测量估计目标位置。仿真结果显示，目标高度估计越准确，加入俯仰角估计后的目标位置估计精度越高。

利益冲突所有作者均声明不存在利益冲突

Conflict of Interests The authors declare that there is no conflict of interests