Discussion :

[helaaa]

Bonjour à tous
Je veux remplacer un vecteur entier par des NaN en cas où j’ai dans ce vecteur 4 valeurs consécutive = NaN, sinon j’interpole entre les valeurs
J’ai fais ceci

Code :

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

1
2
3
4
5
6
7
8
test_nan=isnan(X);
            for k=1:length(X)
                if (test_nan(i)==1) && (test_nan(i+1)==1) && (test_nan(i+2)==1) && (test_nan(i+3)==1)
                    Yi(1:length(X))=NaN;
                else
                    idx=find(~isnan(X));
                    Yi=interp1(idx,X(idx),1:length(X));
                end

l'interpolation marche bien mais le problème c'est qu'on interpole aussi, même en ayant plus de 4 valeurs consécutive en NaN et ça remplace pas la ligne entière par des NaN

[soft001]

Bonjour à tous
Je veux remplacer un vecteur entier par des NaN en cas où j’ai dans ce vecteur 4 valeurs consécutive = NaN,

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
 
clear all
clc
t=[nan 1 2 nan nan nan nan];
t2=isnan(t);
m=0;
 
for i=1: numel(t2)
    if t2(i)==1
        m=m+1;
    else
        m=0;
    end
    if m==4
        t(1:end)=nan
        break
    end
end

[helaaa]

merci pour la réponse
je n'ai pas regardé la discussion car j'ai déjà trouver une solution
là voilà

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
 
            test_nan=isnan(X);
            test_nan=test_nan';
            tester= strfind(test_nan, [1 1 1 1]);
            trouvernan = isempty(tester);
            if trouvernan == 1
                idx=find(~isnan(X));
                Yi=interp1(idx,X(idx),1:length(X));
            else if trouvernan == 0
                    Yi(1:length(X))=NaN;
                end
            end

et merci encore une fois

[soft001]

Votre solution est plus optimisée.

[Jerome Briot]

La condition IF imbriquée n'est pas nécessaire ici :

Code :

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

1
2
3
4
5
6
if trouvernan == 1
    idx=find(~isnan(X));
    Yi=interp1(idx,X(idx),1:length(X));
else
    Yi(1:length(X))=NaN;  
end