Plotting spectrogram for EEG data

조회 수: 27 (최근 30일)
Sparsha Kumari
Sparsha Kumari 2020년 2월 5일
댓글: Walter Roberson 2021년 2월 7일
Greetings!
I am done with the filtering and plotting the power spectrum for the EEG data. Now I want to plot the spectrogram for this data. I have tried using the 'spectrogram' and 'stft' function and played with the various input options but I do not get the plot the way I want it to be. Could someone please show me some light in this regard.
Thanks!
P.S. : I need to get it as shown in figure 1 (Reference: William E. DeCoteau et. al. www.pnas.orgcgidoi10.1073pnas.0700818104). Instead I get it as shown in the figure 2.
  댓글 수: 3
Sparsha Kumari
Sparsha Kumari 2020년 2월 5일
Thanks a lot for your reply! I get it as shown in figure 2. I have used the 'stft' function as follows and got the plot,
stft(data, fs);
Sudon Biswas
Sudon Biswas 2021년 2월 6일
Could you plz provide me spectrogram and magnitude squared coherence (MSC) Matlab code for ECG signal? The code must be for single signal example: 1. For raw signal 2. Noisy Signal 3. For noise free signal

댓글을 달려면 로그인하십시오.

답변 (1개)

Daniel M
Daniel M 2020년 4월 24일
You should plot the absolute value of your output from spectrogram or stft. Maybe even the log, or 10*log10 too (but I don't recall exactly, you'll have to look up the documentation). Basically, when you call those functions with outputs, try and make the plot match the one you get when you call the functions without outputs.
Your spectral and temporal resolution are very different from the first figure. That paper is doing a time-frequency analysis over ~ 4 seconds and 25 Hz. You are looking at ~150 ms and 500 Hz. You can look at this tutorial (with video) on how time window and frequency resolution trade-off with each other.
  댓글 수: 5
Sudon Biswas
Sudon Biswas 2021년 2월 6일
@Walter Roberson
Walter Roberson
Walter Roberson 2021년 2월 7일
Coherence is only defined between two signals, not one signal by itself. https://en.wikipedia.org/wiki/Coherence_(signal_processing)

댓글을 달려면 로그인하십시오.

카테고리

Help CenterFile Exchange에서 EEG/MEG/ECoG에 대해 자세히 알아보기

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by