Transformée et série de Fourier
Bonjour a tous,
j’espère que vous allez pouvoir m'aider :)
Je dois utiliser la transformée de Fourier discrète (TDF ou FFT) pour évaluer les coefficient du développement en série de Fourier d'un signal.
La série de Fourier, ça va, la transformée de Fourier aussi, par contre je viens a peine de commencer Matlab ...
donc j'ai le signal suivant :
Code:
1 2 3 4 5
| f1 = 3 ;
t = 0:0.01:1;
s1 = 1 + 2 * sin (2 * pi * f1 * t);
subplot(3, 1, 1), plot(t, s1); |
jusque la ok, par contre après, je dois faire un vecteur x avec les coefficients et la j'y arrive pas.
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
|
N = numel(t);
x = zeros(1, N);
support = zeros(1, N);
for k = 0:N - 1
x(k + 1) = 0;
for n = 0:N - 1
x(k + 1) = x(k + 1) + s1(n + 1) * exp(-1i * 2 * pi * k * n / N);
end
support(k + 1) = k / N;
end
% Module
subplot(3, 1, 2), plot(support, abs(x));
% Argument
subplot(3, 1, 3), plot(support, angle(x)); |
Voila ... (on se moque pas s'il vous plait :aie:) j'arrive pas après au moment de calculer les coefficients.
Je n'arrive pas a reproduire cette formule : http://upload.wikimedia.org/wikipedi...633d33bdea.png
Merci pour votre aide !