基于小波变换的频谱分析方法.docxVIP

基于小波变换的频谱分析方法 一、 采样频率选取 频带分析通常采用快速傅立叶变换（fft）方法，但高速傅立叶变换信号频率的计算精度由采样频率和采样点决定。假设采样频率为fs, 采样点数为N, 则频率分辨率为Δf=fs/N。对于待检测信号如果采样频率选取过高, 将产生过饱和现象;如果采样频率选取过低, 将产生欠采样现象。这两种情况就使得采样频率的选取受到限制。实际应用中, 采用增加采样点数来提高频率分辨率, 其缺点在于增加运算量, 使实时性降低;此外, 有些实际信号的变化比较快, 而增加采样点数必然增加采样时间, 这样就导致频谱检测跟踪不上信号的变化, 从而限制了使用增加采样点数来提高频率分辨率的应用场合。 本文提出了基于小波变换的频谱分析方法, 通过小波变换方法来提高频率分辨率, 对短数据仍然具有较高的频率分辨率, 从而满足信号频谱检测准确度需求。 二、 小波变换方法用于提高光谱分辨率 1. 基于尺度的频谱检测 小波变换的实质是使信号通过一组带通滤波器, 这些滤波器的通带特性决定着信号频谱检测的精度。对于不同的信号希望有一个适合该信号的带通滤波器。根据小波变换的理论, 小波滤波器的中心频率和通带带宽与尺度α成反比, 这样就可以通过选择变换尺度α的值来选择滤波器的特性完成最佳的频谱检测。从而实现较短数据下的高分辨率的频谱检测, 使用Morlet小波, 其基本小波函数为 Ψ(t) =e-t/Teuω0t(1) 其对应不同尺度的频率特性为 Ψ(αω)=π/T????√?e[?(αω?ω0)2/(4/T)](2)Ψ(αω)=π/Τ?e[-(αω-ω0)2/(4/Τ)](2) 式 (1) 和式 (2) 中,ω0是最大谱估计频率, 为采样频率fs的一半;α是尺度因子;T为采样时间间隔。 由式 (1) 和式 (2) 可以得出, 滤波器的中心频率为ω0/α=fs/2α, 其带宽随中心频率变化而变化。中心频率高, 带宽宽;中心频率低, 带宽窄。但带宽和中心频率的比值不变, 这是常数Q滤波器。缓慢改变尺度就会使得滤波器的中心频率缓慢变化, 两个邻近的滤波器的通带有共同的部分, 信号通过这共同部分时, 都可以得到正确的频谱检测结果。通过改变尺度α的值就可以精确控制滤波器特性, 从而提高信号频谱检测的精度。实际中, 只需要在感兴趣的频率段改变尺度就可以了。 假设单一正弦信号通过连续的一组滤波器, 这些滤波器的通带覆盖着0～fs/2的频率段, 相邻的滤波器之间具有共同的通带频率段。尺度α的变化控制这些滤波器的特性, 如果尺度α=1, 则频率分辨率仅仅决定于采样点数和采样频率, 这就和经典的FFT有同样的频率分辨率;如果使用尺度α=0.5, 在同样采样点数下, 则频率分辨率提高一倍, 从而得到精确的频谱检测。使用尺度α为一个变化的量就可以实现不同频率段具有不同的频率分辨率。实际情况下, 我们通过改变尺度α设置滤波器的带宽为滤波器中心频率的十分之一, 每一个尺度对应一个滤波器, 采样后的信号每经过一个滤波器都会得到一个检测结果, 信号经过其中一个或两个相邻的滤波器后完全通过, 从而得到正确的频谱检测结果。大部分滤波器将滤掉信号, 只得到噪声。这通过比较数据可以判断出哪一组滤波器得到信号。对于信号幅度大于噪声幅度的信号检测, 通过简单的比较数据大小就可以得到正确的结果;对于信号幅度小于噪声幅度的信号检测, 只有通过多次检测才能判断出信号所在的频率段, 此种情况, 不在本论文考虑之中。 2. 调整采样频率 在信号频谱检测中, 总是希望信号频率正好落在带通滤波器通带的中心处。这样, 信号几乎完全通过, 而噪声达到最大限度的抑制。在实际信号频谱检测中, 仅仅改变尺度的值不能保证信号频率正好落在带通滤波器通带的中心处, 从而很难达到对检测精度要求很高的信号频谱检测, 文献中达到10‰左右的检测精度。为了达到1‰的检测准确度, 必须将信号尽可能的通过滤波器的中心处, 也就是信号频率尽可能的接近滤波器通带的中心频率。我们采取调整采样频率来确保信号通过滤波器通带的中心处。这样, 需要进行4～5次左右的采样频率的调整, 一般情况下三次就满足要求。调整采样频率的过程如下:首先, 选择一个适应范围较广的采样频率 (例如, 1kHz采样可以适用于100Hz～500Hz信号频率) 得到一次精度较差的检测结果 (为了提高检测时间, 可以将尺度设置为1, 和FFT有同样的频率分辨率, 也可以直接使用FFT做第一次的频率检测) ;再紧接着以上一次的估计结果调整采样频率;一般选择采样频率为上一次检测结果的16倍, 重复这个过程, 直到最后两次频谱估计结果一样或很接近, 完成频率检测, 得到准确度较高的检测结果。 这样调整采样频率的优点有三点, 第一, 保证了信号频谱检测的准确度;第二, 初