2013-10-14 40 views
0

我試圖重新以下結果:譜圖(重建從數據條結果)

Spectrogram

從以下數據:

https://www.dropbox.com/s/mi3szqgzgku29rn/FS40.dat

的時間以毫秒爲單位(頻率爲40000Hz)

文章狀態說明他們用的是複合Morlet wav ELET創建譜圖:通過與所述形式

瓦特(t的復Morlet小波卷積的信號從 從1-88赫茲的數據的時間 - 頻率頻譜圖計算 「從平均電力個LFP估計,對於每個感興趣的頻率f 0,其中σ= m /2πf0,並且i是(f 0)= A exp(2 * i π f 0 * t)

虛構單位。標準化因子是A = 1 /(σ(2π)^ 0.5),並且定義時間和頻率分辨率之間的折中的常數m是7.我只使用頻譜圖函數設法得到一些「好」結果在matlab中。 但是我沒有太多的關於如何使用morlet復小波的想法。 我試圖用'morl'窗口

嘗試使用cwt時得到了不好的結果謝謝。

P.S.

我試圖重新創建本文:通過低閾值尖峯和快速尖峯的interneurons控制可分離諧振電路:

由細胞類型選擇性光遺傳學驅動器驅動不同新皮質振盪的計算建模。

+1

嗯,我想說你被文中的措詞誤導了。圖中顯示的是信號能量的時間 - 頻率圖,但它不是頻譜圖。後一項更爲普遍,每張譜圖都是信號能量的時間 - 頻率圖,但反之亦然。 Matlab中的'spectrogram'函數特別使用短時傅立葉變換而不是小波變換。其次,如果想要估計頻譜能量,使用復小波很奇怪。 –

+1

對於時間頻率分析,特別是Matlab中的神經生理數據,我推薦[FieldTrip](http://fieldtrip.fcdonders.nl/)。 'ft_freqanalysis'包含一個使用Morlet小波的選項。 –

回答