Cześć
piszę do Was ze swoim jakże osobistym problemem, który spędza mi sen z powiek od dobrych dwóch nocy.
Otóż historia wygląda tak: projekt polega na analizie sygnału EMG pod kątem zmęczenia mięśni. Analizy dokonuje się na podstawie wykresu częstotliwościowego więc to co mam trzeba przemielić Fourierem i wypluć jakieś sensowne dane.
Na stronie Link
Ładnie opisane jest właściwie wszystko co mam zrobić. Jedynie nie potrafię wygenerować odpowiedniego Fouriera, który powinien wyglądać mniej więcej tak jak poniżej.
Czyli po lewej surowy zapis EMG, a po prawej już przemielony przez fft.
W Matlabie wygląda to u mnie mniej więcej tak:
Wyniki wyglądają mniej więcej tak:
-surowe dane
-Fourierek
Napiszcie proszę dobrzy ludzie co robię źle...
Będę wdzięczny za pomoc
Pozdrawiam,
Tomek

piszę do Was ze swoim jakże osobistym problemem, który spędza mi sen z powiek od dobrych dwóch nocy.
Otóż historia wygląda tak: projekt polega na analizie sygnału EMG pod kątem zmęczenia mięśni. Analizy dokonuje się na podstawie wykresu częstotliwościowego więc to co mam trzeba przemielić Fourierem i wypluć jakieś sensowne dane.
Na stronie Link
Ładnie opisane jest właściwie wszystko co mam zrobić. Jedynie nie potrafię wygenerować odpowiedniego Fouriera, który powinien wyglądać mniej więcej tak jak poniżej.

Czyli po lewej surowy zapis EMG, a po prawej już przemielony przez fft.
W Matlabie wygląda to u mnie mniej więcej tak:
Code:
clear all;
close all;
clc
%Obróbka po imporcie danych
sampling = 0.0001; %probkowanie
freq= 10000; % czestotliwosc probkowania
dane=importdata('pomiar5.dat');
[dlugosc nic] = size(dane.data); %dlugosc wektora danych
N=dlugosc;
%wyznaczanie następnej potęgi dwójki do długości Fouriera
NFFT=2^nextpow2(N);
czas=0:sampling:dlugosc/freq-sampling;
%import danych
dane1=dane.data(:,2);
%Fourier właściwy
Y=fft(dane1,NFFT)/N;
f=freq/2*linspace(0,1,NFFT/2+1);
figure(1);
plot(czas,dane1);
figure(3);
plot(f,2*abs(Y(1:NFFT/2+1)));
Wyniki wyglądają mniej więcej tak:
-surowe dane

-Fourierek

Napiszcie proszę dobrzy ludzie co robię źle...
Będę wdzięczny za pomoc

Pozdrawiam,
Tomek