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 : Sélectionner tout - Visualiser dans une fenêtre à part
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');