Discussion :

[jena]

Salut
J'ai un vecteur qui contient des NaN. J'aime bien determiner les positions de début de chaque Groupe de NaN et son fin et que est le groupe le plus longue et son debut et son fin
Exemple ce vecteur
x=[ 1 3 5 NaN NaN 6 7 8 9 11 NaN NaN NaN NaN NaN NaN 13 14 15 NaN NaN NaN NaN 18 NaN 15 NaN NaN NaN NaN ]
on a 5 groupes
le groupe le plus longue commence à pos 11 et se termine à pos 16

Code :

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

1
2
3
4
5
6
7
8
 
for i = 1:length(x)
if isnan (x[i])
posdeb  = i;
elseif not (isnan (x[i]))
posfin  = i;
end
end

Mais ce code ne donne rien il manque la mise à jour de pos

[jena]

Voici un autre essai mais pas vraiment de resultat

Code :

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

1
2
3
4
5
6
7
8
9
10
11
pos = [];
 k=1
for i = 1:length(x)
        if isnan (x(i))
            pos(k)  = i; 
            k=k+1;
        elseif not (isnan (x(i)))
            pos(j)  = i;
            k=k+1;
        end
end

x =

Columns 1 through 23

1 3 5 NaN NaN 6 7 8 9 11 NaN NaN NaN NaN NaN NaN 13 14 15 NaN NaN NaN NaN

Columns 24 through 30

18 NaN 15 NaN NaN NaN NaN

et

pos =

Columns 1 through 23

0 0 0 4 26 0 0 0 0 0 11 12 13 14 15 16 0 0 0 20 21 22 23

Columns 24 through 30

0 25 0 27 28 29 30
aide SVP

[Newenda]

tout ça semble bien compliqué et tu dois pouvoir t'en sortir sans boucle.

Voici un code qui joue avec les indexages et la propriété de isan de mettre de 1 ou 0 :
L'idée générale c'est qu'avec ce vecteur, si tu soustrait l'indice (i+1) par l'indice (i).. Au premier NaN tu auras 1-0 = 1, ensuite pour le groupe de NaN 1-1=0 et enfin pour le dernier NaN: 0-1 = -1.

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
clear all
close all

x=[ 1 3 5 NaN NaN 6 7 8 9 11 NaN NaN NaN NaN NaN NaN 13 14 15 NaN NaN NaN NaN 18 NaN 15 NaN NaN NaN NaN ] ;
x2=[x 0] ; % ajout d'une valeur à la fin (tout sauf NaN) pour l’agrandir ton vecteur d'un indice (servira après)
NaNvector=isnan(x2)  ; % remplacer les NaN par des 1 et le reste par des 0

NaNvector_shifted = [NaN NaNvector(1:end-1)]; % créer un vecteur identique au précédent mais décalé d'un indice

diff = NaNvector-NaNvector_shifted ; % faire la différence entre les deux


deb = find(diff==1)  % lorsque le calcul de diff vaut 1: début des NaN
fin = find(diff==-1)-1  % lorsque le calcul de diff vaut -1: fin des NaN + 1 indice

deb te donne les indices de tous les débuts de groupe de NaN
fin te donne les indices de toutes les fins des groupes de NaN

Code :

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

1
2
3
4
5
6
7
8
deb =
 
     4    11    20    25    27
 
 
fin =
 
     5    16    23    25    30

Voilà, je te laisse regarder ça.

[Invité]

Bonjour,

Une solution :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
x = [1 3 5 NaN NaN 6 7 8 9 11 NaN NaN NaN NaN NaN NaN 13 14 15 NaN NaN NaN NaN 18 NaN 15 NaN NaN NaN NaN];
 
start  = strfind([false ~isnan(x)], [false true]); % Positions de début
stop   = strfind([~isnan(x) false], [true false]); % Positions de fin
groups = arrayfun(@(d,f) x(d:f), start, stop, 'UniformOutput', false); % Formation des groupes
groupsSize = cellfun(@numel, groups); % Taille de chaque groupe
 
[maxSize,idx] = max(groupsSize); % Plus grande taille
start(idx) % Position de début 
stop(idx)  % Position de fin
groups{idx} % Groupe le plus long

[Jerome Briot]

Tu peux aussi t'inspirer de ce défi : [Défis][7] Et maintenant, la suite ! et de la solution donnée par Aurélien sur son blog

Code :

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

1
2
3
4
5
ix = isnan(x);
 
index = find(diff([ix,inf]));
 
data = [ix(index) ; [1 index(1:end-1)+1] ; diff([0,index])]

Ce qui donne avec l'exemple que tu as donné :

Code :

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

1
2
3
     0     1     0     1     0     1     0     1     0     1
     1     4     6    11    17    20    24    25    26    27
     3     2     5     6     3     4     1     1     1     4

Les résultats se lisent colonne par colonne. Chaque colonne correspond à un groupe de valeurs successives NaN ou non.

• La valeur de la première ligne à 1 signifie que la colonne en question correspond à un bloc de NaN consécutifs.
• La deuxième ligne donne l'indice dans x du début du groupe
• La troisième ligne donne le nombre d'élément dans le groupe

Donc dans ce cas, on voit bien que le plus grand groupe de NaN correspond à la colonne 4, avec un indice de début à 11 et qu'il comporte 6 éléments.