Discussion :

[diego45]

Bonjour, je viens vous demander votre avis et aide pour améliorer un morceau de mon code :

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
 
 
  for k=1:29
 
        [n,p]=size(deltaspreadint{k,1}(:,:));
 
        for pp=n:-1:20
            for o=1:p
            if isempty(find(deltaspreadint{k,1}(pp,o)==0))==1 && sum(deltaspreadint{k,1}((pp-10:pp-1),o))==0
 
                deltaspreadint{k,1}(pp,o)=0;
            end
            end
        end
    end

Ce code est très long à tourner !

Le but de la manipulation permet d'enlever les valeurs "abérantes" dans le calcul de delta (Xt+1-Xt). En effet lorsque Xt=0 on à un delta énorme, comparément aux autres valeurs.
Ce pronlème vient des valeurs manquantes dans ma série, qui sont remplacés par 0.

Imaginons un vecteur Xt=

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
 
0
0
0
0
0
10
11
12
13

On à alors un vecteur delta :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
 
0
0
0
0
10
1
1
1

Le 10 représente dans ce que je souhaites faire après, un point abérant;
le code que j'ai présenté avant, cherche dans le vecteur si la valeur pp est non nulle si oui, et si les 10 valeurs avant sont nulles alors il considére que ce pp est à remplacé par 0.

Mais c'est très long, avait vous une idée de la façon de traiter ce problème ? merci

[diego45]

J'ai testé de gérer ça directement au moment du calcul du delta :

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
18
 for k=1:29
     B=['Plot',num2str(k)];
     Matrice_delta_int=diff(spreadint2{k,1}.(B)(:,:));
     [n,p]=size(Matrice_delta_int);
 
     for i=n:-1:20
            for j=1:p
            if isempty(find(Matrice_delta_int(i,j)==0))==1 && sum(Matrice_delta_int((i-10:i-1),j))==0
 
               Matrice_delta_int(i,j)=0;
            end
            end
     end
 
 
     deltaspreadint{k,1}= Matrice_delta_int(:,:);
 
 end

Je vais voir si ça permet d'aller plus vite

[diego45]

C'est bon j'ai supprimé en matriciel les valeurs supérieur à 10 écart type

[Gooby]

Bonjour,

Je n'ai pas exactement saisi ce que tu veux faire mais si tu cherches les delta de grandes valeurs sans boucle, tu peux toujours faire comme ceci je suppose:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

idx = find(abs(deltaspreadint(1:end-1)-deltaspreadint(2:end))>X); % valeur X à ne pas dépasser