Discussion :

[Invité]

Représentation d'un ellipsoïde Matlab
Bonjour,

je souhaiterai représenter un ellipsoïde et le déformer par la suite pour lui donner une forme particulière .
je pense utilisé la fonction pré-enregistrer de Matlab ellipsoid() mais je dois d'abord crée les coefficients k3 et K2 sachant que comme ceci varient je pensais les stocker dans deux matrice vides
c'est coefficient doivent multiplier les valeurs du demi petit axe et demi grand axe de mon ellipsoide
je ne sais pas si c'est une bonne idée
j'ai pour le moment effectuer le code suivant:

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
 
l = 1;
x = 0:0.01:1; 
N = length(x);
k3 = [];
k2 = [];
for  i = 0:N
    xtild = (x(i)-0.2*l)/(0.75*l-0.2*l);
    k3 = -(1/6)*(sin(pi*x(i)-pi/2)+1)-(x(i)^2)/8 +1/2;
    if (x(i)>=0 && x(i)<0.2)
         k2 = 1;
    elseif  (x(gi)>=0.75 && x(i)< 1)
         k2 = 1-(sin(pi*xtild-pi/2)+1)/5; 
    else(x(i)>=0.75 && x(i)<1); 
         k2 = 3/5; 
end 
     k3 = a;
     k2 = b;
end

Du coup je voulais savoir si quelqu'un voyait une erreur flagrante sachant que j'ai l'erreur suivant qui s'affiche:

Code :

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

1
2
3
4
5
 
Array indices must be positive integers or logical values.
 
Error in Untitled4 (line 7)
    xtild = (x(i)-0.2*l)/(0.75*l-0.2*l);

Bien cordialement,

Marsu23

[blablu]

Pour ton erreur ligne 7. Dans Matlab, les indices commencent à 1 et pas à 0. Du coup c'est pour i dans [1,N].

[Invité]

merci pour la réponse

[Invité]

Du coup j'ai procédé de la manière suivante

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
l = 1;
x = 0:0.01:1; 
N = length(x);
xtild = zeros(1,N);
k3 = zeros(1,N);
k2 = zeros(1,N)
for  i = 1:N
    xtild(1,i) =(x(i)-0.2*l)/(0.75*l-0.2*l);
    k3(1,i) = -(1/6)*(sin(pi*x(i)-pi/2)+1)-(x(i)^2)/8 +1/2; 
    if (x(i)>=0 && x(i)<0.2)
         k2(1,i) = 1;
    elseif (x(i)>=0.75 && x(i)< 1)
          k2(1,i) = 1-(sin(pi*xtild(1,i)-pi/2)+1)/5; 
    else(x(i)>=0.75 && x(i)<1); 
       k2(1,i) = 3/5; 
    end
    disp(k2);
    disp(k3);
end