CN112269192B - 一种快速自适应的动态北斗监测实时解算去噪方法 - Google Patents
一种快速自适应的动态北斗监测实时解算去噪方法 Download PDFInfo
- Publication number
- CN112269192B CN112269192B CN202011140262.XA CN202011140262A CN112269192B CN 112269192 B CN112269192 B CN 112269192B CN 202011140262 A CN202011140262 A CN 202011140262A CN 112269192 B CN112269192 B CN 112269192B
- Authority
- CN
- China
- Prior art keywords
- moment
- estimated value
- system state
- value
- state
- 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
- 238000000034 method Methods 0.000 title claims abstract description 42
- 238000012544 monitoring process Methods 0.000 title claims abstract description 20
- 239000011159 matrix material Substances 0.000 claims abstract description 53
- 238000001914 filtration Methods 0.000 claims abstract description 39
- 238000005259 measurement Methods 0.000 claims abstract description 14
- 230000008569 process Effects 0.000 claims abstract description 12
- 238000006073 displacement reaction Methods 0.000 claims description 20
- 238000006243 chemical reaction Methods 0.000 claims description 9
- 238000012545 processing Methods 0.000 claims description 8
- 230000009466 transformation Effects 0.000 claims description 4
- 230000007704 transition Effects 0.000 claims description 4
- 238000013459 approach Methods 0.000 abstract description 5
- 230000003044 adaptive effect Effects 0.000 description 7
- 238000004364 calculation method Methods 0.000 description 6
- 230000000694 effects Effects 0.000 description 6
- 238000010586 diagram Methods 0.000 description 5
- 238000010276 construction Methods 0.000 description 3
- 238000001514 detection method Methods 0.000 description 3
- 238000011161 development Methods 0.000 description 3
- 238000012795 verification Methods 0.000 description 3
- 239000000654 additive Substances 0.000 description 2
- 230000000996 additive effect Effects 0.000 description 2
- 238000012937 correction Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000001228 spectrum Methods 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- 230000006978 adaptation Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- JJWKPURADFRFRB-UHFFFAOYSA-N carbonyl sulfide Chemical compound O=C=S JJWKPURADFRFRB-UHFFFAOYSA-N 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000005530 etching Methods 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 238000005065 mining Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/13—Receivers
- G01S19/14—Receivers specially adapted for specific applications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/13—Receivers
- G01S19/21—Interference related issues ; Issues related to cross-correlation, spoofing or other methods of denial of service
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/13—Receivers
- G01S19/35—Constructional details or hardware or software details of the signal processing chain
- G01S19/37—Hardware or software details of the signal processing chain
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Signal Processing (AREA)
- Testing Or Calibration Of Command Recording Devices (AREA)
Abstract
本发明提供一种快速自适应的动态北斗监测实时解算去噪方法,包括:初始条件;令k=1;得到k时刻的***状态中间估计值得到k‑1时刻***状态噪声矩阵Qk‑1;得到k时刻的误差协方差矩阵中间估计值得到测量过程噪声R;得到k时刻的增益系数Kk;得到k时刻误差协方差矩阵最终估计值Pk;得到k时刻的***状态最终估计值xk;令k=k+1,返回前面步骤,直到对所有时刻的测量值完成滤波操作。优点:通过对原始信号进行滤波,使滤波后的值逼近理论值,进行使滤波后得到的测量值能够更为真实的反映边坡状态,从而实现对边坡的准确监测。
Description
技术领域
本发明属于数据去噪技术领域,具体涉及一种快速自适应的动态北斗监测实时解算去噪方法。
背景技术
伴随着我国国民经济的进一步发展和全国基础性工程建设战略的持续推进,我国加强了桥梁道路等交通建设,全国地区的高等级公路和铁路都得到极其迅猛的发展。但是,西部山区地形,滑坡、崩塌等自然灾害频发,具有分布广泛、发生率高以及快速成灾等特征。据不完全统计,在近10年左右的时间里,平均每年有重大崩滑发生。所以,边坡稳定性的研究已成为工程地质学研究的重要任务之一。
边坡的安全性已成为交通建设的一个重要属性,所以有必要对边坡进行有效的实时监测,确保桥梁等设施正常运行。边坡主要存在变形导致位移甚至滑坡的现象,为了避免灾害,如果发现边坡位移速度急剧增大,应立即撤出人员、采掘和运输设备,以确保人员和财产的安全。通过监测变形的规律,可以进一步的研究边坡变形的原因,预报灾害,为以后的边坡设计积累经验,例如,通过边坡稳定性分析的基础资料,预报滑坡。还可以依据观测资料了解和掌握滑坡的形态、规模和发展趋势,以便采取相应的处理措施。
常见的边坡检测有基于北斗GNSS静态测量方法,虽然静态测量时精度高,但是需要较长的时间,不能及时、实时反映边坡变化的情况。目前,采用位移传感器的方法可获取边坡位移数据,进而判断边坡健康情况。但是,此种方式主要存在以下问题:
由于环境中的白噪声、热噪声、宇宙噪声和一些加性噪声的影响,使得位移传感器数据产生抖动,对边坡检测造成极大的影响,降低边坡检测的有效性。
发明内容
针对现有技术存在的缺陷,本发明提供一种快速自适应的动态北斗监测实时解算去噪方法,可有效解决上述问题。
本发明采用的技术方案如下:
本发明提供一种快速自适应的动态北斗监测实时解算去噪方法,包括以下步骤:
步骤1,初始条件:
设置初始时刻***状态最终估计值x0已知;初始时刻误差协方差矩阵最终估计值P0已知,并且,P0=0;
步骤2,令k=1;
步骤3,建立公式(1)所示的位移传感器***的第一状态转换方程:
其中:
为k时刻的***状态中间估计值;
A为常数,由于位移传感器***为一维信号处理的***,因此,取值为1;
xk-1为k-1时刻***状态最终估计值;
获得k-1时刻***状态最终估计值xk-1;根据位移传感器***的第一状态转换方程,计算得到k时刻的***状态中间估计值
步骤4,根据得到的k时刻的***状态中间估计值以及k-1时刻***状态最终估计值xk-1,根据下式计算得到k-1时刻***状态噪声矩阵Qk-1:
步骤5,得到k-1时刻***状态噪声矩阵Qk-1后,再根据k-1时刻***状态噪声矩阵Qk-1以及k-1时刻误差协方差矩阵最终估计值Pk-1,根据下面的位移传感器***的第二状态方程,计算得到k时刻的误差协方差矩阵中间估计值
步骤6,根据下式,计算得到测量过程噪声R:
步骤7,将步骤5得到的k时刻的误差协方差矩阵中间估计值和步骤6得到的测量过程噪声R代入以下的增益系数表达式,得到k时刻的增益系数Kk:
其中:
H为测量值与测量值对应的***状态值的变换矩阵,由于测量的参数与***状态值参数相同,没有转换关系,因此,取其为1;
步骤8,将步骤7得到的k时刻的增益系数Kk、步骤5得到的k时刻的误差协方差矩阵中间估计值代入以下误差协方差矩阵更新方程,得到k时刻误差协方差矩阵最终估计值Pk:
步骤9,将k时刻的测量值yk、步骤3得到的k时刻的***状态中间估计值步骤7得到的k时刻的增益系数Kk代入以下的观测值状态更新方程,得到k时刻的***状态最终估计值xk:
步骤10,令k=k+1,返回步骤3,直到对所有时刻的测量值完成滤波操作。
本发明提供的一种快速自适应的动态北斗监测实时解算去噪方法具有以下优点:本发明提供一种快速自适应的动态北斗监测实时解算去噪方法,是一种快速自适应卡尔曼滤波方法,通过计算相邻两个时刻的差分值,估计此时间段内的信号噪声协方差矩阵,从而不断更新,通过对原始信号进行滤波,使滤波后的值逼近理论值,进行使滤波后得到的测量值能够更为真实的反映边坡状态,从而实现对边坡的准确监测。
附图说明
图1为本发明提供的一种快速自适应的动态北斗监测实时解算去噪方法的流程示意图;
图2为采用传统卡尔曼滤波对测量值的滤波结果图;
图3为采用传统卡尔曼滤波对状态误差Q的估计结果图;
图4为采用本发明快速自适应实时解算去噪方法进行滤波后的滤波结果图;
图5为采用本发明快速自适应实时解算去噪方法进行状态误差Q的估计结果图。
具体实施方式
为了使本发明所解决的技术问题、技术方案及有益效果更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅用以解释本发明,并不用于限定本发明。
数字信号的滤波是通过数字滤波器来实现的。数字滤波器是一种用于过滤时间离散信号的数字***,是通过对抽样数据进行数学处理来达到频域滤波的目的。最开始通过对信号进行求加权平均来实现对加性噪声的去除,后来通过设计模拟滤波器实现对其他噪声的去除,随着数字信号的发展,逐渐采用FFT变换进行去噪处理,即:首先经变换得到频谱,通过分析频谱进而知道有用信号与噪声所在的带宽,然后设计对应的数字滤波器实现对信号的去噪处理。虽然低通高通等滤波器能实现去噪,但是对于随机噪声的处理效果不佳,尤其是对传感器的数据信号处理达不到有效的应用。因此需要提出一种对随机噪声具有优良性能的滤波方法
卡尔曼滤波器在随机序列估计、空间技术、目标跟踪领域具有广泛的应用,1960年卡尔曼还提出能控性的概念。后人称之为卡尔曼滤波,其实质是由量测值重构***的状态向量。以“预测—实测—修正”的顺序递推,卡尔曼滤波是以最小均方误差为估计的最佳准则,来寻求一套递推估计的算法,其基本思想是:采用信号与噪声的状态空间模型,利用前一时刻的估计值和现时刻的观测值来更新对状态变量的估计,求出现时刻的估计值。它适合于实时处理和计算机运算,卡尔曼滤波的实质是由量测值重构***的状态向量。它以“预测—实测—修正”的顺序递推,根据***的量测值来消除随机干扰,再现***的状态。
但在实际应用中,随着量测值数目的增加,由于估计误差的均值和估计误差协方差可能越来越大,使滤波逐渐失去准确估计的作用,这种现象称为卡尔曼滤波发散。引起滤波器误差发散的主要原因有两点:
(1)描述***动力学特性的数学模型和噪声估计模型不准确,不能直接真实地反映物理过程,使得模型与获得的量测值不匹配而导致滤波发散。这种由于模型建立过于粗糙或失真所引起的发散称为滤波发散。
(2)由于卡尔曼滤波是递推过程,随着滤波步数的增加,舍入误差将逐渐积累。如果计算机字长不够长,这种积累误差很有可能使估计误差方差阵失去非负定性甚至失去对称性,使滤波增益矩阵逐渐失去合适的加权作用而导致发散。这种由于计算舍入误差所引起的发散称为计算发散。
针对上述卡尔曼滤波发散的原因,本发明提供一种快速自适应的动态北斗监测实时解算去噪方法,是一种快速自适应卡尔曼滤波方法,通过计算相邻两个时刻的差分值,估计此时间段内的信号噪声协方差矩阵,从而不断更新,通过对原始信号进行滤波,使滤波后的值逼近理论值,进行使滤波后得到的测量值能够更为真实的反映边坡状态,从而实现对边坡的准确监测。
参考图1,本发明提供一种快速自适应的动态北斗监测实时解算去噪方法,包括以下步骤:
步骤1,初始条件:
设置初始时刻***状态最终估计值x0已知;初始时刻误差协方差矩阵最终估计值P0已知,并且,P0=0;
步骤2,令k=1;
步骤3,建立公式(1)所示的位移传感器***的第一状态转换方程:
其中:
为k时刻的***状态中间估计值;
A为常数,由于位移传感器***为一维信号处理的***,因此,取值为1;
xk-1为k-1时刻***状态最终估计值;
获得k-1时刻***状态最终估计值xk-1;根据位移传感器***的第一状态转换方程,计算得到k时刻的***状态中间估计值
步骤4,根据得到的k时刻的***状态中间估计值以及k-1时刻***状态最终估计值xk-1,根据下式计算得到k-1时刻***状态噪声矩阵Qk-1:
步骤5,得到k-1时刻***状态噪声矩阵Qk-1后,再根据k-1时刻***状态噪声矩阵Qk-1以及k-1时刻误差协方差矩阵最终估计值Pk-1,根据下面的位移传感器***的第二状态方程,计算得到k时刻的误差协方差矩阵中间估计值
步骤6,根据下式,计算得到测量过程噪声R:
步骤7,将步骤5得到的k时刻的误差协方差矩阵中间估计值和步骤6得到的测量过程噪声R代入以下的增益系数表达式,得到k时刻的增益系数Kk:
其中:
H为测量值与测量值对应的***状态值的变换矩阵,由于测量的参数与***状态值参数相同,没有转换关系,因此,取其为1;
步骤8,将步骤7得到的k时刻的增益系数Kk、步骤5得到的k时刻的误差协方差矩阵中间估计值代入以下误差协方差矩阵更新方程,得到k时刻误差协方差矩阵最终估计值Pk:
步骤9,将k时刻的测量值yk、步骤3得到的k时刻的***状态中间估计值步骤7得到的k时刻的增益系数Kk代入以下的观测值状态更新方程,得到k时刻的***状态最终估计值xk:
步骤10,令k=k+1,返回步骤3,直到对所有时刻的测量值完成滤波操作。
因此,对于k时刻的测量值yk,通过以上步骤3到步骤9的去噪滤波方法,得到滤波后的值,即:k时刻的***状态最终估计值xk,从而实现对k时刻的测量值yk的滤波。
在本发明提供的快速自适应的动态北斗监测实时解算去噪方法中,k时刻的***状态最终估计值xk和k时刻的***状态中间估计值均为估计值,只是估计精确程度不同,k时刻的***状态中间估计值/>是直接用第一状态转换方程算出来的估计值,而k时刻的***状态最终估计值xk是对测量值进行修正后的估计值,k时刻的***状态最终估计值xk更精确,作为下一个循环中***状态的初始值。
同样的,k时刻的误差协方差矩阵中间估计值和k时刻误差协方差矩阵最终估计值Pk均为估计值,只是估计精确程度不同,k时刻的误差协方差矩阵中间估计值/>是直接用第二状态转换方程算出来的估计值,而k时刻误差协方差矩阵最终估计值Pk是对k时刻的误差协方差矩阵中间估计值/>进行修正后的估计值,k时刻误差协方差矩阵最终估计值Pk更精确,作为下一个循环中***状态的初始值。
本发明中,初始条件,k=1,初始时刻***的精确状态X0已知,可以看作服从协方差为0的正态分布,X0~N(x0,P0),其中,x0为初始时刻***状态最终估计值,已知;P0为初始时刻误差协方差矩阵最终估计值,P0=0;
建立公式(1)所示的位移传感器***的第一状态转换方程时,设Xk-1为k-1时刻***的精确状态,服从于正态分布Xk-1~N(xk-1,Pk-1),以k-1时刻为基础,仅由状态转换方程进行计算,得到的***k时刻的精确状态Xk服从的正态分布为
用Xk表示***k时刻的精确状态,Xk不可直接求得或测量得到,其服从于一个正态分布,整个计算过程就是在求Xk这个分布的估计。xk、和Pk都是这个正态分布的系数的估计值。其中,/>是正态分布的均值,/>是协方差矩阵的估计值,xk是一个较为精确的正态分布的均值,可以将这个xk近似地认为是我们要求的精确值,并用其作为下一步计算的基出。
与传统卡尔曼不同的是,本发明在步骤4中,通过计算前后时刻的***状态估计值的差值的平方作为***状态噪声矩阵,在步骤6中对测量过程噪声R进行估计计算,并以测量过程噪声R为基础,得到了较准确的增益系数Kk,最后再通过增益系数Kk修正误差协方差矩阵中间估计值,通过增益系数Kk修正***状态中间估计值,进而得到较准确的***状态最终估计值和误差协方差矩阵最终估计值。采用本发明提供的快速自适应的动态北斗监测实时解算去噪方法,实现对原始信号进行滤波,使滤波后的值逼近理论值,进行使滤波后得到的测量值能够更为真实的反映边坡状态,从而实现对边坡的准确监测。
下面列举一个具体验证例:
本验证例分析了某桥梁边坡位移传感器的实时数据,该位移传感器每5秒上传一次测量数据,本验证例中位移传感器横坐标X的真实值为-0.45,将采用传统卡尔曼滤波的结果与本发明提出的快速自适应卡尔曼滤波的结果进行对比分析,如图2所示,为采用传统卡尔曼滤波对测量值的滤波结果图;如图3所示,为采用传统卡尔曼滤波对状态误差Q的估计结果图。
由图2可以看出,传统卡尔曼滤波的效果不佳,主要是由于对实时协方差矩阵估计不准确导致的,从图中看出精确值在-0.45附近,而滤波后的结果却大幅振荡。
由图3可以看出,在第二次迭代以后状态误差Q就趋于一个定值,但是实际上状态误差随着时间是变化的,因此导致最后的估计值达不到理想效果。
如图4所示,为采用本发明快速自适应实时解算去噪方法进行滤波后的滤波结果图;如图5所示,为采用本发明快速自适应实时解算去噪方法进行状态误差Q的估计结果图。
由图4可以看出,本发明的快速自适应卡尔曼滤波能精确地估计每一个***状态,即使在外界误差变化没有规律,即误差方差可变的情况下也能具有良好的估计性能,因此较传统卡尔曼相比,本发明滤波性能有较大提升。
由图5可以看出,在第二次迭代以后状态误差Q随时间变化而变化,与实际情况相符,因此最后的估计值达到理想数据滤波效果。
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视本发明的保护范围。
Claims (1)
1.一种快速自适应的动态北斗监测实时解算去噪方法,其特征在于,包括以下步骤:
步骤1,初始条件:
设置初始时刻***状态最终估计值x0已知;初始时刻误差协方差矩阵最终估计值P0已知,并且,P0=0;
步骤2,令k=1;
步骤3,建立公式(1)所示的位移传感器***的第一状态转换方程:
其中:
为k时刻的***状态中间估计值;
A为常数,由于位移传感器***为一维信号处理的***,因此,取值为1;
xk-1为k-1时刻***状态最终估计值;
获得k-1时刻***状态最终估计值xk-1;根据位移传感器***的第一状态转换方程,计算得到k时刻的***状态中间估计值
步骤4,根据得到的k时刻的***状态中间估计值以及k-1时刻***状态最终估计值xk-1,根据下式计算得到k-1时刻***状态噪声矩阵Qk-1:
步骤5,得到k-1时刻***状态噪声矩阵Qk-1后,再根据k-1时刻***状态噪声矩阵Qk-1以及k-1时刻误差协方差矩阵最终估计值Pk-1,根据下面的位移传感器***的第二状态方程,计算得到k时刻的误差协方差矩阵中间估计值
步骤6,根据下式,计算得到测量过程噪声R:
步骤7,将步骤5得到的k时刻的误差协方差矩阵中间估计值和步骤6得到的测量过程噪声R代入以下的增益系数表达式,得到k时刻的增益系数Kk:
其中:
H为测量值与测量值对应的***状态值的变换矩阵,由于测量的参数与***状态值参数相同,没有转换关系,因此,取其为1;
步骤8,将步骤7得到的k时刻的增益系数Kk、步骤5得到的k时刻的误差协方差矩阵中间估计值代入以下误差协方差矩阵更新方程,得到k时刻误差协方差矩阵最终估计值Pk:
步骤9,将k时刻的测量值yk、步骤3得到的k时刻的***状态中间估计值步骤7得到的k时刻的增益系数Kk代入以下的观测值状态更新方程,得到k时刻的***状态最终估计值xk:
步骤10,令k=k+1,返回步骤3,直到对所有时刻的测量值完成滤波操作。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011140262.XA CN112269192B (zh) | 2020-10-22 | 2020-10-22 | 一种快速自适应的动态北斗监测实时解算去噪方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011140262.XA CN112269192B (zh) | 2020-10-22 | 2020-10-22 | 一种快速自适应的动态北斗监测实时解算去噪方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112269192A CN112269192A (zh) | 2021-01-26 |
CN112269192B true CN112269192B (zh) | 2024-02-02 |
Family
ID=74343037
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011140262.XA Active CN112269192B (zh) | 2020-10-22 | 2020-10-22 | 一种快速自适应的动态北斗监测实时解算去噪方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112269192B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113654554B (zh) * | 2021-08-13 | 2024-06-14 | 北京维思韦尔航空电子技术有限公司 | 一种快速自适应的动态惯性导航实时解算去噪方法 |
CN114046785A (zh) * | 2021-11-10 | 2022-02-15 | 广东微电科技有限公司 | 磁探测信号线性噪声滤波方法及***、计算机可读存储介质、磁导航传感器、agv小车 |
CN114912551B (zh) * | 2022-07-18 | 2023-04-07 | 中国铁路设计集团有限公司 | 面向桥梁变形监测的gnss和加速度计实时融合方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109902568A (zh) * | 2019-01-22 | 2019-06-18 | 长江大学 | 一种自适应卡尔曼滤波算法 |
CN110850450A (zh) * | 2019-12-03 | 2020-02-28 | 航天恒星科技有限公司 | 一种卫星钟差参数的自适应估计方法 |
CN111623703A (zh) * | 2020-07-28 | 2020-09-04 | 湖南联智科技股份有限公司 | 一种基于新型卡尔曼滤波的北斗变形监测实时处理方法 |
-
2020
- 2020-10-22 CN CN202011140262.XA patent/CN112269192B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109902568A (zh) * | 2019-01-22 | 2019-06-18 | 长江大学 | 一种自适应卡尔曼滤波算法 |
CN110850450A (zh) * | 2019-12-03 | 2020-02-28 | 航天恒星科技有限公司 | 一种卫星钟差参数的自适应估计方法 |
CN111623703A (zh) * | 2020-07-28 | 2020-09-04 | 湖南联智科技股份有限公司 | 一种基于新型卡尔曼滤波的北斗变形监测实时处理方法 |
Non-Patent Citations (1)
Title |
---|
一种基于卡尔曼滤波处理的北斗卫星无源组合导航自适应野值剔除方法;王光鼎;张升康;杨汝良;;电子与信息学报(08);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN112269192A (zh) | 2021-01-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112269192B (zh) | 一种快速自适应的动态北斗监测实时解算去噪方法 | |
WO2020191980A1 (zh) | 一种无线传感网络数据漂移盲校准方法 | |
CN104459696A (zh) | 一种基于平地相位的sar干涉基线精确估计方法 | |
US7379844B2 (en) | Particle sampling method and sensor fusion and filtering method | |
CN110231636B (zh) | Gps与bds双模卫星导航***的自适应无迹卡尔曼滤波方法 | |
CN112487604B (zh) | 海洋重力仪输出数据长时间非线性漂移补偿方法 | |
CN113589286B (zh) | 基于D-LinkNet的无迹卡尔曼滤波相位解缠方法 | |
CN113820003B (zh) | 一种适用于桥梁振动监测的加速度实时重构动态位移方法 | |
WO2024016369A1 (zh) | 面向桥梁变形监测的gnss和加速度计实时融合算法 | |
CN113552565B (zh) | 针对sar数据高噪声及大梯度变化区域的相位解缠方法 | |
CN115510578B (zh) | 基于InSAR近实时监测的滑坡失稳时间概率预测方法及产品 | |
CN108120452B (zh) | Mems陀螺仪动态数据的滤波方法 | |
CN110553631B (zh) | 一种关于水位流量关系的水位测量系列误差分析方法 | |
CN110501686A (zh) | 基于一种新型自适应高阶无迹卡尔曼滤波的状态估计方法 | |
CN113281749A (zh) | 一种顾及同质性的时序InSAR高相干点选取方法 | |
CN117387884A (zh) | 一种基于多传感器数据融合的桥梁挠度测量方法 | |
CN109471192B (zh) | 一种全自动重力测试仪高精度动态数据处理方法 | |
CN111427096A (zh) | 全张量重力梯度仪数据质量评价和滤波处理方法 | |
CN109557594B (zh) | 基于重力异常时变的重力基准图时变修正方法及*** | |
KR101693310B1 (ko) | 평면파를 이용한 이방성 매질에서의 전파형 역산방법 | |
CN114236480A (zh) | 一种机载平台传感器***误差配准算法 | |
CN105678716A (zh) | 一种地基sar大气干扰相位校正方法及装置 | |
CN112632454A (zh) | 一种基于自适应卡尔曼滤波算法的mems陀螺滤波方法 | |
CN109581494B (zh) | 叠前偏移方法及装置 | |
CN113654554B (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 |