Discussion :

[laurent12]

Bonjour a tous,

Je travaille sur un projet simulink-matlab, dans lequel je recupere le signal de sortie d'un de mes blocs simulink dans un GUI, pour tracer son spectre. Le probleme est que je n'arrive pas a avoir une bonnes echelles en frequence.

Lorsque je tente d'afficher le spectre d'un sinus d'amplitude 1 et de frequence 900kHz, mon spectre a une echelle de 0 a 130, et affiche la valeur positif et negatif de mon spectre. (cf piece jointe)

Voici mon code:

Code :

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

1
2
3
4
5
6
7
8
9
10
 
% Je recupere tous d'abord la sortie de mon bloc simulink
rtoIn = get_param([bdroot '/Measurement_tool/Inputs_settings/Am1'],'RuntimeObject');
InSignal=rtoIn.OutputPort(1).Data;
 
%Je calcul ensuite son spectre
InSpctr=abs(2*fft(InSignal)/length(InSignal));
 
%Je trace mon spectre
plot(squeeze(InSpctr(:,1)))

Si quelqun a une solution a mon probleme,
Merci

[phryte]

Salut.

mon spectre a une echelle de 0 a 130, et affiche la valeur positif et negatif de mon spectre

Pour avoir une échelle semi-log :

Code :

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

1
2
3
4
5
6
7
NFFT = length(input);
Y = fft(input,NFFT);
f = linspace(0,Fs/2,NFFT/2);
Y_plot = abs(Y(1:floor(NFFT/2)));
figure; semilogx(f, 2*Y_plot) 
xlabel('Frequency [Hz]');
ylabel('Amplitude');grid

[laurent12]

Merci pour ta reponse phryte, mais je ne veut pas une echelle semi-log.
Deplus, avec ta methode, et en utilisant 'plot' au lieu de 'semilogx', la frequence qu'affiche mon spectre n'est pas tout a fait exacte, et l'amplitude est de 115, pour un sinus d'amplitude 1.


Je n'arrive donc pas a visualiser mon veritable spectre...

[phryte]

Salut.
Pour plot je fais comme ça (variable en fonction du signal) :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
f = F*(0:256)/512;
Y = fft(y,512);%FFT
figure(2)
plot(f,abs(Y(1:length(f))))
grid
Pyy = Y.* conj(Y) / 512;%power spectrum, 
figure(3)
plot(f,Pyy(1:257))
title('Frequency content of y')
xlabel('frequency (Hz)')
grid

Tu peux aussi tracer le spectre sous simulink.