Ce fichier était déposé au format PDF, alors que son contenu se résumait à une simple image. Si vous le pouvez, merci de remplacer le ou les fichiers concernés par l'image elle-même (éventuellement extraite du fichier original). Veuillez consulter la liste des formats de fichiers pour plus d'informations. Deutsch ∙ English ∙ español ∙ français ∙ हिन्दी ∙ 日本語 ∙ македонски ∙ português ∙ русский ∙ slovenščina ∙ 中文 ∙ +/−

Abaque de Wulff avec un pas de 2 °, pour la lecture de figure de pôles

Wulff net with a step of 2°, used to read pole figures

Voir aussi/see also : image:abaque wulff.png

Auteur/author : Christophe Dang Ngoc Chan (cdang)

Réalisé avec/made with : Scilab

clear;clf;
pas=0.005;

// Cercle horizontal de base

y=[0:pas:1]';
NY=size(y); ny=NY(1,1);

for i=1:ny
   x(i)=sqrt(1-y(i)^2);
end

// Cercles méridiens et parallèles

for j=1:44
   r1=cos(j*%pi/90); // décalage de 2°
   r2=sin(j*%pi/90);
   for i=1:ny
// Cercles méridiens
      Xm(i,j)=x(i)*r1; // coordonnées des points du cercle incliné
      Ym(i,j)=y(i);
      Zm(i,j)=x(i)*r2;
// Cercles parallèles
      Xp(i,j)=y(i)*r2; // coordonnées des points du cercle parallèle
      Yp(i,j)=r1;
      Zp(i,j)=x(i)*r2;
   end
end


// Projection stéréographique

// Pôle Sud S(0,0,-1), pôle P(X1,Y,Z)
// Droite : S + m*SP = (m*X1)i + (m*Y)j + (-1+m*(Z+1))k

for j=1:44
   for i=1:ny
// Cercles méridiens
      mm = 1/(1+Zm(i,j)); // z=0
      Xm1(i,j)=mm*Xm(i,j);
      Ym1(i,j)=mm*Ym(i,j);
// Cercles parallèles
      mp = 1/(1+Zp(i,j)); // z=0
      Xp1(i,j)=mp*Xp(i,j);
      Yp1(i,j)=mp*Yp(i,j);
   end
end

// Tracé

// Méridiens : simplifié près des extrémités (10°)

k=int((1-cos(%pi/18))/pas);
l=int((1-cos(%pi/90))/pas);
if l=0 then l=1; end

// Graduations par 2°

for j=[[1:4],[6:9],[11:14],[16:19],[21:24],[26:29],[31:34],[16:39],[41:44]]

   plot2d(Xm1(1:ny-k,j),Ym1(1:ny-k,j),style=2)
   plot2d(-Xm1(1:ny-k,j),Ym1(1:ny-k,j),style=2)
   plot2d(Xm1(1:ny-k,j),-Ym1(1:ny-k,j),style=2)
   plot2d(-Xm1(1:ny-k,j),-Ym1(1:ny-k,j),style=2)
   
   plot2d(Xp1(:,j),Yp1(:,j),style=2)
   plot2d(-Xp1(:,j),Yp1(:,j),style=2)
   plot2d(Xp1(:,j),-Yp1(:,j),style=2)
   plot2d(-Xp1(:,j),-Yp1(:,j),style=2)

end

// Graduations principales

for j=5:5:44

   xpoly(Xm1(1:ny-k,j),Ym1(1:ny-k,j))
   a=get('hdl');a.thickness=2;
   xpoly(-Xm1(1:ny-k,j),Ym1(1:ny-k,j))
   a=get('hdl');a.thickness=2;
   xpoly(Xm1(1:ny-k,j),-Ym1(1:ny-k,j))
   a=get('hdl');a.thickness=2;
   xpoly(-Xm1(1:ny-k,j),-Ym1(1:ny-k,j))
   a=get('hdl');a.thickness=2;
   plot2d(Xm1(ny-k:ny-l,j),Ym1(ny-k:ny-l,j))
   plot2d(-Xm1(ny-k:ny-l,j),Ym1(ny-k:ny-l,j))
   plot2d(Xm1(ny-k:ny-l,j),-Ym1(ny-k:ny-l,j))
   plot2d(-Xm1(ny-k:ny-l,j),-Ym1(ny-k:ny-l,j))

   xpoly(Xp1(:,j),Yp1(:,j))
   a=get('hdl');a.thickness=2;
   xpoly(-Xp1(:,j),Yp1(:,j))
   a=get('hdl');a.thickness=2;
   xpoly(Xp1(:,j),-Yp1(:,j))
   a=get('hdl');a.thickness=2;
   xpoly(-Xp1(:,j),-Yp1(:,j))
   a=get('hdl');a.thickness=2;

end

// axes et cercle autour

plot2d(x,y)
plot2d(-x,y)
plot2d(x,-y)
plot2d(-x,-y)

xpoly([0,0],[-1,1])
a=get('hdl');a.thickness=2;
xpoly([-1,1],[0,0])
a=get('hdl');a.thickness=2;

axe=get('current_axes');
axe.axes_visible=['off','off'];
axe.isoview='on';
axe.box='off';

Moi, en tant que détenteur des droits d’auteur sur cette œuvre, je la publie sous les licences suivantes :

Vous avez la permission de copier, distribuer et modifier ce document selon les termes de la GNU Free Documentation License version 1.2 ou toute version ultérieure publiée par la Free Software Foundation, sans sections inaltérables, sans texte de première page de couverture et sans texte de dernière page de couverture. Un exemplaire de la licence est inclus dans la section intitulée GNU Free Documentation License.http://www.gnu.org/copyleft/fdl.htmlGFDLGNU Free Documentation Licensetruetrue

	Ce fichier est disponible selon les termes de la licence Creative Commons Attribution – Partage dans les Mêmes Conditions 3.0 Non Transposé.
	Vous êtes libre : de partager – de copier, distribuer et transmettre cette œuvre d’adapter – de modifier cette œuvre Sous les conditions suivantes : paternité – Vous devez donner les informations appropriées concernant l'auteur, fournir un lien vers la licence et indiquer si des modifications ont été faites. Vous pouvez faire cela par tout moyen raisonnable, mais en aucune façon suggérant que l’auteur vous soutient ou approuve l’utilisation que vous en faites. partage à l’identique – Si vous modifiez, transformez ou vous basez sur cet élément, vous devez distribuer votre contribution sous une license identique ou compatible à celle de l’original.
	Ce bandeau de licence a été ajouté à ce fichier dans le cadre de la procédure de mise à jour des licences des images sous GFDL.http://creativecommons.org/licenses/by-sa/3.0/CC BY-SA 3.0Creative Commons Attribution-Share Alike 3.0truetrue

Vous pouvez choisir l’une de ces licences.