{"nbformat":4,"nbformat_minor":0,"metadata":{"colab":{"provenance":[{"file_id":"1zWlxCU-oN6lGKFwGUSYod6q51zY1cXHu","timestamp":1633599075728}],"authorship_tag":"ABX9TyPrFGL6EzSA05s/JW4yUuGf"},"kernelspec":{"name":"python3","display_name":"Python 3"},"language_info":{"name":"python"}},"cells":[{"cell_type":"code","metadata":{"id":"JC3OsAwZL6oY"},"source":["import matplotlib.pyplot as plt #Pour tracer des graphiques\n","import numpy as np #Pour faire divers calculs\n","\n","## Données expérimentales 2021\n","\n","t=np.array([20,40,60,80,100,120,140,160,180,200,220,240]) # valeur du temps en s\n","A1=np.array(['à remplir']) # valeur de l'absorbance expérience 1 à ??°C\n","A2=np.array(['à remplir']) # valeur de l'absorbance expérience 2 à ??°C\n","A3=np.array(['à remplir']) # valeur de l'absorbance expérience 3 à ??°C\n","\n","Tol=.02 #Tolérance du spectro\n","\n","# Régression linéaire avec np.polyfit\n","\n","p1=np.polyfit(t, A1, 1) #Régression linéaire de A1 en fonction de t (équation A1 = p1[0]*t+p1[1])\n","p2=np.polyfit(t, A2, 1) #Régression linéaire de A2 en fonction de t (équation A2 = p2[0]*t+p2[1])\n","p3=np.polyfit(t, A3, 1) #Régression linéaire de A3 en fonction de t (équation A3 = p3[0]*t+p3[1])\n","\n","coef1=np. corrcoef(t,A1,1)  # calcul du coefficient de corrélation\n","r1=coef1[0,1] # attribue le coefficient de corrélation à r\n","r12=r1**2 # calcul coefficient de détermination\n","\n","coef2=np. corrcoef(t,A2,1)  # calcul du coefficient de corrélation\n","r2=coef2[0,1] # attribue le coefficient de corrélation à r\n","r22=r2**2 # calcul coefficient de détermination\n","\n","coef3=np. corrcoef(t,A3,1)  # calcul du coefficient de corrélation\n","r3=coef3[0,1] # attribue le coefficient de corrélation à r\n","r32=r3**2 # calcul coefficient de détermination\n","\n","# Graphes et droites de régression\n","\n","plt.figure(1)\n","plt.xlabel('temps en s') #Légende de l’axe des abscisses\n","plt.ylabel('A') #Légende de l’axe des ordonnées\n","\n","plt.plot(t,A1,'b+',label='A1') #Représente A1 en fonction de t avec des croix bleues avec une étiquette pour la courbe\n","plt.plot(t,A2,'g+',label='A2') #Représente A2 en fonction de t avec des croix vertes avec une étiquette pour la courbe\n","plt.plot(t,A3,'r+',label='A3') #Représente A3 en fonction de t avec des croix rouges avec une étiquette pour la courbe\n","\n","plt.plot(t,np.polyval(p1,t),'b',label='modèle') #Tracé de la droite de régression en bleu\n","plt.plot(t,np.polyval(p2,t),'g',label='modèle') #Tracé de la droite de régression en vert\n","plt.plot(t,np.polyval(p3,t),'r',label='modèle') #Tracé de la droite de régression en rouge\n","\n","plt.title('suivi cinétique de l\\'iodation de la propanone') #Titre du graphe\n","plt.legend() #Affiche l’étiquette de la courbe\n","plt.grid() #Affiche le quadrillage\n","plt.show() #Affiche le graphique\n","\n","plt.figure(2)\n","plt.xlabel('temps en s') #Légende de l’axe des abscisses\n","plt.ylabel('A') #Légende de l’axe des ordonnées\n","\n","plt.errorbar(t, A1, yerr=2*Tol*A1, fmt='.b',label='barre d\\'incertitude élargie') #Ajoute les barres d’incertitude élargie en bleu (ici 2 fois .02*A)\n","plt.errorbar(t, A2, yerr=2*Tol*A2, fmt='.g',label='barre d\\'incertitude élargie') #Ajoute les barres d’incertitude élargie en vert (ici 2 fois .02*A)\n","plt.errorbar(t, A3, yerr=2*Tol*A3, fmt='.r',label='barre d\\'incertitude élargie') #Ajoute les barres d’incertitude élargie en rouge (ici 2 fois .02*A)\n","\n","plt.plot(t,np.polyval(p1,t),'b',label='modèle') #Tracé de la droite de régression en bleu\n","plt.plot(t,np.polyval(p2,t),'g',label='modèle') #Tracé de la droite de régression en vert\n","plt.plot(t,np.polyval(p3,t),'r',label='modèle') #Tracé de la droite de régression en rouge\n","\n","plt.title('suivi cinétique de l\\'iodation de la propanone') #Titre du graphe\n","plt.legend() #Affiche l’étiquette de la courbe\n","plt.grid() #Affiche le quadrillage\n","plt.show() #Affiche le graphique\n","\n","print('pente 1=',format(p1[0],\"#.2e\"),'/s') #Affiche la pente de la droite expérience 1\n","print('y0_1 =', format(p1[1],\"#.3g\")) #Affiche l’ordonnée à l’origine expérience 1\n","print(' r^2 = ' ,format(r12,\"#.4f\"))\t#Affiche le coefficient r**2\n","\n","print('pente 2=',format(p2[0],\"#.2e\"),'/s') #Affiche la pente de la droite expérience 2\n","print('y0_2 =', format(p2[1],\"#.3g\")) #Affiche l’ordonnée à l’origine expérience 2\n","print(' r^2 = ' ,format(r22,\"#.4f\"))\t#Affiche le coefficient r**2\n","\n","print('pente 3=',format(p3[0],\"#.2e\"),'/s') #Affiche la pente de la droite expérience 3\n","print('y0_3 =', format(p3[1],\"#.3g\")) #Affiche l’ordonnée à l’origine expérience 3\n","print(' r^2 = ' ,format(r32,\"#.4f\"))\t#Affiche le coefficient r**2"],"execution_count":null,"outputs":[]}]}