2010-04-21 47 views
3

我正在尝试为我的应用程序构建一个鳞片图视图,以查看是否存在我们可以从小波变换中检索到的相关信息,而不是使用谱图来查看可以通过FFT检索的内容。生成信号的鳞片图

到目前为止,我可以采取一个波形,我可以对其执行正向小波变换。但是我在下一步中迷失了。我如何将这些信息转化为电力/能源信息?我有一组不同频率的波形,但我有,如我所说,没有频率信息。

任何人都可以告诉我下一步是什么把这个转换后的数据转换成一个鳞片?

任何帮助将不胜感激,因为我的谷歌技能让我失望!

回答

4

在合理的假设下,离散小波变换(DWT)将时间序列的功率/能量/方差分解为尺度。这是一种能量保留变换,原始时间序列中包含的总方差包含在平方小波系数(正确归一化)中,就像DFT一样!我认为Percival和Walden的文字Wavelet Methods for Time Series Analysis是这类信息的极好资源。现在,连续小波变换(CWT)是一种冗余变换,并且在DWT中的能量保持特性不再成立。但是,您仍然可以查看平方小波系数以生成与短时傅里叶变换(STFT;又名谱图)类似的“尺度图”。

Torrence和Compo有一些不错的wavelet software来做到这一点,用一篇论文来解释大气时间序列背景下的方法。我相信matlab小波工具箱也可以执行CWT和相关的尺度图。

+0

谢谢。我已经发现了能量保存的事情。所以我通过计算每个频带的功率来显示一些小波数据。它不是很好,但肯定显示了我一些信息:)虽然......我仍然失去了如何做CWT ......花了相当多的时间来研究它。我认为它不仅仅是做一个滑动窗口DWT? – Goz 2010-04-22 12:07:34

+0

两个(主要)事物将CWT与DWT区分开来:(1)小波和(2)用于计算小波系数的平移。 CWT使用连续函数是小波生成函数(Morlet,高斯的一阶导数等),并计算所有尺度下所有翻译的小波系数。因此,1D信号(时间序列)的CWT产生(严重相关的)小波系数的2D图像。 DWT使用离散函数作为小波GF,并产生与原始观测相同数量的小波系数。使用Torrence/Compo的软件! – 2010-04-23 10:06:46

2

我从来没有做过鳞状图,所以我不会声称尝试知道任何东西。然而,我已经找到了matlab的尺度函数的代码,这个函数的评论方式应该可以帮助你更好地理解事情....我希望:)。

http://www.csse.uwa.edu.au/~pk/Research/MatlabFns/FrequencyFilt/scalogram.m

我假设你使用MATLAB为您的DSP的东西,并知道如何解释是m文件。可能是所有dsp的东西的一个强制性程序,我猜。