Exercice 1 :
cm = 8.630e-3 * 40e3
print("cm = " + str(cm) + " m/s")
On cherche l'incertitude type $u(cm)$ qui s'exprime comme le produit de deux grandeurs, on peut écrire :
$$\frac{u(c_m)}{c_m} = \sqrt{\left(\frac{u(\lambda_m)}{\lambda_m}\right)^2 + \left(\frac{u(f_0)}{f_0}\right)^2}$$Avec un intervalle de confiance à 95% : $u(\lambda_m) = \dfrac{\Delta\lambda_m}{2}\quad$ et $\quad u(f_0) = \dfrac{\Delta f_0}{2}$
Donc
$$u(c_m) = c_m \sqrt{\left(\frac{\Delta(\lambda_m)}{2\lambda_m}\right)^2 + \left(\frac{\Delta(f_0)}{2f_0}\right)^2}$$import math
u_cm = cm * math.sqrt((0.009/8.63)**2 + (0.025/40)**2)
print("u_cm = " + str(u_cm) + " m/s")
Avec un intervalle de confiance à 95% : $\Delta c_m = 2\,u(c_m)$ donc :
$$c = 345,2 \pm 0.8\,\text{m/s}$$Exercice 2 :
Si nous appliquons la relation vue en cours à la relation précédente que l'on peut exprimer sous la forme :
$$ q = f(x, y)\quad\quad \text{avec les mesures $x$ et $y$ indépendantes}$$On obtient :
$$u^2(q) = \left(\dfrac{\partial f}{\partial x}u(x)\right)^2 + \left(\dfrac{\partial f}{\partial x}u(y)\right)^2$$Il nous faut maintenant exprimer la dérivée partielle pour chaque variable :
$$\dfrac{\partial f}{\partial x} = \dfrac{2y}{(x+y)^2}$$
$$\dfrac{\partial f}{\partial y} = -\dfrac{2x}{(x+y)^2}$$
Il n'y a plus qu'à écrire notre incertitude type sur $q$:
$$u^2(q) = \left(\dfrac{2y}{(x+y)^2}\right)^2u^2(x) + \left(-\dfrac{2x}{(x+y)^2}\right)^2u^2(y)$$
$$u(q) = \dfrac{2}{(x+y)^2} \sqrt{y^2 u^2(y) + x^2 u^2(x)}\quad\quad (2)$$
D'après (1) on peut écrire
En remplaçant dans (2) on obtient :
Exercice 3 :
mesures = [538.2, 554.3, 545.7, 552.3, 566.4, 537.9, 549.2, 540.3]
def moyenne(mesures):
somme = 0
for elt in mesures:
somme += elt
return somme/len(mesures)
print("moyenne =",moyenne(mesures))
from math import sqrt
def ecart_type(mesures):
m = moyenne(mesures)
somme = 0
for elt in mesures:
somme += (elt-m)**2
return sqrt((1./(len(mesures)-1))*somme)
print("ecart_type =", ecart_type(mesures))
def incertitude_type(Et, n):
return Et/sqrt(n)
deltaLambda = 2.37*incertitude_type(ecart_type(mesures), len(mesures))
print("deltaLambda",deltaLambda)
Exercice 4 :
from sympy import *
init_printing(use_unicode=True) #affichage à la LaTeX
x = Symbol("x")
$h_1$ est définie et dérivable sur $\bigcup_{k\in\mathbb{Z}}=\left] -\dfrac{\pi}{6}+2k\pi\, ;\, \dfrac{7\pi}{6} + 2k\pi \right[$
h1 = log(1+sqrt(2*sin(x)+1))
La dérivée de $h_1$
diff(h1)
simplify(diff(h1))
$h_2$ est définie et dérivable sur $\mathbb{R}$
h2 = sin(log((exp(2*x)+1)/(exp(2*x)+3)))
La dérivée de $h_2$
diff(h2)
simplify(diff(h2))