Discussion :

[adnino]

Bonjour,
j'utilise ce code pour le calcul de la FFT de mon signal mais je ne sais pas si c'est correcte ou pas!
quelqu'un peut s'assurer de mon code ou de me proposer un autre code plus correcte
Un grand merci merci

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
20
21
22
 
tinc=1e-12;
tmax=1e-6;
t=0:tinc:tmax;
sigma=0.007;
c=3e8;
A0=20;
t0=100e-12;
 
S(t)=A0*exp(-(t-t0).^2/(2*sigma)^2)
Fs=1/tinc;
 
NFFT = 2^nextpow2(length(t));
 
Yfft = fft(S(t),NFFT)/tmax;
f = Fs/2*linspace(0,1,NFFT/2+1);
figure(2)
plot(f,abs((Yfft(1:NFFT/2+1))))
title('FFT signal emission')
xlabel('Frequency (Hz)')
xlim([0 5e11])
ylabel('|Y(f)|')

[Invité]

Bonjour,

As-tu pris la peine de tester ce code sur MATLAB ? Si oui tu aurais déjà ta réponse, vu qu'il génère une erreur...

Subscript indices must either be real positive integers or logicals.

t=0:tinc:tmax; contient des valeurs temporelles non entières, qui ne peuvent en aucun cas servir d'indices de S.

Et ce n'est pas par tmax que tu dois diviser la FFT mais le nombre d'échantillons.