traitement du signal EEG

Bonjour à tous,

je suis en train de réaliser un projet en traitement du signal sous Matlab et j'ai rencontré quelques difficultés

en fait j'ai des données sur excel qui ont des enregistrements du signal EEG c'est une matrice d'une seule colonne

j'ai pu récupérer ce fichier et tracer le signal correspondant avec la fréquence d'échantillonnage qui est déja données par l'utilisation de ce programme

Code:

---

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15

close all ; clear all; clc [fileout,pathout] =uigetfile('*.CSV'); chemin = [pathout,fileout] Valeurs=xlsread(chemin); A=size(Valeurs)-1; Fs=200; ts=1/Fs; t1=0:ts:(A/Fs); t1=t1';   plot(t1,Valeurs,'r') grid on title('signal EEG');

---

mais quand je veux faire la transformée de fourrier de signal il marche pas malgré que cette fonction fonctionne avec les signaux comme le sinus, le cosinus ..
et voila le programme

Code:

---

1 2 3 4 5 6 7 8 9

% Transformé du fourier du signal EEG N=length(t1); f=Fs*(0:N-1)./N; f=f(:); X= fftshift(fft(Valeurs)*ts/2);   subplot(211), plot(t1,Valeurs), title('x(t)') subplot(212), plot(f,abs(X)), title('module(X)') grid on

---

Merci pour votre aide

Citation:

---

il marche pas

---

C'est à dire ? bug (lequel) ? resultats bizarres (en quoi) ?

normalement y'a pas d'erreurs mais la transformé de fourrier ne donne pas le résultat exacte

il reste toujours le problème au niveau d'amplitude quand je fais la transformé de fourrier il me donne des amplitudes très grande par rapport au signal représenté dans le domaine temporelle :(

Essaye ce code mounya

Code:

---

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27

close all ; clear all; clc [fileout,pathout] =uigetfile('*.CSV'); chemin = [pathout,fileout] Valeurs=xlsread(chemin); x=size(Valeurs)-1; %your signal   Fs = 512; % recall Nyquist sampling frequency, it must be % at least twice that of sampled signal % but it is better if it is higher than that   L = length (x); % length of signal (number of sample values)   T = 1./Fs; % sampling period t = [0:L-1]*T; % time scale for time domain f = Fs/2*linspace(-1,1,L); % frequency scale for freq. domain   signal=x; % generated signal is anything we like FFTsignal = fft(signal)/L; % fft of signal (freq domain) figure(); % always create new figure for your plots subplot(2,1,1),plot(t,signal); %axis([0 3/f1 -1.1 1.1]); title('Signal in time domain'); subplot(2,1,2),plot(f,abs(fftshift(FFTsignal))); %axis([-2*f1 2*f1 -1.1 1.1]); title('Signal in frequency domain');

---

:roll:

Bonjour,

SVP, si on est dans un cas inverse; on cherche à restituer le signal à partir de ces paramètres extraites de la transformée de Fourier.
Ici je ne veux pas appliquer la FFT-1.
Y a-t-il un algorithme qui mais la somme des sin et des cos sachant que la série de fourier est une somme des sinus et cosinus, tout en sachant les paramètres?

mercii pour votre réponse, cordialement

Vous voulez faire une transformée de fourier inverse, mais sans la transformée de fourier inverse?

pour restituer le signal sans fft inverse vous pouvez utiliser les ondelettes 8O