problème avec la courbe de spectres obtenue à travers un code

Bonjour,
J'essaye de comprendre un code pour l'utiliser dans mon étude.
Dans une partie du code je n'arrive pas comprendre le résultat obtenu sur la courbe: 4 courbes alors que je pense que le code ne dessine que 2.
J'ai déjà utilisé les fonctions semilogy et line mais je ne comprends pas d'où viennent les 2 autres courbes obtenues sur la figure dans ce cas !!!

voilà mon code
pout et pout2 sont mes deux signaux

Code:

---

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17

figure(24) clf ysig=real(pout); yerr=real(pout)-real(pout2); nfft=389; bd=isnan(ysig); gd=find(~bd); bd([1:(min(gd)-1) (max(gd)+1):end])=0; ysig(bd)=interp1(gd,ysig(gd),find(bd)); [Pxs,F]=spectrum(ysig(isfinite(ysig)),nfft,1,[],ceil(nfft/2)); %[Pxs,F]=pmtm(ysig(isfinite(ysig)),4,4096,1); yerr(bd)=interp1(gd,yerr(gd),find(bd)); [Pxe,F]= spectrum(yerr(isfinite(ysig)),nfft,1,[],ceil(nfft/2)); %[Pxe,F]=pmtm(yerr(isfinite(ysig)),4,4096,1);   semilogy(F,Pxs); line(F,Pxe,'color','r');

---