Discussion :

[aigle noir]

Bonjour,

pour faire l'analyse de données issues de capteurs pt 100 je dois les importer sur Matlab, chose qui a été faite

maintenant en scrutant ses données je trouve plein de fausses valeurs dues au système d'acquisition utilisé comme par exemple -200 degrés alors que ma plage de température varie entre 20 et 30 et aussi des NaN donc j'ai voulu faire un programme pour y remédier que je vous joins

Code :

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

1
2
3
4
5
6
7
function  x = replace(x)
 
    for k= 1:length(x)
        if x(k) == -200,Nan  
                x(k) = x(k-1);
        end
    end;

mais ça ne marche pas
pouvez vous m'aider a résoudre ce problème

merci

[Jerome Briot]

Essaie ceci :

Code :

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

1
2
3
4
5
6
7
function  x = replace(x)
 
for k = 1:length(x)
    if x(k)==-200 || isnan(x(k))
        x(k) = x(k-1);
    end
end

[aigle noir]

Tout d'abord merci Dut,

apres j'ai ressis à faire ce code

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
function  x = replace(x)
 
    for k= 2:length(x)
        I = x(k-1);
        if (x(k)>20) && (x(k)<30)  
            I = x(k);    
        else 
            x(k)=I;
 
 
        end
        x = x;
    end

mais le probleme que j'ai est que je dois conserver les valeurs dans le meme vecteur apres traitement (ecraser les valeurs originales), mais la fonction que j'ai faite ne me le permet pas.

une aide pour faire ça?

cordialement

[cedd70]

La solution de Dut ne fonctionnera pas si deux points consécutifs sont faux.
je ne vois pas trop de solution, pourquoi pas faire une moyenne du point précédent et suivant ( mais ça ne marchera quand même pas dans tous les cas)

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
 
function  x = replace(x)
 
    for k= 2:length(x)
        I = x(k-1);
        if ((I<20) || (I>30)|| (isnan(I))
            x(k)= (x(k-1)+x(k+1))/2;
        else 
            x(k)=I;
        end
        x = x;
    end

[Jerome Briot]

La solution de Dut ne fonctionnera pas si deux points consécutifs sont faux.

Code :

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

1
2
3
4
5
6
7
8
9
10
11
>> x = [0 0 -200 5 6 -200 -200 NaN 5 NaN 4]
 
x =
 
     0     0  -200     5     6  -200  -200   NaN     5   NaN     4
 
>> x = replace(x)
 
x =
 
     0     0     0     5     6     6     6     6     5     5     4

Ma solution ne fonctionne pas si le premier élément de x est faux (-200 ou NaN).
Il suffit de prendre ceci en considération avant d'entrer dans la boucle.

[cedd70]

Au temps pour moi, j'ai regardé trop vite.