--- /dev/null
+\r
+\documentclass{article}\r
+\usepackage[a4paper,margin=2cm]{geometry}\r
+\usepackage[T1]{fontenc}\r
+\usepackage[latin1]{inputenc}%\r
+\usepackage[garamond]{mathdesign}\r
+\usepackage{pst-eqdf,pst-node,pst-tools}\r
+\usepackage{array,amsmath}\r
+\usepackage{animate}\r
+\newpsstyle{vecteurA}{arrowinset=0.05,arrowsize=0.1,linecolor={[rgb]{1 0.5 0}}}\r
+\newpsstyle{vecteurB}{arrowinset=0.05,arrowsize=0.1,linecolor={[rgb]{0 0.5 1}}}\r
+\newpsstyle{vecteurC}{arrowinset=0.1,arrowsize=0.2,linecolor={[rgb]{1 0 0}}}\r
+%timeline\r
+\begin{filecontents}{pb2corps.dat}\r
+::0x0\r
+::1\r
+::2\r
+::3\r
+::4\r
+::5\r
+::6\r
+::7\r
+::8\r
+::9\r
+::10\r
+::11\r
+::12\r
+::13\r
+::14\r
+::15\r
+::16\r
+::17\r
+::18\r
+::19\r
+::20\r
+::21\r
+::22\r
+::23\r
+::24\r
+::25\r
+::26\r
+::27\r
+::28\r
+::29\r
+::30\r
+::31\r
+::32\r
+::33\r
+::34\r
+::35\r
+::36\r
+::37\r
+::38\r
+::39\r
+::40\r
+::41\r
+::42\r
+::43\r
+::44\r
+::45\r
+::46\r
+::47\r
+::48\r
+::49\r
+::50\r
+::51\r
+::52\r
+::53\r
+::54\r
+::55\r
+::56\r
+::57\r
+::58\r
+::59\r
+::60\r
+::61\r
+::62\r
+::63\r
+::64\r
+::65\r
+::66\r
+::67\r
+::68\r
+::69\r
+::70\r
+::71\r
+::72\r
+::73\r
+::74\r
+::75\r
+::76\r
+::77\r
+::78\r
+::79\r
+::80\r
+::81\r
+::82\r
+::83\r
+::84\r
+::85\r
+::86\r
+::87\r
+::88\r
+::89\r
+::90\r
+::91\r
+::92\r
+::93\r
+::94\r
+::95\r
+::96\r
+::97\r
+::98\r
+::99\r
+::100\r
+::101\r
+::102\r
+::103\r
+::104\r
+::105\r
+::106\r
+::107\r
+::108\r
+::109\r
+::110\r
+::111\r
+::112\r
+::113\r
+::114\r
+::115\r
+::116\r
+::117\r
+::118\r
+::119\r
+::120\r
+::121\r
+::122\r
+::123\r
+::124\r
+::125\r
+::126\r
+::127\r
+::128\r
+::129\r
+::130\r
+::131\r
+::132\r
+::133\r
+::134\r
+::135\r
+::136\r
+::137\r
+::138\r
+::139\r
+::140\r
+::141\r
+::142\r
+::143\r
+::144\r
+::145\r
+::146\r
+::147\r
+::148\r
+::149\r
+::150\r
+::151\r
+::152\r
+::153\r
+::154\r
+::155\r
+::156\r
+::157\r
+::158\r
+::159\r
+::160\r
+::161\r
+::162\r
+::163\r
+::164\r
+::165\r
+::166\r
+::167\r
+::168\r
+::169\r
+::170\r
+::171\r
+::172\r
+::173\r
+::174\r
+::175\r
+::176\r
+::177\r
+::178\r
+::179\r
+::180\r
+::181\r
+::182\r
+::183\r
+::184\r
+::185\r
+::186\r
+::187\r
+::188\r
+::189\r
+::190\r
+::191\r
+::192\r
+::193\r
+::194\r
+::195\r
+::196\r
+::197\r
+::198\r
+::199\r
+::200\r
+\end{filecontents}\r
+\title{Gravitation : le problème des deux corps avec PSTricks\\ partie 3}\r
+\date{6 juillet 2\,012}\r
+\begin{document}\r
+\maketitle\r
+\section{Mouvements dans le repère du centre de masse}\r
+D'abord un petit rappel des documents précédents :\r
+\r
+On considère un système de deux corps en interaction gravitationnelle $M_1$ de masse $m_1$ et $M_2$ de masse $m_2$ dans le repère galiléen \textit{inertiel} $\mathcal{R}$. Ils sont supposés ponctuels.\r
+\begin{center}\r
+\begin{pspicture}(-3,-1)(5,6)\r
+\psgrid[subgriddiv=0,gridcolor=lightgray,griddots=10,gridlabels=0pt]%\r
+\pstVerb{/x01 -2 def /y01 1 def\r
+ /x02 5 def /y02 4 def\r
+ /xr0 x02 x01 sub def\r
+ /yr0 y02 y01 sub def\r
+ /M1 3 def\r
+ /M2 1 def\r
+ /Mt M1 M2 add def\r
+ /xG0 x01 M1 mul x02 M2 mul add Mt div def\r
+ /yG0 y01 M1 mul y02 M2 mul add Mt div def}%\r
+\pnode(0,0){O}\r
+\pnode(!x01 y01){M1}\r
+\pnode(!x02 y02){M2}\r
+\pnode(!xG0 yG0){C}\r
+\pscircle[fillstyle=solid,fillcolor=red!50](M1){0.21}\r
+\pscircle[fillstyle=solid,fillcolor=blue!50](M2){0.07}\r
+\psdot(C)\r
+\psline[linestyle=dotted](M1)(M2)\r
+\psline[linestyle=dashed](M2)(O)(M1)\r
+\psline{<->}(5,0)(0,0)(0,6)\r
+\uput[dl](O){$O$}\r
+\uput[r](0,5.8){$y$}\r
+\uput[u](4.8,0){$x$}\r
+\uput[l](0,5.8){$\mathcal{R}$}\r
+\uput{0.22}[l](M1){$M_1$}\r
+\uput{0.1}[u](M2){$M_2$}\r
+\uput[u](C){$C$}\r
+\pcline[offset=5pt,linestyle=none](M1)(M2)\r
+\ncput*[nrot=:U]{$r$}\r
+\rput(M1){\psline[style=vecteurC]{->}(!xr0 5 div yr0 5 div)\uput[u](!xr0 10 div yr0 10 div){$\overrightarrow{F}_{2/1}$}}\r
+\rput(M2){\psline[style=vecteurC]{->}(!xr0 5 div neg yr0 5 div neg)\uput[u](!xr0 10 div neg yr0 10 div neg){$\overrightarrow{F}_{1/2}$}}\r
+\end{pspicture}\r
+\end{center}\r
+On note $\overrightarrow{r}=\overrightarrow{M_1M_2}$. $M_2$ subit de la part de $M_1$ une force attractive $\overrightarrow{F}_{1/2}$ et réciproquement $M_1$ subit de la part de $M_2$ une force attractive $\overrightarrow{F}_{2/1}$ telles que :\r
+\[\r
+\overrightarrow{F}_{1/2}=-\mathcal{G}\frac{m_1m_2}{r^3}\overrightarrow{r}\quad \overrightarrow{F}_{2/1}=+\mathcal{G}\frac{m_1m_2}{r^3}\overrightarrow{r}\r
+\]\r
+La loi de Newton appliquée à chacune des particules s'écrit :\r
+\[\r
+m_1\frac{\mathrm{d}^2\overrightarrow{r}}{\mathrm{d}t^2}=\hphantom{-}\mathcal{G}\frac{m_1m_2}{r^3}\overrightarrow{r}\r
+\]\r
+\[\r
+m_2\frac{\mathrm{d}^2\overrightarrow{r}}{\mathrm{d}t^2}=-\mathcal{G}\frac{m_1m_2}{r^3}\overrightarrow{r}\r
+\]\r
+Sachant que $r=\sqrt{(x_{2}-x_{1})^2+(y_{2}-y_{1})^2}$, nous obtenons un système de quatre équations différentielles :\r
+\[\r
+\left\{\r
+\begin{array}[m]{l}\r
+ \ddot{x_1}=\hphantom{-}\mathcal{G}\displaystyle\frac{m_2}{r^3}(x_2-x_1)\\[1em]\r
+ \ddot{y_1}=\hphantom{-}\mathcal{G}\displaystyle\frac{m_2}{r^3}(y_2-y_1)\\[1em]\r
+ \ddot{x_2}=-\mathcal{G}\displaystyle\frac{m_1}{r^3}(x_2-x_1)\\[1em]\r
+ \ddot{y_2}=-\mathcal{G}\displaystyle\frac{m_1}{r^3}(y_2-y_1)\r
+\end{array}\r
+\right.\r
+\]\r
+\r
+% x y x' y'\r
+% y[0] y[1] y[2] y[3]\r
+\def\FictifAlg{%\r
+ y[2]|y[3]|%\r
+ -(M1+M2)*y[0]/(y[0]^2+y[1]^2)^1.5|%\r
+ -(M1+M2)*y[1]/(y[0]^2+y[1]^2)^1.5}\r
+Voici comment placer sur la pile les différentes variables, ainsi que le système d'équations en notation algébrique. Il faut retenir la notation utilisée pour représenter les variables en fonction de leurs positions respectives sur la pile :\r
+\begin{verbatim}\r
+% 0 1 2 3 4 5 6 7\r
+% y[0] y[1] y[2] y[3] y[4] y[5] y[6] y[7]\r
+% x1 y1 x'1 y'1 x2 y2 x'2 y'2\r
+\def\GravAlgIIcorps{%\r
+ y[2]|y[3]|%\r
+ M2*(y[4]-y[0])/((y[4]-y[0])^2+(y[5]-y[1])^2)^1.5|%\r
+ M2*(y[5]-y[1])/((y[4]-y[0])^2+(y[5]-y[1])^2)^1.5|%\r
+ y[6]|y[7]|%\r
+ M1*(y[0]-y[4])/((y[4]-y[0])^2+(y[5]-y[1])^2)^1.5|%\r
+ M1*(y[1]-y[5])/((y[4]-y[0])^2+(y[5]-y[1])^2)^1.5}\r
+\end{verbatim}\r
+Les conditions initiales, positions et vitesses respectives des deux corps, déterminent l'évolution du système.\r
+\begin{verbatim}\r
+\def\InitCond{x01 y01 v0x1 v0y1 x02 y02 v0x2 v0y2}\r
+\end{verbatim}\r
+Elles sont définies dans l'environnement \verb+\begin{pspicture}+ :\r
+\begin{verbatim}\r
+\pstVerb{\r
+ /x01 2 def\r
+ /y01 2 def\r
+ /v0x1 .2 def\r
+ /v0y1 0.25 def\r
+ /x02 -3 def\r
+ /y02 0 def\r
+ /v0x2 -0.25 def\r
+ /v0y2 -0.5 def\r
+ /M1 3 def\r
+ /M2 1 def\r
+ }%\r
+\end{verbatim}\r
+Dans cette même variable \verb+\pstVerb+, il est utile de faire d'autres calculs, comme la période du mouvement afin de bien choisir l'intervalle de temps sur lequel on dessine les trajectoires. Les calculs s'effectuent en \textsf{postscript} et l'affichage du résultat utilise, par exemple, le package \textsf{pst-tools}.\r
+\begin{verbatim}\r
+\rput(-2,0){\psPrintValue[decimals=4]{periode}\hphantom{00000}s}\r
+\end{verbatim}\r
+On désactivera cette ligne par \% au début, une fois la période connue. Dans l'exemple proposé $T=41,0688\,\mathrm{s}$. On étudie le mouvement sur trois périodes environ.\r
+%% 0 1 2 3 4 5 6 7\r
+%% x1 y1 x'1 y'1 x2 y2 x'2 y'2\r
+\def\GravAlgIIcorps{%\r
+ y[2]|y[3]|%\r
+ M2*(y[4]-y[0])/((y[4]-y[0])^2+(y[5]-y[1])^2)^1.5|%\r
+ M2*(y[5]-y[1])/((y[4]-y[0])^2+(y[5]-y[1])^2)^1.5|%\r
+ y[6]|y[7]|%\r
+ M1*(y[0]-y[4])/((y[4]-y[0])^2+(y[5]-y[1])^2)^1.5|%\r
+ M1*(y[1]-y[5])/((y[4]-y[0])^2+(y[5]-y[1])^2)^1.5}\r
+%\newpage\r
+\r
+Dans $\mathcal{R}$ la position de $C$ est déterminée par la relation :\r
+\[\r
+\overrightarrow{OC}=\frac{m_1\overrightarrow{OM_1}+m_2\overrightarrow{OM_2}}{m_1+m_2}\r
+\]\r
+Avec les conventions utilisées les coordonnées de $C$ seront calculées en postscript ainsi :\r
+\begin{verbatim}\r
+ M1 y[0] mul M2 y[4] mul add M1 M2 add div % abscisse de C\r
+ M1 y[1] mul M2 y[5] mul add M1 M2 add div % ordonnée de C\r
+\end{verbatim}\r
+\r
+Dans le repère $\mathcal{R}^*$ lié au centre de masse, les coordonnées des deux points se calculent par :\r
+\r
+\[\r
+\left\{\r
+\begin{array}[m]{l}\r
+x^*_1=x_1-x_C\\\r
+y^*_1=x_1-y_C\\\r
+x^*_2=x_2-x_C\\\r
+y^*_2=y_2-y_C\\\r
+\end{array}\r
+\right.\r
+\]\r
+Ce qui en postscript s'écrira :\r
+\begin{verbatim}\r
+ y[0] % abscisse de M1\r
+ M1 y[0] mul M2 y[4] mul add M1 M2 add div % abscisse de C\r
+ sub\r
+ y[1] % ordonnée de M1\r
+ M1 y[1] mul M2 y[5] mul add M1 M2 add div % ordonnée de C\r
+ sub\r
+\end{verbatim}\r
+\begin{verbatim}\r
+ y[4] % abscisse de M2\r
+ M1 y[0] mul M2 y[4] mul add M1 M2 add div % abscisse de C\r
+ sub\r
+ y[5] % ordonnée de M2\r
+ M1 y[1] mul M2 y[5] mul add M1 M2 add div % ordonnée de C\r
+ sub\r
+\end{verbatim}\r
+\r
+Les deux options \texttt{[plotfuncx]} et \texttt{[plotfuncy]} permettent justement de faire les calculs nécessaires à ce changement de repère. \texttt{Mt} désignant la masse totale calculée dans la variable \verb+\pstVerb+ :\r
+\begin{verbatim}\r
+ /Mt M1 M2 add def\r
+\end{verbatim}\r
+\begin{verbatim}\r
+ plotfuncx=y 0 get\r
+ y 4 get M2 mul\r
+ y 0 get M1 mul add\r
+ Mt div sub ,\r
+ plotfuncy=1 get\r
+ y 1 get\r
+ M1 mul\r
+ y 5 get M2 mul add\r
+ Mt div sub,\r
+\end{verbatim}\r
+Par défaut, c'est la variable de la fonction étudiée, ici le temps $t$, qui est notée \texttt{x}, qui est placée en abscisse. Toutes les grandeurs calculées :\r
+\begin{verbatim}\r
+% 0 1 2 3 4 5 6 7\r
+% y[0] y[1] y[2] y[3] y[4] y[5] y[6] y[7]\r
+% x1 y1 x'1 y'1 x2 y2 x'2 y'2\r
+\end{verbatim}\r
+sont placées(empilées) dans un tableau noté \texttt{y}, c'est donc dans ce tableau qu'il faut \textit{aller à la pêche} avec l'instruction \textsf{postscript} \texttt{get}.\r
+\r
+Ainsi |\texttt{y 0 get}| donne : \texttt{y[0]} soit \texttt{x1}, |\texttt{y 1 get}| donne : \texttt{y[1]} soit \texttt{y1} etc.\r
+\begin{center}\r
+\begin{pspicture}(-6,-8)(3,3)\r
+\pstVerb{\r
+/arccos {\r
+ dup\r
+ dup mul neg 1 add sqrt\r
+ exch\r
+ atan\r
+} def\r
+ /G 1 def\r
+ /x01 2 def\r
+ /y01 2 def\r
+ /v0x1 .2 def\r
+ /v0y1 0.25 def\r
+ /x02 -3 def\r
+ /y02 0 def\r
+ /v0x2 -0.25 def\r
+ /v0y2 -0.5 def\r
+ /M1 3 def\r
+ /M2 1 def\r
+ /Mt M1 M2 add def\r
+ /mu M1 M2 mul M1 M2 add div def % masse réduite\r
+ /K G Mt mul def\r
+ /xG0 M1 x01 mul M2 x02 mul add M1 M2 add div def\r
+ /yG0 M1 y01 mul M2 y02 mul add M1 M2 add div def\r
+ /vG0x M1 v0x1 mul M2 v0x2 mul add M1 M2 add div def\r
+ /vG0y M1 v0y1 mul M2 v0y2 mul add M1 M2 add div def\r
+ /K1 M2 Mt div neg def\r
+ /K2 M1 Mt div def\r
+% conditions initiales pour le point réduit\r
+ /xr0 x02 x01 sub def\r
+ /yr0 y02 y01 sub def\r
+ /theta_0 yr0 xr0 atan def\r
+ /r0 xr0 dup mul yr0 dup mul add sqrt def\r
+ /v0xr v0x2 v0x1 sub def\r
+ /v0yr v0y2 v0y1 sub def\r
+ /v0r_2 v0xr dup mul v0yr dup mul add def\r
+% constante des aires\r
+ /Cste xr0 v0yr mul yr0 v0xr mul sub def % Cste des aires\r
+ /Energie 0.5 mu mul v0r_2 mul % 1/2 mv^2\r
+ G M1 M2 mul mul r0 div sub\r
+ def\r
+ /par Cste dup mul K div def % p\r
+ /exc 2 Cste dup mul mul Energie mul mu div K dup mul div 1 add sqrt def % e\r
+ /a_2 par 1 exc dup mul sub div def % demi-grand axe\r
+ /b_2 a_2 1 exc dup mul sub sqrt mul def % demi-petit axe\r
+ /c_2 a_2 exc mul def % distance focale\r
+ % phase\r
+ /Phi theta_0 par r0 div 1 sub exc div arccos sub def\r
+ /rP par 1 exc add div def\r
+ /rA par 1 exc sub div def\r
+% vitesses à l'apogée et au périgée\r
+ /vA G Mt mul 2 rA div 1 a_2 div sub mul sqrt def\r
+ /vP G Mt mul 2 rP div 1 a_2 div sub mul sqrt def\r
+% positions du périgée et de l'apogée\r
+ /xP rP Phi cos mul def\r
+ /yP rP Phi sin mul def\r
+ /xA rA Phi cos mul neg def\r
+ /yA rA Phi sin mul neg def\r
+ /xW xA xP add 2 div def\r
+ /yW yA yP add 2 div def\r
+% periode\r
+ /periode 6.28 a_2 3 exp G div Mt div sqrt mul def\r
+ /radius {par 1 exc t Phi sub cos mul add div} def\r
+ /xE {radius t cos mul} def\r
+ /yE {radius t sin mul}def\r
+ }%\r
+%% 0 1 2 3 4 5 6 7\r
+%% x1 y1 x'1 y'1 x2 y2 x'2 y'2\r
+\def\InitCond{ x01 y01 v0x1 v0y1 x02 y02 v0x2 v0y2}\r
+\psset{method=rk4}\r
+\psgrid[subgriddiv=0,gridcolor=lightgray,griddots=10,gridlabels=0pt](-6,-8)(3,3)%\r
+% mouvement de M2 par rapport à M1\r
+%\psequadiff[plotpoints=1000,algebraic,\r
+% plotfuncx=y dup 4 get exch 0 get sub ,\r
+% plotfuncy=dup 5 get exch 1 get sub,\r
+% tabname=XY1]{0}{22.7}{\InitCond}{\GravAlgIIcorps}\r
+%\listplot[unit=1,linecolor=green]{XY1 aload pop}\r
+% mouvement de M1 par rapport à M2\r
+%\psequadiff[plotpoints=1000,algebraic,\r
+% plotfuncx=y dup 0 get exch 4 get sub ,\r
+% plotfuncy=dup 1 get exch 5 get sub,\r
+% tabname=XY2]{0}{22.7}{\InitCond}{\GravAlgIIcorps}\r
+%\listplot[unit=1,linecolor=magenta]{XY2 aload pop}\r
+% mouvement de M1 par rapport à G\r
+\psequadiff[plotpoints=1000,algebraic,\r
+ plotfuncx=y 0 get\r
+ y 4 get M2 mul\r
+ y 0 get M1 mul add\r
+ Mt div sub ,\r
+ plotfuncy=1 get\r
+ y 1 get\r
+ M1 mul\r
+ y 5 get M2 mul add\r
+ Mt div sub,\r
+ tabname=XY1G,\r
+ saveData,filename=XY1G.dat\r
+ ]{0}{41.07}{\InitCond}{\GravAlgIIcorps}\r
+\listplot[unit=1,linecolor=red]{XY1G aload pop}\r
+% mouvement de M2 par rapport à G\r
+\psequadiff[plotpoints=1000,algebraic,\r
+ plotfuncx=y 4 get\r
+ y 4 get M2 mul\r
+ y 0 get M1 mul add\r
+ Mt div sub ,\r
+ plotfuncy=5 get\r
+ y 1 get M1 mul\r
+ y 5 get M2 mul add\r
+ Mt div sub,\r
+ tabname=XY2G,\r
+ saveData,filename=XY2G.dat\r
+ ]{0}{41.07}{\InitCond}{\GravAlgIIcorps}\r
+\listplot[unit=1,linecolor=blue]{XY2G aload pop}\r
+% centre de masse\r
+%\psequadiff[plotpoints=1000,algebraic,\r
+% plotfuncx=y dup 0 get M1 mul exch\r
+% 4 get M2 mul add\r
+% Mt div,\r
+% plotfuncy=dup 1 get M1 mul exch\r
+% 5 get M2 mul add\r
+% Mt div,\r
+% tabname=XGYG\r
+% ,saveData,filename=XYG.dat\r
+% ]{0}{120}{\InitCond}{\GravAlgIIcorps}\r
+%\listplot[unit=1,linecolor=cyan]{XGYG aload pop}\r
+\pnode(!x01 xG0 sub y01 yG0 sub){M01}\r
+\pnode(!x02 xG0 sub y02 yG0 sub){M02}\r
+%\pnode(!xG0 yG0){G0}\r
+\psdots(M01)(M02)\r
+%\psdot[linecolor=cyan](G0)\r
+\rput(M01){\psline[unit=5,linecolor=red]{->}(!v0x1 vG0x sub v0y1 vG0y sub)}\r
+\rput(M02){\psline[unit=5,linecolor=blue]{->}(!v0x2 vG0x sub v0y2 vG0y sub)}\r
+%\rput(G0){\psline[unit=5]{->}(!vG0x vG0y)}\r
+\psline{<->}(3,0)(0,0)(0,3)\r
+\rput(-2,2){\psPrintValue[decimals=4]{periode}\hphantom{0000000}s}\r
+\end{pspicture}\r
+\end{center}\r
+%\end{document}\r
+La partie suivante est une animation en \textit{temps réel}.\r
+\r
+Pour faciliter les calculs les positions successives des deux corps ont été enregistrées, avec l'option \texttt{saveData}, dans deux fichiers. Voici pour le premier :\r
+\begin{verbatim}\r
+\psequadiff[plotpoints=1000,algebraic,\r
+ plotfuncx=y 0 get\r
+ y 4 get M2 mul\r
+ y 0 get M1 mul add\r
+ Mt div sub ,\r
+ plotfuncy=1 get\r
+ y 1 get\r
+ M1 mul\r
+ y 5 get M2 mul add\r
+ Mt div sub,\r
+ tabname=XY1G,\r
+ saveData,filename=XY1G.dat\r
+ ]{0}{41.07}{\InitCond}{\GravAlgIIcorps}\r
+\end{verbatim}\r
+Cette option est ensuite désactivée. Pour l'animation, les fichiers sont lus et les données placées dans un tableau.\r
+\begin{verbatim}\r
+\pstVerb{/XY1 [(XYM1.dat) run] def\r
+ /XY2 [(XYM2.dat) run] def\r
+\end{verbatim}\r
+On pourra, éventuellement, tracer un segment joignant les deux corps en activant la ligne :\r
+\begin{verbatim}\r
+ %\psline[linestyle=dashed,linecolor=white](!X1 Y1)(!X2 Y2)\r
+\end{verbatim}\r
+\r
+\newpage\r
+\section{L'animation}\r
+\begin{center}\r
+\pstVerb{/XY1 [(XY1G.dat) run] def\r
+ /XY2 [(XY2G.dat) run] def\r
+ }%\r
+\def\nFrames{200}% 200 images\r
+\begin{animateinline}[controls,loop,timeline=pb2corps.dat,%\r
+ begin={\begin{pspicture}(-6,-8)(3,3)},\r
+ end={\end{pspicture}}]{20}% 20 images/s\r
+\psframe*[linecolor={[cmyk]{1 1 0 0.7}}](-6,-8)(3,3)\r
+\psgrid[subgriddiv=0,gridcolor=white,griddots=10,gridlabels=0pt]%\r
+\listplot[linecolor=gray]{XY1 aload pop}\r
+\listplot[linecolor=gray]{XY2 aload pop}\r
+\psdot[linecolor=yellow,dotstyle=+](0,0)\r
+\newframe\r
+\multiframe{\nFrames}{i=0+10}{% 1 point sur 10\r
+\pstVerb{/X1 XY1 \i\space get def\r
+ /Y1 XY1 \i\space 1 add get def\r
+ /X2 XY2 \i\space get def\r
+ /Y2 XY2 \i\space 1 add get def\r
+ }%\r
+%\psdot(!X1 Y1)\r
+\pscircle*[linecolor=red](!X1 Y1){0.3}\r
+\pscircle*[linecolor=blue!50](!X2 Y2){0.1}\r
+%\psline[linestyle=dashed,linecolor=white](!X1 Y1)(!X2 Y2)\r
+}\r
+\end{animateinline}\r
+\end{center}\r
+\section{Mouvement d'un corps par rapport à l'autre}\r
+Étudions le mouvement de $M_2$ par rapport à $M_1$. Dans un repère lié à $M_1$ les coordonnées de $M_2$ s'expriment par :\r
+\[\r
+\left\{\r
+\begin{array}[m]{l}\r
+x_{2/1}=x_2-x_1\\\r
+y_{2/1}=y_2-y_1\r
+\end{array}\r
+\right.\r
+\]\r
+\begin{verbatim}\r
+% 0 1 2 3 4 5 6 7\r
+% y[0] y[1] y[2] y[3] y[4] y[5] y[6] y[7]\r
+% x1 y1 x'1 y'1 x2 y2 x'2 y'2\r
+ y[4] y[0] sub % x2 - x1\r
+ y[5] y[1] sub % y2 - y1\r
+\end{verbatim}\r
+\r
+\begin{verbatim}\r
+% mouvement de M2 par rapport à M1\r
+\psequadiff[plotpoints=1000,algebraic,\r
+ plotfuncx=y dup 4 get exch 0 get sub ,\r
+ plotfuncy=dup 5 get exch 1 get sub,\r
+ tabname=XY21,\r
+ saveData,filename=XY21.dat\r
+ ]{0}{41.07}{\InitCond}{\GravAlgIIcorps}\r
+\end{verbatim}\r
+\begin{center}\r
+\begin{pspicture}(-7,-9)(3,3)\r
+\pstVerb{\r
+/arccos {\r
+ dup\r
+ dup mul neg 1 add sqrt\r
+ exch\r
+ atan\r
+} def\r
+ /G 1 def\r
+ /x01 2 def\r
+ /y01 2 def\r
+ /v0x1 .2 def\r
+ /v0y1 0.25 def\r
+ /x02 -3 def\r
+ /y02 0 def\r
+ /v0x2 -0.25 def\r
+ /v0y2 -0.5 def\r
+ /M1 3 def\r
+ /M2 1 def\r
+ /Mt M1 M2 add def\r
+ /mu M1 M2 mul M1 M2 add div def % masse réduite\r
+ /K G Mt mul def\r
+ /xG0 M1 x01 mul M2 x02 mul add M1 M2 add div def\r
+ /yG0 M1 y01 mul M2 y02 mul add M1 M2 add div def\r
+ /vG0x M1 v0x1 mul M2 v0x2 mul add M1 M2 add div def\r
+ /vG0y M1 v0y1 mul M2 v0y2 mul add M1 M2 add div def\r
+ /K1 M2 Mt div neg def\r
+ /K2 M1 Mt div def\r
+% conditions initiales pour le point réduit\r
+ /xr0 x02 x01 sub def\r
+ /yr0 y02 y01 sub def\r
+ /theta_0 yr0 xr0 atan def\r
+ /r0 xr0 dup mul yr0 dup mul add sqrt def\r
+ /v0xr v0x2 v0x1 sub def\r
+ /v0yr v0y2 v0y1 sub def\r
+ /v0r_2 v0xr dup mul v0yr dup mul add def\r
+% constante des aires\r
+ /Cste xr0 v0yr mul yr0 v0xr mul sub def % Cste des aires\r
+ /Energie 0.5 mu mul v0r_2 mul % 1/2 mv^2\r
+ G M1 M2 mul mul r0 div sub\r
+ def\r
+ /par Cste dup mul K div def % p\r
+ /exc 2 Cste dup mul mul Energie mul mu div K dup mul div 1 add sqrt def % e\r
+ /a_2 par 1 exc dup mul sub div def % demi-grand axe\r
+ /b_2 a_2 1 exc dup mul sub sqrt mul def % demi-petit axe\r
+ /c_2 a_2 exc mul def % distance focale\r
+ % phase\r
+ /Phi theta_0 par r0 div 1 sub exc div arccos sub def\r
+ /rP par 1 exc add div def\r
+ /rA par 1 exc sub div def\r
+% vitesses à l'apogée et au périgée\r
+ /vA G Mt mul 2 rA div 1 a_2 div sub mul sqrt def\r
+ /vP G Mt mul 2 rP div 1 a_2 div sub mul sqrt def\r
+% positions du périgée et de l'apogée\r
+ /xP rP Phi cos mul def\r
+ /yP rP Phi sin mul def\r
+ /xA rA Phi cos mul neg def\r
+ /yA rA Phi sin mul neg def\r
+ /xW xA xP add 2 div def\r
+ /yW yA yP add 2 div def\r
+% periode\r
+ /periode 6.28 a_2 3 exp G div Mt div sqrt mul def\r
+ /radius {par 1 exc t Phi sub cos mul add div} def\r
+ /xE {radius t cos mul} def\r
+ /yE {radius t sin mul}def\r
+ }%\r
+%% 0 1 2 3 4 5 6 7\r
+%% x1 y1 x'1 y'1 x2 y2 x'2 y'2\r
+\def\InitCond{ x01 y01 v0x1 v0y1 x02 y02 v0x2 v0y2}\r
+\psset{method=rk4}\r
+\psgrid[subgriddiv=0,gridcolor=lightgray,griddots=10,gridlabels=0pt](-7,-9)(3,3)%\r
+% mouvement de M2 par rapport à M1\r
+\psequadiff[plotpoints=1000,algebraic,\r
+ plotfuncx=y dup 4 get exch 0 get sub ,\r
+ plotfuncy=dup 5 get exch 1 get sub,\r
+ tabname=XY21,\r
+ saveData,filename=XY21.dat\r
+ ]{0}{41.07}{\InitCond}{\GravAlgIIcorps}\r
+\listplot[unit=1,linecolor=green]{XY21 aload pop}\r
+% mouvement de M1 par rapport à M2\r
+%\psequadiff[plotpoints=1000,algebraic,\r
+% plotfuncx=y dup 0 get exch 4 get sub ,\r
+% plotfuncy=dup 1 get exch 5 get sub,\r
+% tabname=XY2]{0}{22.7}{\InitCond}{\GravAlgIIcorps}\r
+%\listplot[unit=1,linecolor=magenta]{XY2 aload pop}\r
+\pnode(!x02 x01 sub y02 y01 sub){M021}\r
+\pscircle*[linecolor=red](0,0){0.3}\r
+\pscircle*[linecolor=blue!50](M021){0.1}\r
+\uput[r](M021){$(t_0)$}\r
+\psdot(M021)\uput[l](M021){$M_{2/1}$}\r
+\rput(M021){\psline[unit=5,linecolor=blue]{->}(!v0x2 v0x1 sub v0y2 v0y1 sub)}\r
+\uput{0.3}[d](0,0){$M_1(0)$}\r
+\psline{<->}(3,0)(0,0)(0,3)\r
+\rput(-2,2){\psPrintValue[decimals=4]{periode}\hphantom{0000000}s}\r
+\end{pspicture}\r
+\end{center}\r
+Le mouvement de $M_1$ par rapport à $M_2$ :\r
+\begin{verbatim}\r
+% mouvement de M1 par rapport à M1\r
+\psequadiff[plotpoints=1000,algebraic,\r
+ plotfuncx=y dup 0 get exch 4 get sub ,\r
+ plotfuncy=dup 1 get exch 5 get sub,\r
+ tabname=XY12,\r
+ saveData,filename=XY12.dat\r
+ ]{0}{41.07}{\InitCond}{\GravAlgIIcorps}\r
+\end{verbatim}\r
+\begin{center}\r
+\begin{pspicture}(-3,-2)(8,8)\r
+\pstVerb{\r
+/arccos {\r
+ dup\r
+ dup mul neg 1 add sqrt\r
+ exch\r
+ atan\r
+} def\r
+ /G 1 def\r
+ /x01 2 def\r
+ /y01 2 def\r
+ /v0x1 .2 def\r
+ /v0y1 0.25 def\r
+ /x02 -3 def\r
+ /y02 0 def\r
+ /v0x2 -0.25 def\r
+ /v0y2 -0.5 def\r
+ /M1 3 def\r
+ /M2 1 def\r
+ /Mt M1 M2 add def\r
+ /mu M1 M2 mul M1 M2 add div def % masse réduite\r
+ /K G Mt mul def\r
+ /xG0 M1 x01 mul M2 x02 mul add M1 M2 add div def\r
+ /yG0 M1 y01 mul M2 y02 mul add M1 M2 add div def\r
+ /vG0x M1 v0x1 mul M2 v0x2 mul add M1 M2 add div def\r
+ /vG0y M1 v0y1 mul M2 v0y2 mul add M1 M2 add div def\r
+ /K1 M2 Mt div neg def\r
+ /K2 M1 Mt div def\r
+% conditions initiales pour le point réduit\r
+ /xr0 x02 x01 sub def\r
+ /yr0 y02 y01 sub def\r
+ /theta_0 yr0 xr0 atan def\r
+ /r0 xr0 dup mul yr0 dup mul add sqrt def\r
+ /v0xr v0x2 v0x1 sub def\r
+ /v0yr v0y2 v0y1 sub def\r
+ /v0r_2 v0xr dup mul v0yr dup mul add def\r
+% constante des aires\r
+ /Cste xr0 v0yr mul yr0 v0xr mul sub def % Cste des aires\r
+ /Energie 0.5 mu mul v0r_2 mul % 1/2 mv^2\r
+ G M1 M2 mul mul r0 div sub\r
+ def\r
+ /par Cste dup mul K div def % p\r
+ /exc 2 Cste dup mul mul Energie mul mu div K dup mul div 1 add sqrt def % e\r
+ /a_2 par 1 exc dup mul sub div def % demi-grand axe\r
+ /b_2 a_2 1 exc dup mul sub sqrt mul def % demi-petit axe\r
+ /c_2 a_2 exc mul def % distance focale\r
+ % phase\r
+ /Phi theta_0 par r0 div 1 sub exc div arccos sub def\r
+ /rP par 1 exc add div def\r
+ /rA par 1 exc sub div def\r
+% vitesses à l'apogée et au périgée\r
+ /vA G Mt mul 2 rA div 1 a_2 div sub mul sqrt def\r
+ /vP G Mt mul 2 rP div 1 a_2 div sub mul sqrt def\r
+% positions du périgée et de l'apogée\r
+ /xP rP Phi cos mul def\r
+ /yP rP Phi sin mul def\r
+ /xA rA Phi cos mul neg def\r
+ /yA rA Phi sin mul neg def\r
+ /xW xA xP add 2 div def\r
+ /yW yA yP add 2 div def\r
+% periode\r
+ /periode 6.28 a_2 3 exp G div Mt div sqrt mul def\r
+ /radius {par 1 exc t Phi sub cos mul add div} def\r
+ /xE {radius t cos mul} def\r
+ /yE {radius t sin mul}def\r
+ }%\r
+%% 0 1 2 3 4 5 6 7\r
+%% x1 y1 x'1 y'1 x2 y2 x'2 y'2\r
+\def\InitCond{ x01 y01 v0x1 v0y1 x02 y02 v0x2 v0y2}\r
+\psset{method=rk4}\r
+\psgrid[subgriddiv=0,gridcolor=lightgray,griddots=10,gridlabels=0pt](-3,-2)(8,9)%\r
+% mouvement de M2 par rapport à M1\r
+%\psequadiff[plotpoints=1000,algebraic,\r
+% plotfuncx=y dup 4 get exch 0 get sub ,\r
+% plotfuncy=dup 5 get exch 1 get sub,\r
+% tabname=XY21,\r
+% saveData,filename=XY21.dat\r
+% ]{0}{41.07}{\InitCond}{\GravAlgIIcorps}\r
+%\listplot[unit=1,linecolor=green]{XY21 aload pop}\r
+% mouvement de M1 par rapport à M2\r
+\psequadiff[plotpoints=1000,algebraic,\r
+ plotfuncx=y dup 0 get exch 4 get sub ,\r
+ plotfuncy=dup 1 get exch 5 get sub,\r
+ tabname=XY12,\r
+ saveData,filename=XY12.dat\r
+ ]{0}{41.07}{\InitCond}{\GravAlgIIcorps}\r
+\listplot[unit=1,linecolor=magenta]{XY12 aload pop}\r
+\pnode(!x01 x02 sub y01 y02 sub){M012}\r
+\pscircle*[linecolor=red](M012){0.3}\r
+\pscircle*[linecolor=blue!50](0,0){0.1}\r
+\uput{0.3}[r](M012){$(t_0)$}\r
+\uput[d](0,0){$M_2(0)$}\r
+\psdot(M012)\uput{0.3}[l](M012){$M_{1/2}$}\r
+\rput(M012){\psline[unit=5,linecolor=red]{->}(!v0x1 v0x2 sub v0y1 v0y2 sub)}\r
+%\rput(G0){\psline[unit=5]{->}(!vG0x vG0y)}\r
+\psline{<->}(3,0)(0,0)(0,3)\r
+\rput(-2,2){\psPrintValue[decimals=4]{periode}\hphantom{0000000}s}\r
+\end{pspicture}\r
+\end{center}\r
+\section{Animation du mouvement de $M_1$ par rapport à $M_2$}\r
+\begin{center}\r
+\pstVerb{/XY12 [(XY12.dat) run] def\r
+ }%\r
+\def\nFrames{199}% 199 images au lieu de 200 ?\r
+\begin{animateinline}[controls,loop,timeline=pb2corps.dat,%\r
+ begin={\begin{pspicture}(-3,-2)(8,10)},\r
+ end={\end{pspicture}}]{20}% 20 images/s\r
+\psframe*[linecolor={[cmyk]{1 1 0 0.7}}](-3,-2)(8,10)\r
+\psgrid[subgriddiv=0,gridcolor=white,griddots=10,gridlabels=0pt]%\r
+\listplot[linecolor=gray]{XY12 aload pop}\r
+\pscircle*[linecolor=blue!50](0,0){0.1}\r
+\newframe\r
+\multiframe{\nFrames}{i=0+10}{% 1 point sur 10\r
+\pstVerb{/X1 XY12 \i\space get def\r
+ /Y1 XY12 \i\space 1 add get def\r
+ }%\r
+\pscircle*[linecolor=red](!X1 Y1){0.3}\r
+}\r
+\end{animateinline}\r
+\end{center}\r
+\end{document}
\ No newline at end of file
Projets / pst-eqdf.git / commitdiff
Les fichiers de Syracuse sont mis à
disposition (sauf mention contraire) selon les termes de la