Discussion :

[Naceur_5]

Bonjour,
je veux avec ce petit programme sur MATLAB d'afficher l'enveloppe convexe en se basant sur le théorème de Caratheodory
avec le coefficient Lemda de Caratheodory lemda1=lemda2=lemda3=lemda4=0. 25

Code :

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

1
2
3
4
5
6
7
8
9
10
x1=-9.771, x2=0.04587, x3=0.009142, x4=-7.131
y1=-9.728, y2=-0.5722, y3=0.841, y4=-1.843
z1=-8.754, z2=-1.354, z3=0.6607, z4=0.3261
 
x=[0.25*x1+0.25*x2+0.25*x3+0.25*x4; -0.25*x1-0.25*x2-0.25*x3-0.25*x4]
y=[0.25*y1+0.25*y2+0.25*y3+0.25*y4; -0.25*y1-0.25*y2-0.25*y3-0.25*y4]
z=[0.25*z1+0.25*z2+0.25*z3+0.25*z4; -0.25*z1-0.25*z2-0.25*z3-0.25*z4]
X=[x y z]
K = convhulln(X);
trisurf(K,X(:,1),X(:,2),X(:,3) )

mais j'obtiens toujours l'erreur suivante:

??? Error using ==> cgprechecks at 47
Not enough unique points specified.

Error in ==> convhulln at 42
cgprechecks(x, nargin, cg_opt);

Error in ==> CONVEXCOMBINATION3 at 25
K = convhulln(X);

pourriez-vous m'aidez Svp

et Merci beaucoup

[FLB]

Salut,
tu ne définis que 2 points : convhulln ne peut pas trouver de face à une droite!
Rajoute au moins un 3ème point pour pouvoir trouver les faces formant l'enveloppe convexe de ton nuage de point.