Block FFT Visualisation

Bonjour,

Je voudrais connaitre la FFT d'un signal sous SIMULINK.

Voici comment j'ai effectué mon schéma:

Signal=>block fft=>Scope

La sortie de fft étant un complex, je ne peux pas visualiser le spectre en fréquence de mon signal. Que dois-je faire?

Merci d'avance pour votre réponse

Bonjour.
Une solution :
http://img13.imageshack.us/img13/1668/fft0905.gif

merci pour cet élément de réponse.

Maintenant j'obtiens l'erreur suivante concernat le buffer: "Continous sample times not allowed"

Effectivement, mon signal d'entrée est un signal continue.

Il faudrait donc je le transforme en signal numérique.

Je me trompe?

Citation:

---

"Continous sample times not allowed"

---

Tu as mis une valeur non nulle dans Sample time aux box des signaux d'entrée (0.001 parexemple) ?

Citation:

---

Envoyé par phryte

Tu as mis une valeur non nulle dans Sample time aux box des signaux d'entrée (0.001 parexemple) ?

---

Je pense qu'il y a un malentendu. En réalité, je veux établir le spectre de 2 signaux.

Le premier étant l'entré d'une fonction de transfert de type 1/(1+to*s) (donc continu) et le second étant la sortie.

Comment obtenir les spectres ?

J'ai essayé en mettant:
Signal=>zero-order hold (avec un sample time différent de 0)=>FFT=>Spectrum Scope

Le spectrum scope est "vide" à la fin.

savez vous où est mon erreur?

Je ferais comme cela :
http://img14.imageshack.us/img14/704...uxspectres.gif

Sur un même graphe :
http://img4.imageshack.us/img4/540/fftdouble.gif

ok, merci beacoup pour votre réponse! Malheureusement, je n'ai pas les résultats espérés...
Je posterai une suite dès que j'en aurai le temps.
Merci encore !

Bonjour, en fait l'image schéma.jpg montre ce que je souhaite simuler.
Le signal d'entrée e est donné par:

Code:

---

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

clear; clc; pas=0.0001; tmax=30; freq=5; alpha=1; A=4; B=2; def=pi*182/180; t=0:pas:tmax;   y1= (A/2)*square((freq/2*pi)*t)+A/2; y2= (B/2)*square((freq/2*pi)*t+def)+B/2; y3=(1-exp(-alpha*t)); e=y1+y2.*y3;

---

Voici ce que j'utilise pour étudier le signal sous matlab.

Code:

---

1 2 3 4 5 6 7

plot(t,e) % % E = fft(e,512); % Pee = E.* conj(E) /512; % f = 100000000*(0:256)/512; % plot(f,Pee(1:257))

---

Que ce soit sous Matlab ou Simulink, je n'obtiens pas les résultats désirés. Explication: Lorsque je cherche la FFT d'un signal carré de fréquence, je devrais obtenir une courbe avec les fréquences: f, 3f, 5f,... mais je n'obtiens pas ça.

Lorsque je prends le signal crée sous Simulink et que je l'importe avec une boite to workspace, E (E = fft(e,512);) n'a qu'une composante non nulle. Pourquoi?

Y a t il un problème dans mon programme ou dans mon schéma SIMULINK?

Merci pour vos réponses

Bonjour.
Une solution avec Matlab :

Code:

---

1 2 3 4 5 6 7 8

NFFT=length(e); E = fft(e,NFFT); f = linspace(0,freq/2,NFFT/2); E_plot = abs(E(1:floor(NFFT/2))); figure(2) semilogx(f, 2*E_plot) xlabel('Frequency [Hz]'); ylabel('Amplitude');grid

---

Avec simulink (le gain K est <0) :
http://img196.imageshack.us/img196/2043/fft0916.gif
http://img8.imageshack.us/img8/5245/scop0916.gif

Pour la simulation, j'utilise les paramètres suivants:

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 28 29 30 31 32 33 34 35 36 37 38

Pulse Generator:   Pulse type: Time Based Time: Use simulation Time Amplitude: A Period: 1/freq Pulse Width: 50 Phase Delay: 0   Pulse Generator1:   Pulse type: Time Based Time: Use simulation Time Amplitude: B Period: 1/freq Pulse Width: 50 Phase Delay: def   Ramp:   Slope: 1 Start time: 0 Initial output: 0     Zero-Order Hold:   Sample time: 1   Magnitude FFT:   Output: Magnitude squared FFT length: 256   Vector Scope:   Input domain: Frequency

---

En ce qui concerne l'étude du signal carré suivant: square(t) de fréquence 1/pi, voici la FFT que j'obtiens (cf FFT.JPG).
Je ne retrouve pas les fréquences f, 3f, 5f,... Est ce normal?

Vous m'avez déjà beaucoup aidé. MERCI!

Deux corrections et c'est bon :
For example, generate a 30 Hz square wave:

Code:

---

1 2 3 4

t = 0:.0001:.0625; y = SQUARE(2*pi*30*t);, plot(t,y)   f = linspace(0,1/pas/2,NFFT/2);

---