analyse spectrale ...

S

sweetsalt

Guest
Salut à tous!
Je suis invité à effectuer une analyse spectrale en utilisant des techniques de Fourier de la 64 suivante points.which sont obtenues par échantillonnage d'un signal continu à 4000Hz.

-2.2084441e-001
1.0841281e-001
-2.2566847e-002
-7.1269979e-001
-7.7570468e-001
1.5490627e 000
1.9458210e-001
-8.9878564e-001
2.5642118e-001
1.2323886e 000
-2.5090429e-001
-6.4754173e-001
-2.8157727e-001
1.8196999e 000
-3.4228194e-001
1.0446433e-000
5.8907411e-001
4.8584173e-001
7.0878880e-001
2.1880127e-001
1.2519184e 000
2.9170207e-001
1.2549903e 000
1.2138000e 000
-3.1823385e-001
-4.6785664e-002
-2.3964629e-001
1.7661788e-000
1.4964292e-000
-3.6095019e-001
1.1122538e-000
2.1902688e-002
6.1575414e-001
1.3917852e 000
-3.2921212e-001
7.7902325e-001
6.3796904e-001
1.2476316e-000
9.7494022e-001
-7.5831465e-001
1.0848348e-000
6.0460012e-001
1.4021865e-000
1.3785847e 000
-6.5360020e-001
-4.9165649e-001
9.8713534e-001
1.0617638e-000
1.0888472e-000
6.9453009e-001
6.1975115e-001
8.5906752e-001
9.9296393e-002
1.8264694e-001
-6.5977471e-001
1.0669556e-000
1.5649163e-001
3.4083024e-001
8.9142667e-001
6.5523161e-001
1.1048760e-001
-8.1966173e-001
-2.1416913e-001
1.6509091e-000

Maintenant, le problème est quand je fais, de la tige (abs (fft (x)), ce que je dois faire doit être symétrique autour de son centre, car elle un signal réel.
signal réel a seulement la fréquence de 0 à pi (les commentaires nécessaires??)
si les valeurs de 0 à pi doit répéter de pi à 2 pi,

Mais dans mon cas FFT n'est pas symétrique, la première valeur est à l'origine problème, comme j'ai 64 FFT point, 31 valeurs autour de valeurs sont symétriques 33e, mais la valeur à 1 est comme nulle part dans le signal de plus, y est tellement thanksggiving? ??

si quelqu'un peut faire l'analyse spectrale du signal donné pour moi, ce serait une grande aide.

je fais

YF = fft (x, 64);
y = yf (1:33);% il suffit de regarder la première moitié
m = abs (y);% m = ampleur des sinusoïdes
p = déballer (angle (y)); p% = phase de sinusoïdes, déballer () Copes% avec un degré sauts 360
w = 0: pi/32: pi;subplot (2,1,1), la tige (W / pi, m); grandeurs complot%
ylabel ('Abs. Magnitude »), grille;
xlabel ('omega \ / \ pi');
subplot (2,1,2), la tige (W / pi, p * 180/pi); phase complot% en degrés
ylabel («Phase [degrés]»), grille;
xlabel ('omega \ / \ pi');plz aidez-moi d'urgence, le délai est 2morow.
Merci d'avance

 
Salut
Il n'y a pas proplem dans le u Obtenez des résultats ...il est correct et semetrical.
FFT NE PAS RETOURNER fréquence de 0 à 2 * PI ..Mais à partir de 0 à (63/64) * 2 * PI à l'étape 2 * pi/64.
La première composante est la DC et si vous ne trouverez nulle part ailleurs dans le signal.
Le spectre du signal échantillonné est des répliques de la continu, si bien
la composante à 33 est la valeur à pi.
la composante à 64 est la valeur à 63/64 2Pi
la valeur à 2 pi est connue et égale à celle à 33

Encore une fois u ne faut pas oublier que le spectre du signal échantillonné n'est pas de 0 à pi ...il est entended de-inf à inf en répétant spectre de-pi à pi (qui est semetrical)

 

Welcome to EDABoard.com

Sponsor

Back
Top