Mon algo pour Fourier est-il corrêcte ?

[bnoch]

Bonjour, voici algorithme que j'utilise pour transformer un signal périodique en séries de coefficients puis en signal de nouveau.

Malheureusement j'obtiens des résultats un peu aléatoire au niveau de la jonction bout-à-bout d'une seule période générée.
(voir post section informatique "Coefficients de Fourier" pour le code).

Algo :

A partir d'une seule période d'un signal périodique formé de différentes harmoniques de phase et amplitude quelconques.
(fréquence d'échantillonnage Fe = Nombre de samples N)

-Calcul d'un niveau identité: Somme des amplitudes fournie divisé par leurs nombre.
-Pour chaque harmonique :
-Calcul de sa fréquence F : Fe / N * rang de l'harmonique.
-Pour chaque sample, préparation des coefficients a et b :
-Calcul de la phase Ph : F * (2Pi / Fe) * index du sample.
-Addition à a de Valeur du sample * Cos(Ph).
-Addition à b de Valeur du sample * Sin(Ph).
-a et b sont multipliés par 2/Fe.
-l'amplitude de cette harmonique est la racine carrée de la somme des carrés de a et b.
-la phase est l'arc tangente de b/a.

A partir de cela pour reformer le signal:
-Valeur d'un sample v d'indice t :
-Pour chaque composante :
-Calcul de la phase Ph : F * 2Pi / Fe * t
-Addition à v de a * Cos(Ph).
-Addition à v de b * Sin(Ph).