MONTAGE SEUWIN

Expérience 1

In [1]:
# Importation des bibliothèques

import numpy as np
import matplotlib.pyplot as plt
from scipy import stats
In [2]:
#--------------------------------------------------------------------------------------------------------------------------------
# Déclaration grandeurs
#--------------------------------------------------------------------------------------------------------------------------------

# Points présentation

X_pres1 = 1.5
Y_pres1 = 2.5

X_pres2 = 3.2
Y_pres2 = 4.2


X_pres = np.array([X_pres1,X_pres2])
Y_pres = np.array([Y_pres1,Y_pres2])


#Ensemble des points

X = np.array([0,1,2,3,4,5,6,7,8,9,10,X_pres1,X_pres2])
Y = np.array([1,2,3,4,5,6,7,8,9,10,11,Y_pres1,Y_pres2])


# Incertitudes

u_X = 0.1
u_Y = 0.5
In [3]:
#--------------------------------------------------------------------------------------------------------------------------------
# Ajustement
#--------------------------------------------------------------------------------------------------------------------------------

reg = stats.linregress(X,Y)
yfit = reg[0]*X + reg[1]
u_b = reg[1]*reg[4]/reg[1]

print("le coefficient directeur est" + u" %3.2e "% (reg[0]) +"+/-" + u" %5.2e "% (reg[4]))
print("l'ordonnée à l'origine est" + u" %3.2e "% (reg[1]) +"+/-" + u" %5.2e "% (u_b))
le coefficient directeur est 1.00e+00 +/- 0.00e+00 
l'ordonnée à l'origine est 1.00e+00 +/- 0.00e+00 
In [4]:
#--------------------------------------------------------------------------------------------------------------------------------
# Tracé
#--------------------------------------------------------------------------------------------------------------------------------

plt.figure(figsize=(20,12))
plt.errorbar(X,Y,xerr = u_X, yerr = u_Y, marker = '+', linestyle = ' ',color = 'b', capsize =6, label = 'points expérimentaux')
plt.plot(X,yfit,'g', label = 'ajustement')
plt.errorbar(X_pres,Y_pres, xerr = u_X, yerr = u_Y, marker = '+', linestyle = ' ', color = 'r', capsize = 6, label = 'points présentation')
plt.xlabel('X (unité de X)',fontsize = 20)
plt.ylabel('Y (unité de Y)',fontsize = 20)
plt.title('X en fonction de Y',fontsize = 20)
plt.legend(loc=0, fontsize = 20)
plt.show()