Discussion :

[sky-mars]

Bonjour tout le monde
J'ai réalisé des hélices sous matlab avec 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
20
21
22
23
24
25
26
% Modelisation helice
clear all;
X0=0:2:220%linspace(0,6.1544,6.1544/499);
A=15;
Y0=A*sin(X0/4)+20*sin(X0/8)+10*sin(X0)
% Formule de la rotation de centre O et d'angle alpha
alpha=42*pi/180;
X=X0*cos(alpha)-Y0*sin(alpha);
Y=Y0*cos(alpha)+X0*sin(alpha);
P=[X',Y'];
save fodil_p.txt P -ASCII
R=50;
teta=0;
Xs(1)=R;
Zs(1)=0;
for i=1:(length(X)-1)
    teta=1/R*(X(i+1)-X(i))+teta;
    Xs(i+1)=R*cos(teta);
    Zs(i+1)=R*sin(teta);
end;
%plot(X,Y)
plot3(Xs(1,:),Zs(1,:),Y)
M=[Xs(1,:)',Zs(1,:)',Y']
length (M)
save fodil.txt M -ASCII
grid on

et en fait je souhaiterai déterminer la longueur, donc ce que je fais c'est le calcul d'une norme euclidienne mais je suis pas sur que ma modélisation sous matlab soit OK , pourrai-t-on me le confirmer s'il vous plait merci

Norme

Code :

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

1
2
3
4
L=0
for j=1:length(X0)-1
L=L+sqrt((X0(j+1)-X0(j))^2+(Y0(j+1)-Y0(j))^2)
end;

[nahouto]

Ca m'a l'air d'être pas mal

[Caro-Line]

Et sinon il y a la fonction NORM

Code :

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

help norm

[sky-mars]

Voici les résultats
1) avec ma formule
2) norm appliqué sur la matrice P

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
 
>> L
 
L =
 
  1.2913e+003
 
>> norm(P)
 
ans =
 
  1.3413e+003

[Newenda]

A vérifier qu'en faisant ceci :

Code :

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

for j=1:length(X0)-1

tu ne prives pas ton calcul de la dernière longueur..ce qui pourrait être à l'origine de la légère différence des deux valeurs.

[sky-mars]

Salut
Ouais mais voila X0 est de taille 111 si j'enlève le -1 beh X0(112) ne sera pas définit

[Jerome Briot]

Code :

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

1
2
3
4
L=0
for j=1:length(X0)-1
L=L+sqrt((X0(j+1)-X0(j))^2+(Y0(j+1)-Y0(j))^2)
end;

Pas besoin de la boucle FOR-END :

Code :

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

L = sqrt(sum((X0(2:end)-X0(1:end-1)).^2+(Y0(2:end)-Y0(1:end-1)).^2))