发明名称 基于小波特征值容差阈值随机统计模拟的故障判别方法
摘要 本发明公开了一种基于小波特征值容差阈值随机统计模拟的故障判别方法,它包括利用信号的小波分析判别故障、利用信号容差阈值判别故障、利用小波分析特征值容差阈值判别故障、利用小波分析敏感特征值阈值判别故障等步骤。本发明采用随机统计模拟的方法得到信号小波分析特征值阈值和敏感特征值阈值,并应用特征值阈值和敏感特征值阈值判断信号是否为故障信号,该方法不需要大量故障信号和正常信号样本即可判断,解决了采用神经网络和支持向量机等众多方法因无法获得训练样本而导致的无法进行故障判别的难题,同时该方法节省了大量的时间,方法简单、准确性高。
申请公布号 CN106292622A 申请公布日期 2017.01.04
申请号 CN201610578532.2 申请日期 2016.07.21
申请人 中国人民解放军军械工程学院 发明人 马彦恒;李刚;尹园威;刘新海
分类号 G05B23/02(2006.01)I 主分类号 G05B23/02(2006.01)I
代理机构 石家庄科诚专利事务所 13113 代理人 张红卫
主权项 一种基于小波特征值容差阈值随机统计模拟的故障判别方法,其特征在于它按照如下的步骤顺序进行:一、利用信号的小波分析判别故障采用小波分析方法,将被测信号和样本信号(未发生故障时的正常信号)进行时间或空间频率的局部化分析,通过伸缩平移运算对信号(函数)逐步进行多尺度细化,提取样本信号和被测信号的多种频率,进行显示比对,并计算各层小波分析的弧距,计算高频分量特征值并进行比对;如果样本信号和被测信号小波分析出的特征值差异较大,则可判断被测装备发生了故障;如果特征值差异较小,则需要依次利用下述步骤的信号容差阈值、特征值阈值和敏感特征值阈值来进行综合判别;二、利用信号容差阈值判别故障信号容差均值的算法如公式(1),<maths num="0001"><math><![CDATA[<mrow><msub><mover><mi>&Delta;</mi><mo>&OverBar;</mo></mover><mrow><mi>x</mi><mi>y</mi></mrow></msub><mo>=</mo><mrow><mo>(</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><mo>|</mo><msub><mi>x</mi><mi>i</mi></msub><mo>-</mo><msub><mi>y</mi><mi>i</mi></msub><mo>|</mo><mo>)</mo></mrow><mo>/</mo><mi>n</mi><mo>,</mo><mi>i</mi><mo>=</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>,</mo><mn>3</mn><mo>,</mo><mo>...</mo><mo>,</mo><mi>n</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001055874030000011.GIF" wi="1557" he="126" /></maths>公式(1)中,<img file="FDA0001055874030000012.GIF" wi="68" he="71" />为信号的容差均值,x<sub>i</sub>为样本信号x的第i个数据,y<sub>i</sub>为被测信号y的第i个数据,n为数据采样点数;给定<img file="FDA0001055874030000013.GIF" wi="72" he="70" />一个限定值<img file="FDA0001055874030000014.GIF" wi="90" he="63" /><img file="FDA0001055874030000015.GIF" wi="62" he="63" />为信号的容差阈值,当<img file="FDA0001055874030000016.GIF" wi="191" he="70" />时,则判定被测信号y为故障信号,否则为正常信号,称<img file="FDA0001055874030000017.GIF" wi="62" he="63" />为信号的容差阈值;三、利用小波分析特征值容差阈值判别故障利用小波分析特征值容差阈值按照如下步骤判别故障:步骤①:根据步骤二中输入的信号容差阈值,以样本信号数据作为基准,将步骤二中输入的信号容差阈值作为信号容差均值期望,仿真模拟出若干组被测信号数据;步骤②:对步骤①获得的仿真模拟信号数据每组分别进行小波分析,求出每组仿真模拟信号数据各自的特征参数的值;步骤③:将所有组仿真模拟信号小波分析同一个特征参数的值相加,除以仿真模拟信号组数,求得此特征参数的均值,用此方法求出所有特征参数的均值,特征参数均值的算法如公式(2),<img file="FDA0001055874030000021.GIF" wi="1462" he="126" />公式(2)中,<img file="FDA0001055874030000022.GIF" wi="47" he="59" />为所有组仿真模拟信号小波分析第j个特征参数的均值,r为小波分析特征参数的个数,m为仿真模拟信号数据的组数,<img file="FDA0001055874030000023.GIF" wi="67" he="54" />为第i组仿真模拟信号数据小波分析的第j个特征参数值;步骤④:将步骤③获得的全部特征参数的均值与样本信号小波分析的同一个特征参数的值相减取绝对值并求和,获得特征值容差阈值<img file="FDA0001055874030000024.GIF" wi="89" he="78" />特征值容差阈值<img file="FDA0001055874030000025.GIF" wi="67" he="78" />的算法如公式(3),<img file="FDA0001055874030000026.GIF" wi="1190" he="127" />公式(3)中,<img file="FDA0001055874030000027.GIF" wi="61" he="78" />为小波分析特征值容差阈值,<img file="FDA0001055874030000028.GIF" wi="51" he="55" />为样本信号第j个特征参数的值;步骤⑤:对实测的被测信号进行小波分析,将其小波分析各特征值分别与样本信号小波分析对应的特征值相减取绝对值并求和,得到实际特征值容差<img file="FDA0001055874030000029.GIF" wi="91" he="78" />如果<img file="FDA00010558740300000210.GIF" wi="219" he="85" />则判断被测信号为故障信号,否则为正常信号,实际特征值容差<img file="FDA00010558740300000211.GIF" wi="67" he="82" />的算法如公式(4),<img file="FDA00010558740300000212.GIF" wi="1174" he="127" />公式(4)中,<img file="FDA00010558740300000213.GIF" wi="69" he="86" />为实际特征值容差,<img file="FDA00010558740300000214.GIF" wi="49" he="63" />为实测的被测信号小波分析第j个特征参数的值;四、利用小波分析敏感特征值阈值判别故障利用小波分析敏感特征值阈值对故障进行判别,具体步骤如下:步骤(a):将步骤②中获得的每组仿真模拟信号数据的小波分析特征值按照小波分析的层进行分组,每一层小波分析对应一组特征值;步骤(b):以样本信号数据小波分析每一层对应的特征值为基准,计算所有仿真模拟信号数据小波分析同一层对应的特征值与这些基准特征值的容差均值,算法如公式(5),<img file="FDA0001055874030000031.GIF" wi="1255" he="136" />公式(5)中,<img file="FDA0001055874030000032.GIF" wi="97" he="86" />为所有仿真模拟信号数据与样本信号数据第h层小波分析特征值容差均值,h为小波分析的层号,r<sub>h</sub>为第h层小波分析的特征值个数,<img file="FDA0001055874030000033.GIF" wi="93" he="70" />为第i组仿真模拟信号数据第h层小波分析的第j个特征值,<img file="FDA0001055874030000034.GIF" wi="73" he="54" />为样本信号数据第h层小波分析第j个特征值,k为总共进行的小波分析的层数;步骤(c):从第1层到第k层,对计算出的<img file="FDA0001055874030000035.GIF" wi="96" he="86" />值进行大小排列,选取值最大的前p个对应的小波分析层的特征值组作为敏感特征值组;步骤(d):将所有组仿真模拟信号小波分析中,步骤(c)所选出的p个敏感特征值组里的同一个特征参数的值相加,除以仿真模拟信号的组数,求得此敏感特征值的均值,用此方法求出这p个特征值组里所有特征值的均值,敏感特征值均值的算法如公式(6),<img file="FDA0001055874030000036.GIF" wi="1510" he="126" />公式(6)中,<img file="FDA0001055874030000037.GIF" wi="78" he="63" />为所有组仿真模拟信号p个敏感特征值组中第j个敏感特征值的均值,r<sub>p</sub>为p个敏感特征值组中所有敏感特征参数的个数,<img file="FDA0001055874030000038.GIF" wi="102" he="55" />为第i组仿真模拟信号数据小波分析的p个敏感特征值组中第j个特征值;步骤(e):将步骤(d)获得的全部敏感特征值的均值与样本信号对应小波分析敏感特征值相减取绝对值并求和,获得敏感特征值容差阈值<img file="FDA0001055874030000041.GIF" wi="114" he="87" />敏感特征值容差阈值<img file="FDA0001055874030000042.GIF" wi="91" he="90" />的算法如公式(7),<img file="FDA0001055874030000043.GIF" wi="1174" he="135" />公式(7)中,<img file="FDA0001055874030000044.GIF" wi="87" he="86" />为小波分析敏感特征值容差阈值;步骤(f):对实测的被测信号进行小波分析,将其p个敏感特征值组中的所有敏感特征值分别与样本信号小波分析的同一个敏感特征参数的值相减取绝对值并求和,得到实际敏感特征值容差<img file="FDA0001055874030000045.GIF" wi="123" he="87" />如果<img file="FDA0001055874030000046.GIF" wi="275" he="87" />则判断被测信号为故障信号,否则为正常信号,实际敏感特征值容差<img file="FDA0001055874030000047.GIF" wi="99" he="87" />的算法如公式(8),<img file="FDA0001055874030000048.GIF" wi="1245" he="135" />公式(8)中,<img file="FDA0001055874030000049.GIF" wi="97" he="83" />为实际敏感特征值容差,<img file="FDA00010558740300000410.GIF" wi="77" he="70" />为实测的被测信号小波分析p个敏感特征值组中第j个敏感特征参数的值。
地址 050011 河北省石家庄市和平西路97号