Discussion :

[ocelote]

Bonjour,

j'ai suivi ce tutoriel pour déterminer un plan a partir de plusieurs points.

J'ai essayé de suivre au mieux les indications, mais mon resultat est peu probant :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
 
    X = Pts(:,1)
    Y = Pts(:,2)
    Z = Pts(:,3)
    U = ones(length(X), 1);
    M = [U X Y X.^2 X.*Y Y.^2]
    K = M\Z;
    Z;
    Z_modele = M*K;
    residus = Z-Z_modele; 
    variance_d_origine = var(Z);
    variance_expliquee = var(Z_modele);
    variance_residuelle = var(residus);
    coefficient_correlation_r2 = variance_expliquee/variance_d_origine

Ou X contient tous les X de mes points (x1,x2,...xn), Y et Z egalement sur le meme principe.

Je souhaiterais dessiner mon plan calculé mais je ne sais pas comment m'y prendre.

Merci

[davcha]

Une autre façon de faire consiste à calculer les valeurs singulières de la matrice formée par tes points, centrés en zéro.

En gros :

Code :

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

1
2
3
4
5
6
7
M=Pts - mean(Pts)
[u s v]=svd(M)
n=v(:,3) % normale à l'hyperplan.
n/=n(3)
[X Y]=meshgrid(-10:10)
Z=-n(1)*X-n(2)*Y
mesh(X,Y,Z)

[ocelote]

Merci pour ton aide, mais il y a un probleme a la ligne : n/=n(3) et je n'ai pas suffisamment compris pour corriger moi meme

[davcha]

n/=n(3) c'est équivalent à n=n/n(3).

[ocelote]

je vois, mais j'ai une erreur de dimension en utilisant mean(). Je dois soustraire la valeur moyenne de X, Y et Z pour chaque point ?

Par exemple : X(1) - moy(X) ?

[davcha]

oui.