Auto-corrélation sous Matlab

Bonjour
je dois faire le traitement d'un signal notamment en utilisant Fourier et compagnie. En fait, j'applique la fonction fft à l'auto-correlation de pression. En gros si P représente un vecteur de pressions, mon programme sur une version récente de Matlab donne

P=P-mean(P)
Auto_corr = xcorr(P)
Spectre = fft(Auto_corr)
Puissance_spect = sqrt(Spectre.*Spectre)

et ensuite je plotte la puissance spectale comme il faut.

Mais le problème, ce que je dois développer sur une version de Matlab qui n'a pas la fonction xcorr, qui a en revanche fft.

Ma question : comment faire l'équivalent du programme ci-dessus, avec une version sans xcorr, i.e. quoi faire pour avoir le même résultat ?

Merci d'avance !