Discussion :

[kurul1]

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

[Jerome Briot]

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

[kurul1]

Justement, je ne sais pas ?

[ol9245]

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 :

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

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

Pour le résoudre :

Code :

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

1
2
3
4
5
6
p(end) = p(end)-y ;
c = root(p) ;
 
% plot
hold on
plot(c, repmat(y, size(c)), 'o')