Courbes paramétriques et équations différentielles pour la physique (Mat307-ex237)Bernard.Parisse@univ-grenoble-alpes.fr2017, 2022 |
![](mat237001.png)
Remarque : la version HTML de ce cours est interactive,
elle contient de nombreuses
commandes Xcas que le lecteur peut exécuter avec ou sans modifications
depuis un navigateur compatible
(l’interactivité est optimisée pour Firefox). Deux versions
sont proposées, une utilisant www.mathjax.org
pour
un rendu plus fidèle des formules, mais qui nécessite un
temps de chargement plus long, l’autre sans. Ces fichiers HTML
ont été générés avec hevea.inria.fr
de Luc Maranget,
et le fork de Yannick Chevallier pour le support mathjax.
Table des matières
- Chapitre 1 Présentation du module
- Chapitre 2 Prérequis
- Chapitre 3 Courbes paramétriques et polaires
- Chapitre 4 Propriétés métriques des courbes.
- Chapitre 5 Formes différentielles et intégrales curvilignes
- Chapitre 6 Équations et systèmes différentiels.
- 6.1 Introduction et représentation graphique.
- 6.2 Existence et unicité
- 6.3 Quelques méthodes de résolution explicite.
- 6.3.1 Équations à variables séparables
- 6.3.2 Équations linéaires
- 6.3.3 Équations linéaires à coefficients constants
- 6.3.4 Systèmes différentiels linéaires à coefficients constants d’ordre 1.
- 6.3.5 Systèmes et équations
- 6.3.6 Allure des courbes en dimension 2.
- 6.3.7 Systèmes d’ordre plus grand que 1
- 6.3.8 Intégrales premières.
- 6.3.9 Le modèle proie-prédateur
- 6.3.10 Quelques autres méthodes
- 6.4 Comportement asymptotique des solutions
- Chapitre 7 Introduction au calcul variationnel
- Chapitre 8 Utilisation de la calculatrice.
Index
|
|
Chapitre 1 Présentation du module
Ce cours commence par l’étude des particularités des courbes paramétrées, en distinguant propriétés cinématiques (dépendant du paramétrage comme la vitesse, l’accélération) et propriétés géométriques d’une courbe (c’est-à-dire intrinsèques, indépendamment du paramétrage, par exemple la longueur, la courbure). La géométrie d’une courbe peut sembler un objet d’étude intéressant uniquement le mathématicien. Mais ce n’est pas le cas : les lois de Képler par exemple traduisent sous forme de contrainte géométrique les lois de la mécanique céleste, et la relativité, c’est s’extraire d’un système de coordonnées ou d’un paramétrage pour trouver des équations covariantes par changement de système et des propriétés invariantes.
On abordera aussi un peu calcul différentiel et intégral dans le plan, en particulier l’intégrale curviligne, notion utile pour calculer des quantités comme le travail d’une force, mais aussi comme outil de calcul pratique pour trouver une aire, un centre d’inertie...
Dans la deuxième moitié du cours, on s’intéresse aux équations différentielles. Ce thème est très lié au précédent puisque les solutions d’une équation différentielle décrivent une courbe, et certaines propriétés de ces équations (constantes du mouvement par exemple) imposent des contraintes géométriques aux solutions. On rappellera quelques méthodes de résolution explicite d’équations et on parlera un peu des propriétés qualitatives ou asymptotiques de solutions, particulièrement importantes lorsqu’on ne sait pas résoudre explicitement une équation différentielles.
On terminera par une petite introduction au calcul variationnel, une théorie très utile pour mettre en équation de phénomènes physiques de manière covariante (les équations d’Euler-Lagrange). Un prolongement de ce chapitre (non étudié ici) est le théorème de Noether qui permet de déduire des constantes du mouvement en fonction de symétries du problème (par exemple quantité de mouvement et invariance par translation d’espace, énergie pour le temps, moment cinétique pour les rotations).
Le programme en mots-clefs :
-
Étude des courbes paramétrées :
paramétrisation, cas des coniques, symétries, périodicité, branches infinies, étude locale (tangente, point singulier, convexité), propriétés métriques (longueur d’arc, courbure, accélération normale et tangentielle). Courbes en polaires, cas des coniques.
Lien avec la physique : la cinématique. - Intégration le long d’un chemin d’une forme différentielle
Lien avec la physique : travail d’une force, force dérivant d’un potentiel, forme différentielle fermée et thermodynamique. - Équations et systèmes différentiels :
champ des tangentes (principe de résolution approchée), existence et unicité d’une courbe intégrale passant par une condition initiale donnée, méthodes de résolution analytiques dans le cas à variables séparables et dans le cas linéaire (cas des coefficients constants, cas non constant, cas des systèmes, comportement asymptotique : régime permanent et transitoire), intégrales premières (constantes du mouvement), équations autonomes (points d’équilibre, linéarisation au voisinage).
Lien avec la physique : circuits RLC, radioactivité, cinématique, ... - introduction au calcul des variations : équations d’Euler-Lagrange. Lien avec la physique : principe de moindre action, formulation lagrangienne des équations de la mécanique.
Le but du module n’est pas de vous apprendre à dériver ou intégrer mécaniquement de grosses expressions mais de comprendre et mettre en oeuvre les notions au programme. Les calculs fastidieux ou techniques se feront directement sur ordinateur (en TP) ou sur calculatrices (en TD) avec un logiciel de calcul formel. Les calculs simples seront faits si possible en calcul mental (qui est un analogue de l’exercice physique pour le cerveau) sinon papier-crayon, à défaut sur calculatrice. On prendra l’habitude de vérifier systématiquement les calculs par divers moyens (calcul vérifié à la machine, cohérence avec la représentation graphique, cohérence du tableau de variations, signe et ordre de grandeur d’une aire, position d’un centre d’inertie, etc.). Les calculatrices pourront être autorisées au DS en fonction de la technicité des calculs du sujet et seront autorisées à l’examen final. Les étudiants n’ayant pas de calculatrices formelles peuvent en emprunter une pour le semestre.
L’évaluation se fait sur :
- 1/2 : DS semaine des partiels, TP et autre évaluation en TD
- 1/2 : l’examen final (avec règle de max)
Chapitre 2 Prérequis
Cette section liste quelques compétences dont l’expŕience montre
qu’elles ne sont pas suffisamment maitrisés par les étudiants, à
revoir avant le début des TD, cf. l’exercice 0 de la feuille d’exercices
et le lien suivant
maths.ac-orleans-tours.fr/ressources_lycee/activites_mentales/articles/fiches_de_calcul_mental_au_lycee/
2.1 Géométrie analytique
Droites du plan : savoir passer de l’équation cartésienne
à des équations
paramétriques et réciproquement,
et de même vers/de droite donnée par deux points, par un point
et un vecteur directeur. Représentation graphique. Coefficient directeur,
ordonnée à l’origine. Coefficient directeur de la tangente en un
graphe de fonction régulière .
Distance entre 2 points et (Pythagore),
distance d’un point à une droite d’équation :
Cercles : lien entre équation cartésienne et centre-rayon et
équation paramétrique trigonométrique.
Exemple d’exercice : calculer la tangente en un point du cercle ou
d’un graphe, l’intersection de la tangente avec l’axe des , etc.
2.2 Trigonométrie
Valeurs remarquables de (en ).
Passer de à et réciproquement.
Formules de duplication
Développement et linéarisation en utilisant les nombres complexes (ci-dessous désigne la partie réelle)
2.3 Complexes
Partie réelle, imaginaire et module, argument.
Conjugué .
2.3.1 Complexe et géométrie analytique
Affixe d’un point du plan, d’un vecteur.
Rotation de centre l’origine et d’angle revient à ajouter à l’argument. Translation revient à ajouter l’affixe du vecteur. Conjugué correspond à une symétrie par rapport à l’axe .
Chapitre 3 Courbes paramétriques et polaires
3.1 Introduction
Le graphe d’une fonction ( un intervalle) est un exemple de courbe du plan, mais il n’est pas assez général pour représenter tous les types de courbe du plan, par exemple un segment de droite vertical, ou un cercle, car deux points distincts d’un graphe doivent avoir des abscisses différentes. D’autre part, il apparait naturellement d’autres types de courbes que les graphes de fonction, par exemple la trajectoire d’un mobile dans le plan dont les coordonnées dépendent du temps (selon une équation différentielle ou un système différentiel), ce sont les courbes paramétriques, ou des courbes vérifiant une équation cartésienne également appelées courbes implicites, par exemple en géométrie le cercle , ou en cinématique des courbes de niveau de l’énergie totale dans le plan position-impulsion.
Dans cette section, on va étudier les courbes en paramétriques, donnée par un couple de fonctions définies pour dans un intervalle ou une réunion d’intervalles et à valeurs dans . (Ceci ne restreint pas trop la généralité, on peut montrer sous des hypothèses assez générales que l’allure locale d’une courbe implicite est identique à celle d’une courbe paramétrique, sauf en certains points dits singuliers, c’est le théorème des fonctions implicites).
Exemples :
-
le graphe d’une fonction est une courbe paramétré
d’équation ( est le temps).
On aurait aussi pu choisir
(ce qui revient à changer l’origine des
temps) ou d’autres paramétrages.
Exemple :
- une droite d’équation est le graphe d’une fonction, donc paramétrable comme ci-dessus. Une droite verticale peut se paramétrer par
- le cercle de centre l’origine et de rayon 1 peut se paramétrer
par . On peut bien sur
le paramétrer par , mais dans ce cas on parcourt plusieurs
fois le cercle (périodicité). On peut aussi paramétrer tout le
cercle sauf un point avec le paramétrage rationnel
. Ce
paramétrage permet de calculer plus facilement des points
du cercle, mais contrairement au paramétrage trigonométrique,
il n’est pas “uniforme”.
Exemple
ou avec les nombres complexes
3.2 Représentation graphique
La plupart des calculatrices graphiques
et de nombreux logiciels de maths permettent de
représenter graphiquement un arc de courbe en donnant des valeurs
extrêmes et (souvent notées tmin
et
tmax
) et un pas (tstep
). Le logiciel évalue
la valeur de et en , , , ... puis relie les points de la courbe obtenue par des segments
(parfois avec des autres arcs de courbes). La plupart du temps
cela donne une bonne idée de la courbe, mais parfois on peut
manquer un détail intéressant (valeur de trop
grande), ou un morceau de courbe (mauvaises
valeurs de et ). Par exemple
Il peut être nécessaire d’ajuster le cadrage graphique
à l’affichage (xmin
, xmax
, ymin
, ymax
) ou
de l’affiner avec un menu de zoom
. Attention, sur
certaines calculatrices,
les opérations de changement de cadrage graphique
provoquent un nouveau calcul complet qui peut durer une
dizaine de secondes.
Mise en oeuvre :
-
avec Xcas: on utilise la commande
plotparam
dans le menuGraphe, Courbes
. Le cadrage graphique est calculé automatiquement et peut être modifié par les touches menus à droite du graphe. On peut spécifier le pas avec l’argument optionneltstep=
. - HP Prime : touche Apps puis Paramétrique, entrer les
formules donnant et en fonction de , puis touche
Plot. Pour régler la discrétisation, faire shift-Plot (Plot
setup).
Ou touche Apps puis Géométrie et utiliser les mêmes commandes qu’avec Xcas. - CAS sur Casio Graph 90+e/35eii: commande
plotparam
du catalogue (F4, puis touche log 2 fois). - sur les calculatrices TI89/92/voyage 200 : il faut sélectionner le mode
de tracé paramétrique avec la touche MODE,
puis sur
Y=
, l’écran de définition de apparait. Si les réglages graphiques ne sont pas satisfaisants, la toucheWINDOW
permet d’y accéder, puis la toucheGRAPH
lance le tracé. - TI Nspire : lancer Graphes depuis Home puis touche Menu, 3. Graph puis 3. Parametric. On peut aussi régler le zoom depuis Menu.
Exemples : essayez de tracer quelques courbes en paramétriques
3.3 Paramétrage
On adoptera souvent la convention d’appeler temps le paramétre . Mais cela ne signifie pas que le paramétrage est réellement le temps mesuré en secondes. On peut très bien paramétrer une courbe avec un paramètre autre, qui peut être un multiple constant ou variable du temps (c’est d’ailleurs conforme au principe de la relativité). Le paramétrage n’est jamais unique, on peut changer de paramétrage pourvu que la fonction donnant le nouveau en fonction de l’ancien paramétrage soit une bijection (qui peut même renverser le sens de déroulement du temps c’est-à-dire le sens de parcours de la courbe). On utilisera d’ailleurs plus loin un paramétrage par la longueur, où la courbe est parcourue à vitesse constante égale à 1.
Le choix d’un paramétrage est ce qui fait la différence entre la cinématique (on prend le temps comme paramètre) et la géométrie (où on cherche à décrire les propriétés intrinséques de la courbe indépendamment du paramétrage). Ainsi, l’équation cartésienne d’une courbe est une propriété géométrique, indépendante du choix de paramétrage choisi pour l’obtenir.
On observe aussi que l’opération inverse, trouver un paramétrage à partir d’une équation cartésienne de courbe n’est pas possible de manière explicite, sauf dans quelques cas particuliers. C’est pour cette raison qu’il est beaucoup plus difficile (et couteux en temps) d’obtenir une représentation graphique d’une courbe donnée par son équation cartésienne.
3.4 Étude analytique d’une courbe en paramétrique
3.4.1 Rappel sur les graphes de fonction
Pour tracer le graphe d’une fonction , on commence par déterminer le domaine de définition, on restreint éventuellement l’intervalle d’étude (parité, périodicité). Puis on calcule les limites aux bornes du domaine de définition :
-
une limite infinie en fini correspond à une asymptote
verticale d’équation .
Exemple
- une limite finie lorsque correspond
à une asymptote horizontale d’équation .
Exemple
- si tend vers lorsque ,
on a une branche infinie. S’il s’agit d’une asymptote , alors
tend vers lorsque tend vers ce qui
permet de trouver en calculant ,
puis . Si est fini, on a une
asymptote, sinon une branche parabolique de direction .
On remarque que si ,
comme tend vers l’infini, on ne peut pas avoir d’asymptote, on
parle de branche parabolique de direction . Si tend vers
l’infini, on a une branche parabolique de direction .
Exemple
- Lorsqu’il y a une asymptote d’équation , on peut raffiner en cherchant la position de la courbe par rapport à l’asymptote en étudiant le signe de
On calcule la dérivée première de pour déterminer le sens de variations, les points d’annulation correspondent à des tangentes horizontales. On peut étudier la convexité de (signe de ), les points d’inflexion de la courbe se produisent lorsque s’annule. On trace alors la courbe en faisant apparaitre les points particuliers et les asymptotes.
Pour une courbe en paramétrique, le plan général est analogue, mais l’étude est un peu plus compliquée puisqu’on a deux fonctions et au lieu d’une seule .
3.4.2 Domaine et périodicité
On supposera dans toute la suite que les fonctions et sont continument dérivables au moins 2 fois, sauf peut-être en un nombre fini de réels d’un intervalle de .
On commence par déterminer le domaine de définition de
et de , et on essaie de le réduire si possible, soit par
périodicité (par exemple pour le cercle ci-dessus, )
soit par l’existence de symétries si les fonctions
et sont paires ou impaires.
Si et
sont paires, alors on parcourt deux fois le même arc de courbe
sur et , on peut restreindre le domaine d’étude
à .
Si est pair et impair, alors
, il y a une symétrie par rapport à
l’axe des , on se restreint à .
Par exemple, pour
Dans le cas périodique, on peut tester des symétries correspondant
à des demi (voire quart) de période.
3.4.3 Branches infinies
On s’intéresse ensuite aux bornes du domaine de définition et aux points où ou/et ne sont pas définis. Si et admettent une limite finie, on peut prolonger la courbe. Si les limites existent mais ne sont pas finies, on a une branche infinie ( ou ). Si l’une des deux valeurs tend vers l’infini, l’autre restant finie, on a une asymptote (horizontale si tend vers l’infini, verticale si tend vers l’infini), on peut déterminer la position de l’arc de courbe par rapport à l’asymptote en cherchant le signe de ou lorsque tend vers la valeur particulière (limite à droite et limite à gauche). Enfin si et tendent vers l’infini tous les deux, on cherche la limite de , Si , on a une branche parabolique de direction asymptotique , on cherche alors la limite , si cette limite est finie et vaut on a une asymptote oblique (on peut déterminer la position en cherchant le signe de .
Exemples :
On peut utiliser la commande limit
dans Xcas pour étudier
une asymptote, par exemple dans le premier cas, pour étudier la
branche infinie pour 1
On définit les fonctions
puis on calcule les limites lorsque
Si est fini non nul et fini, on en déduit que est asymptote à la courbe. Il y a une
autre asymptote pour ( fini, tend vers l’infini)
Autre exemple :
3.4.4 Étude locale
On se place en une valeur de où et sont continument dérivables au moins deux fois. On notera la dérivation par rapport au paramètre par le signe ’ (en physique on utilise aussi le point). On a alors un développement de Taylor à l’ordre 2 du vecteur où et sont compris entre et . Si le vecteur vitesse est non nul, on en déduit un équivalent Lorsque est proche de 0, le vecteur est équivalent à un vecteur colinéaire à (supposé non nul), qui est donc vecteur tangent à la courbe en .
On appelle point singulier un point où la vitesse s’annulle.
On verra dans la suite comment étudier la tangente en un point
singulier d’une courbe. Génériquement, une courbe n’a pas
de points singuliers, car il faut annuler simultanément les
deux dérivées, or on n’a qu’un seul paramètre libre . Par contre
une famille de courbes dépendant d’un paramètre
(par exemple )
possède en général un nombre discret de valeurs du
paramètre pour lesquelles la courbe admet un point singulier.
Dans l’exemple, , les deux
dérivées s’annulent si
(en , ) ou (en ).
Commandes Xcas :
x:=t^2-m*t; y:=m/(1+t^2)+t;
solve([diff(x,t),diff(y,t)],[m,t]);
supposons(m=[-2.0,-5,5,0.1]);
plotparam([x,y],t=((-3) .. 3),affichage=arrow_line);
Remarque : en cinématique, si la vitesse et l’accélération sont nulles en un point et que les équations ne dépendent pas explicitement du temps, on reste indéfiniment en ce point qui est un point d’équilibre, la notion de tangente à la courbe n’a alors pas de sens. On peut aussi suivre une trajectoire qui se rapproche de plus en plus d’un point d’équilibre (la limite de est alors ce point, pour si l’équilibre est stable ou si l’équilibre est instable).
Pour faire une étude locale plus précise dans le cas d’un point régulier, ou pour déterminer la tangente en un point singulier, il faut poursuivre le développement de Taylor à un ordre plus grand. À l’ordre 2, si et sont 3 fois continument dérivables, il existe tels que : Si les vecteurs vitesse et accélération ne sont pas colinéaires, alors forme une base, et dans cette base a pour coordonnées un terme d’ordre 3 en puissances de , l’arc de courbe est à l’ordre 2 identique à un arc de parabole. On parle de point birégulier. Si est une base directe, l’arc est convexe (la vitesse “tourne” dans le sens trigonométrique), sinon il est concave. On peut tester cela en calculant le déterminant des coordonnées de ou le sens de variations de , la pente de la tangente
Remarques :
- pour un graphe de fonction, , on retrouve le critère usuel .
- Attention, la convexité/concavité dépend du sens de paramétrage, cela n’a pas de sens de parler de convexité d’une courbe géométrique, seulement de changement de convexité/point d’inflexion.
Exemple : point d’inflexion en de
La courbe admet deux autres points d’inflexion ( et
) qu’on peut déterminer avec les commandes Xcas suivantes :
Note : on utilise comme paramètre
x
au lieu de t
pour pouvoir
utiliser la notation '
pour dériver (si on utilise t
comme paramètre, il faut utiliser diff(.,t)
pour calculer
la dérivée par rapport à ). L’instruction fsolve
effectue une résolution numérique, pour tenter une résolution
exacte, utiliser solve
, mais on risque alors de manquer
certaines solutions.On observe que la convexité est (presque) une propriété géométrique, en effet si on change de paramétrage on dérive par rapport à puis : on retrouve en facteur qui est positif si on parcourt la courbe dans le même sens ou négatif sinon.
La convexité décrit qualitativement la géométrie de la courbe (orientée) à l’ordre 1. On verra plus loin que le rayon de courbure décrit quantitativement la géométrie de la courbe à l’ordre 2 (comme la tangente décrit la géométrie de la courbe à l’ordre 1).
Dans le cas d’un point singulier (), si l’accélération , alors la tangente est portée par . L’étude complète de la nature d’un point singulier ou de la convexité d’un point régulier tel que est colinéaire à nécessite de faire un développement de Taylor en jusqu’au premier ordre , s’il existe, tel que :
- les dérivées d’ordre 1, .., de s’annulent
- la dérivée d’ordre est non nulle, on la note
- les dérivées d’ordre sont colinéaires à la dérivée d’ordre (ce qui inclus le cas où elles sont nulles)
- la dérivée d’ordre est non colinéaire à , on la note .
Dans la base , les composantes de sont alors respectivement équivalentes à et où . On en déduit que la tangente à la courbe est portée par .
- si est pair, on a un rebroussement de première espèce si est impair (cas générique d’un point singulier, ) ou de deuxième espèce si est pair. On ne peut pas régulariser le point singulier par changement de paramétrage.
- Si est impair, on peut reparamétriser la courbe pour rendre le point non singulier (prendre ) mais au risque de perdre de la régularité,
- Si est impair et impair on a un point d’inflexion géométrique (changement de sens de convexité).
Exemples de points singuliers en avec dans l’ordre rebroussement
de 1ère puis 2ième espèce, méplat et inflexion :
Les deux derniers cas peuvent être reparamétrés (au prix
de la perte de dérivabilité seconde) en posant .
Pour faire l’étude d’un point singulier avec Xcas, on peut utiliser
la fonction series
sur et (ici c’est inutile,
le développement de Taylor est déjà fait).
Remarque : il peut arriver dans des cas pathologiques que toutes les dérivées de s’annulent en un point sans que la fonction soit nulle (par exemple si et contiennent un facteur en , on parle de fonction plate). Il peut aussi arriver que toutes les dérivées soit colinéaires à la première dérivée non nulle, si on se déplace sur une droite ou si la tangeance est plate.
3.5 Plan d’étude d’une courbe
- On détermine et on restreint le domaine de définition (périodicité, symétries).
- On étudie les branches infinies (point exclus du domaine, ) : asymptotes horizontales, verticales, directions asymptotiques, asymptotes obliques.
- Recherche de et , on étudie l’annulation conjointe des deux (points singuliers).
- Signe de et , double tableau de variations faisant apparaitre et mise en évidence des tangentes horizontales et verticales
- Pour préciser le tracé, on peut chercher la convexité en étudiant le signe de .
- Tracé des points remarquables et des asymptotes et on les relie entre eux en suivant les sens de variations du tableau de variations.
3.6 Courbes en polaires
Une courbe en polaire est essentiellement donnée par la distance au centre d’un point de la courbe en fonction de l’angle entre la direction et le vecteur : On autorise toutefois des valeurs négatives pour , si c’est le cas, on prend alors le symétrique par rapport à l’origine du point situé à distance et d’angle .
Représentation graphique : avec Xcas
(ou CAS sur Casio Graph 90+e/35eii), on utilise
la commande plotpolar
(menu Graphes, Courbes dans Xcas ou
Repr. graphique sur Casio), sur calculatrices
graphiques, sélectionner le mode de tracé en polaire (touche Mode
sur TI89/92/V200) ou
l’application Polaire ou Géométrie sur les HP Prime.
Par exemple
Remarque : une courbe en polaires est toujours parcourue dans le sens trigonométrique.
C’est un cas particulier de courbe en paramétriques puisque mais on préfère souvent faire l’étude directement sur la fonction . Le plan d’étude est calqué sur celui d’une courbe en paramétrique, mais on n’a qu’une seule fonction à étudier.
- domaine de définition de , recherche de périodicités et symétries ( ou ajout d’une demi ou d’un quart de période). Si la période n’est pas un multiple de , cela correspond à obtenir un arc de la courbe par rotation à partir d’un autre arc de la courbe.
- branches infinies pour (non infini)
où n’est pas défini. La branche a pour direction
asymptotique la droite faisant un angle avec
l’axe des . On calcule alors la limite si elle existe de
, c’est l’ordonnée dans
le repère obtenu par rotation d’angle , si la limite
est finie et vaut
on a une asymptote (d’équation dans le repère tourné).
Exemple . n’est pas défini pour , donc . Pour , on calcule
La tangente est donc l’image par la rotation de centre et d’angle de la droite - si la fonction n’est pas périodique, il y a lieu d’étudier l’existence de limites de en , si la limite est nulle on s’approche en spiralant de l’origine, si elle est finie, il y a un cercle asymptote, si elle est infinie une spirale.
- comme ,
la vitesse (si le temps est ) est donnée par
où est une base orthonormée directe.
Donc si ou , le point est régulier et l’angle de la tangente avec vérifie (si et , la tangente est portée par ). Si , la tangente est portée par .2 - On ne peut avoir un point singulier que pour . On ne fait
pas leur étude comme en paramétriques, en effet la tangente est
toujours portée par , si change de signe la courbe
a la même allure que pour un point régulier, si ne change pas
de signe on a un rebroussement de première espèce (puisqu’on
traverse la tangente lorsque augmente)
Exemple :
- Convexité : pour avoir un point d’inflexion, il faut que
On peut le montrer de différentes manières :
- En calculant le déterminant de {vitesse,accélération } par rapport à dans le repère , on a
- En calculant la dérivée de l’angle fait avec l’axe qui vaut
- avec Xcas en se ramenant en paramétriques
où on a noté l’angle pour pouvoir dériver avec
'
et et les deux coordonnées.
- de même on calcule la courbure définie en section 4.2
3.7 Coniques
Les coniques sont des courbes implicites dont l’équation
cartésienne est du second degré
Exemples:
On va voir qu’elles sont de trois types : ellipses, hyperbole, parabole3 et on va les paramétriser, à partir de leur équation cartésienne ou à partir de leurs éléments géométriques (le calcul des éléments géométrique à partir de l’équation cartésienne fait intervenir l’étude des formes quadratiques, il ne sera pas abordé dans ce cours). Les coniques sont des courbes importantes en géométrie, ce qui a un intérêt en optique (parabole), mais aussi en cinématique (première loi de Kepler : l’orbite décrite par une planète est une ellipse dont le Soleil occupe un foyer).
3.7.1 Ellipse
Exemple : ouvrir un niveau de géométrie 2d dans Xcas,
choisir le mode ellipse cliquer 2 points (ce sont les foyers)
puis un 3ème point (point de l’ellipse), passer en mode
pointeur et faire bouger l’un des points, observer la forme
de l’ellipse qui en résulte. Ou dans une ligne de commande
normale taper la commande ellipse()
avec
en arguments les 2 points foyers et un point de l’ellipse
ou l’équation cartésienne de l’ellipse, par exemple
ellipse(-1,1,3+i)
trace l’ellipse de foyers
et passant par le point .
On note la distance entre les deux foyers, qui doit être plus petite que pour que l’ellipse soit non vide. L’excentricité de l’ellipse est définie par . Si , on obtient un cercle de centre et de rayon . Si , on va voir qu’il s’agit d’un cercle contracté selon l’axe perpendiculaire à dans un rapport de . On va également calculer l’équation en coordonnées polaires de (c’est sous cette forme que l’on montre que la Terre décrit une ellipse dont le Soleil occupe un foyer).
Soit le milieu de et , on se place dans le repère orthonormé dont le premier axe contient et donc les coordonnées de sont et celles de sont . Soit un point de l’ellipse, on a d’une part : et d’autre part : donc : en additionnant avec et en appliquant , on en déduit : En prenant le carré, on a : d’où : finalement : qui est bien la contraction selon de rapport du cercle de centre et de rayon (appelé grand cercle de l’ellipse).
En coordonnées paramétriques, on peut utiliser le paramétrage suivant :
En coordonnées polaires, on note la distance de à , et l’angle entre l’axe et . L’abscisse de est donc : que l’on combine avec (1) pour obtenir : donc :
Remarques :
- La première loi de Képler dit que l’orbite d’une planète autour du Soleil est une ellipse dont le Soleil occupe un des foyers. La troisième loi de Képler donne la relation suivante entre le demi-grand axe , la période de révolution et le produit de la masse du Soleil4 par la constante de gravitation :
- Attention, . Et dans le cas de l’orbite
de la Terre autour du Soleil, aucun de ces deux paramétrages n’est
le temps . Le paramétrage par le temps se déduit de la loi
des aires
Il nécessite de résoudre
une équation , cf. l’équation du temps dans le cours :
www-fourier.ujf-grenoble.fr/~parisse/climat/orbite.html
L’excentricité de la Terre est faible, en ce moment, elle est responsable d’une petite différence de durée des saisons, elle varie très lentement au cours des siècles sous l’action des autres planètes du système solaire, cette lente variation est une des raisons des glaciations du quartenaire. - On peut aussi définir géométriquement l’ellipse par un foyer et une directrice , c’est l’ensemble des tels que . Dans le repère d’origine , a pour équation (dans le repère d’origine , c’est ). En effet .
Exemple : faites varier la valeur de l’excentricité ci-dessous, que
voit-on pour E=0.0, E un peu inférieur à 1 (par exemple 0.8) et un peu supérieur
à 1 (par exemple 1.3)
3.7.2 Parabole
Si est un point et une droite ne passant pas par , la parabole de foyer et directrice est l’ensemble des points équidistants de et . En choisissant un repère tel que la droite ait pour équation et en prenant , appartient à la parabole si donc en passant au carré : La parabole est donc (ici) un graphe de fonction, donc un cas particulier de courbe paramétrique. On peut trouver son équation en polaire, en prenant comme origine donc l’équation de la droite devient et on a : cf. l’exercice sur les coniques données par foyer et directrice, qui traite aussi le cas des hyperboles. On peut aussi faire à titre d’exercice l’étude de la courbe en polaire : lorsque et .
Un intérêt majeur de la parabole en optique est que les rayons incidents perpendiculaires à la directrice se réfléchissent en passant par le foyer (on peut même montrer que cela caractérise une parabole). Illustration-démonstration avec Xcas dans un niveau de géométrie taper les commandes
P:=plotfunc(x^2/2+1/2,x=-5..5); supposons(a=[-1.4,-5,5,0.1]); D:=line(x=a,color=red); M:=inter_unique(P,D); T:=tangent(P,M); R:=symetrie(T,D,color=red); trace(R);
puis faire varier en cliquant sur les flèches. Pour tester en
ligne, commencez par initialiser la trace en exécutant
puis faites varier en cliquant sur le bouton +
ou -
:
inter_unique(R,line(x=0))
elle est indépendante de et est le foyer. On peut montrer qu’une courbe ayant cette propriété est une parabole.
3.7.3 Hyperbole
Une hyperbole de foyers et est définie comme l’ensemble des points tels que : où est une constante telle que , avec une excentricité .
En physique, les hyperboles interviennent dans les trajectoires non périodiques en mécanique céleste, mais aussi comme courbes de déphasage constant entre deux sources situées aux deux foyers (les figures d’interférence font apparaitre des hyperboles).
On peut faire un calcul analogue à celui de l’ellipse, on en déduit que l’équation cartésienne de l’hyperbole dans le repère centré au milieu des foyers, d’axe l’axe des foyers est donc : On peut paramétrer les deux branches de l’hyperbole par et en polaires
Exercice : faire l’étude de la courbe paramétrée et montrer que l’hyperbole admet deux asymptotes d’équation .
3.7.4 Paramétrisation rationnelle
Si on connait un point d’une conique, on peut effectuer
un changement d’origine en ce point, l’équation cartésienne
devient
On suppose que 5.
On cherche alors l’intersection de la conique avec la droite
(de pente ), on va voir que la droite coupe en général la
conique en deux points, l’origine et un autre point dont on
calcule les coordonnées en fonction de 6.
Graphiquement, par exemple
Si on cherche les points où le dénominateur en s’annule, on doit calculer (pour et en supposant que la fraction est irréductible7) le discriminant8 de l’équation du second degré Il y a trois cas possibles:
- si , il n’y a pas de racine, le paramétrage est défini pour tout et les limites en de sont nulles (car puisque ), la conique est bornée, c’est une ellipse.
- si , il y a une racine double, qui engendre une étude de branche infinie en , on obtient une parabole (deux branches selon que tend vers par la droite ou la gauche). Il n’y a pas d’asymptote, on a bien qui tend vers , mais ne converge pas (le numérateur a une racine simple qui ne compense pas la racine double au dénominateur)
- si , il y a deux racines distinctes , donc deux valeurs de où il faut faire une étude de branche infinie, on a alors une hyperbole avec 4 branches infinies et deux asymptotes parallèles à (en effet le rapport tend bien vers et a une limite car la racine au dénominateur de et est simple donc il y a simplification avec le numérateur)
Exercice : paramétrer et faire l’étude des coniques :
Remarque :
on a vu que les ellipses, paraboles, hyperboles admettent une
équation réduite du second degré. On en déduit facilement que
leur équation dans un repère quelconque est toujours
du second degré. Réciproquement, pour une équation cartésienne
on a calculé une paramétrisation rationnelle, mais pas
démontré que c’était forcément une conique. Pour faire
cela, l’outil adapté est l’étude des formes quadratiques. On peut
toutefois le faire à la main en dimension 2, en faisant une rotation
pour annuler le coefficient de . Par exemple
on voit que l’angle de la rotation à effectuer vérifie
- 1
- on notant et les deux fonctions pour pouvoir utiliser et dans droite
- 2
- Si , cela se lit sur l’expression de la vitesse qui est non nulle, mais c’est encore vrai si et non identiquement nul, pour le voir, on observe que a pour direction qui tend vers lorsque tend vers 0.
- 3
- En toute rigueur il faut ajouter deux autres cas ; l’ensemble vide et les paires éventuellement confondues de droites
- 4
- en négligeant la masse de la planète devant celle du Soleil
- 5
- Si , le polynôme est homogène et se factorise, on obtient l’origine ou la réunion de deux droites
- 6
- Cette méthode fonctionne pour les coniques, mais ne fonctionne malheureusement pas pour n’importe quelle équation cartésienne
- 7
- sinon, on aura deux droites parce que le polynôme se factorise en produit de deux facteurs de degré 1 dont
- 8
- On peut aussi voir ce discriminant comme le déterminant de la matrice de la forme quadratique associée
Chapitre 4 Propriétés métriques des courbes.
4.1 Longueur d’arc
La longueur d’un morceau de courbe régulier parcouru pendant un petit intervalle de temps est égal au premier ordre à la longueur du segment tangent parcouru, ou encore au produit de la norme de la vitesse instantanée par On remarque que cette quantité est invariante par changement de paramétrage, si alors On en déduit
En coordonnées polaires :
Remarque : il est très rare que l’on puisse effectuer
le calcul explicite d’une primitive de ,
il faut alors se contenter d’une valeur approchée de l’intégrale
lorsque et ont des valeurs numériques, calculée
par des méthodes numériques qui généralisent la méthode
des rectangles (cf. le cours de mat406).
Ce calcul se fait avec Xcas (ou une calculatrice formelle) en donnant
une valeur approchée à l’une des bornes.
Il y a quelques exceptions par exemple la longueur d’un arc de
parabole se calcule avec une formule explicite (essayez la
commande int(sqrt(1+4t^2),t,a,b)
ou
La cycloïde1
admet aussi une formule simple pour sa longueur
Par contre, la
longueur d’un arc d’ellipse ne se calcule pas avec les fonctions
usuelles (pour pouvoir le faire, il faut introduire des fonctions spéciales
adaptées, appelées intégrales elliptiques) :
4.2 Courbure, repère de Frenet, accélération normale et tangentielle.
Si on choisit , la longueur d’arc, comme nouveau paramètre de temps, la longueur parcourue est égale au temps, donc la vitesse instantanée par rapport à est de norme 1. On peut aussi le voir en notant : où est la norme de la vitesse avec comme paramètre, donc est bien égal à 1.
Calculons maintenant l’accélération avec ce nouveau paramètre . Comme la vitesse est de norme constante égale à 1, donc de carré 1, en dérivant par rapport à , on vérifie que l’accélération est perpendiculaire à la vitesse pour ce paramétrage par la longueur d’arc . L’accélération par rapport à est donc portée par la normale à la trajectoire, et sa mesure algébrique est appelé courbure (signée), notée , la valeur absolue de l’inverse de est appelé le rayon de courbure (la direction de l’accélération pointe vers le centre de courbure). Si on se déplace sur un cercle de centre et de rayon à vitesse 1, alors , la vitesse est donnée par donc de norme 1, et l’accélération par , sa norme vaut et sa direction pointe vers le centre du cercle. Donc la courbe est, à l’ordre 2 au point considéré, identique à un cercle de rayon .
Revenons au paramètrage initial . Dérivons par rapport à la vitesse , on obtient : L’accélération se décompose donc en deux parties
- le premier terme colinéaire au vecteur tangent est l’accélération tangentielle, de norme ,
- le second terme perpendiculaire au vecteur tangent est l’accélération normale, dont la norme est , où est le rayon de courbure
Autre formule de calcul du rayon de courbure : l’accélération normale vaut donc
Exemple : calcul du cercle osculateur en un point d’une parabole .
C:=cercle_osculateur(G,M)
. Pour visualiser que les cercles
sont inclus les uns dans les autres, initialiser
Remarques :
- La courbure est aussi la dérivée par rapport à l’abscisse curviligne de l’angle fait par la tangente avec une direction fixe, par exemple l’axe . En effet , dont la dérivée est le produit de par le vecteur normal .
- Génériquement, une courbe reste du même coté de sa tangente (car le terme suivant dans le développement est d’ordre 2, de signe constant en 0), les exceptions sont les points d’inflexion. Par contre, génériquement une courbe traverse son cercle osculateur (en y rentrant ou en en sortant), car le terme suivant dans le développement de la différence entre les points des deux courbes est d’ordre 3 et change donc de signe en 0. Les exceptions (tangeance courbe-cercle osculateur d’ordre 3 au lieu de 2) sont appelés sommets d’une courbe, par exemple le sommet d’une parabole.
- on peut calculer les coordonnées du centre du cercle osculateur de manière algébrique (i.e. sans introduire de racines carrées) à partir des coordonnées paramétriques de et de ses dérivées
- la courbe décrite par les lorsque parcourt
la courbe étudiée est appelée
développée
de la courbe .
La vitesse de vaut
on en déduit que la tangente à la développée en a
pour direction la normale (si admet un
point critique, par exemple en un sommet de la courbe,
la développée admet génériquement
un point de rebroussement2).
L’enveloppe3
des normales à une courbe est donc sa développée.
Exemple : développée de l’ellipse
Initialiser la trace avec
puis faire varier :Not evaledOn observe 4 sommets pour l’ellipse, situés sur les grands et petits axes, et donc 4 points de rebroussements pour la développée. Ouvrir un niveau de géométrie 2d dans Xcas, taper une commande par ligne
G:=plotparam([2*cos(t),sin(t)],t=0..2*pi);
M:=element(G);
T:=tangent(M);
N:=perpendiculaire(M,T);
trace(N)
passer en mode pointeur (menu mode du niveau de géométrie) et faire bouger le point le long d’un quart de l’ellipse, ceci trace un faisceau de normales à l’ellipse, dont on voit apparaitre l’enveloppe (limite entre la région couverte et non couverte par des points du faisceau de normales), cette enveloppe est la développée de l’ellipse (vous pouvez utiliser le menu M à droite du dessin pour effacer les traces). Avec Xcas version 1.1.1-18 ou ultérieure on peut tracer la développée avec la commandedeveloppee(G)
.
Cf. aussi les animations de l’article Développée de wikipedia.
Exercice : calculer le repère de Frenet pour une ellipse , par exemple puis le rayon de courbure, puis la développée (on obtient une courbe image par affinité d’une astroïde). Donner une équation paramétrique simple de .
Vérification avec Xcas (version à jour) :
- De plus, comme est normé,
la longueur d’arc de courbe de la développée est donnée par :
Conséquence : si on enroule un fil sur la développée , que ce fil est tendu et que son extrémité coïncide, avant de commencer à le dérouler, avec un point de la courbe alors dans la suite du déroulé, l’extrémité parcoura la courbe (on dit que est une développante de ). - si la courbure est de signe constant, les cercles osculateurs sont inclus les uns dans les autres, c’est une conséquence de la remarque précédente, de l’inégalité triangulaire et du fait que la distance entre deux centres de cercles osculateurs est plus petite que la longueur d’arc sur la développée. Comme les points de la courbe sont sur des cercles osculateurs, il en résulte que la courbe entre ou sort définitivement du cercle osculateur au point de contact.
- La développée peut servir à calculer une caustique en
optique. On envoie des rayons lumineux parallèle à une
direction fixée vers un miroir ayant la forme de la courbe
la caustique
est l’enveloppe des rayons lumineux réfléchis
(que l’on observe par une plus grande intensité lumineuse).
On peut montrer que la caustique est la développée
de l’anticaustique de par rapport à une droite perpendiculaire
aux rayons lumineux (pour déterminer
l’anticaustique d’une courbe par rapport
à une droite, on prend un point de la courbe, on le projète sur
la droite puis on prend le symétrique du projeté par rapport à la
tangente à la courbe au point choisi,
l’anticaustique est le lieu de ces symétriques).
Cf. dans Xcas la session exemple du
menu
Exemple, geometrie, caustique
. Les développées peuvent aussi servir dans le calcul de caustiques par réfraction :http://www.mathcurve.com/courbes2d/caustic/caustic.htm
- On peut faire une étude analogue pour une courbe dans l’espace, dans ce cas la dérivée de par rapport à l’abscisse curviligne fait intervenir une composante sur le troisième vecteur du repère direct ( comme binormal), appelé torsion.
- Équation intrinsèque d’une courbe :
Il s’agit de trouver une courbe vérifiant une relation entre la courbure (ou rayon de courbure) et l’abscisse curviligne, par exemple la relation avec fixé.4 Pour trouver une telle courbe, on la paramètre par l’abscisse curviligne , donc la vitesse est de norme 1 et caractérisée par l’angle fait avec une direction fixe, on a alors et on en tire puis . Dans l’exemple, on a , donc (en choisissant la direction fixe pour annuler la constante d’intégration), puis : puis en choisissant l’origine du repère : en posant on a aussi Pour , le tracé est obtenu par la commande
- 1
- qui tire son nom de la trajectoire d’un point fixé à un cercle roulant sans glisser sur une droite, par exemple l’extrémité d’un rayon sur une roue de vélo.
- 2
- On peut montrer qu’une courbe convexe admet au moins 4 sommets (théorème des quatre sommets), sa développée admet donc au moins 4 points de rebroussements.
- 3
- L’enveloppe d’une famille de droites est une courbe dont l’ensemble des tangentes est la famille de droite
- 4
- Ce type de courbe, appelé spirale d’Euler ou de Fresnel ou clothoïde, est utilisée pour faire des raccordements de chemin de fer (ou de route) entre une portion rectiligne, où l’accélération normale est nulle, et un arc de cercle, où l’accélération normale est constante, en effet si est constant alors l’accélération normale varie linéairement en fonction de l’abscisse curviligne donc du temps à vitesse constante. C’est plus agréable pour les passagers qui passent d’une accélération nulle à une accélération constante progressivement, mais aussi pour créer une pente progressive latérale sur les rails pour compenser la force centrifuge par la gravité et éviter une usure prématurée du rail.
Chapitre 5 Formes différentielles et intégrales curvilignes
Il s’agit ici de calculer des intégrales le long d’un arc de courbe. Cela intervient directement par exemple pour calculer le travail d’une force au cours d’un déplacement le long d’une courbe ou la quantité de chaleur/travail pendant un cycle en thermodynamique (le long d’une courbe dans le plan défini par deux coordonnées indépendantes comme par exemple pression-température ou pression-volume) ou indirectement en transformant un calcul d’intégrale double à l’intérieur d’un domaine en intégrale curviligne sur le bord du domaine (calcul d’aire, de centre d’inertie, de moment d’inertie...). Dans les cas favorables, on a un analogue des primitives, on peut calculer un potentiel et faire la différence de potentiel entre les deux extrémités du chemin pour calculer l’intégrale curviligne.
On va d’abord définir ce qu’on peut intégrer le long d’une courbe, à savoir une forme différentielle (aussi appelée 1-forme), puis on donnera quelques résultats sur les formes fermées et exactes (c’est le cas favorable, il correspond aux forces conservatives en mécanique ou aux différentielles totales de fonctions d’état en thermodynamique).
5.1 Forme différentielle
Soit une fonction de deux variables continument dérivable. On s’intéresse aux variations de lorsqu’on se déplace dans le plan depuis le point dans une direction donnée à la vitesse . On a alors une formule équivalente à celle de la dérivée d’une fonction d’une variable :
En effet : donc
Exemples :
- la dérivée de selon la direction (axe des ) est et selon la direction (axe des ) est .
- Soit et la distance de au point . Alors existe en tout point différent de et est le vecteur unitaire porté par .
Remarque : Différentielle et gradient, forme différentielle et champ de vecteurs
La différentielle a les mêmes composantes que
le gradient de (noté , gradient(V,[x,y])
avec Xcas),
mais ce ne sont pas les mêmes objets :
en un point donné
est une application linéaire (qui a un sens indépendamment
de la définition d’un produit scalaire)
alors que est un vecteur (dont la relation
avec la dérivée directionnelle dépend du produit scalaire),
on a pour tout vecteur la relation
On a la même relation entre le travail d’une force (que l’on peut voir comme une forme linéaire qui s’applique sur les vecteurs déplacement) et la force correspondante qui est un champ de vecteurs . Le champ de vecteurs correspond à la généralisation d’une différentielle, appelée forme différentielle en utilisant le produit scalaire : pour tout vecteur , on a Dans le cas d’une force conservative, est l’opposé de la différentielle du potentiel de la force (le signe moins vient de la convention utilisée en physique ).
On parle parfois de vecteur covariant pour la différentielle (et vecteur contravariant pour le gradient).
Applications :
-
Tangente à une courbe de niveau :
le vecteur tangent en un point d’une courbe de niveau de est dans le noyau de l’application linéaire en ce point (puisque est constant le long de le courbe), ou de manière équivalente est orthogonal à la courbe de niveau. - Calcul du gradient en coordonnées polaires.
le repère est orthonormé, pour connaitre les coordonnées de dans ce repère il suffit de calculer la dérivée directionnelle de dans les directions et car : Or la dérivée directionnelle selon est la dérivée partielle de par rapport à , et la dérivée directionnelle selon est la dérivée partielle de par rapport à divisée par (il faut diviser par parce qu’on se déplace tangentiellement au cercle de rayon ), donc - Tangente et la normale en un point d’une ellipse :
Ce sont la bissectrice extérieure et intérieure issues de du triangle défini par et les foyers. Cela vient du fait que le gradient de la distance est le vecteur unitaire porté par et que l’ellipse est courbe de niveau de la somme des distances aux foyers. On retrouve ainsi qu’un rayon lumineux issu d’un foyer se réfléchit sur l’ellipse en passant par l’autre foyer, et la longueur du chemin parcouru est constante (donc des rayons lumineux issus d’un foyer arrivent en phase à l’autre foyer quel que soit la direction initiale) - De même on peut déterminer géométriquement la tangente et la normale à une parabole, ce sont les bissectrices issues de de la droite (où est le foyer) et de la perpendiculaire en à la directrice de la parabole, on retrouve ainsi que les rayons lumineux perpendiculaires à la directrice se réflechissent sur la parabole en passant par le foyer) et à une hyperbole (comme pour une ellipse).
On note donc [resp. ] la différentielle de [resp. ]1 on a :
Une forme différentielle est la généralisation de la différentielle d’une fonction, elle s’écrit sous la forme où et sont des fonctions des deux variables , mais pas forcément les dérivées partielles d’une fonction .
La définition géométrique d’une forme différentielle est la donnée en tout point du plan (ou d’un domaine ouvert du plan) d’une application linéaire de à valeur dans 2 (ou en tout point de l’espace d’une application linéraire de à valeurs dans pour une courbe de ). Si on prend la base canonique de , une application linéaire de dans est caractérisée par sa matrice qui possède une ligne et deux colonnes et a donc deux coefficients et , une forme différentielle équivaut donc bien à la donnée d’un couple de fonction .
5.2 Intégrale curviligne
Exemple: on prend et on calcule l’intégrale curviligne le long de l’arc de parabole pour , on obtient En paramétrant par avec on retrouve le même résultat.
La valeur de l’intégrale est bien définie indépendamment du paramétrage, en effet si on change de paramétrage avec une bijection envoyant sur , on a (en utilisant la linéarité de à la deuxième ligne) :
Attention à l’orientation, si on change d’orientation, on change le signe, par exemple si on parcourt l’arc de parabole de vers , en utilisant le paramétrage , on obtient l’opposé :
Remarque : le travail d’une force le long d’un arc de courbe est donné par l’intégrale curviligne de la forme différentielle .
L’intégrale curviligne d’une forme différentielle reliant deux points et d’un arc de courbe se calcule en choisissant un paramétrage de , si est une courbe paramétriques, on prendra en général le paramétrage définissant , si est une courbe on peut prendre , si est une courbe en polaires , on peut prendre ().
Pour certaines formes différentielles, on peut faire comme en dimension 1, trouver une primitive, voir la section ci-dessous.
Calcul avec Xcas : sur les versions les plus récentes de Xcas (1.9.0-996),
la commande integrate permet de calculer une intégrale curviligne.
On définit la courbe et son paramétrage, par exemple
puis on passe à integrate
le champ de vecteur associé à la forme différentielle,
exprimé en fonction des variables ,
la liste des variables , la courbe , et les valeurs de début et fin du paramètre
Remarque : l’intégrale curvligne se généralise en dimension supérieure à 2.
5.3 Forme différentielle exacte
Voyons maintenant à quelle condition il existe un analogue du calcul avec une primitive. On a: En effet, on applique la définition (7) où et :
Pour une force qui dérive d’un potentiel, on a donc montré que le travail de la force se calcule en faisant la différence de potentiel entre les deux extrémités. Cette propriété, analogue au calcul d’intégrale classique en utilisant une primitive n’est pas automatique, car elle implique que l’intégrale curviligne ne dépend pas du chemin choisi pour relier les deux points. Or en thermodynamique, la chaleur est modélisée par une forme différentielle, mais la chaleur échangée dépend du chemin suivi (c’est vrai aussi en mécanique pour le travail de forces non conservatives comme les forces de frottement). En mathématiques, on parle de forme différentielle exacte ou non exacte.
Si on choisit comme chemin un segment entre deux points et d’ordonnées identiques et d’abscisses et (i.e. ), alors en divisant par et en faisant tendre vers 0, on a De même . Réciproquement, si et alors donc
Si est deux fois continument différentiable alors . D’où une condition nécessaire pour que soit exacte :
Une forme exacte est toujours fermée, mais la réciproque n’est pas
toujours vraie, une forme fermée n’est pas forcément exacte, cela
dépend où elle est définie. Si elle est définie dans un
domaine ouvert de sans trou ( tout entier,
un rectangle, un disque, etc.),
on peut montrer qu’une forme
fermée est une forme exacte, on se fixe un point et
on définit comme pour
un chemin quelconque reliant à , on montre
que le résultat ne dépend
pas du choix du chemin en appliquant le théorème
de Stokes (voir section suivante). Sinon, il existe des
contre-exemples, comme sur le cercle unité
La forme est fermée :
mais elle n’est pas exacte :
Pour trouver le potentiel
dont une forme différentielle fermée
est la différentielle, on résoud d’abord par exemple
en intégrant par rapport à , étant considéré
comme un paramètre, on obtient à une constante d’intégration
près, cette constante d’intégration en peut dépendre
de , c’est donc une fonction ,
on remplace dans et on intègre
en pour trouver la valeur de (à une constante près).
Cette opération est executée par la commande potential()
de Xcas.
Exemple : Cette forme est bien fermée : La forme est définie dans tout entier, donc est exacte, on intègre par rapport à , on trouve une constante d’intégration, qui est donc constante par rapport à donc est une fonction dépendant de . On détermine ensuite en dérivant par rapport à qui se simplifie en donc une constante, que l’on peut prendre nulle :
Si une forme n’est pas fermée, elle n’est pas exacte, et on ne peut pas calculer une intégrale curviligne par différence de potentiel, il faut utiliser la définition et un paramétrage de , ce qui n’est pas forcément plus couteux en calcul et a des applications au calcul d’intégrale double si est un chemin fermé délimitant le domaine d’intégration.
Les formes exactes ont une autre application (anticipant sur le chapitre suivant) : la recherche d’intégrales premières d’équations différentielles. En effet si le long d’un arc paramétrable par alors et est le graphe d’une solution de cette équation différentielle.
Si n’est pas exacte, il n’y a pas de potentiel mais il peut arriver qu’en multipliant la forme par une fonction, on trouve une nouvelle forme qui elle est fermée, on parle alors de facteur intégrant. (On limite en général la recherche à des fonctions ne dépendant que de ou de ).
La notion de facteur intégrant ne se limite pas à la résolution d’équations différentielles. Par exemple en thermodynamique, la forme chaleur n’est pas fermée, mais en divisant par la température on obtient une forme fermée dont le potentiel est l’entropie.
5.4 Intégrale curviligne et intégrales doubles.
Terminons ce chapitre par le lien entre intégrale curviligne sur un lacet (chemin fermé) et intégrale double à l’intérieur du lacet. C’est évidemment surtout intéressant pour les formes non exactes, car si est un lacet et une forme exacte, alors . On a le théorème de Stokes, aussi appelé en dimension 2 formule de Green-Riemann :
Idée de la preuve : on commence par le cas où est un rectangle , on peut alors calculer on compare avec les intégrales curvilignes sur les segments verticaux orientés et . De même pour et les segments horizontaux.
Pour le cas d’un domaine d’intégration plus général, on approche par une réunion disjointe de petits rectangles.
Application : pour calculer l’aire d’un domaine de frontière , il suffit de calculer l’une des intégrales curvilignes : Par exemple, l’aire à l’intérieur de l’ellipse vaut
On peut aussi calculer des moments d’inertie
ou la position d’un centre de gravité
en se ramenant à une intégrale curviligne.
Exemple : Calculer la position du centre d’inertie d’un quart
de cercle .
On a donc délimité par , réunion de , et .
Pour trouver la position du centre
d’inertie en (en c’est identique), on doit calculer
et on divise par l’aire du quart de cercle, on trouve
donc , on peut visualiser avec
la commande
Remarque
En dimension 2,
si on calcule l’intégrale curviligne du produit scalaire
entre un champ
et le vecteur normal sortant
au contour d’intégration,
alors où
désigne un paramétrage de la courbe par la longueur d’arc
dans le sens trigonométrique,
donc on intégre par rapport à la fonction ,
la forme différentielle est ,
le théorème de Stokes nous donne
l’égalité de l’intégrale de contour avec
l’intégrale double à l’intérieur
La formule analogue en dimension 3 et la relation entre la divergence
du champ électrique et charge permet de montrer que le flux du champ
à travers une surface est proportionnel à la charge intérieure.
5.5 Formulaire courbe et intégrale curviligne
5.5.1 Courbes paramétriques
- Recherche d’asymptote : si ou tend vers l’infini, si un seul tend vers l’infini, asymptote horizontale ou verticale, sinon , si et asymptote oblique
- Point singulier : si . Si l’accération est non nulle, elle dirige la tangente, si de plus la dérivée troisième est non colinéaire on a un rebroussement de première espèce. En général, discuter sur la première dérivée ième non nulle et la première dérivée ième non colin’éaire.
- Convexité : pente de la tangente , chercher le signe de ou de
- longueur , la vitesse,
- courbure
- repère de Frenet formé de , et tel que soit orthonormé direct,
- rayon de courbure
- accélération tangentielle et normale avec et
- cercle osculateur de centre , parcout la développée de la courbe.
5.5.2 Coniques
-
Définition par directrice , foyer et excentricité ,
,
Équation polaire (origine en , d’équation ) - Définition par équation du second degré , le signe de discrimine entre ellipse (-), parabole (nul) ou hyperbole (+).
- Ellipse: , , , .
Équation cartésienne réduite .
Paramétrique réduite
Équation polaire (origine en )
Tangente et normale: bissectrices extérieure et intérieure au point des directions issues des foyers. - Hyperbole , , ,
Équation cartésienne réduite .
Paramétrique réduite
Équation polaire (origine en )
Tangente et normale: bissectrices intérieure et extérieure au point des directions issues des foyers. - Tangente et normale à la parabole : bissectrices de la direction issue du foyer et de la perpendiculaire à la directrice.
5.5.3 Courbes en polaires
- vitesse , d’affixe dirige la tangente, sauf si (point singulier), dans ce cas la tangente fait un angle de avec l’axe des et on a un rebroussement de première espèce si ne change pas de signe.
- recherche d’asymptote si pour , si existe, asymptote dans le repère tourné de
- convexité : étude du signe de ou de
- courbure
5.5.4 Intégrales curvilignes
- forme différentielle
- est exacte si , , et le gradient de .
- est fermée si .
- exact entraine fermé, la réciproque est vraie si le domaine de définition n’a pas de trous.
- calcul d’une intégrale curviligne :
- si est exacte, . Les lignes de niveau de ont leur tangente orthogonale au gradient de et sont solutions de .
- Si est un lacet orienté dans le sens
trigonométrique
- 1
- Géométriquement, [resp. ] est la forme linéaire constante (i.e. indépendante du point du plan choisi) qui a tout vecteur de associe sa première [resp. deuxième] coordonnée :
- 2
- Pour être complet, on suppose de plus que cette application linéaire qui dépend du point du plan en dépend de manière au moins continue et presque toujours de manière continument différentiable
Chapitre 6 Équations et systèmes différentiels.
6.1 Introduction et représentation graphique.
On s’intéresse à l’équation différentielle où et . Si , c’est une équation différentielle, si c’est un système différentiel.
Exemple : en dimension , . On sait résoudre cette équation, les solutions sont de la forme . Si on trace la courbe représentative de ces solutions (appelée courbe intégrale), on observe que par tout point du plan, il passe une solution unique. La tangente à une courbe intégrale a pour pente donc pour vecteur directeur le vecteur de composantes .
C’est vrai de manière plus générale, le vecteur directeur
de la tangente à une courbe intégrale est . Si
on représente dans le plan selon un quadrillage régulier
les vecteurs , une courbe intégrale doit être
tangente à ces vecteurs chaque fois qu’elle passe en un point
du quadrillage, (et à peu près tangente si elle passe à
proximité). Un tel quadrillage est appelé champ des tangentes
(commande plotfield
en Xcas, mode également disponible
sur certaines calculatrices).
Exercice : tracer le champ des tangentes et quelques solutions
pour quelques exemples de fonction , avec Xcas créer
une figure 2d, puis choisir le mode Champ des tangentes
du menu Geo, Graphe, entrer la fonction, puis cliquer en quelques
points pour faire tracer la solution passant par ces points.
Par exemple pour
L’équation (3) est d’ordre 1, or certaines équations différentielles se présentent naturellement comme des équations d’ordre 2, par exemple l’équation fondementale de la dynamique (accélération=somme des forces divisée par la masse). Mais on peut facilement se ramener à un système différentiel d’ordre 1, en augmentant la dimension de . Par exemple, si on pose , où est la position et la vitesse, alors l’équation devient un système d’ordre 1 où est la force, qui dépend de la position (champ électrique, gravitation...) et éventuellement de la vitesse (force de frottement, champ magnétique...). On utilise aussi assez fréquemment où est la position, et la quantité de mouvement (qui dépend de la vitesse, linéairement en mécanique classique).
Représentation graphique : comme précédemment,
on peut se placer dans l’espace
des (si est en dimension 1), mais il est souvent
plus difficile d’observer des phénomènes sur un graphe
en 3-d que dans le plan, on préfère ne pas représenter explicitement le
temps , mais uniquement , on est donc naturellement
ramené à représenter une solution (une courbe intégrale)
par une courbe paramétrique en (ou en position
impulsion). On a encore la notion de champ des tangentes
si ne dépend pas explicitement du temps
(on dit que le système est autonome), dans ce
cas une courbe intégrale a pour tangente en
de direction portée par le vecteur .
Exemple : . La commande
permet d’en représenter le champ des tangentes
et d’avoir une idée approximative de l’allure des solutions
(normalize
affiche des vecteur tangents de norme 1, si
on n’utilise pas cette option, la taille des vecteurs tangents donne
la “vitesse” de déplacement).
On sait résoudre ce système différentiel, soit en appliquant
une technique matricielle présentée ci-dessous,
soit en se ramenant à une équation
linéaire d’ordre 2 à coefficients constants:
donc , étant déterminés
par les conditions initiales sur .
Une équation donnée sous la forme (3) est appelée une équation résolue en , car on a exprimé la dérivée en fonction de et de . Il existe (plus fréquemment en mathématiques) d’autres formes d’équations différentielles (non résolues) où le premier travail de résolution peut consister à exprimer en fonction de et (ce qui n’est pas toujours possible explicitement).
Exemple : en dimension 1, , on sait résoudre exactement
cette équation à variables séparables, les solutions sont de la
forme .
On observe que contrairement à
où passe une solution et une seule par chaque point du plan,
ici toutes les solutions valent 0 en : il passe une
infinité de solutions par le point et il n’en passe aucune
par . Ce phénomène de non unicité/non
existence vient de la mise sous forme résolue qui
fait apparaitre une singularité de en .
On présente dans la suite de cette section des résultats qualitatifs sur les équations sous forme résolue lorsqu’on ne sait pas les résoudre, ainsi que quelques méthodes explicites pour certaines équations différentielles que l’on sait résoudre.
6.2 Existence et unicité
Il s’agit ici de préciser dans quelles conditions le résultat intuitif suivant est vrai : étant donné une condition initiale , il y a une et une seule évolution possible, donc une solution unique de l’équation ou du système (3).
On a le :
Remarques
- Attention, l’existence d’une solution ne signifie absolument pas que l’on sait calculer explicitement .
- L’existence et l’unicité d’une solution permet d’affirmer le caractère déterministe de l’équation. Mais la méconnaissance précise de la condition initiale peut au cours du temps provoquer une erreur tellement grande sur que celle-ci devient imprédictible.
- Le théorème ne dit rien sur la taille de l’intervalle
d’existence de la solution (en temps). Certaines solutions
peuvent exploser en temps fini, par exemple
desolve(y'=y^2 and y(0)=1)
Bien entendu, si on modélise une quantité physique par , dire que explose (en temps fini ou infini du reste) signifie que les approximations utilisées pour la modélisation ne sont plus valable bien avant.
On admettra ce théorème, voici quelques idées heuristiques de la preuve. L’équation peut se réécrire sous la forme intégrale équivalente Si est assez proche de , on peut approcher l’intégrale par C’est exactement ce qu’on fait en suivant le champ des tangentes pour approcher une courbe intégrale graphiquement, et si on discrétise le temps avec un pas petit, cette méthode d’approximation est appelée méthode d’Euler. On peut bien sur utiliser d’autres approximations (meilleures) de l’intégrale pour avoir une meilleure approximation de la solution, et les méthodes dites de Runge-Kutta utilisent cette idée. D’un point de vue théorique, la preuve repose plutôt sur ce qu’on appelle le théorème du point fixe, on met la valeur approchée de trouvée dans l’équation intégrale pour avoir une nouvelle valeur approchée de , on recommence, ainsi de suite, et on montre que le processus converge (il s’agit mathématiquement parlant d’une suite récurrente de fonctions, la preuve rigoureuse de la convergence nécessite des outils mathématiques de niveau L3-M1 de maths, c’est l’analogue des suites récurrentes de réels qui permettent de résoudre numériquement des équations comme abordées en mat406).
Conséquence du théorème 13 : deux courbes intégrales de la même équation différentielle ne peuvent se couper dans . Donc si on connait une courbe intégrale de et qu’on prend une condition initiale en-dehors de cette courbe, la courbe intégrale unique passant par cette condition initiale restera du même coté de . Si on connait deux courbes intégrales de , une courbe intégrale passant par une condition initiale entre les deux courbes restera entre les deux courbes.
Exemple : (équation logistique).
Exercice : toujours pour que se passe-t-il pour une condition initiale ?
6.3 Quelques méthodes de résolution explicite.
6.3.1 Équations à variables séparables
Si on peut factoriser en termes ne dépendant que de ou ne dépendant que de , on dit que l’équation est à variable séparable Cette équation admet des solutions constantes lorsque . Si , par le théorème de Cauchy-Lipschitz ne s’annule nulle part sur son domaine de définition. On peut donc diviser par et intégrer : On obtient une équation implicite de la forme où est une primitive de , de et une constante arbitraire. Dans les cas favorables, on peut exprimer en fonction de (par exemple si l’équation est linéaire sans second membre, on a donc est le log que l’on sait inverser). Dans les cas moins favorables, on peut exprimer et en fonction d’un paramètre : la courbe intégrale est une courbe paramétrée. Dans les cas défavorables, on reste sous forme implicite.
Exercice : résoudre explicitement l’équation
et retrouver les résultats qualitatifs de la section précédente.
6.3.2 Équations linéaires
On commence par résoudre l’équation sans second membre (aussi appelée homogène) sur un intervalle ouvert sur lequel . L’ensemble des solutions est un espace vectoriel (car l’équation est linéaire) et de dimension l’ordre de l’équation : pour le prouver on peut appliquer le théorème de Cauchy-Lipschitz au système d’ordre 1 équivalent, ce système est tel que est un vecteur de , on a ensuite un isomorphisme entre les solutions et la condition initiale.
Si l’ordre est 1, on a une équation à variables séparables et la solution est une exponentielle : Exemple : , on a
Si l’ordre est plus grand que 1, on n’a en général pas de solution explicitable avec les fonctions usuelles et des primitives2, pour certaines équations importantes en physique, des fonctions spéciales ont été créées pour exprimer les solutions, par exemple les fonctions de Bessel. Il existe quelques cas particuliers où le calcul explicite est possible, dont le cas où les coefficients sont constants (section suivante). Si on connait une solution d’une équation linéaire, alors en posant , la fonction vérifie une équation linéaire d’ordre un de moins, ainsi si on connait une solution d’une équation linéaire d’ordre 2, on peut la résoudre complètement.
Le calcul d’une solution particulière d’une équation linéaire avec second membre se fait en faisant varier les constantes d’intégration : on prend la forme générale de la solution de l’équation homogène, on remplace les constantes d’intégration par des fonctions inconnues, on remplace dans l’équation avec second membre et on résoud en les fonctions inconnues, la méthode détaillée dans le cas des coefficients constants s’applique à l’identique. La solution générale est la somme d’une solution particulière et de la solution générale de l’équation sans second membre.
Exemple : , solution générale de l’équation homogène , variation de la constante on remplace dans et on obtient , donc et , d’où la solution générale .
6.3.3 Équations linéaires à coefficients constants
On peut chercher des solutions de l’équation sans second membre sous la forme d’exponentielles , doit alors vérifier une équation polynomiale appelée équation caractéristique, de degré le degré de l’équation différentielle. Plus précisément, si on remplace dans alors
En effet, on a le bon nombre d’éléments (), il suffit donc de montrer qu’il s’agit d’une famille libre.
Pour cela, il suffit de faire tendre vers l’infini si toutes les parties réelles des sont distinctes. Si certaines sont égales, on peut faire tendre vers l’infini sur une demi-droite du plan complexe, on pose où est un petit angle choisi pour que les parties réelles de soient toutes distinctes (ceci revient à choisir tel que les projections des racines sur l’axe tourné de soient toutes distinctes, on brise ainsi une symétrie sur les racines de ).
On peut aussi faire une récurrence. Au rang c’est évident. Si et si vérifient : on factorise et on dérive, on a on est ramené à l’identité précédente au rang donc par récurrence, et si , puis avec la relation du départ.
Dans le cas général, on a :
En effet, on a : Supposons que est solution de l’équation, alors est encore solution si : et on reconnait une équation différentielle linéaire à coefficients constants dont l’équation caractéristique est , on a donc montré la :
Pour prouver le théorème, il faut encore se convaincre que la famille est une base, ce qui revient à prouver l’indépendance linéaire de ces fonctions. Cela se fait comme précédemment. Si toutes les parties réelles des racines sont distinctes, en faisant tendre vers l’infini on se ramène à un équivalent polynômial nul. Si certaines parties réelles sont égales, on peut à nouveau faire tendre vers l’infini dans le plan complexe en tournant d’un petit angle. On peut aussi adapter la récurrence, sur le nombre de racines. S’il y en a une, on a un polynôme. Sinon, on factorise , et on dérive la multiplicité de pour appliquer le résultat au cran , on a alors un système triangulaire sur le groupe d’inconnues de la même exponentielle. On peut aussi se ramener à des polynômes en posant comme ci-dessus et en faisant tendre vers l’infini.
Si est à coefficients réels et admet une racine non réelle alors est encore racine, on peut réécrire la base d’exponentielles complexes en une base de fonctions réelles en utilisant des fonctions trigonométriques. En effet les combinaisons linéaires réelles de et ont des coefficients conjugués :
Exemples :
- , équation caractéristique , deux racines distinctes , donc
- , équation caractéristique , a une racine double , donc
- , équation caractéristique , deux racines conjuguées , donc
On peut trouver une solution particulière de l’équation avec second membre par la méthode de variation des constantes, qui s’applique d’ailleurs également lorsque l’équation est linéaire à coefficients variables. Si la solution générale est engendrée par , on pose : On pose et ainsi de suite jusqu’à la dérivée d’ordre de , ces équations et l’équation différentielle donnent alors un système linéaire en les . Ce système a pour déterminant la matrice de -ième ligne la dérivée d’ordre de . Si on calcule la dérivée de , il faut faire porter la dérivée sur la dernière ligne pour avoir une contribution non nulle, on peut alors utiliser l’équation différentielle pour montrer que vérifie une équation linéaire d’ordre 1 donc il est soit identiquement nul soit jamais nul. On montre alors par récurrence que l’indépendance linéaire de entraine alors que le déterminant est non nul. En effet s’il est nul au rang , alors l’une des colonnes est combinaison linéaire des autres, par exemple la dernière, on a donc : en dérivant on en déduit que on est ramené à un système linéaire homogène en inconnues (les ) à qui on applique l’hypothèse de récurrence, on en déduit que les sont nuls donc les sont des constantes ce qui contredit l’indépendance linéaire des .
Pour des second membre combinaison linéaire de termes avec polynôme, il est plus simple de chercher directement une solution particulière combinaison linéaire de où est de même degré que si n’est pas racine de , ou de degré le degré de plus la multiplicité de comme racine de . On peut aussi utiliser la transformation de Laplace et son inverse.
6.3.4 Systèmes différentiels linéaires à coefficients constants d’ordre 1.
Il s’agit donc de systèmes de la forme où , est une matrice carrée de taille indépendante du temps, et .
On commence par résoudre l’équation homogène .
Si la matrice est diagonalisable, alors où
=diag est diagonale et inversible, le système devient :
donc en posant , on a (puisque est indépendant du temps) :
donc , puis la solution générale
Le calcul avec Xcas se fait en utilisant la commande desolve
,
par exemple
desolve(y'=[[1,2],[2,1]]*y)
ou avec conditions initiales
desolve(y'=[[1,2],[2,1]]*y and y(0)=[1,2])
On peut aussi utiliser la fonction
exp
avec comme argument (on généralise ainsi
la notation
de la dimension 1), multiplié par la condition initiale :exp([[1,2],[2,1]]*t)*[1,2]
Les calculs intermédiaires pour diagonaliser la matrice sont exécutés par les commandes
eigenvals
, eigenvects
,
jordan
.On peut ensuite calculer une solution particulière par la méthode
de variation des constantes, ou encore en résolvant
composante par composante (ou par transformation
de Laplace). Avec Xcas,
il suffit d’ajouter le second membre dans la commande desolve
desolve(y'=[[1,2],[2,1]]*y+[x,x+1])
Si la matrice n’est pas diagonalisable (ce qui entraine qu’elle a au moins une valeur propre de multiplicité plus grande que 1), on peut alors la trigonaliser, on se ramene à résoudre un système triangulaire, ce qui revient à résoudre pour chaque composante une équation différentielle linéaire d’ordre 1 avec un éventuel second membre.
6.3.5 Systèmes et équations
Il y a un lien entre systèmes différentiels linéaires et équations linéaires. En effet une équation d’ordre peut s’écrire comme un système différentiel d’ordre 1, on peut calculer le polynôme caractéristique de la matrice on retrouve alors l’équation caractéristique. Inversement, toute matrice admet un polynôme annulateur tel que 3, le polynôme caractéristique de est un polynôme annulateur (théorème de Cayley-Hamilton). Les composantes des solutions du système différentiel sont des solutions de l’équation différentielle dont l’équation caractéristique est . En effet :
Exemple en dimension 2. Soit Si alors on en déduit puis . Supposons donc , alors (on peut vérifier que ) donc si alors et s’en déduit avec (on peut du reste partir de cette relation pour établir l’équation d’ordre 2 vérifiée par ). On peut ainsi résoudre tous les systèmes de dimension 2, même si la matrice n’est pas diagonalisable.
Exercice : Résoudre de cette manière le système
desolve(y'=[[1,2],[2,1]]*y and y(0)=[1,2])
Autre exemple : système d’ordre 2 se ramenant à une équation d’ordre 2 à coefficients complexes. Les équations pour une particule chargée dans un champ magnétique constant porté par l’axe et un champ électrique constant perpendiculaire (donc dans le plan ), avec vitesse initiale nulle ou contenue dans le plan donnent une trajectoire plane Si on pose alors vérifie l’équation Le polynôme caractéristique de cette équation possède deux racines distinctes 0 et (mais pas le conjugué, l’équation n’est pas à coefficients réels!) donc la solution homogène est Le champ électrique joue ici le rôle de second membre, comme 0 est solution de l’équation caractéristique, la forme de la solution particulière est , en remplaçant on obtient donc La forme générale des solutions est un cercle si parcouru une infinité de fois, qui se déforme sous l’effet du champ électrique en une sorte de spirale de ressort, pour une vitesse initialle nulle, on obtient une cycloïde.
6.3.6 Allure des courbes en dimension 2.
Si on se place dans le repère propre (en prenant les vecteurs propres comme vecteurs de base), et si a deux valeurs propres distinctes ( est alors diagonalisable), alors chaque coordonnée suit une exponentielle, dans ce repère avec . Si et sont réels, l’une des exponentielles domine l’autre lorsque et c’est l’inverse lorsque , la courbe a pour asymptote les directions propres si et sont de signes contraires, ou tend vers 0 et possède une branche parabolique de direction asymptotique l’une des directions propres sinon. Si et sont complexes conjugués de partie réelle non nulle, on a une spirale qui tend vers 0 d’un coté et vers l’infini de l’autre (selon le signe de la partie réelle). Si est symétrique, alors et sont réels, ce cas ne peut pas se produire, de plus on peut choisir un repère propre orthonormé, les courbes ressemblent à des hyperboles. Ce sont des hyperboles si trace (la somme des valeurs propres vaut 0 donc le produit des coordonnées dans le repère propre vaut une constante), ces hyperboles sont équilatères si est symétrique. Quelques exemples :
Remarque :pour un système différentiel à coefficients non constants, il n’existe pas de méthode générale de résolution. Il arrive que dans certains cas particuliers, on puisse résoudre le système, par exemple si on trouve une matrice de passage indépendante du temps ramenant le système à un système diagonal ou triangulaire : un exemple avec Ou si commute avec , on peut prendre comme solution.
6.3.7 Systèmes d’ordre plus grand que 1
On se ramène à un système d’ordre 1. Par exemple deux ressorts couplés on pose , on a On délègue le calcul des valeurs propres à la machine :
6.3.8 Intégrales premières.
Lorsqu’on ne sait pas résoudre explicitement une équation ou un système différentiel, il peut arriver qu’on connaisse une ou des constantes du mouvement en cinématique, appelées aussi intégrales premières.
C’est le cas par exemple de l’énergie totale (mécanique plus cinétique) pour des forces conservatives. En dimension un, la connaissance de l’intégrale première énergie totale permet de ramener l’équation fondamentale de la dynamique d’ordre 2 à une équation du premier ordre à variables séparables : soit donc on peut ainsi calculer le temps en fonction de et tracer le graphe de en fonction de puis le graphe de en fonction de par symétrie par rapport à la première bissectrice.
Exemple : calcul de la période d’un pendule, on repère
une masse reliée à un fil de longueur à un point fixe
par l’angle formé avec la verticale (orienté vers le
bas), de sorte que
l’énergie potentielle de la masse est
on a donc
puis
Si on lache sans vitesse initiale la masse avec un angle
alors
donc
puis
Pour des raisons de symétrie, la période du pendule est donc
L’expression à intégrer n’admet pas de primitive avec les fonctions
usuelles, on l’appelle intégrale elliptique (il y a un lien avec
la longueur d’un arc d’ellipse). On peut calculer une valeur numérique
approchée de cette intégrale si est donné.
Pour de petites valeurs de , on peut approcher
par et calculer l’intégrale
qui ne dépend pas de . On observe que cette approximation
est encore assez bonne pour (erreur<4%).
En dimension plus grande, l’existence d’intégrales premières peut permettre de connaitre la forme de la courbe intégrale et même parfois de résoudre complètement l’équation (cas du problème à deux corps ci-dessous).
Autre exemple, la découverte d’un facteur intégrant pour la forme différentielle donne une intégrale première pour l’équation , en effet est nul sur une courbe intégrale, donc est constant, les courbes intégrales sont donc les courbes de niveau de . Une équation à variables séparables est un cas particulier, avec ne dépendant que de et de .
Pour un système autonome, est une intégrale première si grad, en effet
Problème à deux corps Cas d’un point de soumis à une force centrale comme la gravité ou la force coulombienne : on montre
- la conservation du moment cinétique (vérification immédiate en dérivant). Ceci entraine que le mouvement est dans un plan orthogonal à et la loi des aires (où est l’angle formé par avec une direction fixe du plan) : (ceci est vrai dès que la force est centrale, indépendamment de la norme de la force)
- la conservation du vecteur excentricité défini par : En effet où est orthonormé direct et .
Si on prend l’axe des porté par , en faisant le produit scalaire avec : on obtient en appliquant les propriétés du produit mixte et la définition de : la courbe intégrale est donc une conique d’excentricité ayant l’origine pour foyer et parcourue selon la loi des aires (l’aire balayée par le segment origine-point mobile est proportionnelle au temps).
6.3.9 Le modèle proie-prédateur
C’est un système autonome en dimension 2 pour lequel on sait calculer une intégrale première. Il se présente sous la forme avec des constantes positives, l’effectif des proies, celui des prédateurs, correspond à la reproduction naturelle des proies, à la mortalité par rencontre d’un prédateur, à la mortalité naturelle des prédateurs et à la natalité dépendant du nombre de proies. On peut déterminer les points d’équilibre et leur stabilité comme pour n’importe quel système autonome (exercice), on trouve qui est instable et , les valeurs propres du linéarisé sont 2 imaginaires purs conjugués, donc on ne peut pas conclure sur la stabilité à ce stade.
On peut déterminer une intégrale première en faisant apparaitre des dérivées logarthmiques donc en posant on a d’où : donc : est une intégrale première du mouvement, qui se passe donc sur les courbes de niveau de en ou de en . On observe que ces courbes de niveau sont fermées, impliquant un mouvement périodique, si on exprime en fonction de par le théorème des fonctions implicites donc sur toute la courbe à l’exception des deux points où la tangente est verticale alors on peut calculer la période du mouvement en appliquant : donc
6.3.10 Quelques autres méthodes
On peut encore citer : changement de fonction, changement de variables, équation homogène, équations de Bernoulli, de Clairault, de Ricatti, développements en séries entières..., certaines de ces méthodes sont implémentées par les logiciels de calcul formel.
6.4 Comportement asymptotique des solutions
Les équations de la physique sont souvent des équations autonomes sans second membre (pas de dépendance explicite en temps) ou avec un second membre qui est le seul terme de l’équation dépendant du temps (il s’agit d’un forçage extérieur). Dans le premier cas, les solutions doivent rester bornées (par exemple en énergie), donc ne peuvent pas tendre vers l’infini. Dans le second cas, une question naturelle est alors la suivante : le système atteint-il un équilibre, peut-on décomposer la solution en deux parties : un régime permanent et un régime transitoire ?
On a déjà fait une étude de comportement asymptotique pour l’équation , la solution se comporte comme un point déquilibre instable, si on en dévie même légèrement, on s’en éloigne définitivement, alors que se comporte comme un point déquilibre stable. Nous allons généraliser cette étude, pour les équations linéaires à coefficients constants (avec ou sans second membre, perturbation dépendant du temps), les équations autonomes sans second membre, et dans le cas de systèmes différentiels linéaires à coefficients constants.
6.4.1 Équations linéaires à coefficients constants d’ordre 1 et 2
Pour les équations homogènes d’ordre 1 , la solution générale est , le comportement asymptotique lorsque dépend du signe de , si la limite est 0 et la solution décroit exponentiellement vite. Donc si , quelle que soit la condition initiale, toutes les solutions de l’équation avec second membre ont le même comportement asymptotique, celui d’une solution particulière de l’équation :on a donc un régime transitoire exponentiellement décroissant et un régime permanent.
Pour les équations homogènes d’ordre 2 , la solution générale est si et sont les deux racines simples de ou si l’équation caractéristique admet une racine double. Le comportement à l’infini dépend du signe de la partie réelle de et . Il faut que les deux parties réelles soient strictement négatives pour que la solution tende vers 0, à vitesse exponentielle, si l’une au moins des parties réelles est positive ou nulle, alors il n’y a pas convergence vers 0. Plus précisément
- Si , il y a deux racines complexes conjuguées distinctes de partie réelle , donc la solution décroit exponentiellement vers 0 si , comme , avec des oscillations périodiques en , de période (régime oscillatoire amorti enveloppé par ). Si , la solution ne tend pas vers 0, reste de taille bornée, elle est périodique de période (régime oscillatoire)
- Si , on a deux racines réelles distinctes, qui sont toutes les deux strictement négatives si l’opposé de leur somme et leur produit sont positifs : (régime amorti équivalent à la plus grande des deux exponentielles)
- Si , on a une racine double , il y a convergence vers 0 si .
- Dans tous les autres cas, la partie réelle d’une des racines est positive ou nulle et il n’y a pas de convergence vers 0 de la solution générale. Si on a deux racines imaginaires pures conjuguées, la solution est périodique, sinon la solution tend vers l’infini pour une condition initiale générique.
Exemples
- , deux racines imaginaires pures conjuguées, solution générale périodique
- , deux racines complexes conjuguées de partie réelle négative, il y a convergence exponentielle vers 0 avec des oscillations, la solution générale est .
- , deux racines réelles, une positive, une négative. La solution générale est , elle tend génériquement vers l’infini (sauf condition initiale annulant ).
- ...
On peut généraliser à un ordre quelconque. Si toutes les racines de l’équation caractéristique sont de partie réelle négative, la solution générale de l’équation homogène tend vers 0 à l’infini, elle est appelée régime transitoire. Quelle que soit la condition initiale, on tend vers la solution particulière appelée régime permanent.
6.4.2 Forçage périodique
Il arrive souvent qu’un système physique soit soumis à un forçage extérieur périodique, par exemple pour la température à échelle fine, l’alternance jour-nuit, ou à grande échelle, l’alternance des saisons, ou circuit RCL soumis à un courant périodique. Il est donc utile de déterminer les caractéristiques de la solution en régime permanent.
Exemple : ordre 1 On sait qu’une solution particulière est donnée par , on remplace et on obtient L’amplitude de la solution particulière est donc l’amplitude du second membre divisée par le module , et l’exponentielle est en retard d’un déphasage donné par l’argument de soit . La solution particulière suit donc le second membre, avec un déphasage compris entre 0 et un quart de période, selon la valeur de . Si le système a une forte inertie intrinsèque ( grand pour avoir une solution de l’homogène associée en exponentielle décroissant lentement), on s’approche du quart de période, c’est pour cette raison que la température près de la mer atteint son maximum en été environ 2 mois après le solstice, alors que dans les terres, c’est plutot 3 semaines après (le maximum d’un quart de période étant presque réalisé par la banquise qui atteint son minimum d’extend presque 3 mois après le solstice).
À l’ordre 2, on peut faire la même étude pour , cette fois l’amplitude est divisée par Si (pas de frottements) et si est solution de l’équation caractéristique, la solution particulière est en , il y a résonance (c’est pour éviter d’entrer en résonance avec une fréquence propre d’un pont qu’on ne doit pas le traverser à plusieurs en marchant au même pas cadencé).
6.4.3 Équation autonome sans second membre
Il s’agit d’une équation de la forme où on suppose continument dérivable. Les solutions stationnaires sont données par les racines de (les telles que ). Pour toute condition initiale entre deux racines consécutive de , la solution va rester entre ces deux racines consécutives. Comme ne s’annule pas entre deux racines consécutives, est de signe constant donc la solution est monotone majorée ou minorée donc convergente, donc converge, donc converge vers 0 (sinon ne convergerait pas), donc tend vers une des racines lorsque 4. Si , on tend vers la plus grande des racines lorsque , sinon vers la plus petite. Si la condition initiale est au-delà de la plus grande racine ou en-deça de la plus petite racine, on tend soit vers l’infini, soit vers la racines.
On peut préciser la vitesse de convergence. Si , ( linéaire) la convergence vers se fait comme pour . Dans le cas général, si , ce résultat est encore valable, heuristiquement : où est une fonction qui tend vers 0 lorsque tend vers , donc : d’où le résultat (pour une justification plus rigoureuse il faut appliquer le théorème des fonctions implicites pour déterminer et vérifier que s’intègre).
Si , la solution est appelée équilibre stable : pour toute condition initiale situé entre et la solution tend vers lorsque et la convergence se fait à vitesse exponentielle, comme .
Exemple : pour l’équation logistique , , il y a deux équilibres et , avec et donc un équilibre stable en 1, et un équilibre instable en 0.
6.4.4 Systèmes linéaires
Cas linéaire
L’évolution du système est gouvernée par les valeurs propres
de la matrice du système, exactement comme pour les équations
linéaires où ce sont les racines de l’équation
caractéristique.
La solution générale tend vers 0 si toutes les valeurs propres
ont une partie réelle strictement négative. S’il y a des paires
de valeurs propres conjuguées de partie réelle négative,
des phénomènes cycliques amortis apparaissent.
Si les valeurs propres sont négatives ou nulles mais distinctes, la solution
reste bornée (avec des composantes qui peuvent être
périodiques).
Si une des valeurs propres a une partie réelle strictement positive,
alors pour une condition initiale générique, la solution tend vers
l’infini.
Exemples
-
[Y]:=desolve(y'=A*y and y(0)=[1,0,0])
pour
A:=[[1,2,3],[4,5,6],[7,8,9]]
, puisplot(Y[0],x=0..4)
la solution tend vers l’infini, à vitesse exponentielle comme on peut le voir avecplot(ln(Y[0]),x=0..4)
. En effet 16.12... est valeur propre de (eigenvalues(approx(A))
). On observe le même comportement en remplaçant par (ceci diffère de la dimension 1, où en changeant le sens du temps une solution divergente devient convergente). On peut représenter le graphe de la courbe décrite dans l’espace par exemple avecplotparam(Y,x=0..2)
[Y]:=desolve(y'=A*y and y(0)=[1,0])
pourA:=[[-3,1],[1,-5]]
, la courbe dans le plan est obtenue parplotparam(Y,x=0..10)
, en faisant plusieurs zoom out, on voit la courbe partir de la condition initiale le point et aboutir (presque) en l’origine. Les valeurs propres sont en effet .- Même chose avec
A:=[[-1,2],[-2,-1]];
. La courbe part toujours du point pour aboutir presque en l’origine, cette fois en spiralant (car les valeurs propres sont complexes conjuguées) - Pour
A:=[[0,2],[-2,0]];
, les valeurs propres sont imaginaires pures, la courbe est un cercle décrite de manière périodique.
Cas autonome
On ne sait pas intégrer un système sans plus
de précision sur (ce n’est plus une équation à
variables séparables et il n’y a pas d’ordre dans ,
donc pas de monotonie des solutions à attendre).
On ne peut donc
espérer un résultat général que si la condition initiale
est proche d’un point d’équilibre (une solution de ).
Dans la plupart des cas, on peut conclure sur la stabilité
ou l’instabilité du point déquilibre en fonction de la partie
réelle des valeurs propres de , un peu comme en dimension 1.
Si toutes les valeurs propres ont des parties strictement négative
on peut montrer que
le système revient à l’équilibre exponentiellement vite, si
l’une des parties réelles est strictement positive,
pour une condition initiale générique, le système s’en
éloigne, et s’il y a des parties réelles nulles, on ne peut pas conclure/
6.4.5 Forçage près d’un point d’équilibre de système.
Si on ajoute un terme dépendant du temps , on ne sait plus résoudre l’équation ni décrire son comportement qualitatif en toute généralité. Si la condition initiale est proche d’un équilibre stable, et si la perturbation est “petite” (en tenant compte de l’échelle de temps des exponentielles du système linéarisé) on peut alors linéariser et espérer que la solution se comporte comme la solution de au moins pendant un certain intervalle de temps. Si est petit et constant, le point d’équilibre est déplacé au premier ordre de
Exemple : modèle d’évolution température puis température-CO2.
Le modèle le plus simple ne tient compte que des radiations venues du Soleil et réémises par la Terre, considérée comme un corps noir où modélise l’inertie thermique, est la constante solaire (environ 1364/4) et est relié à la constante de Stefan-Boltzmann (5.67e-8 S.I.). On a alors un équilibre pour , et cet équilibre est stable.
Si on perturbe par un effet de serre additionnel du CO2, on modélise l’évolution de la température de la Terre par où est la température d’équiibre de la Terre et la concentration en ppm de gaz carbonique, modélise la capacité calorifique de la Terre (on peut estimer ), la constante de Stefan-Boltzmann (5.67e-8 S.I.). Par exemple avec un taux de CO2 stabilisé à 450ppm, le nouvel équilibre est donné à l’ordre 1 par
Le taux de CO2 de l’atmosphère peut être considéré comme un forçage extérieur (dépendant de scénarios d’émissions de CO2) mais il dépend aussi de la température de l’océan, on peut donc modéliser l’évolution conjointe des deux variables par un système différentiel autonome auquel on ajoute une composante dépendant du temps (émissions anthropiques). Par exemple un système 2 par 2 avec un second membre constant dans un scénario avec émissions de CO2 constantes. où est négatif (l’océan absorbe l’excédent de CO2 émis par rapport à la valeur avec laquelle il est en équilibre, on peut estimer par les observations : émissions 4.5 ppm par an, hausse de CO2 2ppm/an, donc ), et représente la perturbation anthropique (par exemple si stabilisation des émissions de CO2 à ce niveau). Dans un modèle simplifié ne dépend que de , la dérivée a des coefficients négatifs sur la diagonale et un coefficient nul sous la diagonale, donc les valeurs propres de sont négatives, le climat est stable. On atteint alors un nouvel équilibre avec une température et un taux de CO2 donnés par La valeur de la constante de couplage entre CO2 et T affecte évidemment le calcul de donc des valeurs à l’équilibre. Ici avec nos estimations :
On peut raffiner ce modèle en ajoutant par exemple la glace et ses interactions avec la température (si la température monte, la glace fond, si la glace fond, l’albédo de la Terre diminue ce qui va faire monter la température), ce qui amène à un système différentiel en dimension 3 où est une fonction décroissante, est positif, et représente la perturbation anthropique (la puissance deux tiers appliquée à la masse de glace sert à passer d’un volume à une surface pour représenter l’effet de la variation de volume de glace sur l’albédo).
- 1
- En toute rigueur, il faut prouver que la solution maximale est bien définie sur tout entier. Soit l’intervalle maximal de définition de la solution. Si , alors en intégrant qui est borné sur on obtient une valeur finie pour la limite en de , on peut alors prolonger autour de en appliquant le théorème de Cauchy-Lipschitz en , ce qui est contradictoire avec l’hypothèse de maximalité. Donc et de même
- 2
- On peut d’ailleurs démontrer que certaines équations ne sont pas résolubles de cette manière, cf. la théorie de Galois différentielle
- 3
- Cela vient du fait que les puissances de forment une famille d’un espace vectoriel de dimension finie , donc la famille est liée à partir de éléments, en fait on peut montrer que c’est le cas si on considère .
- 4
- On peut prouver l’existence globale de la solution exactement comme pour l’exemple de la section 6.2
Chapitre 7 Introduction au calcul variationnel
La recherche de minimas/maximas est une des application du calcul différentiel : en dimension 1, la dérivée s’annule lorsque la fonction est maximale ou minimale, en dimension plus grande c’est le gradient qui s’annule. Le calcul variationnel est une généralisation du principe précédent lorsque l’inconnue n’est pas l’endroit où l’extrêmum est atteint (un réel ou un point), mais une fonction . Par exemple, si on recherche le plus court chemin entre 2 points de l’espace, ou entre 2 points situé sur une sphère ou une surface : dans ce cas l’inconnue est le chemin, que l’on peut représenter par une courbe paramétrée. On obtient alors une équation différentielle qui permet de déterminer le chemin, de même que l’équation ou permettait de trouver la position d’un extrêmum. Réciproquement, certaines équations différentielles de la physique peuvent se mettre sous la forme minimiser une fonction dépendant d’un chemin, le chemin étant la courbe intégrale de l’équation différentielle. C’est le cas par exemple des équations de la dynamique en mécanique classique aussi bien qu’en relativité. Un des intérêts d’une formulation variationnelle de ces équations, c’est que ce type de formulation est plus intrinsèque (plus géométrique) elle ne dépend pas des coordonnées.
Dans le cas général on se donne :
- une fonction deux fois continument dérivable dépendant de la position dans un système de coordonnées, de la dérivée de par rapport au temps et du temps appelé lagrangien : En coordonnées cartésiennes , en coordonnées polaires , etc. En coordonnées cartésienne, est la vitesse1, alors qu’en polaires la vitesse n’est pas .
- et deux points
et on cherche parmi les courbes paramétrées deux fois continument dérivables d’origine et extrémité le(s) chemin(s) réalisant le minimum (s’il existe) de l’action2 : En coordonnées cartésiennes, est une courbe paramétrique (mais en coordonnées polaires, n’est pas une courbe en polaires).
Exemples :
- longueur minimale dans le plan, en coordonnées cartésiennes et En relativité générale, la longueur dépend du point (la structure de l’espace(-temps) est déformé par les masses), dépend de
- lagrangien de la mécanique classique ( ou 3, dans un référentiel galiléen en coordonnées cartésiennes) : énergie cinétique moins3 énergie potentielle :
- lagrangien en relativité restreinte : le premier terme est proportionnel à l’opposé du temps propre.
- si on considère une particule chargée de charge dans un champ électrique alors (dans un champ électro-magnétrique dépend de la vitesse, ), où est le potentiel scalaire (et où est le potentiel vecteur, on a rot). Les champs et peuvent dépendre du temps ou non.
Sur les exemples, on obtient
- pour la longueur minimale dans le plan, on a qui est la première composante du vecteur tangent, de même pour la deuxième composante, donc le long de la courbe le vecteur tangent a sa dérivée nulle, donc est constant. Une courbe réalisant un extrêmum de la distance entre deux points dans le plan est donc portée par une droite, c’est le segment reliant ces deux points.
- Pour le deuxième exemple c’est-à-dire la composante sur de la force et de la quantité de mouvement, donc l’équation d’Euler-Lagrange donne l’équation fondamentale de la dynamique.
- Pour le troisième exemple, on a c’est-à-dire la composante sur de la force et de la quantité de mouvement en relativité restreinte, on retrouve donc l’équation fondamentale de la dynamique.
Démonstration (idée) :
On fait varier le chemin en ajoutant à
un vecteur avec , on obtient
une action , on calcule la dérivée en
de , elle doit s’annuler pour avoir un extrêmum, et
ce quel que soit la valeur de la fonction
telle que . Prenons pour commencer
uniquement sur la première composante
, on a :
on dérive par rapport à sous le signe intégrale (on peut
intervertir dérivée et intégrale car sont
deux fois continument dérivables). Comme intervient dans deux
composantes de , il y a deux dérivées partielles qui
interviennent :
On intègre par parties le deuxième terme (), le terme tout intégré est nul
car , d’où :
Comme le résultat doit être nul
pour toute valeur de , on en déduit la première
équation d’Euler-Lagrange (en prenant
si la régularité
est suffisante, ou sinon en raisonnant par l’absurde :
si l’équation n’est pas vérifiée
en un point, alors on prend non nulle seulement
au voisinage de ce point et nulle ailleurs, et on choisit
de même signe que l’autre facteur, l’intégrale est alors
strictement positive, absurde).
Un des intérêts de cette écriture des équations de la
mécanique, c’est de pouvoir effectuer un changement de coordonnées
plus facilement, car la propriété de rendre l’action extrêmale pour un
chemin est indépendant du choix des coordonnées.
Exemple : si , on peut utiliser les coordonnées polaires
, on a alors
Si le potentiel dépend seulement de (en dimension 2), alors ne
dépend pas de (seulement de ) donc
on a donc une intégrale première, qui est le moment cinétique
. L’autre équation est
et s’exprime uniquement en fonction de
tout se passe comme si on était en dimension 1 avec un potentiel
effectif .
Exercice : Calculer le lagrangien en coordonnées sphériques
et donner les équations d’Euler-Lagrange si le potentiel
est radial ().
Solution abrégée
ne dépend pas explicitement de , il y a donc une constante
du mouvement
c’est le moment cinétique par rapport à .
ne dépend pas de mais en dépend, donc
n’est pas conservé :
Toutefois, pour des raisons de symétrie, les moments par rapport à
et sont aussi conservés, on a donc d’autres constantes
du mouvement. On peut continuer de deux manières, soit choisir
le repère pour avoir à la condition initiale, alors
reste nul pendant tout le mouvement qui se passe dans
le plan constant, on est ramené à un lagrangien en
coordonnées polaires, qui ne dépend plus de .
Ou bien on montre que
est constant.
Plus généralement, si ne dépend pas explicitement du temps, alors le hamiltonien défini par : est une intégrale première, en effet
Exercice : calculer pour le lagrangien de la mécanique classique et de la relativité restreinte.
Exemple : On cherche la forme d’un toboggan qui permette de se rendre le plus rapidement possible d’un point (origine du repère) à un point situé à une altitude plus basse sous l’action de la gravité (en négligeant les frottements). Si cette courbe est un graphe de fonction alors la vitesse est donnée par . D’autre part . Donc on en déduit : donc le temps à minimiser est Pour se ramener au problème précédent, on change de notations, devient un “temps virtuel” et est la dérivée de par rapport à ce temps virtuel, il faut minimiser le lagrangien ne dépend pas explicitement de , donc le hamiltonient correspondant est conservé, donc indépendant de donc en revenant à la notation pour l’abscisse on a Après simplification, on obtient l’équation différentielle : soit Comme et , on en déduit que est négatif : Il s’agit d’une équation à variables séparables. En posant on obtient une équation indépendante de : Donc puis (pour trouver la constante d’intégration, on observe que pour ) : Si on pose , on a , la solution est donc une cycloïde renversée. On peut aussi le vérifier directement en remplaçant dans l’équation et par les équations paramétriques de la cycloïde renversée on trouve pour le membre de droite :
Pour aller plus loin :
- le théorème de Noether permet de déduire une constante du mouvement lorsque le lagrangien est invariant par un groupe continu de transformations (par exemple translation ou rotation dans l’espace).
- les équations de la mécanique peuvent se réécrire en terme du hamiltonien, en utilisant les coordonnées de l’espace des phases, sous la forme On obtient un système symplectique qui possède des propriétés géométriques et permet de calculer des solutions numériques de meilleure qualité en en tenant compte (intégration symplectique). La formulation hamiltonienne de la mécanique est aussi très liée à la mécanique quantique (principe de correspondance).
7.1 Résumé équations différentielles et calcul variationnel
7.1.1 Équations et systèmes différentiels
- Cauchy-Lipschitz : avec continument dérivable, alors par toute condition initiale passe une courbe intégrale maximale et une seule.
- ordre 1, à variables séparables , la solution
s’obtient sous forme implicite
en intégrant .
Cas particulier les équations autonomes . - ordre 1, linéaire , solution sans second membre puis variation de la constante pour solution particulière
- Équation linéaire à coefficients constants
Base de solutions de l’équation sans second membre les où racine du polynôme caractéristique. Si est racine multiplicité , il faut ajouter à la base .
Avec second membre : si du type avec polynôme et réel ou complexe, il existe une solution du même type avec de même degré en général (ajouter au degré la multiplicité de comme racine du polynôme caractéristique). Pour un second membre général, utiliser la méthode de variation des constantes poser (si l’ordre est , poser pour ), remplacer dans l’équation différentielle puis résoudre le système en les . - Systèmes différentiels linéaires à coefficients
constants .
Sans second membre, si la matrice est diagonalisable
diag
alors .
Avec second membre, faire varier les constantes .
7.1.2 Comportement asymptotique des solutions
-
Équations linéaires à coefficients constants
Sans second membre : la solution tend vers 0 en l’infini si toutes les racines de l’équation caractéristique ont une partie réelle strictement négative. La décroissance est exponentielle (racine réelle de partie réelle maximale) ou exponentielle oscillante (paire de racines conjuguées de partie réelle maximale)
Avec second membre : un régime permanent s’établit, (régime transitoire exponentiellement décroissant apériodique ou oscillant). Si le second membre est périodique, le régime permanent est périodique de même période sauf s’il y a résonance avec une fréquence propre de l’équation ( racine du polynôme caractéristique). - Équation autonome de solutions stationnaires avec . Si une condition initiale est entre 2 solutions stationnaires consécutives, elle y reste et tend vers l’une en et vers l’autre en . Si , est un équilibre stable, si la condition initiale est proche de , alors tend vers pour
- Systèmes différentels linéaires à coefficients
constants de matrice
Sans second membre : la solution tend vers 0 en l’infini si toutes les valeurs propres de ont une partie réelle strictement négative. Décroissance exponentielle apériodique (valeur propre réelle de partie réelle maximale) ou oscillante (paire de valeurs propres conjuguées de partie réelle maximale).
Avec second membre : régime permanent, de même période si le second membre est périodique, sauf s’il y a résonance avec une fréquence propre du système ( valeur propre de ). - Système autonome : si est point d’équilibre (solution de ), et si les valeurs propres du linéarisé sont de partie réelle strictement négative, alors on montre que l’équilibre est stable.
7.1.3 Calcul variationnel
- lagrangien, système de coordonnées, pas forcément cartésien on peut par exemple prendre (coordonnées polaires) ou (coordonnées sphériques)
- action le long d’un chemin (courbe paramétrée) d’extrémités et
- condition nécessaire pour minimiser l’action, la courbe paramétrée doit vérfier les équations d’Euler-Lagrange
- Si ne dépend pas explicitement d’une des coordonnées alors est une constante du mouvement.
- Si ne dépend pas explicitement du temps alors le hamiltonien est une constante du mouvement
- En mécanique classique pour des forces conservatives
, énergie
cinétique, potentiel.
Si ne dépend pas d’une coordonnée cartésienne, la composante correspondante de la quantité de mouvement est conservée, si ne dépend que de , c’est le moment cinétique qui est conservé.
Le hamiltonien est l’énergie totale.
En relativité restreinte, est proportionnel au temps propre .
- 1
- on appelle espace des configurations l’espace des positions, vitesses, en coordonnées cartésiennes, espace que l’on a déjà utilisé implicitement pour passer d’une équation différentielle d’ordre 2 à un système différentiel d’ordre 1
- 2
- On pourrait bien sur avoir une action dépendant de dérivées d’ordre supérieur de
- 3
- Le signe moins vient de la convention adoptée en physique pour le lien entre potentiel et force
Chapitre 8 Utilisation de la calculatrice.
On présente ici un résumé des commandes utiles dans ce module pour les calculatrices empruntables au DLST (Casio Graph 90+e), ainsi que quelques pistes de documentation pour les étudiants possédant d’autres calculatrices (Numworks, HP Prime, TI Nspire, 89/92/Voyage 200, Casio Classpad, etc.).
8.1 Emprunt ou achat
Le DLST propose un emprunt de Casio Graph 90 aux étudiants pour le semestre ou l’année. Quelques informations pour les étudiants souhaitant acheter leur propre calculatrice parmi les modèles compatibles Xcas : la Casio Graph 35eii est la moins chère en monochrome (se trouve à moins de 50 euros à la rentrée avec l’offre de remboursement différé de Casio), la Casio Graph 90+e la moins chère en couleurs (peut se trouver 60 euros à cette rentrée), vient ensuite la Numworks (83 euros, calculatrice plus rapide). La HP Prime est la plus performante mais presque deux fois plus chère. En occasion, les TI Nspire CX et Numworks N0110 non verrouillées sont de bonnes affaires.
8.2 Reset
Aucun logiciel de calcul formel n’est exempt de bugs, surtout s’il n’a pas été beaucoup testé (cas des calculatrices couleurs récentes). Cela peut se traduire par un calcul qui n’en finit pas ou un plantage de la calculatrice. Localisez l’emplacement du bouton Reset qui permet de réinitialiser votre calculatrice : Casio Graph 90+e/35eii/Numworks/TI Nspire activable avec une pointe de stylo, HP Prime, activable avec une pointe de trombone ou une épingle.
8.3 Casio Graph 90+e/35eii
Si vous ne voyez pas l’icone de Xcas depuis MENU, installez-le en suivant la documentation de CAS Le lien précédent explique plus en détails comment utiliser Xcas sur cette calculatrice, on en donne ici un petit résumé.
8.3.1 Lancement
Taper sur la touche MENU, puis déplacez le curseur jusqu’à voir l’icone de Xcas et la légende Khicas ou Khicas90, puis tapez EXE. Khicas est l’ancienne version, plus testée, Khicas90 est la nouvelle version, plus complète.
Pour éteindre la calculatrice, tapez sur la touche MENU puis shift ON.
8.3.2 Interruption, nettoyage.
Pour interrompre un long calcul, tapez sur AC/ON.
En cas de crash (affichage SYSTEM ERROR...), essayez de taper sur la touche MENU puis ouvrez une autre application. Assez souvent, cela permet de sauvegarder la session en cours, on tape alors MENU, on relance CAS et on peut reprendre le travail sans perdre de données ni devoir réinitialiser la calculatrice.
En cas de crash persistant, il est conseillé
d’effacer les données de session de Xcas. Tapez sur la touche
MENU, puis Memoire, puis F2 (Memoire de stockage), effacez
le fichier session.xw
.
8.3.3 Remarques
Pour saisir une commande, cherchez-la dans les menus rapides (F1, F2, shift ou alpha F1 à F6) ou depuis le catalogue (touche F4). On peut aussi taper la commande en toutes lettres en bloquant le clavier en minuscules (F5). Avec Khicas90, vous pouvez taper shift INS pour saisir un caractère non présent sur le clavier.
L’unité d’angle par défaut est le radian, on peut le vérifier dans la ligne d’état (RAD). Si ce n’est pas le cas, faire F6, 7, EXIT. Vérifiez aussi que vous êtes en mode Xcas, si Python apparait dans la ligne d’état, faire shift, SETUP, 2, 1.
Pour obtenir une approximation numérique d’une expression, on
peut utiliser
la commande approx(
, raccourci clavier touches shift F2 2.
8.3.4 Courbes
Pour tracer la représentation graphique d’une courbe, on peut utiliser l’application intégrée de Casio (cf. le manuel du constructeur), mais il n’est alors pas possible de faire l’étude analytique de la courbe On conseille donc de rester dans CAS et d’utiliser
-
pour un graphe de fonction la commande
plot
, - pour une courbe en paramétriques la commande
plotparam
, - pour une courbe en polaires la commande
plotpolar
.
Ces commandes se trouvent dans le menu rapide shift F4 ou depuis le menu complet dans le sous-menu Graphiques : touches F4 puis 7 sélectionner la commande avec le curseur haut ou bas. Taper ensuite F2 pour recopier un exemple en ligne de commande, vous pouvez ensuite modifier cet exemple et l’exécuter.
La touche permet de saisir ou selon le réglage du menu shift-SETUP.
Par exemple
X:=sin(2t);
Y:=cos(3t);
plotparam([X,Y],t,0,2*pi)
.
Pour saisir X
, taper ALPHA puis la touche + (légende X en rouge) ;
pour saisir :=
, taper shift INS, , etc.
Pour définir une variable, on peut omettre le :
, i.e. écrire =
à la place de :=
Utilisez les touches de déplacement du Pad pour bouger le graphique, les touches + et - pour faire un zoom in ou out et EXIT pour quitter.
Pour faire les calculs nécessaires à l’étude de la courbe,
vous pouvez utiliser les fonctions de calcul formel (factor
,
simplify
, solve
, diff
, limit
,
int
, etc.) qui se trouvent dans les menus rapides F1/F2 ou
par le menu complet (F4, puis Algebre
ou Analyse ou Resoudre) ou via des raccourcis claviers
(* pour factoriser, + pour développer).
Si vous tracez une courbe, les variables X1, X2, Y1, Y2
contiendront automatiquement les composantes de la vitesse et de l’accélération.
Par exemple X1
contient diff(X,t)
.
Vous pouvez aussi avec le raccourci (F2, 2) calculer la dérivée de
par rapport à .
X1=>*
factorise la dérivée
(touche pour saisir =>
),
solve(X1=0,t)
cherche les zéros de la dérivée, etc.
Pour étudier une courbe en polaire, on peut utiliser les complexes, par
exemple F:=R(x)*exp(i*x) stocke dans F
l’expression de l’affixe du point, on peut ensuite calculer F(x=)
pour avoir l’affixe du point d’angle ou
F’(x=) pour l’affixe du vecteur directeur
de la tangente en ce point.
Les calculs d’intégrales (F4 Analyse integrate) se font par défaut en
cherchant une primitive ce qui peut être long ou/et ne pas aboutir,
vous pouvez forcer le calcul approché d’une intégrale
définie en mettant une des bornes sous forme d’un nombre approché
par exemple 1.0
au lieu de 1
.
Sur les Graph 90, la commande integrate
permet de calculer une intégrale curviligne
(cf. la section 5.2).
8.3.5 Équations différentielles
Pour résoudre une équation ou un système différentiel linéaire,
ouvrez le catalogue F4 puis Resoudre puis la commande
desolve(
, puis ses arguments comme dans Xcas, par
exemple desolve(Y'=x*Y-x,x,Y)
ou avec condition
initiale desolve([Y'=x*Y-x,Y(0)=2],x,Y)
.
Pour saisir le caractère prime de dérivation, taper F2 1.
Pour les systèmes différentiels linéaires à coefficients
constants, la syntaxe est similaire
par exemple desolve(Y'=[[1,-1],[2,4]]*Y,x,Y)
ou avec condition
initiale desolve([Y'=[[1,-1],[2,4]]*Y,Y(0)=[1,2]],x,Y)
.
Pour saisir une matrice, vous pouvez l’entrer directement comme liste de listes en ligne de commande ou utiliser le menu F6, Editer matrice (ou le raccourci shift-Mat).
Pour vérifier des étapes intermédiaires des calculs,
vous pouvez utiliser, depuis les menus rapides shift-Mat ou shift-3 ou depuis
le sous-menu Algebre (bi)lineaire
, les commandes det
(déterminant), charpoly
(polynôme caractéristique),
eigenvals
(valeurs propres), eigenvects
(vecteurs
propres), jordan
(diagonalisation). Par exemple
A:=[[1,-1],[2,4]]
factor(det(A-x*idn(2)))
eigenvals(A)
eigenvects(A)
P,D:=jordan(A); P*D*inv(P)
Pour tracer un champ des tangentes, choisissez plotfield
dans le sous-menu Graphique
. Par exemple
plotfield(sin(t*Y),[t=-3..3,Y=-3..3],plotode=[[0,1],[0,-1]])
trace le champ des tangentes pour l’équation différentielle
avec les solutions passant par les conditions initiales
et .
8.4 Numworks N0110/N0115/N0120, TI Nspire CX et CX II
Vous pouvez installer CAS sur les calculatrices :
- Installation Numworks : Page d’installation et kit de connexion Numworks
- TI Nspire CX et CX2, avec l’OS 4.5.0 ou
4.5.3 ou 4.5.4 (Nspire CX) ou bien 5.2 ou 5.3 (Nspire CX 2). Tapez sur
la touche Home puis 4 et 5 pour connaitre votre numéro de version.
Attention ne faites pas de mise à jour sur le site de TI, cela
rendrait votre calculatrice irréversiblement incompatible avec
Xcas.
Installation : Si vous devez installer une version de l’OS compatible, suivez les instructions du début du tutoriel d’installation de ndless pour CX ou CX 2.
Ensuite téléchargez et désarchivez le fichier CX ou CX II. Envoyez le dossierndless
de l’archive avec le logiciel de connexion sur votre calculatrice, puis installez Ndless en lançant dans le dossierndless
l’installeur correspondant à votre OS, puis lancezkhicas
.
Le mode de fonctionnement est très proche de celui des Casio Graph 90, il y a quelques différences pour les raccourcis clavier (les menus rapides se lancent en tapant shift suivi de 1 à 9, ou 0 ou les parenthèses, le point décimal,...).
8.5 HP Prime
Pour un guide plus complet, vous pouvez consulter pour le calcul formel
http://www-fourier.ujf-grenoble.fr/~parisse/calc/hprime.pdf
pour le reste le guide du constructeur
https://tiplanet.org/forum/archives_voir.php?id=19716
Emulateur gratuit disponible sur hpcalc.org, choisir “HP Prime Virtual Calculator Emulator” pour son système
d’exploitation.
8.5.1 Nettoyage.
Pour effacer les variables de CAS, tapez la commande restart
(disponible depuis la touche boite à outils, catalog).
En cas d’anomalies persistantes après la commande restart
,
vous pouvez tout effacer. Pour cela éteindre la calculatrice. Appuyez
simultanément sur les touches F, O et C, laissez ces 3 touches
enfoncées et appuyer sur ON, relachez ON en laissant F, O et C
enfoncées. Relachez F, O et C quelques secondes plus tard.
Sélectionnez 4. FLS Utilities
, puis FORMAT C:
,
tapez ensuite sur la touche Esc puis sur 9 RESET
.
8.5.2 Remarques
Vérifiez que les angles sont en radians (touche Settings).
Pour accéder à une commande depuis CAS
,
on peut utiliser l’arborescence des
menus depuis la touche représentant une boite à outils,
ou la liste des commandes
(Ctlg
) ou taper la commande en toutes lettres. La touche Help
affiche de l’aide sur la commande saisie.
Pour taper un caractère alphabétique, taper auparavant sur
la touche ALPHA
orange. En tapant deux fois sur ALPHA
on verrouille le mode alphabétique (taper une troisième fois pour
déverouiller). Pour obtenir un caractère majuscule dans CAS, faire
précéder la touche du caractère souhaité de la touche
Shift
.
Attention pour faire du calcul formel, tapez sur la touche CAS,
et utilisez des noms de
variables en minuscules, n’utilisez pas les variables
A
à Z
ni .
Laissez le mode de calcul par défaut en mode exact (touche CAS Settings),
si vous voulez
une approximation numérique, il suffit de saisir une des données
avec un point décimal (par exemple sin(1.0)
au lieu
de sin(1)
).
Pour commencer un nouvel exercice, vous pouvez faire le ménage en
utilisant la commande restart
dans l’écran CAS
.
8.5.3 Courbes.
Appuyer sur la touche Apps, sélectionner sur l’écran tactile
Paramétrique ou Polaire (ou Fonction).
Entrez X1(T)
et Y1(T)
ou
R1(), vous pouvez vous aider de la touche
1
pour taper la variable et de la touche à sa gauche pour entrer des
fractions, racines carrées ... Pour voir le graphe, appuyer sur la
touche Plot
. Pour changer le zoom ou le cadrage, utiliser les doigts sur
l’écran tactile comme sur un smartphone ou une tablette.
Pour revenir à la définition de fonction
appuyer sur la touche Symb
. Pour changer la discrétisation,
appuyer sur Shift
puis Plot_Setup
.
Pour faire les calculs nécessaires à l’étude d’une courbe,
tapez sur la touche CAS. Les
fonctions de calcul formel (factor
,
simplify
, solve
, diff
, limit
,
int
, etc.) se trouvent dans le menu
qu’on ouvre en tapant sur la touche boite à outils (B en orange),
puis appuyer sur le mot CAS en bas de l’écran tactile.
Par exemple, si vous avez entré un graphe en paramétriques dans
, vous pouvez stocker la dérivée de dans la
variable v1
en tapant
diff(X1(t),t)
(pour saisir X
en majuscules,
il faut taper la touche
alpha puis la touche shift puis la touche X orange)
puis Sto sur l’écran tactile (ou touches
Shift
puis EEX_Sto
)
puis v1
, puis factor(v1)
ou solve(v1=0,t)
(attention à bien préciser la variable recherchée dans solve
si ce n’est pas x
),
de même pour dy1
, etc.
Pour une courbe en polaires, on peut utiliser les complexes, par
exemple r:=R1(x)*exp(i*x) stocke dans r
l’expression de l’affixe du point, on peut ensuite calculer r(x=)
pour avoir l’affixe du point d’angle ou
r’(x=) pour l’affixe du vecteur directeur
de la tangente en ce point.
Les calculs d’intégrales se font par défaut en
cherchant une primitive ce qui peut être long ou/et ne pas aboutir,
vous pouvez forcer le calcul approché d’une intégrale
définie en mettant une des bornes sous forme d’un nombre approché
par exemple 1.0
au lieu de 1
.
Pour visualiser des courbes implicites, on peut utiliser l’App
Graphique Avancé ou l’App Géométrie (Tracé avancé)
ou directement dans CAS, instruction
plotimplicit
. Pour les courbes
de niveau, utiliser dans CAS ou l’App Géométrie la
commande plotcontour
.
8.5.4 Équations différentielles
Pour résoudre une équation ou un système différentiel linéaire,
appuyer sur la touche CAS si vous n’êtes pas dans la vue CAS,
puis touche boite à outils, menu CAS, Résoudre,
Equation différentielle, ou taper directement la commande
desolve(
, puis ses arguments comme dans Xcas, par
exemple desolve(y'=x*y-x,x,y)
ou avec condition
initiale desolve(y'=x*y-x and y(0)=2,x,y)
.
Pour saisir le caractère prime de dérivation, vous pouvez taper
sur shift () et supprimer un des 2 primes.
Pour les systèmes différentiels linéaires à coefficients
constants, la syntaxe est similaire
par exemple desolve(y'=[[1,-1],[2,4]]*y)
ou avec condition
initiale desolve(y'=[[1,-1],[2,4]]*y and y(0)=[1,2])
.
Pour avoir des détails, vous pouvez utiliser les commandes det
(déterminant), charpoly
(polynôme caractéristique),
eigenvals
(valeurs propres), eigenvects
(vecteurs
propres), jordan
(diagonalisation). Par exemple
a:=[[1,-1],[2,4]]
factor(det(a-x*identity(a)))
eigenvals(a)
eigenvects(a)
p,d:=jordan(a); p*d*inv(p)
Pour tracer un champ des tangentes, vérifiez que x
et y
ne sont pas affectées (purge(x,y)
si nécessaire),
ouvrir l’application
Géométrie, puis dans la vue par défaut Plot, menu Cmds,
Tracé, puis champ de directions.
Pour tracer la solution approchée d’une équation différentielle
avec condition initiale, déplacez le curseur dans la
fenêtre Plot
au point de coordonnées
la condition initiale et tapez Entree. Attention, le calcul
est à nouveau effectué à chaque rafraichissement de l’écran
graphique, ce calcul peut devenir très long dès qu’on trace
plus de 2 courbes. Vous pouvez effacer une courbe intégrale
depuis Symb
.
On peut aussi choisir ODE au lieu de champ de directions
mais faut alors saisir
les arguments de la commande plotode
, pour vous aider à le faire,
tapez sur la touche Help
, recopiez un exemple (touche menu Ex)
et adaptez-le.
8.6 TI 89, 92, Voyage 200
Guide complet du constructeur
https://tiplanet.org/forum/archives_voir.php?id=609
8.6.1 Nettoyage.
Pour nettoyer, depuis HOME, taper sur la touche F6 Cleanup,
vous pouvez effacer soit les variables à 1 lettre (a
à z
)
soit tout effacer (NewProb
)
8.6.2 Remarques
Vérifiez que les angles sont en radians (touche MODE).
Pour saisir une commande dans l’écran HOME
,
on peut utiliser les menus F2, etc.
ou la liste des commandes par ordre alphabétique
touche CATALOG
, la syntaxe des commandes s’y affiche, ou
taper la commande en toutes lettres.
Pour taper un caractère alphabétique sur les TI89, taper auparavant sur
la touche ALPHA
violette. En tapant deux fois sur ALPHA
on verrouille le mode alphabétique (taper une troisième fois pour
déverouiller).
Attention, les calculatrices TI font la distinction entre le -
binaire de la soustraction et le -
unaire de changement de signe,
ceci peut engendrer des erreurs.
Il est conseillé de conserver le MODE
de calcul AUTO
,
les calculs se feront par défaut en mode exact, si vous voulez
une approximation numérique, il suffit de saisir une des données
avec un point décimal (par exemple sin(1.0)
au lieu
de sin(1)
).
Pour commencer un nouvel exercice, vous pouvez faire le ménage en
utilisant le menu F6
de HOME
.
8.6.3 Courbes
Touche MODE, sélectionnez le type de tracé paramétrique ou
polaire (ou fonction). Puis touche verte, F1 (Y=), définir les
fonctions xt1
et yt1
(attention les TI font la
différence entre le -
unaire opposé et le -
binaire soustraction et il faut explicitement fermer les parenthèses ouvertes).
Puis touche verte, F2 (WINDOW) pour définir la discrétisation
et le cadrage graphique, puis touche verte, F3 (GRAPH) pour avoir le
tracé. On peut ensuite zoomer avec le menu F2.
Pour faire les calculs nécessaires à l’étude d’une courbe,
tapez sur la touche HOME.
Les fonctions de calcul formel (factor
,
solve
, d
, limit
(par exemple
limit(1/x,x,0,1)
pour la limite de en ),
, etc.) se trouvent dans les menus F2 et F3.
Par exemple, si vous avez entré un graphe en paramétriques dans
, vous pouvez stocker la dérivée de dans la
variable v1
en tapant
(xt1(t),t)
touche Sto
puis
v1
, puis factor(v1)
ou solve(v1=0,t)
, de même
pour dy1
, etc.
Pour une courbe en polaires, on peut utiliser les complexes, par
exemple r:=r1(x)*exp(i*x) stocke dans r
l’expression de l’affixe du point, on peut ensuite calculer r|x=
pour avoir l’affixe du point d’angle ou
(r,x)|x= pour l’affixe du vecteur directeur
de la tangente en ce point.
Les calculs d’intégrales se font par défaut en commençant
par chercher une primitive ce qui peut être long ou/et ne pas aboutir,
vous pouvez forcer le calcul approché d’une intégrale
définie en mettant une des bornes sous forme d’un nombre approché
par exemple 1.0
au lieu de 1
.
Pour visualiser des courbes de niveau ou implicite, touche MODE puis
3d, Enter puis touche verte, Y=
, puis menu F1, puis 9, puis Style.
8.6.4 Équations différentielles
Pour résoudre une équation différentielle, vous pouvez utiliser
la commande deSolve
(menu F3-Calc de HOME), par exemple
deSolve(y'=x*y-x,x,y)
. Pour taper le prime de dérivation,
utiliser la touche 2nd puis = sur les 89 ou B sur les 92/V200.
S’il y a une condition initiale, on utilise and
, par exemple
deSolve(y'=x*y-x and y(0)=2,x,y)
.
Il n’y a pas de commande permettant de résoudre directement un
système différentiel. Mais on peut calculer le polynôme
caractéristique d’une matrice avec la commande det
, par
exemple
[[1,-1][2,4]]
sto a
det(a-x*identity(2))
sto p
solve(p=0,x)
Les commandes eigVc(a)
et eigVl(a)
renvoient les
valeurs approchées des vecteurs propres et valeurs propres d’une
matrice numérique a
, mais elles ne fonctionnent pas pour
une matrice à coefficients symboliques ou si on veut un résultat
exact. Pour chaque valeur propre
renvoyée (ici on prend 2) par la commande solve
,
on peut utiliser la commande
rref(a-2*identity(a))
pour réduire sous forme
échelonnée le système linéaire à résoudre pour trouver
les vecteurs propres. La résolution du système est alors triviale.
Pour reprśenter un champ de tangentes, taper sur la touche MODE
et sélectionner DIFF EQUATION pour Graph. Ensuite, touche verte
Y= pour entrer l’équation différentielle (t0
est l’instant de la
condition initiale, 0 par défaut, yi1
est la valeur de y1
à
l’instant t0
), touche verte WINDOW
pour régler la fenêtre, touche verte GRAPH pour tracer
le champ des tangentes et la solution de l’équation différentielle
passant par la condition initiale donnée. On peut tracer plusieurs
solutions simultanément en donnant une liste pour yi1
,
par exemple {1,2}
tracera les deux solutions passant par les
deux conditions initiales et .
8.7 Autres calculatrices
D’autres modèles graphiques formels sont partiellement documentés sur
http://www-fourier.ujf-grenoble.fr/~parisse/calc/calc.html
-
Casio Classpad, voir aussi
https://tiplanet.org/forum/archives_cat.php
ouhttp://www.casio-education.fr
- TI Nspire (CX) CAS, voir aussi
https://tiplanet.org/forum/archives_list.php?cat=Manuels+Nspire
- Tracé de courbes : lancer Graphes depuis le menu de départ puis touche Menu, 3. Graph puis 3. Parametric. On peut aussi régler le zoom depuis Menu.
- Représentation d’équations différentielles :
lancer Graphes depuis le menu de départ
puis touche Menu, 7. Equa diffs. Donner ensuite
y1'
en fonction dex
ety1
. Si vous ne précisez pas de condition initiale, le champ des tangentes est tracé, sinon la solution passant par la condition initiale. - Calcul formel : revenir à l’écran de départ (touche ON/home sur les CX), puis A Calculate. Utiliser la touche Menu pour saisir des commandes de calcul formel ou taper directement au clavier, les commandes sont les mêmes que pour les TI89/92 décrites ci-dessus.
- HP40G, HP49G, HP50G, voir
http://www-fourier.ujf-grenoble.fr/~parisse/#calc
La Casio Graph 100 dispose aussi d’un module de calcul formel
(un peu moins performant), vous pouvez télécharger son manuel
depuis
www.support.casio-europe.com/fr/download/manuals/calc/GRAPH100_MAN1_Fr.pdf
Modèles graphiques sur lesquels on peut ajouter un logiciel de calcul formel (moins performant) :
-
Casio Graph 75, 85, 95 (fx9860G*) : on peut installer Algebra Eigenmath
https://tiplanet.org/forum/archives_voir.php?id=448883
- Casio Graph 35+USB (fx9750Gii, ceci ne concerne pas les
Graph 35+eii où on peut installer KhiCAS) :
il faut d’abord transformer la calculatrice en Casio Graph 75 ou 85 ce qui
est un peu compliqué mais bien expliqué ici :
https://tiplanet.org/forum/archives_voir.php?id=4489
puis on peut installer Eigenmath fx (item précédent). - TI Nspire (non CAS) : on peut installer Ndless puis CAS
(prononcer Khicas, le portage de Xcas pour TI Nspire) :
http://www-fourier.ujf-grenoble.fr/~parisse/install_fr.html#ti