Cours 1/13 (mercredi 4/9)
graphe fonction, courbe paramétrée
Présentation du module.
1er chapitre Courbes paramétrées:
Motivation: graphes de fonction insuffisant, cinématique
Définition de courbe paramétrée sur un intervalle ou une réunion d'intervalles de R. Cas particulier t=x pour les fonctions.
Rappel pour les fonctions: 1/ domaine de définition, 2/ domaine d'étude: restriction par parité ou/et périodicité, 3/ étude aux bornes de l'intervalle, asymptote, 4/ recherche des variations avec f', 5/ raffinement optionnel: convexité, 6/ tableau de variations 7/ tracé. Exemple avec sqrt(x^4+1)/x
Vérification ou calcul à la machine, tracé à la machine. Discrétisation: attention! plotfunc(x+0.01/(x-sqrt(2)),x=-2..2,xstep=0.1) et xstep=0.01
Généralisation aux courbes paramétrées:
1/ Domaine,
2/ domaine d'étude, restriction du domaine d'étude par parité (4 cas), par périodicité (par ex. pour x=cos(t), y=sin(3t))
3/ Etude des branches infinies si x ou y -> infini ou les 2. Asymptote horizontale/verticale si un seul tend vers l'infini. Recherche d'asymptote oblique si x et y tendent vers l'infini.
mat307 2024/25
Modérateur : xcasadmin
Re: mat307 2024/25
cours 2/13 vendredi 6 septembre
Exemple x(t)=2t+1/(2t+1), y(t)=t^2-1/(2t+1), recherche des asymptotes obliques en t=-+inf (branche para) et t=-1/2 x+y=-3/4
Etude locale en t0
La tangente est portée par la vitesse si elle est non nulle/point régulier.
Si x'=0 ou y'=0 tangente verticale ou horizontale.
Def point singulier. Si l'accélération est non nulle la tangente est portée par l'accélération pour un point singulier et on a un rebroussement.
Exemple ci-dessus accélération (-8,10)
Double tableau de variations et tracé de l'exemple. Il y a un point d'inflexion t>0 (apres la tangente verticale).
Raffinement: étude de la convexité: Si (vitesse,accélération) est un repère direct la courbe est convexe,
si indirect concave, si nul inflexion analytique, si nul en changeant de signe, inflexion géométrique.
Critère équivalent m=y'/x' pente de la tangente en un point régulier, si le signe de m' est +, la pente croit, convexe, si - concave, si nul inflexion analytique.
Exemple: les calculs sont vite pénibles, on les fait à la machine
Exemple x(t)=2t+1/(2t+1), y(t)=t^2-1/(2t+1), recherche des asymptotes obliques en t=-+inf (branche para) et t=-1/2 x+y=-3/4
Etude locale en t0
La tangente est portée par la vitesse si elle est non nulle/point régulier.
Si x'=0 ou y'=0 tangente verticale ou horizontale.
Def point singulier. Si l'accélération est non nulle la tangente est portée par l'accélération pour un point singulier et on a un rebroussement.
Exemple ci-dessus accélération (-8,10)
Double tableau de variations et tracé de l'exemple. Il y a un point d'inflexion t>0 (apres la tangente verticale).
Raffinement: étude de la convexité: Si (vitesse,accélération) est un repère direct la courbe est convexe,
si indirect concave, si nul inflexion analytique, si nul en changeant de signe, inflexion géométrique.
Critère équivalent m=y'/x' pente de la tangente en un point régulier, si le signe de m' est +, la pente croit, convexe, si - concave, si nul inflexion analytique.
Exemple: les calculs sont vite pénibles, on les fait à la machine
Code : Tout sélectionner
X:=2t+1/(2t+1); Y=t^2-1/(2t+1);
X1:=diff(X,t); Y1:=diff(Y,t);
factor(X1*diff(Y,t,2)-diff(X,t,2)*Y1);
X(t=1/2); Y(t=1/2). X1(t=1/2); Y1(t=1/2);
Re: mat307 2024/25
cours 3/13 du 12 septembre
Position par rapport à la tangente: étude en fonction de la parité de p et q, premières dérivées formant un repère. Exemple X:=2t+1/(2t+1); Y:=t^2-1/(2t+1) en t=-1: D.L. en posant t=-1+h, rebroussement de 1ère espèce.
Courbes en polaires:
définition, attention r(theta) peut être négatif
plan d'étude comme en paramétriques: domaine de définition, d'étude avec les symétries/périodicité.
Exemple avec 1/(1+2*cos(theta)), symétrie Ox et domaine 0,pi sauf en 2*pi/3
Branches infinies: si theta=+/-inf est dans le domaine d'étude on peut avoir un cercle ou un point asymptote. En un theta0 fini, si r tend vers l'infini, calcul de limite r*sin(theta-theta0)=l, si l est fini asymptote Y=l dans le repère tourné de theta0. Si l=infini on a une branche parabolique de direction asymptotique faisant un angle theta0.
Exemple avec 1/(1+2*cos(theta))
Tableau de variations
courbes en polaire
Position par rapport à la tangente: étude en fonction de la parité de p et q, premières dérivées formant un repère. Exemple X:=2t+1/(2t+1); Y:=t^2-1/(2t+1) en t=-1: D.L. en posant t=-1+h, rebroussement de 1ère espèce.
Courbes en polaires:
définition, attention r(theta) peut être négatif
plan d'étude comme en paramétriques: domaine de définition, d'étude avec les symétries/périodicité.
Exemple avec 1/(1+2*cos(theta)), symétrie Ox et domaine 0,pi sauf en 2*pi/3
Branches infinies: si theta=+/-inf est dans le domaine d'étude on peut avoir un cercle ou un point asymptote. En un theta0 fini, si r tend vers l'infini, calcul de limite r*sin(theta-theta0)=l, si l est fini asymptote Y=l dans le repère tourné de theta0. Si l=infini on a une branche parabolique de direction asymptotique faisant un angle theta0.
Exemple avec 1/(1+2*cos(theta))
Tableau de variations
courbes en polaire
Re: mat307 2024/25
Cours du 19/9 (4/13)
Courbes en polaires (fin): étude locale
Calcul de la "vitesse" r'*e_r+r*e_theta
Attention vitesse cinématique, multiplier par d theta/dt
Point régulier tan(V)=r/r', V=(e_r,tangente), pi/2 si r'=0
Point singulier si r=r'=0, en l'origine uniquement. Tangente portée par e_r. Rebroussement génériquement, cas général selon la parité de la 1ère dérivée de r qui s'annule.
Inflexion analytique: 1/r+(1/r)''=0
Exemple 1/(1+2*cos(theta)))
Tutoriel utilisation des calculatrices (1ère partie): allumer/éteindre/touche MENU et lancement de Xcas. Menus rapides F1 à F3, accès à l'aide, menu complet avec F4, touche F5 pour bloquer en mode alphabétique minuscule, touche ALPHA pour taper un caractère alphabétique majuscule. Exemple du tracé en polaire de 1/(1+2*cos(theta)).
Complément sur les coniques, définies comme r(theta)=A/(1+e*cos(theta)), e=excentricité, e=2 exemple étudié précédemment, e=0 cercle.
Si e<1, alors défini partout, tangente verticales en theta=0 et pi, tracé.
Equation cartésienne : c'est un cercle en x,Y si on pose Y=y/sqrt(1-e^2), demi-grand axe a=A/(1-e^2), dont le centre est décalé en abscisse de l'origine de a*e.
Lien avec la physique : lois de Kepler. Pour la Terre, e=0.0167 et theta=0 vers le 3 ou le 4 janvier. Au début de l'hiver (Nord), la Terre reçoit plus d'énergie du Soleil qu'au début de l'été, environ 6% (4*e).
En climatologie: un mot sur l'explication astronomique des glaciations du quaternaire (théorie de Milankovitch)
N.B.: les résultats de ce complément ne sont pas exigibles en examen.
Courbes en polaires (fin): étude locale
Calcul de la "vitesse" r'*e_r+r*e_theta
Attention vitesse cinématique, multiplier par d theta/dt
Point régulier tan(V)=r/r', V=(e_r,tangente), pi/2 si r'=0
Point singulier si r=r'=0, en l'origine uniquement. Tangente portée par e_r. Rebroussement génériquement, cas général selon la parité de la 1ère dérivée de r qui s'annule.
Inflexion analytique: 1/r+(1/r)''=0
Exemple 1/(1+2*cos(theta)))
Tutoriel utilisation des calculatrices (1ère partie): allumer/éteindre/touche MENU et lancement de Xcas. Menus rapides F1 à F3, accès à l'aide, menu complet avec F4, touche F5 pour bloquer en mode alphabétique minuscule, touche ALPHA pour taper un caractère alphabétique majuscule. Exemple du tracé en polaire de 1/(1+2*cos(theta)).
Complément sur les coniques, définies comme r(theta)=A/(1+e*cos(theta)), e=excentricité, e=2 exemple étudié précédemment, e=0 cercle.
Si e<1, alors défini partout, tangente verticales en theta=0 et pi, tracé.
Equation cartésienne : c'est un cercle en x,Y si on pose Y=y/sqrt(1-e^2), demi-grand axe a=A/(1-e^2), dont le centre est décalé en abscisse de l'origine de a*e.
Lien avec la physique : lois de Kepler. Pour la Terre, e=0.0167 et theta=0 vers le 3 ou le 4 janvier. Au début de l'hiver (Nord), la Terre reçoit plus d'énergie du Soleil qu'au début de l'été, environ 6% (4*e).
En climatologie: un mot sur l'explication astronomique des glaciations du quaternaire (théorie de Milankovitch)
N.B.: les résultats de ce complément ne sont pas exigibles en examen.
Re: mat307 2024/25
Cours 5 sur 13 du 26 septembre
Propriétés métriques des courbes (début):
Longueur d'arc paramétré (C1 par morceaux)=intégrale de ||vitesse instantanée|| * dt
=intégrale sqrt(x'(t)^2+y'(t)^2)*dt ou intégrale sqrt(r'^2+r^2) *d theta
Exemple: segment, arc de cercle, arc de parabole (t,t^2/2)
Remarque: il n'y a pas souvent de formule explicite -> calcul approché à la machine en donnant une des bornes sous forme d'un nombre approché (1.0 au lieu de 1 par exemple): integrate(sqrt(1+t^2),t,0,1.0)
Remarque: ne dépend pas du paramétrage choisi
Paramétrage par s, abscisse curviligne. Vecteur vitesse=vitesse scalaire*vecteur tangent, vecteur tangent=dM/ds.
Calcul de l'accélération d2M/dt^2, formule a=dv/dt*T+v^2*dT/ds
dT/ds est orthogonal a T, dT/ds= kappa * N, kappa la courbure et N le vecteur normal tel que (M,T,N) soit orthonormé direct, c'est le repère de Frenet.
Formules pour le calcul de la courbure et repère de Frenet en paramétriques et en polaire
Tutoriel d'utilisation de la calculatrice: vérification de résultats de l'étude analytique de x(t)=2t+1/(2t+1), y(t)=t^2-1/(2t+1). En utilisant plotparam, la calculatrice définit des variables x0,x1,x2,y0,y1,y2 contenant x(t) et ses dérivées et y(t) et ses dérivées, ce qui permet facilement de vérifier les points singuliers ou les points d'inflexion solve(x1*y2-x2*y1=0)
Commandes utiles: tabvar (menu F4 Analyse), courbure (F4 Courbes)
Propriétés métriques des courbes (début):
Longueur d'arc paramétré (C1 par morceaux)=intégrale de ||vitesse instantanée|| * dt
=intégrale sqrt(x'(t)^2+y'(t)^2)*dt ou intégrale sqrt(r'^2+r^2) *d theta
Exemple: segment, arc de cercle, arc de parabole (t,t^2/2)
Remarque: il n'y a pas souvent de formule explicite -> calcul approché à la machine en donnant une des bornes sous forme d'un nombre approché (1.0 au lieu de 1 par exemple): integrate(sqrt(1+t^2),t,0,1.0)
Remarque: ne dépend pas du paramétrage choisi
Paramétrage par s, abscisse curviligne. Vecteur vitesse=vitesse scalaire*vecteur tangent, vecteur tangent=dM/ds.
Calcul de l'accélération d2M/dt^2, formule a=dv/dt*T+v^2*dT/ds
dT/ds est orthogonal a T, dT/ds= kappa * N, kappa la courbure et N le vecteur normal tel que (M,T,N) soit orthonormé direct, c'est le repère de Frenet.
Formules pour le calcul de la courbure et repère de Frenet en paramétriques et en polaire
Tutoriel d'utilisation de la calculatrice: vérification de résultats de l'étude analytique de x(t)=2t+1/(2t+1), y(t)=t^2-1/(2t+1). En utilisant plotparam, la calculatrice définit des variables x0,x1,x2,y0,y1,y2 contenant x(t) et ses dérivées et y(t) et ses dérivées, ce qui permet facilement de vérifier les points singuliers ou les points d'inflexion solve(x1*y2-x2*y1=0)
Commandes utiles: tabvar (menu F4 Analyse), courbure (F4 Courbes)
Re: mat307 2024/25
Cours 6 sur 13 du 5/10:
Rappel repère de Frenet, courbure. Formule kappa=d/ds(angle de T avec un axe fixe), formule dN/ds
Recherche d'un cercle au plus proche de la courbe: centre sur la normale, même accélération normale, -> Omega=M+1/kappa*N
Exemple de calcul complet courbure/centre cercle osculateur pour la parabole (t,t^2/2).
Illustration machine du fait que génériquement la courbe traverse le cercle osculateur.
Les centres des cercles osculateurs décrivent une autre courbe, la développée.
Calcul de la tangente en un point de la développée, c'est la normale à la courbe de départ
Illustration du fait que la développée est l'enveloppe des normales à la courbe.
https://www-fourier.univ-grenoble-alpes ... evpara.xws
Si kappa est extrémal, la développée admet un point singulier au centre du cercle osculateur correspondant.
Définition de sommet d'une courbe, point ou la courbure est extrémale, c'est l'analogue à l'ordre 2 de point d'inflexion a l'ordre 1. Exemple de la parabole.
Remarque: équation intrinsèque d'une courbe: recherche d'une courbe vérifiant une relation entre le rayon de courbure et la longueur d'arc s.
Motivation: on veut parcourir une courbe à vitesse constante avec une accélération normale augmentant linéairement entre 0 et une valeur fixée (celle qu'on a sur un arc de cercle) pour faire un raccord entre une ligne droite et un arc de cercle.
v^2/R=s donc d/ds(theta)=cste*s où theta est l'angle de la tangente avec une direction fixe => on intègre puis T=(cos(theta),sin(theta))=dM/ds
Tracé à la machine:
plotparam(integrate(exp(i*s^2),s,0,t),t,0,10)
2ème partie: formes différentielles et intégrales curvilignes (intersection non vide avec le cours de mat304)
Dérivée directionnelle d'une fonction, dV(w)=dérivée directionnelle de V selon w. Calcul en coordonnées dV((a,b))=a*partial_x V+b*partial_y V. Exemple V(x,y)=x [respectivement y], dx(w)=1ère [resp. dy(w)=2ème] composante de w.
On a alors dV=partial_x V dx+ partial_y V y
Rappel repère de Frenet, courbure. Formule kappa=d/ds(angle de T avec un axe fixe), formule dN/ds
Recherche d'un cercle au plus proche de la courbe: centre sur la normale, même accélération normale, -> Omega=M+1/kappa*N
Exemple de calcul complet courbure/centre cercle osculateur pour la parabole (t,t^2/2).
Illustration machine du fait que génériquement la courbe traverse le cercle osculateur.
Les centres des cercles osculateurs décrivent une autre courbe, la développée.
Calcul de la tangente en un point de la développée, c'est la normale à la courbe de départ
Illustration du fait que la développée est l'enveloppe des normales à la courbe.
https://www-fourier.univ-grenoble-alpes ... evpara.xws
Si kappa est extrémal, la développée admet un point singulier au centre du cercle osculateur correspondant.
Définition de sommet d'une courbe, point ou la courbure est extrémale, c'est l'analogue à l'ordre 2 de point d'inflexion a l'ordre 1. Exemple de la parabole.
Remarque: équation intrinsèque d'une courbe: recherche d'une courbe vérifiant une relation entre le rayon de courbure et la longueur d'arc s.
Motivation: on veut parcourir une courbe à vitesse constante avec une accélération normale augmentant linéairement entre 0 et une valeur fixée (celle qu'on a sur un arc de cercle) pour faire un raccord entre une ligne droite et un arc de cercle.
v^2/R=s donc d/ds(theta)=cste*s où theta est l'angle de la tangente avec une direction fixe => on intègre puis T=(cos(theta),sin(theta))=dM/ds
Tracé à la machine:
plotparam(integrate(exp(i*s^2),s,0,t),t,0,10)
2ème partie: formes différentielles et intégrales curvilignes (intersection non vide avec le cours de mat304)
Dérivée directionnelle d'une fonction, dV(w)=dérivée directionnelle de V selon w. Calcul en coordonnées dV((a,b))=a*partial_x V+b*partial_y V. Exemple V(x,y)=x [respectivement y], dx(w)=1ère [resp. dy(w)=2ème] composante de w.
On a alors dV=partial_x V dx+ partial_y V y
Re: mat307 2024/25
Exemples
calcul de gradient: V(M)=d(O,M), gradient=e_r
calcul de gradient en coordonnées polaire et lien avec la dérivée directionnelle
Prop: Les courbes de niveau de V sont orthogonales au gradient.
Exemple: cercle.
Cas de l'ellipse: le gradient est porté par la bissectrice intérieure des demi-droites issues des 2 foyers.
Intégrale curviligne d'une forme différentielle omega le long d'un arc paramétré gamma(t), t \in [a,b]
Exemple arc de parabole gamma(t)=(t,t^2), t \in [0,1] et omega=y*dx
Autre paramétrage (u^2,u^4), u \in [0,1]. Même résultat.
Ceci se généralise, justifiant la définition, il suffit de faire un changement de variables dans l'intégrale.
Attention, si on change de sens de parcours, on change de signe.
Peut-on calculer l'intégrale par différence entre une "primitive" aux 2 extrémités comme en dimension 1?
Si omage=dV, alors oui
On appelle forme diff exacte une forme pour lequel c'est possible.
Si omega est exacte alors omega=dV, idée de preuve en intégrant sur le segment (x,y)->(x+h,y) on trouve M=derivee partielle de V selon x
Remarque: convention de signe opposée à la physique avec F=-grad(V)
Condition nécessaire pour être exacte (si V est C^2): partial M/partial y=partial N/partial x. On parle de forme fermée.
On a exacte -> fermée, mais la réciproque n'est pas toujours vraie, s'il y a des trous dans le domaine de définition de omega (admis pour l'instant).
Si fermée, recherche de potentiel en intégrant une des 2 équations puis en remplaçant dans l'autre.
Exemples la prochaine fois.
calcul de gradient: V(M)=d(O,M), gradient=e_r
calcul de gradient en coordonnées polaire et lien avec la dérivée directionnelle
Prop: Les courbes de niveau de V sont orthogonales au gradient.
Exemple: cercle.
Cas de l'ellipse: le gradient est porté par la bissectrice intérieure des demi-droites issues des 2 foyers.
Intégrale curviligne d'une forme différentielle omega le long d'un arc paramétré gamma(t), t \in [a,b]
Exemple arc de parabole gamma(t)=(t,t^2), t \in [0,1] et omega=y*dx
Autre paramétrage (u^2,u^4), u \in [0,1]. Même résultat.
Ceci se généralise, justifiant la définition, il suffit de faire un changement de variables dans l'intégrale.
Attention, si on change de sens de parcours, on change de signe.
Peut-on calculer l'intégrale par différence entre une "primitive" aux 2 extrémités comme en dimension 1?
Si omage=dV, alors oui
On appelle forme diff exacte une forme pour lequel c'est possible.
Si omega est exacte alors omega=dV, idée de preuve en intégrant sur le segment (x,y)->(x+h,y) on trouve M=derivee partielle de V selon x
Remarque: convention de signe opposée à la physique avec F=-grad(V)
Condition nécessaire pour être exacte (si V est C^2): partial M/partial y=partial N/partial x. On parle de forme fermée.
On a exacte -> fermée, mais la réciproque n'est pas toujours vraie, s'il y a des trous dans le domaine de définition de omega (admis pour l'instant).
Si fermée, recherche de potentiel en intégrant une des 2 équations puis en remplaçant dans l'autre.
Exemples la prochaine fois.
Re: mat307 2024/25
Intégrales curvilignes, potentiel:
Exemple 1: y*dx non fermée donc non exacte, calcul sur l'arc de parabole déjà fait, calcul sur le segment reliant (0,0) a (1,1), résultat différent
Exemple 2: y*dx+x*dy
Fermée, exacte; calcul du potentiel V
calcul de 2 manières de l'intégrale curviligne de omega sur l'arc de parabole reliant (0,0) a (1,1): par V(B)-V(A) ou en paramétrant avec (t,t^2) ou (t,t), vérification que c'est la même chose
Exemple 3: (y*dx-x*dy)/(x^2+y^2) fermée mais pas exacte à cause de (0,0). Calcul de l'intégrale sur le cercle unité
Intégrales curvilignes sur un lacet et intégrale double
Théorème de Stokes/Green-Riemann.
Idée de démonstration: rectangle puis réunion de rectangles et passage à la limite.
Application au calcul d'intégrales doubles: recherche de M et N, ils ne sont pas uniques, en général on prend un des deux nul.
calcul d'aire avec x*dy ou -y*dx ou 1/2*(-y*dx+x*dy), exemple aire de l'ellipse,
centre d'inertie d'un quart de cercle.
Tout début des équations et systèmes différentiels: ordre, dimension, problème de Cauchy, existence d'une solution unique passant par une condition initiale donnée si la fonction f(t,y) est régulière.
Exemple 1: y*dx non fermée donc non exacte, calcul sur l'arc de parabole déjà fait, calcul sur le segment reliant (0,0) a (1,1), résultat différent
Exemple 2: y*dx+x*dy
Fermée, exacte; calcul du potentiel V
calcul de 2 manières de l'intégrale curviligne de omega sur l'arc de parabole reliant (0,0) a (1,1): par V(B)-V(A) ou en paramétrant avec (t,t^2) ou (t,t), vérification que c'est la même chose
Exemple 3: (y*dx-x*dy)/(x^2+y^2) fermée mais pas exacte à cause de (0,0). Calcul de l'intégrale sur le cercle unité
Intégrales curvilignes sur un lacet et intégrale double
Théorème de Stokes/Green-Riemann.
Idée de démonstration: rectangle puis réunion de rectangles et passage à la limite.
Application au calcul d'intégrales doubles: recherche de M et N, ils ne sont pas uniques, en général on prend un des deux nul.
calcul d'aire avec x*dy ou -y*dx ou 1/2*(-y*dx+x*dy), exemple aire de l'ellipse,
centre d'inertie d'un quart de cercle.
Tout début des équations et systèmes différentiels: ordre, dimension, problème de Cauchy, existence d'une solution unique passant par une condition initiale donnée si la fonction f(t,y) est régulière.
Re: mat307 2024/25
7/11, cours 9 sur 13
Théorème de Cauchy-Lipschitz (admis): si f continument dérivable sur RxR^n il existe une solution maximale sur un intervalle ouvert en temps I passant par une condition initiale donnée.
Valable aussi sur ]t0,t1[x]y0,y1[
Deux solutions ne peuvent pas se croiser.
Ne dit rien sur l'intervalle en temps, qui peut être R, ]-inf,t1[, ]t0,+inf[, ]t0,t1[. Tout dépend si y tend vers l'infini en temps fini.
Exemple 1: y'=t*y^2, on a la solution nulle, les autres sont non nulles, donc on peut diviser par y^2 pour trouver les autres solutions.
Exemple 2: (Contrôle d'une solution par des solutions connues, pas de croisement), par exemple y'=y*(1-y), si 0<y(t0)<1 la solution y reste, y est croissante et majorée -> y converge en +infini, y'->0 donc y tend vers 1. allure du graphe d'une solution
Exemple 3: t*y'=y, attention parait régulier, mais il faut exprimer y' en fonction de t et y, y'=y/t et CL ne s'applique pas en 0 on a une infinité de solutions passant par (0,0) (y=C*t) et pas de solution pour y(0)!=0.
Illustration machine
session
Idée très vague de preuve: équation intégrale équivalente -> méthodes d'approximation numérique généralisant le champ des tangentes, suite de fonctions dont il faut montrer la convergence.
2/ Méthodes de résolution explicite
-> Variables séparables y'=g(t)*h(y),
solutions particulières stationnaires: y=racine de h,
sinon
H(y)=G(t)+Cte, G primitive de g et H de 1/h
Cas favorable: on peut calculer y en fonction de t
cas général: équation implicite entre t et y, à tracer dans le plan des (t,y)
Exemple : résolution de y'=t*y^2,
3/ Cas particulier équations linéaires du 1er ordre homogène ( sans second membre):
-> : y=K*exp(H(t)), K reel.
Exemple : y'=t*y
4/ Equations linéaires du 1er ordre avec second membre: méthode de variation de la constant, exemple y'=t*y-t
Théorème de Cauchy-Lipschitz (admis): si f continument dérivable sur RxR^n il existe une solution maximale sur un intervalle ouvert en temps I passant par une condition initiale donnée.
Valable aussi sur ]t0,t1[x]y0,y1[
Deux solutions ne peuvent pas se croiser.
Ne dit rien sur l'intervalle en temps, qui peut être R, ]-inf,t1[, ]t0,+inf[, ]t0,t1[. Tout dépend si y tend vers l'infini en temps fini.
Exemple 1: y'=t*y^2, on a la solution nulle, les autres sont non nulles, donc on peut diviser par y^2 pour trouver les autres solutions.
Exemple 2: (Contrôle d'une solution par des solutions connues, pas de croisement), par exemple y'=y*(1-y), si 0<y(t0)<1 la solution y reste, y est croissante et majorée -> y converge en +infini, y'->0 donc y tend vers 1. allure du graphe d'une solution
Exemple 3: t*y'=y, attention parait régulier, mais il faut exprimer y' en fonction de t et y, y'=y/t et CL ne s'applique pas en 0 on a une infinité de solutions passant par (0,0) (y=C*t) et pas de solution pour y(0)!=0.
Illustration machine
session
Idée très vague de preuve: équation intégrale équivalente -> méthodes d'approximation numérique généralisant le champ des tangentes, suite de fonctions dont il faut montrer la convergence.
2/ Méthodes de résolution explicite
-> Variables séparables y'=g(t)*h(y),
solutions particulières stationnaires: y=racine de h,
sinon
H(y)=G(t)+Cte, G primitive de g et H de 1/h
Cas favorable: on peut calculer y en fonction de t
cas général: équation implicite entre t et y, à tracer dans le plan des (t,y)
Exemple : résolution de y'=t*y^2,
3/ Cas particulier équations linéaires du 1er ordre homogène ( sans second membre):
-> : y=K*exp(H(t)), K reel.
Exemple : y'=t*y
4/ Equations linéaires du 1er ordre avec second membre: méthode de variation de la constant, exemple y'=t*y-t
Re: mat307 2024/25
Cours 10/13 (14/11)
Equation linéaire d'ordre n homogène -> espace vectoriel de solutions, on montre que c'est de dimension n en appliquant Cauchy-Lipschitz sur le système associé. Pas de méthode générale de résolution si n>1.
Si on connait une solution s(t), on peut se ramener à résoudre une équation d'ordre n-1 par une méthode de type variation de la constante, en posant y=K(t)s(t), z=K' va vérifier une équation d'ordre n-1. Exemple t^2/2*y''=y avec s(t)=t^2 => t*z'+4*z=0 => z=Ct^(-4) => K=C1 t^(-3) +C2
Avec second membre: solution particulière+solution générale homogène, comme dans le cas à coeffs constants.
Linéaire à coefficient constant,
Homogène: équation caractéristique P(r)=0.
Prop: Si P a n racines simples r_1, .., r_j, ..., r_n l'espace vectoriel des solutions est engendré par les e^{r_j*t}
Exemple 1: y''+3y'-4y=0,
Exemple 2: y''+2y'+5y=0, et passage au réel
Exemple 3: y''+4y'+4y=0, pas assez de solutions
Idée de la démonstration
Prop: Si r racine multiple, remplacer "e^{r*t} m fois" par e^{r*t},t*e^{r*t}, ..., t^{m-1}*e^{r*t}
Exemple: y''+4y'+4y=0, vérification, généralisation: cf poly.
Résolution avec second membre: solution particulière+solution générale.
Méthode de variation des constantes pour trouver une solution particulière,
exemple dans un cas d'ordre 2: y''+3y'-4y=exp(t)
y=lambda_1*exp(-4t)+lambda_2*exp(t)
résolution de système linéaire en les lambda_j' en imposant lambda_1'*exp(-4t)+lambda_2'*exp(t)=0, puis en remplaçant on a lambda_1'*-4exp(-4t)+lambda_2'*exp(t)=second membre
Cas général: cf poly.
Les calculs sont assez pénibles, si le second membre est du type G(t)*exp(r*t), solution particulière de la même forme exp(r*t)*A(t) avec A de même degré que G + multiplicité de r dans le polynôme caractéristique, ce qui donne des calculs plus simples. Cas de l'exemple.
Equation linéaire d'ordre n homogène -> espace vectoriel de solutions, on montre que c'est de dimension n en appliquant Cauchy-Lipschitz sur le système associé. Pas de méthode générale de résolution si n>1.
Si on connait une solution s(t), on peut se ramener à résoudre une équation d'ordre n-1 par une méthode de type variation de la constante, en posant y=K(t)s(t), z=K' va vérifier une équation d'ordre n-1. Exemple t^2/2*y''=y avec s(t)=t^2 => t*z'+4*z=0 => z=Ct^(-4) => K=C1 t^(-3) +C2
Avec second membre: solution particulière+solution générale homogène, comme dans le cas à coeffs constants.
Linéaire à coefficient constant,
Homogène: équation caractéristique P(r)=0.
Prop: Si P a n racines simples r_1, .., r_j, ..., r_n l'espace vectoriel des solutions est engendré par les e^{r_j*t}
Exemple 1: y''+3y'-4y=0,
Exemple 2: y''+2y'+5y=0, et passage au réel
Exemple 3: y''+4y'+4y=0, pas assez de solutions
Idée de la démonstration
Prop: Si r racine multiple, remplacer "e^{r*t} m fois" par e^{r*t},t*e^{r*t}, ..., t^{m-1}*e^{r*t}
Exemple: y''+4y'+4y=0, vérification, généralisation: cf poly.
Résolution avec second membre: solution particulière+solution générale.
Méthode de variation des constantes pour trouver une solution particulière,
exemple dans un cas d'ordre 2: y''+3y'-4y=exp(t)
y=lambda_1*exp(-4t)+lambda_2*exp(t)
résolution de système linéaire en les lambda_j' en imposant lambda_1'*exp(-4t)+lambda_2'*exp(t)=0, puis en remplaçant on a lambda_1'*-4exp(-4t)+lambda_2'*exp(t)=second membre
Cas général: cf poly.
Les calculs sont assez pénibles, si le second membre est du type G(t)*exp(r*t), solution particulière de la même forme exp(r*t)*A(t) avec A de même degré que G + multiplicité de r dans le polynôme caractéristique, ce qui donne des calculs plus simples. Cas de l'exemple.
Re: mat307 2024/25
Cours 11 sur 13
EDO linéaires à coefficients constants.
Principe de superposition : application au calcul de solution particulière pour un second membre trigonométrique en utilisant les complexes.
Exemples: y'+y=cos(w*t), calcul de l'amplitude et du déphasage de la solution particulière.
y''+2y'+5y=cos(w*t), idem
y''+y=cos(w*t), discussion en fonction de la valeur de w
Comportement asymptotique des solutions d'équations différentielles
cas 1: linéaires à coeff constants en dimension 1
Cas homogène:
Si toutes les racines de l'équation caractéristique sont de partie réelle < 0 alors convergence vers 0
Si <=0 et si les racines de partie réelle sont simples, les solutions sont bornées.
Exemples: y''+y=0, y''+y'+y=0, y''+3y'+2y=0, y''+4y'+4y=0 allure des solutions
Régime permanent: toutes les solutions s'approchent d'une solution particulière si parties réelles des racines<0
Cas d'une partie réelle nulle (simple): fonction bornée sauf si i*omega est solution de l'équation caractéristique -> résonance.
EDO linéaires à coefficients constants.
Principe de superposition : application au calcul de solution particulière pour un second membre trigonométrique en utilisant les complexes.
Exemples: y'+y=cos(w*t), calcul de l'amplitude et du déphasage de la solution particulière.
y''+2y'+5y=cos(w*t), idem
y''+y=cos(w*t), discussion en fonction de la valeur de w
Comportement asymptotique des solutions d'équations différentielles
cas 1: linéaires à coeff constants en dimension 1
Cas homogène:
Si toutes les racines de l'équation caractéristique sont de partie réelle < 0 alors convergence vers 0
Si <=0 et si les racines de partie réelle sont simples, les solutions sont bornées.
Exemples: y''+y=0, y''+y'+y=0, y''+3y'+2y=0, y''+4y'+4y=0 allure des solutions
Régime permanent: toutes les solutions s'approchent d'une solution particulière si parties réelles des racines<0
Cas d'une partie réelle nulle (simple): fonction bornée sauf si i*omega est solution de l'équation caractéristique -> résonance.
Re: mat307 2024/25
Cours du 28 novembre (12 sur 13)
Comportement des solutions d'équations autonomes y'=f(y): équilibre stable si f passe de positif à négatif en y=r, instable dans le cas contraire. Si f'(r)<0 stable, si f'(r)>0 instable.
Exemple y'=y0^4-y^4 (cf. l'exam de janvier 2023)
Remarque: linéarisé de l'équation en un point d'équilibre, permet d'avoir une idée du comportement de y'=f(y)+g(t) près d'un point d'équilibre de y'=f(y) pendant un court intervalle de temps.
Autres méthodes de résolution explicite d'équations différentielles: constantes du mouvement:
* y'=f(x,y). S'il existe g telle que omega=g(dy-fdx) est exacte, les courbes integrales sont courbes de niveau du potentiel de omega. Exemple y'=-y/(x+y). Si omega n'est pas fermée, on peut chercher à la rendre fermée en utilisant un facteur intégrant.
* système conservatif en dimension 1 permet de remplacer m*x''=somme des forces en
1/2*m*x'^2+V(x)=E, équation à variables séparables,
Exemple: calcul de la période du pendule sans approximation, fait intervenir int(1/sqrt(cos(theta)-cos(theta0)),theta=0..theta0), on trouve numériquement une période 4% plus grande avec un angle initial theta0=pi/4
session Xcas
Systèmes différentiels linéaires à coefficients constants:
Y'=A*Y homogène, Y dans R^d, ou Y'=A*Y+b(t) avec second membre
Dans le cas diagonalisable
A=P*D*P^(-1), D diagonale, Y'=A*Y+b, Y=P*Z, Z vérifie un système diagonal Z'=D*Z+P^-1*b, résolution, puis Y=P*Z.
Exemple: A:=[[1,-1],[2,4]] résolution avec second membre b=[exp(t),0]
Si n=2, A=[[a,b],[c,d]], Y=[x,y], on peut se ramener à une équation d'ordre 2 en x si b!=0 par substitution. On observe que son équation caractéristique est le polynome caractéristique de la matrice A. Ceci permet de traiter les cas où A n'est pas diagonalisable.
Comportement des solutions d'équations autonomes y'=f(y): équilibre stable si f passe de positif à négatif en y=r, instable dans le cas contraire. Si f'(r)<0 stable, si f'(r)>0 instable.
Exemple y'=y0^4-y^4 (cf. l'exam de janvier 2023)
Remarque: linéarisé de l'équation en un point d'équilibre, permet d'avoir une idée du comportement de y'=f(y)+g(t) près d'un point d'équilibre de y'=f(y) pendant un court intervalle de temps.
Autres méthodes de résolution explicite d'équations différentielles: constantes du mouvement:
* y'=f(x,y). S'il existe g telle que omega=g(dy-fdx) est exacte, les courbes integrales sont courbes de niveau du potentiel de omega. Exemple y'=-y/(x+y). Si omega n'est pas fermée, on peut chercher à la rendre fermée en utilisant un facteur intégrant.
* système conservatif en dimension 1 permet de remplacer m*x''=somme des forces en
1/2*m*x'^2+V(x)=E, équation à variables séparables,
Exemple: calcul de la période du pendule sans approximation, fait intervenir int(1/sqrt(cos(theta)-cos(theta0)),theta=0..theta0), on trouve numériquement une période 4% plus grande avec un angle initial theta0=pi/4
session Xcas
Systèmes différentiels linéaires à coefficients constants:
Y'=A*Y homogène, Y dans R^d, ou Y'=A*Y+b(t) avec second membre
Dans le cas diagonalisable
A=P*D*P^(-1), D diagonale, Y'=A*Y+b, Y=P*Z, Z vérifie un système diagonal Z'=D*Z+P^-1*b, résolution, puis Y=P*Z.
Exemple: A:=[[1,-1],[2,4]] résolution avec second membre b=[exp(t),0]
Si n=2, A=[[a,b],[c,d]], Y=[x,y], on peut se ramener à une équation d'ordre 2 en x si b!=0 par substitution. On observe que son équation caractéristique est le polynome caractéristique de la matrice A. Ceci permet de traiter les cas où A n'est pas diagonalisable.
Re: mat307 2024/25
jeudi 5 décembre, cours 13 sur 13
Systèmes différentiel linéaires:
Vérification de calculs avec la calculatrice: shift-2 (Mat) saisie de matrice, recherche valeurs/vecteurs propres, desolve pour un système.
Allure des courbes en dimension 2 (en lien avec l'exo 3.6 du TP).
Comportement à l'infini des systèmes homogènes: comme pour les équas diff en remplaçant racine de l'équation caractéristique par valeur propre, si parties réelles < 0 convergent vers 0, si <=0 et diagonalisable, solutions bornées.
Introduction au calcul variationnel:
Lagrangien L(x,x',t) (attention j'ai noté x point en cours, pas x' mais sur ce forum je ne sais pas faire un x point...), action S, exemple: chemin le plus court en dimension 2, lagrangien de la mécanique classique=énergie cinétique-potentiel
Equations d'Euler-Lagrange pour minimiser S sur les chemins d'extrémités fixées vus comme courbes paramétrées par le temps avec paramètres fixé aux extrémités.
Exemple:
chemin le plus court: on trouve que la dérivée du vecteur tangent est nulle
lagrangien de la mécanique classique en cartésien : les équations d'Euler-Lagrange redonnent m*accélération=forces
Exemple en polaire pour le problème à 2 corps, conservation du moment cinétique
Cas où L ne dépend pas explicitement du temps: conservation du hamiltonien associé au lagrangien L
Systèmes différentiel linéaires:
Vérification de calculs avec la calculatrice: shift-2 (Mat) saisie de matrice, recherche valeurs/vecteurs propres, desolve pour un système.
Allure des courbes en dimension 2 (en lien avec l'exo 3.6 du TP).
Comportement à l'infini des systèmes homogènes: comme pour les équas diff en remplaçant racine de l'équation caractéristique par valeur propre, si parties réelles < 0 convergent vers 0, si <=0 et diagonalisable, solutions bornées.
Introduction au calcul variationnel:
Lagrangien L(x,x',t) (attention j'ai noté x point en cours, pas x' mais sur ce forum je ne sais pas faire un x point...), action S, exemple: chemin le plus court en dimension 2, lagrangien de la mécanique classique=énergie cinétique-potentiel
Equations d'Euler-Lagrange pour minimiser S sur les chemins d'extrémités fixées vus comme courbes paramétrées par le temps avec paramètres fixé aux extrémités.
Exemple:
chemin le plus court: on trouve que la dérivée du vecteur tangent est nulle
lagrangien de la mécanique classique en cartésien : les équations d'Euler-Lagrange redonnent m*accélération=forces
Exemple en polaire pour le problème à 2 corps, conservation du moment cinétique
Cas où L ne dépend pas explicitement du temps: conservation du hamiltonien associé au lagrangien L