Régression sphérique

[robinho81]

bonjour
je travaille sur une méthode de régression sphérique qui peut donner les coordonnes de la sphère en fonction des covariants selon la méthode:
http://www.cs.bsu.edu/homepages/kerryj/kjones/circles.pdf
j'ai pu coder la méthode en C++ mais je trouve toujours des écarts par rapport à la théorique
je pense que j'ai fait des erreurs sur l'estimation des covariants sxx2,sxy2,sxz2 utilisé dans l'expression analytique des coordonnées:
ax=((Sxx2 )+Sxy2 )+Sxz2 ) )(Sy2 Sz2-Syz2 )+(Syx2 )+Syy2 )+Syz2 ) )(Sxz Syz-Sxz Sz^2 )+(Szx2 ....
est ce qu'il y'a quelqu'un qui peut m'aider

[Cliffe]

Faut appliquer les formules non ?

[CENTER] [/CENTER]

En C++ :

Code: Tout sélectionner

// n = ...
// Calculer m_x = moyenne(x)
// Calculer m_y = moyenne(y)
// Calculer m_z = moyenne(z)
// Calculer les produits

double S(std::vector & v1, std::vector & v2, double m_v) {
    double sum = 0.;
    for (int i = 0; i < n; ++i)
        sum += v1[i] * v2[i]
    return sum / n - m_v;
}


[JeanJ]

Pour information : une méthode de régression sphérique facile à programmer est exposée dans l'article : "Régressions coniques, quadriques, circulaire, sphérique,..." (sphérique p. 17-18), par le lien :
http://www.scribd.com/JJacquelin/documents

[robinho81]

Faut appliquer les formules non ?

[CENTER] [/CENTER]

En C++ :

Code: Tout sélectionner

// n = ...
// Calculer m_x = moyenne(x)
// Calculer m_y = moyenne(y)
// Calculer m_z = moyenne(z)
// Calculer les produits

double S(std::vector & v1, std::vector & v2, double m_v) {
    double sum = 0.;
    for (int i = 0; i < n; ++i)
        sum += v1[i] * v2[i]
    return sum / n - m_v;
}


merci mais je veux connaitre l'expression de Sxx2

[Cliffe]

Tu prend .

[robinho81]

mais je veux aussi connaitre les expressions Sxx2, Sxy2

[robinho81]

sxx2=1/n somme(x^3)-moy(x)^3

[JeanJ]

Salut robinho81
Tous ces termes sont écrits explicitement page 18 du document déjà cité :
http://fr.scribd.com/doc/14819165/Regressions-coniques-quadriques-circulaire-spherique