Discussion :

[membreComplexe12]

Salut tous,

j'utilise actuellement matlab pour essayer de fitter des courbes "brut" (qui sont donc un peu bruité) afin de pouvoir associer une fonction mathématique à ma courbe et afin d'éliminer le bruit grace à ma courbe mathématique.

Voici en PJ mes courbes et les données associées à la première courbe.

j'ai fais pas mal de tentatives (matlab, fonction polyfit) avec des modèles exponentiels/logarithmiques...
mais ça ne fonctionne pas tout à fait comme je voudrais

Je viens vous voir car les personnes du forum matlab m'ont renvoyées vers vous pour avoir des informations complémentaires :
http://www.developpez.net/forums/d12...ifiee-polyfit/

ps: j'aimerai avoir une fonction monotone de préférence (sinon j'utiliserai un polynome et ça marcherai bien)

je vous remercie d'avance pour votre aide

[Alexis.M]

Si tu veux modéliser une courbe, c'est bien de savoir si justement il y a un modèle utilisable derrière. Donc première question:

D'où viennent ces données, quel phénomène les produit ?

Sinon à vue de nez la fonction:

a * log(1+exp(b*(x-x0))) + c

fait bien l'affaire pour la partie linéaire en tout cas ensuite ça colle moins bien( voir la pièce jointe)

D'où vient cette formule ? log(1+exp(b*x))∕b est parfois appelé la fonction de cout logistique et c'est une approximation lisse de max(0,x).
Donc en ajoutant quelques coefficients pour gérer l'échelle le décalage horizontal et le décalage vertical on obtient la formule donnée au début.

[Alexis.M]

Tu peux faire mieux:

a*log(1+exp(b*(x-x0)) + c + d*x

avec ça j'obtiens un fit presque parfait. Tu "fittes" d'abord sans le terme d*x
(la pente résiduelle etant faible) et tu réutilise les paramètre trouvés pour fitter la formule avec le d*x.

Explication des paramètres:

x0: donne la localisation du "coude"
b: donne la le caractère lisse ou non de la transition au niveau du coude.
a*b: donne la contribution à la pente loin du coude
c: le décalage vertical
d: la pente résiduelle.

donc dans la partie la plus pentue asymptotiquement la pente est a*b*d
est dans la partie résiduelle la pente est d.

Avec cela en tête il est facile d'estimer les paramètres de départ.

[membreComplexe12]

Tu peux faire mieux:

a*log(1+exp(b*(x-x0)) + c + d*x

avec ça j'obtiens un fit presque parfait. Tu "fittes" d'abord sans le terme d*x
(la pente résiduelle etant faible) et tu réutilise les paramètre trouvés pour fitter la formule avec le d*x.

Explication des paramètres:

x0: donne la localisation du "coude"
b: donne la le caractère lisse ou non de la transition au niveau du coude.
a*b: donne la contribution à la pente loin du coude
c: le décalage vertical
d: la pente résiduelle.

donc dans la partie la plus pentue asymptotiquement la pente est a*b*d
est dans la partie résiduelle la pente est d.

Avec cela en tête il est facile d'estimer les paramètres de départ.

ouahou

c nikel ça
tu as effectué ceci sur matlab ?

es ce possible de me montrer ton bout de code ?

[membreComplexe12]

merci beaucoup Alexis d'avoir pris le temps de répondre.

Si tu veux modéliser une courbe, c'est bien de savoir si justement il y a un modèle utilisable derrière. Donc première question:
D'où viennent ces données, quel phénomène les produit ?

oui je vois ce que tu veux dire. Ici c'est un cas un peu complexe car justement on a pas de modèle connu qui permet de réprésenter ceci physiquement. (perte de non linéarité d'un ressort dû à une variation interne de la structure).
=> du coup je veux avoir un modèle simple mais et voir comment évolue les coefficient identifiés en fonction de la variation interne de la structure
=> je veux qqch de monotone car dans la réalité de phénomène est monotone

Sinon à vue de nez la fonction:
a * log(1+exp(b*(x-x0))) + c
fait bien l'affaire pour la partie linéaire en tout cas ensuite ça colle moins bien( voir la pièce jointe)

tu as oublié la PJ

Sinon je pense qu'une relation

Code :

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

1-exp(-X/A)

ou
doit coller pas

Code :

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

1-exp(-A/X)

mal mais je ne sais pas trop comment les identifier avec polyfit de matlab....

D'où vient cette formule ? log(1+exp(b*x))∕b est parfois appelé la fonction de cout logistique et c'est une approximation lisse de max(0,x).
Donc en ajoutant quelques coefficients pour gérer l'échelle le décalage horizontal et le décalage vertical on obtient la formule donnée au début.

d'acccord, je comprends le principe, je regarderai.
merci

[Alexis.M]

Attention nous avons eu des collisions de messages

Non j'ai fais ça avec python/scipy/matplotlib

[membreComplexe12]

Attention nous avons eu des collisions de messages

oui, en effet

Non j'ai fais ça avec python/scipy/matplotlib

ah OK. Peux tu me montrer quand meme le bout de code ?

=> quand tu fais l'identification c'est semi automatiquement alors ?

=> tu identifie sans le terme d*x de manière automatique et ensuite tu identifie avec le terme d*x en donnant "à la main" la position de la cassure ?

[Alexis.M]

si cela t'intéresse:

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
 
import numpy as np
from scipy.optimize import leastsq
 
x,y = np.load("data.txt").T
 
def myfunc(x,x0,a,b,c):
    return b* np.log(1+np.exp(a*(x-x0))) + c
 
def myfunc2(x,x0,a,b,c,d):
    return b* np.log(1+np.exp(a*(x-x0))) + c + d*x
 
def err(params):
    return myfunc(x,*params) - y
 
def err2(params):
    return myfunc2(x,*params) - y
 
x0, a, b, c  = leastsq(err, [0.004, 1000, 10.,300.0])[0]
x0, a, b, c, d = leastsq(err2, [x0,a,b,c,0.0])[0]

[Alexis.M]

L'emplacement du coude est déjà estimer au premier tour, donc tu peux juste mettre la pente à 0.0 pour le deuxième tour.

[membreComplexe12]

merci Alexis.

=> je vais regarder merci beaucoup pour ton aide.


ps: j'ai trouvé aussi d'autre choses :
http://fr.wikipedia.org/wiki/%C3%89q...chaelis-Menten
et
je vais regarder du coté des exponentielles ssaturantes (1-exp) exp(A/x) ...