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;