Interpolation log - log

Bonjour,

je cherche à faire une interpolation entre deux points représentés dans un espace log-log (échelle logarithme sur les deux axes). Ok, pas de problème je cherche dans mon grimoire d'algo et je ressorts cette formule :

Code:

---

1 2	y = 10^{ log(y0)*log(x1/x) + log(y1)*log(x/x0) / log(x1/x0) }

---

avec (x0, y0) le premier point, (x1, y1) le second point et (x, y) le point à interpoler. Le problème c'est que ça ne marche pas, exemple avec :

Code:

---

1 2 3

x0 = 1.0 x1 = 100.0 y0 = 100.0 y1 = 2.0

---

Si je choisis x=10, j'obtiens y=141.4!!!!, la valeur devrait être y1<y<y0.

Quelqu'un à des infos sur l'interpolation log-log ?

Merci

Bon, j'ai résolu mon problème, avec mon petit crayon et ma petite feuille j'ai retrouvé l'équation originale (une erreur sur une parenthèse). Si ça peut aider d'autres personnes voila l'équation :

Code:

---

1 2	y = 10^{ ( log(y0)*log(x1/x) + log(y1)*log(x/x0) ) / log(x1/x0) }

---

Ce calcul à besoin de deux *, quatre / et d'une +, le même résultat peut être obtenu avec une multiplication en moins :

Code:

---

1 2	y = 10^ { log(y0) + log(y1/y0) * log(x/x0) / log(x1/x0) }

---

et cette fois-ci le résultat est correct y=14.1.

A+