发明名称 一种用于识别振动声调制机理的信号解调方法
摘要 本发明属于无损检测领域,公开了一种用于识别振动声调制机理的信号解调方法,包括:设定所需要WPD解调的等级j;对原始信号进行小波包解调;将解调后的2j个小波包系数进行重构;通过小波包能量比选取特征信息;将步骤四中保留下来的子带信号全部累加,并对累加后的信号进行连续小波变换,得到其特征信息的时频分布图;根据调频信号和调幅信号在时频图中表现的不同特征,判断调制信号的调制类型;根据调制类型确定调制度。本发明所述方法能够有效地、自动地获取含噪声信号的特征信息,抗干扰能力强;通过连续小波变换时频分析,能够非常清晰地识别信号的调制类型,并进一步识别调制度。
申请公布号 CN103472141A 申请公布日期 2013.12.25
申请号 CN201310406433.2 申请日期 2013.09.09
申请人 北京工业大学 发明人 焦敬品;谢辉;宋国荣;吴斌;何存富
分类号 G01N29/36(2006.01)I;G01N29/44(2006.01)I 主分类号 G01N29/36(2006.01)I
代理机构 北京思海天达知识产权代理有限公司 11203 代理人 张慧
主权项 1.一种用于识别振动声调制机理的信号解调方法,其特征在于能够有效识别调制信号的调制类型和调制度,抗干扰能力强,包括以下步骤:步骤一,根据原始信号的特征和解调的目标,设定所需要WPD解调的等级j,4≤j≤8;解调等级越高降噪效果越好,但对电脑硬件要求也越高;步骤二,将原始信号通过小波包解调,小波包函数<img file="FDA00003791534800011.GIF" wi="138" he="82" />定义式为:<maths num="0001"><![CDATA[<math><mrow><msubsup><mi>w</mi><mrow><mi>j</mi><mo>,</mo><mi>k</mi></mrow><mi>m</mi></msubsup><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><msup><mn>2</mn><mrow><mi>j</mi><mo>/</mo><mn>2</mn></mrow></msup><msup><mi>w</mi><mi>m</mi></msup><mrow><mo>(</mo><msup><mn>2</mn><mi>j</mi></msup><mi>t</mi><mo>-</mo><mi>k</mi><mo>)</mo></mrow></mrow></math>]]></maths>式中,j=1,2,…,为尺度参数,k=1,2,…,N/2<sup>j</sup>,为平移因子,N为原始信号长度,m=0,1,…2<sup>j</sup>-1,为振荡因子;步骤三,将解调后的2j个小波包系数<img file="FDA00003791534800013.GIF" wi="502" he="145" />进行重构,得到长度与原始信号相同的信号;步骤四,通过小波包能量比选取特征信息,方法如下:首先设定能量比阈值,计算出各个子带信号的能量比,保留高于能量比阈值的子带信号,滤去低于能量比阈值的子带信号;各个子带信号的能量比<img file="FDA00003791534800014.GIF" wi="58" he="79" />为:<maths num="0002"><![CDATA[<math><mrow><msubsup><mi>L</mi><mi>j</mi><mi>m</mi></msubsup><mo>=</mo><mfrac><mrow><mo>&Integral;</mo><mrow><mo>(</mo><mi>x</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>-</mo><msubsup><mi>p</mi><mi>j</mi><mi>m</mi></msubsup><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>)</mo></mrow><mi>dt</mi></mrow><mrow><msup><mrow><mo>&Integral;</mo><mi>x</mi></mrow><mn>2</mn></msup><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mi>dt</mi></mrow></mfrac></mrow></math>]]></maths>步骤五,将步骤四中保留下来的子带信号全部累加,并对累加后的信号进行连续小波变换,得到其特征信息的时频分布图;连续小波变换把原来的一维信号变为二维信号,其定义式为:<maths num="0003"><![CDATA[<math><mrow><mrow><mo>(</mo><msub><mi>W</mi><mi>&Psi;</mi></msub><mi>f</mi><mo>)</mo></mrow><mrow><mo>(</mo><mi>a</mi><mo>,</mo><mi>b</mi><mo>)</mo></mrow><mo>=</mo><msup><mrow><mo>|</mo><mi>a</mi><mo>|</mo></mrow><mrow><mo>-</mo><mn>1</mn><mo>/</mo><mn>2</mn></mrow></msup><munderover><mo>&Integral;</mo><mrow><mo>-</mo><mo>&infin;</mo></mrow><mo>&infin;</mo></munderover><mi>f</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mover><mrow><mi>&Psi;</mi><mrow><mo>(</mo><mfrac><mrow><mi>t</mi><mo>-</mo><mi>b</mi></mrow><mi>a</mi></mfrac><mo>)</mo></mrow></mrow><mo>&OverBar;</mo></mover><mi>dt</mi></mrow></math>]]></maths>式中,a为伸缩因子,b为平移因子,Ψ(x)称为基本小波或母小波,<img file="FDA00003791534800017.GIF" wi="213" he="147" />是母小波经平移和伸缩后产生的一族函数,称为小波基函数,<img file="FDA00003791534800018.GIF" wi="202" he="156" />是<img file="FDA00003791534800019.GIF" wi="208" he="149" />的共轭运算,W<sub>Ψ</sub>f可认为是原始信号f(t)与一族小波基函数的内积;步骤六,根据调频信号和调幅信号在视频图中表现的不同特征,判断调制信号的调制类型;判断方法如下:如果时频图中信号的能量沿时间轴分布不均匀,调制类型为幅度调制;否则,调制类型为频率调制;步骤七,根据调制类型确定调制度。
地址 100124 北京市朝阳区平乐园100号