一种舰艇水平姿态摇摆周期的计算方法技术

本发明专利技术涉及一种舰艇水平姿态摇摆周期的计算方法，摇摆周期选择峰值周期，包括如下步骤：步骤1、面向用户实际需求的水平姿态摇摆周期边界值确定；步骤2、对有效的水平姿态角数分别据进行曲线拟合和提取极值点，其中水平姿态角数包括舰艇动态纵摇角数据和舰艇动态横摇角；步骤3、舰艇水平姿态摇摆周期量化计算；步骤4、面向用户实际需求的舰艇水平姿态摇摆周期定量显示。本计算方法能够为用户提供在线计算的舰艇最近时刻摇摆周期及历时摇摆周期。算的舰艇最近时刻摇摆周期及历时摇摆周期。算的舰艇最近时刻摇摆周期及历时摇摆周期。

【技术实现步骤摘要】
一种舰艇水平姿态摇摆周期的计算方法


[0001]本专利技术属于船舰艇装备
，具体涉及一种舰艇水平姿态摇摆周期的计算方法。

技术介绍

[0002]目前，海军舰艇装备的惯性导航类设备(包括基于静电陀螺仪、液浮陀螺仪、挠性陀螺仪、激光陀螺仪、光纤陀螺仪等的惯性导航设备、航姿设备、罗经设备等)均输出纵摇角、横摇角信息，但并未计算纵摇角、横摇角对应的往复周期值。但舰艇官兵在固定翼舰载机、舰载直升机出动回收等场景下，亟需了解舰艇摇摆周期，以对舰载机着舰等动作及时做出决策，以提高任务安全性等任务特性。

技术实现思路

[0003]本专利技术的目的是在于克服现有技术的不足之处，提供一种舰艇水平姿态摇摆周期的计算方法。
[0004]本专利技术的上述目的通过如下技术方案来实现：
[0005]一种舰艇水平姿态摇摆周期的计算方法，其特征在于：摇摆周期选择峰值周期，包括如下步骤：
[0006]步骤1、面向用户实际需求的水平姿态摇摆周期边界值确定；
[0007]步骤2、对有效的水平姿态角数据分别进行曲线拟合和提取极值点，其中水平姿态角数据包括舰艇动态纵摇角数据和舰艇动态横摇角数据；
[0008]步骤3、舰艇水平姿态摇摆周期量化计算；
[0009]步骤4、面向用户实际需求的舰艇水平姿态摇摆周期定量显示。
[0010]进一步的：步骤1中进行水平姿态摇摆周期边界值确定，边界条件包括：惯导摇摆角采样频率不低于10Hz；惯导摇摆角最大测量误差小于2角分；舰船摇摆角度不小于30角分。
[0011]进一步的：步骤2包括：
[0012]2.1、取1s数据采用3次多项式进行拟合，每次拟合后采用滑动窗方式滑动0.1s；设拟合多项式最高次为n，则多项式可表达为：
[0013][0014]其中，y为姿态角测量值，t为时间，a
i
为待拟合参数；
[0015]将惯导输出时刻t
k
及对应摇摆角y
k
带入上式，可得关于多项式系数的方程：
[0016][0017]式中，时刻t
k
是相对于拟合数据段开始时刻的时间间隔，N为待拟合的数据总个数；记
[0018][0019]根据最小二乘算法解得系数为：
[0020]x＝(A
T
A)
‑1A
T
y
…………………………………
(10)
[0021]2.2、定义拟合误差为离散时间处拟合曲线与原曲线最大差值，拟合误差Δy计算公式如下，设定本方案可接受最大拟合误差为1分；
[0022]Δy＝max(Y)
…………………………………
(11)
[0023]式中，Y＝{|y
i
‑
y(t
i
)||i＝1，
……
，N}表示各拟合时刻拟合误差的集合；
[0024]2.3、计算拟合曲线导数和导数判别式，导数计算公式如下：
[0025][0026]判别式计算公式如下：
[0027]Δ＝(2a2)2‑
12a3a1…………………………………
(13)
[0028]2.4、计算极值点：
[0029]当判别式大于0时,计算两个极值点，计算公式如下：
[0030][0031]当a3大于0时，较小极值点为极大值点，较大极值点为极小值点，当a3小于0，较小极值点为极小值点，较大极值点为极大值点；
[0032]当极值点处于区间[0.3,0.7]中，补偿拟合初始时间后，比同类极值点大0.4s时判定该极值点为有效极值点。
[0033]进一步的：步骤3中舰艇水平姿态摇摆周期量化计算的规则为：连续两个极大值点中间包含一个极小值点，且两个极大值点与极小值点的摇摆角幅值差都大于30角分，定义两个极大值点时间间隔作为当前时刻摇摆周期；
[0034]若没有搜索到有效极大值点，则在极大值点列中向前滑动一个极大值点再次进行判断，若距离当前时刻1min内无有效极大值点对，则认为当前舰船平稳，无明显摇摆。
[0035]进一步的：步骤4中舰艇水平姿态摇摆周期定量显示包括：一、30s内舰艇水平姿态平均摇摆周期；二、60s内舰艇水平姿态平均摇摆周期；三、5min内舰艇水平姿态平均摇摆周期。
[0036]本专利技术具有的优点和积极效果：
[0037]本专利技术舰艇摇摆周期的计算方法，主要以舰艇动态纵摇角、横摇角参数为输入数据，通过滑窗拟合等方法建立摇摆周期计算方案，运用极值判断等方法进行赋值，最终实现舰艇摇摆周期动态确定，能够为用户提供在线计算的舰艇最近时刻摇摆周期及历时摇摆周期。
附图说明
[0038]图1是本专利技术舰船摇摆峰值计算流程图。
具体实施方式
[0039]以下结合附图并通过实施例对本专利技术的结构作进一步说明。需要说明的是本实施例是叙述性的，而不是限定性的。
[0040]一种舰艇水平姿态摇摆周期的计算方法，主要包括：(1)面向用户实际需求的水平姿态摇摆周期边界值确定；(2)对有效的水平姿态角数分别据进行曲线拟合和提取极值点，其中水平姿态角数包括舰艇动态纵摇角数据和舰艇动态横摇角；(3)舰艇水平姿态摇摆周期量化计算；(4)面向用户实际需求的舰艇水平姿态摇摆周期定量显示。具体的：
[0041]在海浪等外界激励作用下，舰船绕平衡位置作往复摇摆，由于外界激励作用的复杂性，使得舰船的摇摆运动曲线呈现复杂的形态，具有丰富的频谱成分，因此，舰船运动周期的定义也有多种方式，如峰值周期、谱峰周期、谱平均周期等。峰值周期定义为两次摇摆角峰值对应的时间间隔，体现摇摆运动的短时特性，反映上一个摇摆的周期，具有一定的延时特性；谱峰周期定义为一段时间内摇摆频谱成分中幅值最大的周期；谱平均周期是摇摆频谱的一种综合评价，反映一段时间内摇摆能量主要分布的周期。
[0042]本专利技术主要以反映短时摇摆周期特性为目的，选择峰值周期作为舰船摇摆周期运动的计算方式。
[0043]上述面向用户实际需求的水平姿态摇摆周期边界值确定为：在舰艇摇摆幅度过小时，计算摇摆周期的需求不存在。因此，本专利技术的应用条件如下：
[0044]1、惯导摇摆角采样频率不低于10Hz；
[0045]2、惯导摇摆角最大测量误差小于2角分；
[0046]3、舰船摇摆角度不小于30角分，低于30角分时认为舰船处于平稳状态。
[0047]上述对有效的水平姿态角数分别据进行曲线拟合和提取极值点的计算过程如下：
[0048]1、取1s数据采用3次多项式进行拟合，每次拟合后采用滑动窗方式滑动0.1s；设拟合多项式最高次为n，则多项式可表达为：
[0049][0050]其中，y为姿态角测量值，t为时间，ai为待拟合参数。
[0051]将惯导输出时刻t
k
及对应摇摆角y
k
带入上式，可得关于多项式系数的方程：<本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种舰艇水平姿态摇摆周期的计算方法，其特征在于：摇摆周期选择峰值周期，包括如下步骤：步骤1、面向用户实际需求的水平姿态摇摆周期边界值确定；步骤2、对有效的水平姿态角数据分别进行曲线拟合和提取极值点，其中水平姿态角数据包括舰艇动态纵摇角数据和舰艇动态横摇角数据；步骤3、舰艇水平姿态摇摆周期量化计算；步骤4、面向用户实际需求的舰艇水平姿态摇摆周期定量显示。2.根据权利要求1所述的舰艇水平姿态摇摆周期的计算方法，其特征在于：步骤1中进行水平姿态摇摆周期边界值确定，边界条件包括：惯导摇摆角采样频率不低于10Hz；惯导摇摆角最大测量误差小于2角分；舰船摇摆角度不小于30角分。3.根据权利要求2所述的舰艇水平姿态摇摆周期的计算方法，其特征在于：步骤2包括：2.1、取1s数据采用3次多项式进行拟合，每次拟合后采用滑动窗方式滑动0.1s；设拟合多项式最高次为n，则多项式可表达为：其中，y为姿态角测量值，t为时间，a
i
为待拟合参数；将惯导输出时刻t
k
及对应摇摆角y
k
带入上式，可得关于多项式系数的方程：式中，时刻t
k
是相对于拟合数据段开始时刻的时间间隔，N为待拟合的数据总个数；记根据最小二乘算法解得系数为：x＝(A
T
A)
‑1A
T
y
…………………………………
(3)2.2、定义拟合误差为离散时间处拟合曲线与原曲线最大差值，拟合误差Δy计算公式如下，设定本方案可接受最大拟合误...

【专利技术属性】
技术研发人员：邓福建，常红军，韩旭，
申请(专利权)人：中国船舶重工集团公司第七零七研究所，
类型：发明
国别省市：