基于SINS/LDV组合的连续高程测量方法及系统技术方案

本发明专利技术涉及一种基于SINS/LDV组合的连续高程测量方法及系统，属于大地测量领域。所述系统包括捷联惯导系统、激光多普勒测速仪、UPS电源以及导航计算机，其中，捷联惯导系统用于敏感载体的角运动和线运动；捷联惯导系统分别与激光多普勒测速仪和导航计算机连接；激光多普勒测速仪与导航计算机连接，UPS电源分别与捷联惯导系统和激光多普勒测速仪连接；为了提高高程测量精度，以捷联惯导系统的脉冲输出为基准去触发激光多普勒测速仪的输出，预防捷联惯导系统与激光多普勒测速仪内部晶振不一致导致的数据丢失现象。导致的数据丢失现象。导致的数据丢失现象。

【技术实现步骤摘要】
基于SINS/LDV组合的连续高程测量方法及系统


[0001]本专利技术涉及一种基于捷联惯导系统(Strapdown Inertial Navigation System,SINS)/激光多普勒测速仪(Laser Doppler Velocimeter,LDV,通常简称测速仪)组合的连续高程测量方法及系统，属于大地测量领域。

技术介绍

[0002]高程测量在地质勘探、地形绘制、重力测量等方面具有重要作用。在野外进行作业时，通常需要进行连续高程测量，而目前的高程测量方法如几何水准测量、三角高程测量等方式无法实现连续高程测量，全球卫星导航系统(Global Navigation Satellite System,GNSS)虽然可以实现连续高程测量，但是在丛林、山谷以及高楼密集的环境中GNSS信号较差导致测量精度急剧下降。
[0003]为了实现全自主式连续高程测量，可以采用惯导系统与激光多普勒测速仪组合进行测量，惯导系统可以选用旋转惯导系统或捷联惯导系统，两套系统与激光多普勒测速仪组合的高程测量都可以避免环境的影响。其中旋转惯导系统/激光多普勒测速仪组合连续高程测量对惯导系统的陀螺和加速度计精度要求较低，同时存在以下问题：
[0004]第一：旋转惯导系统中包含转位机构，增加了惯导系统设计和安装的复杂度，同时也增加了测量系统的不稳定性；
[0005]第二：转位机构在安装时存在水平倾角，其补偿过程增加了测量方法的复杂度。
[0006]因此，可以采用包含更高精度陀螺和加速度计的捷联惯导系统/激光多普勒测速仪组合的方式进行连续高程测量，一方面降低了系统设计和测量方法的复杂度，另一方面又提高了测量方法的精度。

技术实现思路

[0007]针对当前连续高程测量方法的缺陷，本专利技术的目的是提出一种基于SINS/LDV组合的连续高程测量方法及系统。
[0008]为了实现上述技术目的，本专利技术的技术方案如下：
[0009]一种基于SINS/LDV组合的连续高程测量方法，基于由捷联惯导系统、激光多普勒测速仪、UPS电源以及导航计算机组成的连续高程测量系统，该方法分为以下步骤：
[0010]S1：将捷联惯导系统分别与激光多普勒测速仪和导航计算机连接；
[0011]将激光多普勒测速仪与导航计算机连接，同时将UPS电源分别与捷联惯导系统和激光多普勒测速仪连接；
[0012]S2：在将搭建好的系统安装于载体(例如实验车)时，由于安装误差的存在，测速仪坐标系m系与载体坐标系b系之间不能完全重合，需要对其进行标定；其中影响激光多普勒测速仪速度投影的安装误差有测速仪坐标系m系与载体坐标系b系之间的俯仰安装角误差θ、航向安装角误差ψ以及比例因子K，其标定方法如下：
[0013]S2.1：打开激光多普勒测速仪和捷联惯导系统进行数据采样；
[0014]S2.2：首先设θ＝0、ψ＝0以及K＝1，起点处地标点A的坐标为(X
A
，Y
A
，Z
A
)；
[0015]S2.3：经过1～2min直线运动后载体行驶到地标点B，其坐标为(X
B
，Y
B
，Z
B
)，点A到点B位移大小为L1；
[0016]激光多普勒测速仪经航迹推算后的位置为点C，其坐标为(X
C
，Y
C
，Z
C
)，点A到点C位移大小为L2，具体航迹推算方法参见参考文献《One
‑
dimensional reference
‑
beam LDV for accurate altitude estimation in a land vehicle》(Rong Huang,Qi Wang,Xiaoming Nie,et al,《Applied Optics》,2020.11)；
[0017]则标定的结果为：
[0018][0019]S2.4：通过标定后的俯仰安装角误差θ、航向安装角误差ψ计算得到测速仪坐标系m系和载体坐标系b系之间的安装误差矩阵
[0020][0021]通过计算的测速仪坐标系m系和载体坐标系b系之间的安装误差矩阵以及比例因子K，得到某一导航时刻激光多普勒测速仪速度在b系下的投影：
[0022][0023]式中，l为导航信息的更新迭代时刻，为某一导航时刻激光多普勒测速仪测量的速度在测速仪坐标系m系下投影，由激光多普勒测速仪提供。
[0024]S3：完成系统的安装和激光多普勒测速仪安装误差标定后，接着进入高程测量；
[0025]进入高程测量后，导航计算机对测量过程中捷联惯导系统的输出进行姿态解算，得到载体坐标系b系与导航坐标系n系之间实时的3
×
3姿态矩阵
[0026]通过实时姿态矩阵进一步得到某一导航时刻测速仪速度在导航坐标系n系下的投影
[0027][0028]惯导系统由于天向无阻尼会导致系统误差迅速发散，实时姿态矩阵的误差也会迅速增大，因此测速仪速度的误差也会增大；
[0029]为了得到稳定的高精度实时姿态矩阵必须借助激光多普勒测速仪对惯导系统的天向进行阻尼，通常选择卡尔曼滤波器进行不同传感器之间的信息融合并对系统误差
进行估计与反馈，从而抑制惯导系统误差的发散，卡尔曼滤波器的构建步骤参见参考文献《基于二维激光多普勒测速仪的车载组合导航系统》(陈红江，聂晓明，王梦成，《红外与激光工程》，2018.12)；
[0030]至此可以稳定输出高精度的实时姿态矩阵
[0031]S4:将(4)式展开，某一导航时刻激光多普勒测速仪速度在导航坐标系n系下各个方向上的投影为：
[0032][0033]上式中，v
E_LDV(l)
、v
N_LDV(l)
、v
U_LDV(l)
为某一导航时刻激光多普勒测速仪在导航坐标系n系的东、北、天三个方向上的速度投影，为某一导航时刻激光多普勒测速仪在载体坐标系b系的x、y、z轴三个方向上的速度投影；
[0034]将激光多普勒测速仪速度进行如下计算，得到任意导航时刻载体的位置信息：
[0035][0036]上式中，T为激光多普勒测速仪的速度更新周期，L
(l)
、λ
(l)
、h
(l)
分别为某一导航时刻航迹推算得到的载体的纬度、经度和高程信息，L
(l
‑
1)
、λ
(l
‑
1)
、h
(l
‑
1)
分别为上一导航时刻航迹推算得到的载体的纬度、经度和高程信息，其中起始点位置L0、λ0、h0由地标点给出，R
E
、R
N
分别为载体所在位置的卯酉圈半径和子午圈半径；
[0037]最后，经过(6)式中第三式，可以连续测得载体行驶过程中的高程信息。<本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于SINS/LDV组合的连续高程测量方法，基于由捷联惯导系统、激光多普勒测速仪、UPS电源以及导航计算机组成的连续高程测量系统，其特征在于，该方法分为以下步骤：S1：将捷联惯导系统分别与激光多普勒测速仪和导航计算机连接；将激光多普勒测速仪与导航计算机连接，同时将UPS电源分别与捷联惯导系统和激光多普勒测速仪连接；S2：在将搭建好的系统安装于载体时，由于安装误差的存在，测速仪坐标系m系与载体坐标系b系之间不能完全重合，需要对其进行标定；其中影响激光多普勒测速仪速度投影的安装误差有测速仪坐标系m系与载体坐标系b系之间的俯仰安装角误差θ、航向安装角误差ψ以及比例因子K，其标定方法如下：S2.1：打开激光多普勒测速仪和捷联惯导系统进行数据采样；S2.2：首先设θ＝0、ψ＝0以及K＝1，起点处地标点A的坐标为(X
A
，Y
A
，Z
A
)；S2.3：经过1～2min直线运动后载体行驶到地标点B，其坐标为(X
B
，Y
B
，Z
B
)，点A到点B位移大小为L1；激光多普勒测速仪经航迹推算后的位置为点C，其坐标为(X
C
，Y
C
，Z
C
)，点A到点C位移大小为L2；则标定的结果为：S2.4：通过标定后的俯仰安装角误差θ、航向安装角误差ψ计算得到测速仪坐标系m系和载体坐标系b系之间的安装误差矩阵载体坐标系b系之间的安装误差矩阵通过计算的测速仪坐标系m系和载体坐标系b系之间的安装误差矩阵以及比例因子K，得到某一导航时刻激光多普勒测速仪速度在b系下的投影：式中，l为导航信息的更新迭代时刻，为某一导航时刻激光多普勒测速仪测量的速度在测速仪坐标系m系下投影，由激光多普勒测速仪提供；S3：完成系统的安装和激光多普勒测速仪安装误差标定后，接着进入高程测量；进入高程测量后，导航计算机对测量过程中捷联惯导系统的输出进行姿态解算，得到载体坐标系b系与导航坐标系n系之间实时的3
×
3姿态矩阵
通过实时姿态矩阵进一步得到某一导航时刻测速仪速度在导航坐标系n系下的投影投影至此可以稳定输出高精度的实时姿态矩阵S4:将(4)式展开，某一导航时刻激光多普勒测速仪速度在导航坐标系n系下各个方向上的投影为：上式中，v
E_LDV(l)
、v
N_LDV(l)
、v
U_LDV(l)
为某一导航时刻激光多普勒...

【专利技术属性】
技术研发人员：魏国，杨泽坤，高春峰，周健，聂晓明，王林，于旭东，周文健，刘万青，侯承志，
申请(专利权)人：中国人民解放军国防科技大学，
类型：发明
国别省市：