How to extract mean frequency from continuous wavelet transform (CWT)

조회 수: 3 (최근 30일)
Brent
Brent 2017년 3월 22일
댓글: Christoph Mohl 2023년 11월 27일
Hello,
I have used the cwt function in Matlab 2017a to create a scalogram of an electromyography (EMG) signal. My goal is to compare the time varying mean frequencies of different EMG signals. The issue is that I'm not sure how to use the complex double output (wt) to calculate a usable mean frequency 1D waveform.

채택된 답변

Santhana Raj
Santhana Raj 2017년 3월 22일
Hi,
If your question is to how to get useful information from the complex matrix of 'wt', then use the command 'abs' to get the absolute magnitude of the complex variable.
If your question is to how to extract mean data from a vector, then compute the average of the signal. You can do the same for different scales in your wt or overall average value across all scales.
Hope this helps.
Raj
  댓글 수: 4
Brent
Brent 2017년 11월 27일
Hi Jungyeon. Below is the code I used to calculate IMNF. It's based on the paper above and uses the [wt, f, coi] outputs from the cwt function. you'll likely need to filter the IMNF as the outcome is quite noisy.
wta=abs(wt);
s_max=length(f);
s=1:s_max;
w0 = f(1);
for i=1:size(wta,2)
w=w0./(s(2:end));
IMNFi(i)=trapz(w',w'.*wta(2:end,i))./trapz(w',wta(2:end,i));
end
Christoph Mohl
Christoph Mohl 2023년 11월 27일
Hi Brent,
even though it´s been a while since you posted the code above, I hope you can still help me with one question.
Why do you divide w0 by the second to the last scale?
Kind regards,
Christoph
w=w0./(s(2:end));

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

추가 답변 (0개)

카테고리

Help CenterFile Exchange에서 Time-Frequency Analysis에 대해 자세히 알아보기

Community Treasure Hunt

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

Start Hunting!

Translated by