Discussion :

[azassma]

Bonsoir,

Je calcule l'amplitude spectrale d'un signal avec LabVIEW en utilisant la FFT(r,teta) et quant j'utilise la fonction FFT de matlab je n'obtient pas le bon résultat?

Code :

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

1
2
3
4
5
6
7
8
9
10
11
 
fid = fopen('fichier.txt', 'r');
x1 = fscanf(fid, '%f ')  ; 
fclose(fid);
 
n=length(x1);  Te=1/100000;  tmax=n*Te;  tmin=0;  t=tmin:Te:tmax-Te;
 
fe=1/Te;  fmax=fe*n;  fmin=0;  df=(fmax-fmin)/n;  f=fmin:df:fmax-df;
 
y=(fft(x1));
figure,plot(f,y),grid,title('spectre du signal')

S'il vous plait est ce que vous pouvez m'aider ? Merci beaucoup

[Gooby]

Bonjour,

Je pense que tu dois regarder le module de la FFT pour avoir ton spectre. Dans tout les cas, je te conseille de faire un tour sur ce sujet qui pourra peut être t'aider dans ton travail sur le spectre d'un signal.

[azassma]

Merci pour votre réponse
J'ai utilisé le module mais ça ne donne pas le bon résultat comme dans LabVIEW
Voila le code que j'ai utilisé

Code :

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

1
2
3
4
5
6
7
8
 
 
Y2=fft(y1);
 
for i=1:n
    amplitudeSpectrale(i) = (imag(Y2(i)))^2+(real(Y2(i)))^2;
    amplitudeSpectrale(i) = sqrt(amplitudeSpectrale(i));
end

[Gooby]

Pour le module d'un nombre complexe tu peux te simplifier la vie en utilisant abs(monNombre)

Tu peux également jeter un oeil sur la fonction proposée par WinJerome pour le Calcul et tracé de la FFT, peut être y trouveras-tu ce que tu désires.

PS: peux tu nous donner le résultat que tu obtiens?

[azassma]

Merci Gooby pour votre réponse, j'ai utilisé la fonction abs mais ça ne donne pas le bon résultat

voila le résultat obtenu à l'aide de LabVIEW, sachant que l'axe des abscisses est une fréquence



et voici celui obtenu avec MATLAB

[Gooby]

Premier point: c'est étrange que l'axe des abscisses de ton grahique LabView soit 'Time' !

Second point, l'axe de symétrie du graphique Matlab se trouve au milieu de ton axe des abscisses, tu peux représenter uniquement la moitié de ton graph (ce qui représente une réalité physique puisque les fréquences négatives n'existent pas)

Je pense que tu devrais essayer de zoomer sur ta zone basse fréquence et regarder un peu où se trouvent tes pics, il se peut que ton graphique ressemble plus à celui de LabView que tu ne le penses.

[azassma]

Merci mr Gooby, effectivement quand je fais un zoom j’obtiens une bonne représentation.

Pour le premier point l'axe des abscisses c'est la fréquence qui étais représentée, juste dans le LabVIEW je n'ai pas modifié le titre "time"
Pour le deuxième point les fréquences sont toutes positives donc je n'ai pas compris pourquoi existe-t-il de symétrie.

[jyber]

Cette symétrie (dite hermitienne) est une propriété de la transformée de Fourier pour les signaux réels. La seconde moitié du résultat est un miroir de la première. Attention avec cette symétrie! Le dernier nombre n'est pas le symérique du premier, mais du second. Le premier joue un rôle particulier (amplitude à la fréquence zéro, soit composante continue), ainsi que le nombre médian à l'indice N/2 (c'est l'amplitude à la fréquence de Nyquist). C'est une source d'erreur fréquente quand on manipule des FFT complètes.

Pour un signal réel, on peut se contenter de la "couper en deux".

[azassma]

Merci Jyber j'ai bien compris