Deux exemples récapitulatifs

Les deux exemples d'activités ci-dessous récapitulent l'essentiel des procédures abordées dans ce support de cours.

Les ressources sur le tir au basket

1
#!/usr/bin/env python3
2
# -*- coding: utf-8 -*-
3
"""
4
Created on Wed Oct 14 14:29:37 2020
5
6
@author: cedric
7
"""
8
9
import matplotlib.pyplot as plt
10
import numpy as np
11
from scipy.optimize import curve_fit        #permet de chercher une modélisation
12
13
# Lecture des donnees du fichier txt
14
t, x, y = np.loadtxt('tirbasket.txt', unpack = True , usecols=(0, 1, 2), delimiter = '\t', skiprows = 2)
15
16
# Creation du graphique avec les points de mesure
17
plt.axis('equal')
18
plt.title('Trajectoire du ballon')
19
plt.xlabel('Abscisse x (m)')
20
plt.ylabel('Altitude y (m)')
21
plt.scatter(x, y, color = 'blue', label = 'Trajectoire', marker = '+')
22
#plt.show()
23
24
25
#Définition du modèle que l'on souhaite vérifier
26
def FonctionModele(x,a,b,c):
27
    return a*x**2 + b*x + c
28
29
#Calcul du modèle
30
parametre,covariance= curve_fit(FonctionModele,x,y)
31
32
a=parametre[0]
33
b=parametre[1]
34
c=parametre[2]
35
36
37
#Tracé du modèle
38
xth=np.linspace(1,3.5,200)
39
yth=a*xth**2 + b*xth+c
40
plt.plot(xth, yth, color = 'green', label = 'Modele', marker = '')
41
42
#affichage de l'équation
43
plt.text(2,3,'y= {:.3f} $x^2$+{:.3f} x + {:.3f}'.format(a,b,c))
44
45
plt.show()
46

Ci-dessous, le fichier contenant les pointages exploités dans le programme.

Les ressources pour la loi de Newton

Description de l'expérience analysée

Fichiers ressources

1
#!/usr/bin/env python3
2
# -*- coding: utf-8 -*-
3
"""
4
Created on Wed Oct 14 14:29:37 2020
5
6
@author: cedric
7
"""
8
9
import matplotlib.pyplot as plt
10
import numpy as np
11
from scipy.optimize import curve_fit        #permet de chercher une modélisation
12
13
14
# Lecture des donnees du fichier txt
15
t, T = np.loadtxt('newtonlaw2.txt', unpack = True , usecols=(0, 1), delimiter = ';', skiprows = 0)
16
17
# Creation du graphique avec les points de mesure
18
#plt.axis('equal')                               #Pour choisir sir les 2 axes ont le même échelle
19
plt.title('Temperature en fonction du temps')
20
plt.xlabel('temps t (s)')
21
plt.ylabel('température T(°C) ')
22
plt.scatter(t, T, color = 'blue', label = 'Points expérimentaux', marker = '+')
23
plt.legend()
24
#plt.show()
25
26
27
#Définition du modèle que l'on souhaite vérifier
28
def FonctionModele(t,Ta,Ti,tau):
29
        return Ta + (Ti-Ta)*np.exp(-t/tau)
30
    #math.exp(-t2/tau)
31
32
#Calcul du modèle
33
parametre,covariance= curve_fit(FonctionModele,t,T,p0=(30,60,2000))
34
35
Ta=parametre[0]
36
Ti=parametre[1]
37
tau=parametre[2]
38
39
40
#Tracé du modèle
41
tth=np.linspace(0,7000,5000)
42
Tth=Ta + (Ti - Ta)*np.exp(-tth/tau)
43
plt.plot(tth, Tth, color = 'green', label = 'Modele', marker = '')
44
plt.legend()
45
46
#affichage de l'équation
47
plt.text(1500,50,'T(t) = {: .1f} + ({: .1f} - {: .1f}) . exp(-t/{: .1f})'.format(Ta,Ti,Ta,tau))
48
49
plt.show()
50

Ci-dessous, le fichier contenant les mesures exploitées dans le programme.