Difficultés avec fzero

Bonjour,

J'ai quelques difficultés avec le M-File suivant :

Code:

---

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16

n=40000; D=[]; d=0.1892719; for k=0:20 format long e sol= fzero(@(x) myfun(x,d),0); d=sol+d; D=[D;d]; end d=2.4; for k=20:40 format long e sol= fzero(@(x) myfun(x,d),3); d=sol+d; D=[D;d]; end

---

La première boucle for fonctionne. Elle renvoie des valeurs positives.

Le problème est au niveau de la deuxième boucle for. Elle renvoie des valeurs positives et négatives par alternance. Or je voudrais qu'elle renvoie seulement des valeurs positives correspondant au zéro de la fonction au voisinage de 3.

Comment se fait-il que fzero renvoie des valeurs négatives alors que je lui précise de renvoyer des valeurs proches de 3 ?

Je vous remercie par avance pour votre aide.
Cordialement.

Salut

je te conseille de tracer sur un graphique ta fonction, sur une plage de x englobant les solutions qu'elle te retourne, et 3
tu devrais voir visuellement d'ou vient ton problème
(si ce n'est pas le cas, montre nous un de ces graphiques)

sinon tu peux aussi te retourner sur des fonctions d'optimisation te permettant d'ajouter des contraintes, comme fminbnd

Bonjour,

Merci pour votre réponse. Mettant rendue compte d'une erreur, j'ai modifié le programme. J'ai cependant toujours une erreur avec ce programme qui utilise la fonction fzero.

Pour voir d'où venait le problème, j'ai écrit deux programmes utilisant une fonction dont le zéro est facile à calculer en faisant un dl en dt=0.

La fonction est la suivante : f(dt)=dt*103000(exp(3*(t+dt))-1)-30.

Un programme calcule les valeurs de dt en lui donnant l'expression de dt et l'autre calcule les valeurs de dt en utilisant la fonction fzero.

Les deux programmes sont :

Programme calculant dt avec l'expression de dt :

Code:

---

1 2 3 4 5 6 7 8

dt=0; for k=0:10 a=103000*exp(3*(dt))*3; b=103000*(exp(3*(dt))-1); c=-30; delta=b*b-4*a*c; dt=dt+(-b+(delta)^(1/2))/(2*a) end

---

Programme calculant dt avec la fonction fzero :

Code:

---

1 2 3 4 5 6 7 8 9 10 11 12 13 14

H=[]; t=0; for k=0:10 format long e dt= fzero(@(dt) myfunpa(dt,t),0)+dt; t=dt; H=[H;dt]; end   function f = myfunpa(dt,t) %UNTITLED2 Summary of this function goes here % Detailed explanation goes here f=dt*(103000*(exp(3*(dt+t))-1))-30; end

---

Je n'ai cependant pas les mêmes résultats dans les deux programmes.

Quelqu'un pourrait-il me dire d'où vient mon erreur ?

Je vous remercie par avance,

Cordialement.

Même sans effectuer de boucle comme ci-dessous, les deux programmes ne donnent pas le même résultat.

Programme calculant dt avec l'expression de dt :

Code:

---

1 2 3 4 5

a=103000*exp(3*(dt))*3; b=103000*(exp(3*(dt))-1); c=-30; delta=b*b-4*a*c; dt=(-b+(delta)^(1/2))/(2*a)

---

Programme calculant dt avec la fonction fzero :

Code:

---

1 2 3 4 5 6 7 8 9

format long e dt= fzero(@(dt) myfunpa(dt),0);   La fonction est écrite dans le fichier à part ci-dessous : function f = myfunpa(dt,t) %UNTITLED2 Summary of this function goes here % Detailed explanation goes here f=dt*103000*(3*dt)-30; end

---

Je ne comprends vraiment pas d'où vient le problème.
La fonction fzero s'utilise-t-elle que sous certaines conditions ?

Vous remerciant par avance pour votre réponse.

salut

as tu suivi mon conseil ? a savoir de tracer ta fonction ?
a mon avis ca t'aiderai à y voir plus clair

fzero cherche une solution en partant d'un point initial, mais il ne trouvera pas la plus proche du point initiale :
en effet il la recherche en suivant des gradients, un peu comme si tu mettais une bille sur ta courbe au point x0 et qu'elle suivait la gravité

de manière générale quand on fait de l'optimisation sur des problèmes simples (cad avec 1 ou 2 variables, donc graphiquement représentables) il est bien de tracer à la fin une courbe de ta fonction avec le point de départ x0 et la solution x (et dans ton cas les solution x)

Bonjour,

J'ai tracé la fonction : f(dt)=dt*103000*(3*dt)-30.

Autour de 0, il y a deux zéros :

1) Pour f(0.0098)<0 et pour f(0.0099)>0.

2) Pour f(-0.0099)>0 et f(-0.0098)<0

La fonction fzero donne le deuxième zéro. Or je voudrais le premier zéro.

J'ai essayé la commande "Zero found in the interval: [-2,2]"
pour spécifier l'intervalle dans lequel la solution doit être trouvée. Cependant avec cette commande l'intervalle doit être symétrique.

Est-il possible de modifier cette commande pour pouvoir inscrire un intervalle non symétrique ?

Merci d'avance.

Salut,
Tu peux symétriser l’intervalle toi même : si tu cherches pour x près de a, cherche pour x-a proche de 0.

Bonjour,

Merci pour votre réponse.

En fait, la commande : "Zero found in the interval: [-2,2]" ne fonctionne pas.
Il s'agit d'une ligne que matlab inscrit lors de l'exécution.

Exite-t-il une autre ligne de commande qui permette de spécifier l'intervalle dans lequel le zéro doit être trouvé ?

Merci d'avance.

non

il faut pour cela changer de solveur (et du coup redéfinir le problème)
soit fminbnd soit fmincon

Fabien

edit : sinon tu peux truander, en donnant une valeur fixe non nulle à ta fonction en dehors de l'interval souhaité (dans la définition de ta fonction)