Fourier transformé d'un signal RMN sujet

Bonjour,

Je travaille avec les valeurs omega d'un ensemble de 5000 spins, qui se promène entre deux états différents, sur une période dt.
Mon Workspace:
Nspins=5000
Ntimes=4096
dt=2.0000e-4 s
sw_h=5000

J'ai réussi à obtenir le signal en fonction du temps pour un spin, comme il m'est demandé.
Cependant on me demande par la suite d'obtenir le spectre de fréquence avec la transformée de Fourier pour l'ensemble des 5000 spins (que l'ensemble des pics se trouvent sur le même spectre.
Voià le programme que j'ai réalisé pour le début:

>> Omega1988=OmegaArray(1988,:);
>> X=cumsum(Omega1988);
>> Phi=X*dt;
>> S=exp(i*Phi);
>> figure(1)
>> subplot(2,1,1)
>> plot(t,Phi)
>> subplot(2,1,2)
>> plot(t,S)
Warning: Imaginary parts of complex X and/or Y arguments ignored
>> Y=cumsum(OmegaArray,2);
>> Phi2=Y*dt;
>> S2=exp(i*Phi2);
>> figure(2)
>> subplot(3,1,1)
>> plot(t,Phi2)
>> subplot(3,1,2)
>> plot(t,S2)
Warning: Imaginary parts of complex X and/or Y arguments ignored

A la fin de ces étapes je visualise correctement le spectre pour un seul spin, et sur la figure(2) le spectre pour les 5000 spins. Comment par la suite effectuer la transformée de Fourier avec Matlab:
>> nu1=sw_h*t;
>> nu2=nu1-(5000/(2+1));
>> S2=S2*0.5;
>> I2=fft(S2,5000);
>> I2=fftshift(I);
>> figure(3)
>> subplot(2,1,1)
>> plot(nu2,I2)
J'ai réalisé cette programmation mais le résultat est loin d'être satisfaisant. :/

Je vous remercie d'avance pour votre aide.