用于AUV载重力仪的垂直加速度改正方法技术

本发明专利技术公开了一种用于AUV载重力仪的垂直加速度改正方法，包括：利用低通滤波消除AUV载重力仪的高频垂直扰动加速度，利用均值滤波抑制AUV载重力仪的低频垂直扰动加速度，利用博弈滤波改正AUV载重力仪的垂直扰动加速度与重力异常在频谱上有重合的噪声；该方法能够解决AUV搭载重力仪进行海洋重力测量时，由于自身震动、海浪、海流等因素导致重力仪摇摆从而附加额外的加速度，致使测量精度低的问题，能够提高海洋重力测量精度，具有较强的动态响应能力和对外界干扰的鲁棒性。

【技术实现步骤摘要】
用于AUV载重力仪的垂直加速度改正方法
本专利技术涉及重力仪的加速度改正方法
，特别涉及一种用于AUV载重力仪的垂直加速度改正方法。
技术介绍
自主式水下航行器(AutonomousUnderwaterVehicle，AUV)，如常规潜艇、核潜艇和水下自主作业机器人等已广泛应用于海洋油气源勘探和国防等领域。要实现水下长时间作业，必须依赖于自主水下导航定位系统。通过AUV搭载的相对重力仪采集的原始重力异常信号经过建模与实时滤波处理之后，需要进行改正才能得到反映AUV位置信息的实时重力异常，进而与已存储的重力数据库进行匹配得到AUV的位置信息，辅助于导航系统。AUV航行时会受到海浪、海流的作用而产生摇摆。重力测量是在AUV不断地摇晃条件下进行的，纵摇导致垂直扰动加速度，横摇导致水平扰动加速度。因此，要实现高精度的海洋重力测量，增强海洋重力测量系统动态响应能力和对外界干扰的鲁棒性，对AUV载重力仪的垂直加速度改正十分必要。传统垂直扰动加速度改正的方法利用数字滤波技术确定AUV垂直加速度用于垂直扰动加速度改正，该方法误差大，精度难以保证。
技术实现思路
本专利技术的目的在于提供一种用于AUV载重力仪的垂直加速度改正方法，该方法能够提高海洋重力测量精度。本专利技术是通过以下技术方案实现的：本专利技术的第一方面提供一种用于AUV载重力仪的垂直加速度改正方法，包括:利用低通滤波消除AUV载重力仪的高频垂直扰动加速度；利用均值滤波抑制所述AUV载重力仪的低频垂直扰动加速度；利用博弈滤波改正所述AUV载重力仪的垂直扰动加速度与重力异常在频谱上有重合的噪声。进一步地，所述利用低通滤波消除AUV载重力仪的高频垂直扰动加速度中所述的低通滤波采用巴特沃思滤波器。进一步地，所述利用低通滤波消除AUV载重力仪的高频垂直扰动加速度中，周期大于120s的扰动信号衰减不大于1dB，周期小于10s的扰动信号衰减不小于100dB。进一步地，所述利用均值滤波抑制AUV载重力仪的低频垂直扰动加速度，包括：建立AUV载重力仪在水中的垂直方向的运动方程：其中，z为垂直方向的位移，A为运行幅度，t为运动时间，T为运动周期，为初相位；AUV载重力仪在水中的垂直加速度为：AUV在一个周期内垂直加速度的平均值为：进一步地，所述高频垂直扰动加速度包括AUV运行振动造成的高频噪声。进一步地，所述高频垂直扰动加速度包括海浪和外部颠簸造成的高频噪声。进一步地，所述低频垂直扰动加速度包括AUV自身上下航行的低频噪声。进一步地，所述利用博弈滤波改正AUV载重力仪的垂直扰动加速度与重力异常在频谱上有重合的噪声包括：利用AUV博弈滤波找到适当系统状态的线性组合估计使得代价函数J′最小化，所述代价函数J′如下所示：式中，ωk为过程噪声，υk为量测噪声，Zk为系统状态的线性组合，系统状态的线性组合估计，X0为系统状态初始值，为系统状态初始值估计，P0、Qk、Rk和Sk均为基于特定问题选择的对称正定矩阵。进一步地，所述P0、Qk、Rk和Sk分别代表稳定性问题、噪声问题、发散问题、滤波精度问题选择的对称正定矩阵。本专利技术的用于AUV载重力仪的垂直加速度改正方法利用低通滤波消除AUV载重力仪的高频垂直扰动加速度，利用均值滤波抑制AUV载重力仪的低频垂直扰动加速度，利用AUV博弈滤波进一步改正AUV载重力仪的垂直扰动加速度与重力异常在频谱上有重合的噪声，该方法能够提高海洋重力测量精度。附图说明为了更清楚地说明本专利技术实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图做简单地介绍，显而易见地，下面描述中的附图是本专利技术的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。图1为本实施例提供的用于AUV载重力仪的垂直加速度改正方法的流程图；图2为受高频垂直扰动加速度干扰的重力异常仿真测量值；图3为图2重力异常仿真测量值仿真测量值在巴特沃思滤波器滤波后输出波形；图4为均值滤波误差示意图；图5为AUV载重力仪水下重力测量仿真数据在AUV博弈滤波改正后的曲线。具体实施方式为使本专利技术实施例的目的、技术方案和优点更加清楚，下面将结合本专利技术实施例中的附图，对本专利技术实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本专利技术一部分实施例，而不是全部的实施例。基于本专利技术中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本专利技术保护的范围。应当理解，本文中使用的术语“包括”和“具有”以及他们的任何变形，意图在于覆盖不排他的包含，例如，包含了一系列步骤或单元的过程、方法、系统、产品或设备不必限于清楚地列出的那些步骤或单元，而是可包括没有清楚地列出的或对于这些过程、方法、产品或设备固有的其它步骤或单元。应当理解，本文中使用的术语“和/或”仅仅是一种描述关联对象的关联关系，表示可以存在三种关系，例如，A和/或B，可以表示：单独存在A，同时存在A和B，单独存在B这三种情况。AUV航行时会受到海浪、海流的作用而产生摇摆。重力测量是在AUV不断地摇晃条件下进行的，纵摇导致垂直扰动加速度，横摇导致水平扰动加速度。传统垂直扰动加速度改正的方法采用GNSS差分技术确定垂直扰动加速度并给予补偿,但需要在地面架设GNSS参考基站，成本高且限制了作业范围。另外还有利用数字滤波技术确定AUV垂直加速度用于垂直扰动加速度改正，该方法误差大，精度难以保证。本专利技术基于以上问题提出一种用于AUV载重力仪的垂直加速度改正方法。下面结合具体实施例对本专利技术提供的用于AUV载重力仪的垂直加速度改正方法进行详细介绍。图1为本实施例提供的用于AUV载重力仪的垂直加速度改正方法的流程图，本实施例提供一种用于AUV载重力仪的垂直加速度改正方法，包括：S11、利用低通滤波消除AUV载重力仪的高频垂直扰动加速度；具体地，将高频垂直扰动加速度利用低通滤波器进行改正，其中，低通滤波是一种过滤方式，规则为低频信号能正常通过，而超过设定临界值的高频信号则被阻隔、减弱，但是阻隔、减弱的幅度则会依据不同的频率以及不同的滤波程序(目的)而改变。本实施例的低通滤波可采用巴特沃思滤波器，周期大于120s的信号(重力异常的主要成分)衰减不大于1dB，周期小于10s的信号(高频垂直扰动加速度)衰减不小于100dB，使得整体幅频特性单调递减。图2为受高频垂直扰动加速度干扰的重力异常仿真测量值，图3为图2重力异常仿真测量值仿真测量值在巴特沃思滤波器滤波后输出波形，该重力异常仿真测量值中假设高频垂直扰动加速度的周期为2π，重力异常原始信号的周期为100π；且海浪起伏引起的干扰加速度的幅值大于重力异常信号的幅值，取重力异常原始信号幅值为1mGal，搭载重力仪本文档来自技高网...

【技术保护点】
1.一种用于AUV载重力仪的垂直加速度改正方法，其特征在于，包括：/n利用低通滤波消除AUV载重力仪的高频垂直扰动加速度；/n利用均值滤波抑制所述AUV载重力仪的低频垂直扰动加速度；/n利用博弈滤波改正所述AUV载重力仪的垂直扰动加速度与重力异常在频谱上有重合的噪声。/n

【技术特征摘要】
1.一种用于AUV载重力仪的垂直加速度改正方法，其特征在于，包括：
利用低通滤波消除AUV载重力仪的高频垂直扰动加速度；
利用均值滤波抑制所述AUV载重力仪的低频垂直扰动加速度；
利用博弈滤波改正所述AUV载重力仪的垂直扰动加速度与重力异常在频谱上有重合的噪声。


2.根据权利要求1所述的用于AUV载重力仪的垂直加速度改正方法，其特征在于，
所述利用低通滤波消除AUV载重力仪的高频垂直扰动加速度中所述的低通滤波采用巴特沃思滤波器。


3.根据权利要求2所述的用于AUV载重力仪的垂直加速度改正方法，其特征在于，所述利用低通滤波消除AUV载重力仪的高频垂直扰动加速度中，周期大于120s的扰动信号衰减不大于1dB，周期小于10s的扰动信号衰减不小于100dB。


4.根据权利要求2所述的用于AUV载重力仪的垂直加速度改正方法，其特征在于，所述利用均值滤波抑制所述AUV载重力仪的低频垂直扰动加速度，包括：
建立AUV载重力仪在水中的垂直方向的运动方程：



其中，z为垂直方向的位移，A为运行幅度，t为运动时间，T为运动周期，为初相位；
AUV载重力仪在水中的垂直加速度为：



AUV在一个周期内垂直加速度的平均值为：

【专利技术属性】
技术研发人员：李超超，刘伟铭，张晓栋，徐飞，曾凯，赵毛毛，时小华，吕志斌，
申请(专利权)人：宁夏大学，
类型：发明
国别省市：宁夏;64