Discussion :

[braniac]

Bonjour,
Je suis un novice de MatLab et j'aimerai avoir un coup de pouce concernant un programme de traitement de signal que j'essaie de réaliser.

En fait, je dois implémenter un programme qui permet de récupérer des valeurs de pressions échantillonnées à 1kHz à partir d'un fichier .txt et les filtrer pour éliminer le bruit qui vient se superposer sur le signal d'origine. Ensuite je dois les ploter en fonction du temps et de la fréquence.

Pour l'importation du fichier tout va bien j'ai réussi à importer toutes les valeurs et à les indexer
J'ai aussi créé une base temporelle de 0.001*30099, avec N=30099 le nombre d'échantillons de mon signal. et une base fréquentielle. Après je fais une fft de mon signal et c'est à partir de là ou commence mon problème car une fois je fais la fft je me retrouve avec des pics répétés et je vois pas comment je pourrais procéder pour éliminer ceux qui représentent les harmoniques.

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

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
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
Pression=smooth(Rampe(:,3));
 
Fs=1000;                                    %Fréquence de sampling
N=30099;                                    %Nombre d'échantillons
ind=0.001;                                  %Pas temporel de 0.001s
 
 
t(:,1)=0:(1/Fs):(N-1)/Fs;                   %Base temporelle
f(:,1)=(0:N-1)*Fs/N;                        %Changement de la base fréquentielle
 
X=fft(Pression,30099);                      %Serie de Fourier discrete
 
 
 
ty=abs(X);
 
 
 
%Changement des indices temporelles - incrémentation de 0.001 à chaque pas
for i=1:1:N
 
        Rampe(i,1)=ind;                     %Introduction du Pas temporel
        ind=ind+0.001;                      %Incrémentation
 
        if(i==N)                            %Sortir de la boucle lorsqu'on atteint le pas temporel maximal 
            break; 
        end;
end;
 
 
 
figure(1)
plot(t,Pression);
title('P2 - Temporelle');
 
figure(2)
semilogx(f(1:length(f)/2),ty(1:length(f)/2));
title('FFT de P2 sans filtrage - Frequentielle');
 
 
%Elimination du bruit et des harmoniques
 for i=1:1:30099
     if(ty(i)<6.765e5)          
          ty(i)=0;
     end;
 end;
 
 %Recomposition du signal temporel à partir du spectre du sgnal
 Ry=ifft(ty);

Je vous remercie d'avance.
Cdlt.

[Bébar]

Salut,

j'aurais plusieurs remarques à faire sur ton code :

*MATLAB te propose un outils assez puissant, la manipulation de matrices. La plupart du temps, tu n'as pas besoin de boucles et les oppérations se font dans des matrices. Je ne saurais trop dire si tu gagnes en temps de calcul (je crois), mais en tout cas tu gagne en lisibilité et ça fait plus classe (paie l'argumentation). Peut importe, ça dit, j'avoue utiliser des boucles quand je ne sais pas comment faire d'autre.
Par exemple toute ta boucle peut être remplacée par

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
 
pas = 0.001;
ind = 0:pas:(N-1)*pas %le tableau "pas" contient N élements linéairement répartis de 0 à (N-1)*pas, équivalent à linespace(,,).
%ou
ind = pas:pas:N*pas %pour ne pas commencer à 0;

*la fft, pour des raisons obscure (faut regarder comment est codée la fft) te fera gagner en vitesse si tu la fait sur une puissance de 2. Par exemple, pour toi, ça sera mieux sur 32768 (=2^15) points que sur 30099 points. Du coup, j'entend par là qu'il faut rajouter des points à ton signal... or tu ne peux pas en inventer d'autres. Mais tu peux rajouter des 0 (c'est du "zero-padding"). En fait c'est même bien de le faire :
-si deux signaux ont des fréquences proches, tu risque de ne pas bien les différencier si le nombre d'échantillons est faible (ici deux sinus de fréquences 1000 et 1100Hz avec 50 échantillons

(http://fr.mathworks.com/help/signal/...o-padding.html)

-si deux signaux ont des fréquences éloignés, tu peux avoir une fausse impression sur leur amplitudes (ici deux sinus d'amplitude '1' )




pour faire du zero-padding, (avec une fréquence de 2) tu peux entrer le code suivant :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
 
n_pad = nextpow2(length(Pression))+2; %tu utilise la puissance de 2 suivante (nextpow2) + 2 par exemple pour un encore plus important zero-padding.
%tu peux mettre juste nextpow2(...) ou rajouter +1, +3... à toi de voir
X=fft(Pression,2^n_pad); %cela va completer Pression avec plein de '0', pour avoir un signal de longueur 2^n_pad

*un autre point à voir est le fenêtrage. Comme il est tard, je laisse wikipedia expliquer : https://fr.wikipedia.org/wiki/Fen%C3%AAtrage . Si tu as une question, un truc pas compris : demande
à noter que le fenêtrage se fait donc avant la fft. De plus tu peux essayer deux fenêtre "blackman" et "hanning", matlab a deux fonctions du même nom pour les créer.

Une fois tout ça fait tu devrais pouvoir mieux observer ton signal.

Si tu peux envoyer ton fichier txt sur le forum, on pourra mieux y regarder

bien à toi,