From: Juergen Gilg Date: Sat, 23 Jun 2012 17:09:02 +0000 (+0200) Subject: Ajouté L'animation du satellite ... X-Git-Url: https://melusine.eu.org/syracuse/G/git/?a=commitdiff_plain;h=35b2d821966345a0f542207db4cf8b9cb97fd112;p=pst-eqdf.git Ajouté L'animation du satellite ... --- diff --git a/gravitation/LISTE.txt b/gravitation/LISTE.txt index c229127..6f7257e 100644 --- a/gravitation/LISTE.txt +++ b/gravitation/LISTE.txt @@ -1,20 +1,34 @@ -== Satellite +== Satellite + eqdf-grav_01.tex + eqdf-grav_01.pdf + gravitation_01.tex -gravitation_01.pdf + +gravitation_01.pdf + gravitation_02.tex + gravitation_02.pdf -== Exemples avec Distiller + +== Animation du satellite +animation_satellite.tex +animation_satellite.pdf +== Exemples avec Adobe Distiller gravitation_01_distiller.tex + gravitation_01_distiller.pdf -== Les fichiers du Paquet -pst-eqdf.sty +== Les fichiers du paquet pst-eqdf.tex + +pst-eqdf.sty + + pstricks-add.tex == Les figures g1.png g2.png g3.png h1.png -h2.png \ No newline at end of file +h2.png diff --git a/gravitation/animation_satellite.pdf b/gravitation/animation_satellite.pdf new file mode 100644 index 0000000..503a600 Binary files /dev/null and b/gravitation/animation_satellite.pdf differ diff --git a/gravitation/animation_satellite.tex b/gravitation/animation_satellite.tex new file mode 100644 index 0000000..931d118 --- /dev/null +++ b/gravitation/animation_satellite.tex @@ -0,0 +1,156 @@ +\documentclass{article} +\usepackage[a4paper,margin=2cm]{geometry} +\usepackage[latin1]{inputenc} +\usepackage{pstricks,pst-eqdf} +\usepackage{animate} +\newpsstyle{vecteurA}{arrowinset=0.05,arrowsize=0.125,linecolor={[rgb]{1 0.5 0}}} +\newpsstyle{vecteurB}{arrowinset=0.05,arrowsize=0.1,linecolor={[rgb]{0 0.7 1}}} +\newpsstyle{vecteurC}{arrowinset=0.1,arrowsize=0.2,linecolor={[rgb]{1 0.5 0}}} +%%%%%%%%%%%%%%%%%% +\title{Animation du mouvement d'un satellite} +\date{23 juin 2\,012} +%timeline +\begin{filecontents}{satellite.dat} +::0x0 +::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 +::27 +::28 +::29 +::30 +::31 +::32 +::33 +::34 +::35 +::36 +::37 +::38 +::39 +::40 +::41 +::42 +::43 +::44 +::45 +::46 +::47 +::48 +::49 +::50 +::51 +::52 +\end{filecontents} +\newcommand{\Jupiter}{% + \psclip{\pscircle[fillstyle=solid,fillcolor=yellow]{1.5}} + \psset{fillstyle=solid,fillcolor={[cmyk]{0 0.2 0.4 0}},linestyle=none} + \psframe(-2,.6)(2,1) + \psframe(-2,0.45)(2,0.52) + \psframe(-2,.05)(2,.3) + \psframe(-2,-0.35)(2,-0.2) + \psframe(-2,-0.9)(2,-.5) + \endpsclip + \psellipse[fillstyle=solid,fillcolor={[cmyk]{0 0.4 0.6 0}},linestyle=none](-0.3,-0.6)(0.35,0.2) + \pscircle{1.5}} +\begin{document} +\def\eqsatellite{% +y[2]|y[3]|-GM*y[0]/((sqrt(y[0]^2+y[1]^2))^3)|-GM*y[1]/((sqrt(y[0]^2+y[1]^2))^3)}% +\begin{center} +%\begin{pspicture}(-10,-10)(6,6) +\def\nFrames{50}% 50 images +\begin{animateinline}[controls,loop,timeline=satellite.dat,% + begin={\begin{pspicture}(-10,-10)(6,6)}, + end={\end{pspicture}}]{5}% 5 images/s +%\uput[r](0,3){$y$}\uput[u](3,0){$x$} +%\psset{unit=2} +\psframe*[linecolor={[rgb]{0 0 0.5}}](-10,-10)(6,6) +\pstVerb{/tabXiYi [(XiYi.dat) run] def}% +\pstVerb{/GM 1 def + /theta0 30 def + /r0 2 def + /x0 r0 theta0 cos mul def + /y0 r0 theta0 sin mul def + /v0 0.85 def + /v0x v0 theta0 sin mul neg def + /v0y v0 theta0 cos mul def + /Lc r0 v0 mul def % moment cinetique + /par Lc dup mul GM div def % paramètre de l'ellipse +% excentricité + /exc 1 0.5 v0 4 exp mul r0 dup mul mul GM r0 mul v0 dup mul mul sub GM dup mul div 2 mul add sqrt def +%%%%%%%%%%%%%% + /a_2 par 1 exc dup mul sub div def % demi-grand axe + /b_2 par 1 exc dup mul sub sqrt div def % demi-petit axe + /periode 2 3.1416 dup mul a_2 3 exp mul GM div sqrt mul def}% +\psequadiff[method=rk4, + plotpoints=1000, + algebraic, + whichabs=0, + whichord=1, + tabname=XiYi, +% saveData,filename=XiYi.dat +]{0}{43}{x0 y0 v0x v0y}{\eqsatellite}% +\listplot[linecolor=gray,unit=2]{XiYi aload pop} +\psequadiff[method=rk4, + plotpoints=1000, + algebraic, + whichabs=2, + whichord=3, + tabname=vxvy, +% saveData,filename=vxvy.dat +]{0}{43}{x0 y0 v0x v0y}{\eqsatellite}% +%\listplot[unit=2,linecolor=red]{vxvy aload pop} +% on dessine la vitesse un point sur 50 +%\pscircle[fillcolor=gray!70,fillstyle=solid,unit=2](0,0){0.4} +\Jupiter% +%\psgrid[subgriddiv=2,gridcolor=lightgray,gridlabels=8pt,unit=2](-5,-5)(3,3) +%\psline[arrowinset=0.05,arrowsize=0.1,unit=2]{<->}(3,0)(0,0)(0,3) +\newframe +\multiframe{\nFrames}{i=0+40}{ +\psset{unit=2}% +\pstVerb{/tabXYpartiel { + 2 2 \i\space {/I exch def + tabXiYi I get + tabXiYi I 1 add get + } for + } def + /vX vxvy \i\space get def + /vY vxvy \i\space 1 add get def + /xi tabXiYi \i\space get def + /yi tabXiYi \i\space 1 add get def + /ri3 xi dup mul yi dup mul add sqrt 3 exp def + /Fx xi ri3 div neg def + /Fy yi ri3 div neg def}% +\listplot[linecolor=red,linewidth=0.05]{tabXYpartiel} +\rput(!xi yi){\pscircle[fillstyle=solid,fillcolor=gray!50]{0.075}\psline[style=vecteurA]{->}(! vX 2 mul vY 2 mul)\psline[style=vecteurB]{->}(! Fx 5 mul Fy 5 mul)}} +\newframe +\listplot[linecolor=gray,unit=2]{XiYi aload pop} +\listplot[linecolor=red,unit=2,linewidth=0.05]{XiYi aload pop} +\psset{unit=2}% +\rput(!x0 y0){\pscircle[fillstyle=solid,fillcolor=gray!50]{0.075}\psline[style=vecteurA]{->}(! v0x 2 mul v0y 2 mul)\psline[style=vecteurB]{->}(! x0 r0 3 exp div 5 mul neg y0 r0 3 exp div 5 mul neg)} +\end{animateinline} +\end{center} +\end{document}