2010-04-21 49 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的東西的一個強制性程序,我猜。