Image:Abaque wulff.png

	This Math image should be recreated using vector graphics as an SVG file. This has several advantages; see Commons:Media for cleanup for more information. If an SVG form of this image is already available, please upload it. After uploading an SVG, replace this template with template {{Vector version available|new image name.svg}} in this image.
العربية | Български | Català | Česky | Dansk | Deutsch | English | Esperanto | Español | Français | 한국어 | Italiano | Magyar | Lietuvių | Nederlands | 日本語 | Polski | Português | Română | Русский | Suomi | Svenska | Türkçe | Українська | Tiếng Việt | मराठी | ‪中文(繁體)‬ | ‪中文(简体)‬ | +/-

[edit] Summary

Abaque de Wulff (pas de 10°) pour la lecture de figure de pôles

Wulff net (step of 10°) used to read pole figures

Voir aussi/see also Image:Abaque wulff 2 degres.pdf, Image:Abaque wulff simplifie.png

Auteur/author : Christophe Dang Ngoc Chan (cdang)

Réalisé avec/made with : Scilab

clear;clf;
pas=0.01;

// Cercle 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

for j=1:8
   r1=cos(j*%pi/18); // décalage de 10°
   r2=sin(j*%pi/18);
   for i=1:ny
      Xm(i,j)=x(i)*r1; // coordonnées des points du cercle incliné
      Ym(i,j)=y(i);
      Zm(i,j)=x(i)*r2;
   end
end

// Cercles parallèles

for j=1:8
   r1=cos(j*%pi/18); // décalage de 10°
   r2=sin(j*%pi/18);
   for i=1:ny
      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:8
   for i=1:ny
      m = 1/(1+Zm(i,j)); // z=0
      Xm1(i,j)=m*Xm(i,j);
      Ym1(i,j)=m*Ym(i,j);
   end
end

for j=1:8
   for i=1:ny
      m = 1/(1+Zp(i,j)); // z=0
      Xp1(i,j)=m*Xp(i,j);
      Yp1(i,j)=m*Yp(i,j);
   end
end

// Tracé

angle=['10','20','30','40','50','60','70','80'];

for j=1:8
   plot2d(Xm1(:,j),Ym1(:,j),style=2)
   plot2d(-Xm1(:,j),Ym1(:,j),style=2)
   plot2d(Xm1(:,j),-Ym1(:,j),style=2)
   plot2d(-Xm1(:,j),-Ym1(:,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

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

xpoly([0,0],[-1,1])
xpoly([-1,1],[0,0])

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

[edit] Licensing

I, the copyright holder of this work, have published or hereby publish it under the following license: Permission is granted to copy, distribute and/or modify this document under the terms of the GNU Free Documentation license, Version 1.2 or any later version published by the Free Software Foundation; with no Invariant Sections, no Front-Cover Texts, and no Back-Cover Texts. A copy of the license is included in the section entitled "GNU Free Documentation license". Aragonés | العربية | Asturianu | Български | বাংলা | ইমার ঠার/বিষ্ণুপ্রিয়া মণিপুরী | Brezhoneg | Bosanski | Català | Cebuano | Česky | Dansk | Deutsch | Ελληνικά | English | Esperanto | Español | Eesti | Euskara | فارسی | Suomi | Français | Gaeilge | Galego | עברית | Hrvatski | Magyar | Bahasa Indonesia | Ido | Íslenska | Italiano | 日本語 | ქართული | ភាសាខ្មែរ | 한국어 | Kurdî / كوردی | Latina | Lëtzebuergesch | Lietuvių | Bahasa Melayu | Nnapulitano | Nederlands | ‪Norsk (nynorsk)‬ | ‪Norsk (bokmål)‬ | Occitan | Polski | Português | Română | Русский | Slovenčina | Slovenščina | Shqip | Српски / Srpski | Svenska | తెలుగు | ไทย | Türkçe | Українська | اردو | Tiếng Việt | Volapük | Yorùbá | ‪中文(中国大陆)‬ | ‪中文(台灣)‬ | +/-

La page ci-dessous contient cette image :

• Figure de pôles