CN104048676A - 基于改进粒子滤波的mems陀螺随机误差补偿方法 - Google Patents
基于改进粒子滤波的mems陀螺随机误差补偿方法 Download PDFInfo
- Publication number
- CN104048676A CN104048676A CN201410298926.3A CN201410298926A CN104048676A CN 104048676 A CN104048676 A CN 104048676A CN 201410298926 A CN201410298926 A CN 201410298926A CN 104048676 A CN104048676 A CN 104048676A
- Authority
- CN
- China
- Prior art keywords
- particle
- mems
- weights
- particle filter
- data
- 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.)
- Granted
Links
- 239000002245 particle Substances 0.000 title claims abstract description 101
- 238000000034 method Methods 0.000 title claims abstract description 35
- 238000005070 sampling Methods 0.000 claims abstract description 12
- 230000006870 function Effects 0.000 claims description 15
- 238000001914 filtration Methods 0.000 claims description 14
- 239000011159 matrix material Substances 0.000 claims description 14
- 238000005259 measurement Methods 0.000 claims description 9
- 238000012731 temporal analysis Methods 0.000 claims description 9
- 238000000700 time series analysis Methods 0.000 claims description 9
- 241001123248 Arma Species 0.000 claims description 7
- 238000012952 Resampling Methods 0.000 claims description 6
- 238000007689 inspection Methods 0.000 claims description 6
- 238000005457 optimization Methods 0.000 claims description 4
- 238000012545 processing Methods 0.000 claims description 4
- 238000012549 training Methods 0.000 claims description 4
- 210000000262 cochlear duct Anatomy 0.000 claims description 2
- 230000000694 effects Effects 0.000 abstract description 6
- 238000009826 distribution Methods 0.000 abstract description 3
- 238000012300 Sequence Analysis Methods 0.000 abstract 1
- 238000007781 pre-processing Methods 0.000 abstract 1
- 230000008569 process Effects 0.000 description 9
- 238000013528 artificial neural network Methods 0.000 description 5
- 238000011160 research Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 3
- 238000010923 batch production Methods 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 238000000342 Monte Carlo simulation Methods 0.000 description 1
- XUIMIQQOPSSXEZ-UHFFFAOYSA-N Silicon Chemical compound [Si] XUIMIQQOPSSXEZ-UHFFFAOYSA-N 0.000 description 1
- 239000000853 adhesive Substances 0.000 description 1
- 230000001070 adhesive effect Effects 0.000 description 1
- 238000005311 autocorrelation function Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005265 energy consumption Methods 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000010365 information processing Effects 0.000 description 1
- 230000002401 inhibitory effect Effects 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000002715 modification method Methods 0.000 description 1
- 238000003825 pressing Methods 0.000 description 1
- 230000000644 propagated effect Effects 0.000 description 1
- 229910052710 silicon Inorganic materials 0.000 description 1
- 239000010703 silicon Substances 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000013179 statistical model Methods 0.000 description 1
- 238000005309 stochastic process Methods 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C25/00—Manufacturing, calibrating, cleaning, or repairing instruments or devices referred to in the other groups of this subclass
- G01C25/005—Manufacturing, calibrating, cleaning, or repairing instruments or devices referred to in the other groups of this subclass initial alignment, calibration or starting-up of inertial devices
Landscapes
- Engineering & Computer Science (AREA)
- Manufacturing & Machinery (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Gyroscopes (AREA)
Abstract
本发明属于陀螺导航领域,具体涉及一种基于改进粒子滤波的MEMS陀螺随机误差补偿方法。包括:数据的采集、预处理与检验:对处理后的MEMS陀螺的输出数据进行时间序列分析建模:利用小波网络算法优化粒子滤波,改善粒子滤波性能;对时间序列模型进行粒子滤波,对MEMS陀螺仪误差补偿,提高精度。M本发明采用一种改进的粒子滤波算法,将小波网络算法与粒子滤波典型采样算法结合,增大位于概率分布尾部的粒子权值,使较高权值的粒子分解为若干个较小权值的粒子,提高粒子样本的多样性,减小误差,提高滤波效果。将基于小波网络的粒子滤波算法应用到MEMS陀螺仪的误差补偿中,可以有效减少随机误差,提高MEMS陀螺仪精度。
Description
技术领域
本发明属于陀螺导航领域,具体涉及一种基于改进粒子滤波的MEMS陀螺随机误差补偿方法。
背景技术
早在20世纪50年底,粒子滤波在统计学和理论物理领域得到了广泛的应用。粒子滤波是一种基于蒙特卡罗方法和递推贝叶斯的统计滤波方法,它对系统的过程噪声和量测噪声没有任何限制,突破了Kalman滤波理论的框架。粒子滤波通过预测和更新来自于系统概率密度函数的采样样本,来近似非线性系统的贝叶斯估计,是现代信号与信息处理科学、统计模拟理论之间的交叉科学,具有重要的研究意义和现实价值。但是随着多次递归更新后粒子样本的多样性丧失,使得滤波效果受到影响,将小波网络算法与粒子滤波典型采样算法结合,利用小波网络调整粒子的权值,提高粒子样本的多样性,改善滤波效果。
随着微电子技术的发展,MEMS(Micro Electro Mechanical System,MEMS)系统得到迅速的发展。由于其独特的制造工艺及微型化、易于集成化、易于批量生产等特点,使其在导航制导、汽车、飞行器、机器人、生物医学、通信等领域得到广泛的应用。目前,国内外很多研究机构和高校都在致力于MEMS陀螺的研究,成为各国研究的重点内容。MEMS陀螺采用类集成电路的硅加工工艺,器件尺寸小,且重量轻,适合批量生产;性能稳定且抗干扰能力强;可靠性也比较高且易集成、功耗低。目前,基于MEMS技术的惯性器件得到了广泛的应用,但是MEMS惯性器件在高精度的惯性传感器中的应用却极大地受到了限制。但是由于MEMS陀螺精度比传统陀螺低,使其应用受到限制。如何改善MEMS陀螺精度成为MEMS陀螺的研究重点。
由于MEMS惯性器件的精度受到误差、漂移的影响,而MEMS陀螺易受制造工艺及使用环境的影响,使得MEMS陀螺的误差产生机理非常复杂,很难对其建立准确的模型;因此对MEMS陀螺仪进行误差建模与补偿是十分必要的。
发明内容
本发明的目的在于提供一种提高滤波效果提高陀螺仪精度的基于改进粒子滤波的MEMS陀螺随机误差补偿方法。
本发明的目的是这样实现的:
基于改进粒子滤波的MEMS陀螺随机误差补偿方法,包括:
(1)数据的采集、预处理与检验:
采集MEMS陀螺的输出数据,把安装在转台上惯导系统,采样频率设为100Hz,采样时间为20ms,进行陀螺输出数据采集,保存采集的数据,从采集的MEMS陀螺仪输出数据中选取前10000个数据,去除陀螺的确定性误差,得到包含噪声的陀螺随机漂移数据,对MEMS陀螺的输出数据进行预处理,得到平稳、正态、零均值的MEMS陀螺输出的时间序列;
(2)对处理后的MEMS陀螺的输出数据进行时间序列分析建模:
分析MEMS陀螺仪的输出数据,利用AIC准则确定模型为ARMA(2,1):
xk由yk一阶差分得到,最终得到的MEMS陀螺随机漂移的时间序列模型为:
为模型的输出,即估计的时间序列,ak为白噪声;
(3)利用小波网络算法优化粒子滤波,改善粒子滤波性能:
(3.1)初始化,k=0时,采样得到N个粒子
(3.2)对i=1,2,…,N计算重要性权值:
式中,表示k时刻第i个粒子的权值,
将得到的粒子权值用矩阵表示并按降序排列,;将权值矩阵W分为高权值矩阵 与低权值矩阵
(3.3)权值分裂:将高权值的粒子分裂成中阶权值的粒子并取代低权值的粒子;
(3.4)权值调整:利用小波网络调整低权值的粒子;
低权值的粒子作为小波网络的输入数据,粒子的状态值作为小波网络的初始权值,任意时刻的量测值为小波网络的期望输出;样本的学习函数为系统的量测方程,训练网络,小波网络输出的为调整后的粒子权值,将得到的新的粒子权值归一化处理;
(3.5)重采样:若则进行重采样。其中Neff为样本的有效抽样尺度,Nthreshold为提前设定的阈值;
(3.6)输出:
状态估计:
方差估计:
(3.7)判断是否结束,若是则退出,否则返回步骤(3.4);
(4)对时间序列模型进行粒子滤波,对MEMS陀螺仪误差补偿,提高精度:
基于所建模型,利用步骤(3)中改进的滤波算法对MEMS陀螺随机误差进行有效地补偿:
为随机时间序列;θj为自回归系数和滑动平均系数;ak为零均值的白噪声;
确定状态向量为观测向量为建立粒子滤波的系统和量测方程为:
X(k+1)=AX(k)+BW(k)
Z(k)=CX(k)+V(k)
对MEMS陀螺仪的随机漂移的时间序列模型进行滤波处理,粒子数N选为100。
本发明的有益效果:
MEMS陀螺仪的误差由确定性误差和随机漂移两部分组成。其中确定性误差可以通过实验的方法进行补偿;而随机漂移具有很大的随机性,且随时间的增加而变化,无线性规律可循,只能用统计模型来近似描述。在现有的陀螺仪随机漂移建模中,对于陀螺随机误差通常采用的是Allan方差方法和时间序列分析的建模方法。对于Allan方差方法是建立在统计学基础上的方法,仅适用于平稳的随机信号且对数据样本要求较多;而时间序列分析方法可以建立精确的MEMS陀螺仪随机漂移的模型,是提高陀螺仪精度的一种有效方法。虽然Kalman滤波可以获得最优的状态估计,但对于非线性、非高斯分布的状态模型,其滤波和预测精度很难保证。但是粒子滤波完全突破了Kalman滤波理论的框架,是一种基于蒙特卡罗思想的非线性、非高斯系统的滤波方,对系统的过程噪声和量测噪声没有任何限制。但是随着多次递归更新后粒子样本的多样性丧失,具有较高权值的粒子被多次统计计算,使得滤波效果受到影响。本发明采用一种改进的粒子滤波算法,将小波网络算法与粒子滤波典型采样算法结合,增大位于概率分布尾部的粒子权值,使较高权值的粒子分解为若干个较小权值的粒子,提高粒子样本的多样性,减小误差,提高滤波效果。将基于小波网络的粒子滤波算法应用到MEMS陀螺仪的误差补偿中,可以有效减少随机误差,提高MEMS陀螺仪精度。
附图说明
图1是基于改进粒子滤波的MEMS陀螺仪随机误差的滤波算法流程图;
图2是时间序列分析建模流程图;
图3是小波网络结构图;
图4是基于小波网络改进的粒子滤波算法的流程图;
图5是MEMS陀螺原始数据输出示意图;
图6是基于ARMA模型的粒子滤波图。
具体实施方式
下面结合附图对本发明做进一步描述。
针对MEMS陀螺仪精度较低、随机漂移较大,本发明提出了一种可以有效对MEMS陀螺仪随机误差补偿的方法。该方法的步骤如下:
步骤1:数据采集、预处理与检验。采集MEMS陀螺的输出信号,得到观测数据,并对数据进行预处理并检验,得到平稳、正态、零均值的时间序列。
步骤2:对处理后的MEMS陀螺的输出数据进行时间序列建模,确定模型结构及参数,并检验所建模型的适应性。
步骤3:利用小波网络算法优化粒子滤波,改善粒子滤波性能。
步骤4:对所建的时间序列模型进行粒子滤波,对MEMS陀螺仪误差补偿,提高精度。
本发明描述的方法是一种MEMS陀螺随机误差补偿方法,该发明采用了时间序列分析方法对MEMS陀螺输出数据建模,采用改进的粒子滤波方法对陀螺仪输出数据进行滤波处理,可以有效减小随机误差,明显降低陀螺仪的随机漂移,对误差进行有效的补偿,提高陀螺仪的精度。本发明设计方案如图1所示,步骤如下:
步骤1:数据采集、预处理与检验。
首先采集MEMS陀螺的输出数据。实际上MEMS陀螺输出是连续的,而时间序列分析法建模需要离散的数据;因此要以一定的采样频率采集数据。本发明中采样频率设为100Hz,采样时间为20ms,得到离散的数据序列yk。把惯导系统安装在转台上,接通电源预热15分钟,并对串口接收程序设置,利用所编写的导航系统界面进行陀螺输出数据采集。并将采集的数据保存在文件夹中。从采集的MEMS陀螺仪输出数据中选取前5000个数据,通过陀螺确定性误差补偿实验,去除陀螺的确定性误差,得到包含噪声的陀螺随机漂移数据。其次对MEMS陀螺的输出数据进行预处理。通过对MEMS陀螺输出数据分析,陀螺仪原始数据中的常值分量通过求均值来提取,存在的趋势项利用差分方程法去除。分析陀螺仪输出数据可知,作一阶差分即可去除趋势项。最后对数据进行检验。确定处理后的数据满足建模要求。
步骤2:对预处理后的MEMS陀螺输出数据进行时间序列分析建模。
时间序列分析法建模包括AR模型、MA模型和ARMA模型。通常可以把一个平稳的随机过程看成不相关的白噪声通过上述三种中的一种模型而产生的过程,ARMA(p,q)一般形式如下:
式中,xk为随机时间序列;p、q为自回归阶数和滑动平均阶数;θj(j=1,2,…,q)为自回归系数和滑动平均系数;ak为零均值的白噪声。
上式中若q=0,则变为p阶自回归模型,记为AR(p);若p=0,则变为q阶滑动平均模型,记为MA(q)。
分析MEMS陀螺仪的输出数据,可知自相关函数具有明显的拖尾现象,因此模型可以选为AR(p)模型或ARMA(p,q)模型。利用AIC准则确定模型为ARMA(2,1),形式如下:
xk由yk一阶差分得到,所以最终的MEMS陀螺随机漂移的时间序列模型为:
式中,为模型的输出,即估计的时间序列,ak为白噪声。
步骤3:利用小波网络算法优化粒子滤波算法。本发明提出的小波网络优化粒子滤波算法,调整粒子的权值,使其进入高权值区域,提高粒子样本的多样性,减小误差,避免粒子退化,可以有效提高滤波精度。具体实现过程如下所示:
第一步:初始化。k=0时,从重要性函数采样得到N个粒子重要性密度函数取转移先验:
式中,表示k时刻第i个粒子;zk为系统量测值;为k时刻的重要密度函数;为k时刻先验概率密度函数。
第二步:对i=1,2,…,N计算重要性权值。其重要性权值为:
式中,表示k时刻第i个粒子的权值。
粒子的权值矩阵按降序排列,以权值ωp为界,将权值矩阵W分为权值较高的矩阵 与权值较低的矩阵
第三步:权值分裂。根据上一步中分离出的高权值的粒子和低权值的粒子,将高权值的粒子分裂成较低权值的粒子并取代权值很低的粒子。
第四步:权值调整。在经过上一步中权值分裂步骤后,权值矩阵W变为其中,为按降序排列的粒子权值。取矩阵中粒子权权值较低的部分,利用小波网络调整权值较小的粒子。
小波网络以BP神经网络为拓扑结构,小波基函数作为隐含层激励函数,信号前向传播的同时误差反向传播的网络。小波网络的学习算法具有更稳定的性能,能收敛到全局最小点;且尺度因子可以根据误差进行修正,以使网络达到最好的性能。本发明中采用梯度修正法修正小波网络的权值,从而得到最佳的粒子权值。小波网络权值修正过程如下:
(1)计算网络预测误差e:
e=Zk-Yk,k=1,2…m (6)
式中,Zk为网络期望输出,Yk为网络预测输出。
(2)根据预测误差e修正小波网络权值
式中,是根据网络预测误差e计算得到的。
(8)
式中,η为学习率。
本发明中将权值较小的粒子作为小波网络的输入数据,粒子的状态值作为神经网络的初始权值,任意时刻的量测值为小波网络的期望输出,样本的学习函数为系统的量测方程。由图3知网络的隐含层输出公式为:
式中,h(j)为隐含层第j个节点输出值;为粒子的权值;wij为输入层和隐含层的连接权值;bj为小波基函数hj的平移因子;aj为小波基函数hj的伸缩因子;hj为小波基函数。
小波网络输出层计算公式为:
式中,wjk为隐含层到输出层权值;h(j)为第j个隐含层节点的输出;l为隐含层节点数;m为输出层节点数。
小波网络的学习步长设为0.05,小波基函数为Morlet小波函数训练小波网络,当训练误差满足要求时,小波网络输出为调整后的粒子权值。
将得到的新的粒子权值归一化处理:
式中,表示k时刻第i个粒子的权值;N表示粒子总数;为神经网络调整后的粒子权值。
第五步:重采样。若则进行重采样,即将原来的带权样本映射为等权样本其中Neff为样本的有效抽样尺度,Nthreshold为设定的阈值。
第六步:输出。
状态估计:
方差估计:
第七步:判断是否结束,若是则退出本算法,否则返回第四步。
步骤4:对所建的时间序列模型进行粒子滤波,对MEMS陀螺仪误差补偿。
基于所建模型,利用步骤三中改进的滤波算法来对MEMS陀螺随机误差进行有效地补偿。模型表达式为:
式中,为随机时间序列;θj为自回归系数和滑动平均系数;ak为零均值的白噪声。
由此可确定状态向量为观测向量为建立粒子滤波的系统和量测方程为:
X(k+1)=AX(k)+BW(k)
(10)
Z(k)=CX(k)+V(k)
(11)
其中, C=[1 0 0],W(k)=[ak ak-1]T。
根据上述状态方程和量测方程采用步骤4中改进的粒子滤波算法对MEMS陀螺仪的随机漂移的时间序列进行滤波处理。粒子数N选为100,进行仿真,如图6所示,可以看出步骤3中改进的粒子滤波算法对MEMS陀螺随机漂移起到了很好的抑制作用,有效地改善了陀螺仪的精度。
Claims (1)
1.基于改进粒子滤波的MEMS陀螺随机误差补偿方法,其特征在于:
(1)数据的采集、预处理与检验:
采集MEMS陀螺的输出数据,把安装在转台上惯导系统,采样频率设为100Hz,采样时间为20ms,进行陀螺输出数据采集,保存采集的数据,从采集的MEMS陀螺仪输出数据中选取前10000个数据,去除陀螺的确定性误差,得到包含噪声的陀螺随机漂移数据,对MEMS陀螺的输出数据进行预处理,得到平稳、正态、零均值的MEMS陀螺输出的时间序列;
(2)对处理后的MEMS陀螺的输出数据进行时间序列分析建模:
分析MEMS陀螺仪的输出数据,利用AIC准则确定模型为ARMA(2,1):
xk由yk一阶差分得到,最终得到的MEMS陀螺随机漂移的时间序列模型为:
为模型的输出,即估计的时间序列,ak为白噪声;
(3)利用小波网络算法优化粒子滤波,改善粒子滤波性能:
(3.1)初始化,k=0时,采样得到N个粒子
(3.2)对i=1,2,…,N计算重要性权值:
式中,表示k时刻第i个粒子的权值,
将得到的粒子权值用矩阵表示并按降序排列,;将权值矩阵W分为高权值矩阵 与低权值矩阵
(3.3)权值分裂:将高权值的粒子分裂成中阶权值的粒子并取代低权值的粒子;
(3.4)权值调整:利用小波网络调整低权值的粒子;
低权值的粒子作为小波网络的输入数据,粒子的状态值作为小波网络的初始权值,任意时刻的量测值为小波网络的期望输出;样本的学习函数为系统的量测方程,训练网络,小波网络输出的为调整后的粒子权值,将得到的新的粒子权值归一化处理;
(3.5)重采样:若则进行重采样。其中Neff为样本的有效抽样尺度,Nthreshold为提前设定的阈值;
(3.6)输出:
状态估计:
方差估计:
(3.7)判断是否结束,若是则退出,否则返回步骤(3.4);
(4)对时间序列模型进行粒子滤波,对MEMS陀螺仪误差补偿,提高精度:
基于所建模型,利用步骤(3)中改进的滤波算法对MEMS陀螺随机误差进行有效地补偿:
为随机时间序列;θj为自回归系数和滑动平均系数;ak为零均值的白噪声;
确定状态向量为观测向量为建立粒子滤波的系统和量测方程为:
X(k+1)=AX(k)+BW(k)
Z(k)=CX(k)+V(k)
对MEMS陀螺仪的随机漂移的时间序列模型进行滤波处理,粒子数N选为100。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410298926.3A CN104048676B (zh) | 2014-06-26 | 2014-06-26 | 基于改进粒子滤波的mems陀螺随机误差补偿方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410298926.3A CN104048676B (zh) | 2014-06-26 | 2014-06-26 | 基于改进粒子滤波的mems陀螺随机误差补偿方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104048676A true CN104048676A (zh) | 2014-09-17 |
CN104048676B CN104048676B (zh) | 2017-02-15 |
Family
ID=51501834
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410298926.3A Expired - Fee Related CN104048676B (zh) | 2014-06-26 | 2014-06-26 | 基于改进粒子滤波的mems陀螺随机误差补偿方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104048676B (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104251712A (zh) * | 2014-10-09 | 2014-12-31 | 哈尔滨工程大学 | 基于小波多尺度分析的mems陀螺随机误差补偿方法 |
CN104820757A (zh) * | 2015-05-18 | 2015-08-05 | 中国电子科技集团公司第二十研究所 | 基于物理模型的mems陀螺温漂特性神经网络建模方法 |
CN107330149A (zh) * | 2017-05-27 | 2017-11-07 | 哈尔滨工业大学 | 基于arma和bpnn组合模型的mimu陀螺随机漂移预测方法 |
CN108444471A (zh) * | 2018-05-03 | 2018-08-24 | 天津大学 | 一种基于粒子滤波与小波变换的加速度计信号去噪方法 |
CN109379762A (zh) * | 2018-09-10 | 2019-02-22 | 中国联合网络通信集团有限公司 | 一种随机信号数据统计方法和系统 |
CN109540808A (zh) * | 2018-11-02 | 2019-03-29 | 湖南文理学院 | 一种变压器检测系统和故障诊断方法 |
CN112212862A (zh) * | 2020-09-24 | 2021-01-12 | 天津理工大学 | 一种改进粒子滤波的gps/ins组合导航方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101813493A (zh) * | 2010-04-09 | 2010-08-25 | 南京航空航天大学 | 一种基于粒子滤波的惯性导航系统初始对准方法 |
US20110167893A1 (en) * | 2008-09-17 | 2011-07-14 | Bae Systems Information And Electronic Systems Integration Inc. | Method and apparatus for in-flight calibration of gyroscope using magnetometer reference |
WO2012159939A1 (fr) * | 2011-05-20 | 2012-11-29 | Sagem Defense Securite | Procede de calibration d'un ensemble inertiel comportant une phase dynamique entre deux phases statiques |
CN103217176A (zh) * | 2013-05-06 | 2013-07-24 | 南京智真电子科技有限公司 | 一种消除mems陀螺零漂的方法 |
CN103674062A (zh) * | 2013-12-05 | 2014-03-26 | 广东电网公司电力科学研究院 | 基于Allan方差与ARMA模型分析提高陀螺仪测量精度的方法 |
-
2014
- 2014-06-26 CN CN201410298926.3A patent/CN104048676B/zh not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20110167893A1 (en) * | 2008-09-17 | 2011-07-14 | Bae Systems Information And Electronic Systems Integration Inc. | Method and apparatus for in-flight calibration of gyroscope using magnetometer reference |
CN101813493A (zh) * | 2010-04-09 | 2010-08-25 | 南京航空航天大学 | 一种基于粒子滤波的惯性导航系统初始对准方法 |
WO2012159939A1 (fr) * | 2011-05-20 | 2012-11-29 | Sagem Defense Securite | Procede de calibration d'un ensemble inertiel comportant une phase dynamique entre deux phases statiques |
CN103217176A (zh) * | 2013-05-06 | 2013-07-24 | 南京智真电子科技有限公司 | 一种消除mems陀螺零漂的方法 |
CN103674062A (zh) * | 2013-12-05 | 2014-03-26 | 广东电网公司电力科学研究院 | 基于Allan方差与ARMA模型分析提高陀螺仪测量精度的方法 |
Non-Patent Citations (3)
Title |
---|
刘宇等: "基于提升小波的MEMS陀螺去噪方法研究", 《压电与声光》 * |
李杰等: "基于时间序列分析的Kalman滤波方法在MEMS陀螺仪随机漂移误差补偿中的应用研究", 《传感技术学报》 * |
钱华明等: "基于Kalman滤波的MEMS陀螺仪滤波算法", 《哈尔滨工业大学学报》 * |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104251712A (zh) * | 2014-10-09 | 2014-12-31 | 哈尔滨工程大学 | 基于小波多尺度分析的mems陀螺随机误差补偿方法 |
CN104251712B (zh) * | 2014-10-09 | 2017-10-31 | 哈尔滨工程大学 | 基于小波多尺度分析的mems陀螺随机误差补偿方法 |
CN104820757A (zh) * | 2015-05-18 | 2015-08-05 | 中国电子科技集团公司第二十研究所 | 基于物理模型的mems陀螺温漂特性神经网络建模方法 |
CN104820757B (zh) * | 2015-05-18 | 2018-02-06 | 中国电子科技集团公司第二十研究所 | 基于物理模型的mems陀螺温漂特性神经网络建模方法 |
CN107330149A (zh) * | 2017-05-27 | 2017-11-07 | 哈尔滨工业大学 | 基于arma和bpnn组合模型的mimu陀螺随机漂移预测方法 |
CN107330149B (zh) * | 2017-05-27 | 2022-10-28 | 哈尔滨工业大学 | 基于arma和bpnn组合模型的mimu陀螺随机漂移预测方法 |
CN108444471A (zh) * | 2018-05-03 | 2018-08-24 | 天津大学 | 一种基于粒子滤波与小波变换的加速度计信号去噪方法 |
CN109379762A (zh) * | 2018-09-10 | 2019-02-22 | 中国联合网络通信集团有限公司 | 一种随机信号数据统计方法和系统 |
CN109379762B (zh) * | 2018-09-10 | 2021-09-14 | 中国联合网络通信集团有限公司 | 一种随机信号数据统计方法和系统 |
CN109540808A (zh) * | 2018-11-02 | 2019-03-29 | 湖南文理学院 | 一种变压器检测系统和故障诊断方法 |
CN112212862A (zh) * | 2020-09-24 | 2021-01-12 | 天津理工大学 | 一种改进粒子滤波的gps/ins组合导航方法 |
Also Published As
Publication number | Publication date |
---|---|
CN104048676B (zh) | 2017-02-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104048676A (zh) | 基于改进粒子滤波的mems陀螺随机误差补偿方法 | |
CN105843073B (zh) | 一种基于气动力不确定降阶的机翼结构气动弹性稳定性分析方法 | |
Wang et al. | Central difference particle filter applied to transfer alignment for SINS on missiles | |
CN107084714B (zh) | 一种基于RoboCup3D的多机器人协作目标定位方法 | |
CN106683122A (zh) | 一种基于高斯混合模型和变分贝叶斯的粒子滤波方法 | |
CN103900610A (zh) | 基于灰色小波神经网络的mems陀螺随机误差预测方法 | |
CN103217175A (zh) | 一种自适应容积卡尔曼滤波方法 | |
CN110728357B (zh) | 一种基于循环神经网络的imu数据去噪方法 | |
Slivinski et al. | A hybrid particle–ensemble Kalman filter for Lagrangian data assimilation | |
CN104112079A (zh) | 一种模糊自适应变分贝叶斯无迹卡尔曼滤波方法 | |
CN103973263B (zh) | 一种逼近滤波方法 | |
CN104506162A (zh) | 基于ls-svr建模的高阶粒子滤波器的故障预示方法 | |
CN103604430A (zh) | 一种基于边缘化ckf重力辅助导航的方法 | |
CN104251712A (zh) | 基于小波多尺度分析的mems陀螺随机误差补偿方法 | |
Xu et al. | A novel robust filter for outliers and time-varying delay on an SINS/USBL integrated navigation model | |
CN115270239A (zh) | 基于动力特性和智能算法响应面法的桥梁可靠性预测方法 | |
CN103776449A (zh) | 一种提高鲁棒性的动基座初始对准方法 | |
CN105678015A (zh) | 一种高超声速三维机翼的非概率可靠性气动结构耦合优化设计方法 | |
CN103793614B (zh) | 一种突变滤波方法 | |
CN109471192A (zh) | 一种全自动重力测试仪高精度动态数据处理方法 | |
CN110186483B (zh) | 提高惯性制导航天器落点精度的方法 | |
CN106646452A (zh) | 一种基于摄动多高斯拟合的空间目标跟踪方法 | |
CN104320108A (zh) | 基于ahcif的集中式测量值加权融合方法 | |
Rezaie et al. | Shrinked (1− α) ensemble Kalman filter and α Gaussian mixture filter | |
CN104268597A (zh) | 基于ahcif的集中式测量值扩维融合方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170215 |