Trouver les zéros

:lun: Aloha :lun: 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:(:cry:

Je vous joins le graph pour que vous comprenniez mon dilemme

Merci d'avance pour votre aide;)

Code:

---

indice_des_0 = find(X==0)

---

un petit exemple :

Code:

---

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

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:

---

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

---

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

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 :koi:

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 ?

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:

---

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 :D

Citation:

---

Envoyé par chris05

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 :D

---

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 !

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

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

Code:

---

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:

---

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 ^^

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

Code:

---

1 2	help sign;

---

;)

... décidément je passe pour quoi ??? :oops: ... 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:

---

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:cry: ç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

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).

Citation:

---

Envoyé par Caro-Line

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

Citation:

---

Envoyé par Caro-Line

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:P

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

Code:

---

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:

---

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 !! :mur:

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.

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

Code:

---

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...:calim2::calim2:

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:

---

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

---

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

Code:

---

1 2 3 4 5

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

---

Citation:

---

Envoyé par Caro-Line

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

Code:

---

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:

---

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... :arf::arf::arf: une explication??.... ou un fusil...:pan:

Alors en ligne de commande :

Code:

---

1 2 3

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

---

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

Code:

---

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)

---

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