宝典小波相干性分析.doc

综 述 小波相干分析及其应用 摘 要：将小波变换与相干分析相结合构成的小波相干分析，探测Fourier相干无法探测的特征信息，小波相干分析不仅能提供傅立叶分析类似的谱图，还能捕捉信号之间短时相互作用，因此小波相干分析在临床上的应用越来越广泛。本文主要介绍小波相干分析方法以及在生活中的应用。 关键词：小波分析；相干分析；小波相干；脑电信号；肌电信号 1 引言 随着科技的进步，信号处理在我们的生活中的作用越来越明显。在临床方面，脑电信号和肌电信号的分析，不仅有助于医师诊断病人的身体状况，而且还可以帮助医师进行康复工作。但因为生理信号是一种非常复杂的信号，信号本身非常微弱，稳定性较差，随机性很强，因而传统的Fourier相干在分析这些信号时存在一定的局限性[1-2]。小波分析方法对非平稳信号的特殊处理能力，使其在脑电和肌电信号的分析和处理中显示出极大的优越性。因此与相干分析相结合构成小波相干分析，既能够获取待分析信号的幅值和相位信息，又能够衡量相干性随时间的变化规律[3-4] 。 2 相干分析 对于两个复随机信号x和y，相干性系数定义为功率谱密度(power-spectrum density，PSD) 和互谱密度(cross-spectrum density，CSD ) 的函数，计算公式如下: (1) 公式(1) 中,Pxx(f)和Pyy(f)分别表示信号x和信号y的PSD,Pxy(f)表示信号x和y之间的CSD，PSD是频率f的实函数，而CSD是f的复函数。Cohxy表示信号x和信号y在频率f处的相干性系数，式中0≤Cohxy≤1，且Cohxy=0，x和y不相干;Cohxy=1，x和y完全相干。 相干性系数反映的是两信号之间的同步性相似性，或两信号的变化规律是否具有线性关系，该理论在地球物理雷达通信等方面都有着重要的应用，近年来也越来越多地应用于医学信号，如EEG和EMG。 当公式( 1) 中的信号x和y 分别为EEG中两个通道信号时，即可实现EEG信号的相干性分析，按照经典的频谱分析方法，设计步骤如下： (1) 对记录到的EEG时域信号进行傅立叶变换( FFT) ,得到F(x)和F(y) ； (2) F(x)与F(y)的乘积作为CSD,F(x) F(y)分别与其共轭相乘作为PSD； (3) 用CSD除以两信号的PSD进行归一化处理，如公式(1) 所示，得到Cohxy; 归一化处理是为了使相干性系数与两信号的震荡幅度相独立，从而保证相干性分析在动态功率谱变化中的有效性。 3 小波变换 一般认为，实际信号中不同频率成分的分量具有不同的时变特性，通常，慢变信号具有较低频率成分的频谱，变化激烈的信号具有较高频率成分的频谱。小波变换是由法国科学家莫莱特（Morlet）在1980年分析地震信号时提出的，在小波变换中，Morlet引入了多尺度分析的概念，可以由粗及细地逐步观察信号。 小波母函数的定义如下： （2） 其中，a, b∈R, a≠0, 分别是尺度参数和时间参数，母函数可以为实函数或者复函数。小波变换的实质是将信号与一个在时域和频域上均具有局域化性质的平移伸缩小波权函数进行卷积，从而将信号分解成位于不同时间和频率上的各个成分。小波变换的定义如下： （3） 式中*表示复共轭，(S,ψab)即表示小波系数[5]。小波变换克服了短时傅立叶变换窗函数固定尺度的缺陷，高频处采用短时窗以提高时间分辨率，低频处采用长时窗以提高频率分辨率。在小波分解中，随着分解尺度的增加，小波逐渐向低频方向聚焦。Morlet小波由于具有良好的时域与频域局部化特性，因此在信号的时频分析应用中经常被采用。复Morlet 小波是高斯窗口的复正弦函数，其表达式如下： （4） 复Morlet小波变换的定义如下： （5） 小波分析不仅具有完美的数学内涵，而且具有重要的应用价值，特别适用于信号的瞬态分析、图像边沿检测、图像去噪处理、模式识别、数据压缩、分形信号分析等方面。Akey[6]、Blinowska[7]、Meste[8]等人利用小波相干的时间—尺度能量分布分析方法对心血管音、脑电、晚电位等生物医学信号进行分析，所得尺度谱的分辨率比一般谱图的分辨率要高。在晚电位分析方面，目前仪器中采用的手段多是累加平均，人们希望能发展逐拍的动态检测，但是由于噪声(主要是肌电)干扰，小幅度的肌电与之很难区分。为了研究小波分析对心室晚电位动态分析的有效性，Tuteur[9]