Problème d'affichage d'une fft

[Alays]

bonjour je veux afficher 8000 échantillons d'un cosinus donc je fais un cos de fréquence 900hz donc la fréquence d’échantillonnage est >900*2 donc je prends large 3000 hz et je fais:

Code :

1
2
x=fft(cos(2*pi*900*(0:7999)/8000*3000));
plot(abs(x))

le problème c'est que je ne vois pas une barre correspondante à 900 donc je ne vois vraiment pas pourquoi

[Winjerome]

Bonsoir,

En remaniant un peu l'expression de ton signal, on a :

Code :

1
2
cos(2*pi*(900*3000)*(0:7999)/8000)
cos(2*pi*2700000*(0:7999)/8000)

Autrement-dit un signal de fréquence 2700000 Hz

Commence par générer ton vecteur temps à part.

[Alays]

je l'ai fais à part et c'est pareil:
pourtant mon vecteur temps c'est bien ça fs*(0:n-1)/n

Code :

1
2
t=3000*(0:7999)/8000
cos(2*pi*900*t)

[Winjerome]

Heu, je voulais dire, en réfléchissant en même temps sur ce que je t'ai dit
Ta formule est fausse. Fais une analyse des unités.

[Alays]

ben je dirais t=0:1/3000:7999/3000; 1 échantillon toutes 1/3000 s autrement je ne vois vraiment pas

[Alays]

alors c'est bien ça ? je ne suis vraiment pas sur

[tachmou]

Salut

Code :

1
2
3
4
5
6
7
8
9
10
11
12
t_initial = 0;
t_final = 2; %sec
dt = 1/3000;
frequence_echantillonage = 1/dt;
t = t_initial:dt:t_final;
f = 900; %Hz
w = 2*pi*f; %rad/sec
signal = cos(w*t);
x=fft(signal);
Pxx = x.*conj(x); %puissance
f = frequence_echantillonage*(0:(size(x,2)-1))/size(x,2);
plot(f,Pxx);

Pas testé, mais il vaut mieux faire les étapes pas à pas. C'est plus clair et ça va plus vite.

Cordialement,

[Winjerome]

@tachmou : Ce code est faux, la mise ne correspondance fréquences/amplitude est fausse.

@Alays : As-tu retracé ta fft depuis ? Tu en conclus quoi ?

Le sujet sur la fft a été abordé d’innombrables fois, surtout récemment. Une recherche ne sera pas de trop.

[tachmou]

Salut

Euh, qu'est-ce que tu appelles la mise en correspondance amplitude/fréquence ? Je viens de tester mon code, et j'obtiens bien les fréquences qui ressortent...

Cordialement,