发明名称 基于内禀波形匹配的端点效应抑制方法
摘要 该方法公开了一种基于内禀波形匹配的端点效应抑制方法,属于信号分析、故障诊断领域,特别是基于波形匹配的端点效应抑制方法。该发明为对波形两端分别进行延拓,首先对接收波形进行消去中心趋势,然后选取末端的一段波形,与消除中心趋势的波形中其它段波形相匹配;如果匹配成功,则选取匹配成功的该段原始波形对接收波形进行延拓,若匹配不成功则用末端的一段波形进行延拓。从而具有准确性高、稳定性好的效果。
申请公布号 CN104182617B 申请公布日期 2017.04.12
申请号 CN201410370393.5 申请日期 2014.07.30
申请人 电子科技大学 发明人 罗大双;王科盛
分类号 G06F17/50(2006.01)I 主分类号 G06F17/50(2006.01)I
代理机构 电子科技大学专利中心 51203 代理人 张杨
主权项 基于内禀波形匹配的端点效应抑制方法,该方法包括对接收波形左右两端的延拓计算,首先对左端波形的延拓计算过程包括:步骤1、接收到信号,设该信号波形中所有极大值点序列构成的集合为M={M<sub>1</sub>,M<sub>2</sub>,…,M<sub>i</sub>,M<sub>i+1</sub>,…,M<sub>m</sub>},设所有极大值点的相邻极小值点序列构成的集合为N={N<sub>1</sub>,N<sub>2</sub>,…,N<sub>i</sub>,N<sub>i+1</sub>,…,N<sub>m</sub>};设信号的起始点为S<sub>1</sub>(x<sub>s1</sub>,y<sub>s1</sub>),第一个极值点为极大值点,则第一段子波为s<sub>1</sub>为以M<sub>1</sub>为起点向右包含a个极大值点和a个极小值点的波段;则第i段子波起始点为S<sub>i</sub>(x<sub>si</sub>,y<sub>si</sub>),第i段子波为s<sub>i</sub>为以M<sub>i</sub>为起点向右包含a个极大值点和a个极小值点的波段;设极大值点的坐标为M<sub>i</sub>=(x<sub>mi</sub>,y<sub>mi</sub>),极小值点的坐标为N<sub>i</sub>=(x<sub>ni</sub>,y<sub>ni</sub>);根据相似波形的时间尺度比例,计算第i段子波的起始点S<sub>i</sub>的坐标x<sub>si</sub>,由于得到的x<sub>si</sub>点可能不在采样点上,这时采用四舍五入取整;所述第i段子波的起始点S<sub>i</sub>的坐标x<sub>si</sub>计算公式为:<maths num="0001"><math><![CDATA[<mrow><msub><mi>x</mi><mrow><mi>s</mi><mi>i</mi></mrow></msub><mo>=</mo><mi>r</mi><mi>o</mi><mi>u</mi><mi>n</mi><mi>d</mi><mrow><mo>(</mo><mrow><msub><mi>x</mi><mrow><mi>m</mi><mi>i</mi></mrow></msub><mo>-</mo><mfrac><mrow><mo>(</mo><msub><mi>x</mi><mrow><mi>m</mi><mn>1</mn></mrow></msub><mo>-</mo><msub><mi>x</mi><mrow><mi>s</mi><mn>1</mn></mrow></msub><mo>)</mo><mo>(</mo><msub><mi>x</mi><mrow><mi>m</mi><mrow><mo>(</mo><mi>i</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow></mrow></msub><mo>-</mo><msub><mi>x</mi><mrow><mi>m</mi><mi>i</mi></mrow></msub><mo>)</mo></mrow><mrow><msub><mi>x</mi><mrow><mi>m</mi><mn>2</mn></mrow></msub><mo>-</mo><msub><mi>x</mi><mrow><mi>m</mi><mn>1</mn></mrow></msub></mrow></mfrac></mrow><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001147660700000011.GIF" wi="854" he="202" /></maths>其中,x<sub>si</sub>是点S<sub>i</sub>的坐标,round(·)表示四舍五入取整;步骤2、用曲线l<sub>m1</sub>连接s<sub>1</sub>的所有极大值点,用曲线l<sub>n1</sub>连接s<sub>1</sub>的所有极小值点,计算出l<sub>m1</sub>、l<sub>n1</sub>中心线l<sub>1</sub>,则有<img file="FDA0001147660700000012.GIF" wi="251" he="124" />用曲线l<sub>mi</sub>连接s<sub>i</sub>的所有极大值点,用曲线l<sub>ni</sub>连接s<sub>i</sub>的所有极小值点,计算出l<sub>mi</sub>、l<sub>ni</sub>的中心线l<sub>i</sub>,则有<img file="FDA0001147660700000013.GIF" wi="252" he="119" />将s<sub>1</sub>减去l<sub>1</sub>得到一个消去中心趋势的特征子波,记为<img file="FDA0001147660700000014.GIF" wi="89" he="75" />将s<sub>i</sub>减去l<sub>i</sub>得到一个消去中心趋势的匹配子波<img file="FDA0001147660700000015.GIF" wi="89" he="70" />步骤3、求取特征子波<img file="FDA0001147660700000016.GIF" wi="61" he="77" />与各匹配子波<img file="FDA0001147660700000017.GIF" wi="63" he="71" />的匹配度,找出其中最大的匹配度并将其与事先设定的阈值比较,若大于或等于阈值则匹配成功,转到步骤4,若匹配不成功则转到步骤5;步骤4、在输入的信号波段中将最佳匹配子波的起始点S<sub>i</sub>左边的包含k个极大值和k个极小值的波段作为延拓样本,消去该样本的中心趋势,结果记为s<sup>α</sup>,将s<sup>α</sup>加上步骤2中的l<sub>1</sub>,如果l<sub>1</sub>长度不够,则将l<sub>1</sub>延长,得到延拓波形S;将波形S延拓到信号端点S<sub>1</sub>左边;步骤5、如果ρ&lt;β,则表示匹配不成功,此时直接用s<sub>1</sub>进行镜像延拓,如果数据长度不够,则重复延拓该段波形,即完成了数据左端的内禀波形延拓;步骤6、按照步骤1到步骤5的相同方法对输入波形右端进行延拓。
地址 611731 四川省成都市高新区(西区)西源大道2006号