Discussion :

[IFREMER]

Aloha Matlabeuzs et Matlabeurs

J'ai un petit soucis... (et oui encore ).

J'ai une série de mesure, avec laquel je fais la moyenne et je retranche cette moyenne à mes mesures. Je trace cette soustraction avec un plot et maintenant je dois chercher les abscisses de la courbes où celle ci passe par 0...
Je n'ai pas la moindre idée pour procéder... je peux pas passer par l'équation de la courbe (trop compliquée..) et dans ma série de valeur ayant subis la soustraction il n'y a aucune valeur à 0...
Donc je suis un peu perdu

Je vous joins le graph pour que vous comprenniez mon dilemme

Merci d'avance pour votre aide

[Newenda]

Code :

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

indice_des_0 = find(X==0)

un petit exemple :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
>> X=[5 0 2 4 8 0 4 5]
 
X =
 
     5     0     2     4     8     0     4     5
 
>> indice_des_0 = find(X==0)
 
indice_des_0 =
 
     2     6

EDIT : désolé je n'avais pas lu qu'aucune valeur = 0, dans ce cas il faut récupérer les valeurs, il me semble que y'a une aide sur la FAQ

[chris05]

Le probleme est que le vecteur solution de IFREMER ne contient pas forcement les zeros ...

Je ne crois pas qu'il est possible d'identifier les zeros d'une fonction a partir d'une figure (à vérifier avec les pro de Matlab présents sur ce forum (Duf, Caro?? )).

Néanmoins si tu possèdes l'equation de ta fonction, tu peux toujours utiliser la commande solve :

Code :

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

1
2
 
zeros = solve ('function = 0', 'variable');

[IFREMER]

Aie aie, je n'ai pas l'équation de la courbe malheursement... j'espère qu'avec le graph c'est possible...

[Caro-Line]

Dans un 1er temps ce serait bien que tu nous donnes une idée de tes données.

Après la question c'est aussi de savoir quelle précision tu veux et quelle méthode tu dois utiliser.

Je ne comprends pas dans ton 1er post tu dis que l'équation de la courbe est trop complexe et dans ton dernier que tu ne l'as pas

Personnellement à part essayer de lire les valeurs sur le graphe (absolument pas scientifique ni précis) ou résoudre l'équation, je ne vois pas trop d'autres solutions (l'utilisation d'INTERP1 dans ce sens me parait difficile mais tu peux regarder quand même).
Mais bon je n'y connais pas grand-chose dans ce domaine-là.

Peut-être tenter un DIFF aussi pour voir où ont lieu les changements de signe puis procéder par dichotomie ?

[chris05]

Si tu ne possede pas l'equation de ta fonction, tu peux toujours essayer d'en trouver une se rapprochant de ta serie de mesures avec une formule d'interpolation matlab (par exemple polyfit) puis calculer les zeros de cette fonction.

Code :

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

1
2
 
f=polyfit(x,y,n)

Cette commande interpole ta serie de mesure sous forme d'une fonction polynomiale d'ordre n. Vu la complexité de ta serie, utilise au minimum n=3.

Il existe d'autre fonction de ce type mais je ne les ai pas en tête.

Le probleme est que tu auras des erreurs numériques, les zeros obtenus ne correspondront peut etre pas exactement aux zeros de ta serie de mesures. Cela te donnera néanmoins une bonne idée de leur localisation.

Sinon je ne vois que la recherche des zeros 'manuelle' en vue de ton graphique, mais cela va etre long et pas bien bon pour tes yeux

[Newenda]

Le probleme est que tu auras des erreurs numériques, les zeros obtenus ne correspondront peut etre pas exactement aux zeros de ta serie de mesures. Cela te donnera néanmoins une bonne idée de leur localisation.

Sinon je ne vois que la recherche des zeros 'manuelle' en vue de ton graphique, mais cela va etre long et pas bien bon pour tes yeux

De toute façon si il n'y a aucune valeur en 0, il y aura forcément une erreur, ça on ne peut pas y échapper, il ne s'agira que d'une estimation, quelque soit la méthode.
Je ne sais pas si la recherche "manuelle" est mieux adapté ici que ta première idée ! Autant rester sur un régression, c'est quand même plus propre et permet en outre d'avoir un modèle de ton étude : ce qui est extrêmement intéressant !

[tug83]

Si tu as la Spline Toolbox, utilise fnzeros
sinon tu as cette soluce :
How do I find all the zeros of a function?

[IFREMER]

Excusez moi si je n'ai pas été clair voici quelques une de mes valeurs:

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
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
 
>>mois 
 
mois=
 
    1
    2
    3
    4
    5
    6
    7
    8
    9
   10
   11
   12
   13
 
>> sstr
 
sstr =
 
   -0.5486
   -0.6258
   -0.6165
   -0.3113
    0.2365
    0.2242
    0.3509
    0.8701
    1.1467
    1.2606
    0.8623
    0.1576
   -0.8050

Et je n'ai pas l'équation de la courbe car celle ci est trop compliquée à trouver mais je vais essayer les méthodes qu'on me propose merci bcp.

EDIT:
Recoucou,

Ayant réfléchis en mangeant ^^, j'ai un peu travaillé sur la solution que me proposer caroline et cela m'ira
J'ai donc fait un interp1 en interpolant mes données à 0.001 près et j'ai des résultats assez convaincants (pas tout à fait correctes mais presque à 0.01)

Je vous remercie encore pour le coup de main et je mets le code en ligne dès que j'ai terminé ^^

Code :

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

1
2
3
xi = 1:.001:108;
yi = interp1(mois,sstr,xi);
plot (xi,yi)

Bon par contre je suis obliger de chercher les solutions à la main car aucune valeur n'est à zero... mais bon c'est déjà ça ^^

[IFREMER]

Petite question qui pourrais m'aider dans la fonction que je suis en train de réaliser pour mon pbm:

est ce qu'il existe une fonction matlab qui detecte le signe (+ ou - ) d'un nombre?

Merci

[duf42]

Code :

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

1
2
 
help sign;

[IFREMER]

... décidément je passe pour quoi ??? ... enfin bref

Donc voic le code que j'ai mis en place c'est le mieux que j'ai réussis à faire donc si vous avez des retouches pour qu'il soit encore plus performant merci à vous ^^

Code :

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

1
2
3
4
5
6
7
8
9
xi = 1:.001:108
yi = interp1(mois,sstr,xi);
Xnul=[];
Ynul=[];
for a = 1:107000;
if sign(yi(a)) + sign(yi(a+1))==0;
Xnul=[Xnul;xi(a)];
Ynul=[Ynul;yi(a)];
end
end

EDIT:
Euh ça marche pas il me fait un décalage... l'indéxation est bonne pourtant non?? mm si yi et xi sont des matrices colonnes

[Caro-Line]

Qu'entends-tu par décalage ?

Sinon ton code tel qu'il est devrait planter car le compteur a devrait aller jusqu'à length(yi)-1 (107000 dans ton cas) et non pas length(yi).

[IFREMER]

Qu'entends-tu par décalage ?

Sur le graph le changement de signe s'effectue trois valeurs avant celle affichée... Si besoin je peux vous envoyer le fichier mat pour que vous regardiez

--> Réglé j'avais dut laisser une valeur sans faire attention

Sinon ton code tel qu'il est devrait planter car le compteur a devrait aller jusqu'à length(yi)-1 (107000 dans ton cas) et non pas length(yi).

Oui dsl je l'avais remarqué mais oublié de le changer sur le post

Sinon j'ai essayer de la mettre en fonction pour regarder avec le DEBUG mais là ça fonctionne plus du tout....

Code :

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

1
2
3
4
5
??? Error using ==> interp1
NaN is not an appropriate value for X.

Error in ==> detection at 4
yi = interp1(mois,sstr,xi);

Voici la fonction que j'ai faite :

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
function[]=detection()
load MaxetMin.mat
xi = 1:.001:108;
yi = interp1(mois,sstr,xi);
Xnul=[];
Ynul=[];
plot(xi,yi)
for a = 1:107000;
 
    if sign(yi(a)) + sign(yi(a+1))==0;
        xi(a)
        yi(a)
        Xnul=[Xnul;xi(a)];
        Ynul=[Ynul;yi(a)];
    end 
 
end

Je comprend plus rien !!

[Caro-Line]

Pour le message d'erreur c'est étrange il te signale que tu as des NaN (Not A Number) dans ton vecteur mois...
Tu devrais le vérifier.

[IFREMER]

Je viens de vérifier, j'ai aucun NaN dans mes matrices....

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
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
 
>>mois
 
mois =
 
     1
     2
     3
     4
     5
     6
     7
     8
     9
    10
    11
    12
    13
    14
    15
    16
    17
    18
    19
    20
    21
    22
    23
    24
    25
    26
    27
    28
    29
    30
    31
    32
    33
    34
    35
    36
    37
    38
    39
    40
    41
    42
    43
    44
    45
    46
    47
    48
    49
    50
    51
    52
    53
    54
    55
    56
    57
    58
    59
    60
    61
    62
    63
    64
    65
    66
    67
    68
    69
    70
    71
    72
    73
    74
    75
    76
    77
    78
    79
    80
    81
    82
    83
    84
    85
    86
    87
    88
    89
    90
    91
    92
    93
    94
    95
    96
    97
    98
    99
   100
   101
   102
   103
   104
   105
   106
   107
   108
 
>>sstr
 
sstr =
 
   -0.5486
   -0.6258
   -0.6165
   -0.3113
    0.2365
    0.2242
    0.3509
    0.8701
    1.1467
    1.2606
    0.8623
    0.1576
   -0.8050
   -1.2233
   -1.1546
   -1.0485
   -0.6807
   -0.4645
   -0.1124
    0.5952
    1.0818
    1.0535
    0.6184
    0.1269
   -0.8370
   -1.4343
   -1.4317
   -1.3190
   -0.8238
   -0.6650
   -0.2506
    0.3245
    0.6119
    0.8988
    0.6313
   -0.4930
   -1.1365
   -1.5806
   -1.2757
   -1.1309
   -0.6641
   -0.3049
    0.0158
    0.3981
    1.0065
    0.9795
    0.8788
    0.2215
   -0.2801
   -1.0067
   -0.8466
   -0.6783
   -0.6609
   -0.4435
   -0.2466
    0.2959
    0.8417
    1.1158
    0.5450
   -0.2592
   -0.7573
   -1.2999
   -1.2755
   -0.7971
   -0.7414
   -0.5260
   -0.1849
    0.6993
    1.2744
    1.3727
    0.8536
    0.2579
   -0.4618
   -0.9230
   -1.0133
   -0.4296
   -0.5178
   -0.3073
    0.0197
    0.8132
    1.3018
    1.3595
    0.8055
    0.3128
   -0.1252
   -0.3387
   -0.2154
    0.0168
    0.7801
    0.7554
    0.6667
    1.1430
    1.4477
    1.4770
    0.8460
    0.3681
   -0.4969
   -0.9319
   -1.0935
   -0.6348
    0.0404
   -0.1569
    0.0091
    0.6737
    1.2661
    1.3825
    1.0872
    0.2080

et ce que je comprend pas c'est que quand je réalise cette série d'instruction directement sur MatLab ça fonctionne...

EDIT: Sinon pour être plus précis dans la valeur des coordonnées du zéro j'ai ajouté une moyenne pour avoir encore un peu moin d'erreur (en réalisant les commandes sur la "fenêtre de commande" ):

Code :

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

1
2
3
4
5
6
7
8
9
10
xi = 1:.001:108;
yi = interp1(mois,sstr,xi);
Xnul=[];
Ynul=[];
for a = 1:107000;
if sign(yi(a)) + sign(yi(a+1))==0;
b =(xi(a)+xi(a+1))/2;
c=(yi(a)+yi(a+1))/2;
Xnul=[Xnul;b]; Ynul=[Ynul;c];
end
end

[Caro-Line]

Pour être sur il faudrait vérifier juste après avoir chargé ton fichier MaxetMin.mat :

Code :

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

1
2
3
4
5
function ...
load MaxetMin.mat
if any(isnan(mois))
   error('il y a des NaN dans mois !!');
end

[IFREMER]

Pour être sur il faudrait vérifier juste après avoir chargé ton fichier MaxetMin.mat :

Code :

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

1
2
3
4
5
function ...
load MaxetMin.mat
if any(isnan(mois))
   error('il y a des NaN dans mois !!');
end

Je viens de le faire et en effet il me mets qu'il y a des NaN pourtant quand je regarde avec le debug étape par étape au moment où il charge mois je vérifie les valeurs et il n'y a pas de NaN... je montre

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
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
>> detection()
2   load MaxetMin.mat
3   if any(isnan(mois))
K>> mois

mois =

     1
     2
     3
     4
     5
     6
     7
     8
     9
    10
    11
    12
    13
    14
    15
    16
    17
    18
    19
    20
    21
    22
    23
    24
    25
    26
    27
    28
    29
    30
    31
    32
    33
    34
    35
    36
    37
    38
    39
    40
    41
    42
    43
    44
    45
    46
    47
    48
    49
    50
    51
    52
    53
    54
    55
    56
    57
    58
    59
    60
    61
    62
    63
    64
    65
    66
    67
    68
    69
    70
    71
    72
    73
    74
    75
    76
    77
    78
    79
    80
    81
    82
    83
    84
    85
    86
    87
    88
    89
    90
    91
    92
    93
    94
    95
    96
    97
    98
    99
   100
   101
   102
   103
   104
   105
   106
   107
   108

K>> %---Pas de NaN dans mois
4      error('il y a des NaN dans mois !!');

Je comprend pas... une explication??.... ou un fusil...

[Caro-Line]

Alors en ligne de commande :

Code :

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

1
2
3
load MaxetMin.mat
idx=find(any(isnan(mois))
mois(idx)

[IFREMER]

J'ai donc rajouté la ligne de code voila le résultat :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
detection
 
idx =
 
     1
 
??? Error using ==> mois
Too many input arguments.
 
Error in ==> detection at 4
mois(idx)

ça veut dire qu'il y a un NaN caché ??? O.o