MPSI

On aura une approche calculatoire. Introduction et objectifs Une équation différentielle est une équation dans laquelle une fonction inconnue et une ou plusieurs de ses dérivés apparait. On appelle équation différentielle linéaire une équation différentielle de la forme an(x)y(n)+…+a0(x)y=b(x)a_n(x)y^{(n)}+\ldots+a_0(x)y=b(x)an​(x)y(n)+…+a0​(x)y=b(x) ces équations différentielles sont mieux connues, on verra ici des solutions à des équations du premier (n=1n=1n=1) et deuxième ordre à coefficients (ana_nan​) constants et second membre (b(x)b(x)b(x)) particulier (de la forme P(x)ekxP(x)e^{kx}P(x)ekx).

On aura une approche technique du sujet. Théorème fondamental du calcul intégral Soient f∈C(I,C)f\in\mathcal{C}(I,\mathbb{C})f∈C(I,C) et (a,b)∈I2(a,b)\in I^2(a,b)∈I2. Alors la fonction c↦∫axf(t)dtc\mapsto \int_a^x f(t)dtc↦∫ax​f(t)dt est une primitive de fff sur III, souvent noté FFF. Pour toute primitive FFF de fff sur III, ∫abf(t)dt=∫abf=F]ab=F(b)−F(a)\int_a^b f(t)dt=\int_a^b f=F]^b_a=F(b)-F(a)∫ab​f(t)dt=∫ab​f=F]ab​=F(b)−F(a). Primitives usuelles Ces primitives usuelles sont à connaitre. Dans ce qui suit, uuu et vvv désignent deux fonctions d’une variable réelle xxx, et c,kc,kc,k deux constantes réelle.

Questions préliminaires Script : def subdivision1(tmax:float, n:int) -> list: Dt = tmax/n # Pas de temps return [k*Dt for k in range(n+1)] print(subdivision1(1., 5)) Renvoie : [0.0, 0.2, 0.4, 0.6, 0.8, 1.0] Script : def subdivision2(tmax:float, Dt:float) -> list: n = int(tmax/Dt) return [k*Dt for k in range(n+1)] print(subdivision2(0.3, 0.07)) Renvoie : [0.0, 0.07, 0.14, 0.21000000000000002, 0.28] (Ignorez l’erreur de python) def resolution(e, U:float, tau:float, tmax:float, Dt:float) -> (list, list): t = subdivision2(tmax, Dt) u = [U] uk = U r = Dt/tau for k in range(len(t)-1): uk = uk + r * (e(t[k])-uk) u.