Discussion :

[conanran]

Bonjour,
J'ai des données (X,Y), j'ai tracé la courbe Y=f(X) : (voir courbe rouge de la figure).
Je voudrai savoir comment je peux trouver la forme de la fonction à partir de cette courbe. j'ai essayé d'utiliser la fonction fil mais j'ai pas réussi de trouver le modèle général ainsi que les bonne entrée de la fonction fit dans mon cas.

Code :

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

1
2
3
f=fit(X,Y,'smoothingspline')
 
plot(f,X,Y)

Merci pour votre aide.

[thecrazydonut]

Bonjour,

f est un objet contenant tous les informations concernant le fit dans le champ p.

Code :

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

1
2
3
4
5
6
7
8
f.p
ans = 
      form: 'pp'
    breaks: [1x168 double]
     coefs: [167x4 double]
    pieces: 167
     order: 4
       dim: 1

En revanche, tu peux voir que la forme n'est pas triviale puisqu'il s'agit d'un fit polynomial par partie (piecewise polynomial fitting). Il y a donc autant d’équations décrivant ta courbe que de segments. "L'équation générale" est donc un ensemble d'équations. Si tu veux une équation toute belle et simple de type polynomiale, il faut abandonner l'argument 'smoothingspline' et effectuer un fit plus simple, par exemple avec 'poly3'.

[conanran]

Merci pour votre réponse
J'ai utilisé 'poly5' et il a l'air pas mal.
En effet, en jouant sur un autre paramètre je trouve presque la même allure avec un petit décalage (courbe bleu). Je voudrais savoir si possible de fixer les mêmes coefs de l'équation générale trouvés pour la courbe rouge et identifier que le coef qui joue sur la position de la courbe.
Merci

[thecrazydonut]

En effet, en jouant sur un autre paramètre je trouve presque la même allure avec un petit décalage (courbe bleu). Je voudrais savoir si possible de fixer les mêmes coefs de l'équation générale trouvés pour la courbe rouge et identifier que le coef qui joue sur la position de la courbe.

J'ai pas tout compris. Que représente la courbe bleue ? C'est certain que ce n'est pas le fit avec poly5....
Deuxième phrase, ce que je comprends c'est que tu veux décaler la courbe obtenue avec le fit horizontalement (fit avec les 4 premiers paramètres déjà déterminés, seul le 5e étant libre). Je ne crois pas qu'il existe une option fit avec contrainte dans les fonctions de base de matlab de cette forme, mais elle peut se coder facilement (avec fminsearch par exemple). Mais bon, peut-être que pour tes besoins tu peux juste faire une faire ça manuellement avec appréciation visuelle.

[SpaceX]

Salut,

Quand tu visualise ta courbe, tu peux utiliser la fonction intégré à Matlab : "Basic Fitting" dans l'onglet Tools.

Ensuite tu coche "Show Equations" pour visualiser l'équation de la courbe qu'il a réussi a approximé.