Discussion :

[zizo08]

Bonjour,
J'ai une question et j'ai besoin svp de votre aide.
J'ai une courbe et je veux déterminer les coordonnées du point ou la fonction atteint 95% de sa valeur finale ( la valeur ou la fonction devient constante).

est ce que je peux implémenter par matlab un tel problème.
pour être franche j'en ai aucune idée comment pourrais je faire ça.

merci

[ricololo]

Salut,
donc tu as un vecteur x et un vecteur y=f(x) qui tend vers une valeur contante.
Tu peux peut-être procéder ainsi:

Code :

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

1
2
 
l=find(y==0.95*y(end));

Les coordonnées de ton point sont donc (x(l);y(l))

[zizo08]

merci pour votre idée.
est ce qu'on ne peut pas l'indiquer dans la figure: je trace f(x) et sur la même
figure indiqué ce point.
Est ce possible ça en matlab?

[ricololo]

Si, tu peux le faire, c'est même très facile.
Tu m'as l'air de débuter en matlab, faudrait peut être que tu lises quelques tutoriels !!
Je t'aide tout de même (même si ce n'est pas dans la politique du forum de 'pondre' les réponses!):

Code :

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

1
2
 
plot(x,y,'b',x(l),y(l),'ro')

La courbe apparaît en bleu et le point en rouge!!

[zizo08]

merci beaucoup pour ton aide.
je tout à fait débutante en matlab.
mais ce n'est pas ce que je cherche, j'ai mal exprimé ma question:la fonction tend à atteindre une valeur constante, je veux limiter la zone sur laquelle ou elle va devenir constante par deux axes par exemple verticales ( je veux le faire grahiquement).

merci une autre fois.

[Caro-Line]

je veux le faire grahiquement

Les images étant plus parlantes que les mots, si tu nous faisais justement une image de ce que tu veux obtenir ?
Parce que personnellement je n'ai rien compris à ce que tu veux.

Ton problème de départ (déterminer les coordonnées (x,y) telles que f(x,y)=95% de f(xfin,yfin)) est purement algorithmique.

[zizo08]

vous avez raison.
je veux déterminer le point algorithmiquement et en suite encadé la zone 'constatnte' qui commence à partir de ce point par de axes verticales.

espérons que je viens d'exprimer cette fois ci bien momn problème.

[Caro-Line]

je veux déterminer le point algorithmiquement

As-tu essayé la méthode de ricololo ?
(bien qu'à mon avis il faille surement rajouter un petit epsilon car il se peut qu'on n'ait jamais exactement la valeur, voir la Pourquoi 0.3-0.2-0.1 est-il différent de 0 ? vers la fin)

et en suite encadé la zone 'constatnte' qui commence à partir de ce point par de axes verticales.

Tu ne veux vraiment pas faire un dessin (même tout moche sous Paint), je t'assure que ce serait plus clair...
De plus il doit y avoir un problème de vocabulaire car ta courbe est déjà dans des axes et tu ne peux pas avoir des axes dans des axes

[ricololo]

Tentative de solution à ton problème (pas sûr de l'avoir bien compris):
si on considère la fonction f(x)=100*(1-exp(-x/10)) qui tend vers 100 en plus l'infini (et déjà bien avant!!):

Code :

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

1
2
3
4
5
6
 
x=[0:.01:100];
y=100*(1-exp(-x/10));
l=find(y<.95001*y(end) & y>.94999*y(end));
plot(x,y) % en bleu
plot([x(l),x(l)],[0,y(l)],'r') % en rouge

c'est ce dont tu as besoin? si c'est le cas, il existe surement une façon plus élégante de tracer cet 'axe' vertical!

[zizo08]

merci beaucoup ricololo.
c'est bon maintenent ça mrache.

grand merci.
bonne journée

[Jerome Briot]

Code :

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

1
2
3
4
5
6
 
x=[0:.01:100];
y=100*(1-exp(-x/10));
l=find(y<.95001*y(end) & y>.94999*y(end));
plot(x,y) % en bleu
plot([x(l),x(l)],[0,y(l)],'r') % en rouge

Et encore une fois, le FIND ne sert à rien ici :

Code :

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

1
2
3
4
5
6
7
 
x=[0:.01:100];
y=100*(1-exp(-x/10));
l=y<.95001*y(end) & y>.94999*y(end);
plot(x,y) % en bleu
hold on
plot([x(l),x(l)],[0,y(l)],'r') % en rouge

[ricololo]

BONJOUR Mr Dut,
excuse moi pour cette absurde FIND!

[Jerome Briot]

Mr Dut,

Mr Dut... c'est nouveau ça

[Jerome Briot]

Mr Dut est de retour... avec la fonction FIND

Code :

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

1
2
3
4
5
6
7
figure
x=0:.01:100;
y=100*(1-exp(-x/10));
idx = find((y-95)>=0,1,'first');
plot(x,y) % en bleu
hold on
plot([x(idx),x(idx)],[0,y(idx)],'r') % en rouge

D'ailleurs si on compare les deux codes :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
>> y(idx)
 
ans =
 
   95.0013
 
>> y(l)
 
ans =
 
   94.9963

[zizo08]

Merci beaucoup Mr.Dut.
votre code est parfait.

meri

[zizo08]

bonjour,

je vais poser une autre question dans le même context.
si j'ai un vecteur x et je vais le dessiner plot(x), est ce que je peux utiliser la même stratégie pour obtenir sa valeur limite càd ja fais comme suit:*

Code :

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

1
2
3
4
5
6
7
8
9
 
 
a=0:12:840;
 
idx = find((x-95)>=0,1,'first');
%
plot(a,x,'r')
hold on
plot([a(idx),a(idx)],[0,x(idx)],'r')

merci

[ricololo]

Salut,
bien sûr! Ce qu'on t'as proposé avec des vaecteurs x et y, tu peux le faire avec des vecteurs a et x!!
C'est bien ta question? sinon, j'ai mal compris.

[zizo08]

même si x est construit à partir d'une boucle?

Code :

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

1
2
 
x=[x x(i)]

je pose cette question parce que je n'ai pas eu le bon résultat que je cherche.
si joit la figure de vecteur x dont je parle et je veux déterminer sa valeur limite.


merci pour l'aide

[Jerome Briot]

Tu devrais vraiment lire ce tutoriel pour te familiariser avec les matrices : Introduction à la gestion des matrices

[zizo08]

j'arrive à comprendre le problème.
la valeur de idx est un réel et si je fais x(idx), le point ne s'affiche pas car idx doit être un entier ou une valeur logique.
que dois je faire?

merci