本发明涉及导航技术领域，特别涉及一种观测时滞系统的组合导航数据解算方法、装置及导航设备。
背景技术：
随着导航技术的不断发展，单一导航系统已难以满足用户对导航精度的要求，因此，一种将至少两种导航系统组合在一起的组合导航技术应用而生。
相关技术中，组合导航技术被应用在全球卫星导航系统/惯性导航系统(globalnavigationsatellitesystem/inertialnavigationsystem，gnss/ins)组合导航系统中。相较于单一的导航系统，gnss/ins组合导航系统在保持ins中信息全面性和自主性的同时，利用gnss克服了ins中导航误差随时间累积的问题，在导航精度和导航性能上明显优于单一导航系统。
然而在实时差分(real-timekinematic，rtk)模式下，gnss/ins组合导航系统存在观测时滞。gnss数据通过电台或网络传输时存在较大传输延迟，而ins数据却不存在传输延迟，导致同一时刻采集的gnss数据和ins数据的接收时刻存在较大差异，致使进行导航数据组合解算时，gnss数据和ins数据无法进行时间同步，严重影响组合导航的导航精度。
技术实现要素：
为了解决相关技术中rtk模式下，gnss数据通过电台或网络传输时存在较大传输延迟，而ins数据却不存在传输延迟，导致同一时刻采集的gnss数据和ins数据的接收时刻存在较大差异，致使进行导航数据组合解算时，gnss数据和ins数据无法进行时间同步，严重影响组合导航的导航精度的问题，本发明实施例提供了一种观测时滞系统的组合导航数据解算方法、装置及导航设备。所述技术方案如下：
第一方面，提供了一种观测时滞系统的组合导航数据解算方法，用于gnss/ins组合导航系统中，该方法包括：
存储k时刻的预测状态误差协方差矩阵以及ins解算得到的所述k时刻的第一导航数据，所述k时刻是gnss采样时刻，所述预测状态误差协方差矩阵用于指示所述k时刻的预测状态参数估计值的误差；
在k’时刻接收到所述k时刻采集的gnss数据时，根据所述gnss数据解算得到的第二导航数据、所述第一导航数据和所述预测误差协方差矩阵进行组合更新解算，组合更新解算结果包括所述k时刻的状态参数估计值和状态误差协方差矩阵；
根据累积状态转移矩阵和所述组合更新解算结果计算j时刻的预测状态参数估计值和预测状态误差协方差矩阵，所述j时刻为完成组合更新解算的时刻，所述累积状态转移矩阵用于指示所述k时刻到所述j时刻的状态变化；
当k-时刻的组合更新解算发生在所述k时刻与所述k’时刻之间时，根据所述k-时刻的组合更新解算结果对所述j时刻的预测状态参数估计值进行修正，所述k-时刻为所述k时刻之前的gnss采样时刻；
根据修正后的所述j时刻的预测状态参数估计值和所述j时刻的预测状态误差协方差矩阵修正所述j时刻的ins误差。
该实现方式中，导航设备在gnss数据采集时刻对ins解算的ins数据进行存储，并在接收到gnss采集的gnss数据时对两者进行组合更新解算，从而将组合更新解算结果转换至解算完成时刻，实现gnss数据和ins数据的时间同步；同时，当k-时刻的组合更新解算发生在k时刻与k’时刻之间，即gnss数据的传输延迟大于gnss的单位采样间隔时，导航设备考虑k-时刻的组合更新解算结果对状态参数估计值的影响，对解算完成时刻的预测状态参数估计值进行修正，提高了解算完成时刻的预测状态参数估计值的准确度，进而提高了组合导航的导航精度。
在一种可能的实现方式中，所述根据所述k-时刻的组合更新解算结果对所述j时刻的预测状态参数估计值进行修正，包括：
根据所述k时刻的状态参数估计值计算k-’时刻的第一预测状态参数估计值，所述k-’时刻为获取到所述k-时刻所采集的gnss数据的时刻；
根据所述k-时刻的组合更新解算结果中包含的状态参数估计值计算所述k-’时刻的第二预测状态参数估计值；
根据所述第一预测状态参数估计值和所述第二预测状态参数估计值对所述j时刻的预测状态参数估计值进行修正。
在一种可能的实现方式中，所述根据所述第一预测状态参数估计值和所述第二预测状态参数估计值对所述j时刻的预测状态参数估计值进行修正，包括：
将所述第一预测状态参数估计值和所述第二预测状态参数估计值的差值确定为修正值；
根据所述修正值对所述j时刻的预测状态参数估计值进行修正；
其中，所述修
该实现方式中，通过将k时刻的状态参数估计值和k-时刻的状态参数估计值状态转移至同一时刻，得到对应同一时刻的两个预测状态参数估计值，从而根据两个预测状态参数估计值的差值确定出k-时刻的组合更新解算结果对状态参数估计值所产生的影响，并应用该差值对j时刻的预计状态参数估计值进行修正，提高了解算完成时刻的预测状态参数估计值的准确度，进而提高了组合导航的导航精度。
在一种可能的实现方式中，所述根据所述gnss数据解算得到的第二导航数据、所述第一导航数据和所述预测误差协方差矩阵进行组合更新解算，包括：
根据第一公式
根据第二公式
其中，
在一种可能的实现方式中，所述根据累积状态转移矩阵和所述组合更新解算结果计算j时刻的预测状态参数估计值和预测状态误差协方差矩阵，包括：
根据第三公式
根据第四公式
其中，
第二方面，提供了一种观测时滞系统的组合导航数据解算装置，该装置包括至少一个单元，该至少一个单元用于实现上述第一方面或第一方面任意一种可能的实现方式所提供的观测时滞系统的组合导航数据解算方法。
第三方面，本发明实施例提供了一种导航设备，该导航设备包括gnss/ins组合导航系统、处理器以及存储器，所述存储器用于存储一个或一个以上的指令，所述指令被指示由所述处理器执行，所述处理器用于实现上述第一方面或第一方面中任意一种可能的实现方式所提供的观测时滞系统的组合导航数据解算方法。
第四方面，本发明实施例提供一种计算机可读存储介质，该计算机可读存储介质中存储有用于实现上述第一方面或第一方面中任意一种可能的实现方式所提供的观测时滞系统的组合导航数据解算方法的可执行程序。
附图说明
为了更清楚地说明本发明实施例中的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍。
图1示出了本发明各个实施例提供的观测时滞系统的组合导航数据解算方法的原理示意图；
图2a示出了本发明一个实施例提供的观测时滞系统的组合导航数据解算方法的流程图；
图2b是图2a所示组合导航数据解算方法的实施示意图；
图3a示出了本发明另一个实施例提供的观测时滞系统的组合导航数据解算方法的流程图；
图3b是图3a所示组合导航数据解算方法的实施示意图；
图4示出了本发明一个实施例提供的导航设备的结构框图；
图5示出了本发明一个实施例提供的观测时滞系统的组合导航数据解算装置的框图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚，下面将结合附图对本发明实施方式作进一步地详细描述。
ins是一种利用角速度传感器和加速度传感器进行导航的导航系统。工作状态下，ins对加速度传感器采集到的加速度数据进行积分，确定物体当前的行进速度；对角速度传感器采集到的角速度数据进行积分，确定物体当前的姿态，从而根据实时的行进速度和姿态绘制出物体的行进轨迹，并根据该行进轨迹进行导航。相较于其他类型导航系统，ins采集数据的频率较高(通常为100hz或200hz)，所采集的数据全面且具有自主性；但是在进行长时间导航的情况下，连续积分所产生的累计误差将影响ins的导航准确度。而gnss利用卫星进行导航时不存在累计误差，因此，将gnss和ins进行组合，并利用gnss采集的gnss数据对同一时刻ins的ins数据进行修正(即组合更新解算过程)，能够消除ins长时间导航时产生的累计误差。
但是在rtk模式下，gnss需要通过电台或网络接收gnss基站差分信号(即接收gnss数据)，而gnss基站差分信号在传输过程中存在传输延迟，且该传输延迟与gnss基站差分信号的传输距离有关。当gnss基站差分信号的传输距离较长时，传输延迟将大于gnss的单位采样时间间隔(即传输延迟超过1s)，导致利用gnss数据对ins数据进行修正时产生严重的时间异步问题，影响组合导航的导航精度。
而本发明各个实施例提供的导航数据解算方法中，如图1所示，在gnss112的采样时刻，ins111通过惯性导航单元(inertialmeasurementunit，imu)采集的imu数据经过imu数据解算后存储在存储器120中(通常为寄存器)。当接收到采集的gnss数据时，处理器即对存储器120中存储的ins数据(imu数据经过解算后得到)以及gnss数据(经过gnss数据解算得到)进行组合更新解算，实现了gnss数据与ins数据的同步化。完成组合更新解算后，处理器对组合更新解算结果进行状态转移，并在gnss数据的传输延迟较长的情况下，利用上一次组合更新解算结果对状态转移后的结果进行修正，最终输出最优的gnss/ins实时组合导航结果。在ins修正阶段，处理器即利用解算完成时刻的预测状态参数估计值对ins111进行修正，从而消除ins111在长时间导航过程中的累计误差和传感器误差。
本发明各个实施例适用于安装有gnss/ins组合导航系统的导航设备中，该导航设备可以是具有导航功能的移动终端或安装在车辆中的车载导航设备。gnss/ins中的gnss可以是全球定位系统(globalpositioningsystem，gps)、格洛纳斯导航系统(glonass)、伽利略导航系统(galileo)或北斗导航系统，且ins由多种imu构成，该imu包括加速度传感器和角速度传感器等等。下面采用示意性的实施例进行说明。
请参考图2a，其示出了本发明一个实施例提供的观测时滞系统的组合导航数据解算方法的流程图。本实施例以该组合导航数据解算方法用于安装有gnss/ins组合导航系统的导航设备中来举例说明，该方法包括：
步骤201，存储k时刻的预测状态误差协方差矩阵以及ins解算得到的k时刻的第一导航数据，k时刻是gnss采样时刻。
可选的，当在k时刻检测到采集gnss数据时，ins对同为k时刻采集的ins数据进行解算，并对解算得到的第一导航数据进行存储。其中，解算的方法为机械编排解算，第一导航数据用于指示k时刻的行进速度、行进姿态以及行进轨迹。
在存储第一导航数据的同时，导航设备对k时刻对应的预测误差协方差矩阵进行存储，该预测状态误差协方差矩阵根据ins上一采样时刻量测后的状态误差协方差矩阵迭代计算得到，用于指示k时刻的预测状态参数估计值的误差。
比如，如图2b所示，gnss的采样间隔大于ins的采样间隔，gnss在t0、t1和t2时刻进行数据采样时，导航设备存储t0、t1和t2时刻ins采样并解算得到的第一导航数据。
步骤202，在k’时刻接收在k时刻采集的gnss数据时，根据gnss数据解算得到的第二导航数据、第一导航数据和预测误差协方差矩阵进行组合更新解算。
由于gnss数据在传输过程中存在传输延迟，因此导航设备需要经过一段时间后才能获取k时刻采集的gnss数据。比如，如图2b所示，导航设备在t0’时刻获取gnss在t0时刻采集到的gnss数据，在t1’时刻获取gnss在t1时刻采集到的gnss数据。
获取到gnss数据后，导航设备对该gnss数据进行解算，得到对应的第二导航数据，并进一步根据存储的k时刻的第一导航数据、预测状态误差协方差矩阵以及该第二导航数据进行组合更新解算，即利用k时刻的gnss数据进行修正。其中，组合更新解算结果包括k时刻的状态参数估计值和状态误差协方差矩阵。
步骤203，根据累积状态转移矩阵和组合更新解算结果计算j时刻的预测状态参数估计值和预测状态误差协方差矩阵，j时刻为完成组合更新解算的时刻。
在j时刻，导航设备完成对k时刻的组合更新解算，并根据k时刻与j时刻之间各个采集时刻的状态转移矩阵迭代计算出k时刻到j时刻的累积状态转移矩阵。
进一步的，导航设备利用该累积状态转移矩阵对k时刻的组合更新解算结果进行状态转移，从而得到j时刻的预测状态参数估计值和预测状态误差协方差矩阵。
步骤204，当k-时刻的组合更新解算发生在k时刻与k’时刻之间时，根据k-时刻的组合更新解算结果对j时刻的预测状态参数估计值进行修正，k-时刻为k时刻之前的gnss采样时刻。
当gnss数据的传输延迟较短时，gnss数据的接收以及组合更新解算过程在一个gnss采样间隔内完成；但是当gnss数据的传输延迟较长时，gnss数据的接收以及组合更新解算过程可能会跨历多个gnss采样间隔。
比如，如图2b所示，导航设备在t1时刻之后的t0’时刻接收到t0时刻采集的gnss数据，导航设备在t2时刻之后的t1’时刻接收到t1时刻采集的gnss数据。
当gnss的采样时刻(k时刻)到gnss数据接收时刻(k’时刻)之间发生了组合更新解算(k-时刻的组合更新解算)时，导航设备在进行状态转移过程中，需要考虑组合更新解算对j时刻的预测状态参数估计值的影响，并进一步对j时刻的预测状态参数估计值进行修正，从而提高j时刻的预测状态参数估计值的准确性。比如，如图2b所示，导航设备使用t0时刻的组合更新解算结果对t1’时刻的预测状态参数估计值进行修正。
步骤205，根据修正后的j时刻的预测状态参数估计值和j时刻的预测状态误差协方差矩阵修正j时刻的ins误差。
通过上述步骤201至204完成gnss采集时刻到当前时刻(即j时刻)的状态转移后，导航设备进一步根据当前时刻的预测状态参数估计值和预测状态误差协方差矩阵修正j时刻的ins误差，从而消除长时间导航过程中产生的累计误差和传感器误差。
综上所述，本实施例中，导航设备在gnss数据采集时刻对ins解算的ins数据进行存储，并在接收到gnss采集的gnss数据时对两者进行组合更新解算，从而将组合更新解算结果转换至解算完成时刻，实现gnss数据和ins数据的时间同步；同时，当k-时刻的组合更新解算发生在k时刻与k’时刻之间，即gnss数据的传输延迟大于gnss的单位采样间隔时，导航设备考虑k-时刻的组合更新解算结果对状态参数估计值的影响，对解算完成时刻的预测状态参数估计值进行修正，提高了解算完成时刻的预测状态参数估计值的准确度，进而提高了组合导航的导航精度。
请参考图3a，其示出了本发明另一个实施例提供的观测时滞系统的组合导航数据解算方法的流程图。本实施例以该组合导航数据解算方法用于安装有gnss/ins组合导航系统的导航设备中来举例说明，该方法包括：
步骤301，存储k时刻的预测状态误差协方差矩阵以及ins解算得到的k时刻的第一导航数据，k时刻是gnss采样时刻。
在ins中，当前采样时刻的状态与上一采样时刻的状态有关，且当前采样时刻的状态可以根据上一采样时刻的状态和状态转移矩阵计算得到。其中，根据上一采样时刻的状态和状态转移矩阵计算当前采样时刻的状态的公式如下。
其中，
需要说明的是，k-1时刻和k时刻为相邻的ins采集时刻。
相似的，导航设备可以根据k-1时刻的状态误差协方差矩阵计算得到k时刻的预测状态误差协方差矩阵并进行存储。其中，根据计算k时刻的预测状态误差协方差矩阵的公式如下：
其中，
需要说明的是，导航设备在存储第一导航数据和预测状态误差协方差矩阵的同时，从k时刻开始累计状态转移矩阵和系统噪声协方差矩阵，供后续进行状态转移时使用。
步骤302，在k’时刻接收到k时刻采集的gnss数据时，根据gnss数据解算得到的第二导航数据、第一导航数据和预测误差协方差矩阵进行组合更新解算。
当在k’时刻接收到k时刻采集的gnss数据时，导航设备即通过gnss对接收到的gnss数据进行解算，得到指示k时刻行进方向和行进速度的第二导航数据。
完成gnss数据解算后，导航设备获取上述步骤301存储的k时刻的第一导航数据以及预测状态误差协方差矩阵，并利用解算得到的第二导航数据对其进行修正，即根据第二导航数据、第一导航数据和预测误差协方差矩阵进行组合更新解算(卡尔曼量测更新过程)，得到k时刻量测更新后的状态参数估计值以及量测更新后的状态误差协方差矩阵。
在一种可能的实施方式中，本步骤可以包括如下步骤。
一、根据第一公式
其中，
可选的，在计算k时刻的状态参数估计值和状态误差协方差矩阵前，导航设备根据存储的k时刻的预计状态误差协方差矩阵和观测矩阵计算卡尔曼滤波增益矩阵。
其中，卡尔曼滤波增益矩阵
计算得到卡尔曼滤波增益矩阵后，导航设备即通过第一公式计算量测更新后k时刻的状态参数估计值。
二、根据第二公式
在计算k时刻状态参数估计值的同时，导航设备根据在k时刻预测状态误差协方差矩阵的基础上，通过第二公式计算k时刻的状态误差协方差矩阵。其中，第二公式中，ι为单位矩阵，(ι-kkhk)t为ι-kkhk的转置矩阵，rk为k时刻的观测噪声方差矩阵，
需要说明的是，根据第二公式计算得到的状态误差协方差矩阵
通过上述步骤，导航设备完成组合更新解算后，即获得k时刻的状态参数估计值以及相应的状态误差协方差矩阵。
步骤303，根据累积状态转移矩阵和组合更新解算结果计算j时刻的预测状态参数估计值和预测状态误差协方差矩阵，j时刻为完成组合更新解算的时刻。
经过上述步骤302可以得到k时刻的状态参数估计值以及相应的状态误差协方差矩阵，为了进一步根据k时刻的状态得到j时刻(即组合更新解算完成时刻)的状态，导航设备需要进一步根据k时刻与j时刻之间的状态转移矩阵对k时刻的状态进行状态转移。
导航设备在进行状态转移前，需要计算k时刻到j时刻的累积状态转移矩阵。在一种可能的实施方式中，导航设备将k时刻至j时刻之间的状态转移矩阵进行迭代，即可计算得到k时刻至j时刻的累积状态转移矩阵φj/k。
进一步的，导航设备根据累计状态转移矩阵从k时刻迭代至j时刻，可以得到
其中，φj/k为累积状态转移矩阵，
相应的，导航设备根据k时刻的状态参数估计值以及k时刻至j时刻的累积状态转移矩阵，通过第三公式即可计算得到j时刻的预测状态参数估计值；
导航设备根据k时刻的状态误差协方差矩阵以及k时刻至j时刻的累积状态转移矩阵，通过第四公式即可计算得到j时刻的预测状态误差协方差矩阵。
步骤304，根据k时刻的状态参数估计值计算k-’时刻的第一预测状态参数估计值，k-’时刻为获取到k-时刻所采集的gnss数据的时刻。
当gnss数据的传输延迟较长时，gnss数据的接收以及组合更新解算过程可能会跨历多个gnss采样间隔。相应的，当gnss数据的接收以及组合更新解算过程跨历多个gnss采样间隔时，k-时刻的组合更新解算即发生在k时刻与k’时刻之间。
比如，如图3b所示，导航设备在t2时刻之后的t1’时刻接收到t1时刻采集的gnss数据(即t1时刻gnss数据的接收以及组合更新解算过程跨历两个gnss采样间隔)，且t0时刻的组合更新解算发生在t1时刻与t1’时刻之间。
考虑到k-时刻的组合更新解算对j时刻带来的影响，导航设备需要根据k-时刻的组合更新解算对j时刻的预测状态参数估计值进行修正，从而提高j时刻预测状态参数估计值的准确性。
在一种可能的实施方式中，如图3b所示，导航设备可以根据转移区间t1和转移区间t2内状态参数估计值的变换情况确定t0时刻的组合更新解算产生的影响。
如图3b所示，以k时刻为t1时刻，k’时刻为t1’时刻，k-时刻为t0时刻，k-’时刻为t0’时刻为例，导航设备可以根据t1时刻的状态参数估计值以及t0’时刻至t1时刻的状态转移矩阵计算t0’时刻的第一预测状态参数估计值。
具体的，第一预测状态参数估计值
步骤305，根据k-时刻的组合更新解算结果中包含的状态参数估计值计算k-’时刻的第二预测状态参数估计值。
与上述步骤304相似的，导航设备可以根据t0时刻的状态参数估计值以及t0时刻至t0’时刻的状态转移矩阵计算t0’时刻的第二预测状态参数估计值。
具体的，第二预测状态参数估计值
步骤306，根据第一预测状态参数估计值和第二预测状态参数估计值对j时刻的预测状态参数估计值进行修正。
通过上述步骤304和305，导航设备根据k时刻的状态参数估计值和k-时刻的状态参数估计值得到同一时刻(k-’时刻)的两个预测结果；进一步的，通过分析两个预测结果的差异，导航设备即可对j时刻的预测状态参数估计值进行修正。
在一种可能的实施方式中，导航设备可以将第一预测状态参数估计值和第二预测状态参数估计值的差值作为j时刻的修正值，并利用该修正值进行修正。可选的，本步骤可以包括如下步骤。
一、将第一预测状态参数估计值和第二预测状态参数估计值的差值确定为修正值。
结合上述步骤304和步骤305，导航设备计算得到的修
二、根据修正值对j时刻的预测状态参数估计值进行修正。
在一种可能的实施方式中，导航设备可以根据上述步骤计算得到的修正值以及t0’时刻至t1’时刻的状态转移矩阵进行状态参数修正。比如，导航设备可以通过如下公式对j时刻的预测状态参数估计值进行修正：
其中，
步骤307，根据修正后的j时刻的预测状态参数估计值和j时刻的预测状态误差协方差矩阵修正j时刻的ins误差。
通过上述步骤301至306完成gnss采集时刻到当前时刻(即j时刻)的状态转移以及状态参数修正后，导航设备进一步根据当前时刻的预测状态参数估计值和预测状态误差协方差矩阵修正j时刻的ins误差，从而消除长时间导航过程中产生的累计误差和传感器误差。
本实施例中，通过将k时刻的状态参数估计值和k-时刻的状态参数估计值状态转移至同一时刻，得到对应同一时刻的两个预测状态参数估计值，从而根据两个预测状态参数估计值的差值确定出k-时刻的组合更新解算结果对状态参数估计值所产生的影响，并应用该差值对j时刻的预计状态参数估计值进行修正，提高了解算完成时刻的预测状态参数估计值的准确度，进而提高了组合导航的导航精度。
以下为本发明实施例的装置实施例，对于装置实施例中未详细阐述的部分，可以参考上述方法实施例中的细节。
请参考图4，其示出了本发明一个实施例提供的观测时滞系统的组合导航设备的结构框图。该导航设备包括：gnss/ins组合导航系统410、处理器420以及存储器430。
处理器420通过总线分别与gnss/ins组合导航系统410和存储器430相连。
处理器420包括一个或一个以上处理核心。处理器420通过运行操作系统或应用程序模块，实现导航数据的解算及修正功能。
可选地，存储器430可存储操作系统432、至少一个功能所需的应用程序模块434。可选的，应用程序模块434包括存储模块434a、解算模块434b、转换模块434c和修正模块434d。其中，存储模块434a用于实现有关数据存储的步骤；解算模块434b用于实现有关导航数据解算的步骤；转换模块434c用于实现有关状态转移的步骤；修正模块434d用于实现有关数据修正的步骤。
此外，存储器430可以由任何类型的易失性或非易失性存储设备或者它们的组合实现，如静态随机存取存储器(sram)，电可擦除可编程只读存储器(eeprom)，可擦除可编程只读存储器(eprom)，可编程只读存储器(prom)，只读存储器(rom)，磁存储器，快闪存储器，磁盘或光盘。
本领域技术人员可以理解，图4中所示出的结构并不构成对导航设备的限定，导航设备可以包括比图示更多或更少的部件或组合某些部件，或者不同的部件布置。
请参考图5，其示出了本发明一个实施例提供的观测时滞系统的组合导航数据解算装置的框图。该装置可以通过软件或硬件实现成为导航设备的全部或一部分。该装置包括：存储单元510、解算单元520、转换单元530和修正单元540。
存储单元510用于实现上述图2a或图3a实施例中的步骤201和301等步骤，以及其它隐含的由导航设备执行的存储步骤；
解算单元520用于实现上述图2a或图3a实施例中的步骤202和302等步骤，以及其它隐含的由导航设备执行的导航数据解算步骤；
转换单元530用于实现上述图2a或图3a实施例中的步骤203和303等步骤，以及其它隐含的由导航设备执行的状态转移步骤；
修正单元540用于实现上述图2a或图3a实施例中的步骤204、205、304、305、306和307等步骤，以及其它隐含的由导航设备执行的数据修正步骤。
其中，存储单元510可以由图4中的处理器420执行存储器430中的存储模块434a实现；解算单元520可以由图4中的处理器420执行存储器430中的解算模块434b实现；转换单元530可以由图4中的处理器420执行存储器430中的转换模块434c实现；修正单元540可以由图4中的处理器420执行存储器430中的修正模块434d实现。
本领域普通技术人员可以理解实现上述实施例的全部或部分步骤可以通过硬件来完成，也可以通过程序来指令相关的硬件完成，所述的程序可以存储于一种计算机可读存储介质中，上述提到的存储介质可以是只读存储器，磁盘或光盘等。
以上所述仅为本发明的较佳实施例，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。