Discussion :

[membreComplexe12]

Bonjour,

J'ai une fonction périodique f(t), j'ai sur un tableau pour une periode (de 1secondes) les valeurs de f(t).

[voir piece jointe]

Je veux calculer les coefficients de Fourier mais je ne sais pas trop comment faire.

[voir en piece jointe l'expression des coeff de Fourier]

1°) pour le coefficients a0:
c'est facile il suffit de faire ==> a0=1/T * trapz(temps,f)

2°) pour les deux autres je cherche une méthode plus rapide:

==> Je pensais rajouter deux lignes:
ligneA: j'ecris le resultat de cos(n.2pi/T)
ligneB: j'ecris f*cos(n.2pi/T)

ensuite je fais an=2/T.trapz(ligneB,t)

---------------------------------------------

le probleme c'est que pour chaques coeff je dois rajouter 2lignes puis faire ensuite l'integration mais je me dit qu'il y a peut etre une methode plus rapide
que je ne connais pas?

[Jean Dumoncel]

Bonjour,

je ne comprends pas bien ce que tu as fait, peux-tu nous montrer ton code?

[membreComplexe12]

Disons que mon tableau est la matrice A. je fais une autre matrice B qui aura les valeurs des cosinus*fonction. et C qui sera les resultats des An

je ne traite dans ce cas que les coeff An car c'est la meme chose pour les Bn

Code :

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

1
2
3
4
5
6
7
8
9
 
T=1;
%si on s'interesse aux dix premiers harmoniques n=1:10
for n=1:10
B(n,:)=A(2,:)*cos(n*(2*pi)/T);
C(n)=2/T*trapz(A(1,:),B(n,:));
end
 
a0=1/T*trapz(A(1,:),A(2,:))

en faite je sais comment trouver les coefficients de fourier à partir de deux lignes (temps et valeur du signal) en faisant des boucles mais je voulais juste savoir si il existe une fonction toute prete sous matlab

et je voulais savoir si pour la fonction trapz on peut integrer dedans une multiplication et non directement nos deux vecteurs tout prets.

[Jean Dumoncel]

en faite je sais comment trouver les coefficients de fourier à partir de deux lignes (temps et valeur du signal) en faisant des boucles mais je voulais juste savoir si il existe une fonction toute prete sous matlab

Si elle existe, je ne la connais pas.

et je voulais savoir si pour la fonction trapz on peut integrer dedans une multiplication et non directement nos deux vecteurs tout prets.

Je ne sais pas si ça va répondre à ta question, mais tu peux par exemple faire :

Code :

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

1
2
3
for n=1:10
C(n)=2/T*trapz(A(1,:),A(2,:)*cos(n*(2*pi)/T));
end

[membreComplexe12]

Si elle existe, je ne la connais pas.

Code :

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

1
2
3
for n=1:10
C(n)=2/T*trapz(A(1,:),A(2,:)*cos(n*(2*pi)/T));
end

a oui en effet


merci beaucoup

[membreComplexe12]

je ne pense pas qu'une fonction existe pour faire se que je veux je vais donc le faire manuellement, merci quand meme pour l'aide