Discussion :

[hollowdeadoss]

Bonjour tout le monde,

J'ai un souci à comprendre le code suivant:

Code :

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

[x,resnorm] = lsqcurvefit(@(x,xdata) myfunction(x,4),x0,[],ydata,[-Inf,-Inf,0,0],[],options_lsqcurvefit) ;

Ce qui me gêne c'est surtout ce @(x,xdata) suivi de "myfuntcion" où xdata n'est pas appelé.
De plus un peu plus loin on voit que xdata est égal à [].
Du coup j'avoue ne pas saisir le fonctionnement.
Enfin je comprends le principe de cette fonction (moindre carré pour les problèmes non linéaire) mais là vraiment je bloque sur ce @(x,xdata) ...

Si quelqu'un a une idée il ne faut pas hésiter

(j'essaye de me servir de ce code pour un problème vraiment similaire!).

Cordialement,

Hollow

[Jean Dumoncel]

Bonjour,

Le @(x,xdata) suivi de "myfuntcion" permet de précsiser que lsqcurvefit s'applique sur la fonction myfuntcion avec comme variable d'entrées : x et xdata.

Par contre, je trouve étonnant la façon d'appeler myfonction :

Code :

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

myfunction(x,4)

Regarde les exemples de la doc, il faut plutôt mettre :

Code :

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

myfunction(x,xdata,4)

(et du coup xdata doit apparaître dans la fonction myfonction)

De plus comment l'algorithme peut s'exécuter si xdata est un vecteur vide?? c'est sur xdata et ydata que s'effectue l'identification de x... tu l'as testé ce code? Il donne un résultat?

Je te conseille de bien regarder les exemples de la doc, si tu as des questions, n'hésite pas.

[tesca]

Le xdata peut etre nul (pas sûr; si tu l'utilise pas, il faut qu'il soit au moins de même longueur que ydata). En fait lsqcurvefit fait une approximation des moindres carré entre ydata et la sortie de myfunction. Le seul parametre vraiment obligatoire que doit prendre myfunction est le x.

myfunction est un truc du genre

Code :

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

function y=myfonction(x,xdata,autre variable)

"autre variable n'est possible que si tu passe par le @(x,xdata)"

par exemple tu cherche les paramétre de tes points ydata en Ae-t/tau

Code :

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

y=x(1)*exp(-xdata/x(2))

Mais si tu veut faire des "fitting originale" tu n'a pas besoin du xdata, il suffit que y soit de meme longueur que ydata pour faire l'approximation

[Jean Dumoncel]

Oui effectivement, j'ai écrit un peu vite.

[hollowdeadoss]

Merci bien!

En fait dans la fonction "myfunction" les arguments sont justes les paramètres "x" et l'ordre d'un polynôme (ici égal à 4).
La fonction n'est pas de moi, mais j'ai testé et ça tourne sans problème.

Mais j'avoue que j'y comprenais pas grand chose avec ces xdata en vecteur vide, et le myfunction où le xdata n'intervient pas.

Merci pour votre aide,

Hollow