信号检测与估计论文

  • 格式:doc
  • 大小:321.64 KB
  • 文档页数:7

利用自相关-小波阈值变换联合去噪和伪哈密顿量的变尺度Duffing振子来检测微弱周期信号的方法

崔妮婷

【摘要】:Duffing振子微弱周期信号检测方法只对低率参数信号比较敏感,并且系统临界阈值会随待测信号频率改变。文章在分析Duffing非线性动力学系统运动特性基础上,针对Duffing振子微弱信号检测方法存在的上述问题,提出基于自相关-小波阈值变换联合去噪和伪哈密顿量的变尺度Duffing振子弱信号检测方法。该方法对待测信号二次采样,把高频信号变换为一个确定低频参数信号,这样就唯一确定了系统临界阈值;采用自相关和小波变换联合去噪方法,最大程度改善了待测信号信噪比;然后构造Duffing振子伪哈密顿量,能够在只有较少数据量时定量、准确判断有无待测信号。

【关键词】:Duffing振子 尺度变换 联合去噪 伪哈密顿量

【Abstract】:Duffing oscillator signal detection method is more sensitive for low rate parameter

signals,and the system threshold will change with signal frequency. This article based on the

analysis on nonlinear dynamics of Duffing system movement characteristics, in view of the

problems mentioned the above, put forward the method of the scale transformation Duffing

oscillator based on autocorrelation jointed with wavelet threshold transformation and

Pseudo-Hamiltonian.This method samples the signal twice to transform the high frequency signals

to a low frequency parameters, thus the unique state transition threshold was determined; Using

autocorrelation jointed with wavelet threshold transformation to improve the SNR to greatest

degree; And then construct the Duffing oscillator Hamiltonian, which can be used to make

quantitative, accurate judgment for measuring signal when there is only little amount of data .

【Key words】:Duffing oscillator scale transformation united denoising Pseudo-Hamiltonian

高频工程信号经过二次采样变换为低频参数信号,此低频参数信号通过自相关器,抑制部分噪声,同时去掉了待测信号相位信息,对相关信号进行小波阈值变换,去除自相关后的残留噪声。去噪后信号输入Duffing振子检测系统,并设置系统初始状态参数。然后计算Duffing系统平均伪哈密顿量,设定检测系统伪哈密顿量门限阈值。最后由系统输出的伪哈密顿量判断系统状态及是否存在待测信号。

信号检测基本原理框图如下:

尺度变换平均伪哈密顿量T待测信号PH自相关器小波阈值变换Duffing系统门限判决

1.基于Duffing方程微弱信号检测原理

选用连续动力学系统中Duffing振子作为研究对象,Duffing方程标准形式为:

3cosxkxxxrwt (1)

式中:k为阻尼比,r为策动力振幅,w为策动力力角频率。

Duffing振子有如下特点:当k固定时,调整内置策动力幅值r大小,当r =

0时,相点(,)xx最终停在两焦点之一;当0crr时相点在两焦点附近做周期运动,系统处于倍周期分岔状态;crr时为从倍周期分岔到混沌态的过渡态;当cdrrr时系统进入混沌状态;drr时为从混沌状态向周期状态转变的临界混沌态。这种系统从混沌态到大尺度周期状态的跃迁特性被用来检测噪声环境中的微弱信号。实验中设置信号测系统中策动力rcoswt幅值大小为r = rd,此时系统处于临界混沌状态,若加入与策动力同频率的微弱信号,系统会跳变到大尺度周期状态,然后调节内置策动力幅值为rx时系统恰好进入周期状态,则待测信号幅值大小为dxrr。 为提高系统检测模型普适性,实现对任意频率微弱正弦信号的检测,对式(1)变形,令tw,得:

2232cos()xkwxwxwxrww (2)

其中dxdxxwddt,令dxydt,dyyd,则(2)式中加入待测微弱信号并写成状态方程形式为:

35[cos()()]xwyywxxkyrwi (3)

上式中()cos(())sihwwz为外部引入的湮没在噪声中的微弱正弦信号,h为待测正弦信号幅度,w为待测信号与Duffing系统内置周期策动力信号频率差,为待测信号初相位,sz为待测信号中混有的噪声。

对于cos()cosrwhww 令 w,得

cos()cos()cos(())rwhwwFw (4)

其中: 22()2cos()Frrh

(5)

sin()()arctancos()hwrhw (6)

这里 ()F为总策动力幅值,()为总策动力初相角。因为 h << r,所以()的影响可以忽略不计。当0w时,只存在相位差时系统特性如下:

若()dFr系统可能发生相变,此时解得的范围为:

arccosarccos22ddhhrr (7)

当 不在此范围时,()dFr,系统始终处于混沌状态,待测信号不能被检测出来。

2.高频参数待测信号尺度变换

Duffing系统只有在低频参数条件下有较好动态特性和检测效果[10],且Duffing振子检测信号时,不同频率待测信号对应的相变阈值也不同[11],如果每次检测过程都要搜索相变阈值,将大大增加检测复杂度。所以,在处理实际工程信号时,文章对待测信号进行二次采样,即引入变尺度系数R,对待测信号进行频率/时间尺度变换。这种变换并不改变参与计算数据的数值,只是将其在时间尺度上进行压缩或放大。若待测信号角频率为w,其采样频率为fs,则数值计算的步长为1/sdtf。对检测系统引入变尺度系数R相当于将信号的时间间隔增大了R倍,相应的信号角频率被压缩R倍后变为/wR,此时数值计算步长变为/sdtRdtRf。经过二次采样后的信号在通过Duffing系统时可被识别,对该信号进行尺度恢复即可得到原采集信号中的特征信号。

3.自相关去噪

设已知频率待测信号为:()()()xtstnt,()st是周期信号,()nt是均值为零的高斯白噪声,信号自相关输出为:

221()lim()()()()()()TTxxTsssnnsnnRxtxtdtTRRRR (8)

根据互相关函数性质,由于信号()st与噪声()nt不相关,并且噪声平均值为0,得到()()0snnsRR,则有

()()()xxssnnRRR (9)

随着增大,()nnR趋近于零,则对充分大的,可得()()xxssRR。

工程应用中,当待测信号为()()cosntxthwt时,有

2()cos()()2xxaRwnt (10)

式(10)中,()nt是相关信号中混有的噪声。

可见,实际中由于积分时间不可能无限长,噪声只能得到一定程度抑制,另外,自相关器输出信号相位为0,由式(7)可知,由于丢失了信号相位信息,反而方便了检测。

4.小波阈值去噪

小波阈值消噪过程中,信号经过小波变换后,可以认为由信号产生的小波系数包含有信号的重要信息,其幅值大,但数目较少,而噪声对应的小波系数幅值小。因此,通过在不同尺度上选取一合适阈值,并将小于该阈值的小波系数置零,而保留大于该阈值的小波系数,从而使信号中的噪声得到有效抑制。最后进行逆小波变换,得到去噪后的重构信号。

假设离散含噪待测信号可由式(11)表示为

()()()xisini (1)iM (11)

式中12(,,,,,)iMxxxxx为含噪信号,12(,,,,,)iMsssss为不含噪的真实信号,12(,,,,,)iMnnnnn并且n为噪声,为噪声标准差,小波阈值降噪具体步骤如下:

(1)选择合适的小波函数对含噪观测信号x进行 L尺度小波分解,以获得相应的尺度系数12(,,,,,)llllliMAxAxAxAxAx与小波系数12(,,,,,)llllliMWxWxWxWxWx,1,2,,lL。

(2)由信号x,所对应的各尺度上小波系数lWx在某些特定位置上有较大值,这些点对应于原始信号s的奇变位置和重要信息,对于白噪声n,它所对应的各尺度上小波系数在lWx每一尺度上的分布是均匀的,并随着尺度的增加系数lWx的幅值有所减少,因此,消噪办法是寻找一个合适的小波系数的变化函数,把由噪声引起的小波系数设为0,而对特定位置的较大系数lWx保留@@,计算求得相应的小波系数 12ˆˆˆˆˆ(,,,,,)llllliMWxWxWxWxWx。

(3)利用尺度系数12(,,,,,)llllliMAxAxAxAxAx与经阈值函数变化后的小波系数12ˆˆˆˆˆ(,,,,,)llllliMWxWxWxWxWx重构信号,可得消噪后真实信号s的估计12(,,,,,)iMsssss。