FFT of EMG Signal

조회 수: 35 (최근 30일)
Karl Rueth
Karl Rueth 2017년 10월 22일
댓글: ANADI BISWAS 2021년 12월 28일
I am trying to obtain the FFT of an EMG Signal, I have already done some other type of processing, including; removing DC Offset, Rectification and Signal Envelope, however I am stuck in the FFT part ? Any help would be appreciated ..
Attached the Code and Photo of current situation can be found.
  댓글 수: 3
Karl Rueth
Karl Rueth 2017년 10월 23일
For some reason, the final result, isn't the required signal.
dpb
dpb 2017년 10월 23일
What "final result", there's no FFT uncommented in the sample code; the result will be what it shall be; you can't really "require" anything.
You don't even say what your desired result is that you're looking for, specifically, what are we, mind-readers?
Look at the example at
doc fft
for example of doing PSD via FFT for pointers.

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

답변 (3개)

Star Strider
Star Strider 2017년 10월 22일
See the documentation on fft (link).

Luc Dransfeld
Luc Dransfeld 2018년 1월 8일
If you follow the FFT documentation from Matlab, you should be able to get there. Maybe though, you should set your plot window differently, as motion artifacts can contain a power that is proportionately way bigger that the rest of the EMG-frequencies resulting that you will only visibly see the motion artifacts. You can always use matlab function: periodogram, which basically gets you the logarithm from the power, expressed in dB.

Trishna Sah
Trishna Sah 2019년 1월 23일
편집: Trishna Sah 2019년 1월 23일
This is using your variables in the fft doc by Matlab
%fft plot for emg4
Fs=1000;
L=length(emg4);
t=L/Fs;
x=emg4;
Y = fft(x);
P2 = abs(Y/L);
P1 = P2(1:L/2+1);
P1(2:end-1) = 2*P1(2:end-1);
f = Fs*(0:(L/2))/L;
plot(f,P1)
title('Frequency Spectrum Plot')
xlabel('f (Hz)')
ylabel('|P1(f)|')
f = Fs*(0:(L/2))/L;
plot(f,P1)
title('Single-Sided Amplitude Spectrum of X(t)')
xlabel('f (Hz)')
ylabel('|P1(f)|')
  댓글 수: 5
dpb
dpb 2020년 11월 15일
Ayup.
ANADI BISWAS
ANADI BISWAS 2021년 12월 28일
Units of x axis and y axis after abs(fft)?

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

카테고리

Help CenterFile Exchange에서 Spectral Measurements에 대해 자세히 알아보기

Community Treasure Hunt

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

Start Hunting!

Translated by