on a définit f comme étant
janvier T1 sur t=0 à t=1/12
T2 sur t=1/12 à 2/12
..
decembre T12 sur t=11/12 à 12/12

j'y ai pas prêté attention, mais cette fonction n'est pas jolie, car il y a un saut entre la valeur T1 en t=0 et la valeur T12
en t=1, et les series de fourier n'aiment pas ce genre de saut qui casse la périodicité.

la manière de m'en sortir, est de redéfinir le signal température comme étant
T1 sur t=0 à t=1/24
T2 sur t= 1/12 - 1/24 à t= 1/12+ 1/24
...
T12 sur t= 11/12-1/24 à t=11/12+1/24
T1 sur t= 1-1/24 à 1

ma fonction est périodique de période = 1. la fonction calculée dans la procédure, va etre proche de celle la. Donc je comprends que le décallage soit de 1/24, un demi mois.

maintenant, pourquoi tes fonctions sont décalées de 1 mois, je ne peux pas le dire. est ce du au fait que tu omis un mois?

sous matlab, tu construira un vecteur de 13 points, ou le 13 ieme est egal au premier,
Tu fais la tfd, qui te donne c_n pour n allant de -6 à 6 (donc 13 valeurs)
argument , module de ce machin.
Je ne sais plus s il y a une fonction qui de donne la somme de dirichlet directement.