Re: Fit

Top Page

Reply to this message
Author: Patrick Dupre
Date:  
To: Edgar Bonet
CC: guilde
Subject: Re: Fit
> « complexes » comme les nombres complexes, ou « complexes » comme
> compliquées ?
>

Cela peut etre les deux. J'ai reussi a ne pas utiliser des nombres complexes
pour l'intant !
Je dois eventuellement faire des diagonalisations.

> > Nonlinear squares.
>
> Tu veux dire « nonlinear least-squares » (moindres carrés).
>
> > Je peux ustiliser des outils simples, genre gnuplot.


gnuplot funcionne tres bien pour une fonction simple
f(x) = Ampl * (1 - (cos (x-x0))**2 * (1 - epsil))
Il faudrait que je puisse faire une sous-routine. Peut-etre cela
est possible.
Je n'ai pas vraiment vu le moyen d'ecrire une procedure.
Si on pouvait appeler du c depuis gnuplot, cela m'irait aussi.
Mais cela n'a pas l'air possible alors que l'inverse et possible.

> Et il ne convient pas ?
>
> > Cependant GSL n'a pas l'air de fonctionner (<gsl/gsl_multifit_nlin.h>).
>
> Ça veut dire quoi « n'a pas l'air » ?
>
> > Je suis etonne de ne pas trouver tout cela dans GSL !
>
> Trouver quoi ? Moi j'ai l'impression qu'on y trouve des choses [1].

Oui, sur le papier c'est bien mais en fait je n'ai pas reussi
a faire tourner
gsl_multifit_fsolver * gsl_multifit_fsolver_alloc (const gsl_multifit_fsolver_type * T, size_t n, size_t p)

Ca plante a l'allocation memoire.
Je n'ai pas vu d'exemple.
Il y a peut-ete une mailing liste GSL ou je pourrais avoir des infos.
gsl_multifit_dffsolver_alloc semble fonctionner mais cela necessite de
calculer les derives : pas forcement possible

Comme suggerer par Jereome
scipy.optimize.curve_fit semble interesante, dommage cela n'a pas ete porte en
perl.
Il y a bien
perl-Algorithm-CurveFit
mais j'ai compris que l'appel de la fonction ne pouvait etre qu'une function simple
(comme sous gnuplot).

Il va falloir que j'abandonne perl pour python !




===========================================================================
 Patrick DUPRÉ                                 | | email: pdupre@???
 Laboratoire de Physico-Chimie de l'Atmosphère | |
 Université du Littoral-Côte d'Opale           | |
 Tel.  (33)-(0)3 28 23 76 12                   | | Fax: 03 28 65 82 44
 189A, avenue Maurice Schumann                 | | 59140 Dunkerque, France
===========================================================================




>
> « complexes » comme les nombres complexes, ou « complexes » comme
> compliquées ?
>
> > Nonlinear squares.
>
> Tu veux dire « nonlinear least-squares » (moindres carrés).
>
> > Je peux ustiliser des outils simples, genre gnuplot.
>
> Et il ne convient pas ?
>
> > Cependant GSL n'a pas l'air de fonctionner (<gsl/gsl_multifit_nlin.h>).
>
> Ça veut dire quoi « n'a pas l'air » ?
>
> > Je suis etonne de ne pas trouver tout cela dans GSL !
>
> Trouver quoi ? Moi j'ai l'impression qu'on y trouve des choses [1].
>
> À+,
>
> Edgar.
>
> [1] https://www.gnu.org/software/gsl/manual/html_node/Nonlinear-Least_002dSquares-Fitting.html
>
>