Discussion :

[michastro]

Bonjour,
Je suis débutant sur Mathlab et je cherche juste à tracer la fonction 0.7sin(50HZ)+sin(120Hz) puis sa FFT.
Il y a visiblement un soucis dans ma partie abscisse pour la FFT. Voici mon source:

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
 
Fe = 3000;            % Frequence echantillonage
T = 0.1;             % Largeur fenetre en s
t = (0:1/Fe:T);        % vecteur temps de 0 à Largeur
S = 0.7*sin(2*pi*50*t) + sin(2*pi*120*t);
% plot(1000*t,S);
subplot(121)
plot(1000*t,S); % *1000 ca on veut des ms
title('Signal 0.7sin(50HZ)+sin(120Hz)');
xlabel('t (milliseconds)')
ylabel('X(t)')
Y=fft(S);
P=abs(Y);
F=linspace(-Fe*T,Fe*T,Fe*T+1);
subplot(122);
plot(F,P);

Merci de votre aide

[phryte]

Bonjour,

Quelques petites modifications :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
clear
Fe = 3000;            % Frequence echantillonage
T = 0.1;             % Largeur fenetre en s
t = (0:1/Fe:T);        % vecteur temps de 0 à Largeur
S = 0.7*sin(2*pi*50*t) + sin(2*pi*120*t);
% plot(1000*t,S);
subplot(121)
plot(1000*t,S); % *1000 ca on veut des ms
title('Signal 0.7sin(50HZ)+sin(120Hz)');
xlabel('t (milliseconds)')
ylabel('X(t)')
grid
Y=fft(S);
P=abs(Y);
%F=linspace(-Fe*T,Fe*T,Fe*T+1);
F=linspace(0,Fe,length(S));
subplot(122);
plot(F(1:20),P(1:20));
grid