Approximation polynomiale d'un nuage de points

Bonjour à tous

J'ai un nuage de point pour lequel je veux calculer l'équation de la courbe passant par tous les points.

J'ai utilisé la fonction polyfit pour me calculer les coefficients du polynôme. J'ai utilisé un degré n - 1 ou n est mon nombre de points.

jusque là pas de souci et avec la fonction polyval j'ai construit la courbe que j'ai affiché en même temps que mes points pour vérifier.

Maintenant ce que j'aimerais faire, c'est l'inverse. c'est à dire obtenir la valeur de x pour une valeur y donnée. (polyval elle me donnant y en fonction de x).

Est-ce possible ?

Merci d'avance pour votre aide

Il "suffit" de calculer les racines du polynôme P(x) - y = 0 avec la fonction ROOTS, non ?

Justement, je ne sais pas ?

si p est le polynome obtenu avec polyfit alors le code suivant te donne le polynome dont t'a parlé dut dans son post.

Code:

---

p(end) = p(end)-y ;

---

Pour le résoudre :

Code:

---

1 2 3 4 5 6

p(end) = p(end)-y ; c = root(p) ;   % plot hold on plot(c, repmat(y, size(c)), 'o')

---