TSpe-TP9.ipynb

{"metadata":{"kernelspec":{"name":"python3","display_name":"Python 3 (ipykernel)","language":"python"},"language_info":{"name":"python","version":"3.7.12","mimetype":"text/x-python","codemirror_mode":{"name":"ipython","version":3},"pygments_lexer":"ipython3","nbconvert_exporter":"python","file_extension":".py"}},"nbformat_minor":4,"nbformat":4,"cells":[{"cell_type":"markdown","source":"# Mouvement de rotation  \n\nLes séparateurs décimaux en langage python sont des points.  \nCependant nous autres français avons l'habitude d'utiliser des virgules ... Tous nos logiciels sont donc adaptés et utilisent des virgues comme séparateurs décimaux.  \n\nLa première étape consiste donc à remplacer toutes les virgules par des points, dans notre fichier de pointage.\n\n> Le programme ci-dessous remplacerait les **?** par des **!** - Modifier ce programme pour qu'il remplace les **,** par des **.** puis l'exécuter.","metadata":{}},{"cell_type":"code","source":"import matplotlib.pyplot as plt\nimport numpy as np\n\n# Remplacement des , par des . dans le fichier de données\nfmaster = open(\"pointage.txt\", \"rt\") #fichier maitre\nfout = open(\"pointagepoint.txt\", \"wt\") #fichier modifié\nfor line in fmaster:\n    fout.write(line.replace('?', '!'))\nfmaster.close()\nfout.close()","metadata":{"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"markdown","source":"Le programme ci-dessous permet alors de tracer la trajectoire de la pastille, ainsi que le vecteur vitesse et le vecteur accélération en chaque point.\n\n> Exécuter le programme.","metadata":{}},{"cell_type":"code","source":"# Lecture des données du fichier\nt,x,y = np.loadtxt('pointagepoint.txt', unpack=True, skiprows = 3)\n\n# Calcul de la vitesse, à partir de la position\nvx, vy = [], []\nfor i in range(1, len(x)-1) :\n    vx.append((x[i+1]-x[i-1])/(t[i+1]-t[i-1]))\n    vy.append((y[i+1]-y[i-1])/(t[i+1]-t[i-1]))\n\n# Calcul de l'accélération, à partir de la vitesse\nax, ay = [], []\nfor i in range(1, len(vx)-1) :\n    ax.append((vx[i+1]-vx[i-1])/(t[i+1]-t[i-1]))\n    ay.append((vy[i+1]-vy[i-1])/(t[i+1]-t[i-1]))\n\n# Creation du graphique\nplt.axis('equal')\nplt.title('Trajectoire de la pastille')\nplt.xlabel('x (m)')\nplt.ylabel('y (m)')\nplt.plot(x, y, color = 'green', marker = '+')\n\n# Tracé des vecteurs vitesse\nfor i in range(len(vx)) :\n    plt.quiver(x[i+1], y[i+1], vx[i], vy[i], angles = \"xy\", scale_units = \"xy\", scale = 5, color = \"red\")\n\n# Tracé des vecteurs accélération\nfor i in range(len(ax)) :\n    plt.quiver(x[i+2], y[i+2], ax[i], ay[i], angles = \"xy\", scale_units = \"xy\", scale = 50, color = \"blue\")\n\n# Fin de la creation du graphique\nplt.show()","metadata":{"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"code","source":"","metadata":{},"execution_count":null,"outputs":[]}]}