inverser une fonction

Bonjour tout le monde,

J'ai une fonction polynomiale dite "ep" de degré 9. Je cherche à l'inverser de la manière suivante:

Code:

---

ep_inverse = finverse(ep);

---

Cela fonctionne bien!

Ensuite, je veux calculer une valeur en utilisant "ep_inverse" et la ca explose... :aie:

Les erreurs sont les suivantes:

Citation:

---

??? Error using ==> sym.sym>checkindex at 2729
Index must be a positive integer or logical.

Error in ==> sym.sym>privformatscalar at 2679
checkindex(x);

Error in ==> sym.sym>privformat at 2663
s = privformatscalar(x);

Error in ==> sym.sym>sym.subsref at 1379
[inds{k},refs{k}] = privformat(inds{k});

Error in ==> Polynomium_v5 at 151
abs_point_cle_1 = ep_inverse(ord_point_cle_1);

---

Est-ce que quelqu'un pourrait m'aider à résoudre ce problème svp?

Merci d'avance! :ccool:

il ne faut pas utiliser la syntaxe g(5) mais la fonction subs ou matlabFunction
avec

Code:

---

1 2 3 4 5 6 7 8 9 10 11 12

>> syms x >> f = x^2 + 3*x +4   f =   x^2 + 3*x + 4   >> g = finverse(f)   g =   (4*x - 7)^(1/2)/2 - 3/2

---

tu peux faire

Code:

---

1 2 3 4 5

>> subs(g,x,5)   ans =   0.3028

---

ou

Code:

---

1 2 3 4 5 6 7 8 9 10 11 12 13

>> x=5   x =   5   >> subs(g)   ans =   0.3028   >>

---

ou encore

Code:

---

1 2 3 4 5 6 7 8 9 10 11

>> gg = matlabFunction(g)   gg =   @(x)sqrt(x.*4.0-7.0).*(1.0./2.0)-3.0./2.0   >> gg(5)   ans =   0.3028

---

voir aussi sym2poly puis polyval

Bonjour et merci de ta réponse le fab! :ccool:

Malheureusement, je pense que j'ai un problème pour appliquer tes conseils. Avec la syntaxe suivante, j'ai des problèmes car abs_point_cle_1 et abs_point_cle_2 ne sont pas représentés dans un tableur dans mon explorateur de variable mais par un "petit cube"...

Code:

---

1 2 3 4 5 6 7 8 9 10 11 12 13 14

ord_point_cle_1 = (E1+E2)/2; ord_point_cle_2 = (E2+E3)/2; if degre_polyfit == 7 ep = p(1)*x^6 + p(2)*x^5 + p(3)*x^4 + p(4)*x^3 + p(5)*x^2 + p(6)*x^1 + p(7); elseif degre_polyfit == 8 ep = p(1)*x^7 + p(2)*x^6 + p(3)*x^5 + p(4)*x^4 + p(5)*x^3 + p(6)*x^2 + p(7)*x^1 + p(8); end   ep_inverse = finverse(ep); syms x; abs_point_cle_1 = subs(ep_inverse, x, ord_point_cle_1); abs_point_cle_2 = subs(ep_inverse, x, ord_point_cle_2); sol(1, 1) = abs_point_cle_1; sol(2, 1) = abs_point_cle_2;

---

le petit cube signifie qu'il s'agit d'un objet symbolique
E1 et E2 sont des objets symboliques ?? et p ??
car là c'est comme si tu avais une autre variable que x dans f

sinon, plutôt que d'écrire f à la main (visiblement après un polyfit) utilise plutôt poly2sym

tu devrais également à mon avis déclarer ton objet symbolique x avant de déclarer f puis d'inverser f. d'ailleurs je suis étonne que tu n'ai pas une erreur en déclarant f fonction de x avant d'avoir déclarer x via syms !!!
tu n'aurais pas déjà déclarer x d'une autre façon ??

si tu t'en sort pas, donnes moi un exemple complet (avec tes données)

essayes ca :

Code:

---

1 2 3 4 5 6 7 8 9 10 11 12

ord_point_cle_1 = (E1+E2)/2; ord_point_cle_2 = (E2+E3)/2;   syms x; ep = poly2sym(p,x)   ep_inverse = finverse(ep);   abs_point_cle_1 = subs(ep_inverse, x, ord_point_cle_1); abs_point_cle_2 = subs(ep_inverse, x, ord_point_cle_2); sol(1, 1) = abs_point_cle_1; sol(2, 1) = abs_point_cle_2;

---