CN113052297B - 基于卷积神经网络融合ekf的拖缆姿态解算方法及系统 - Google Patents
基于卷积神经网络融合ekf的拖缆姿态解算方法及系统 Download PDFInfo
- Publication number
- CN113052297B CN113052297B CN202110253991.4A CN202110253991A CN113052297B CN 113052297 B CN113052297 B CN 113052297B CN 202110253991 A CN202110253991 A CN 202110253991A CN 113052297 B CN113052297 B CN 113052297B
- Authority
- CN
- China
- Prior art keywords
- layer
- neural network
- data
- convolution
- ekf
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000013528 artificial neural network Methods 0.000 title claims abstract description 42
- 230000004927 fusion Effects 0.000 title claims abstract description 42
- 238000004364 calculation method Methods 0.000 title claims abstract description 25
- 238000013527 convolutional neural network Methods 0.000 claims abstract description 85
- 238000001914 filtration Methods 0.000 claims abstract description 32
- 238000000034 method Methods 0.000 claims abstract description 26
- 238000005259 measurement Methods 0.000 claims abstract description 10
- 238000012549 training Methods 0.000 claims abstract description 7
- 238000011176 pooling Methods 0.000 claims description 96
- 239000011159 matrix material Substances 0.000 claims description 59
- 230000006870 function Effects 0.000 claims description 49
- 230000004913 activation Effects 0.000 claims description 36
- 210000002569 neuron Anatomy 0.000 claims description 21
- 239000013598 vector Substances 0.000 claims description 12
- 238000005070 sampling Methods 0.000 claims description 6
- 230000007704 transition Effects 0.000 claims description 3
- 238000010276 construction Methods 0.000 claims description 2
- 238000005457 optimization Methods 0.000 claims description 2
- 238000007500 overflow downdraw method Methods 0.000 abstract description 2
- 239000010410 layer Substances 0.000 description 166
- 238000010586 diagram Methods 0.000 description 7
- 230000008569 process Effects 0.000 description 6
- 238000013461 design Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 238000003062 neural network model Methods 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- RZVHIXYEVGDQDX-UHFFFAOYSA-N 9,10-anthraquinone Chemical compound C1=CC=C2C(=O)C3=CC=CC=C3C(=O)C2=C1 RZVHIXYEVGDQDX-UHFFFAOYSA-N 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000000295 complement effect Effects 0.000 description 1
- 239000012792 core layer Substances 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000008676 import Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/04—Architecture, e.g. interconnection topology
- G06N3/045—Combinations of networks
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C21/00—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
- G01C21/10—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration
- G01C21/12—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration executed aboard the object being navigated; Dead reckoning
- G01C21/16—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration executed aboard the object being navigated; Dead reckoning by integrating acceleration or speed, i.e. inertial navigation
- G01C21/18—Stabilised platforms, e.g. by gyroscope
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/08—Learning methods
-
- H—ELECTRICITY
- H03—ELECTRONIC CIRCUITRY
- H03H—IMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
- H03H21/00—Adaptive networks
- H03H21/0012—Digital adaptive filters
- H03H21/0025—Particular filtering methods
- H03H21/0029—Particular filtering methods based on statistics
- H03H21/003—KALMAN filters
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- General Physics & Mathematics (AREA)
- Biophysics (AREA)
- Mathematical Physics (AREA)
- Data Mining & Analysis (AREA)
- Evolutionary Computation (AREA)
- General Health & Medical Sciences (AREA)
- Molecular Biology (AREA)
- Computing Systems (AREA)
- General Engineering & Computer Science (AREA)
- Biomedical Technology (AREA)
- Computational Linguistics (AREA)
- Software Systems (AREA)
- Artificial Intelligence (AREA)
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- Automation & Control Theory (AREA)
- Probability & Statistics with Applications (AREA)
- Navigation (AREA)
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
Abstract
本发明公开了一种基于卷积神经网络融合EKF的拖缆姿态解算方法及系统。该融合方法包括:通过多个传感器测得原始姿态数据;将原始姿态的九轴数据解算为角度数据;构建卷积神经网络与EKF融合模型,将上述角度数据作为模型的输入;使用原始角度数据训练卷积神经网络与EKF融合模型,让模型可以根据原始数据自主调整EKF的参数,找到最优的滤波参数模型;将滤波后的姿态角度数据作为卷积神经网络与EKF融合模型的输出进行输出。本发明的融合解算方法,提高了姿态角度数据的测量精度。
Description
技术领域
本发明属于神经网络融合扩展卡尔曼滤波的技术领域,具体涉及一种基于卷积神经网络融合EKF的拖缆姿态解算方法及系统。
背景技术
扩展卡尔曼滤波(Extended Kalman Filter,EKF)是标准卡尔曼滤波在非线性情形下的一种扩展形式,它是一种高效率的递归滤波器(自回归滤波器)。EKF的基本思想是利用泰勒级数展开将非线性系统线性化,然后采用卡尔曼滤波框架对信号进行滤波,因此它是一种次优滤波。
在许多实际问题,遇到的系统模型都是非线性的,因此,扩展卡尔曼滤波器应用广泛。扩展卡尔曼滤波即先进行模型的局部线性化近似,再利用卡尔曼滤波公式得到结果。EKF在函数的非线性不是很剧烈的情形下,能够具有很不错的滤波效果。但是EKF也有它的不足之处:引入了线性化误差,对非线性强的系统,容易导致滤波结果下降。神经网络具有很强的容错性以及自学习、自组织及自适应能力,能够模拟复杂的非线性映射。神经网络的这些特性和强大的非线性处理能力,恰好弥补了EKF的不足之处。而针对神经网络与卡尔曼滤波融合方法的研究也已经成为许多工业应用领域最关键的问题之一。
相关技术中公开了一种小波神经网络与EKF的姿态解算方法,采用小波神经网络结合卡尔曼滤波来进行姿态解算,最终得到校准的姿态信息。
相关技术中公开了一种基于BP神经网络融合卡尔曼滤波算法的无人机目标跟踪方法,该控制方法使用融合BP神经网络与卡尔曼滤波算法预测被遮挡目标中心的位置坐标,发送给飞行模式下的无人机飞行控制系统实现目标跟踪。
上述相关技术具有以下缺陷:(1)扩展卡尔曼滤波器只有在系统噪声和测量噪声已知的情况下才能达到最优的结果,并且EKF方法在泰勒展开后,需要忽略二阶及以上高阶项,在复杂模型中误差相对较大;(2)基于小波神经网络融合EKF的姿态解算,姿态检测的误差较大,数据融合精度有待进一步的提高。而对于一般的神经网络算法,在实际应用中,随着传感器数量的增加,会产生由于参数量增大而引起计算量变大的问题,在神经网络的收敛性和时间特性上无法达到最优。
发明内容
本发明旨在至少解决现有技术或相关技术中存在的技术问题之一。
为此,本新型的目的在于提出一种基于卷积神经网络融合EKF的拖缆姿态解算方法。为了实现上述目的,本发明所要解决的技术问题在于提供一种基于卷积神经网络融合EKF的拖缆姿态解算方法及系统。
采用如下的方案:
一种基于卷积神经网络融合EKF的拖缆姿态解算方法,包括:
通过多个传感器测得原始姿态数据;
将原始姿态数据解算为原始角度数据;
构建卷积神经网络与EKF融合模型,将原始角度数据作为卷积神经网络与EKF融合模型的输入;
使用原始角度数据训练卷积神经网络与EKF融合的模型,让卷积神经网络根据原始角度数据自主调整EKF的滤波参数,找到最优的滤波参数,将最优的滤波参数输入至EKF滤波后输出角度数据;
将滤波后的角度数据作为卷积神经网络与EKF融合模型的输出。
进一步地,所述多个传感器包括:陀螺仪、加速度计和磁力计,所述原始姿态数据解算后的原始角度数据包括:俯仰角、翻滚角和航向角。
进一步地,所述EKF通过上一时刻的角度数据X(k)预测更新下一时刻的角度数据X(k-1);包括:
建立离散线性系统的状态方程为:
X(k)=F(X(k-1))+ΓW(k-1) (1)
式中,X(k)是系统在k时刻的状态,n维状态序列;F是n×n维状态转移矩阵;W(k)是系统的白噪声,r维状态序列;Γ为n×r维噪声驱动矩阵;
建立离散线性系统的观测方程为:
Z(k)=H·X(k)+V(k) (2)
式中,Z(k)是对应X(k)的观测信号,是m维观测序列;V(k)是观测的测量噪声,m维状态序列;H是单位矩阵,m×n维观测矩阵;
按以下步骤进行系统状态的预测、更新:
①状态一步预测:
②协方差矩阵一步预测:
P(k|k-1)=F·P(k-1|k-1)·FT+ГQГT (4)
③卡尔曼滤波增益的更新:
K(k)=P(k|k-1)·HT/(H·P(k|k-1)·HT+R) (5)
式中,K(k)为卡尔曼增益;R为观测噪声的协方差矩阵;
④状态的更新:
式中Z(k)为观测矩阵;
⑤协方的更新:
P(k|k)=(E-K(k)·H)·P(k|k-1) (7)
式中E为单位矩阵。
进一步地,构建卷积神经网络与EKF融合模型,将原始姿态数据解算的角度数据作为卷积神经网络的输入,将经过卷积层、池化层、全连接层和第一激活函数后输出的姿态数据作为所述卷积神经网络的输出,具体包括:
构建卷积神经网络的输入层,所述输入层的输入数据为所述多个传感器测得的原始姿态数据解算后的角度数据,所述原始角度数据为m×n矩阵,m为角度种类,n为时间序列采样数;
构建卷积神经网络的第一卷积层,并设置所述第一卷积层的过滤器大小和数量,通过所述第一卷积层的过滤器对所述输入层的输入数据进行卷积运算;
构建卷积神经网络的第一池化层,并设置所述第一池化层的过滤器大小和数量,采用均值池化对所述第一卷积层输出的数据进行运算;
构建卷积神经网络的第二卷积层,设置所述第二卷积层的过滤器大小和数量,通过所述第二卷积层的过滤器对所述第一池化层的输出数据进行卷积运算;
构建卷积神经网络的第二池化层,设置所述第二池化层的过滤器大小和数量,采用均值池化对所述第二卷积层输出的数据进行运算;
构建卷积神经网络的第三卷积层,设置所述第三卷积层的过滤器大小和数量,通过所述第三卷积层的过滤器对所述第二池化层的输出数据进行卷积运算;
构建卷积神经网络的第三池化层,设置所述第三池化层的过滤器大小和数量,采用均值池化对所述第三卷积层输出的数据进行运算;
构建卷积神经网络的输出层,通过所述输出层的全连接层和所述第一激活函数将所述第三池化层输出的数据转化为一维向量进行输出。
进一步地,所述第一激活函数和第二激活函数均为tanh,所述tanh的函数表达式为:
其中,z为所述第一激活函数和第二激活函数的输入;以及所述z的函数表达式为:
一种基于卷积神经网络融合EKF的拖缆姿态解算系统,包括:
采集单元,通过多个传感器测得原始姿态数据;
EKF数据构建单元,用于将原始姿态数据解算为原始角度数据;构建扩展卡尔曼滤波所需的数据,
卷积神经网络融合EKF模型单元,根据原始角度数据自主调整EKF的滤波参数,找到最优的滤波参数,将最优的滤波参数输入至EKF滤波后输出角度数据;
获取单元,用于获取激活函数;
激活函数优化单元,用于卷积神经网络融合EKF模型进行优化以输出姿态角度数据,以提高输出姿态角度数据的精度。
进一步地,构建卷积神经网络与EKF融合模型,将原始姿态数据解算的角度数据作为卷积神经网络的输入,将经过卷积层、池化层、全连接层和第一激活函数后输出的姿态数据作为所述卷积神经网络的输出,具体包括:
构建卷积神经网络的输入层,所述输入层的输入数据为所述多个传感器测得的原始姿态数据解算后的角度数据,所述原始角度数据为m×n矩阵,m为角度种类,n为时间序列采样数;
构建卷积神经网络的第一卷积层,并设置所述第一卷积层的过滤器大小和数量,通过所述第一卷积层的过滤器对所述输入层的输入数据进行卷积运算;
构建卷积神经网络的第一池化层,并设置所述第一池化层的过滤器大小和数量,采用均值池化对所述第一卷积层输出的数据进行运算;
构建卷积神经网络的第二卷积层,设置所述第二卷积层的过滤器大小和数量,通过所述第二卷积层的过滤器对所述第一池化层的输出数据进行卷积运算;
构建卷积神经网络的第二池化层,设置所述第二池化层的过滤器大小和数量,采用均值池化对所述第二卷积层输出的数据进行运算;
构建卷积神经网络的第三卷积层,设置所述第三卷积层的过滤器大小和数量,通过所述第三卷积层的过滤器对所述第二池化层的输出数据进行卷积运算;
构建卷积神经网络的第三池化层,设置所述第三池化层的过滤器大小和数量,采用均值池化对所述第三卷积层输出的数据进行运算;
构建卷积神经网络的输出层,通过所述输出层的全连接层和所述第一激活函数将所述第三池化层输出的数据转化为一维向量进行输出。
进一步地,所述第一激活函数和第二激活函数均为tanh,所述tanh的函数表达式为:
其中,z为所述第一激活函数和第二激活函数的输入;以及所述z的函数表达式为:
本发明的有益效果:
(1)本发明提供的基于卷积神经网络融合EKF的拖缆姿态解算方法,通过将卷积神经网络和扩展卡尔曼滤波器进行有效结合,并通过预设的通用核结构对激活函数优化,最终提高了姿态角度数据的测量精度;与传统的数据解算方法相比,本发明采用卷积神经网络(Convolutional Neural Networks,简写为CNN),通过对已知数据集的反复学习训练,逐步调整改变神经元连接权重,达到处理信息、模拟输入输出之间关系的目的,实现姿态数据的解算,解决了传统姿态解算存在的精度低、数据误差较大等缺点。
(2)本发明提供的基于卷积神经网络融合EKF的拖缆姿态解算方法,与传统的拖缆姿态解算方法相比,有较好的鲁棒适应性,可以根据不同的传感器数据,形成多样化的数学模型。
(3)本发明提供的基于卷积神经网络融合EKF的拖缆姿态解算方法,对拖缆姿态进行解算,其包含非线性关系,卷积层作为很好的非线性函数逼近器,可以很好的确定姿态传感器数据特征,卷积层的计算方式可以很大提高计算效率,降低了训练层数和训练难度,实现了在较少参数设计下,达到更优的神经网络模型,并且可以减少部分的计算量。
本发明的附加方面和优点将在下面的描述中变得明显,或通过本发明的实践了解到。
附图说明
图1示出了本发明的一个实施例的基于卷积神经网络融合EKF的拖缆姿态解算方法的示意流程图;
图2示出了本发明的另一个实施例的基于卷积神经网络融合EKF的拖缆姿态解算方法的总体设计示意图;
图3示出了本发明的一个实施例的卡尔曼滤波算法的基本原理;
图4示出了本发明的一个实施例的卷积神经网络的结构示意图;
图5示出了本发明的一个实施例的激活函数的示意图(a)和(b)分别示出了ReLU函数、tanh函数的激活函数的曲线示意图。
具体实施方式
为了能够更清楚地理解本发明的上述目的、特征和优点,下面结合附图和具体实施方式对本发明进行进一步的详细描述。需要说明的是,在不冲突的情况下,本申请的实施例及实施例中的特征可以相互结合。
在下面的描述中阐述了很多具体细节以便于充分理解本发明,但是,本发明还可以采用其他不同于在此描述的其他方式来实施,因此,本发明的保护范围并不限于下面公开的具体实施例的限制。
图1示出了本发明的一个实施例的基于卷积神经网络融合EKF的拖缆姿态解算方法的示意流程图。如图1所示,一种基于卷积神经网络融合EKF的拖缆姿态解算方法,包括:
步骤102,通过多个传感器测得原始姿态九轴数据;
步骤104,将原始姿态九轴数据解算为原始角度数据;
步骤106,构建卷积神经网络与EKF融合模型,将上述原始角度数据作为模型的输入;
步骤108,使用原始角度数据训练卷积神经网络与EKF融合模型,让卷积神经网络根据原始数据自主调整EKF的参数,找到最优的滤波参数输入对EKF角度数据进行滤波;
步骤110,将滤波后的姿态角度数据作为卷积神经网络与EKF融合模型的输出进行输出。
图2示出了本发明的另一个实施例的基于卷积神经网络融合EKF的拖缆姿态解算方法的总体设计示意图。如图2所示,一种基于卷积神经网络融合EKF的拖缆姿态解算方法主要内容是扩展卡尔曼滤波(EKF)与卷积神经网络融合来解算姿态数据。
本发明采用基于扩展卡尔曼滤波进行运算。EKF是一种建立系统状态方程、观测方程进行递归推算,通过上一时刻的角度数据X(k)预测更新下一时刻的角度数据X(k-1)。
建立离散线性系统的状态方程为:
X(k)=F(X(k-1))+ΓW(k-1) (1)
式中,X(k)是系统在k时刻的状态,n维状态序列;F是n×n维状态转移矩阵;W(k)是系统的白噪声,r维状态序列;Γ为n×r维噪声驱动矩阵。
建立离散线性系统的观测方程为:
Z(k)=H·X(k)+V(k) (2)
式中,Z(k)是对应X(k)的观测信号,是m维观测序列;V(k)是观测的测量噪声,m维状态序列;H是单位矩阵,m×n维观测矩阵。
按以下步骤进行系统状态的预测、更新:
A状态一步预测:
B协方差矩阵一步预测:
P(k|k-1)=F·P(k-1|k-1)·FT+ΓQΓT (4)
C卡尔曼滤波增益的更新:
K(k)=P(k|k-1)·HT/(H·P(k|k-1)·HT+R) (5)
式中,K(k)为卡尔曼增益;R为观测噪声的协方差矩阵。
D状态的更新:
式中Z(k)为观测矩阵。
E协方的更新:
P(k|k)=(E-K(k)·H)·P(k|k-1) (7)
式中E为单位矩阵。
利用扩展卡尔曼滤波,只需给定初值的X和P,就可以通过递推计算得到k时刻的状态估计值如图3所示,观测向量是陀螺仪原始数据解算出来的三轴数据。状态向量包括原始加速度计、磁力计数、陀螺仪数据结算分别得到的三轴数据,通过四元数互补滤波姿态解算的角度值,包括通过四元数法求出来的欧拉角,根据欧拉角可以求出翻滚角、俯仰角、航向角,但这里是欧拉角,还没有求出翻滚角、俯仰角以及航向角。通过四元数法解出的角度数据含有噪声,通过EKF调整参数进行噪声滤除,观测向量就是陀螺仪数据可以帮助修正噪声。
卡尔曼滤波的更新过程分为两步:分别是时间更新过程和量测更新过程。式(3)是利用k-1时刻的状态估计来预测k时刻状态的方法,式(4)对预测的结果进行了定量的描述。上述两个公式的运算仅与动态特性的信息相关。其它的式子属于量测更新过程,利用量测向量Z(k)来计算预测值的修正量。所有这些方程的目的都是合理地利用观测Z(k),来对状态值的更新进行修正。
在本实施例中,激活函数为tanh,tanh的函数表达式为:
其中,z为激活函数的输入;以及z的函数表达式为:
在本发明的一个实施例中,构建卷积神经网络,将原始姿态数据作为所述卷积神经网络的输入,将经过卷积层、池化层、全连接层和第一激活函数后输出的姿态数据作为卷积神经网络的输出,具体包括:构建卷积神经网络的输入层,输入层的输入数据为原始姿态九轴数据解算后的角度数据,角度数据为m×n矩阵,m为数据的种类数,n为时间序列采样数;构建卷积神经网络的第一卷积层,并设置第一卷积层的过滤器大小和数量,通过所述第一卷积层的过滤器对所述输入层的输入数据进行卷积运算;构建卷积神经网络的第一池化层,并设置所述第一池化层的过滤器大小和数量,采用均值池化对所述第一卷积层输出的数据进行运算;构建卷积神经网络的第二卷积层,设置所述第二卷积层的过滤器大小和数量,通过所述第二卷积层的过滤器对所述第一池化层的输出数据进行卷积运算;构建卷积神经网络的第二池化层,设置所述第二池化层的过滤器大小和数量,采用均值池化对所述第二卷积层输出的数据进行运算;构建卷积神经网络的第三卷积层,设置所述第三卷积层的过滤器大小和数量,通过所述第三卷积层的过滤器对所述第二池化层的输出数据进行卷积运算;构建卷积神经网络的第三池化层,设置所述第三池化层的过滤器大小和数量,采用均值池化对所述第三卷积层输出的数据进行运算;构建卷积神经网络的输出层,通过所述输出层的全连接层和所述第一激活函数将所述第三池化层输出的数据转化为一维向量进行输出。
在本实施例中,通过构建卷积神经网络的核心层数为6层,包括第一卷积层、第一池化层、第二卷积层、第二池化层、第三卷积层和第三池化层,再通过将原始姿态九轴数据解算后的角度数据作为卷积神经网络的输入,最后构建卷积神经网络的输出层,通过对卷积池化后的特征数据经过该输出层的全连接层和第一激活函数,将二维特征数据转化为一维的向量进行输出。
图4示出了本发明的一个实施例的卷积神经网络的结构示意图;如图2所示,通过加速度计、磁力计和陀螺仪产生原始姿态数据;通过EKF,将原始姿态数据解算为角度数据;进一步地,构建卷积神经网络(Convolutional Neural Network,简称CNN)与EKF融合的神经网络模型,通过使用卷积神经网络根据输入的原始角度数据输出EKF的调节参数,被EKF直接使用,以将所述原始角度数据作为所述神经网络的输入,将经过卷积层、池化层、全连接层和第一激活函数后输出的姿态数据作为所述神经网络的输出进行输出;最后输出包括俯仰角、横滚角和航向角的姿态角度数据,通过将卷积神经网络和EKF进行有效结合,最终提高了姿态角度数据的测量精度。
进一步地,图4示出了本发明的一个实施例的卷积神经网络的结构示意图。如图4所示,卷积神经网络的层数共有5层,具体包括:层1、层2、层3、层4和层5。层1为卷积神经网络的输入层,层2包括:第一卷积层和第一池化层,层3包括:第二卷积层和第二池化层,层4包括:第三卷积层和第三池化层,层5为卷积神经网络的输出层,该卷积神经网络的输出层包括全连接层,通过输出层的全连接层将第三池化层输出的数据转化为一维向量进行输出,最后通过输出层进行数据输出。
进一步地,图4所示的卷积神经网络过程如下:
(1)输入层:输入数据为M,具体为多个传感器产生地原始姿态数据,定义为m×n矩阵,m取值3,为原始角度数(包括翻滚角、俯仰角、航向角),n取值13,为时间序列采样数,所以输入数据大小为(3×13)矩阵。
(2)输入层到第一卷积层:过滤器C1(f1×f1×1)对输入数据M进行卷积运算。设置过滤器数量为32,每个大小为3×3。则单个过滤器卷积运算输出矩阵的元素计算表达式为:
其中,C1ij为过滤器内第i行第j个元素,Mc1ij为输入数据对应框内第i行第j个元素,n1为卷积输出矩阵的一个元素。通过一定的卷积步长s,重复公式(12),输出新的二维矩阵。
矩阵大小计算公式为:
其中,n为输入数据矩阵大小;此处为13,p为对输入数据填充(padding)大小,在卷积神经网络进行卷积或者池化操作时,操作过程中过滤器可能不能将某个方向上的数据刚好处理完,造成部分原始数据的丢失,所以需要对输入数据的矩阵边缘进行补“0”操作,但是在本网络结构中不存在上述部分数据尚未利用的情况。所以不需要进行填充操作,即此处p的值为0;f为过滤器大小,此处为3;s为卷积步长,此处为1;c1为输出矩阵大小。
所以,c1=13-3+1=11,输出矩阵大小为(11×11×32)。
(3)第一卷积层到第一池化层:由于在相邻时间段内的姿态数据是类似或近乎不变,因此采用均值池化,在缩减数据量的同时保留该时间段内主要数据特征。在此CNN结构中本池化层过滤器大小为(2×2),则输出矩阵元素的计算表达式为:
其中,P1ij为池化过滤器内第i行第j个元素,Mp1ij为输入数据对应框内第i行第j个元素。p1为池化输出矩阵的一个元素。通过一定的步长s,重复公式(12),输出新的二维矩阵。池化层通道数为32,由公式(11)可得池化后输出矩阵大小为:
[(11-2+1)×(11-2+1)],即(10×10×32)。
(4)第一池化层到第二卷积层:以第一池化层的输出作为该层输入,每个过滤器(3×3)对其进行卷积。其矩阵元素和大小计算与第一卷积层一致。本层过滤器个数为64,由公式(13)可得输出矩阵大小为:
[(10-3+1)×(10-3+1)],即(8×8×64)。
(5)第二卷积层到第二池化层:与第一池化层相同,也为均值池化,其矩阵元素和大小计算与第一池化层一致。本层过滤器个数为64,设置过滤器大小为2,步长为1;由公式(13)计算,池化后输出矩阵大小为:
[(8-3+1)×(8-3+1)],即(6×6×64)。
(6)第二池化层到第三卷积层:以第二池化层的输出作为该层输入,每个过滤器(3×3)对其进行卷积。其矩阵元素和大小计算与第一卷积层一致。本层过滤器个数为128,由公式(13)可得输出矩阵大小为:
[(6-3+1)×(6-3+1)],即(4×4×128)。
(7)第三卷积层到第三池化层:与第一池化层相同,也为均值池化,其矩阵元素和大小计算与第一池化层一致。本层过滤器个数为128,设置过滤器大小为3,步长为1;由公式(13)计算,池化后输出矩阵大小为:
[(4-3+1)×(4-3+1)],即(2×2×128)。
(8)全连接层:将上述输出的数据集(2×2×128)转化为一维向量,即(512×1),并经过第一激活函数进行输出。
进一步地,图5(a)和(b)示出了ReLU函数、tanh函数的激活函数的曲线示意图。由于本实施例中将原始姿态九轴数据解算后的角度数据作为输入时,其数据有正有负,而图5中只有tanh函数具有有正有负的函数特点,因此本实施例中的所述第一激活函数为tanh。
以上所述仅为本发明的优选实施例而已,并不用于限制本本发明,对于本领域的技术人员来说,本本发明可以有各种更改和变化。凡在本本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本本发明的保护范围之内。
Claims (2)
1.一种基于卷积神经网络融合EKF的拖缆姿态解算方法,其特征在于,包括:
通过多个传感器测得原始姿态数据;
将原始姿态数据解算为原始角度数据;
构建卷积神经网络与EKF融合模型,将原始角度数据作为卷积神经网络与EKF融合模型的输入;
使用原始角度数据训练卷积神经网络与EKF融合的模型,让卷积神经网络根据原始角度数据自主调整EKF的滤波参数,找到最优的滤波参数,将最优的滤波参数输入至EKF滤波后输出角度数据;
将滤波后的角度数据作为卷积神经网络与EKF融合模型的输出;所述多个传感器包括:陀螺仪、加速度计和磁力计,所述原始姿态数据解算后的原始角度数据包括:俯仰角、翻滚角和航向角;所述EKF通过上一时刻的角度数据X(k)预测更新下一时刻的角度数据X(k-1);包括:
建立离散线性系统的状态方程为:
X(k)=F(X(k-1))+ГW(k-1) (1)
式中,X(k)是系统在k时刻的状态,n维状态序列;F是n×n维状态转移矩阵;W(k)是系统的白噪声,r维状态序列;Γ为n×r维噪声驱动矩阵;
建立离散线性系统的观测方程为:
Z(k)=H·X(k)+V(k) (2)
式中,Z(k)是对应X(k)的观测信号,是m维观测序列;V(k)是观测的测量噪声,m维状态序列;H是单位矩阵,m×n维观测矩阵;
按以下步骤进行系统状态的预测、更新:
①状态一步预测:
②协方差矩阵一步预测:
P(k|k-1)=F·P(k-1|k-1)·FT+ΓQΓT (4)
③卡尔曼滤波增益的更新:
K(k)=P(k|k-1)·HT/(H·P(k|k-1)·HT+R) (5)
式中,K(k)为卡尔曼增益;R为观测噪声的协方差矩阵;
④状态的更新:
式中Z(k)为观测矩阵;
⑤协方的更新:
P(k|k)=(E-K(k)·H)·P(k|k-1) (7)
式中E为单位矩阵;构建卷积神经网络与EKF融合模型,将原始姿态数据解算的角度数据作为卷积神经网络的输入,将经过卷积层、池化层、全连接层和第一激活函数后输出的姿态数据作为所述卷积神经网络的输出,具体包括:
构建卷积神经网络的输入层,所述输入层的输入数据为所述多个传感器测得的原始姿态数据解算后的角度数据,所述原始角度数据为m×n矩阵,m为角度种类,n为时间序列采样数;
构建卷积神经网络的第一卷积层,并设置所述第一卷积层的过滤器大小和数量,通过所述第一卷积层的过滤器对所述输入层的输入数据进行卷积运算;
构建卷积神经网络的第一池化层,并设置所述第一池化层的过滤器大小和数量,采用均值池化对所述第一卷积层输出的数据进行运算;
构建卷积神经网络的第二卷积层,设置所述第二卷积层的过滤器大小和数量,通过所述第二卷积层的过滤器对所述第一池化层的输出数据进行卷积运算;
构建卷积神经网络的第二池化层,设置所述第二池化层的过滤器大小和数量,采用均值池化对所述第二卷积层输出的数据进行运算;
构建卷积神经网络的第三卷积层,设置所述第三卷积层的过滤器大小和数量,通过所述第三卷积层的过滤器对所述第二池化层的输出数据进行卷积运算;
构建卷积神经网络的第三池化层,设置所述第三池化层的过滤器大小和数量,采用均值池化对所述第三卷积层输出的数据进行运算;
构建卷积神经网络的输出层,通过所述输出层的全连接层和所述第一激活函数将所述第三池化层输出的数据转化为一维向量进行输出;所述第一激活函数和第二激活函数均为tanh,所述tanh的函数表达式为:
其中,z为所述第一激活函数和第二激活函数的输入;以及所述z的函数表达式为:
2.一种基于卷积神经网络融合EKF的拖缆姿态解算系统,其特征在于,包括:
采集单元,通过多个传感器测得原始姿态数据;
EKF数据构建单元,用于将原始姿态数据解算为原始角度数据;构建扩展卡尔曼滤波所需的数据,
卷积神经网络融合EKF模型单元,根据原始角度数据自主调整EKF的滤波参数,找到最优的滤波参数,将最优的滤波参数输入至EKF滤波后输出角度数据;
获取单元,用于获取激活函数;
激活函数优化单元,用于卷积神经网络融合EKF模型进行优化以输出姿态角度数据,以提高输出姿态角度数据的精度;构建卷积神经网络与EKF融合模型,将原始姿态数据解算的角度数据作为卷积神经网络的输入,将经过卷积层、池化层、全连接层和第一激活函数后输出的姿态数据作为所述卷积神经网络的输出,具体包括:
构建卷积神经网络的输入层,所述输入层的输入数据为所述多个传感器测得的原始姿态数据解算后的角度数据,所述原始角度数据为m×n矩阵,m为角度种类,n为时间序列采样数;
构建卷积神经网络的第一卷积层,并设置所述第一卷积层的过滤器大小和数量,通过所述第一卷积层的过滤器对所述输入层的输入数据进行卷积运算;
构建卷积神经网络的第一池化层,并设置所述第一池化层的过滤器大小和数量,采用均值池化对所述第一卷积层输出的数据进行运算;
构建卷积神经网络的第二卷积层,设置所述第二卷积层的过滤器大小和数量,通过所述第二卷积层的过滤器对所述第一池化层的输出数据进行卷积运算;
构建卷积神经网络的第二池化层,设置所述第二池化层的过滤器大小和数量,采用均值池化对所述第二卷积层输出的数据进行运算;
构建卷积神经网络的第三卷积层,设置所述第三卷积层的过滤器大小和数量,通过所述第三卷积层的过滤器对所述第二池化层的输出数据进行卷积运算;
构建卷积神经网络的第三池化层,设置所述第三池化层的过滤器大小和数量,采用均值池化对所述第三卷积层输出的数据进行运算;
构建卷积神经网络的输出层,通过所述输出层的全连接层和所述第一激活函数将所述第三池化层输出的数据转化为一维向量进行输出;所述第一激活函数和第二激活函数均为tanh,所述tanh的函数表达式为:
其中,z为所述第一激活函数和第二激活函数的输入;以及所述z的函数表达式为:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110253991.4A CN113052297B (zh) | 2021-03-04 | 2021-03-04 | 基于卷积神经网络融合ekf的拖缆姿态解算方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110253991.4A CN113052297B (zh) | 2021-03-04 | 2021-03-04 | 基于卷积神经网络融合ekf的拖缆姿态解算方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113052297A CN113052297A (zh) | 2021-06-29 |
CN113052297B true CN113052297B (zh) | 2022-11-22 |
Family
ID=76510568
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110253991.4A Active CN113052297B (zh) | 2021-03-04 | 2021-03-04 | 基于卷积神经网络融合ekf的拖缆姿态解算方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113052297B (zh) |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111103275A (zh) * | 2019-12-24 | 2020-05-05 | 电子科技大学 | Pat先验信息辅助的基于cnn和自适应ekf的动态fmt重建方法 |
WO2020230123A1 (en) * | 2019-05-12 | 2020-11-19 | Makesense Digital Health Technologies Ltd. | A system and a method for health and diet management and nutritional monitoring |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106525042B (zh) * | 2016-09-27 | 2019-10-18 | 哈尔滨工程大学 | 一种基于蚁群与扩展卡尔曼滤波相结合的多auv协同定位方法 |
CN108981694A (zh) * | 2018-07-18 | 2018-12-11 | 兰州交通大学 | 基于小波神经网络与ekf的姿态解算方法及系统 |
CN109409431B (zh) * | 2018-10-29 | 2020-10-09 | 吉林大学 | 基于神经网络的多传感器姿态数据融合方法及系统 |
CN110595466B (zh) * | 2019-09-18 | 2020-11-03 | 电子科技大学 | 轻量级的基于深度学习的惯性辅助视觉里程计实现方法 |
CN111027692A (zh) * | 2019-11-04 | 2020-04-17 | 中国电子科技集团公司第十四研究所 | 一种目标运动态势预测方法及装置 |
CN111693047B (zh) * | 2020-05-08 | 2022-07-05 | 中国航空工业集团公司西安航空计算技术研究所 | 一种高动态场景下的微小型无人机视觉导航方法 |
-
2021
- 2021-03-04 CN CN202110253991.4A patent/CN113052297B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2020230123A1 (en) * | 2019-05-12 | 2020-11-19 | Makesense Digital Health Technologies Ltd. | A system and a method for health and diet management and nutritional monitoring |
CN111103275A (zh) * | 2019-12-24 | 2020-05-05 | 电子科技大学 | Pat先验信息辅助的基于cnn和自适应ekf的动态fmt重建方法 |
Also Published As
Publication number | Publication date |
---|---|
CN113052297A (zh) | 2021-06-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108803349B (zh) | 非线性多智能体系统的最优一致性控制方法及系统 | |
CN110275436B (zh) | 一种多单臂机械手的rbf神经网络自适应控制方法 | |
CN110647042B (zh) | 一种基于数据驱动的机器人鲁棒学习预测控制方法 | |
CN101871782B (zh) | 基于set2fnn的gps/mems-ins组合导航系统定位误差预测方法 | |
CN110162046B (zh) | 基于事件触发型模型预测控制的无人车路径跟随方法 | |
CN107193210B (zh) | 一种非线性系统的自适应学习预设性能控制方法 | |
CN113900380B (zh) | 一种异构集群系统鲁棒输出编队跟踪控制方法及系统 | |
CN109472418A (zh) | 基于卡尔曼滤波的机动目标状态预测优化方法 | |
CN111399375A (zh) | 一种基于非线性系统的神经网络预测控制器 | |
CN111948937B (zh) | 多智能体系统的多梯度递归强化学习模糊控制方法及系统 | |
CN108168577A (zh) | 基于bp神经网络的mems陀螺随机误差补偿方法 | |
Hwang et al. | A fuzzy CMAC learning approach to image based visual servoing system | |
CN113052297B (zh) | 基于卷积神经网络融合ekf的拖缆姿态解算方法及系统 | |
CN113341696A (zh) | 一种运载火箭姿态控制参数智能整定方法 | |
CN109857127A (zh) | 训练神经网络模型以及飞行器姿态解算的方法、装置 | |
CN112631320A (zh) | 一种无人机自适应控制方法及系统 | |
CN110912535A (zh) | 一种新型无先导卡尔曼滤波方法 | |
CN114186477A (zh) | 一种基于Elman神经网络的轨道预测算法 | |
CN113848701B (zh) | 不确定性有向网络的分布式平均跟踪方法 | |
CN115344047A (zh) | 基于神经网络模型的机器人切换式预测控制轨迹跟踪方法 | |
CN112489075B (zh) | 基于特征函数的序惯式多传感器融合滤波方法 | |
Chen et al. | Neural network compensator-based robust iterative learning control scheme for mobile robots nonlinear systems with disturbances and uncertain parameters | |
CN111216146B (zh) | 一种适用于网络化机器人系统的二部一致性量化控制方法 | |
CN111582461B (zh) | 神经网络训练方法、装置、终端设备和可读存储介质 | |
CN114185273A (zh) | 一种饱和受限下分布式前置时间一致性控制器的设计方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |