一种零相位升沉测量方法及系统技术方案

本发明专利技术提供了一种零相位升沉测量方法及系统，包括设置延迟时间和边界条件确定升沉计算区间；采集升沉计算区间内惯性导航系统测量的等效垂向加速度，对其进行滤波处理得到垂向加速度正向滤波输出数列；垂向加速度正向滤波输出数列进行倒序后进行高通滤波，得到垂向加速度逆向滤波输出数列；垂向加速度逆向滤波输出数列进行倒序后得到按时间正向排列的数列即为等效垂向加速度的滤波结果；对滤波结果进行二次积分，得到升沉计算区间内对应各点的升沉测量结果

【技术实现步骤摘要】
一种零相位升沉测量方法及系统


[0001]本专利技术属于惯性测量算法
，具体涉及一种零相位升沉测量方法及系统
。

技术介绍

[0002]水面舰船
、
水下航行器等水中运载体在航行过程中，受到波浪
、
风力因素的影响，会出现周期性的往复运动，称为升沉运动
。
舰船升沉准确测量在舰船减摇
、
舰载机起降
、
水下精确测量等应用领域均有重要作用
。
但由于水面对电磁波的反射作用，以及水下地形的复杂性等因素，通过
GNSS、
深度计等外界参考信息难以对升沉运动进行准确的测量
。
惯性导航系统在测量舰船的升沉运动上具有独到的优势，虽然惯导系统存在周期性振荡误差，但与升沉周期相比属于极低频分量，因此可以通过高通滤波器对惯导系统输出的位置信息进行滤波，但任何实时滤波算法都会产生相位差，导致滤波结果存在时间延迟，通过设计低延迟的高通滤波器，能够将特定周期内的信号测量延迟控制在
0.1s
以内，测量误差能够控制在厘米级，可以满足一般舰船控制系统精度要求
。
对于某些应用场景，如水下测量
、
建图等领域，精度需要达到毫米级，实时测量方法难以满足要求，但这些领域对于实时性要求不高，因此可以通过后处理的方式，实现零相位延迟的高精度升沉测量
。

技术实现思路

[0003]本专利技术的目的在于提供一种零相位升沉测量方法及系统，采用基于时间延迟的方法，构建零相位滤波器，保证升沉测量结果的准确性，并针对区间延迟计算边界不重合的问题，设计了边界条件计算方法
。
[0004]为了实现上述目的，本专利技术提供了一种零相位升沉测量方法，包括如下步骤
[0005]依据使用要求设置计算窗口时间宽度，将数据序列按计算窗口时间宽度分段；
[0006]确定延迟时间
、
升沉计算区间及左右边界条件计算区间；
[0007]采集升沉计算区间
、
左右边界条件计算区间内的惯性导航系统测量的等效垂向加速度，对其进行滤波处理得到垂向加速度正向滤波输出数列；
[0008]垂向加速度正向滤波输出数列进行倒序后进行高通滤波，得到垂向加速度逆向滤波输出数列；
[0009]垂向加速度逆向滤波输出数列进行倒序后得到按时间正向排列的数列即为等效垂向加速度的滤波结果；
[0010]对滤波结果进行二次积分，得到升沉计算区间内对应各点的升沉测量结果
。
[0011]进一步地，所述确定延迟时间
、
升沉计算区间及左右边界条件计算区间的方法为
[0012]令所述计算窗口时间宽度为
Δ
T
；
[0013]延迟时间设定为
Δ
T/2
；
[0014]升沉计算区间为
[T
k
‑1‑
Δ
T2,T
k
‑1+
Δ
T2]，
T
k
表示当前时刻，
T
k
＝
T
k
‑1+
Δ
T
；
[0015]左右边界条件计算区间分别为
[T
k
‑2,T
k
‑1‑
Δ
T/2]和
[T
k
‑1+
Δ
T/2,T
k
]，
T
k
‑1＝
T
k
‑2+
Δ
T。
进一步地，所述垂向加速度正向滤波输出数列为
[0016]y1(n)
＝
h(n)*a(n)
[0017][0018]其中，
a(n)
为等效垂向加速度计测量值，
n
为离散化的测量时间序号，
n
＝
1,2,...
，
h(n)
为高通滤波器，
*
表示卷积计算，
N
为等效垂向加速度数据序列总长度
。
[0019]进一步地，所述升沉测量结果计算公式如下
[0020][0021]其中，
y(n)
是等效垂向加速度计测量值的滤波结果
。
[0022]进一步地，还包括验证等效垂向加速度计测量值
a(n)
与滤波结果
y(n)
之间是否无相位差的方法，
[0023]对
a(n)、h(n)
和
y(n)
进行傅里叶变换，得到
A(e
j
ω
)、H(e
j
ω
)
和
Y(e
j
ω
)
；
[0024]依据滤波过程计算
[0025]Y1(e
j
ω
)
＝
H(e
j
ω
)A(e
j
ω
)
[0026]Y2(e
j
ω
)
＝
e
‑
j
ω
(N
‑
1)
Y1(e
‑
j
ω
)
[0027]Y3(e
j
ω
)
＝
H(e
j
ω
)Y2(e
j
ω
)
[0028]Y(e
j
ω
)
＝
e
‑
j
ω
(N
‑
1)
Y3(e
‑
j
ω
)
[0029]其中，
ω
为频率角速度，
N
为等效垂向加速度数据序列总长度；
[0030]推导得到
[0031]Y(e
j
ω
)
＝
H(e
‑
j
ω
)H(e
j
ω
)A(e
j
ω
)
[0032]确定
Y(e
j
ω
)
与
A(e
j
ω
)
之间无相位差
。
[0033]进一步地，所述零相位升沉测量方法用于水面舰船或者水下航行器
。
[0034]本专利技术还提供了一种零相位升沉测量系统，包括
[0035]数据分段模块，用于根据使用要求对数据序列分段，以及根据延迟时间和边界确定升沉计算区间
、
边界条件计算区间；
[0036]本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.
一种零相位升沉测量方法，其特征在于，包括如下步骤依据使用要求设置计算窗口时间宽度，将数据序列按计算窗口时间宽度分段；确定延迟时间
、
升沉计算区间及左右边界条件计算区间；采集升沉计算区间
、
左右边界条件计算区间内的惯性导航系统测量的等效垂向加速度，对其进行滤波处理得到垂向加速度正向滤波输出数列；垂向加速度正向滤波输出数列进行倒序后进行高通滤波，得到垂向加速度逆向滤波输出数列；垂向加速度逆向滤波输出数列进行倒序后得到按时间正向排列的数列即为等效垂向加速度的滤波结果；对滤波结果进行二次积分，得到升沉计算区间内对应各点的升沉测量结果
。2.
根据权利要求1所述的零相位升沉测量方法，其特征在于，所述确定延迟时间
、
升沉计算区间及左右边界条件计算区间的方法为令所述计算窗口时间宽度为
Δ
T
；延迟时间设定为
Δ
T/2
；升沉计算区间为
[T
k
‑1‑
Δ
T/2,T
k
‑1+
Δ
T/2]
，
T
k
表示当前时刻，
T
k
＝
T
k
‑1+
Δ
T
；左右边界条件计算区间分别为
[T
k
‑2,T
k
‑1‑
Δ
T/2]
和
[T
k
‑1+
Δ
T/2,T
k
]
，
T
k
‑1＝
T
k
‑2+
Δ
T。3.
根据权利要求2所述的零相位升沉测量方法，其特征在于，所述垂向加速度正向滤波输出数列为
y1(n)
＝
h(n)*a(n)
其中，
a(n)
为等效垂向加速度计测量值，
n
为离散化的测量时间序号，
n
＝
1,2,...
，
h(n)
为高通滤波器，
*
表示卷积计算，
N
为等效垂向加速度数据序列总长度
。4.
根据权利要求3所述的零相位升沉测量方法，其特征在于，所述升沉测量结果计算公式如下其中，
y(n)
是等效垂向加速度计测量值的滤波结果
。5.
根据权利要求4所述的零相位升沉测量方法，其特征在...

【专利技术属性】
技术研发人员：周祖洋，庄广琛，郭玉胜，宋雅兰，邓继权，
申请(专利权)人：北京自动化控制设备研究所，
类型：发明
国别省市：