jeudi 29 septembre 2011

Anamorphoses cylindriques : partie 3

On connaît la célèbre image au format .eps (tiger.eps) et je ne pensais pas possible d'en faire l'anamorphose, vectorielle bien sûr. Et bien, grâce à Jean-Michel Sarlat c'est fait !!! Jean-Michel a mis au point un script permettant de convertir un fichier .eps en commandes PStricks et partant de là de lui appliquer les formules de l'anamorphose cylindrique. Bonne nouvelle, tous les types d'anamorphose à venir : coniques, pyramidale, sphérique pourront bénéficier de cette considérable avancée !

\begin{center}
\begin{pspicture}[showgrid](-6,-8)(6,4)
\psset{Xv=0,Yv=-10}
\input{tiger.tex}
\input{ACtiger.tex}
\multido{\n=-2.00+0.50}{9}{%
    \pnode(! \n\space -2.00){A}
    \pnode(! \n\space 2.00){B}
    \psline(A)(B)
    \pslineAC(A)(B)
    }
\multido{\N=-2.00+0.50}{9}{%
    \pnode(!-2.0 \N){A}
    \pnode(!2.0 \N){B}
    \pslineAC[linestyle=solid](A)(B)
    \psline(A)(B)
    }
\pscircle[doubleline=true,linecolor=black]{3}
\end{pspicture}
\end{center}

La vérification expérimentale réalisée par Juergen Gilg :
montre que le codage n'est pas encore parfait !
Les fichiers sont ici :

Anamorphose cylindrique 3D

Il est facile d'établir les formules pour l'anamorphose en 3 D. L'anamorphose en 3D consiste à obtenir une image (virtuelle) en 3D , à partir d'une image(qui joue le rôle d'objet pour le miroir) plate.
C'est le package pst-solides3d avec l'option [transform=...] qui est utilisé.


 
 Les fichiers sont dans ce répertoire :

mardi 27 septembre 2011

Anamorphose cylindrique : partie 2

La construction d'une anamorphose cylindrique se présente maintenant sous forme d'un package : pst-anamorphosis.
Les commandes spécifiques sont calquées sur celles de PStricks suivies de AC, exemples : \pscircleAC, \pscurveAC, \psccurveAC, \psbezierAC etc. on peut aussi créer l'anamorphose d'un texte avec la macro :
 \pstextAC[fillstyle=solid,fillcolor=green!50,fontsize=50,Yoffset=-0.5]{PStricks}
qui comprend différentes options.
Le détail des commandes est dans la documentation. Voici quelques images, puis les fichiers.



 L'anamorphose d'une fleur de lys. Cette fleur de lys est l'un des ornements du package de Patricks Fradin psvectorian.


http://manuel.luque.perso.neuf.fr/Anamorphoses/pst-anamorphosis/
ou
http://melusine.eu.org/syracuse/G/pstricks/Anamorphoses

samedi 24 septembre 2011

Anamorphose cylindrique : partie 1

Cette première partie rappelle le principe optique du phénomène et en déduit les formules mathématiques nécessaires au calcul de l'image. Quelques commandes PStricks dédiées à l'anamorphose cylindrique permettent de dessiner l'image anamorphosée.
Une deuxième partie enrichira la base d'images et essaiera de construire des macros PStricks plus évoluées.

Les fichiers :
http://manuel.luque.perso.neuf.fr/Anamorphoses/anamorphose_cylindrique.tex

Une version plus évoluée comprenant des commandes nouvelles, que je décrirai plus tard, et un fichier .pro à placer dans le même répertoire :
http://manuel.luque.perso.neuf.fr/Anamorphoses/anacylinder.pro
http://manuel.luque.perso.neuf.fr/Anamorphoses/anamorphose_miroir_cylindrique.tex
http://manuel.luque.perso.neuf.fr/Anamorphoses/anamorphose_miroir_cylindrique.pdf

jeudi 22 septembre 2011

Lignes du champ magnétique créé par deux fils rectilignes infinis parallèles

La macro \PStricks \psfieldlinestwowires[a=2,x0=-6,y0=0,I1=1,I2=1,points=2500,H=0.05,drawarrows=false]
permet de tracer une ligne de champ dans un plan perpendiculaire aux deux fils. La représentation en 3D est aussi facilement réalisable. On peut choisir
- l'intensité de chaque courant, elles peuvent être différentes ainsi que le signe du courant.
- la distance entre les deux fils ;
- le point de départ pour le tracé de la ligne de champ ;
- le nombre de points et le pas du tracé pour ajuster au mieux la ligne ;
- la possibilité de flècher les lignes.
Voici quelques exemples :





Les fichiers .tex et .pdf sont téléchargeables ici :
http://melusine.eu.org/syracuse/G/pstricks/Champ_magnetique_2_fils/
ou
Champ_magnetique_2fils

mercredi 21 septembre 2011

Enjoliver un texte

Enjoliver un texte par des motifs abstraits constitués de courbes, volutes et arabesques ou végétaux (feuilles de houx par exemple) ou bien d'oiseaux, c'est tout un art et le faire grâce à PStricks c'est désormais possible. C'est ce qu'a réalisé Patrick Fradin à partir de la version libre des ornements de :
http://www.vectorian.net/download-free-vector-ornaments.html
 Un exemple réalisé par Patrick :
Les fichiers du package sont accessibles ici :
http://melusine.eu.org/syracuse/pstricks/vectorian/

lundi 19 septembre 2011

Transformation de Joukowski

 Les documents sur la transformation de Joukowski sont très nombreux. En voici deux sur lesquels je me suis beaucoup appuyé :
Le cours des arts&métiers :
 http://192.93.173.165/activites/cours/cours_en_ligne/cours_rey_2_TOMEII.pdf
 et le cours de l'université de Gênes :
http://www.diam.unige.it/~irro/
qui contient l'applet :
http://www.diam.unige.it/~irro/java/conformi1_1.html

Une  macro PStricks  : \psTransformJoukowski[radius=1,xc=0.1,yc=0.2] Ce sont les valeurs par défaut pour étudier la transformation de Joukowski : cercle -> profil aile d'avion (pisciforme).

 Les fichiers .tex et .pdf  :
http://melusine.eu.org/syracuse/G/pstricks/Transformation_de_Joukowski

dimanche 18 septembre 2011

Écoulement autour d'un cylindre avec circulation

C'est une première ébauche de l'écoulement d'un fluide autour d'un cylindre avec circulation. Seules les lignes d'écoulement et les équipotentielles sont tracées. Le calcul et le tracé des équipotentielles demandent beaucoup de temps. Elle se tracent toutes à la fois et non pas l'une après l'autre ce qui peut être énervant d'attendre devant l'écran vide. Mais plus tard, on pourra envisager de tracer les lignes équipotentielles les unes après les autres.
Il y a un problème non résolu lorsque k>1, la ligne d'écoulement d'arrêt ne boucle pas entièrement.
La suite sera l'écoulement autour d'un profil aérodynamique par la transformation de Joukowski.
Voici les images qui illustrent les trois cas, puis les fichiers.
Ces dessins sont très fortement inspirés par ceux de l'excellent site :
Les fichiers .tex et .pdf :
Dans le sous-dossier `macrosPstricks' du dossier:
ou bien ici :

mercredi 14 septembre 2011

Projet PStricks sous Git

Jean-Michel Sarlat a pris l'initiative de créer une page de "Projets Syracuse sous Git" :
http://melusine.eu.org/syracuse/G/
Un projet PStricks est en cours d'élaboration avec Juergen Gilg et l'aide de Jean-Michel :
http://melusine.eu.org/syracuse/G/pstricks/
Vous êtes invités à y participer, en français, anglais, allemand etc.
Jean-Michel a commencé à écrire un tutoriel pour installer git sous Windows:
http://melusine.eu.org/syracuse/G/git-windows/

Vos courbes aux couleurs de l'arc-en-ciel !

Cette commande \parametricplotHSB, est dérivée de \parametricplot du package pst-plot de Timothy Van Zand. Elle en reprend l’essentiel.
Elle s’emploie comme celle de \parametricplot, mais ne supporte pas les options de style linestyle.
Par défaut le codage HSB=true est activé, on alors une courbe qui déploie toutes les couleurs de l’arc-en-ciel sur sa longueur.
On peut choisir la couleur de début et la couleur de la fin de la courbe avec les paramètres : HueBegin=0,HueEnd=0.5, par exemple ; les valeurs de H étant choisies entre 0 et 1.
Avec l’option HSB=false, les options de couleurs classiques redeviennent opérantes.
Le nombre de points se fixe avec le paramètre : plotpoints=1000
L’option [algebraic] est possible, si l'on préfère.
Les fichiers peuvent être téléchargés ici : http://manuel.luque.perso.neuf.fr/plotHSB2011/
Le code :
\documentclass[a4paper]{article}
\usepackage[latin1]{inputenc}
\usepackage[T1]{fontenc}
\usepackage[height=24cm,width=17cm]{geometry}
\usepackage{pstricks,pst-plot,pst-slpe}

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Essai de Manuel Luque 19 février 2003
% transformé par Denis Girou le 25 février 2003
% révision le 14 septembre 2011
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\makeatletter
\pst@addfams{pst-plothsb}

\define@key[psset]{pst-plothsb}{HueBegin}{% Between 0 and 1
\edef\PstParametricplotHSB@HueBegin{#1}}

\define@key[psset]{pst-plothsb}{HueEnd}{% Between 0 and 1
\edef\PstParametricplotHSB@HueEnd{#1}}

\newif\ifPst@HSB
\define@key[psset]{pst-plothsb}{HSB}[true]{\@nameuse{Pst@HSB#1}}

% Default values
\psset{HueBegin=0,HueEnd=1,HSB=true}

\def\parametricplotHSB{\pst@object{parametricplotHSB}}
\def\parametricplotHSB@i{\@ifnextchar[{\parametricplotHSB@do}{\parametricplotHSB@do[]}}
\def\parametricplotHSB@do[#1]#2#3#4{{%
\psset{#1}%
\begin@ClosedObj
\addto@pscode{%
/t  #2 def
/dt #3 t sub \psk@plotpoints\space div def
/t t dt sub def
/Counter 0 def
\psk@plotpoints {
  /t t dt add def
  /Counter Counter 1 add def
/F@pstplot \ifPst@algebraic (#4)
                    tx@AlgToPs begin AlgToPs end cvx
                 \else { #4 }
          \fi
   def
\ifPst@algebraic
   F@pstplot
    \else  #4
\fi
  \pst@number\psyunit mul exch
  \pst@number\psxunit mul exch
  1 Counter eq
    {moveto}                    % First point
    {\ifPst@HSB                 % Other points than the first one
       /PointY exch def
       /PointX exch def
       Counter \psk@plotpoints\space div
           \PstParametricplotHSB@HueEnd\space
             \PstParametricplotHSB@HueBegin\space sub mul
           \PstParametricplotHSB@HueBegin\space add
         1 1 sethsbcolor
       PointX PointY lineto
       stroke
       PointX PointY moveto
     \else
       lineto
     \fi} ifelse
    } repeat}% fin du code ps
\end@ClosedObj}} % fin de la commande PSTricks

\makeatother

% La façon de procéder a été copiée sur pst-plot

\title{\psframebox[fillstyle=slopes]{\huge Vos courbes aux couleurs de l'arc-en-ciel !}}
\author{ Denis GIROU, Manuel LUQUE}
\date{19-25 février 2003 et 14 septembre 2011}

\begin{document}

\maketitle

Cette commande \verb+\parametricplotHSB+, est dérivée de
\verb+\parametricplot+ du package \texttt{pst-plot} de Timothy Van Zand.
Elle en reprend l'essentiel.

Elle s'emploie comme celle de \verb+\parametricplot+, mais ne
supporte pas les options de style \verb+linestyle+. Par défaut le codage
\texttt{HSB=true} est activé, on alors une courbe qui déploie toutes les
couleurs de l'arc-en-ciel sur sa longueur.

Une nouvelle option : on peut choisir la couleur de début et la
couleur de la fin de la courbe avec les paramètres :
\verb+HueBegin=0,HueEnd=0.5+, par exemple ; les valeurs de
\texttt{H} étant choisies entre 0 et 1.

Avec l'option \texttt{HSB=false}, les options de couleurs classiques
redeviennent opérantes.

Le nombre de points se fixe avec le paramètre : \texttt{plotpoints=1000}

L'option \texttt{[algebraic]} est possible.

\begin{verbatim}
\begin{pspicture}(0,-4.5)(10,4.5)
  \psgrid(0,-4)(10,4)
  \psset{plotpoints=360}
  \parametricplotHSB[linewidth=1mm,HSB=false,linecolor=red,algebraic]{0}{6.28}
     {10*t/6.28|4*sin(t)}
  \parametricplotHSB[linewidth=1mm,algebraic]{0}{6.28}{10*t/6.28|4*cos(t)}
\end{pspicture}
\end{verbatim}

\begin{center}
\begin{pspicture}(0,-5)(10,5)
  \psgrid(0,-4)(10,4)
  \psset{plotpoints=360}%
  \parametricplotHSB[linewidth=1mm,HSB=false,linecolor=red,algebraic]{0}{6.28}
     {10*t/6.28|4*sin(t)}
  \parametricplotHSB[linewidth=1mm,algebraic]{0}{6.28}{10*t/6.28|4*cos(t)}
\end{pspicture}
\end{center}

\newpage
\begin{verbatim}
\begin{pspicture}(0,-5)(10,5)
  \psgrid(0,-4)(10,4)
  \psset{plotpoints=300}%
  \parametricplotHSB[linewidth=1mm,HueBegin=0,HueEnd=0.5]{0}{360}
    {t t 36 div t sin 4 mul}
  \parametricplotHSB[linewidth=1mm,HueBegin=0.5,HueEnd=0.7]{0}{360}
    {t t 36 div t cos 4 mul}
  \parametricplotHSB[linewidth=1mm,HueBegin=0.8,HueEnd=1]{0}{360}
    {t t 36 div t sin 2 mul}
\end{pspicture}
\end{verbatim}
\begin{center}
\begin{pspicture}(0,-5)(10,5)
  \psgrid(0,-4)(10,4)
  \psset{plotpoints=300}%
  \parametricplotHSB[linewidth=1mm,HueBegin=0,HueEnd=0.5]{0}{360}
    {t t 36 div t sin 4 mul}
  \parametricplotHSB[linewidth=1mm,HueBegin=0.5,HueEnd=0.7]{0}{360}
    {t t 36 div t cos 4 mul}
  \parametricplotHSB[linewidth=1mm,HueBegin=0.8,HueEnd=1]{0}{360}
    {t t 36 div t sin 2 mul}
\end{pspicture}
\end{center}

\newpage
\begin{verbatim}
\begin{pspicture}(-5,-5)(5,5)
  \psframe*[linecolor=lightgray](-5,-4)(5,4)
  \psgrid[gridlabels=0,subgriddiv=0](-5,-4)(5,4)
  \multido{\nxDiv=-5+0.2}{50}{%
    \psline(\nxDiv,-.1)(\nxDiv,0.1)}
  \multido{\nyDiv=-4.0+0.2}{40}{%
    \psline(-0.1,\nyDiv)(0.1,\nyDiv)}
  \psset{linewidth=1mm}%
  \parametricplotHSB[plotpoints=1000,linecolor=blue,HSB=false,yunit=0.5]{-5}{5}{%
    t
    /temps t 2e-3 mul def
    /frequence2 100 def
    frequence2 360 mul temps mul cos
    1 mul 3 add
    }
  \parametricplotHSB[plotpoints=10000,yunit=1]{-5}{5}{%
    t
    /temps t 2e-3 mul def
    /frequence1 1200 def
    /frequence2 100 def
    frequence2 360 mul temps mul cos
    1 mul 3 add
    frequence1 360 mul temps mul cos
    4 mul
    mul
    0.1 mul
    }
\end{pspicture}
\end{verbatim}
\begin{center}
\begin{pspicture}(-5,-5)(5,5)
  \psframe*[linecolor=lightgray](-5,-4)(5,4)
  \psgrid[gridlabels=0,subgriddiv=0](-5,-4)(5,4)
  \multido{\nxDiv=-5+0.2}{50}{%
    \psline(\nxDiv,-.1)(\nxDiv,0.1)}
  \multido{\nyDiv=-4.0+0.2}{40}{%
    \psline(-0.1,\nyDiv)(0.1,\nyDiv)}
  \psset{linewidth=1mm}%
  \parametricplotHSB[plotpoints=1000,linecolor=blue,HSB=false,yunit=0.5]{-5}{5}{%
    t
    /temps t 2e-3 mul def
    /frequence2 100 def
    frequence2 360 mul temps mul cos
    1 mul 3 add
    }
  \parametricplotHSB[plotpoints=1000,linecolor=red,HSB=false,yunit=0.5]{-5}{5}{%
    t
    /temps t 2e-3 mul def
    /frequence2 100 def
    frequence2 360 mul temps mul cos
    -1 mul 3 sub
    }
  \parametricplotHSB[plotpoints=10000,yunit=1]{-5}{5}{%
    t
    /temps t 2e-3 mul def
    /frequence1 1200 def
    /frequence2 100 def
    frequence2 360 mul temps mul cos
    1 mul 3 add
    frequence1 360 mul temps mul cos
    4 mul
    mul
    0.1 mul
    }
\end{pspicture}
\end{center}
\newpage
\begin{verbatim}
\begin{pspicture}(-2,-2)(2,2)
  \psgrid(-2,-2)(2,2)
  \psset{plotpoints=1000}%
  \parametricplotHSB[linewidth=1,algebraic=true]{0}{6.28}
    {1.5*cos(t)|1.5*sin(t)}
\pscircle{2}\pscircle{1}
\end{pspicture}
\end{verbatim}
\begin{center}
\begin{pspicture}(-2,-2)(2,2)
  \psgrid(-2,-2)(2,2)
  \psset{plotpoints=1000}%
  \parametricplotHSB[linewidth=1,algebraic=true]{0}{6.28}
    {1.5*cos(t)|1.5*sin(t)}
\pscircle{2}\pscircle{1}
\end{pspicture}
\end{center}
\end{document}

vendredi 9 septembre 2011

Figures d'interférences avec deux fentes d'Young

C'est une commande pour dessiner les figures d'interférences lumineuses obtenues avec deux fentes d'Young.

Elle s’utilise ainsi : \psYoung(x,y){a}{D}{l}{λ}.
(x,y) centre la figure de Young au point de coordonnées (x,y) du \begin{pspicture}.
Où a est la distance entre les deux fentes
a s’exprime en m, ex : 5e-5.
l est la largeur des fentes s’exprime en m, ex : 1e-5.
D est la distance à l’écran, exprimée en m.
λ est la longueur d’onde en nm(nanomètres).
Exemples : λ = 700 nm pour le rouge, λ = 546 nm pour le vert et λ = 436 nm pour le bleu.
La figure est vue en grandeur réelle sur l’écran et la couleur s'ajuste automatiquement.
Pour obtenir la figure d’interférences en lumière blanche, il suffit de laisser vide l’argument de la longueur d’onde λ.
 \psYoung(0,0){100e-6}{2}{10e-6}{700}
 \psYoung(0,0){100e-6}{2}{10e-6}{546}
\psYoung(0,0){100e-6}{2}{10e-6}{436}
\psYoung(0,0){50e-6}{2}{10e-6}{}
Les fichiers :

Une animation :

mercredi 7 septembre 2011

Construction de Huygens du rayon réfracté

C'est une commande qui s'utilise ainsi :

\Huygens{1.5}{60} % {indice}{angle d'incidence}

Elle a été mise au point avec Juergen Gilg. Voici le résultat de la macro ci-dessus :


En modifiant les paramètres :  

\Huygens{1.33}{45}


Le code :

\documentclass[12pts]{article}
\usepackage{pstricks,pst-node}
\usepackage[latin1]{inputenc}
\usepackage[T1]{fontenc}
\usepackage[a4paper,left=1.5cm,right=1.5cm]{geometry}
\SpecialCoor

\newcommand{\Huygens}[2]{
% #1  indice
% #2  angle d'incidence
\psset{unit=5}
\begin{pspicture}*(-1.5,-1.2)(1.5,1.2)
\pstVerb{% les dimensions
/arcsin { % JP Vignault
   dup 1 eq {
      90
   } {
      dup
      dup mul neg 1 add sqrt
      atan
      dup 90 lt
         {}
         {360 sub}
      ifelse
   } ifelse
} def
    /iN #1 def % indice du milieu réfractant
    /R1 1 def
    /R2 1 iN div def
    /i1 #2 def
    /i2 i1 sin iN div arcsin def
% les points pour la construction
    /xA R1 i1 sin div def
    /yA 0 def
    /xB R2 i2 sin mul def
    /yB R2 i2 cos mul neg def
    /xC R1 i1 sin mul def
    /yC R1 i1 cos mul neg  def
    /IE xA i1 cos mul def
    /xE IE i1 cos mul def
    /yE IE i1 sin mul def
% le rayon des ondelettes
    /Rx {neg xA add i2 sin mul} def
    }%
\psframe*[linecolor=cyan!20](-1.5,-1.2)(1.5,0)
\psline(-1.5,0)(1.5,0)
\psarc[linestyle=dashed](0,0){!R1}{180}{360}
\psarc[linestyle=dashed](0,0){!R2}{180}{360}
\psline[linestyle=dotted](0,0.5)(0,-0.5)
\psarc[linestyle=dotted]{->}(0,0){0.25}{90}{!i1 90 add}
\uput{0.25}[!i1 2 div 90 add](0,0){$i_1$}
\psarc[linestyle=dotted]{->}(0,0){0.25}{270}{!i2 270 add}
\uput{0.25}[!i2 2 div 270 add](0,0){$i_2$}
\pnode(!xA yA){A}
\pnode(!xB yB){B}
\pnode(!xC yC){C}
\pnode(!xE yE){E}
\pnode(! i1 90 add dup cos 1.38 mul exch sin 1.38 mul){S}
\pnode(! i1 270 add dup cos R1 mul exch sin R1 mul){I1}
\pnode(! i2 270 add dup cos R2 mul exch sin R2 mul){I2}
\psline[linecolor=red!80](S)(0,0)
\pcline[nodesepB=0.8,linewidth=2\pslinewidth,nodesepA=0.4,linecolor=red]{->}(S)(0,0)
\psline[linestyle=dashed](0,0)(I1)
\pcline[nodesepB=-0.5,linecolor=red!80](0,0)(I2)
\pcline[nodesepB=0.45,linewidth=2\pslinewidth,linecolor=red]{->}(0,0)(I2)
\pcline[linestyle=dashed](!R2 230 cos mul R2 230 sin mul)(0,0)
\naput[nrot=:U,npos=0.5,labelsep=0.01]{$r_2=\frac{1}{n_2}$}
\uput[230](!R2 230 cos mul R2 230 sin mul){$(\mathcal{C}_2)$}
\pcline[linestyle=dashed,nodesepB=-0.5](A)(I2)
\pcline[nodesepB=-0.5](A)(E)
\pcline[linestyle=dashed,nodesepB=-0.2](0,0)(E)
\pcline[nodesepB=-0.4,linestyle=dashed](I1)(A)
\pcline[linestyle=dashed](!R1 200 cos mul R1 200 sin mul)(0,0)
\naput[nrot=:U,npos=0.3,labelsep=0.01]{$r_1=\frac{1}{n_1}$}
\uput[230](!R1 200 cos mul R1 200 sin mul){$(\mathcal{C}_1)$}
\psdots[dotstyle=o](0,0)(I2)(I1)(A)(B)(C)(E)
\rput{!i1}(0,0){\psline(0.05,0)(0.05,0.05)(0,0.05)}
\rput{!i1}(I1){\psline(0.05,0)(0.05,0.05)(0,0.05)}
\rput{!i2}(I2){\psline(0.05,0)(0.05,0.05)(0,0.05)}
\rput(-1.2,0.1){(1)}
\rput(-1.2,-0.1){$n_2>n_1$}
\rput(-1.2,-0.2){(2)}
\uput{0.05}[!i1 2 div](0,0){$I$}
\uput[r](I1){$I_1$}
\uput[r](I2){$I_2$}
\uput[u](A){$A$}
\end{pspicture}}

\begin{document}
\begin{center}
\Huygens{1.5}{60} % {indice}{angle d'incidence}
\end{center}
\begin{enumerate}
\item On dessine les demi-cercles $(\mathcal{C}_1)$ et $(\mathcal{C}_2)$ de centre $I$ et de rayons respectifs $r_1=\frac{1}{n_1}$ et $r_2=\frac{1}{n_2}$.
\item On trace le prolongement du rayon incident qui coupe $(\mathcal{C}_1)$ en $I_1$.
\item On construit la tangente au cercle $(\mathcal{C}_1)$ en $I_1$ qui coupe le plan d'incidence en $A$ ;
\item Du point $A$ on mène la tangente au cercle $(\mathcal{C}_2)$. Le point de contact $I_2$ permet de tracer le rayon réfracté $II_2$.
\end{enumerate}
\newpage
\begin{center}
\Huygens{1.33}{45}
\end{center}
\end{document}
 

On peut le télécharger ici :  Huygens

jeudi 1 septembre 2011

Les problèmes du disque tournant et du cycliste crotté d'après Henri Bouasse

Le problème du cycliste crotté  a été traité par Henri Bouasse en 1923 dans son livre "Dynamique générale" paru aux éditions Delagrave, au paragraphe 133 (pages 192 à 194). Il a récemment fait l'objet de deux articles dans Europhysics News. Bruno Macke qui m'a informé de la parution de ces deux articles, en a profité pour me signaler un oubli de ma part et une erreur dans le texte de Henri Bouasse que j'avais mis en ligne :
J'ai apporté les corrections au texte initial.
Tous les fichiers sont ici : disquetournant

Pour les animations .gif, j'ai d'abord créé les fichiers au format .eps grâce au package pst-eps, puis en mode commande DOS (sous windows) je les ai convertis au format .jpg :

FOR %G IN (*.eps) DO convert "%G" -quality 100 "%~nG.jpg"

L'animation s'obtient avec la commande suivante :

convert -delay 20 -loop 0 *.JPG animation.gif

Il faut, bien évidemment, avoir installé ImageMagick.