APP下载

弹上电气系统大基数时序台阶量自动判读方法

2022-10-12刘洪政普承恩赖玉敏崔允红

弹箭与制导学报 2022年4期
关键词:测试数据幅值差分

刘洪政,普承恩,胡 洲,赖玉敏,曹 海,崔允红,黄 松

(四川航天系统工程研究所,成都 610100)

0 引言

随着武器装备的发展,弹上电气系统试验测试任务呈现出数据量大、测试项目杂、参数种类多的特点,测试数据判读是保证电气系统质量性能的重要一环,有助于问题的排查,对故障进行定位。目前在弹上电气系统系统测试工作方面,主要依靠人工进行数据判读及复查,不仅工作量大、效率低,而且容易造成人为判读遗漏和偏差。因此,迫切需要在弹上电气系统数据判读过程中引入自动判读,提高判读效率和判读质量。

针对人工进行测试数据判读的不足,近年来有学者进行了相关的测试数据自动判读研究,取得了显著成果。陈二雷等针对制导弹药的测试数据人工判读时间长,误判概率大等不足,根据测试数据的特点对各物理量参数进行分类并设计了判读规则,依据该判读规则可对物理参量是否合格进行自动判断,最后将判读结果生成测试报告,提高了测试效率和判读的准确性。为提高数据判读的实时性、完整性、准确性,韩璐分析试飞测试数据的特点,建立通用判读规则,设计试飞数据自动判读系统,具备快速判读、详细判读、异常诊断的功能,很大程度提高对飞行结果的综合分析评估能力。陈策等为了提高运载火箭遥测参数起始电平的判读效率和覆盖性,更好地满足高密度测试的需求,基于LabVIEW平台和ACCESS数据库设计开发了一款遥测参数起始电平自动判读软件,该软件融合了数据库模型、自动判读、包络分析、实时显示、故障报警等多个功能,可实现遥测参数全波道数据的自动测试。惠国娟等从测试数据的复杂性出发,对测试数据的判读比对和统计分析方法进行了研究,提高了产品全寿命周期产生测试数据分析和判读能力。此外,曹永梅、吴婧等针对卫星控制系统测试数据的特点和判读难点,开展了测试数据判读规则、判读方法和判读系统研究,提高了测试效率和判读的准确性。

以上多种判读方法的研究,提高了测试数据的判读效率和判读准确性。然而,针对时序台阶量尤其是大基数时序台阶量的研究依旧较少。周辉峰等提出了中值滤波算法实现了台阶参数和脉冲参数的滤波平滑处理,并采用双边多点阈值判断方法与符号判读方法相结合的台阶识别算法,实现了台阶的自动准确识别,但存在两点不足:第一,该方法对于数据的边界点未作处理,当台阶量出现在边界点处会导致阈值矩阵下标溢出,无法准确识别出台阶量;第二,该方法的阈值取值需要依靠经验获取,因此在大基数台阶量测试数据判读中,当不同量级的台阶幅值同时存在时,取值不合适则会导致台阶量漏判。

因此针对弹上电气系统大基数时序台阶量测试数据开展自动判读方法研究,提出了时序差分矩阵判别和台阶幅值映射两步方法,分别实现台阶量的准确识别和大基数的台阶自动判读,以提高大基数时序台阶量的判读效率和判读准确性。

1 大基数时序台阶量的特点分析

时序台阶量主要由弹上电气系统控制时序决定。实际工程测试中发现,由于某些指令的特殊性,时序台阶量中往往存在控制指令基数大而指令之间数值相差小的台阶量,导致不同量级的指令混合在一起,这类数据通常具有时序性、基数大、跳变复杂的特点,将具有上述特点的台阶量称为大基数时序台阶量。大基数时序台阶量示意如图1所示。

图1 大基数台阶量示意图

由图1可得,大基数时序台阶量测试数据具有以下一些特点:

1)台阶指令是严格按照时序产生的,不会在时间上产生跳跃,因而具有时序性的特点。

2)部分台阶指令幅值较大,而部分台阶的指令幅值较小,由于小量级与大量级台阶同时存在,台阶数量多,且在时序上紧邻,若依靠人工进行大基数数据判读,则很难在数据图中发现量级较小的台阶指令。

3)不同台阶指令的持续时间不同,一些指令持续时间较长,而一些指令持续时间较短,导致了台阶测试数据的跳变复杂,人工判读时,对于持续时间较长的台阶容易判断,而持续时间较短的台阶则容易漏判。

2 系统总体方案设计

大基数时序台阶参数判读的原理是首先准确寻找测试原数据中的每一个台阶量,并对每一个台阶进行标记和保存,然后将台阶量与指令库进行比较,并将比较结果按照幅值映射规则进行映射,将映射结果保存和显示,最后制作判读报表,输出判读结果。

在上述过程中,关键在于大基数时序台阶量的准确识别和台阶幅值正确性的自动判读,根据大基数时序台阶参数判读的原理,将自动判读过程分为3个模块:第1部分为台阶量识别模块;第2部分为幅值映射模块;第3部分为判读结果显示模块。大基数时序台阶量自动判读总体方案如图2所示。

在台阶量识别模块中,主要完成原始测试数据的台阶量识别,保存和记录台阶量的时间、幅值等信息,为后续工作做准备;在台阶幅值映射模块中,主要对所识别到的台阶量,完成台阶参数的正确性判读以及幅值映射任务;在判读结果可视化显示中,主要完成判读结果报表的制作,幅值映射结果的可视化显示任务。

图2 大基数时序台阶量自动判读总体方案

3 判读算法设计

3.1 大基数时序台阶识别算法

为准确识别出每一个时序台阶参数,针对大基数时序台阶参数的特点提出时序差分矩阵判断法,该方法可检测到数据中的所有跳点,并结合判断规则对跳点进行判断,能准确识别出台阶量。

假设一组电气系统的测试数据按照时序排列为{,,,…,},其中为测试数据的长度。大基数台阶识别具体步骤如下:

1)查找跳点

由于台阶量的第一个周期数值相对于前一个周期数值具有跳变的特征,因此要实现台阶量的识别,应该首先查找原始数据中的跳点。查找跳点的方法为:

按照时序遍历原始测试数据,将第(+1)个数据与第个数据进行比较,即+1-,若|+1-|>,则将+1标记为跳跃点,其中=1,2,…,-1;为判断数据+1是否为跳点的条件阈值。

为了避免对大基数台阶量的漏判,阈值的确定应当根据预先装订的指令库进行确定,指令库由一组各不相同的指令值按照时序排列组成,设指令库为:

(1)

式中:为指令的个数;代表产生指令所对应的时刻。则有:

(2)

式中:=2,3,…,;=1,2,…,-1。

2)创建核窗口

当在测试数据中查找到跳点+1后,创建一个核窗口。创建方法是以+1为中心,取其前后对称周期的数据赋值给核窗口,为后续计算时序差分矩阵做准备。则核窗口为:

{,,…,}=
{-(-3)2,…,,+1,+2,…,+(+1)2}

(3)

式中:为核窗口大小,为保证核窗口中心对称,为奇数。的取值决定了后续台阶量识别的准确度其值需要根据台阶量的特点来确定;的大小决定了差分矩阵的大小,若取值太大,将导致差分矩阵维数增大,从而使自动判读速度减慢,因此取5~9。

当核窗口的大小确定下来之后,将源数据按照时序赋值给核窗口,然后进行台阶量判断,判断完成之后,继续寻找下一个跳点,直至遍历完所有测试数据为止。假设检测到+1为跳点,则核窗口工作原理如图3所示:

图3 大小为m的核窗口工作原理

3)核窗口越界处理

从图3中可看出,核窗口的工作原理是将核窗口中心点之后的元素与中心点之前的元素分别进行差分。需要注意的是,当<(+1)2-1以及>-(-1)2-1时,即跳点+1位于测试数据的边界处时,核窗口与源数据无法对齐,将导致核窗口的赋值越界,以=5为例,核窗口越界示意图如图4所示。

图4 m=5的核窗口越界示意图

为了避免越界,核窗口取值时需进行如下处理:

(4)

4)求解时序差分矩阵

时序差分矩阵是根据核窗口大小来进行计算的,计算原理是核窗口中跳点之后的数据与跳点之前的数据分别差分,假设核窗口的时序差分矩阵为,则时序差分矩阵计算如下:

当(+1)2-1<<-(-1)2-1时:

(5)

当<(+1)2-1时,核窗口左越界:

(6)

当>-(-1)2-1时,核窗口右越界:

(7)

5)时序差分矩阵映射

时序差分矩阵是由一系列核窗口数据之间的偏差量组成,为了直观反映出偏差量的正负情况,引入映射函数,将时序差分矩阵中的元素通过映射关系映射到指定域中,映射函数设计如下:

假设自变量为,为阈值条件,映射域为[-1,0,1],则有:

(8)

结合阈值条件将时序差分矩阵通过映射函数映射为差分映射矩阵,其中的元素根据式(8)可表示为:

=()

(9)

式中:=1,2,…,;=1,2,…,(-);为时序差分矩阵的行数。

6)求解映射矩阵判读规则

为便于差分映射矩阵的判读,根据矩阵中的映射值情况,将其中的所有元素求和并设置相应判断规则:

(10)

规则的作用原理是:若差分映射矩阵中值为1的元素个数等于所有元素总个数,则可判定为上跳台阶,=1;若差分映射矩阵中值为-1的元素个数等于所有元素总个数,则可判定为下跳台阶,=-1;其他情况则可判定为不是台阶,=0。完成了台阶量的识别,台阶参数识别算法原理如图5所示。

图5 台阶参数识别算法原理

3.2 台阶幅值映射

台阶量识别完成后,应对台阶的正确性进行判断,台阶正确性判断包含两个部分,一部分是幅值正确性判断,另二部分是时间正确性判断。为了使输出结果直观易判,提出台阶幅值映射方法,将大基数时序台阶参数判读结果按照一定规则进行映射,幅值映射的原理如下:

将识别到的各个台阶参数幅值记为集合,则幅值映射的基本原理就是将中的元素按照对应法则映射到集合当中去,集合称之为映射域,记为:→。映射规则设计如下:

(11)

式中:为系统测试时间;为实测台阶值;为台阶的序号;为第个台阶出现的时间;Δ为出现第个台阶的时间偏差;为第个台阶的期望值。

式(11)的具体原理是:根据预先装订的指令库查找当前实测的台阶值,若在指令库中未找到当前台阶值,则认为台阶幅值异常,将台阶幅值映射为-1;若在指令库中找到当前台阶值且当前时间在期望的时间范围内,则认为台阶幅值正常,将台阶幅值映射为;若在指令库中找到当前台阶值但当前时间不在期望的时间范围内,则认为台阶时间异常,将台阶幅值映射为0。因此,台阶量判读结果为正常时,(,)=;台阶量判读结果为台阶幅值异常时(,)=-1;台阶量判读结果为台阶时间异常时(,)=0。台阶幅值映射原理如图6所示。

图6 台阶幅值映射原理

4 判读算法验证

为对判读算法的有效性进行验证,对某导弹电气系统配电器指令码测试数据开展了自动判读试验验证工作,分别包括台阶识别算法验证和台阶幅值映射算法验证两部分,预先装订的配电器指令如表1所示。

表1 指令表

4.1 台阶识别算法验证

取时序差分矩阵的阈值为=0.8,根据配电器指令码测试数据的特点将核窗口大小=5,Δ=0.1 s。台阶自动识别结果如图7所示。

图7 时序台阶指令识别结果

图7(a)原始测试数据时序台阶指令识别结果;图7(b)为上跳台阶识别结果以及放大图,当原始数据由32 807到35 623变化时,台阶识别标识由0到1突变,在指令无变化区域,标识为0;图7(c)为下跳台阶识别结果以及放大图,当原始数据由35 623到32 836变化时,台阶识别标识由0到-1突变,在指令无变化区域,表识为0;图7(d)为非台阶跳点识别结果以及放大图,当原始数据由35 637到65 535变化以及由65 535到35 637变化时,台阶识别标识为0未发生变化。通过图片数据分析可以得出,所提出的时序差分矩阵法能够准确的识别出了大基数测试数据中的所有台阶量以及台阶类型,并且自动判别出了非台阶量,避免了异常跳点对台阶识别的影响。

4.2 台阶幅值映射验证

在台阶识别的基础上,进行大基数时序台阶量幅值映射验证,分为正常测试数据验证和异常测试数据验证,验证结果如图8所示。

图8 台阶幅值映射结果

表2 异常数据判读报表

正常数据的幅值映射结果如图8(a)所示,21个台阶和21个非零时序指令按照时间顺序依次完成动作,且幅值依次增大。异常数据的幅值映射结果如图8(b)所示,指令随时间不是依次增大,存在2个-1指令(指令异常)和3个0指令(指令正常,时间异常),具体异常情况在后续报表中详细说明。

图8(a)为正常指令映射结果,与表1中的指令数据一致,故不再列出判读报表。图8(b)的判读报表如表2所示,报表内容包括指令序号、时间、指令值、相对时间、台阶类型、映射值、判读结论等。表2中台阶1和台阶13均为指令异常,指令序号显示为“无”,代表此台阶指令不属于指令库中的指令值,因此判定为指令异常;台阶2、台阶6、台阶14的幅值分别对应了指令库中的指令2、指令5、指令12,均属于指令库中的指令,但相应的时间不在指令库对应的时间范围内,因此判定为时间异常。

通过台阶幅值的映射,将大基数指令自动映射到指定域中,使配电器测试数据中的每个指令的正确性都能直观显示,不会由于指令基数较大而影响判读,可解决容易误判和漏判的问题,大幅提高判读效率。

将判读算法用于某型号电气系统综合试验测试。结果表明,在20个时序指令条件下,自动判读时间不大于10 s,人工判读时间不小于5 min,判读效率显著提升。

5 结论

针对弹上电气系统测试数据中产生的大基数时序台阶量的特点,提出了时序差分矩阵判别和幅值映射两步判读方法,该方法可自动准确识别台阶量并使台阶参数按照映射法则自动映射到映射域,使判读结果直观清晰,解决了人工对时序台阶量容易误判和漏判的问题。为验证所提判读方法的正确性,进行了数据判读试验。试验结果表明,该方法有效可行,有效提高了判读效率。

猜你喜欢

测试数据幅值差分
一类分数阶q-差分方程正解的存在性与不存在性(英文)
室温下7050铝合金循环变形研究
一个求非线性差分方程所有多项式解的算法(英)
测试数据管理系统设计与实现
一类caputo分数阶差分方程依赖于参数的正解存在和不存在性
可靠性步进电机细分驱动技术研究
基于烟花爆炸优化算法的测试数据生成方法
基于差分隐私的数据匿名化隐私保护方法
平地机作业负载谱分析
基于自适应粒子群优化算法的测试数据扩增方法