Discussion :

[xe4b4ct]

Bonjour à tous,

Je vous écrits car j'essaye de résoudre un problème depuis si longtemps que j'en ai perdu mes cheveux.

voici un code sous matlab pour fitter une solutions à une équation de diffusion simple :

Code :

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

1
2
3
4
5
6
7
c=@(d,xdata)(10^6*m0./(d(4)*sqrt(4*pi*d(2).*xdata)).*exp(-((d(1)-d(3).*xdata).^2)./(4*d(2).*xdata)));
 
    xdata=[1:dt:length(J)-1];
    ydata=J(2:length(J),4)';
 
    option=optimset('TolFun',1e-15,'TolX',1e-15,'MaxIter',10000);
    [d,resnorm]=lsqcurvefit(c,d0,xdata,ydata,[x0-2 0.95*k0 v0-0.1 0.95*s0],[x0+2 0.95*k0 v0+0.1 1.05*s0],option);

et voici le denrier message d'erreur que matlab m'affiche sans se lasser.

??? Error using ==> optim\private\snls
Equal upper and lower bounds not permitted.

Error in ==> optim\private\lsqncommon at 213
[x,FVAL,LAMBDA,JACOB,EXITFLAG,OUTPUT,msg]=...

Error in ==> lsqcurvefit at 165
[x,Resnorm,FVAL,EXITFLAG,OUTPUT,LAMBDA,JACOB] = ...

Quelqu'un peut m'aider?

Merci par avance

[Caro-Line]

Je n'ai pas l'Optimization Toolbox et ne suis pas spécialiste mais si je traduis juste le message d'erreur en anglais il te dit qu'il n'accepte pas que tes bornes inférieures et supérieurs soient égales.
En regardant l'aide je pense qu'il doit s'agir de ces 2 paramètres :

Code :

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

[x0-2 0.95*k0 v0-0.1 0.95*s0],[x0+2 0.95*k0 v0+0.1 1.05*s0]

qu'il faudrait donc que tu revois.
Apparemment les 2èmes termes de chaque vecteur sont égaux, je suppose que c'est de là que vient l'erreur.

[xe4b4ct]

Merci à toi Caro_line,

voilà le genre d'erreur toute bête qui m'énerve au plus au point.

J'aurai du suivre les cours d'anglais,

Encore merci,

xe4b4ct