kinematics_calculus.tex

<< Sujets Mathématiques\ de la 3D graphique

Sommaire

Mécanique 2 : Dynamique Linéaire et Rotationnelle >>

Chapitre 11

Mécanique 1 : Cinématique Linéaire et Calcul

Ne confondez jamais le mouvement avec l’action.

— Ernest Hemingway (1899–1961)

« Mesdames et messieurs, puis-je attirer votre attention sur la piste centrale. Contemplez devant vous deux manuels ordinaires, l’un intitulé Physique Générale et l’autre Calcul Intégral. Leurs 2 500 pages combinées pèsent plus de 10 kg. Pourtant, dans ce chapitre et le suivant, vos courageux auteurs-cascadeurs vont tenter le spectacle le plus vertigineux et impossible de mysticisme et de subterfuge : réduire ces deux massifs ouvrages en à peine 150 pages ! »

Comme tout bon numéro de cirque, celui-ci est précédé de beaucoup d’annonces pour calibrer vos attentes. La différence ici est que l’objectif de notre préambule est d’abaisser vos attentes.

11.1Vue d’ensemble et autres remarques visant à réduire vos attentes

D’accord, il n’y a aucun moyen de couvrir toute la physique et le calcul en deux chapitres. Comme tout bon politicien le sait, le secret pour communiquer efficacement une matière complexe en peu de temps est d’utiliser des mensonges, tant par omission que par commission. Parlons de chacun de ces types de mensonges afin que vous sachiez ce qui vous attend réellement.

11.1.1Ce qui est laissé de côté

À peu près tout — commençons par ce que nous laissons de côté en physique. Mettre le mot « physique » sur ce chapitre serait encore plus insultant pour les personnes qui font de la vraie physique que ce chapitre ne l’est déjà. Nous nous intéressons uniquement à la mécanique, et à une mécanique très simple de corps rigides qui plus est. Parmi les sujets traditionnellement abordés dans un manuel de physique de première année qui ne sont pas traités dans ce livre, on trouve :

Une remarque s’impose concernant l’énergie et le travail, car même dans le contexte limité de la mécanique, le concept fondamental d’énergie joue un rôle central dans les présentations traditionnelles. De nombreux problèmes sont plus faciles à résoudre en utilisant la conservation de l’énergie qu’en considérant les forces et en appliquant les lois de Newton. (En fait, il existe une alternative à la dynamique newtonienne que nous étudions dans ce livre. Elle est connue sous le nom de dynamique lagrangienne et se concentre sur l’énergie plutôt que sur les forces. Utilisés correctement, les deux systèmes produisent les mêmes résultats, mais la dynamique lagrangienne peut résoudre certains problèmes de manière plus élégante et est particulièrement habile pour traiter la friction, comparée à la dynamique newtonienne.) Cependant, au moment de la rédaction, les simulations numériques de base à usage général sont basées sur la dynamique newtonienne, et l’énergie n’y joue pas de rôle direct. Cela ne signifie pas qu’une compréhension de l’énergie soit inutile ; en fait, la violation de la loi de conservation de l’énergie est au cœur de nombreux problèmes de simulation ! Ainsi, l’énergie apparaît souvent davantage comme un moyen de comprendre le comportement (défectueux) d’une simulation numérique, même si elle n’apparaît pas directement dans le code de simulation.

Parlons maintenant des façons dont ce livre va irriter les professeurs de calcul. Nous pensons qu’une compréhension basique du calcul est vraiment importante pour saisir pleinement de nombreux concepts de physique. Inversement, la physique fournit certains des meilleurs exemples pour expliquer le calcul. Le calcul et la physique sont souvent enseignés séparément, généralement avec le calcul en premier. C’est notre avis que cela rend le calcul plus difficile à apprendre, car cela prive l’étudiant des exemples les plus intuitifs — les problèmes de physique pour lesquels le calcul a été inventé ! Nous espérons qu’entrelacer le calcul avec la physique facilitera votre apprentissage du calcul.

Nos besoins en calcul sont extrêmement modestes dans ce livre, et nous avons laissé de côté encore plus de calcul que de physique. Après avoir lu ce chapitre, vous devriez savoir :

Bien sûr, nous savons qu’un grand nombre de lecteurs possèdent déjà ces connaissances. Prenez un moment pour vous classer dans l’une des catégories suivantes :

  1. Je ne sais absolument rien sur les dérivées ou les intégrales.

  2. Je connais l’idée de base des dérivées et des intégrales, mais je serais probablement incapable de résoudre un problème de calcul de première année avec un crayon et du papier.

  3. J’ai étudié du calcul.

Le niveau 2 de connaissance du calcul est suffisant pour ce livre, et notre objectif est de faire passer en catégorie 2 tous ceux qui sont actuellement en catégorie 1. Si vous êtes en catégorie 3, nos discussions sur le calcul seront une révision (espérons-le divertissante). Nous n’avons pas la prétention de pouvoir faire passer en catégorie 3 quelqu’un qui n’y est pas déjà.

11.1.2Quelques mensonges utiles sur notre univers

L’univers est communément considéré comme discret à la fois dans l’espace et dans le temps. Non seulement la matière est divisée en morceaux discrets appelés atomes, mais il existe des preuves que le tissu même de l’espace et du temps est également divisé en pièces discrètes. Il existe une différence d’opinion quant à savoir si c’est vraiment ainsi ou si cela apparaît seulement de cette façon parce que la seule façon dont nous pouvons interagir avec l’espace est d’y lancer des particules, mais notre avis est que si ça ressemble à un canard, marche comme un canard, cancane comme un canard, a des pattes palmées et un bec, alors c’est une bonne hypothèse de travail que ça a bon goût dans des rouleaux de printemps avec un peu de sauce sombre.

Pendant longtemps, la simple idée que l’univers pourrait ne pas être continu n’avait même pas effleuré l’esprit de quiconque, jusqu’à ce que les Grecs anciens aient l’idée saugrenue et totalement injustifiée que les choses pourraient être faites d’atomes. Le fait que cela se soit avéré exact plus tard est considéré par beaucoup comme relevant davantage de la chance que du bon jugement. Franchement, qui l’aurait cru ? Après tout, les objets du quotidien, comme le bureau sur lequel l’un des auteurs pose ses poignets pendant qu’il tape cette phrase, donnent toute l’apparence d’avoir des surfaces lisses et continues. Mais qui s’en soucie ? Penser au bureau comme ayant une surface lisse et continue est une illusion inoffensive mais utile qui permet à l’auteur de poser ses poignets confortablement sans se préoccuper du tout de l’énergie de liaison atomique et de la théorie de l’incertitude quantique.

Non seulement cette astuce consistant à penser le monde comme continu est une rationalisation psychologique pratique, mais c’est aussi de bonnes mathématiques. Il s’avère que les mathématiques des choses continues sont beaucoup moins maniables que les mathématiques des choses discrètes. C’est pourquoi les personnes qui réfléchissaient au fonctionnement du monde au XVe siècle étaient heureuses d’inventer une mathématique pour un univers continu ; expérimentalement, c’était une bonne approximation de la réalité, et théoriquement les mathématiques fonctionnaient bien. Sir Isaac Newton a ainsi pu découvrir de nombreux résultats fondamentaux sur les mathématiques continues, que nous appelons « calcul », et son application à l’exploration d’un univers continu, que nous appelons « physique ».

Maintenant, nous faisons principalement cela pour pouvoir modéliser un monde de jeu dans un ordinateur, qui est intrinsèquement discret. Il y a une certaine dissonance cognitive à programmer une simulation discrète d’un modèle continu d’un univers discret, mais nous allons essayer de ne pas nous en préoccuper. Il suffit de dire que nous contrôlons entièrement l’univers discret à l’intérieur de notre jeu, ce qui signifie que nous pouvons choisir le type de physique qui s’applique dans cet univers. Tout ce dont nous avons besoin, c’est que les lois physiques ressemblent suffisamment à celles auxquelles nous sommes habitués pour que le joueur éprouve une suspension consentie de l’incrédulité, et espérons-le dise : « Wow ! Cool ! » et veuille dépenser plus d’argent. Pour presque tous les jeux, cela signifie un univers newtonien douillet sans les détails désagréables de la mécanique quantique ou de la relativité. Malheureusement, cela signifie aussi qu’il y a deux trolls désagréables qui rôdent sous le pont, répondant aux noms de chaos et d’instabilité, mais nous ferons de notre mieux pour les apaiser.

Pour l’instant, nous nous intéressons au mouvement d’un petit objet appelé « particule ». À tout moment donné, nous connaissons sa position et sa vitesse.1 La particule a une masse. Nous ne nous préoccupons pas de l’orientation de la particule (pour l’instant), et donc nous ne la considérons pas comme tournant sur elle-même. La particule n’a pas de taille non plus. Nous différerons l’ajout de ces éléments à plus tard, quand nous passerons des particules aux corps rigides.

Nous étudions la mécanique classique, également connue sous le nom de mécanique newtonienne, qui comporte plusieurs hypothèses simplificatrices incorrectes en général mais vraies dans la vie quotidienne dans la plupart des aspects qui nous importent. Nous pouvons donc nous assurer qu’elles sont vraies dans notre monde informatique, si nous le souhaitons. Ces hypothèses sont :

Les deux premières sont brisées par la relativité, et les deux dernières par la mécanique quantique. Heureusement, ces deux sujets ne sont pas nécessaires pour les jeux vidéo, car vos auteurs n’en ont qu’une compréhension superficielle.

Nous allons commencer notre incursion dans le domaine de la mécanique en apprenant la cinématique, qui est l’étude des équations décrivant le mouvement d’une particule dans des situations simples mais courantes. Lors de l’étude de la cinématique, nous ne nous préoccupons pas des causes du mouvement — c’est le sujet de la dynamique, qui sera couverte dans le Chapitre 12. Pour l’instant, « ce n’est pas à nous de demander pourquoi », il s’agit juste de faire les mathématiques pour obtenir des équations prédisant la position, la vitesse et l’accélération de la particule à tout moment tt, ou de mourir. Enfin, oublions cette dernière partie.

Parce que nous traitons nos objets comme des particules et que nous ne suivons que leur position, nous ne considérerons pas leur orientation ni les effets rotationnels avant le Chapitre 12. Lorsque la rotation est ignorée, toutes les idées de la cinématique linéaire s’étendent en 3D de manière simple, et pour l’instant nous nous limiterons à la 2D (et à la 1D). C’est pratique, car les auteurs ne savent pas comment concevoir ces petites choses origami qui s’aplatissent et surgissent quand on ouvre le livre, et l’éditeur ne nous laisserait pas faire même si nous étions assez compulsifs pour apprendre.

11.2Grandeurs de base et unités

La mécanique traite des relations entre trois grandeurs fondamentales de la nature : la longueur, le temps et la masse. La longueur est une grandeur que vous connaissez sans aucun doute ; nous la mesurons en centimètres, en pouces, en mètres, en pieds, en kilomètres, en miles et en unités astronomiques.2 Le temps est une autre grandeur que nous mesurons très aisément, en fait la plupart d’entre nous ont probablement appris à lire l’heure avant d’apprendre à mesurer des distances.3 Les unités utilisées pour mesurer le temps sont la seconde, la minute, le jour, la semaine, la quinzaine,4 et ainsi de suite. Le mois et l’année ne sont souvent pas de bonnes unités à utiliser pour le temps car les différents mois et années ont des durées différentes.

La grandeur masse n’est pas aussi intuitive que la longueur et le temps. La mesure de la masse d’un objet est souvent perçue comme mesurant la « quantité de matière » dans l’objet. Ce n’est pas une mauvaise définition (ou du moins, pas complètement terrible), mais elle n’est pas tout à fait correcte non plus [1]. Une définition plus précise pourrait être que la masse est une mesure de l’inertie, c’est-à-dire de la résistance d’un objet à être accéléré. Plus un objet est massif, plus il faut de force pour le mettre en mouvement, arrêter son mouvement ou modifier son mouvement.

La masse est souvent confondue avec le poids, d’autant plus que les unités utilisées pour mesurer la masse servent aussi à mesurer le poids : le gramme, la livre, le kilogramme, la tonne, etc. La masse d’un objet est une propriété intrinsèque de l’objet, tandis que le poids est un phénomène local qui dépend de l’intensité de l’attraction gravitationnelle exercée par un objet massif voisin. Votre masse sera la même que vous soyez à Chicago, sur la Lune, près de Jupiter ou à des années-lumière du corps céleste le plus proche, mais dans chacun de ces cas votre poids sera très différent. Dans ce livre et dans la plupart des jeux vidéo, nos préoccupations se limitent à un patch relativement petit sur une Terre plate, et nous approximons la gravité par une force constante vers le bas. Il ne sera pas trop dommageable de confondre masse et poids, car la gravité sera pour nous une constante. (Mais nous n’avons pas pu résister à quelques exercices intéressants sur la Station Spatiale Internationale.)

Dans de nombreuses situations, nous pouvons discuter de la relation entre les grandeurs fondamentales sans nous préoccuper des unités de mesure que nous utilisons. Dans de telles situations, nous trouverons utile de désigner la longueur, le temps et la masse par LL, TT et MM respectivement. Un cas important est la définition des grandeurs dérivées. Nous avons dit que la longueur, le temps et la masse sont les grandeurs fondamentales — mais qu’en est-il des autres grandeurs, comme l’aire, le volume, la densité, la vitesse, la fréquence, la force, la pression, l’énergie, la puissance ou n’importe laquelle des nombreuses grandeurs mesurables en physique ? Nous n’attribuons à aucune d’elles sa propre lettre majuscule, car chacune peut être définie en termes des grandeurs fondamentales.

Par exemple, nous pourrions exprimer une mesure de surface en nombre de « pieds carrés ». Nous avons créé une unité en termes d’une autre unité. En physique, nous disons qu’une mesure de surface a l’unité « longueur au carré », soit L2L^{2}. Qu’en est-il de la vitesse ? Nous mesurons la vitesse en miles par heure ou en mètres par seconde. Ainsi la vitesse est le rapport d’une distance par unité de temps, soit L/TL/T.

Un dernier exemple est la fréquence. Vous savez probablement que la fréquence mesure combien de fois quelque chose se produit dans un intervalle de temps donné (à quelle « fréquence » cela arrive). Par exemple, un adulte en bonne santé a un rythme cardiaque moyen d’environ 70 battements par minute (BPM). Le moteur d’une voiture peut tourner à 5 000 tours par minute (RPM). La norme télévisuelle NTSC est définie à 29,97 images par seconde (FPS). Notez que dans chacun de ces cas, nous comptons combien de fois quelque chose se produit dans une durée donnée. Nous pouvons donc écrire la fréquence en unités génériques comme 1/T1/T ou T1T^{- 1}, que vous pouvez lire comme « par unité de temps ». L’une des mesures de fréquence les plus importantes est le Hertz, abrégé Hz, qui signifie « par seconde ». Lorsque vous exprimez une fréquence en Hz, vous décrivez le nombre d’événements, d’oscillations, de battements de cœur, d’images vidéo ou de quoi que ce soit d’autre par seconde. Par définition, 1Hz=1s11\ {Hz} = 1\ s^{- 1}.

Le Tableau 11.1 résume plusieurs grandeurs mesurées en physique, leur relation avec les grandeurs fondamentales, et quelques unités courantes utilisées pour les mesurer.

Grandeur Notation Unité SI Autres unités
Longueur LL mm cmcm , kmkm , inin , ftft , mimi , année-lumière, furlong
Temps TT ss minmin , hrhr , msms
Masse MM kgkg gg , slug, lblb (livre-masse)
Vitesse L/TL/T m/sm/s ft/sft/s , m/hrm/hr , km/hrkm/hr
Accélération L/T2L/T^{2} m/s2{m/s}^{2} ft/s2{ft/s}^{2} , (m/hr)/s(m/hr)/s , (km/hr)/s(km/hr)/s
Force ML/T2ML/T^{2} NN  (Newton) = kgm/s2kg \cdot m/s^{2} lblb (livre-force), poundal
Aire L2L^{2} m2m^{2} mm2mm^{2} , cm2cm^{2} , km2km^{2} , in2in^{2} , ft2ft^{2} , mi2mi^{2} , acre, hectare
Volume L3L^{3} m3m^{3} mm3mm^{3} , cm3cm^{3} , L (litre), in3in^{3} , ft3ft^{3} , cuillère à café, oz liq. (once liquide), tasse, pinte, quart, gallon
Pression Force/Aire = (ML/T2)/L2(ML/T^{2})/L^{2} = M/(LT2)M/(LT^{2}) PaPa  (Pascal) = N/m2N/m^{2} = kg/(ms2)kg/(m \cdot s^{2}) psipsi  ( lbs/in2lbs/in^{2} ), millibar, pouce de mercure, atmatm  (atmosphère)
Énergie Force ×\times Longueur = (ML/T2)L(ML/T^{2}) \cdot L = ML2/T2ML^{2}/T^{2} JJ  (Joule) = Nm2N \cdot m^{2} = kgms2m\frac{kg \cdot m}{s^{2}} \cdot m = kgm2s2\frac{kg \cdot m^{2}}{s^{2}} kWhrkW \cdot hr  (kilowattheure), pied-livre, erg, calorie, BTUBTU  (unité thermique britannique), tonne de TNT
Puissance Énergie / Temps = (ML2/T2)/T(ML^{2}/T^{2})/T = ML2/T3ML^{2}/T^{3} WW  (Watt) = J/sJ/s = kgm2s2s1\frac{kg \cdot m^{2}}{s^{2}} \cdot s^{- 1} = kgm2s3\frac{kg \cdot m^{2}}{s^{3}} hphp  (cheval-vapeur)
Fréquence 1/T=T11/T = T^{- 1} Hz=1/s=s1Hz = 1/s = s^{- 1} = « par seconde » KHz=1000HzKHz = 1 000\ Hz , MHz=1000000HzMHz = 1 000 000\ Hz , « par minute », « par an »

Tableau 11.1Grandeurs physiques sélectionnées et unités de mesure courantes

Bien entendu, toute mesure réelle n’a de sens qu’avec des unités spécifiques attachées. Une façon de s’assurer que vos calculs ont toujours un sens est de transporter les unités à tout moment et de les traiter comme des variables algébriques. Par exemple, si vous calculez une pression et que votre réponse ressort avec les unités m/s, vous savez que vous avez fait une erreur ; la pression a des unités de force par unité de surface, soit ML/(T2L2)ML/(T^{2}L^{2}). En revanche, si vous résolvez un problème et obtenez une réponse en livres par pouce carré (psi), mais que vous cherchez une valeur en Pascals, votre réponse est probablement correcte mais doit juste être convertie dans les unités souhaitées. Ce type de raisonnement est connu sous le nom d’analyse dimensionnelle. Transporter les unités et les traiter comme des variables algébriques met souvent en évidence les erreurs causées par des unités de mesure différentes, et facilite également la conversion d’unités.

Puisque la conversion d’unités est une compétence importante, revoyons-la brièvement. Le concept de base est que pour convertir une mesure d’un ensemble d’unités à un autre, on multiplie cette mesure par une fraction bien choisie qui a la valeur 1. Prenons un exemple simple : combien de pieds font 14,57 mètres ? En cherchant le facteur de conversion,5 nous voyons que 1m3,28083ft1{\ m} \approx 3,28083{\ ft} . Cela signifie que 1m/3,28083ft11{\ m}/3,28083{\ ft} \approx 1 . Prenons donc notre mesure et multiplions-la par une « valeur spéciale » de 1 :

14,57m=14,57m×114,57m×3,28083ft1m47,80ft.\begin{matrix} {14,57{\ m} = 14,57{\ m} \times 1 \approx 14,57{\ m} \times \frac{3,28083{\ ft}}{1{\ m}} \approx 47,80{\ ft}.} \\ \end{matrix}

Notre facteur de conversion nous indique que le numérateur et le dénominateur de la fraction de l’Équation (11.1) sont égaux : 3,28083 pieds égalent 1 mètre. Comme le numérateur et le dénominateur sont égaux, la « valeur » de cette fraction est 1. (Au sens physique, bien sûr, la fraction n’égale pas numériquement 1.) Et nous savons que multiplier quelque chose par 1 ne change pas sa valeur. Parce que nous traitons les unités comme des variables algébriques, le m à gauche s’annule avec le m au bas de la fraction.

Bien sûr, appliquer un seul facteur de conversion simple n’est pas très difficile, mais considérons un exemple plus compliqué. Convertissons 188 km/h en ft/s. Cette fois nous devons multiplier par « 1 » plusieurs fois :

188kmhr×1hr3600s×1000m1km×3,28083ft1m171fts.188\frac{km}{hr} \times \frac{1{\ hr}}{3600{\ s}} \times \frac{1000{\ m}}{1{\ km}} \times \frac{3,28083{\ ft}}{1{\ m}} \approx 171\frac{ft}{s}.

11.3Vitesse moyenne

Nous commençons notre étude de la cinématique en examinant de plus près le concept simple de vitesse. Comment mesure-t-on la vitesse ? La méthode la plus courante est de mesurer le temps qu’il faut pour parcourir une distance fixe. Par exemple, dans une course, nous disons que le coureur le plus rapide est celui qui franchit la ligne d’arrivée dans le temps le plus court.

Considérons la fable de la tortue et du lièvre. Dans l’histoire, ils décident de faire une course, et le lièvre, après avoir pris une avance rapide, devient trop confiant et se laisse distraire. Il s’arrête pendant la course pour faire une sieste, sentir les fleurs, ou quelque autre forme de flânerie. Pendant ce temps, la tortue avance patiemment, finit par dépasser le lièvre et franchit la ligne d’arrivée en premier. Or, c’est un livre de mathématiques et non un livre de développement personnel, aussi veuillez ignorer les leçons morales sur la concentration et la persévérance que l’histoire contient, et considérez plutôt ce qu’elle nous enseigne sur la vitesse moyenne. Examinez la Figure 11.1, qui représente la position de chaque animal au fil du temps.

image

Figure 11.1 Graphe de la position en fonction du temps lors de la course entre la tortue et le lièvre

Voici le déroulement de la course. Le coup de feu part au temps t0t_{0}, et le lièvre sprinte jusqu’au temps t1t_{1}. À ce stade, son orgueil lui fait ralentir, jusqu’au temps t2t_{2} où une jolie lapine passe en sens inverse. (Sa position au fil du temps n’est pas représentée dans le diagramme.) À ce stade, un autre trait masculin tragique pousse le lièvre à faire demi-tour et à marcher avec elle, et il entreprend de la courtiser. À t3t_{3}, il réalise que ses avances ne mènent nulle part, et il commence à arpenter la piste de long en large, abattu, jusqu’au temps t4t_{4}. À ce moment, il décide de faire une sieste. Pendant ce temps, la tortue progresse lentement mais sûrement, et au temps t5t_{5}, elle rattrape le lièvre endormi. La tortue avance et franchit le fil d’arrivée à t6t_{6}. Peu après, le lièvre, peut-être réveillé par le bruit de la foule célébrant la victoire de la tortue, se réveille à t7t_{7} et se précipite frénétiquement vers l’arrivée. À t8t_{8}, le lièvre franchit la ligne d’arrivée, où il est humilié par tous ses pairs, et par la jolie lapine aussi.

Pour mesurer la vitesse moyenne de l’un ou l’autre animal pendant n’importe quel intervalle de temps, on divise le déplacement de l’animal par la durée de l’intervalle. Nous allons nous concentrer sur le lièvre, et nous noterons la position du lièvre xx, ou plus explicitement x(t)x(t), pour souligner le fait que la position du lièvre varie en fonction du temps. Il est courant d’utiliser la lettre grecque majuscule delta (« Δ\Delta ») comme préfixe pour signifier « variation de ». Par exemple, Δx\Delta x signifierait « la variation de la position du lièvre », c’est-à-dire un déplacement du lièvre. De même, Δt\Delta t signifie « la variation du temps actuel », ou simplement « temps écoulé entre deux instants ». En utilisant cette notation, la vitesse moyenne du lièvre de tat_{a} à tbt_{b} est donnée par l’équation

Définition de la vitesse moyenne

vitessemoyenne=déplacementtempsécoulé=ΔxΔt=x(tb)x(ta)tbta.{vitesse\ moyenne} = \frac{déplacement}{temps\ écoulé} = \frac{\Delta x}{\Delta t} = \frac{x(t_{b}) - x(t_{a})}{t_{b} - t_{a}}.

C’est la définition de la vitesse moyenne. Quelles que soient les unités spécifiques utilisées, la vitesse décrit toujours le rapport d’une longueur divisée par un temps, ou pour utiliser la notation discutée dans la Section 11.2, la vitesse est une grandeur avec les unités L/TL/T.

Si nous traçons une droite passant par deux points quelconques sur le graphe de la position du lièvre, la pente de cette droite mesure la vitesse moyenne du lièvre sur l’intervalle entre les deux points. Par exemple, considérons la vitesse moyenne du lièvre au moment où il décélère du temps t1t_{1} à t2t_{2}, comme indiqué dans la Figure 11.2. La pente de la droite est le rapport Δx/Δt\Delta x/\Delta t. Cette pente est également égale à la tangente de l’angle marqué α\alpha, bien que pour l’instant les valeurs Δx\Delta x et Δt\Delta t soient celles que nous aurons sous la main, donc nous n’aurons pas besoin de faire de trigonométrie.

image

Figure 11.2Détermination graphique de la vitesse moyenne

En revenant à la Figure 11.1, remarquez que la vitesse moyenne du lièvre de t2t_{2} à t3t_{3} est négative. En effet, la vitesse est définie comme le rapport du déplacement net sur le temps. Comparez cela à la rapidité, qui est la distance totale divisée par le temps et ne peut pas être négative. Le signe du déplacement et de la vitesse sont sensibles à la direction du déplacement, alors que la distance et la rapidité sont intrinsèquement non négatives. Nous avons déjà évoqué ces distinctions bien plus tôt dans Section 2.2. Bien sûr, il est évident que la vitesse moyenne est négative entre t2t_{2} et t3t_{3}, puisque le lièvre reculait pendant tout l’intervalle. Mais la vitesse moyenne peut aussi être négative sur un intervalle même dans des situations où une progression en avant est réalisée pendant une partie de l’intervalle, comme le plus grand intervalle entre t2t_{2} et t4t_{4}. C’est un cas de « un pas en avant, deux pas en arrière ».

La vitesse moyenne peut aussi être nulle, comme illustré pendant la sieste du lièvre de t4t_{4} à t7t_{7}. En fait, la vitesse moyenne sera nulle chaque fois qu’un objet commence et finit au même endroit, même s’il était en mouvement pendant tout l’intervalle ! (« Deux pas en avant, deux pas en arrière. ») Deux tels intervalles sont illustrés dans la Figure 11.3.

image

Figure 11.3 Deux intervalles pendant lesquels le lièvre n’a aucun déplacement net, et donc sa vitesse moyenne est nulle

Et bien sûr, la leçon finale de la fable est que la vitesse moyenne de la tortue est supérieure à la vitesse moyenne du lièvre, du moins de t0t_{0} à t7t_{7}, quand la tortue franchit la ligne d’arrivée. Cela est vrai malgré le fait que la rapidité moyenne du lièvre était plus élevée, puisqu’il a certainement parcouru une plus grande distance avec toutes les distractions féminines et les allées et venues.

Une dernière chose à souligner. Si nous supposons que le lièvre a tiré les leçons de son expérience et félicité la tortue (après tout, n’attribuons pas au pauvre animal tous les traits de personnalité négatifs !), alors à t=t8t = t_{8} ils se trouvaient au même endroit. Cela signifie que leurs déplacements nets de t0t_{0} à t8t_{8} sont identiques, et donc qu’ils ont la même vitesse moyenne pendant cet intervalle.

11.4Vitesse instantanée et la dérivée

Nous avons vu comment la physique définit et mesure la vitesse moyenne d’un objet sur un intervalle, c’est-à-dire entre deux valeurs de temps qui diffèrent d’une quantité finie Δt\Delta t. Cependant, il est souvent utile de pouvoir parler de la vitesse instantanée d’un objet, c’est-à-dire la vitesse de l’objet pour une valeur de tt, un seul instant dans le temps. On voit que ce n’est pas une question triviale car les méthodes habituelles pour mesurer la vitesse, telles que

vitessemoyenne=déplacementtempsécoulé=ΔxΔt=x(tb)x(ta)tbta,{vitesse\ moyenne} = \frac{déplacement}{temps\ écoulé} = \frac{\Delta x}{\Delta t} = \frac{x(t_{b}) - x(t_{a})}{t_{b} - t_{a}},

ne fonctionnent pas quand on considère un seul instant dans le temps. Que sont tat_{a} et tbt_{b} quand on ne regarde qu’une seule valeur de temps ? En un seul instant, le déplacement et le temps écoulé sont tous deux nuls ; quel est donc le sens du rapport Δx/Δt\Delta x/\Delta t ? Cette section introduit un outil fondamental du calcul connu sous le nom de dérivée. La dérivée a été inventée par Newton pour étudier précisément les questions de cinématique que nous posons dans ce chapitre. Cependant, son applicabilité s’étend à pratiquement tous les problèmes où une grandeur varie en fonction d’une autre grandeur. (Dans le cas de la vitesse, nous nous intéressons à la façon dont la position varie en fonction du temps.)

Étant donné le vaste éventail de problèmes auxquels la dérivée peut être appliquée, Newton n’était pas le seul à l’avoir étudiée. Les applications primitives du calcul intégral pour calculer des volumes remontent à l’Égypte antique. Dès le Ve siècle, les Grecs exploraient les éléments constitutifs du calcul comme les infinitésimaux et la méthode d’exhaustion. Newton partage généralement le mérite avec le mathématicien allemand Gottfried Leibniz6 (1646–1716) pour avoir inventé le calcul au XVIIe siècle, bien que des écrits persans et indiens contiennent des exemples de concepts de calcul utilisés. De nombreux autres penseurs ont apporté des contributions significatives, dont Fermat, Pascal et Descartes.7 Il est assez intéressant de constater que nombre des premières applications du calcul portaient sur les intégrales, même si la plupart des cours de calcul couvrent la dérivée « plus facile » avant l’intégrale « plus difficile ».

Nous suivons d’abord les traces de Newton et commençons par l’exemple physique de la vitesse, que nous considérons comme le meilleur exemple pour acquérir de l’intuition sur le fonctionnement de la dérivée. Ensuite, nous considérons plusieurs autres exemples où la dérivée peut être utilisée, en passant du physique à l’abstrait.

11.4.1Arguments de limite et la définition de la dérivée

Retour à la question posée : comment mesurer la vitesse instantanée ? Observons d’abord une situation particulière pour laquelle c’est facile : si un objet se déplace à vitesse constante sur un intervalle, alors la vitesse est la même à chaque instant de l’intervalle. C’est la définition même de la vitesse constante. Dans ce cas, la vitesse moyenne sur l’intervalle doit être la même que la vitesse instantanée pour tout point de cet intervalle. Dans un graphe comme la Figure 11.1, il est facile de déterminer quand l’objet se déplace à vitesse constante car le graphe est une droite. En fait, presque toute la Figure 11.1 est composée de segments de droite,8 donc déterminer la vitesse instantanée est aussi simple que de choisir deux points quelconques sur un intervalle rectiligne (les extrémités de l’intervalle semblent un bon choix, mais deux points quelconques feront l’affaire) et de déterminer la vitesse moyenne entre ces extrémités.

Mais considérons l’intervalle de t1t_{1} à t2t_{2}, pendant lequel la confiance excessive du lièvre l’amène à décélérer progressivement. Sur cet intervalle, le graphe de la position du lièvre est une courbe, ce qui signifie que la pente de la droite, et donc la vitesse du lièvre, change continuellement. Dans cette situation, mesurer la vitesse instantanée demande un peu plus de finesse.

Pour plus de clarté dans cet exemple, attribuons des nombres particuliers. Pour garder ces nombres ronds (et aussi pour rester dans le thème des courses), permettez le choix fantaisiste de mesurer le temps en minutes et la distance en furlongs.9 Nous posons t1=1mint_{1} = 1\ {min} et t2=3mint_{2} = 3\ {min}, la durée totale est donc de 2 minutes. Disons que pendant cet intervalle, le lièvre parcourt de x(1)=4furx(1) = 4\ {fur} à x(3)=8furx(3) = 8\ {fur} .10 À des fins d’illustration, nous visons la réponse à la question : quelle est la vitesse instantanée du lièvre à t=2,5mint = 2,5\ {min} ? Tout ceci est représenté dans la Figure 11.4.

image

Figure 11.4 Quelle est la vitesse du lièvre à t=2,5mint = 2,5\ {min} ?

Il n’est pas immédiatement évident comment mesurer ou calculer la vitesse au moment exact t=2,5t = 2,5, mais observons que nous pouvons obtenir une bonne approximation en calculant la vitesse moyenne d’un très petit intervalle près de t=2,5t = 2,5. Pour un intervalle suffisamment petit, le graphe ressemble presque à un segment de droite, et la vitesse est presque constante, donc la vitesse instantanée à tout instant donné dans l’intervalle ne sera pas très éloignée de la vitesse moyenne sur tout l’intervalle.

Dans la Figure 11.5, nous fixons l’extrémité gauche d’un segment à t=2,5t = 2,5 et déplaçons l’extrémité droite de plus en plus près. Comme vous pouvez le voir, plus l’intervalle est court, plus le graphe ressemble à une droite, et meilleure est notre approximation. En pensant graphiquement, à mesure que le deuxième point se rapproche de plus en plus de t=2,5t = 2,5, la pente de la droite entre les deux points convergera vers la pente de la droite tangente à la courbe en ce point. Une droite tangente est l’équivalent graphique de la vitesse instantanée, car elle mesure la pente de la courbe exactement en ce seul point.

image

Figure 11.5 Approximation de la vitesse instantanée comme vitesse moyenne d’un intervalle de plus en plus petit

Effectuons cette expérience avec de vrais nombres et voyons si nous pouvons approcher la vitesse instantanée du lièvre. Pour ce faire, nous aurons besoin de connaître la position du lièvre à tout instant donné, il est donc temps de vous dire que la position du lièvre est donnée par la fonction11

x(t)=t2+6t1.x(t) = - t^{2} + 6t - 1.

Le Tableau 11.2 montre des calculs tabulés de vitesse moyenne sur des intervalles avec une extrémité droite t+Δtt + \Delta t qui se rapproche de plus en plus de t=2,5t = 2,5.

La colonne la plus à droite, qui représente la vitesse moyenne, semble converger vers une vitesse de 1 furlong/minute. Mais dans quelle mesure sommes-nous certains que c’est la valeur correcte ? Bien qu’il n’existe aucun calcul qui produira une vitesse résultante exactement égale à 1 furlong/minute, à toutes fins pratiques, nous pouvons atteindre n’importe quel degré de précision souhaité en utilisant cette technique d’approximation et en choisissant Δt\Delta t suffisamment petit. (Nous ignorons les problèmes liés à la précision de la représentation en virgule flottante des nombres dans un ordinateur.)

tt Δt\Delta t t+Δtt + \Delta t x(t)x(t) x(t+Δt)x(t + \Delta t) x(t+Δt)x(t)x(t + \Delta t) - x(t) x(t+Δt)x(t)Δt\frac{x(t + \Delta t) - x(t)}{\Delta t}
[6pt] 2.500 0.500 3.000 7.750 8.0000 0.2500 0.5000
2.500 0.100 2.600 7.750 7.8400 0.0900 0.9000
2.500 0.050 2.550 7.750 7.7975 0.0475 0.9500
2.500 0.010 2.510 7.750 7.7599 0.0099 0.9900
2.500 0.005 2.505 7.750 7.7549 0.0049 0.9950
2.500 0.001 2.501 7.750 7.7509 0.0009 0.9990

Tableau 11.2Calcul de la vitesse moyenne pour des intervalles de durées variées

C’est un argument puissant. Nous avons essentiellement attribué une valeur à une expression que nous ne pouvons pas évaluer directement. Bien qu’il soit mathématiquement illégal de substituer Δt=0\Delta t = 0 dans l’expression, nous pouvons soutenir que pour des valeurs de plus en plus petites de Δt\Delta t, nous convergeons vers une valeur particulière. Dans le jargon du calcul, cette valeur de 1 furlong/minute est une valeur limite, ce qui signifie que lorsque nous prenons des valeurs positives de plus en plus petites pour Δt\Delta t, le résultat de notre calcul approche 1, sans le dépasser (ni l’atteindre exactement).

Les arguments de convergence tels que celui-ci sont définis rigoureusement en calcul en utilisant un outil formalisé connu sous le nom de limite. La notation mathématique pour cela est

v(t)=limΔt0x(t+Δt)x(t)Δt.\begin{matrix} {v(t) = \lim\limits_{\Delta t \rightarrow 0}\frac{x(t + \Delta t) - x(t)}{\Delta t}.} \\ \end{matrix}

La notation « \rightarrow » se lit généralement « approche » ou « tend vers ». Donc le côté droit de l’Équation (11.2) pourrait se lire

Lalimitedex(t+Δt)x(t)ΔtquandΔttendverszéro,{''La\ limite\ de\ }\frac{x(t + \Delta t) - x(t)}{\Delta t}\ {quand}\ \Delta t\ {tend\ vers\ zéro,\text{''}}

ou

LalimitequandΔttendverszérodex(t+Δt)x(t)Δt.{''La\ limite\ quand}\ \Delta t\ {tend\ vers\ zéro\ de\ }\frac{x(t + \Delta t) - x(t)}{\Delta t}.\text{''}

En général, une expression de la forme limak[truc]\lim\limits_{a \rightarrow k}\lbrack{truc}\rbrack est interprétée comme signifiant « La valeur vers laquelle [truc] converge, quand aa se rapproche de plus en plus de kk. »

C’est une idée importante, car elle définit ce que nous entendons par vitesse instantanée.

La vitesse instantanée à un temps tt donné peut être interprétée comme la vitesse moyenne d’un intervalle contenant tt, dans la limite où la durée de l’intervalle tend vers zéro.

Nous n’aurons pas grand besoin d’explorer toute la puissance des limites ni de nous enliser dans les détails ; c’est le domaine mathématique de l’analyse, et nous entraînerait un peu trop loin de nos objectifs actuels, plutôt limités.12 Nous passons sous silence certains détails importants13 pour pouvoir nous concentrer sur un cas particulier : l’utilisation des limites pour définir la dérivée.

La dérivée mesure le taux de variation d’une fonction. Rappelons que « fonction » n’est qu’un terme savant pour désigner toute formule, calcul ou procédure qui prend une entrée et produit une sortie. La dérivée quantifie le taux auquel la sortie de la fonction changera en réponse à une variation de l’entrée. Si xx désigne la valeur d’une fonction à un temps spécifique tt, la dérivée de cette fonction en tt est le rapport dx/dtdx/dt. Le symbole dxdx représente la variation de la sortie produite par une très petite variation de l’entrée, représentée par dtdt. Nous en dirons plus sur ces « petites variations » dans un instant.

Pour l’instant, nous sommes sur un circuit imaginaire où des lapins et des tortues font la course et où des leçons morales sont enseignées par métaphore. Nous avons une fonction avec une entrée tt, le nombre de minutes écoulées depuis le début de la course, et une sortie xx, la distance du lièvre sur le circuit. La règle que nous utilisons pour évaluer notre fonction est l’expression x(t)=t2+6t1x(t) = - t^{2} + 6t - 1. La dérivée de cette fonction nous indique le taux de variation de la position du lièvre par rapport au temps et constitue la définition de la vitesse instantanée. Juste avant, nous avons défini la vitesse instantanée comme la vitesse moyenne prise sur des intervalles de plus en plus petits, mais c’est essentiellement la même chose que la définition de la dérivée. Nous avons simplement formulé les choses la première fois en utilisant une terminologie spécifique à la position et à la vitesse.

Lorsque nous calculons une dérivée, nous n’obtenons pas un seul nombre. Attendre que la réponse à « Quelle est la vitesse du lièvre ? » soit un seul nombre n’a de sens que si la vitesse est la même partout. Dans un tel cas trivial, nous n’avons pas besoin de dérivées, nous pouvons utiliser la vitesse moyenne. La situation intéressante se produit quand la vitesse varie dans le temps. Lorsque nous calculons la dérivée d’une fonction de position dans de tels cas, nous obtenons une fonction de vitesse, ce qui nous permet de calculer la vitesse instantanée à tout moment.

Les trois paragraphes précédents expriment les concepts les plus importants de cette section, aussi permettez-nous de les répéter.

Une dérivée mesure un taux de variation. Puisque la vitesse est le taux de variation de la position par rapport au temps, la dérivée de la fonction de position est la fonction de vitesse.

Les prochaines sections discutent des mathématiques des dérivées de manière un peu plus détaillée, et nous revenons à la cinématique dans la Section 11.5. Ce matériel est destiné à ceux qui n’ont pas14 suivi le calcul de première année. Si vous avez déjà une formation en calcul, vous pouvez passer directement à la Section 11.5, sauf si vous souhaitez un rappel.

La Section 11.4.2 liste plusieurs exemples de dérivées pour vous donner une meilleure compréhension de ce que signifie mesurer un taux de variation, et aussi pour appuyer notre affirmation selon laquelle la dérivée a une applicabilité très large. La Section 11.4.3 donne la définition mathématique formelle de la dérivée15 et montre comment utiliser cette définition pour résoudre des problèmes. Nous allons enfin déterminer la vitesse du lièvre à t=2,5t = 2,5. La Section 11.4.4 liste diverses notations alternatives couramment utilisées pour les dérivées, et enfin, la Section 11.4.5 liste juste assez de règles sur les dérivées pour satisfaire les besoins en calcul différentiel très modestes de ce livre.

11.4.2Exemples de dérivées

La vitesse est peut-être l’introduction la plus facile à la dérivée, mais ce n’est certainement pas le seul exemple. Examinons d’autres exemples pour vous donner une idée du large éventail de problèmes auxquels la dérivée est appliquée.

Les types d’exemples les plus simples consistent à considérer d’autres grandeurs qui varient avec le temps. Par exemple, si R(t)R(t) est la lecture d’un pluviomètre à un temps donné tt, alors la dérivée, notée R(t)R^{\prime}(t), décrit l’intensité de la pluie au temps tt. Peut-être que P(t)P(t) est la lecture d’une soupape de pression sur un réservoir contenant un type de gaz. En supposant que la lecture de pression est proportionnelle à la masse du gaz dans la chambre,16 le taux de variation P(t)P^{\prime}(t) indique la vitesse à laquelle le gaz entre ou sort de la chambre au temps tt.

Il existe aussi des exemples physiques pour lesquels la variable indépendante n’est pas le temps. Le cas prototypique est une fonction y(x)y(x) qui donne la hauteur d’une surface au-dessus d’un point de référence à la position horizontale xx. Par exemple, xx peut être la distance le long de notre circuit métaphorique et yy mesure la hauteur en ce point au-dessus ou en dessous de l’altitude au point de départ. La dérivée y(x)y^{\prime}(x) de cette fonction est la pente de la surface en xx, où les pentes positives signifient que les coureurs montent, et les valeurs négatives indiquent une portion descendante de la course. Cet exemple n’est pas vraiment nouveau, car nous avons examiné des graphes de fonctions et considéré comment la dérivée mesure la pente du graphe en 2D.

Devenons maintenant un peu plus abstraits, tout en gardant une dimension physique comme variable indépendante. Disons que pour un mur d’escalade populaire, nous connaissons une fonction S(y)S(y) qui décrit, pour une hauteur donnée yy, quel pourcentage des grimpeurs est capable d’atteindre cette hauteur ou davantage. Si nous supposons que les grimpeurs commencent à y=0y = 0, alors S(0)=100%S(0) = 100\%. Clairement S(y)S(y) est une fonction non croissante qui finit par descendre jusqu’à 0% à une hauteur maximale ymaxy_{max} que personne n’a jamais atteinte.

Considérons maintenant l’interprétation de la dérivée S(y)S^{\prime}(y). Bien sûr, S(y)0S^{\prime}(y) \leq 0, puisque S(y)S(y) est non croissante. Une grande valeur négative de S(y)S^{\prime}(y) indique que la hauteur yy est une zone où les grimpeurs ont tendance à lâcher17. Peut-être que le mur à cette hauteur est une zone difficile. S(y)S^{\prime}(y) proche de zéro indique que moins de grimpeurs abandonnent à la hauteur yy. Peut-être y a-t-il un plateau que les grimpeurs peuvent atteindre, et là ils se reposent. On pourrait s’attendre à ce que S(y)S^{\prime}(y) diminue juste après ce plateau, car les grimpeurs sont plus reposés. En fait, S(y)S^{\prime}(y) pourrait aussi se rapprocher de zéro juste avant le plateau, car à mesure que les grimpeurs commencent à approcher de ce jalon, ils poussent un peu plus fort et sont plus réticents à abandonner.18

image

Figure 11.6Bonheur en fonction du salaire

Un dernier exemple. La Figure 11.6 montre le bonheur en fonction du salaire. Dans ce cas, la dérivée est essentiellement la même chose que ce que les économistes appelleraient l’« utilité marginale ». C’est le rapport d’unités supplémentaires de bonheur par unité supplémentaire de revenu. Selon cette figure, l’utilité marginale du revenu diminue, ce qui est bien sûr la célèbre loi des rendements décroissants. Selon nos recherches,19 elle devient même négative après un certain point, où les problèmes associés à un revenu élevé commencent à l’emporter sur les avantages psychologiques. L’expression économique « utilité marginale négative » se traduit en langage courant par « arrêtez de faire ça ».

11.4.3Calcul des dérivées à partir de la définition

Nous sommes maintenant prêts pour la définition20 officielle de la dérivée que l’on trouve dans la plupart des manuels de mathématiques, et pour voir comment nous pouvons calculer des dérivées en utilisant la définition. Une dérivée peut être comprise comme la valeur limite de Δx/Δt\Delta x/\Delta t, le rapport de la variation de la sortie divisée par la variation de l’entrée, prise quand nous rendons Δt\Delta t infinitésimalement petit. Répétons cette description en notation mathématique. C’est une équation que nous avons donnée plus tôt dans le chapitre, mais cette fois nous la mettons dans un grand cadre, car c’est ce que font les manuels de mathématiques avec les équations qui sont des définitions.

La Définition d’une Dérivée

dxdt=limΔt0ΔxΔt=limΔt0x(t+Δt)x(t)Δt.\begin{matrix} {\frac{dx}{dt} = \lim\limits_{\Delta t \rightarrow 0}\frac{\Delta x}{\Delta t} = \lim\limits_{\Delta t \rightarrow 0}\frac{x(t + \Delta t) - x(t)}{\Delta t}.} \\ \end{matrix}

Ici, la notation pour la dérivée dx/dtdx/dt est connue sous le nom de notation de Leibniz. Les symboles dxdx et dtdt sont connus sous le nom d’infinitésimaux. Contrairement à Δx\Delta x et Δt\Delta t, qui sont des variables représentant des variations finies de valeur, dxdx et dtdt sont des symboles représentant « une variation infinitésimalement petite ». Pourquoi est-il si important d’utiliser une très petite variation ? Pourquoi ne peut-on pas simplement prendre le rapport Δx/Δt\Delta x/\Delta t directement ? Parce que le taux de variation varie continuellement. Même sur un très petit intervalle de Δt=0,0001\Delta t = 0,0001, il n’est pas constant. C’est pourquoi un argument de limite est utilisé, pour rendre l’intervalle aussi petit que possible — infinitésimalement petit.

Dans certaines circonstances, les infinitésimaux peuvent être manipulés comme des variables algébriques (et vous pouvez aussi y attacher des unités de mesure et effectuer une analyse dimensionnelle pour vérifier votre travail). Le fait que de telles manipulations soient souvent correctes est ce qui donne à la notation de Leibniz son attrait intuitif. Cependant, parce qu’ils sont des valeurs infiniment petites, ils nécessitent une manipulation spéciale, similaire au symbole \infty, et ne doivent donc pas être utilisés à la légère. La plupart du temps, nous interprétons la notation dxdt\frac{dx}{dt} non pas comme un rapport de deux variables, mais comme un seul symbole signifiant « la dérivée de xx par rapport à tt ». C’est la procédure la plus sûre et elle évite tout risque d’utilisation intempestive précédemment mentionnée. Nous aurons plus à dire sur Leibniz et d’autres notations plus loin, mais d’abord, calculons enfin une dérivée et répondons à la question brûlante : à quelle vitesse le lièvre voyageait-il à t=2,5t = 2,5 ?

Différencier une fonction simple en utilisant la définition Équation (11.3) est un rite de passage important, et nous sommes fiers de vous aider à franchir ce seuil. La procédure typique est la suivante :

  1. Substituer x(t)x(t) et x(t+Δt)x(t + \Delta t) dans la définition. (Dans notre cas, x(t)=t2+6t1x(t) = - t^{2} + 6t - 1 ).

  2. Effectuer des manipulations algébriques jusqu’à ce qu’il soit légal de substituer Δt=0\Delta t = 0. (Cela revient souvent à éliminer Δt\Delta t du dénominateur.)

  3. Substituer Δt=0\Delta t = 0, ce qui évalue l’expression « à la limite », supprimant la notation de limite.

  4. Simplifier le résultat.

L’application de cette procédure à notre cas donne

v(t)=dxdt=limΔt0x(t+Δt)x(t)Δt=limΔt0[(t+Δt)2+6(t+Δt)1](t2+6t1)Δt=limΔt0(t22t(Δt)(Δt)2+6t+6(Δt)1)+(t26t+1)Δt=limΔt02t(Δt)(Δt)2+6(Δt)Δt=limΔt0Δt(2tΔt+6)Δt=limΔt02tΔt+6.\begin{matrix} {v(t) = \frac{dx}{dt}} & {= \lim\limits_{\Delta t \rightarrow 0}\frac{x(t + \Delta t) - x(t)}{\Delta t}} \\ & {= \lim\limits_{\Delta t \rightarrow 0}\frac{\lbrack - (t + \Delta t)^{2} + 6(t + \Delta t) - 1\rbrack - ( - t^{2} + 6t - 1)}{\Delta t}} \\ & {= \lim\limits_{\Delta t \rightarrow 0}\frac{( - t^{2} - 2t(\Delta t) - (\Delta t)^{2} + 6t + 6(\Delta t) - 1) + (t^{2} - 6t + 1)}{\Delta t}} \\ & {= \lim\limits_{\Delta t \rightarrow 0}\frac{- 2t(\Delta t) - (\Delta t)^{2} + 6(\Delta t)}{\Delta t}} \\ & {= \lim\limits_{\Delta t \rightarrow 0}\frac{\Delta t\left( - 2t - \Delta t + 6 \right)}{\Delta t}} \\ & {= \lim\limits_{\Delta t \rightarrow 0}{- 2t - \Delta t + 6}.} \\ \end{matrix}

Nous sommes maintenant à l’étape 3. Prendre la limite de l’ Équation (11.4) est maintenant facile ; nous substituons simplement Δt=0\Delta t = 0. Cette substitution n’était pas légale auparavant car il y avait un Δt\Delta t au dénominateur :

v(t)=dxdt=limΔt02tΔt+6=2t(0)+6=2t+6.\begin{matrix} {v(t) = \frac{dx}{dt}} & {= \lim\limits_{\Delta t \rightarrow 0}{- 2t - \Delta t + 6}} \\ & {= - 2t - (0) + 6} \\ & {= - 2t + 6.} \\ \end{matrix}

image image image
v(2.0)=2(2.0)+6=2.0v(2.0) = - 2(2.0) + 6 = 2.0 v(2.4)=2(2.4)+6=1.2v(2.4) = - 2(2.4) + 6 = 1.2 v(2.5)=2(2.5)+6=1.0v(2.5) = - 2(2.5) + 6 = 1.0
image image
v(2.6)=2(2.6)+6=0.8v(2.6) = - 2(2.6) + 6 = 0.8 v(3.0)=2(3.0)+6=0.0v(3.0) = - 2(3.0) + 6 = 0.0

Figure 11.7La vitesse du lièvre et la droite tangente correspondante à des instants sélectionnés

Enfin ! L’Équation (11.5) est la fonction de vitesse que nous cherchions. Elle nous permet d’insérer n’importe quelle valeur de tt et de calculer la vitesse instantanée du lièvre à cet instant. En substituant t=2,5t = 2,5, nous obtenons la réponse à notre question :

v(t)=2t+6,v(2.5)=2(2.5)+6=1.\begin{matrix} {v(t)} & {= - 2t + 6,} \\ {v(2.5)} & {= - 2(2.5) + 6 = 1.} \\ \end{matrix}

Ainsi, la vitesse instantanée du lièvre à t=2,5t = 2,5 était précisément 1 furlong par minute, exactement comme nos arguments précédents l’avaient prédit. Mais maintenant nous pouvons l’affirmer avec confiance.

La Figure 11.7 montre ce point et plusieurs autres le long de l’intervalle étudié. Pour chaque point, nous avons calculé la vitesse instantanée en ce point selon l’Équation (11.5) et avons tracé la droite tangente avec la même pente.

Il est très instructif de comparer côte à côte les graphes de position et de vitesse. La Figure 11.8 compare la position et la vitesse de nos coureurs légendaires.

image

Figure 11.8 Comparaison de la position et de la vitesse

Plusieurs observations intéressantes peuvent être faites à propos de la Figure 11.8.

À ce stade, il convient de reconnaître quelques façons dont notre explication de la dérivée diffère de la plupart des manuels de calcul. Notre approche a été de se concentrer sur un exemple spécifique, celui de la vitesse instantanée. Cela a conduit à quelques différences cosmétiques, comme la notation. Mais il y avait aussi de nombreux points subtils que nous passons sous silence. Par exemple, nous n’avons pas pris la peine de définir les fonctions continues, ni de donner des définitions rigoureuses pour quand la dérivée est définie et quand elle ne l’est pas. Nous avons discuté de l’idée derrière ce qu’est une limite, mais n’avons pas fourni de définition formelle ni considéré les limites approchées par la gauche et par la droite, et les critères d’existence d’une limite bien définie. Nous pensons que commencer par le meilleur exemple intuitif est toujours la façon optimale d’enseigner quelque chose, même si cela signifie « mentir » au lecteur un court moment. Si nous écrivions un manuel de calcul, nous reviendrions à ce stade pour corriger certains de nos mensonges, en passant en revue les points plus subtils et en donnant des définitions plus précises.

Cependant, comme ceci n’est pas un manuel de calcul, nous vous avertirons seulement que ce qui précède est la vue d’ensemble, mais n’est pas suffisant pour gérer de nombreux cas limites lorsque les fonctions font des choses étranges comme partir à l’infini ou présenter des « sauts » ou des « lacunes ». Heureusement, de tels cas limites ne surviennent pas trop souvent pour des fonctions qui modélisent des phénomènes physiques, et ces détails ne deviendront pas un problème pour nous dans le contexte de la physique.

Nous avons bien la place, cependant, de mentionner des notations alternatives pour la dérivée que vous êtes susceptible de rencontrer.

11.4.4Notations pour la dérivée

Plusieurs notations différentes pour les dérivées sont couramment utilisées. Indiquons quelques façons dont d’autres textes pourraient se présenter différemment de ce que nous avons présenté ici. Tout d’abord, il y a une question triviale de dénomination. La plupart des manuels de calcul définissent la dérivée en termes très généraux, où la variable de sortie est nommée yy, le symbole xx désigne la variable d’entrée plutôt que la variable de sortie, et la fonction est simplement nommée ff. En d’autres termes, la fonction à différencier est y=f(x)y = f(x). De plus, beaucoup assignent la « quantité de pas » décroissante à la variable hh plutôt que d’utiliser la notation Δ\Delta, ce qui présente des avantages lors de la résolution des équations résultant du calcul des dérivées à partir de la définition.22 Avec ces variables, ils définiraient la dérivée comme

Définition d’une dérivée en utilisant les variables de la plupart des manuels de calcul

dydx=limh0y(x+h)y(x)h.\begin{matrix} {\frac{dy}{dx} = \lim\limits_{h \rightarrow 0}\frac{y(x + h) - y(x)}{h}.} \\ \end{matrix}

Les différences entre les Équations (11.3) et (11.6) sont clairement cosmétiques.

Une variante de la notation de Leibniz que nous préférons dans ce livre est de préfixer une expression par d/dtd/dt pour signifier « la dérivée par rapport à tt de cette chose à droite ». Par exemple

ddt(t2+5t)\frac{d}{dt}\,(t^{2} + 5t)

peut se lire comme « la dérivée par rapport à tt de t2+5tt^{2} + 5t ». C’est une notation très descriptive et intuitive. Si nous appelons l’expression à droite xx, et interprétons la juxtaposition de symboles comme une multiplication, nous pouvons ramener le xx en haut de la fraction pour retrouver notre notation originale, comme dans

ddt(t2+5t)=ddtx=dxdt.\frac{d}{dt}\,(t^{2} + 5t) = \frac{d}{dt}\, x = \frac{dx}{dt}.

Il est important d’interpréter ces manipulations comme des manipulations notionnelles plutôt que comme ayant une vraie signification mathématique. La notation est attrayante parce que de telles manipulations algébriques avec les infinitésimaux fonctionnent souvent. Mais nous réitérons notre mise en garde contre l’attribution d’une trop grande signification mathématique à de telles opérations.

Une autre notation courante consiste à désigner la dérivée d’une fonction f(x)f(x) avec un prime : f(x)f^{\prime}(x). C’est ce qu’on appelle la notation prime ou notation de Lagrange. Elle est utilisée quand la variable indépendante par rapport à laquelle on différencie est implicite ou comprise par le contexte. En utilisant cette notation, nous définirions la vitesse comme la dérivée de la fonction de position par v(t)=x(t)v(t) = x^{\prime}(t).

Une dernière notation, inventée par Newton et principalement utilisée quand la variable indépendante est le temps (comme dans les équations physiques inventées par Newton), est la notation point. Une dérivée est indiquée en plaçant un point au-dessus de la variable ; par exemple, v(t)=x˙(t)v(t) = \overset{˙}{x}(t).

Voici un résumé des différentes notations pour la dérivée que vous rencontrerez, en utilisant la vitesse et la position comme exemple :

v(t)=dxdt=ddtx(t)=x(t)=x˙(t).v(t) = \frac{dx}{dt} = \frac{d}{dt}x(t) = x^{\prime}(t) = \overset{˙}{x}(t).

11.4.5Quelques règles et raccourcis de différenciation

Revenons maintenant au calcul des dérivées. En pratique, il est rarement nécessaire de revenir à la définition de la dérivée pour différencier une expression. Au lieu de cela, il existe des règles simplificatrices qui permettent de décomposer des fonctions complexes en morceaux plus petits qui peuvent ensuite être différenciés. Il existe aussi des fonctions spéciales, comme lnx\ln x et tanx\tan x, pour lesquelles le travail difficile d’appliquer la définition a déjà été fait et consigné dans ces tableaux qui garnissent l’intérieur des couvertures des manuels de calcul. Pour différencier des expressions contenant de telles fonctions, on se réfère simplement au tableau (bien que nous allions faire un peu de ce « travail difficile » nous-mêmes pour le sinus et le cosinus).

Dans ce livre, nos préoccupations se limitent aux dérivées d’un très petit ensemble de fonctions, qui peuvent heureusement être différenciées avec quelques règles simples. Malheureusement, nous n’avons pas la place ici pour développer les justifications mathématiques derrière ces règles, aussi allons-nous simplement accompagner chaque règle d’une brève explication de son utilisation, et d’un argument intuitif (mathématiquement non rigoureux) pour vous aider à vous convaincre qu’elle fonctionne.

Notre première règle, connue sous le nom de règle des constantes, stipule que la dérivée d’une fonction constante est zéro. Une fonction constante est une fonction qui produit toujours la même valeur. Par exemple, x(t)=3x(t) = 3 est une fonction constante. Vous pouvez substituer n’importe quelle valeur de tt, et cette fonction produit la valeur 3. Puisque la dérivée mesure la rapidité avec laquelle la sortie d’une fonction change en réponse aux changements de l’entrée tt, dans le cas d’une fonction constante, la sortie ne change jamais, et donc la dérivée est x(t)=0x^{\prime}(t) = 0.

La Règle des Constantes

ddtk=0,kest toute constante.\frac{d}{dt}\, k = 0,\quad k\ \text{est\ toute\ constante.}

La règle suivante, parfois appelée règle de la somme, dit que la différenciation est un opérateur linéaire. La signification de « linéaire » est essentiellement identique à notre définition donnée dans le Chapitre 5, mais revoyons-la dans le contexte de la dérivée. Dire que la dérivée est un opérateur linéaire signifie deux choses. Premièrement, pour prendre la dérivée d’une somme, on peut prendre la dérivée de chaque morceau individuellement, et additionner les résultats. C’est intuitif — le taux de variation d’une somme est le taux de variation total de toutes les parties additionnées. Par exemple, considérons un homme qui se déplace sur un train. Sa position dans l’espace monde peut être décrite comme la somme de la position du train, plus la position de l’homme dans l’espace corps du train.23 De même, sa vitesse par rapport au sol est la somme de la vitesse du train par rapport au sol, plus sa vitesse par rapport au train.

Dérivée d’une Somme

ddt[f(t)+g(t)]=ddtf(t)+ddtg(t).\begin{matrix} {\frac{d}{dt}\,\left\lbrack f(t) + g(t) \right\rbrack = \frac{d}{dt}\, f(t) + \frac{d}{dt}\, g(t).} \\ \end{matrix}

La deuxième propriété de la linéarité est que si nous multiplions une fonction par une constante, la dérivée de cette fonction est multipliée par cette même constante. Une façon simple de voir que cela doit être vrai est de considérer les conversions d’unités. Revenons à notre fonction favorite qui donne le déplacement d’un lièvre en fonction du temps, mesuré en furlongs. En prenant la dérivée de cette fonction par rapport au temps, on obtient une vitesse en furlongs par minute. Si quelqu’un ne veut pas utiliser les furlongs, on peut passer des furlongs aux mètres en multipliant la fonction de position originale par un facteur de 201,168. Cela doit multiplier la dérivée par le même facteur, sinon le lièvre changerait soudainement de vitesse juste parce que nous avons changé d’unité.

Dérivée d’une Fonction Multipliée par une Constante

ddt[kf(t)]=k[ddtf(t)],kest toute constante.\begin{matrix} {\frac{d}{dt}\,\left\lbrack kf(t) \right\rbrack = k\left\lbrack \frac{d}{dt}\, f(t) \right\rbrack,\quad k\ \text{est\ toute\ constante}.} \\ \end{matrix}

Si nous combinons les Équations (11.7) et (11.8), nous pouvons énoncer la règle de linéarité de manière plus générale.

La Règle de la Somme

ddt[af(t)+bg(t)]=a[ddtf(t)]+b[ddtg(t)].\begin{matrix} {\frac{d}{dt}\,\left\lbrack af(t) + bg(t) \right\rbrack = a\left\lbrack \frac{d}{dt}\, f(t) \right\rbrack + b\left\lbrack \frac{d}{dt}\, g(t) \right\rbrack.} \\ \end{matrix}

La propriété linéaire de la dérivée est très importante car elle nous permet de décomposer de nombreuses fonctions courantes en morceaux plus petits et plus faciles.

L’une des fonctions les plus importantes et les plus courantes qui doivent être différenciées est aussi la plus facile : le polynôme. En utilisant la propriété linéaire de la dérivée, nous pouvons décomposer, par exemple, un polynôme de degré 4 avec facilité :

x(t)=c4t4+c3t3+c2t2+c1t+c0,dxdt=ddt[c4t4+c3t3+c2t2+c1t+c0]=c4[ddtt4]+c3[ddtt3]+c2[ddtt2]+c1[ddtt]+[ddtc0].\begin{matrix} {x(t)} & {= c_{4}t^{4} + c_{3}t^{3} + c_{2}t^{2} + c_{1}t + c_{0},} \\ \frac{dx}{dt} & {= \frac{d}{dt}\lbrack c_{4}t^{4} + c_{3}t^{3} + c_{2}t^{2} + c_{1}t + c_{0}\rbrack} \\ & {= c_{4}\left\lbrack \frac{d}{dt}\, t^{4} \right\rbrack + c_{3}\left\lbrack \frac{d}{dt}\, t^{3} \right\rbrack + c_{2}\left\lbrack \frac{d}{dt}\, t^{2} \right\rbrack + c_{1}\left\lbrack \frac{d}{dt}\, t \right\rbrack + \left\lbrack \frac{d}{dt}\, c_{0} \right\rbrack.} \\ \end{matrix}

La dernière dérivée ddtc0\frac{d}{dt}c_{0} est nulle par la règle des constantes, puisque c0c_{0} ne varie pas. Il nous reste quatre dérivées simples, dont chacune peut être insérée dans la définition de la dérivée, Équation (11.3), sans trop de difficulté. Résoudre chacune de ces quatre individuellement est considérablement plus facile qu’insérer le polynôme original dans l’Équation (11.3). Si vous faites cet exercice (comme tout étudiant de première année en calcul), vous remarquez deux choses. Premièrement, la fastidiosité algébrique augmente à mesure que la puissance de tt s’élève. Deuxièmement, un motif assez évident se révèle, connu sous le nom de règle de la puissance.

La Règle de la Puissance

ddttn=ntn1,nest un entier.\frac{d}{dt}\, t^{n} = nt^{n - 1},\quad n\ \text{est\ un\ entier}.

Cette règle nous donne les réponses aux quatre dérivées nécessaires ci-dessus :

ddtt4=4t3,ddtt3=3t2,ddtt2=2t1=2t,ddtt=1t0=1.\begin{matrix} {\frac{d}{dt}\, t^{4}} & {= 4t^{3},} & {\frac{d}{dt}\, t^{3}} & {= 3t^{2},} \\ {\frac{d}{dt}\, t^{2}} & {= 2t^{1} = 2t,} & {\frac{d}{dt}\, t} & {= 1t^{0} = 1.} \\ \end{matrix}

Notez que dans la dernière équation, nous avons utilisé l’identité t0=1t^{0} = 1. Cependant, même sans cette identité,24 il devrait être très clair que ddtt\frac{d}{dt}t doit être égal à l’unité. Rappelons que la dérivée répond à la question : « Quel est le taux de variation de la sortie, par rapport au taux de variation de l’entrée ? » Dans le cas de ddtt\frac{d}{dt}t, la « sortie » et l’« entrée » sont toutes deux la variable tt, et donc leurs taux de variation sont égaux. Ainsi, le rapport qui définit la dérivée est égal à un.

Un dernier commentaire avant d’insérer ces résultats dans l’ Équation (11.9) pour différencier notre polynôme. En utilisant l’identité t0=1t^{0} = 1, la règle de la puissance est mise en harmonie avec la règle des constantes :

Dérivée d’une constante, en utilisant la règle de la puissance

ddtk=ddt(kt0)Enutilisantt0=1,=k[ddtt0]Propriétélinéairedeladérivée,=k[0(t1)]Règledelapuissancepourn=0,=0.\begin{matrix} {\frac{d}{dt}\, k} & {= \frac{d}{dt}\,(kt^{0})} & & {{En\ utilisant\ }t^{0} = 1,} \\ & {= k\left\lbrack \frac{d}{dt}\, t^{0} \right\rbrack} & & {{Propriété\ linéaire\ de\ la\ dérivée},} \\ & {= k\lbrack 0(t^{- 1})\rbrack} & & {{Règle\ de\ la\ puissance\ pour\ }n = 0,} \\ & {= 0.} & & \\ \end{matrix}

Revenons à notre polynôme de degré 4. Avec la règle de la somme et la règle de la puissance à notre disposition, nous pouvons l’expédier rapidement :

x(t)=c4t4+c3t3+c2t2+c1t+c0,dxdt=4c4t3+3c3t2+2c2t+c1.\begin{matrix} {x(t)} & {= c_{4}t^{4} + c_{3}t^{3} + c_{2}t^{2} + c_{1}t + c_{0},} \\ \frac{dx}{dt} & {= 4c_{4}t^{3} + 3c_{3}t^{2} + 2c_{2}t + c_{1}.} \\ \end{matrix}

Voici plusieurs autres exemples de l’utilisation de la règle de la puissance. Notez que la règle de la puissance fonctionne aussi pour les exposants négatifs :

ddt(3t54t)=15t44,ddt(t100100+π)=t99,ddt(1t+4t3)=ddt(t1+4t3)=t212t4=1t212t4.\begin{matrix} {\frac{d}{dt}\,(3t^{5} - 4t)} & {= 15t^{4} - 4,} \\ {\frac{d}{dt}\,\left( \frac{t^{100}}{100} + \sqrt{\pi} \right)} & {= t^{99},} \\ {\frac{d}{dt}\,\left( \frac{1}{t} + \frac{4}{t^{3}} \right) = \frac{d}{dt}\,\left( t^{- 1} + 4t^{- 3} \right)} & {= - t^{- 2} - 12t^{- 4} = \frac{- 1}{t^{2}} - \frac{12}{t^{4}}.} \\ \end{matrix}

11.4.6Dérivées de quelques fonctions spéciales avec les séries de Taylor

Cette section examine quelques exemples très spéciaux de différenciation de polynômes. Étant donné une fonction arbitraire f(x)f(x), la série de Taylor de ff est une façon d’exprimer ff comme un polynôme. Chaque terme successif du polynôme est déterminé en prenant une dérivée d’ordre supérieur de la fonction, ce qui est peut-être le point principal des séries de Taylor que vous devriez apprendre quand vous suivrez un vrai cours de calcul, mais pour l’instant nous ne nous intéressons pas à l’origine des séries de Taylor, simplement au fait qu’elles existent. La série de Taylor est un outil très utile dans les jeux vidéo car elle fournit des approximations polynomiales, « faciles » à évaluer par un ordinateur, pour des fonctions qui sont autrement « difficiles » à évaluer. Nous n’avons pas la place pour discuter de beaucoup de choses sur les séries de Taylor en général, mais nous aimerions examiner quelques exemples importants de séries de Taylor. Les séries de Taylor pour les fonctions sinus et cosinus sont

Séries de Taylor pour sin(𝐱)\mathbf{\sin}(\mathbf{x}) et cos(𝐱)\mathbf{\cos}(\mathbf{x})

sinx=xx33!+x55!x77!+x99!+,cosx=1x22!+x44!x66!+x88!+.\begin{matrix} {\sin x} & {= x - \frac{x^{3}}{3!} + \frac{x^{5}}{5!} - \frac{x^{7}}{7!} + \frac{x^{9}}{9!} + \cdots,} \\ {\cos x} & {= 1 - \frac{x^{2}}{2!} + \frac{x^{4}}{4!} - \frac{x^{6}}{6!} + \frac{x^{8}}{8!} + \cdots.} \\ \end{matrix}

Ce motif se poursuit indéfiniment ; en d’autres termes, calculer la valeur exacte de sinx\sin x nécessiterait d’évaluer un nombre infini de termes. Cependant, notez que les dénominateurs des termes croissent très rapidement, ce qui signifie que nous pouvons approximer sinx\sin x simplement en nous arrêtant après un certain nombre de termes et en ignorant le reste.

C’est exactement le processus par lequel les fonctions trigonométriques sont calculées dans un ordinateur. Premièrement, des identités trigonométriques sont utilisées pour ramener l’argument dans une plage restreinte (puisque les fonctions sont périodiques). Cela est fait parce que lorsque la série de Taylor est tronquée, sa précision est la plus grande près d’une valeur particulière de xx, et dans le cas des fonctions trigo, ce point est généralement choisi comme étant x=0x = 0 .25 Ensuite, le polynôme de série de Taylor avec, disons, quatre termes est évalué. Cette approximation est très précise. S’arrêter au terme x7x^{7} est suffisant pour calculer sinx\sin x à environ cinq virgules et demie pour 1<x<+1- 1 < x < + 1.

Toutes ces anecdotes sur les approximations sont intéressantes, mais notre vraie raison d’introduire les séries de Taylor est de les utiliser comme exemples non triviaux de différenciation de polynômes avec la règle de la puissance, et aussi pour apprendre des faits intéressants sur les fonctions sinus, cosinus et exponentielle. Utilisons la règle de la puissance pour différencier le développement en série de Taylor de sin(x)\sin(x). Ce n’est pas si compliqué — nous devons juste différencier chaque terme séparément. Nous ne sommes même pas intimidés par le fait qu’il y a un nombre infini de termes :

Différenciation de la série de Taylor de sin(𝐱)\mathbf{\sin}(\mathbf{x})

ddxsinx=ddx(xx33!+x55!x77!+x99!+)=ddxxddxx33!+ddxx55!ddxx77!+ddxx99!+(Règle de la somme)=13x23!+5x44!7x67!+9x89!+(Règle de la puissance)=1x22!+x44!x66!+x88!+\begin{matrix} {\frac{d}{dx}\,\sin x} & {= \frac{d}{dx}\left( x - \frac{x^{3}}{3!} + \frac{x^{5}}{5!} - \frac{x^{7}}{7!} + \frac{x^{9}}{9!} + \cdots \right)} & & \\ & {= \frac{d}{dx}\, x - \frac{d}{dx}\frac{x^{3}}{3!} + \frac{d}{dx}\frac{x^{5}}{5!} - \frac{d}{dx}\frac{x^{7}}{7!} + \frac{d}{dx}\frac{x^{9}}{9!} + \cdots} & & \text{(Règle\ de\ la\ somme)} \\ & {= 1 - \frac{3x^{2}}{3!} + \frac{5x^{4}}{4!} - \frac{7x^{6}}{7!} + \frac{9x^{8}}{9!} + \cdots} & & \text{(Règle\ de\ la\ puissance)} \\ & {= 1 - \frac{x^{2}}{2!} + \frac{x^{4}}{4!} - \frac{x^{6}}{6!} + \frac{x^{8}}{8!} + \cdots} & & \\ \end{matrix}

Dans la dérivation ci-dessus, nous avons d’abord utilisé la règle de la somme, qui dit que pour différencier tout le polynôme de Taylor, nous pouvons différencier chaque terme individuellement. Ensuite, nous avons appliqué la règle de la puissance à chaque terme, en multipliant à chaque fois par l’exposant et en le décrémentant de un. (Et en se rappelant aussi que ddxx=1\frac{d}{dx}\ x = 1 pour le premier terme.) Pour comprendre la dernière étape, rappelons la définition de l’opérateur factoriel : n!=1×2×3××nn! = 1 \times 2 \times 3 \times \cdots \times n. Ainsi, la constante au numérateur de chaque terme annule le facteur le plus élevé de la factorielle au dénominateur.

L’Équation (11.11) vous semble-t-elle familière ? Elle devrait, car c’est la même que l’Équation (11.10), la série de Taylor pour cosx\cos x. En d’autres termes, nous connaissons maintenant la dérivée de sinx\sin x, et par un processus similaire nous pouvons aussi obtenir la dérivée de cosx\cos x. Énonçons ces faits formellement.26

Dérivées du Sinus et du Cosinus

ddxsinx=cosx,ddxcosx=sinx.\begin{matrix} {\frac{d}{dx}\,\sin x = \cos x,} & & {\frac{d}{dx}\,\cos x = - \sin x.} \\ \end{matrix}

Les dérivées des fonctions sinus et cosinus deviendront utiles dans les sections suivantes.

Examinons maintenant une autre fonction spéciale importante qui jouera un rôle important plus loin dans ce livre, qu’il sera pratique de pouvoir différencier, et qui a aussi une belle série de Taylor bien ordonnée. La fonction dont nous parlons est la fonction exponentielle, notée exe^{x}. La constante mathématique e2,718282e \approx 2,718282 possède de nombreuses propriétés bien connues et intéressantes, et apparaît dans toutes sortes de problèmes, de la finance au traitement du signal. Une grande partie du statut spécial de ee est liée à la nature unique de la fonction exe^{x}. Une manifestation de cette nature unique est que exe^{x} a une série de Taylor si belle :

Série de Taylor de 𝐞𝐱\mathbf{e}^{\mathbf{x}}

ex=1+x+x22!+x33!+x44!+x55!+\begin{matrix} {e^{x} = 1 + x + \frac{x^{2}}{2!} + \frac{x^{3}}{3!} + \frac{x^{4}}{4!} + \frac{x^{5}}{5!} + \cdots} \\ \end{matrix}

En prenant la dérivée, on obtient

ddxex=ddx(1+x+x22!+x33!+x44!+x55!+)=0+1+x1!+x22!+x33!+x44!+=1+x+x22!+x33!+x44!+\begin{matrix} {\frac{d}{dx}\, e^{x}} & {= \frac{d}{dx}\ \left( 1 + x + \frac{x^{2}}{2!} + \frac{x^{3}}{3!} + \frac{x^{4}}{4!} + \frac{x^{5}}{5!} + \cdots \right)} \\ & {= 0 + 1 + \frac{x}{1!} + \frac{x^{2}}{2!} + \frac{x^{3}}{3!} + \frac{x^{4}}{4!} + \cdots} \\ & {= 1 + x + \frac{x^{2}}{2!} + \frac{x^{3}}{3!} + \frac{x^{4}}{4!} + \cdots} \\ \end{matrix}

Mais ce résultat est équivalent à la définition de exe^{x} de l’ Équation (11.12) ; la seule différence entre eux est la question cosmétique du moment où arrêter d’énumérer les termes explicitement et de terminer par « \cdots ». En d’autres termes, la fonction exponentielle est sa propre dérivée : d/dxex=exd/dx\ e^{x} = e^{x}. La fonction exponentielle est la seule fonction pouvant se vanter de cette propriété unique. (Pour être plus précis, tout multiple de la fonction exponentielle, y compris zéro, possède cette qualité.)

La Fonction Exponentielle est sa Propre Dérivée

ddxex=ex.\frac{d}{dx}\, e^{x} = e^{x}.

C’est cette propriété spéciale de la fonction exponentielle qui la rend unique et lui fait apparaître si fréquemment dans les applications. Chaque fois que le taux de variation d’une valeur est proportionnel à la valeur elle-même, la fonction exponentielle surgira presque certainement quelque part dans les mathématiques qui décrivent la dynamique du système.

L’exemple que la plupart d’entre nous connaissent est celui des intérêts composés. Que P(t)P(t) soit la somme d’argent sur votre compte bancaire au temps tt ; supposons que la somme accumule des intérêts. Le taux de variation par intervalle de temps — le montant des intérêts gagnés — est proportionnel au montant d’argent sur votre compte. Plus vous avez d’argent, plus vous gagnez d’intérêts, et plus il croît vite. Ainsi, la fonction exponentielle s’introduit dans la finance avec l’équation P(t)=P0ertP(t) = P_{0}e^{rt}, qui décrit le montant d’argent à tout instant tt, en supposant qu’un montant initial P0P_{0} croît à un taux d’intérêt rr, les intérêts étant composés en continu.

Vous avez peut-être remarqué que la série de Taylor de exe^{x} ressemble frappamment aux représentations en série de sinx\sin x et cosx\cos x. Cette similarité laisse entrevoir une relation profonde et surprenante entre les fonctions exponentielles et les fonctions trigonométriques, que nous explorons dans l’ Exercice 11.

Nous espérons que cette brève rencontre avec les séries de Taylor, bien que quelque peu en dehors de notre fil conducteur principal, a éveillé votre intérêt pour un outil mathématique qui est très pratique, notamment pour son importance fondamentale dans toutes sortes de calculs d’approximation et numériques sur ordinateur. Nous espérons également que c’était un exemple non trivial intéressant de différenciation d’un polynôme. Cela nous a aussi donné l’occasion de discuter des dérivées des fonctions sinus, cosinus et exponentielle ; ces dérivées réapparaîtront dans les sections ultérieures.

11.4.7La règle de la chaîne

La règle de la chaîne est la dernière règle de différenciation que nous abordons ici. Elle nous dit comment déterminer le taux de variation d’une fonction quand l’argument de cette fonction est lui-même une autre fonction que nous savons différencier.

Dans la course entre la tortue et le lièvre, nous n’avons jamais vraiment réfléchi à ce que mesurait exactement notre fonction x(t)x(t), nous avons simplement dit que c’était la « position » du lièvre. Supposons que le circuit était en réalité un parcours sinueux avec des collines, des ponts et même une boucle verticale, et que la fonction que nous avons représentée et précédemment nommée x(t)x(t) mesure en fait la distance linéaire le long de ce chemin sinueux, plutôt que, disons, une position horizontale. Pour éviter les connotations horizontales associées au symbole xx, introduisons la variable ss, qui donne la distance le long du circuit (en furlongs, bien sûr).

Disons que nous avons une fonction y(s)y(s) qui décrit l’altitude du circuit à une distance donnée. La dérivée dy/dsdy/ds nous indique des informations très basiques sur le circuit à cet endroit. Une valeur nulle signifie que le circuit est plat à cet endroit, une valeur positive signifie que les coureurs montent une côte, et une grande valeur positive ou négative indique un endroit où le circuit est très pentu.

Considérons maintenant la fonction composite y(s(t))y(s(t)). Vous devriez pouvoir vous convaincre que cela nous indique l’altitude du lièvre à tout moment tt donné. La dérivée dy/dtdy/dt nous dit à quelle vitesse le lièvre se déplaçait verticalement, à un moment tt donné. C’est très différent de dy/dsdy/ds. Comment pourrait-on calculer dy/dtdy/dt ? Vous pourriez être tenté de dire que pour faire cette détermination, il suffit de savoir où se trouvait le lièvre sur le circuit au temps tt, puis la réponse est la pente du circuit à cet endroit. En symboles mathématiques, vous dites que la vitesse verticale est y(s(t))y^{\prime}(s(t)). Mais ce n’est pas correct. Par exemple, pendant que le lièvre faisait sa sieste (ds/dt=0ds/dt = 0), peu importe la pente du circuit ; puisqu’il ne se déplaçait pas dessus, sa vitesse verticale est nulle ! En fait, à un certain moment de la course il a fait demi-tour et a couru sur le circuit dans la mauvaise direction (ds/dt<0ds/dt < 0), donc sa vitesse verticale dy/dtdy/dt serait opposée à la pente du circuit dy/dsdy/ds. Et évidemment, s’il sprinte rapidement sur un endroit du circuit, sa vitesse verticale sera plus élevée que s’il marchait lentement sur ce même endroit. Mais de même, là où le circuit est plat, peu importe la vitesse à laquelle il court, sa vitesse verticale sera nulle. Nous voyons donc que la vitesse verticale du lièvre est le produit de sa vitesse (mesurée paramétriquement le long du circuit) et de la pente du circuit en ce point.

Cette règle est connue sous le nom de règle de la chaîne. Elle est particulièrement intuitive lorsqu’elle est écrite en notation de Leibniz, car les infinitésimaux dsds semblent « s’annuler ».

La Règle de la Chaîne de Différenciation

dydt=dydsdsdt.\begin{matrix} \frac{dy}{dt} & {= \frac{dy}{ds}\,\frac{ds}{dt}.} \\ \end{matrix}

Voici quelques exemples, utilisant des fonctions que nous savons maintenant différencier :

Exemples de la règle de la chaîne

ddtsin3x=3cos3x,ddtsin(x2)=2xcos(x2),ddtecosx+3x=(sinx+3)ecosx+3x,ddtesin3x+sin(x2)=(3cos3x+2xcos(x2))esin3x+sin(x2).\begin{matrix} {\frac{d}{dt}\,\sin 3x} & {= 3\cos 3x,} \\ {\frac{d}{dt}\,\sin(x^{2})} & {= 2x\cos(x^{2}),} \\ {\frac{d}{dt}\, e^{\cos x + 3x}} & {= ( - \sin x + 3)e^{\cos x + 3x},} \\ {\frac{d}{dt}\, e^{\sin 3x + \sin(x^{2})}} & {= (3\cos 3x + 2x\cos(x^{2}))e^{\sin 3x + \sin(x^{2})}.} \\ \end{matrix}

Nous allons mettre le calcul d’un point de vue purement mathématique de côté pour un moment et revenir à la cinématique. (Après tout, notre but en discutant du calcul était, comme Isaac Newton, d’améliorer notre compréhension de la mécanique.) Cependant, il ne tardera pas à revenir sur le tapis avec la discussion de l’intégrale et du théorème fondamental du calcul.

11.5Accélération

Nous avons beaucoup insisté sur la distinction entre vitesse instantanée et vitesse moyenne, et cette distinction est importante (et l’insistance est justifiée) quand la vitesse change continuellement. Dans de telles situations, nous pourrions vouloir connaître le taux auquel la vitesse change. Heureusement, nous venons d’apprendre la dérivée, dont la raison d’être est d’étudier les taux de variation. Quand nous prenons la dérivée d’une fonction de vitesse v(t)v(t), nous obtenons une nouvelle fonction décrivant à quelle vitesse la vitesse augmente ou diminue à cet instant. Ce taux de variation instantané est une grandeur importante en physique, et elle porte un nom familier : l’accélération.

Dans la conversation courante, le verbe « accélérer » signifie généralement « aller plus vite ». Cependant, en physique, le mot « accélération » a une signification plus générale et peut désigner tout changement de vitesse, pas seulement une augmentation de rapidité. En fait, un corps peut subir une accélération même quand sa rapidité est constante ! Comment est-ce possible ? La vitesse est une valeur vectorielle, ce qui signifie qu’elle a à la fois une magnitude et une direction. Si la direction de la vitesse change, mais pas sa magnitude (sa rapidité), nous disons que le corps éprouve une accélération. Une telle terminologie n’est pas de la simple chicane sur les mots ; l’accélération dans ce cas est une sensation très réelle qui serait ressentie, par exemple, par deux personnes assises à l’arrière d’une voiture qui zigzague et qui se retrouvent pressées l’une contre l’autre d’un côté. Nous en dirons plus sur cette situation particulière dans la Section 11.8.

Nous pouvons apprendre beaucoup sur l’accélération en nous demandant simplement quels types d’unités nous devons utiliser pour la mesurer. Pour la vitesse, nous avons utilisé les unités génériques L/TL/T, longueur unitaire par unité de temps. La vitesse est un taux de variation de position (LL) par unité de temps (TT), ce qui est logique. L’accélération est le taux de variation de la vitesse par unité de temps, et doit donc être exprimée en « vitesse unitaire par unité de temps ». En fait, les unités utilisées pour mesurer l’accélération sont L/T2L/T^{2}. Si vous êtes troublé par l’idée de « temps au carré », pensez-y plutôt comme à (L/T)/T(L/T)/T, ce qui rend plus explicite le fait qu’il s’agit d’une unité de vitesse (L/T)(L/T) par unité de temps.

Par exemple, un objet en chute libre près de la surface de la Terre accélère à un taux d’environ 32ft/s232\ {ft}/s^{2}, ou 9,8m/s29,8\ m/s^{2}. Disons que vous suspendez une bille métallique au bord de la Willis Tower.27 Vous lâchez la bille, et elle commence à accélérer, ajoutant 9,8m/s9,8\ m/s à sa vitesse vers le bas chaque seconde. (Nous ignorons la résistance de l’air.) Après, disons, 2,4 secondes, sa vitesse sera

2,4s×32fts2=76,8fts.2,4\ s \times 32\ \frac{ft}{s^{2}} = 76,8\ \frac{ft}{s}.

Plus généralement, la vitesse à un temps arbitraire tt d’un objet sous accélération constante est donnée par la formule linéaire simple

v(t)=v0+at,\begin{matrix} {v(t) = v_{0} + at,} \\ \end{matrix}

v0v_{0} est la vitesse initiale au temps t=0t = 0, et aa est l’accélération constante. Nous étudions le mouvement des objets en chute libre plus en détail dans la Section 11.6, mais d’abord, examinons une représentation graphique de l’accélération. La Figure 11.9 montre les graphes d’une fonction de position et des fonctions de vitesse et d’accélération correspondantes.

image

Figure 11.9 Graphes de la position, de la vitesse et de l’accélération en fonction du temps.

Vous devriez étudier la Figure 11.9 jusqu’à ce qu’elle vous soit compréhensible. Voici notamment quelques observations dignes d’intérêt :

Les accélérations subies par un objet peuvent varier en fonction du temps, et nous pouvons continuer ce processus de différenciation, ce qui donne encore une autre fonction du temps, que certains appellent la fonction de « secousse » (jerk). Dans ce livre, nous nous en tenons à la fonction de position et à ses deux premières dérivées. De plus, il est très instructif de considérer des situations dans lesquelles l’accélération est constante (ou du moins a une magnitude constante). C’est précisément ce que nous allons faire dans les prochaines sections.

La Section 11.6 considère les objets sous accélération constante, comme les objets en chute libre et les projectiles. Cela fournira un excellent contexte pour introduire l’intégrale, le complément de la dérivée, dans la Section 11.7. Ensuite, la Section 11.8 examine les objets se déplaçant sur un trajet circulaire, qui subissent une accélération de magnitude constante mais de direction changeant continuellement et pointant toujours vers le centre du cercle.

11.6Mouvement sous accélération constante

Examinons maintenant la trajectoire que suit un objet quand il accélère à un taux constant dans le temps. C’est un cas simple, mais fréquent, et important à bien comprendre. En fait, les équations du mouvement que nous présentons dans cette section sont parmi les équations mécaniques les plus importantes à connaître par cœur, surtout pour la programmation de jeux vidéo.

Avant de commencer, considérons un type de mouvement encore plus simple — le mouvement à vitesse constante. Le mouvement à vitesse constante est un cas particulier du mouvement à accélération constante — le cas où l’accélération est constamment nulle. Le mouvement d’une particule à vitesse constante est une équation linéaire intuitive, essentiellement identique à l’ Équation (9.1), l’équation d’un rayon. En une dimension, la position d’une particule en fonction du temps est

x(t)=x0+vt,\begin{matrix} {x(t) = x_{0} + vt,} \\ \end{matrix}

x0x_{0} est la position de la particule au temps t=0t = 0, et vv est la vitesse constante.

Considérons maintenant les objets se déplaçant avec une accélération constante. Nous avons déjà mentionné au moins un exemple important : quand ils sont en chute libre, accélérant en raison de la gravité. (Nous ignorons la résistance de l’air et toutes les autres forces.) Le mouvement en chute libre est souvent appelé mouvement balistique. Nous commençons en une dimension pour garder les choses simples. Notre objectif est une formule x(t)x(t) pour la position d’une particule à un moment donné.

Reprenons notre exemple du lancer illégal de billes depuis la Willis Tower. Établissons un référentiel où xx augmente dans la direction vers le bas, et x0=0x_{0} = 0. En d’autres termes, x(t)x(t) mesure la distance que l’objet a parcourue depuis sa hauteur de lâcher au temps tt. Nous supposons aussi pour l’instant que la vitesse initiale est v0=0ft/sv_{0} = 0\ {ft/s}, ce qui signifie que vous relâchez simplement la bille sans la lancer.

À ce stade, nous ne savons même pas quelle forme devrait prendre x(t)x(t), nous sommes donc un peu bloqués. La « porte d’entrée » vers cette solution semble verrouillée pour nous en ce moment, aussi essayons de contourner par l’arrière, en utilisant une approche similaire à celle que nous avons utilisée plus tôt pour définir la vitesse instantanée. Nous considérerons des façons d’approximer la réponse et observerons ce qui se passe quand les approximations deviennent de plus en plus précises.

6Tranches,Δt=0,40t0v0Δx0.000.000.000.4012.805.120.8025.6010.241.2038.4015.361.6051.2020.482.0064.0025.60Total76.8012Tranches,Δt=0,20t0v0Δx0.000.000.000.206.401.280.4012.802.560.6019.203.840.8025.605.121.0032.006.401.2038.407.681.4044.808.961.6051.2010.241.8057.6011.522.0064.0012.802.2070.4014.08Total84.4824Tranches,Δt=0,10t0v0Δx0.000.000.000.103.200.320.206.400.640.309.600.960.4012.801.280.5016.001.600.6019.201.920.7022.402.240.8025.602.560.9028.802.881.0032.003.201.1035.203.521.2038.403.841.3041.604.161.4044.804.481.5048.004.801.6051.205.121.7054.405.441.8057.605.761.9060.806.082.0064.006.402.1067.206.722.2070.407.042.3073.607.36Total88.32\begin{matrix} \begin{matrix} {{6\ Tranches},\Delta t = 0,40} \\ \\ \begin{matrix} t_{0} & v_{0} & {\Delta x} \\ 0.00 & 0.00 & 0.00 \\ 0.40 & 12.80 & 5.12 \\ 0.80 & 25.60 & 10.24 \\ 1.20 & 38.40 & 15.36 \\ 1.60 & 51.20 & 20.48 \\ 2.00 & 64.00 & 25.60 \\ & {Total} & 76.80 \\ \end{matrix} \\ \\ \\ {{12\ Tranches},\Delta t = 0,20} \\ \\ \begin{matrix} t_{0} & v_{0} & {\Delta x} \\ 0.00 & 0.00 & 0.00 \\ 0.20 & 6.40 & 1.28 \\ 0.40 & 12.80 & 2.56 \\ 0.60 & 19.20 & 3.84 \\ 0.80 & 25.60 & 5.12 \\ 1.00 & 32.00 & 6.40 \\ 1.20 & 38.40 & 7.68 \\ 1.40 & 44.80 & 8.96 \\ 1.60 & 51.20 & 10.24 \\ 1.80 & 57.60 & 11.52 \\ 2.00 & 64.00 & 12.80 \\ 2.20 & 70.40 & 14.08 \\ & {Total} & 84.48 \\ \end{matrix} \\ \end{matrix} & \begin{matrix} {{24\ Tranches},\Delta t = 0,10} \\ \\ \begin{matrix} t_{0} & v_{0} & {\Delta x} \\ 0.00 & 0.00 & 0.00 \\ 0.10 & 3.20 & 0.32 \\ 0.20 & 6.40 & 0.64 \\ 0.30 & 9.60 & 0.96 \\ 0.40 & 12.80 & 1.28 \\ 0.50 & 16.00 & 1.60 \\ 0.60 & 19.20 & 1.92 \\ 0.70 & 22.40 & 2.24 \\ 0.80 & 25.60 & 2.56 \\ 0.90 & 28.80 & 2.88 \\ 1.00 & 32.00 & 3.20 \\ 1.10 & 35.20 & 3.52 \\ 1.20 & 38.40 & 3.84 \\ 1.30 & 41.60 & 4.16 \\ 1.40 & 44.80 & 4.48 \\ 1.50 & 48.00 & 4.80 \\ 1.60 & 51.20 & 5.12 \\ 1.70 & 54.40 & 5.44 \\ 1.80 & 57.60 & 5.76 \\ 1.90 & 60.80 & 6.08 \\ 2.00 & 64.00 & 6.40 \\ 2.10 & 67.20 & 6.72 \\ 2.20 & 70.40 & 7.04 \\ 2.30 & 73.60 & 7.36 \\ & {Total} & 88.32 \\ \end{matrix} \\ \end{matrix} \\ \end{matrix}

Tableau 11.3Valeurs pour différents nombres de tranches

Rendons notre exemple un peu plus spécifique. Plus tôt, nous avons calculé qu’après 2,4 secondes de chute libre, la bille aurait une vitesse de v(2,4)=76,8ft/sv(2,4) = 76,8\ {ft/s}. Cependant, nous n’avons pas dit combien de distance elle avait parcouru pendant ce temps. Essayons de calculer cette distance, qui est x(2,4)x(2,4). Pour ce faire, nous découpons l’intervalle total de 2,4 secondes en un certain nombre de plus petites « tranches » de temps, et approximons la distance que parcourt la bille pendant chaque tranche. Nous pouvons approximer la distance totale parcourue comme la somme des distances parcourues pendant chaque tranche. Pour approximer la distance parcourue par la bille pendant une seule tranche, nous calculons d’abord la vitesse de la bille au début de la tranche en utilisant l’ Équation (11.13). Nous approximons ensuite la distance parcourue pendant la tranche en insérant cette vitesse comme vitesse constante pour la tranche dans l’ Équation (11.14).

Le Tableau 11.3 montre les valeurs tabulées pour 6, 12 et 24 tranches. Pour chaque tranche, t0t_{0} désigne le temps de début de la tranche, v0v_{0} est la vitesse au début de la tranche (calculée selon l’ Équation (11.13) comme v0=t0×32ft/s2v_{0} = t_{0} \times 32\ {ft}/s^{2}), Δt\Delta t est la durée de la tranche, et Δx\Delta x est notre approximation du déplacement pendant la tranche (calculée selon l’Équation (11.14) comme Δx=v0Δt\Delta x = v_{0}\Delta t).

Puisque chaque tranche a une vitesse initiale différente, nous tenons compte du fait que la vitesse change sur l’ensemble de l’intervalle. (En fait, le calcul de la vitesse initiale de la tranche n’est pas une approximation — il est exact.) Cependant, comme nous ignorons le changement de vitesse à l’intérieur d’une tranche, notre réponse n’est qu’une approximation. En prenant de plus en plus de tranches, nous obtenons des approximations de plus en plus précises, bien qu’il soit difficile de voir vers quelle valeur ces approximations convergent. Examinons le problème graphiquement pour voir si nous pouvons avoir un aperçu.

Dans la Figure 11.10, chaque rectangle représente un intervalle de temps dans notre approximation. Notez que la distance parcourue pendant un intervalle est la même que l’aire du rectangle correspondant :

(aire du rectangle)=(largeur du rectangle)×(hauteur du rectangle)=(durée de la tranche)×(vitesse utilisée pour la tranche)=(déplacement pendant la tranche).\begin{matrix} {(\text{aire\ du\ rectangle})} & {= (\text{largeur\ du\ rectangle}) \times (\text{hauteur\ du\ rectangle})} \\ & {= (\text{durée\ de\ la\ tranche}) \times (\text{vitesse\ utilisée\ pour\ la\ tranche})} \\ & {= (\text{déplacement\ pendant\ la\ tranche}).} \\ \end{matrix}

image image
image

Figure 11.10Représentation graphique du Tableau 11.3

Voici maintenant l’observation clé. En augmentant le nombre de tranches, l’aire totale des rectangles se rapproche de plus en plus de l’aire du triangle sous la courbe de vitesse. À la limite, en prenant un nombre infini de rectangles, les deux aires seront égales. Puisque le déplacement total de la bille en chute libre est égal à l’aire totale des rectangles, qui est égale à l’aire sous la courbe, nous sommes conduits à une découverte importante.

La distance parcourue est égale à l’aire sous la courbe de vitesse.

Nous sommes arrivés à cette conclusion en utilisant un argument de limite très similaire à celui que nous avons utilisé pour définir la vitesse instantanée — nous considérons comment une série d’approximations converge dans la limite quand l’erreur d’approximation tend vers zéro.

Notez que nous n’avons fait aucune hypothèse dans cet argument sur v(t)v(t). Dans l’exemple en cours, c’est une simple fonction linéaire, et le graphe est une droite ; cependant, vous devriez pouvoir vous convaincre que cette procédure fonctionnera pour toute fonction de vitesse arbitraire.28 Cet argument de limite est un outil formalisé en calcul connu sous le nom d’intégrale de Riemann, que nous considérerons dans la Section 11.7. Ce sera aussi le moment approprié pour considérer le cas général de tout v(t)v(t). Cependant, puisqu’il y a tant de choses à apprendre de cet exemple spécifique, gardons-le simple aussi longtemps que possible.

Rappelons la question à laquelle nous essayons de répondre : quelle distance parcourt un objet après avoir été lâché à une vitesse initiale nulle, puis accéléré par la gravité pendant 2,4 secondes à un taux constant de 32ft/s232\ {ft}/s^{2} ? Comment cette nouvelle réalisation sur l’équivalence entre la distance parcourue et l’aire sous le graphe de v(t)v(t) nous aide-t-elle ? Dans ce cas particulier, v(t)v(t) est une simple fonction linéaire, et l’aire sous la courbe de t=0t = 0 à t=2,4t = 2,4 est un triangle. C’est une forme dont nous savons facilement calculer l’aire. La base de ce triangle a une longueur 2,4s2,4\ s, et la hauteur est v(2,4)=76,8ft/sv(2,4) = 76,8\ {ft/s}, donc l’aire est

base×hauteur2=2,4s×76,8ft/s2=92,16ft.\frac{\text{base} \times \text{hauteur}}{2} = \frac{2,4\ s \times 76,8\ {ft/s}}{2} = 92,16\ {ft}.

Ainsi, en seulement 2,4 secondes, la bille avait déjà chuté de plus de 92 pieds !

Cela résout le problème spécifique posé, mais soyons plus généraux. Rappelons que l’objectif plus large était une équation cinématique x(t)x(t) qui prédit la position d’un objet en fonction d’une position initiale et d’une vitesse initiale quelconques. Premièrement, remplaçons la constante 2,4 par un temps arbitraire tt. Ensuite, supprimons l’hypothèse que l’objet a initialement une vitesse nulle, et permettons une vitesse initiale arbitraire v0v_{0}. Cela signifie que l’aire sous la courbe v(t)v(t) n’est plus un triangle — c’est un triangle posé sur un rectangle, comme indiqué dans la Figure 11.11.

image

Figure 11.11 Calcul du déplacement au temps tt, avec une vitesse initiale v0v_{0} et une accélération constante aa

Le rectangle a une base tt et une hauteur v0v_{0}, et son aire représente la distance qui serait parcourue s’il n’y avait pas d’accélération. Le triangle au-dessus du rectangle a aussi une base tt, et la hauteur est atat, la différence de v(t)v(t) par rapport à la vitesse initiale résultant de l’accélération au taux aa sur la durée de tt secondes. En additionnant ces deux parties, on obtient le déplacement total, que nous notons Δx\Delta x :

Δx=(Aire du rectangle)+(Aire du triangle)=(Base durectangle)(Hauteur durectangle)+12(Base dutriangle)(Hauteur dutriangle)=(t)(v0)+(1/2)(t)(at)=v0t+(1/2)at2.\begin{matrix} {\Delta x} & {= (\text{Aire\ du\ rectangle}) + (\text{Aire\ du\ triangle})} \\ & {= \begin{pmatrix} \text{Base du} \\ \text{rectangle} \\ \end{pmatrix}\begin{pmatrix} \text{Hauteur du} \\ \text{rectangle} \\ \end{pmatrix} + \frac{1}{2}\begin{pmatrix} \text{Base du} \\ \text{triangle} \\ \end{pmatrix}\begin{pmatrix} \text{Hauteur du} \\ \text{triangle} \\ \end{pmatrix}} \\ & {= (t)(v_{0}) + (1/2)(t)(at)} \\ & {= v_{0}t + (1/2)at^{2}.} \\ \end{matrix}

Nous venons de dériver une équation très utile, mettons-la en valeur pour que les lecteurs qui parcourent le texte la remarquent.

Formule du Déplacement avec Vitesse Initiale et Accélération Constante

Δx=v0t+(1/2)at2.\begin{matrix} {\Delta x = v_{0}t + (1/2)at^{2}.} \\ \end{matrix}

L’Équation (11.15) fait partie d’une poignée d’équations dans ce livre qui valent la peine d’être mémorisées. Elle est très utile pour résoudre des problèmes pratiques qui se posent dans les simulations physiques.29

Il est courant que nous n’ayons besoin que du déplacement Δx\Delta x, et que la position absolue x(t)x(t) n’ait pas d’importance. Cependant, puisque la fonction x(t)x(t) était notre objectif déclaré, nous pouvons facilement l’exprimer en termes de l’Équation (11.15) en ajoutant le déplacement à notre position initiale, que nous notons x0x_{0} :

x(t)=x0+Δx=x0+v0t+(1/2)at2.x(t) = x_{0} + \Delta x = x_{0} + v_{0}t + (1/2)at^{2}.

Travaillons sur quelques exemples pour montrer les types de problèmes qui peuvent être résolus en utilisant l’ Équation (11.15) et ses variantes. Un scénario tentant est de laisser notre bille toucher le sol. La plate-forme d’observation au 103e étage de la Willis Tower est à 1 353 pieds au-dessus du trottoir. Si elle est lâchée de cette hauteur, combien de temps mettra-t-elle pour tomber jusqu’en bas ? En résolvant l’ Équation (11.15) pour tt, nous obtenons

Résolution pour le temps

Δx=v0t+(1/2)at20=(a/2)t2+v0tΔxt=v0±v024(a/2)(Δx)2(a/2)(formule quadratique)t=v0±v02+2aΔxa.\begin{matrix} {\Delta x} & {= v_{0}t + (1/2)at^{2}} & & \\ 0 & {= (a/2)t^{2} + v_{0}t - \Delta x} & & \\ t & {= \frac{- v_{0} \pm \sqrt{v_{0}^{2} - 4(a/2)( - \Delta x)}}{2(a/2)}} & & \text{(formule\ quadratique)} \\ t & {= \frac{- v_{0} \pm \sqrt{v_{0}^{2} + 2a\Delta x}}{a}.} & & \\ \end{matrix}

L’Équation (11.16) est une équation générale très utile. En insérant les valeurs spécifiques à ce problème, nous obtenons

t=v0±v02+2aΔxa=(0)±(0)2+2(32ft/s2)(1353ft)32ft/s2\begin{matrix} t & {= \frac{- v_{0} \pm \sqrt{v_{0}^{2} + 2a\Delta x}}{a}} \\ & {= \frac{- (0) \pm \sqrt{(0)^{2} + 2(32\ {ft/s}^{2})(1\, 353\ {ft})}}{32\ {ft/s}^{2}}} \\ \end{matrix}

=±86592(ft/s)232ft/s2±294,3ft/s32ft/s2±9,197s.\begin{matrix} & {= \pm \frac{\sqrt{86\, 592\ ({ft/s})^{2}}}{32\ {ft/s}^{2}}} \\ & {\approx \pm \frac{294,3\ {ft/s}}{32\ {ft/s}^{2}}} \\ & {\approx \pm 9,197\ s.} \\ \end{matrix}

La racine carrée dans l’ Équation (11.16) introduit la possibilité de deux solutions. Nous utilisons toujours la racine qui donne une valeur positive pour tt.30

Naturellement, une personne qui s’adonne au lancer de billes depuis de grandes hauteurs s’intéresse aux dégâts qu’elle peut causer, aussi la prochaine question logique est : « À quelle vitesse la bille voyage-t-elle quand elle touche le trottoir ? » Pour répondre à cette question, nous insérons le temps de trajet total dans l’ Équation (11.13) :

v(t)=v0+at=0ft/s+(32ft/s2)(9,197s)=294,3ft/s.v(t) = v_{0} + at = 0\ {ft/s} + (32\ {ft/s}^{2})(9,197\ s) = 294,3\ {ft/s}.

Si nous ignorons la résistance de l’air, au moment de l’impact, la bille voyage à une vitesse qui couvre une distance d’environ un terrain de football en une seconde ! On comprend pourquoi ce que nous faisons dans notre imagination est illégal dans la réalité. Continuons quand même.

Supposons maintenant qu’au lieu de simplement lâcher la bille, nous lui donnons une vitesse initiale (nous la lançons vers le haut ou vers le bas). C’était notre libre choix de décider si le haut ou le bas est positif dans ces exemples, et nous avons choisi +x+ x pour la direction vers le bas, ce qui signifie que la vitesse initiale sera négative. Quelle doit être la vitesse initiale pour que la bille reste en l’air quelques secondes de plus, disons un total de 12 secondes ? Une fois encore, nous allons d’abord manipuler l’ Équation (11.15) pour obtenir une solution générale ; cette fois, nous résolvons pour v0v_{0} :

Résolution pour la vitesse initiale

Δx=v0t+(1/2)at2,v0t=Δx+(1/2)at2,v0=Δx/t(1/2)at.\begin{matrix} {\Delta x} & {= v_{0}t + (1/2)at^{2},} \\ {- v_{0}t} & {= - \Delta x + (1/2)at^{2},} \\ v_{0} & {= \Delta x/t - (1/2)at.} \\ \end{matrix}

En insérant maintenant les valeurs spécifiques à notre problème, nous obtenons

v0=Δx/t(1/2)at=(1353ft)/(12,0s)(1/2)(32ft/s2)(12,0s)=112,8ft/s192ft/s=79,2ft/s.\begin{matrix} v_{0} & {= \Delta x/t - (1/2)at} \\ & {= (1\, 353\ {ft})/(12,0\ s) - (1/2)(32\ {ft/s}^{2})(12,0\ s)} \\ & {= 112,8\ {ft/s} - 192\ {ft/s}} \\ & {= - 79,2\ {ft/s}.} \\ \end{matrix}

Notez que le résultat est négatif, indiquant une vitesse vers le haut. Si nous donnons cette vitesse initiale à la bille, nous pourrions nous demander combien de temps il lui faut pour revenir à sa position initiale. En utilisant l’ Équation (11.16) et en posant Δx=0\Delta x = 0, nous obtenons

t=v0±v02+2aΔxa=(79,2ft/s)±(79,2ft/s)2+2(32ft/s2)(0ft)32ft/s2=79,2ft/s±(79,2ft/s)232ft/s2=79,2ft/s±79,2ft/s32ft/s2=0ou4,95s.\begin{matrix} t & {= \frac{- v_{0} \pm \sqrt{v_{0}^{2} + 2a\Delta x}}{a}} \\ & {= \frac{- ( - 79,2\ {ft/s}) \pm \sqrt{( - 79,2\ {ft/s})^{2} + 2(32\ {ft/s}^{2})(0\ {ft})}}{32\ {ft/s}^{2}}} \\ & {= \frac{79,2\ {ft/s} \pm \sqrt{( - 79,2\ {ft/s})^{2}}}{32\ {ft/s}^{2}}} \\ & {= \frac{79,2\ {ft/s} \pm 79,2\ {ft/s}}{32\ {ft/s}^{2}}} \\ & {= 0\ \text{ou}\ 4,95\ s.} \\ \end{matrix}

Il n’est pas surprenant que t=0t = 0 soit une solution ; nous résolvions pour les valeurs de temps quand la bille se trouvait à sa position initiale.

Examinez le graphe de la Figure 11.12, qui représente la position et la vitesse d’un objet se déplaçant sous accélération constante aa avec une vitesse initiale v0v_{0}, où v0v_{0} et aa ont des signes opposés. Faisons trois observations clés. Bien que nous utilisions des termes comme « hauteur », qui sont spécifiques au mouvement de projectile, des affirmations similaires sont vraies chaque fois que les signes de v0v_{0} et aa sont opposés.

image

Figure 11.12Mouvement d’un projectile

La première observation est que le projectile atteint sa hauteur maximale, notée xmaxx_{max}, quand l’accélération a consommé toute la vitesse et que v(t)=0v(t) = 0. Il est facile de résoudre le temps auquel cela se produit en utilisant l’ Équation (11.13), v(t)=v0+atv(t) = v_{0} + at :

Temps pour atteindre l’apogée

v(t)=0,v0+at=0,t=v0/a.\begin{matrix} {v(t)} & {= 0,} \\ {v_{0} + at} & {= 0,} \\ t & {= - v_{0}/a.} \\ \end{matrix}

Nous sommes actuellement en une dimension et ne considérons que la hauteur. Mais si nous sommes en plus d’une dimension, seule la vitesse parallèle à l’accélération doit s’annuler. Il pourrait y avoir une vitesse horizontale, par exemple. Nous discutons du mouvement de projectile en plus d’une dimension dans un moment.

La deuxième observation est que le temps qu’il faut à l’objet pour descendre de son altitude maximale à son altitude initiale, noté tet_{e} dans la Figure 11.12, est le même que le temps nécessaire pour atteindre le maximum. En d’autres termes, le projectile atteint son apogée à te/2t_{e}/2.

La troisième et dernière observation est que la vitesse à t=tet = t_{e}, que nous avons notée vev_{e}, a la même magnitude que la vitesse initiale v0v_{0}, mais de signe opposé.

Avant d’examiner le mouvement de projectile en plus d’une dimension, résumons les formules que nous avons dérivées dans cette section. Les deux premières sont les seules qui valent la peine d’être mémorisées ; les autres peuvent en être déduites.

Résumé des Équations Cinématiques pour l’Accélération Constante

v(t)=v0+at,Δx=v0t+(1/2)at2,x(t)=x0+Δx=x0+v0t+(1/2)at2,v0=Δx/t(1/2)at,t=v0±v02+2aΔxa,a=2Δxv0tt2.\begin{matrix} {v(t)} & {= v_{0} + at,} \\ {\Delta x} & {= v_{0}t + (1/2)at^{2},} \\ {x(t)} & {= x_{0} + \Delta x = x_{0} + v_{0}t + (1/2)at^{2},} \\ v_{0} & {= \Delta x/t - (1/2)at,} \\ t & {= \frac{- v_{0} \pm \sqrt{v_{0}^{2} + 2a\Delta x}}{a},} \\ a & {= 2\frac{\Delta x - v_{0}t}{t^{2}}.} \\ \end{matrix}

L’extension des idées de la section précédente en 2D ou 3D est essentiellement une question de passage à la notation vectorielle ; xx, vv et aa deviennent respectivement 𝐩\mathbf{p}, 𝐯\mathbf{v} et 𝐚\mathbf{a}.31 Bien sûr, le temps tt reste un scalaire :

Équations du mouvement sous accélération constante, en forme vectorielle

𝐯(t)=𝐯0+t𝐚,Δ𝐩=𝐯0t+(t2/2)𝐚,𝐩(t)=𝐩0+Δ𝐩=𝐩0+t𝐯0+(t2/2)𝐚,𝐯0=Δ𝐩/t(1/2)at,𝐚=2Δ𝐩t𝐯0t2.\begin{matrix} {\mathbf{v}(t)} & {= \mathbf{v}_{0} + t\mathbf{a},} \\ {\Delta\mathbf{p}} & {= \mathbf{v}_{0}t + (t^{2}/2)\mathbf{a},} \\ {\mathbf{p}(t)} & {= \mathbf{p}_{0} + \Delta\mathbf{p} = \mathbf{p}_{0} + t\mathbf{v}_{0} + (t^{2}/2)\mathbf{a},} \\ \mathbf{v}_{0} & {= \Delta\mathbf{p}/t - (1/2)at,} \\ \mathbf{a} & {= 2\frac{\Delta\mathbf{p} - t\mathbf{v}_{0}}{t^{2}}.} \\ \end{matrix}

Notez que nous n’avons pas créé une version vectorielle de l’ Équation (11.17) ; nous y reviendrons dans un moment.

Ce changement de notation apparemment trivial cache en réalité deux faits assez profonds. Premièrement, au sens algébrique, la notation vectorielle n’est vraiment qu’un raccourci pour des ensembles d’équations scalaires parallèles pour xx, yy et zz. Le point important est que les trois coordonnées (cartésiennes) sont complètement indépendantes les unes des autres. Par exemple, nous pouvons faire des calculs concernant yy et ignorer complètement les autres dimensions, à condition que l’hypothèse d’accélération constante soit vérifiée pour le mouvement de l’objet. Sans l’indépendance des coordonnées, nous ne pourrions pas effectuer ce changement de notation. Le deuxième fait caché dans cette notation est que, lorsque nous considérons les vecteurs des équations ci-dessus comme des entités géométriques plutôt qu’algébriques, le système de coordonnées particulier utilisé pour décrire ces vecteurs est sans importance. Nous n’avons même pas besoin d’en spécifier un. Bien sûr, c’est un principe de base de la physique : la Nature ne connaît pas le système de coordonnées que vous utilisez.

Nous avons pu passer de 1D à 3D en mettant simplement quelques lettres en gras, grâce à l’indépendance des coordonnées. Cependant, il y a quelques remarques supplémentaires sur le mouvement de projectile en plusieurs dimensions, car il y a des situations où nous devons considérer les effets de toutes les coordonnées en même temps. Une situation a déjà été suggérée par l’absence d’une équation vectorielle correspondant à l’Équation (11.17). Autrement dit, comment pourrait-on résoudre pour le temps tt en connaissant le déplacement Δ𝐩\Delta\mathbf{p}, l’accélération 𝐚\mathbf{a} et la vitesse initiale 𝐯0\mathbf{v}_{0} ? En une dimension, le projectile est « confiné » et ne peut pratiquement pas éviter d’atteindre la cible impliquée par Δx\Delta x.32 Mais en deux dimensions ou plus, la situation est plus compliquée. L’augmentation de complexité liée à l’augmentation des dimensions est analogue au calcul de l’intersection de deux rayons (voir la Section A.8). En 2D, deux rayons quelconques doivent s’intersecter à moins d’être parallèles, tandis qu’en 3D, il existe la possibilité de rayons gauches, qui ne sont pas parallèles mais ne s’intersectent pas.

Par exemple, nous avons calculé plus tôt le temps qu’il faudrait à une bille lâchée d’une grande hauteur pour toucher le trottoir en dessous, ce qui est un problème unidimensionnel. Le problème tridimensionnel correspondant serait d’essayer de faire tomber la bille dans un seau qui peut se déplacer librement sur le trottoir. Disons que le seau est à notre gauche. Notre vitesse initiale devra alors avoir une composante vers la gauche, sinon la bille ne tombera pas dans le seau. Autre signe que le cas multidimensionnel est plus complexe que le 1D : une traduction directe de l’ Équation (11.17) en forme vectorielle donne des opérations insensées comme prendre la racine carrée d’un vecteur et diviser un vecteur par un autre.

La clé pour résoudre ce problème est de réaliser que les changements horizontaux (soit à la position du seau, soit à la vitesse initiale de la bille) n’affectent pas le temps qu’il faut à la bille pour atteindre le trottoir. C’est parce que les coordonnées sont indépendantes les unes des autres. La vitesse et l’accélération horizontales n’interagissent pas avec la vitesse et l’accélération verticales. Pour être précis, passons à notre système de coordonnées 3D standard, où +y+ y pointe vers le haut et xx et zz sont dans le plan horizontal. Le temps qu’il faut à la bille pour atteindre l’altitude du seau dépend uniquement des équations concernant yy ; les coordonnées xx et zz peuvent être ignorées à cette fin. 33 En d’autres termes, calculer le moment où un projectile atteindra une cible est toujours un calcul unidimensionnel — il suffit de choisir quelle direction utiliser. Nous pouvons appliquer l’ Équation (11.17) pour résoudre un temps d’impact tt. Mais cette solution n’est qu’une proposition. Nous savons que si le projectile devait toucher la cible, il le ferait à ce moment. Pour vérifier que nous avons bien touché la cible, nous devons insérer ce temps d’impact présumé dans l’ Équation (11.19) pour voir où se trouvera le projectile à ce moment, et vérifier que la position du projectile est dans les tolérances appropriées.

Parlons un peu plus de ce que signifie exactement « choisir quelle direction utiliser », comme indiqué dans le paragraphe précédent. Dans des cas de mouvement de projectile simples, comme l’exemple de la bille, où la gravité est l’accélération constante, la direction à choisir est évidente : utiliser la direction de la gravité. De plus, parce que les systèmes de coordonnées sont choisis de sorte que le « haut » soit l’un des axes cardinaux, le processus de résolution d’un problème unidimensionnel dans cette direction est une affaire triviale consistant à extraire la coordonnée cartésienne appropriée et à ignorer les autres. En général, cependant, la situation peut être plus compliquée. Mais avant de discuter des détails du cas général, il y a encore quelques choses à dire sur cette situation particulière très importante et fréquente.

image

Figure 11.13Mouvement de projectile

Pour étudier le mouvement de projectile où l’accélération est uniquement due à la gravité, qui est une constante agissant le long d’un axe cardinal, établissons un espace de coordonnées 2D où +y+ y est vers le haut et xx est l’axe horizontal. Sans perte de généralité, nous pouvons faire pivoter notre plan de sorte qu’il contienne la vitesse initiale, et donc toute la trajectoire de la particule. Nous choisissons +x+ x dans la direction horizontale de la vitesse initiale. Nous simplifions aussi les choses en plaçant l’origine à la position initiale de l’objet. Cette notation (ainsi que quelques autres éléments dont nous aurons besoin dans un moment) est illustrée dans la Figure 11.13.

En examinant la notation de la Figure 11.13, nous voyons que nous pouvons exprimer la position de la particule en fonction du temps soit comme 𝐩(t)\mathbf{p}(t), soit en se référant à une coordonnée individuelle avec x(t)x(t) et y(t)y(t). La vitesse instantanée (non représentée dans le diagramme) peut être notée en forme vectorielle comme 𝐯(t)\mathbf{v}(t) ou en notation dérivée comme 𝐩˙(t)\overset{˙}{\mathbf{p}}(t). Les composantes scalaires de la vitesse seront notées en notation dérivée comme x˙(t)\overset{˙}{x}(t) et y˙(t)\overset{˙}{y}(t). La position et la vitesse initiales seront notées en ajoutant un indice 0 (y˙0{\overset{˙}{y}}_{0} est la vitesse verticale initiale). Nous notons l’accélération due à la gravité comme gg ou 𝐠\mathbf{g}.

Énonçons les équations de vitesse et de position en utilisant la notation qui vient d’être décrite :

𝐩˙(t)=𝐯0+t𝐠,x˙(t)=x˙0,y˙(t)=y˙0+gt,𝐩(t)=t𝐯0+(t2/2)𝐠,x(t)=tx˙0,y(t)=ty˙0+(1/2)gt2.\begin{matrix} {\overset{˙}{\mathbf{p}}(t)} & {= \mathbf{v}_{0} + t\mathbf{g},} & {\overset{˙}{x}(t)} & {= {\overset{˙}{x}}_{0},} & {\overset{˙}{y}(t)} & {= {\overset{˙}{y}}_{0} + gt,} \\ {\mathbf{p}(t)} & {= t\mathbf{v}_{0} + (t^{2}/2)\mathbf{g},} & {x(t)} & {= t{\overset{˙}{x}}_{0},} & {y(t)} & {= t{\overset{˙}{y}}_{0} + (1/2)gt^{2}.} \\ \end{matrix}

Les distances étiquetées hh et dd dans la Figure 11.13 sont souvent d’intérêt ; il s’agit respectivement de la hauteur à l’apogée et de la distance horizontale parcourue. Comme discuté plus tôt dans un contexte unidimensionnel, la hauteur maximale est atteinte quand toute la vitesse initiale dans la direction montante a été consommée par la gravité, c’est-à-dire quand y˙(t)=0\overset{˙}{y}(t) = 0. Cela se produit au temps

Temps pour atteindre l’apogée

ta=y˙0/g,t_{a} = - {\overset{˙}{y}}_{0}/g,

et à ce moment, la hauteur est égale à

Altitude à l’apogée

h=y(ta)=tay˙0+(1/2)gta2=(y˙0/g)y˙0+(1/2)g(y˙0/g)2=(y˙02/g)+(1/2)(y˙02/g)=y˙02/2g.\begin{matrix} h & {= y(t_{a}) = t_{a}{\overset{˙}{y}}_{0} + (1/2)gt_{a}^{2}} \\ & {= ( - {\overset{˙}{y}}_{0}/g){\overset{˙}{y}}_{0} + (1/2)g( - {\overset{˙}{y}}_{0}/g)^{2}} \\ & {= ( - {\overset{˙}{y}}_{0}^{2}/g) + (1/2)({\overset{˙}{y}}_{0}^{2}/g)} \\ & {= - {\overset{˙}{y}}_{0}^{2}/2g.} \\ \end{matrix}

Nous avons dit plus tôt que le temps qu’il faut à l’objet pour redescendre à son altitude initiale (que nous avons noté tet_{e}) était le double du temps nécessaire pour atteindre son apogée ; mais à ce moment, nous nous sommes simplement appuyés sur un diagramme. Cette fois, vérifions-le algébriquement :

Temps pour revenir à l’altitude initiale.

y(t)=ty˙0+(1/2)gt2,0=tey˙0+(1/2)gte2,(la position initiale est à l’origine)(1/2)gte2=tey˙0,te=2y˙0/g.(diviser par (1/2)gte)\begin{matrix} {y(t)} & {= t{\overset{˙}{y}}_{0} + (1/2)gt^{2},} & & \\ 0 & {= t_{e}{\overset{˙}{y}}_{0} + (1/2)gt_{e}^{2},} & & \text{(la\ position\ initiale\ est\ à\ l'origine)} \\ {- (1/2)gt_{e}^{2}} & {= t_{e}{\overset{˙}{y}}_{0},} & & \\ t_{e} & {= - 2{\overset{˙}{y}}_{0}/g.} & & {\text{(diviser\ par~}{- (1/2)gt_{e}}\text{)}} \\ \end{matrix}

Comme prévu, le temps de vol tet_{e} est le double du temps nécessaire pour atteindre l’apogée. Calculons maintenant dd, la distance horizontale parcourue :

Distance horizontale parcourue

d=x(te)=tex˙0=2y˙0x˙0/g.d = x(t_{e}) = t_{e}{\overset{˙}{x}}_{0} = - 2{\overset{˙}{y}}_{0}{\overset{˙}{x}}_{0}/g.

Bien sûr, tet_{e} et dd sont basés sur l’hypothèse que nous voulons savoir quand le projectile revient à son altitude initiale. C’est important pour le lancement d’un projectile depuis un plan de sol plat. Si le projectile n’est pas lancé depuis le sol, ou si le sol n’est pas plat, nous devrons considérer où la parabole intersecte le plan du sol.

Nous souhaitons souvent spécifier la vitesse initiale en termes d’angle et de vitesse, plutôt qu’en termes de vitesses le long de chaque axe. Autrement dit, nous souhaitons utiliser des coordonnées polaires plutôt que cartésiennes. Comme indiqué dans la Figure 11.13, nous notons la vitesse de lancement initiale comme s0s_{0} (qui est égale à la magnitude de 𝐯0\mathbf{v}_{0}) et l’angle de lancement comme θ\theta. En convertissant la vitesse initiale des coordonnées cartésiennes aux coordonnées polaires (voir la Section 7.1.3 si vous ne vous souvenez plus comment), nous obtenons

x˙0=s0cosθ,y˙0=s0sinθ.{\overset{˙}{x}}_{0} = s_{0}\cos\theta,{\overset{˙}{y}}_{0} = s_{0}\sin\theta.

En insérant cela dans nos équations cinématiques (11.20) et (11.21), nous obtenons les équations du mouvement pour un projectile en termes d’angle et de vitesse de lancement :

x˙(t)=s0cosθ,y˙(t)=s0sinθ+gt,x(t)=ts0cosθ,y(t)=ts0sinθ+(1/2)gt2.\begin{matrix} {\overset{˙}{x}(t)} & {= s_{0}\cos\theta,} & {\overset{˙}{y}(t)} & {= s_{0}\sin\theta + gt,} \\ {x(t)} & {= ts_{0}\cos\theta,} & {y(t)} & {= ts_{0}\sin\theta + (1/2)gt^{2}.} \\ \end{matrix}

Nous pouvons aussi exprimer tet_{e}, hh et dd en termes de s0s_{0} et θ\theta :

Grandeurs importantes dans le mouvement de projectile, exprimées en termes d’angle et de vitesse de lancement

ta=y˙0/g=(s0sinθ)/g=s0(sinθ)/g,te=2y˙0/g=2(s0sinθ)/g=2s0(sinθ)/g,d=2y˙0x˙0/g=2(s0sinθ)(s0cosθ)/g=2s02(sinθ)(cosθ)/g,h=(1/2)y˙02/g=(1/2)(s0sinθ)2/g=s02(sin2θ)/2g.\begin{matrix} t_{a} & {= - {\overset{˙}{y}}_{0}/g = - (s_{0}\sin\theta)/g} & & {= - s_{0}(\sin\theta)/g,} \\ t_{e} & {= - 2{\overset{˙}{y}}_{0}/g = - 2(s_{0}\sin\theta)/g} & & {= - 2s_{0}(\sin\theta)/g,} \\ d & {= - 2{\overset{˙}{y}}_{0}{\overset{˙}{x}}_{0}/g = - 2(s_{0}\sin\theta)(s_{0}\cos\theta)/g} & & {= - 2s_{0}^{2}(\sin\theta)(\cos\theta)/g,} \\ h & {= - (1/2){\overset{˙}{y}}_{0}^{2}/g = - (1/2)(s_{0}\sin\theta)^{2}/g} & & {= - s_{0}^{2}(\sin^{2}\theta)/\, 2g.} \\ \end{matrix}

Ces équations sont très pratiques car elles capturent directement la relation entre les grandeurs « conviviales » que sont la vitesse de lancement, l’angle de lancement, le temps de vol et la distance de vol.

À ce stade, faisons une observation intéressante sur la relation entre la vitesse initiale s0s_{0} et la distance horizontale parcourue dd. C’est une relation quadratique, ce qui signifie que lorsque nous augmentons s0s_{0} d’un facteur kk, nous augmentons dd d’un facteur k2k^{2}. Il pourrait sembler plus naturel que la relation soit linéaire, ce qui signifierait que dd augmenterait du même facteur kk. Nous pouvons comprendre la relation quadratique en décomposant la vitesse initiale en ses composantes horizontale et verticale, notées précédemment x˙0{\overset{˙}{x}}_{0} et y˙0{\overset{˙}{y}}_{0} respectivement. Il n’est pas difficile de voir qu’augmenter x˙0{\overset{˙}{x}}_{0} augmentera dd du même facteur. Moins évident est le fait que la même chose est vraie pour y˙0{\overset{˙}{y}}_{0}. C’est vrai parce que la durée pendant laquelle l’objet est en l’air est proportionnelle à y˙0{\overset{˙}{y}}_{0}. Donc si nous augmentons la vitesse verticale, nous donnons à l’objet plus de temps pour voyager. Ainsi, tout facteur d’échelle que nous appliquons à ss affectera la distance deux fois : une fois en raison de l’augmentation de la vitesse horizontale due à x˙0{\overset{˙}{x}}_{0}, et une autre fois en raison de l’augmentation du temps de vol due à y˙0{\overset{˙}{y}}_{0}. Cela produit une relation quadratique entre ss et dd.

Revenons maintenant à une question que nous avions mise de côté plus tôt : comment déterminer le point d’impact pour des vecteurs arbitraires Δ𝐩\Delta\mathbf{p}, 𝐚\mathbf{a} et 𝐯0\mathbf{v}_{0} ? Nous avons dit que la clé était de « choisir une direction » et de résoudre un problème unidimensionnel dans cette direction. Si une direction cardinale est choisie, nous jetons simplement les autres coordonnées. Pour une direction arbitraire, nous projetons le problème sur une droite dans cette direction. Toute composante de déplacement, de vitesse ou d’accélération perpendiculaire à cette droite est écartée lors de la projection. Nous avons appris comment projeter sur une droite et mesurer le déplacement dans une direction particulière en utilisant le produit scalaire dans la Section 2.11. Il ne reste qu’à choisir une direction.

En supposant que le projectile atteint la cible, nous obtiendrons la même valeur pour tt quelle que soit la direction choisie. Mais cela ne signifie pas que le choix est sans importance. Par exemple, dans l’exemple de la bille, il serait catastrophique de choisir les directions +x+ x ou +z+ z, car il n’y a pas d’accélération dans l’une ou l’autre de ces directions et l’application de l’ Équation (11.17) entraînerait une division par zéro. Cela suggère la stratégie consistant à utiliser simplement 𝐚\mathbf{a} lui-même comme direction de projection. Pour ce faire, nous calculons le produit scalaire de chaque quantité vectorielle avec 𝐚\mathbf{a}, en effectuant les substitutions Δx=Δ𝐩𝐚\Delta x = \Delta\mathbf{p} \cdot \mathbf{a}, v=𝐯𝐚v = \mathbf{v} \cdot \mathbf{a}, et a=𝐚𝐚a = \mathbf{a} \cdot \mathbf{a}. Ces quantités scalaires peuvent ensuite être insérées dans l’ Équation (11.17). L’Exercice 10 explore cela plus en détail.

11.7L’intégrale

Nous venons de montrer que le déplacement total d’un objet dans un intervalle de temps est égal à l’aire sous le graphe de la vitesse de l’objet. Nous avons utilisé l’exemple de l’accélération constante, qui a un graphe simple, et l’aire était facile à calculer géométriquement. Nous n’avons pas poursuivi dans une plus grande généralité l’argument de limite qui nous a conduits à la surprenante équivalence, car ce cas particulier a de si remarquables applications. Nous sommes maintenant prêts à discuter de cas plus généraux. Le besoin de calculer une « sommation continue », où le taux de croissance est une fonction connue, est un concept courant en ingénierie et en sciences. L’outil du calcul utilisé pour calculer ces sommes est l’intégrale.

Si vous avez déjà étudié le calcul intégral et avez une bonne intuition de l’utilité de l’intégrale, vous pouvez passer directement à la Section 11.8, quand notre attention retourne au sujet de la mécanique. Cependant, si vous n’avez jamais fait de calcul intégral ou si votre intuition sur l’intégrale est un peu fragile, continuez à lire.

Il y a deux façons importantes d’aborder l’intégrale. La première consiste essentiellement à rendre la notion de « sommer beaucoup de petits éléments » un peu plus précise et à introduire une certaine formalisation mathématique. L’autre façon consiste à comparer l’intégrale à la dérivée. Il est important de comprendre les deux interprétations. L’intégrale est un peu plus difficile à appréhender que la dérivée, mais pour des raisons qui deviendront apparentes plus tard, elle joue un rôle beaucoup plus important dans les simulations physiques et de nombreuses autres domaines de la programmation de jeux vidéo. Comprendre ce que fait l’intégrale est très important, même si le vaste assortiment de techniques analytiques sur papier pour calculer des intégrales n’est pas très utile dans notre cas, étant remplacé par des techniques d’intégration numérique.

Transformons notre sommation informelle en notation mathématique, dans laquelle nous calculons l’aire sous la courbe f(x)f(x) dans l’intervalle axba \leq x \leq b. Nous partitionnons cet intervalle en nn tranches, chacune ayant la largeur Δx=(ba)/n\Delta x = (b - a)/n. Le ii-ième rectangle aura une coordonnée gauche xix_{i}, une hauteur égale à f(xi)f(x_{i}), et une aire de f(xi)Δxf(x_{i})\,\Delta x. En utilisant la notation de sommation, nous additionnons tous ces rectangles :

Airei=1nf(xi)Δx.{Aire} \approx \sum\limits_{i = 1}^{n}f(x_{i})\Delta x.

L’erreur de cette approximation diminue à mesure que nous augmentons le nombre de tranches nn, et désormais, à moins que vous ne soyez un nouveau venu, vous savez que nous devons prendre la limite encore une fois.34 En prenant la limite quand nn augmente sans borne et que les tranches deviennent infinitésimalement petites, nous obtenons notre définition de l’intégrale définie.

Intégrale Définie

abf(x)dx=limni=1nf(xi)Δx.\begin{matrix} {\int_{a}^{b}f(x)\ dx = \lim\limits_{n \rightarrow \infty}\sum\limits_{i = 1}^{n}f(x_{i})\Delta x.} \\ \end{matrix}

Dans cette équation,

Δx=(ba)/n,xi=a+iΔx.\begin{matrix} {\Delta x} & {= (b - a)/n,} \\ x_{i} & {= a + i\Delta x.} \\ \end{matrix}

L’Équation (11.22) se lit « L’intégrale de aa à bb de f(x)dxf(x)\ dx. » Certains lisent dxdx comme « par rapport à xx. » La grande similitude dans la notation entre les membres gauche et droit de l’ Équation (11.22) est intentionnelle. Tout comme pour la dérivée, le pas fini Δx\Delta x devient l’infinitésimal dxdx. Le symbole sigma \sum utilisé pour les sommations discrètes est remplacé par le symbole \int, un S allongé que Leibniz voulait faire signifier « sommation ».35 Les aa et bb sont connus comme les « bornes d’intégration » et définissent les points de départ et d’arrivée. La fonction étant intégrée est appelée l’intégrande.

Une intégrale définie comme somme de « tranches verticales » de cette façon est connue sous le nom d’intégrale de Riemann. C’est la définition la plus courante, mais pas la plus générale. En fait, notre définition n’est pas tout à fait aussi générale que la définition typique d’une intégrale de Riemann. Le lecteur attentif peut remarquer que Δx\Delta x est une constante et pourrait être sorti devant la sommation, la rendant Δxi=1nf(xi)\Delta x\sum\limits_{i = 1}^{n}f(x_{i}). Cela fonctionne dans ce cas parce que nous utilisons une partition régulière, et toutes les tranches ont la même largeur. En général, cependant, cette restriction n’est pas nécessaire. La définition traditionnelle de l’intégrale de Riemann prend la limite quand la largeur de la plus grande tranche tend vers zéro. Notre définition est certainement assez puissante pour les fonctions bien comportées avec lesquelles nous travaillons, mais des définitions plus puissantes sont nécessaires pour intégrer des fonctions plus ésotériques. De plus, vous pourriez vous demander pourquoi nous calculons l’aire du rectangle en utilisant la valeur de la fonction au côté gauche du rectangle, plutôt que, par exemple, le point central. Ce serait certainement plus précis. À des fins théoriques de définition de l’intégrale de Riemann, ces choix deviennent identiques dans la limite et nous sommes donc libres de faire le choix que nous voulons. Cependant, lors de l’approximation numérique d’intégrales, il est utile de considérer de telles options.

11.7.1Exemples d’intégrales

À ce stade, nous avons introduit suffisamment de notation et de terminologie pour examiner quelques exemples d’intégrales. Nous aimerions le faire avant d’aller plus loin dans les détails mathématiques. De nombreuses applications de l’intégrale en programmation de jeux vidéo (et dans de nombreuses autres disciplines d’ingénierie) se conçoivent plus directement non pas comme une aire sous une courbe, mais comme un « total cumulatif ». Pensez à un compteur électrique. À tout moment donné, le compteur augmente à un taux déterminé par la quantité d’électricité consommée à cet instant. Le compteur est un total cumulatif continu, et nous disons qu’il intègre le taux de consommation. Quand la climatisation se met en marche, le taux de consommation augmente et le compteur avance plus vite ; la nuit, quand toutes les lumières sont éteintes et les fenêtres ouvertes parce qu’il fait beau dehors, la consommation est la plus faible et le compteur tourne lentement. Le taux de consommation est une fonction qui varie avec le temps et est la fonction étant intégrée. Une intégrale définie de cette fonction entre deux valeurs temporelles aa et bb nous donnerait la quantité totale d’énergie consommée pendant cet intervalle de temps :

Calcul de la consommation d’électricité

(Énergie totaleconsommée)=Heure de débutHeure de fin(Taux de consommationinstantané)dt.\begin{matrix} {\begin{pmatrix} {\text{Énergie totale}} \\ {\text{consommée}} \\ \end{pmatrix} = \int_{\text{Heure de début}}^{\text{Heure de fin}}\begin{pmatrix} {\text{Taux de consommation}} \\ {\text{instantané}} \\ \end{pmatrix}\ dt.} \\ \end{matrix}

Bien que ce ne soit pas essentiel à notre discussion ici, mentionnons quels sont les termes et unités physiques appropriés. Retournant à notre analyse dimensionnelle de la Section 11.2, l’énergie est une grandeur dérivée ; c’est le produit d’une force et d’une longueur. Le Chapitre 12 montre que la force est elle-même une grandeur dérivée qui a les unités ML/T2ML/T^{2} et est mesurée dans le système SI en Newtons (N). Ainsi, l’énergie a des unités abstraites ML2/T2ML^{2}/T^{2}, une combinaison d’unités fondamentales qui laisse vraiment perplexe. Dans le système SI, l’énergie est mesurée en Joules (J), et 1J=1Nm=1kgm2/s21\ J = 1\ N\ m = 1\ {kg}\ m^{2}/s^{2}. Le terme physique approprié pour « taux de transfert d’énergie par unité de temps » est la puissance, et l’unité SI de puissance est le watt (W), qui est égal à un joule par seconde.

En notant la consommation totale d’énergie comme EE et le taux de consommation instantanée comme P(t)P(t), nous pouvons réécrire l’ Équation (11.23) comme

Calcul de la consommation d’électricité, cette fois avec une notation plus digne

E=Heure de débutHeure de finP(t)dt.\begin{matrix} {E = \int_{\text{Heure de début}}^{\text{Heure de fin}}P(t)\ dt.} \\ \end{matrix}

Bien que les détails sur la quantification de l’énergie ne soient pas essentiels à notre discussion, il y a une observation très importante à faire : l’ Équation (11.24) est dimensionnellement cohérente. À gauche, la grandeur mesurée est l’énergie, qui dans le système SI est mesurée en joules. Mais à droite, le taux de consommation est mesuré en watts. Comment est-ce possible ? Rappelons que l’intégrale représente une sommation, et les éléments infinitésimaux étant sommés sont le produit de l’intégrande (dans ce cas, P(t)P(t)) et d’un infinitésimal du domaine d’intégration (dans ce cas, dtdt). En termes d’intégrale de Riemann, le premier détermine la hauteur de chaque tranche, et le second détermine sa largeur. Ici, dtdt représente un pas de temps infinitésimalement petit, mesuré en secondes, donc les unités à droite sont W×s=(J/s)×s=J.W \times s = (J/s) \times s = J. Ainsi, les membres gauche et droit de l’ Équation (11.24) sont mesurés en joules.

Nous pouvons étendre cet exemple en calculant la facture d’électricité, plutôt que simplement la consommation totale. Bien sûr, si le prix de l’énergie est fixe, nous multiplions simplement la consommation par le prix. Mais que faire si le prix variait de moment en moment ? (Cela ne devrait pas être difficile à imaginer de nos jours.) Dans ce cas, nous intégrerions le coût plutôt que l’énergie. Nous déterminons comment calculer le coût d’un seul intervalle de durée dtdt (une tranche « différentielle » de temps) et nous sommons ensuite sur tous les intervalles :

Calcul du coût de l’électricité

(Coûttotal)=Heure de débutHeure de finTauxDeDépense(t)dt=Heure de débutHeure de finTauxDeConsommation(t)Prix(t)dt.\begin{matrix} \begin{pmatrix} {\text{Coût}} \\ {\text{total}} \\ \end{pmatrix} & {= \int_{\text{Heure de début}}^{\text{Heure de fin}}{\text{TauxDeDépense}}(t)\ dt} \\ & {= \int_{\text{Heure de début}}^{\text{Heure de fin}}{\text{TauxDeConsommation}}(t)\,{\text{Prix}}(t)\ dt.} \\ \end{matrix}

Passons à un autre exemple d’intégrale. Imaginez un homme utilisant une machine à coudre avec une pédale qui a une réponse à vitesse variable. S’il appuie légèrement sur la pédale, la machine à coudre fait avancer le tissu lentement, et s’il « appuie à fond », la machine à coudre avance à sa vitesse maximale. Maintenant, imaginez sa fille assise sous la table qui regarde son père coudre. Elle ne peut voir que la pédale, mais pas la machine à coudre ni le tissu. Les seules informations disponibles pour la fille sont l’enfoncement de la pédale, et nous supposons que, sur la base de sa connaissance des machines à coudre et des pédales, elle peut inférer une fonction f(t)f(t) qui décrit la vitesse à laquelle le tissu avance au temps tt. La fille observe la pédale pendant une minute environ, puis son père s’arrête et lui demande : « Quelle longueur de tissu ai-je cousue ? » Disons que cette fille est particulièrement brillante et connaît le calcul intégral, donc elle intègre la fonction f(t)f(t) pour obtenir la longueur totale de tissu qui est passée sous l’aiguille. Comme nous le verrons plus tard, ce type de question est en fait très proche des types de problèmes de mécanique résolus avec des intégrales dans les jeux vidéo !

Une dernière analogie utile : pensez à une dérivée comme à un compteur de vitesse qui vous indique un taux de variation instantané, et à l’intégrale comme à un odomètre décrivant la sommation continue de ce taux de variation. Notez que la lecture du compteur de vitesse ne dépend pas de ce voyage en voiture de l’été dernier, ni même de ce qui s’est passé il y a deux secondes. La lecture du compteur de vitesse n’est affectée que par ce qui se passe à cet instant. L’odomètre, en revanche, est un cumul en cours d’exécution, et toute l’histoire depuis que la voiture a quitté la concession est incluse dans sa lecture. Notre fille sous la table de couture doit faire attention à la pédale tout le temps si elle veut faire une estimation précise de la longueur totale de tissu consommée à tout moment donné.

De nombreux types de problèmes d’ingénierie résolus avec des intégrales sont formulés en termes de sommations continues comme celles-ci : Quel est le déplacement total, quand je connais la fonction de vitesse v(t)v(t) ? Quelle est la quantité totale d’eau dans la baignoire, étant donné l’historique de l’angle de déviation du robinet ? Quelle est la quantité restante de carburant, en fonction du taux de consommation comme fonction du temps ? Pour formuler l’intégrale pour des problèmes comme ceux-ci, nous pouvons d’abord imaginer d’approximer la valeur que nous souhaitons calculer en utilisant une somme finie (\sum) et un pas de taille finie (Δx\Delta x). Nous utilisons ensuite un argument de limite pour remplacer la \sum par un \int, et le Δx\Delta x par un dxdx (voir l’ Équation (11.22)). C’est l’essence de ce que l’on entend par « sommation continue ».

Bien sûr, nous pouvons utiliser l’intégrale pour calculer l’aire sous une courbe, comme les manuels de calcul aiment tant à le souligner. Quand nous faisons glisser une ligne de gauche à droite, la fonction intégrée détermine le taux auquel nous accumulons de l’aire. Là où la fonction a une grande valeur, notre aire totale s’accumule plus rapidement, car les « tranches » dans cette zone sont hautes. Cependant, du point de vue d’un programmeur de jeux vidéo, les manuels de calcul semblent se concentrer sur cette application particulière de l’intégrale de manière disproportionnée par rapport à son application aux problèmes du monde réel.

11.7.2La relation entre la dérivée et l’intégrale

Voyons comment calculer des intégrales maintenant que l’utilité d’une intégrale est (nous l’espérons) solidement ancrée dans votre esprit. En regardant la définition de l’ Équation (11.22), on se demande comment diable on peut évaluer cette limite. Pour la dérivée, nous avions pu manipuler l’expression prise à la limite de sorte que nous pouvions simplement substituer Δt=0\Delta t = 0, mais cela ne semble pas possible dans l’ Équation (11.22). Il s’avère que l’ Équation (11.22) est surtout utile comme moyen de reconnaître quand le problème que vous avez est une intégrale, et est utile pour transformer correctement ce problème en notation intégrale. Elle est aussi utilisée quand nous approximons des intégrales numériquement, où au lieu de ramener la largeur des tranches à zéro, nous nous arrêtons à un Δx\Delta x petit mais fini. Mais cette définition n’est pas utilisée pour résoudre des intégrales à la main.

Explorons l’ Équation (11.15), l’intégrale que nous avons pu résoudre par un simple argument géométrique. Puisque c’est une fonction qui décrit la position en fonction du temps, nous devrions pouvoir prendre sa dérivée et obtenir une fonction décrivant la fonction de vitesse v(t)v(t), puis reprendre la dérivée pour obtenir la fonction d’accélération a(t)a(t). Vérifions que c’est bien le cas :

x(t)=x0+v0t+(1/2)at2,x˙(t)=v(t)=v0+at,(en prenant la dérivée)x¨(t)=v˙(t)=a(t)=a.(en reprenant la dérivée)\begin{matrix} {x(t)} & {= x_{0} + v_{0}t + (1/2)at^{2},} & & \\ {\overset{˙}{x}(t) = v(t)} & {= v_{0} + at,} & & \text{(en\ prenant\ la\ dérivée)} \\ {\overset{¨}{x}(t) = \overset{˙}{v}(t) = a(t)} & {= a.} & & \text{(en\ reprenant\ la\ dérivée)} \\ \end{matrix}

C’est conforme aux attentes. Pas de surprises ici, mais il est rassurant de confirmer que les mathématiques et la physique fonctionnent réellement. Nous savions que la dérivée de la fonction de position est la fonction de vitesse. La question est : pourquoi n’avons-nous pas utilisé cette connaissance plus tôt ? Rappelons que nous connaissions v(t)v(t) et essayions de déterminer ce qu’était x(t)x(t). Nous avons pu obtenir la réponse par un argument graphique, mais il semble qu’il pourrait y avoir une autre voie. Au lieu de chercher une fonction pour calculer l’aire sous la courbe, nous aurions pu chercher une fonction de position dont la dérivée était la fonction de vitesse que nous connaissions déjà. Une telle fonction est connue sous le nom d’antidérivée.

Explorons un peu plus cette idée d’« intégration comme antidérivée ». Pour ce faire, nous n’avons essentiellement qu’à appliquer les règles de différenciation, y compris le petit sous-ensemble que nous avons appris dans la Section 11.4.5, à l’envers.36 Supposons que nous commencions avec la fonction de vitesse v(t)=v0+atv(t) = v_{0} + at, et que nous cherchions un x(t)x(t) dont la dérivée est v(t)v(t). Faites semblant pour l’instant de ne pas connaître déjà la réponse. Pour trouver x(t)x(t), nous décomposons v(t)v(t) en ses termes (en utilisant la règle de la somme à l’envers), puis nous prenons l’antidérivée de chaque terme (en utilisant la règle de la puissance à l’envers). Rappelons que la règle de la puissance de différenciation dit essentiellement : « Multiplier par l’exposant, puis diminuer l’exposant de un. » Donc la règle de la puissance pour l’antidifférenciation est : « Augmenter l’exposant de un, puis diviser par le nouvel exposant. » En appliquant ces deux règles à v0+atv_{0} + at, nous écrivons

v(t)=x˙(t)=v0+at,x(t)=v0t+(1/2)at2.\begin{matrix} {v(t) = \overset{˙}{x}(t)} & {= v_{0} + at,} \\ {x(t)} & {= v_{0}t + (1/2)at^{2}.} \\ \end{matrix}

Mais comparez ce résultat à l’ Équation (11.25) ; vous remarquerez qu’il manque un terme x0x_{0}. Que s’est-il passé ? Il y a une certaine « perte d’information » qui se produit quand nous prenons la dérivée. Si nous savons à quelle vitesse nous allions, nous pouvons toujours déterminer combien de distance nous avons parcouru. Cependant, nous ne pouvons pas savoir où nous avons fini à moins de savoir où nous avons commencé. Ce terme supplémentaire x0x_{0} est le « point de départ » que la dérivée élimine, parce que toute valeur constante a une dérivée nulle. Pour cette raison, il n’est pas tout à fait exact de parler de « l’ » antidérivée de v(t)v(t), puisqu’il n’y a pas de fonction unique dont la dérivée est v(t)v(t), mais une infinité. Toutes les différentes antidérivées sont en réalité des copies les unes des autres, décalées verticalement sur le graphe selon leur valeur particulière de x0x_{0}.

Nous avons dit de manière générale qu’il existe une relation entre l’intégrale (définie) et l’antidérivée. Nous savons donc que dans un certain sens, les opérations d’intégration et de différenciation sont des opérations inverses. Le théorème du calcul qui résume ces relations avec précision porte un nom important : le théorème fondamental du calcul. Le théorème est composé de deux parties. (Les sources ne les énumèrent pas toujours dans le même ordre.)

La première partie montre comment une intégrale définie peut être calculée en utilisant une antidérivée.

Théorème Fondamental du Calcul, Partie 1

Soit f(t)=F(t)f(t) = F^{\prime}(t). (Autrement dit, F(t)F(t) est une antidérivée quelconque de f(t)f(t).) Alors l’intégrale définie abf(x)dx\int_{a}^{b}f(x)\ dx peut être calculée comme

abf(t)dx=F(b)F(a).\begin{matrix} {\int_{a}^{b}f(t)\ dx = F(b) - F(a).} \\ \end{matrix}

L’Équation (11.26) peut sembler un peu mystérieuse en termes abstraits, mais quand nous remplaçons le F(t)F(t) et f(t)f(t) génériques par des notations spécifiques à la position et à la vitesse, la première partie du théorème fondamental du calcul semble énoncer l’évidence :

abv(t)dt=x(b)x(a).\int_{a}^{b}v(t)\ dt = x(b) - x(a).

Cela dit que l’effet cumulatif de la vitesse de l’instant aa à l’instant bb (le déplacement net pendant cet intervalle), est égal à la différence entre la position au temps bb et la position au temps aa.

Remarquez que toute antidérivée fonctionnera — peu importe laquelle. C’est parce que le décalage constant x0x_{0} dans x(t)x(t) s’annule lors de la soustraction x(b)x(a)x(b) - x(a). Pour le voir, considérons la métaphore du compteur électrique. Vous pouvez penser à la lecture numérique brute sur le compteur comme à une antidérivée de votre taux de consommation. Les lectures sur le cadran au début et à la fin du mois correspondent respectivement à F(a)F(a) et F(b)F(b). Notez que la valeur numérique brute de la lecture est en grande partie sans intérêt. Elle pourrait contenir des données influencées par quelqu’un qui vivait dans la maison avant vous. La différence entre les deux lectures, cependant, est très pertinente. Elle correspond à l’intégrale définie, et déterminera le montant de votre facture d’électricité pour le mois.

Ou considérons l’odomètre d’une voiture. Disons que vous vouliez mesurer la longueur d’un trajet particulier. Pour ce faire, au début du trajet vous appuieriez sur le bouton de remise à zéro du compteur de trajet dédié que toutes les voitures ont depuis environ 1980, et à la fin du trajet vous liriez simplement la valeur du compteur de trajet. Puis vous vous réjouiriez de n’avoir pas eu à faire le moindre effort mental ni à utiliser un seul principe du calcul. Mais que faire si le compteur de trajet était cassé et que tout ce dont vous disposiez était le compteur total ? Celui-ci ne peut pas être facilement remis à zéro.37 Dans ce cas, armé des connaissances du calcul que vous avez acquises dans ce livre (ou peut-être juste du bon sens que vous auriez pu acquérir n’importe où), vous soustrairiez la lecture du compteur kilométrique à la fin du trajet de la lecture au début du trajet pour obtenir la distance du trajet. Les lectures réelles du compteur sont F(a)F(a) et F(b)F(b), les valeurs de l’antidérivée. Tout comme avec le compteur électrique, les valeurs brutes ne sont pas utiles38 — seule leur différence compte.

La première partie du théorème fondamental du calcul est très importante car c’est ainsi que nous calculons réellement les intégrales, du moins à la main. Rappelons que nous avons défini l’intégrale définie comme une somme d’un grand nombre de tranches dans la limite quand le nombre de tranches tend vers l’infini et les tranches deviennent infinitésimalement minces. Cette définition ne se prête pas à des manipulations algébriques, comme l’était la définition de la dérivée. La première partie du théorème fondamental du calcul dit que bien que nous puissions formuler des problèmes en utilisant la définition de l’intégrale, nous calculons des intégrales définies en trouvant une antidérivée de la fonction intégrée (à la main, du moins).

La deuxième partie du théorème fondamental du calcul est l’envers de la première partie. La première partie disait que les intégrales définies peuvent être calculées en utilisant des antidérivées ; la deuxième partie montre comment définir une antidérivée en termes d’une intégrale définie.

Théorème Fondamental du Calcul, Partie 2

Soit F(t)F(t) définie par

F(t)=t0tf(u)du.\begin{matrix} {F(t) = \int_{t_{0}}^{t}f(u)\ du.} \\ \end{matrix}

Alors la dérivée de F(t)F(t) est donnée par

F(t)=f(t).F^{\prime}(t) = f(t).

Il peut falloir un certain effort pour déchiffrer cette élégance concise, aussi reformulons-la en français courant. Nous commençons avec une fonction donnée ff. Nous formons ensuite une nouvelle fonction FF, dont la valeur est déterminée en prenant l’intégrale définie de ff depuis un point de départ arbitraire t0t_{0}, et un point d’arrivée tt. Notez que l’argument de FF est utilisé pour définir quand arrêter l’intégration de ff. La variable uu est une variable muette de notation pour l’intégration ; elle n’est pas vue en dehors de l’intégrale. La deuxième partie du théorème fondamental du calcul dit que si nous prenons la dérivée de cette nouvelle fonction FF, le résultat est notre fonction originale ff. En ce sens, l’intégration et la différenciation sont des opérations inverses.

Il peut être difficile de saisir la raison pour laquelle tt se retrouve dans ce qui peut sembler un emplacement étrange, définissant la borne supérieure de l’intégration, mais c’est le point essentiel. Le second théorème dit qu’une fonction définie comme une intégrale telle que l’ Équation (11.27) croîtra à un taux déterminé par l’intégrande. Si nous ajustons légèrement la borne supérieure de l’intégration, le changement dans le résultat de la somme globale sera proportionnel à la valeur de l’intégrande. En pensant à une intégrale comme calcul d’une aire, la borne supérieure de l’intégration, tt, détermine la frontière droite. Si nous repoussons cette frontière un peu vers la droite, l’augmentation de la quantité d’aire dépendra de la hauteur de la fonction en tt.

Réécrivons le théorème en utilisant une notation particulière aux déplacements et aux vitesses :

x(t)=t0tv(u)du,x(t)=v(t).\begin{matrix} {x(t)} & {= \int_{t_{0}}^{t}v(u)\ du,} \\ {x^{\prime}(t)} & {= v(t).} \\ \end{matrix}

Nous voyons maintenant que, pour définir le déplacement x(t)x(t) en termes de v(t)v(t), il n’y a vraiment qu’un seul endroit logique où placer tt. La vitesse avant tt est pertinente pour le déplacement qui s’est produit au temps tt, et l’histoire après tt n’est pas pertinente. Nous utilisons tt pour définir le point d’arrêt de l’intervalle de temps des vitesses à intégrer.

D’où vient t0t_{0} ? C’est un point de départ arbitraire, reflétant un degré d’incertitude (ou de liberté) très similaire à la position de départ inconnue (ou non pertinente) x0x_{0}. Nous pouvons choisir t0t_{0} comme bon nous semble, par rapport à ce que nous voulons que nos mesures soient relatives. La valeur de t0t_{0} définit le point où x(t)=0x(t) = 0. Il est probablement plus précis de dire que x(t)x(t) décrit notre position relative. Relative à quoi ? À l’endroit où nous étions au temps t0t_{0}.

Nous sommes maintenant prêts à clarifier la relation parfois confuse entre l’intégrale définie et l’intégrale indéfinie. L’adjectif « définie » dans « intégrale définie » vient du fait que nous avons spécifié les bornes d’intégration. De ce fait, la « réponse » à une intégrale définie peut être un seul nombre. Quand nous évaluons une intégrale définie, telle que

tstarttendv(t)dt,\int_{t_{start}}^{t_{end}}v(t)\ dt,

le tt est « intégré » et n’apparaît pas dans le résultat. La signification de ce qui précède est « la sommation continue de la vitesse pendant l’intervalle de temps tstartt_{start} à tendt_{end}. » Il n’aurait pas de sens que le résultat contienne tt — de quel tt parlerait-on ? Ainsi, si toutes les autres variables de v(t)v(t) sont connues, et que les bornes tstartt_{start} et tendt_{end} sont connues, nous pouvons réduire la réponse à un simple nombre. Si, cependant, v(t)v(t) contient d’autres quantités inconnues (peut-être une densité variable ρ\rho), ou si les bornes d’intégration elles-mêmes sont des paramètres, alors le résultat sera une fonction en termes de ces variables. Dans tous les cas, dans une intégrale définie, le tt ne fera pas partie du résultat. Si vous êtes programmeur, vous pouvez penser au tt comme à une « variable locale » de l’intégrale définie.

Une intégrale indéfinie, en revanche, étant une antidérivée, aura une « réponse » qui est une fonction, pas un seul nombre. Elle est notée simplement en supprimant les bornes d’intégration, comme dans

v(t)dt.\int v(t)\ dt.

Encore une fois, soulignons que bien que cela puisse sembler très similaire à la notation utilisée pour désigner une intégrale définie, sa signification est en réalité assez différente. Le résultat de l’évaluation de cette intégrale ne devrait pas être un nombre, mais une antidérivée de v(t)v(t) ; c’est-à-dire que nous devrions obtenir une fonction de tt. De plus, un résultat correct aura une constante arbitraire ajoutée, connue sous le nom de constante d’intégration, qui nous rappelle qu’il existe toute une famille de fonctions dont la dérivée est v(t)v(t). Ainsi, la signification de l’intégrale indéfinie ci-dessus est « une fonction qui exprime la sommation continue de la vitesse en fonction du temps, depuis un point de départ inconnu ». Nous avons noté ce décalage constant comme x0x_{0}, mais dans un cadre plus général, il est généralement écrit avec un CC majuscule. Par exemple,

Constante d’intégration

v(t)dt=x(t)+x0(notation déplacement et vitesse),f(t)dt=F(t)+C(notation abstraite courante).\begin{matrix} {\int v(t)\ dt} & {= x(t) + x_{0}} & & {\text{(notation\ déplacement\ et\ vitesse)},} \\ {\int f(t)\ dt} & {= F(t) + C} & & \text{(notation\ abstraite\ courante).} \\ \end{matrix}

Nous n’avons pas besoin d’écrire les bornes d’intégration dans une intégrale indéfinie car elles sont implicites. Comme nous l’avons vu dans la deuxième partie du théorème fondamental du calcul, l’interprétation d’une antidérivée en termes d’intégrale définie est d’utiliser l’argument de l’antidérivée comme borne supérieure de l’intervalle d’intégration. En d’autres termes, une intégrale indéfinie est simplement une intégrale définie avec des bornes d’intégration implicites de la forme dans l’ Équation (11.27). Le degré de liberté dans l’ Équation (11.27) connectant l’ensemble des antidérivées possibles était capturé par la borne inférieure d’intégration inconnue (t0t_{0}). Dans une intégrale indéfinie, nous n’écrivons pas les bornes d’intégration, et à la place l’incertitude est contenue dans la constante d’intégration (x0x_{0} ou CC). Nous pouvons résumer cela (en utilisant les deux schémas de nommage) par

L’intégrale indéfinie

v(t)dt=t0tv(u)du=x(t)+x0,x0=x(t0),f(t)dt=t0tv(u)du=F(t)+C,C=F(t0).\begin{matrix} {\int v(t)\ dt = \int_{t_{0}}^{t}v(u)\ du} & {= x(t) + x_{0},} & x_{0} & {= - x(t_{0}),} \\ {\int f(t)\ dt = \int_{t_{0}}^{t}v(u)\ du} & {= F(t) + C,} & C & {= - F(t_{0}).} \\ \end{matrix}

11.7.3Résumé du calcul infinitésimal

Nous avons achevé notre présentation principale du calcul infinitésimal dans ce livre, à l’exception de quelques petits éléments qui apparaîtront dans des sections ultérieures. Notre objectif était d’amener un lecteur sans aucune connaissance du calcul à comprendre la vue d’ensemble de l’utilité des dérivées et des intégrales. Nous avons survolé à grande vitesse les nombreux détails et techniques qui surgissent dans les situations pratiques — ces détails remplissent des milliers de pages dans les manuels de calcul.

Résumons les points importants que vous devez connaître sur le calcul infinitésimal pour utiliser pleinement le reste de ce livre.

11.8Mouvement circulaire uniforme

Assez de calcul — revenons à la physique. Cette section étudie le mouvement d’une particule se déplaçant en cercle à vitesse constante. Nous étudions le mouvement d’une particule car de nombreux calculs physiques peuvent être simplifiés en représentant un corps rigide comme une masse ponctuelle en son soi-disant centre de masse. Puisqu’un trajet circulaire est intrinsèquement limité à un plan, la Section 11.8.1 commence notre étude en deux dimensions. Après avoir établi les relations fondamentales, la Section 11.8.2 montre comment les appliquer dans un monde où le plan orbital est orienté arbitrairement en trois dimensions.

11.8.1Mouvement circulaire uniforme dans le plan

image

Figure 11.14Mouvement circulaire uniforme

Une particule se déplaçant en cercle à vitesse constante n’a pas une vélocité constante ; si c’était le cas, elle se déplacerait en ligne droite. Puisque la vélocité de l’objet varie dans le temps, il doit être soumis à une certaine accélération. Voyons si nous pouvons déterminer laquelle. Considérons un objet se déplaçant à vitesse constante ss sur un trajet circulaire de rayon rr. Pour faciliter nos calculs, et sans perte de généralité, nous établissons un référentiel bidimensionnel qui se trouve dans le plan de mouvement et dont l’origine est au centre du cercle. Rappelons que la vélocité instantanée 𝐯(t)\mathbf{v}(t) d’une particule est toujours tangente à sa trajectoire, donc le vecteur vélocité en tout point donné sera toujours tangent au cercle en ce point. De plus, d’après la définition de la vitesse, nous savons que 𝐯(t)=s{\parallel {\mathbf{v}(t)} \parallel} = s.

Sur le côté gauche de la Figure 11.14, on voit une particule en mouvement circulaire uniforme pendant un pas de temps fini Δt\Delta t. La figure examine l’état de la particule au temps tt et également à un temps ultérieur t+Δtt + \Delta t.

Considérons la vélocité et l’accélération instantanées, en commençant par une approche géométrique. Examinons les triangles sur le côté droit de la Figure 11.14. Le triangle du haut montre la variation de position sur un certain intervalle de temps Δt\Delta t, résultant de la variation angulaire Δθ\Delta\theta. C’est un triangle isocèle dont les côtés ont pour longueur rr, le rayon du cercle, et dont la base est Δ𝐩\Delta\mathbf{p}, qui est la variation nette de position pendant l’intervalle. Le triangle du bas représente la variation de vélocité sur ce même intervalle, et c’est également un triangle isocèle. Les côtés du triangle du bas ont pour longueur ss, puisque nous supposons que la vélocité a une magnitude constante, et la base est Δ𝐯\Delta\mathbf{v}. Les deux triangles sont semblables, car tous deux sont isocèles avec l’angle au sommet Δθ\Delta\theta, donc on peut écrire

Δ𝐯s=Δ𝐩r.\frac{\parallel {\Delta\mathbf{v}} \parallel}{s} = \frac{\parallel {\Delta\mathbf{p}} \parallel}{r}.

En général, la longueur de Δ𝐩\Delta\mathbf{p} mesure un « raccourci » à travers le cercle, plutôt que la distance réellement parcourue le long du périmètre du cercle, qui est rΔθ=sΔtr\Delta\theta = s\Delta t. Mais considérons ce qui se passe lorsque Δt\Delta t et Δθ\Delta\theta deviennent très petits, comme le montre la Figure 11.15.

image

Figure 11.15Une petite rotation.

Remarquons que lorsque Δθ\Delta\theta devient de plus en plus petit, la longueur de Δ𝐩\Delta\mathbf{p} s’approche de la distance réelle, et à la limite, les deux distances sont égales :

limΔt0Δ𝐩=sΔt.\lim\limits_{\Delta t \rightarrow 0}{\parallel {\Delta\mathbf{p}} \parallel} = s\Delta t.

En substituant ce résultat dans notre relation issue des triangles semblables, nous obtenons

Δ𝐯s=Δ𝐩r,limΔt0Δ𝐯s=sΔtr,limΔt0Δ𝐯Δt=s2r.\begin{matrix} \frac{\parallel {\Delta\mathbf{v}} \parallel}{s} & {= \frac{\parallel {\Delta\mathbf{p}} \parallel}{r},} \\ {\lim\limits_{\Delta t \rightarrow 0}\frac{\parallel {\Delta\mathbf{v}} \parallel}{s}} & {= \frac{s\Delta t}{r},} \\ {\lim\limits_{\Delta t \rightarrow 0}\frac{\parallel {\Delta\mathbf{v}} \parallel}{\Delta t}} & {= \frac{s^{2}}{r}.} \\ \end{matrix}

Le membre gauche de l’Équation (11.28) est une variation de vélocité sur un intervalle lorsque la longueur de l’intervalle tend vers zéro. C’est la définition de l’accélération instantanée ! Ainsi, la magnitude de l’accélération est s2/rs^{2}/r.

Bien sûr, l’accélération est une grandeur vectorielle, et tout ce que nous avons déterminé jusqu’ici est sa magnitude (constante). Quelle est sa direction ? Pour le voir, comparons les vecteurs 𝐩(t)\mathbf{p}(t) et Δ𝐯\Delta\mathbf{v} dans la Figure 11.15. Remarquons qu’ils pointent dans des directions opposées. En fait, à la limite lorsque Δθ\Delta\theta tend vers zéro, ils pointent dans des directions exactement opposées. C’est-à-dire que l’accélération est toujours dirigée vers le centre du cercle, c’est pourquoi elle est appelée accélération centripète (« qui cherche le centre »).

Vélocité et accélération du mouvement circulaire uniforme

Lorsqu’un objet se déplace à vitesse constante ss sur un trajet circulaire de rayon rr, la vélocité 𝐯\mathbf{v} est tangente au cercle. L’accélération à tout instant est dirigée vers le centre du cercle et a pour magnitude

a=s2/r.\begin{matrix} {a = s^{2}/r.} \\ \end{matrix}

En combinant quelques éléments de géométrie élémentaire avec quelques idées du calcul, nous avons obtenu les faits les plus importants sur le mouvement circulaire uniforme. Une combinaison légèrement différente de géométrie et de calcul permettra d’obtenir les équations cinématiques réelles. À cette fin, il sera utile de se référer à θ(t)\theta(t), l’angle que fait le vecteur 𝐩\mathbf{p} avec l’axe +x+x en utilisant les conventions mathématiques traditionnelles, comme le montre la Figure 11.16.

image

Figure 11.16 La position de la particule peut être identifiée par l’angle θ\theta

Précédemment, nous nous préoccupions de Δθ\Delta\theta, la variation de cet angle, mais nous considérons maintenant sa valeur en tant que fonction du temps. Nous notons l’angle initial comme θ(0)=θ0\theta(0) = \theta_{0}. Nous définissons également la fréquence angulaire comme ω=s/r\omega = s/r, qui

est mesurée en radians par seconde.39 Ainsi, nous pouvons exprimer l’angle à tout instant comme

θ(t)=θ0+ωt.\theta(t) = \theta_{0} + \omega t.

Angle en fonction du temps

Nous avons déjà vu l’équation paramétrique d’un cercle dans la Section 9.1, nous savons donc comment exprimer les équations cinématiques de la position de la particule en termes du rayon rr et de l’angle θ(t)\theta(t), comme

Position en fonction du temps

x(t)=rcos(θ(t))=rcos(θ0+ωt),y(t)=rsin(θ(t))=rsin(θ0+ωt).\begin{matrix} {x(t)} & {= r\cos(\theta(t)) = r\cos(\theta_{0} + \omega t),} \\ {y(t)} & {= r\sin(\theta(t)) = r\sin(\theta_{0} + \omega t).} \\ \end{matrix}

Puisque la fonction vélocité est la dérivée de la fonction position, nous pouvons différentier ces équations pour obtenir les équations de vélocité. Heureusement, nous avons appris les dérivées des fonctions sinus et cosinus dans la Section 11.4.6 et la règle de la chaîne dans la Section 11.4.7. La différentiation nous donne

Vélocité en fonction du temps

x˙(t)=ddt(rcos(θ0+ωt))=rωsin(θ0+ωt),y˙(t)=ddt(rsin(θ0+ωt))=rωcos(θ0+ωt).\begin{matrix} {\overset{˙}{x}(t)} & {= \frac{d}{dt}\ (r\cos(\theta_{0} + \omega t)) = - r\omega\sin(\theta_{0} + \omega t),} \\ {\overset{˙}{y}(t)} & {= \frac{d}{dt}\ (r\sin(\theta_{0} + \omega t)) = r\omega\cos(\theta_{0} + \omega t).} \\ \end{matrix}

En différentiant encore une fois pour obtenir l’accélération, nous avons

Accélération en fonction du temps

x¨(t)=ddt(rωsin(θ0+ωt))=rω2cos(θ0+ωt),y¨(t)=ddt(rωcos(θ0+ωt))=rω2sin(θ0+ωt).\begin{matrix} {\overset{¨}{x}(t)} & {= \frac{d}{dt}\ ( - r\omega\sin(\theta_{0} + \omega t)) = - r\omega^{2}\cos(\theta_{0} + \omega t),} \\ {\overset{¨}{y}(t)} & {= \frac{d}{dt}\ (r\omega\cos(\theta_{0} + \omega t)) = - r\omega^{2}\sin(\theta_{0} + \omega t).} \\ \end{matrix}

Ces résultats concordent avec nos conclusions précédentes. En comparant les fonctions d’accélération à la position, nous confirmons qu’elles pointent bien dans des directions opposées. De plus, en se rappelant que ω=s/r\omega = s/r, nous notons que, comme prévu, l’accélération a pour magnitude s2/rs^{2}/r.

Parfois, ω\omega est plus directement accessible que ss. Dans ces situations, il est utile de pouvoir exprimer la magnitude de l’accélération centripète uniquement en termes de ω\omega et rr. En résolvant ω=s/r\omega = s/r pour ss, on obtient s=rωs = r\omega. En substituant cela dans l’Équation (11.29), nous obtenons

Accélération en termes de vitesse angulaire 𝛚\mathbf{\omega} et de rayon 𝐫\mathbf{r}

a=s2/r=(rω)2/r=rω2.\begin{matrix} {a = s^{2}/r = (r\omega)^{2}/r = r\omega^{2}.} \\ \end{matrix}

Travaillons sur un exemple intéressant dont les résultats seront utiles dans les sections suivantes. Nous sommes tous à bord d’une centrifugeuse en rotation en ce moment : la Terre ! La rotation de la Terre crée une force centrifuge apparente, qui tend à nous éloigner du centre de la Terre. Heureusement, la gravité terrestre est suffisamment forte pour nous retenir ici. Sachant que le rayon moyen de la Terre est de 6 371 km, quelle est l’accélération centripète ressentie à l’équateur ?

Pour répondre à cette question, nous utilisons l’Équation (11.30). Le rayon est donné comme r=6371kmr = 6~371\ {km}, et le taux de rotation est ω=2π/jour\omega = 2\pi/{jour}.

Accélération centripète à l’équateur due à la rotation de la Terre

a=rω2=(6371km)(2π/jour)2=(6.371×106m)(2π/(86400s))2(6.371×106m)(5.2885×109s2)0.03369m/s2.\begin{matrix} a & {= r\omega^{2} = (6\, 371\ {km})(2\pi/{jour})^{2} = (6.371 \times 10^{6}\ m)(2\pi/(86\, 400\ s))^{2}} \\ & {\approx (6.371 \times 10^{6}\ m)(5.2885 \times 10^{- 9}\ s^{- 2}) \approx 0.03369\ m/s^{2}.} \\ \end{matrix}

Qu’en est-il de la magnitude de l’accélération centripète aux pôles ? Est-elle la même ? Gardez cette question à l’esprit ; nous y reviendrons dans la Section 12.2.1.

11.8.2Mouvement circulaire uniforme en trois dimensions

Jusqu’à présent, nous avons essentiellement travaillé en deux dimensions, en opérant « dans le plan » sans nous préoccuper de la façon dont ce plan pourrait être orienté en trois dimensions. Considérons maintenant le cas plus général. Nous souhaitons décrire la position, la vélocité et l’accélération de la particule sous forme de vecteurs tridimensionnels, où l’axe de rotation (qui est perpendiculaire au plan contenant le trajet circulaire) est orienté arbitrairement.

Supposons qu’une particule en position 𝐩\mathbf{p} se déplace sur un trajet circulaire autour du point 𝐨\mathbf{o}. Comme il existe de nombreux trajets circulaires différents contenant à la fois 𝐨\mathbf{o} et 𝐩\mathbf{p}, nous devons également spécifier un axe de rotation perpendiculaire au plan. Comme nous l’avons fait dans les chapitres précédents (voir la Section 5.1.3 et la Section 8.4), nous décrivons la direction de l’axe en utilisant un vecteur unitaire 𝐧̂\hat{\mathbf{n}}, et, comme précédemment, le signe de 𝐧̂\hat{\mathbf{n}} nous indique quelle direction est considérée comme rotation positive selon la règle de la main gauche. Le scalaire ω\omega définit le taux de rotation, en radians par unité de temps. La question à laquelle nous voulons répondre est la suivante : quelle est la vélocité 𝐯\mathbf{v} de la particule à cet instant ?

image

Figure 11.17Mouvement circulaire uniforme en trois dimensions

Rappelons ce que nous savons déjà. Premièrement, d’après la relation entre vitesse et fréquence angulaire observée précédemment, nous savons que la vitesse s=𝐯s = {\parallel \mathbf{v} \parallel} doit être ωr\omega r, où rr est le rayon du cercle, c’est-à-dire la distance entre 𝐨\mathbf{o} et 𝐩\mathbf{p}. Deuxièmement, 𝐯\mathbf{v} doit être perpendiculaire à 𝐧̂\hat{\mathbf{n}}, sinon la particule s’écarterait du plan contenant le trajet circulaire, et 𝐯\mathbf{v} doit également être tangente à ce trajet. Ainsi, nous connaissons à la fois la magnitude et la direction de la vélocité 𝐯\mathbf{v}, il nous faut juste un moyen de l’exprimer algébriquement. Pour ce faire, introduisons le vecteur 𝐫=𝐩𝐨\mathbf{r} = \mathbf{p} - \mathbf{o}, le vecteur radial de 𝐨\mathbf{o} vers 𝐩\mathbf{p}. Notons que 𝐫\mathbf{r} se trouve dans le plan de rotation et a une longueur constante, le rayon du trajet circulaire, comme le montre la Figure 11.17.

Or, 𝐯\mathbf{v} est perpendiculaire à la fois à 𝐫\mathbf{r} (puisqu’elle est tangente au trajet) et à 𝐧̂\hat{\mathbf{n}} (puisqu’elle se trouve dans le plan orbital). Vous vous souviendrez peut-être que nous avons un outil permettant de calculer un vecteur perpendiculaire à deux autres vecteurs donnés : le produit vectoriel. Peut-être que 𝐧̂×𝐫=𝐯\hat{\mathbf{n}} \times \mathbf{r} = \mathbf{v} ? La direction est correcte,40 mais considérons la longueur. Rappelons d’après la Section 2.12.2 que la longueur du produit vectoriel est égale au produit des magnitudes des entrées, multiplié par le sinus de l’angle entre les deux vecteurs. Or 𝐧̂\hat{\mathbf{n}} est un vecteur unitaire par hypothèse, et 𝐧̂\hat{\mathbf{n}} et 𝐫\mathbf{r} sont perpendiculaires, donc le sinus de l’angle entre eux est égal à un. Ainsi, la longueur du produit vectoriel 𝐧̂×𝐫\hat{\mathbf{n}} \times \mathbf{r} est simplement 𝐫\parallel \mathbf{r} \parallel. La vitesse correcte est ωr=ω𝐫\omega r = \omega{\parallel \mathbf{r} \parallel}, il nous manque donc juste un facteur ω\omega.

En rassemblant tout cela, nous obtenons la formule de la vélocité d’une particule de vecteur radial 𝐫=𝐩𝐨\mathbf{r} = \mathbf{p} - \mathbf{o} tournant autour de l’axe 𝐧̂\hat{\mathbf{n}} à un taux angulaire de ω\omega radians par unité de temps :

Calcul de la vélocité linéaire d’un point à partir de la vélocité angulaire

𝐯=ω𝐧̂×𝐫.\mathbf{v} = \omega\hat{\mathbf{n}} \times \mathbf{r}.

Comme nous l’avons évoqué dans la Section 8.4, la vélocité angulaire est souvent décrite sous forme de carte exponentielle par un seul vecteur ωω=ω𝐧̂{\omega\omega} = \omega\hat{\mathbf{n}} (noter le ωω\omega\omega en gras pour indiquer une grandeur vectorielle). Dans ce cas, la formule est encore plus simple.

Calcul de la vélocité linéaire d’un point à partir de la vélocité angulaire

𝐯=ωω×𝐫.\begin{matrix} {\mathbf{v} = {\omega\omega} \times \mathbf{r}.} \\ \end{matrix}

Considérons maintenant le problème inverse. Supposons que nous connaissions 𝐩\mathbf{p} et 𝐯\mathbf{v}, et que nous souhaitions mesurer la vélocité angulaire par rapport à 𝐨\mathbf{o}. À nouveau, nous pouvons utiliser le produit vectoriel, mais cette fois, nous avons besoin d’une division pour obtenir la bonne longueur :

Vélocité angulaire d’une particule par rapport à un point arbitraire

ωω=𝐫×𝐯𝐫2.\begin{matrix} {{\omega\omega} = \frac{\mathbf{r} \times \mathbf{v}}{{\parallel \mathbf{r} \parallel}^{2}}.} \\ \end{matrix}

Pour comprendre la division par 𝐫2{\parallel \mathbf{r} \parallel}^{2}, considérons deux points sur un disque rigide qui tourne autour de son centre. Supposons que la vélocité angulaire est mesurée par rapport à ce centre. Un point a le vecteur radial 𝐫\mathbf{r}, et un autre a le vecteur radial k𝐫k\mathbf{r}, qui est dans la même direction depuis le centre, mais à une distance mise à l’échelle par le facteur kk. Ces deux points (en fait, tous les points du disque) devraient avoir la même vélocité angulaire. Ainsi, une division par r\parallel r \parallel est nécessaire pour compenser la variation de 𝐫\mathbf{r} lorsqu’on ajuste le rayon. La division supplémentaire est nécessaire car les points extérieurs ont une vélocité plus élevée ; si on se déplace sur le disque en mettant 𝐫\mathbf{r} à l’échelle par kk, le nouveau point aura une vélocité également mise à l’échelle par kk.

Bien que jusqu’à présent nous ayons supposé que 𝐩\mathbf{p} tourne effectivement autour de 𝐨\mathbf{o}, ce n’est pas forcément le cas. Il pourrait tourner autour d’un autre point, ou se déplacer en ligne droite. Cependant, nous pouvons quand même calculer la vélocité angulaire de 𝐩\mathbf{p} par rapport à 𝐨\mathbf{o}. En essence, ce que l’Équation (11.32) nous dit, c’est quelle serait la vélocité angulaire si 𝐩\mathbf{p} orbitait effectivement autour de 𝐨\mathbf{o}, dans le plan contenant à la fois 𝐫\mathbf{r} et 𝐯\mathbf{v}. L’axe de rotation, qui est parallèle à ωω\omega\omega, est perpendiculaire à ce plan. Il y a en fait une légère complication — 𝐫\mathbf{r} et 𝐯\mathbf{v} pourraient ne pas être perpendiculaires, ce qui serait le cas si la particule orbitait autour de 𝐨\mathbf{o}. Le produit vectoriel dans l’Équation (11.32) rejette essentiellement toute vélocité parallèle à 𝐫\mathbf{r} ; seule la vélocité perpendiculaire à 𝐫\mathbf{r} contribue aux résultats.

Si la particule 𝐩\mathbf{p} orbite effectivement autour de 𝐨\mathbf{o} à vitesse constante, alors la vélocité angulaire calculée par l’Équation (11.32) sera constante. En général, cependant, la vélocité angulaire mesurée par rapport à un point arbitraire n’est pas constante. Par exemple, considérons une particule se déplaçant avec une vélocité linéaire constante. La vélocité angulaire mesurée par rapport à un point stationnaire 𝐨\mathbf{o} croîtra à mesure que la particule s’approche de 𝐨\mathbf{o}, atteindra un maximum au point de plus proche approche, puis diminuera. De plus, même si la particule se déplace bien sur un trajet orbital, la vélocité angulaire ne sera constante que lorsqu’elle est mesurée par rapport au centre de l’orbite.

Un exemple extrêmement important de mouvement orbital en 3D est celui d’une particule attachée à un corps rigide tournant autour d’un axe. Choisissons 𝐨\mathbf{o} à l’intersection de l’axe de rotation et du plan contenant l’orbite circulaire de 𝐩\mathbf{p} ; cela rend 𝐫\mathbf{r} perpendiculaire à l’axe de rotation. Dans ces conditions, la vélocité angulaire orbitale calculée par l’Équation (11.32) est la même pour chaque particule, et c’est également la même que la vélocité angulaire propre du corps rigide. Nous en dirons davantage dans le prochain chapitre.

Nous n’avons pas souvent besoin de calculer la vélocité angulaire d’un point par rapport à un point qui n’est pas le centre de l’orbite. (Cependant, l’Équation (11.31) est fréquemment utilisée pour calculer une vélocité linéaire de point basée sur sa vélocité orbitale.) Pourquoi en parlons-nous alors ? Parce que le calcul est similaire à la façon dont on mesure le couple (voir la Section 12.5) pour une force appliquée dans une direction arbitraire à un emplacement arbitraire.

Exercices

  1. Le Pascal est une unité de mesure de la pression, définie comme un Newton par mètre carré. Un Pascal est égal à combien de psi ? (Le psi correspond à une livre de force par pouce carré.)

  2. La position 1D d’une particule est décrite de façon continue par morceaux

    x(t)={2tt20t<2,02t<4,sin(πt)4t<7,7t7t.\begin{matrix} {x(t) = \left\{ \begin{matrix} {2t - t^{2}} & {0 \leq t < 2,} \\ 0 & {2 \leq t < 4,} \\ {\sin(\pi t)} & {4 \leq t < 7,} \\ {7 - t} & {7 \leq t.} \\ \end{matrix} \right.} \\ \end{matrix}

    Tracez le graphe du mouvement de la particule.

  3. Quelle est la vélocité moyenne de la particule de l’Exercice 2 sur les intervalles suivants ?

    1. (a) t=01t = 0\ldots 1 ?

    2. (b) t=12t = 1\ldots 2 ?

    3. (c) t=02t = 0\ldots 2 ?

    4. (d) t=5.56.5t = 5.5\ldots 6.5 ?

    5. (e) t=09t = 0\ldots 9 ?

  4. Écrivez une fonction continue par morceaux similaire v(t)v(t) décrivant la vélocité de la particule de l’Exercice 2 au temps tt. Dans ce cas, la vélocité n’est pas définie aux « jonctions » entre les morceaux, donc ne vous préoccupez que de ce qui se passe au milieu de chaque morceau. (C’est malheureusement l’un de ces points plus délicats que nous avons dû sauter.)

  5. Quelle est la vélocité instantanée de la particule de l’Exercice 2 aux instants suivants ?

    1. (a) t=0.1t = 0.1

    2. (b) t=1.0t = 1.0

    3. (c) t=1.9t = 1.9

    4. (d) t=4.1t = 4.1

    5. (e) t=5t = 5

    6. (f) t=6.5t = 6.5

    7. (g) t=8t = 8

    8. (h) t=9t = 9

  6. Écrivez une fonction continue par morceaux similaire a(t)a(t) décrivant l’accélération de la particule de l’Exercice 2 au temps tt. À nouveau, ne vous préoccupez pas de ce qui se passe aux points de jonction.

  7. Quelle est l’accélération de la particule aux instants suivants ?

    1. (a) t=0.1t = 0.1

    2. (b) t=1.0t = 1.0

    3. (c) t=1.9t = 1.9

    4. (d) t=4.1t = 4.1

    5. (e) t=5t = 5

    6. (f) t=6.5t = 6.5

    7. (g) t=8t = 8

    8. (h) t=9t = 9

  8. Quelle situation physique est signifiée par un discriminant négatif dans l’Équation (11.16), donnant des solutions complexes ? Et si le discriminant est nul et qu’il n’y a qu’une seule solution ?

  9. Un projectile est lancé avec une vitesse initiale de 150 ft/s, avec un angle d’inclinaison de 40° depuis la position initiale 𝐩0=(0ft,10ft)\mathbf{p}_{0} = (0\ {ft},{10\ ft}).

    1. (a)Quelle est la vélocité initiale sous forme vectorielle ?

    2. (b)À quel moment le projectile atteindra-t-il son apogée ?

    3. (c)Quelles sont les coordonnées du projectile à l’apogée ?

    4. (d)Combien de temps faudra-t-il au projectile pour revenir à une altitude de y=10y = 10 ?

    5. (e)Quel sera le déplacement horizontal à cet instant ?

  10. À la fin de notre discussion sur les projectiles dans la Section 11.6, nous avons posé le problème de résoudre pour le temps d’intersection lorsque l’accélération est un vecteur arbitraire 𝐚\mathbf{a}. Prenez l’Équation (11.18) et faites le produit scalaire des deux membres avec 𝐚\mathbf{a}, puis résolvez pour tt. (Utilisez la formule quadratique, comme précédemment.)

  11. Les exponentielles complexes telles que eixe^{ix} (où ii est le nombre imaginaire tel que i2=1i^{2} = - 1) sont très importantes en équations différentielles, systèmes de contrôle et traitement du signal. Bien qu’il semble étrange de mettre un nombre complexe dans l’exposant, la formule d’Euler donne une interprétation significative. Pour trouver cette interprétation, vous pourriez aller sur wikipedia.com et la chercher (ce qui explique pourquoi la réponse se trouve dans la section des solutions en fin de livre, de toute façon). Mais avant de le faire, développez la série de Taylor de eixe^{ix} et voyez si vous pouvez le découvrir par vous-même. (Ensuite allez en ligne lire sur l’importance surprenante de cette expression.)

  12. La Station spatiale internationale orbite la Terre à environ 340 km au-dessus de la surface terrestre. (L’orbite est en réalité elliptique, mais ignorez cela et supposez qu’elle se déplace en mouvement circulaire uniforme.) Sachant que la vitesse moyenne est d’environ 27 740 km/h, quelle est la période orbitale ? Quelle est l’accélération centripète en m/s2{m/s}^{2} ? Réfléchissez bien !

Bien sûr, les mathématiciens savent comment écrire tous ces nombres. Vous pourrez lire un jour dans un livre de mathématiques comment les écrire tous sous une forme élaborée et élégante, mais il est d’abord utile de savoir de façon approximative de quoi vous essayez de parler.

— Richard Feynman (1918–1988) extrait de
The Feynman Lectures on Physics

Références

[1] Robert Resnick et David Halliday.   Physics, troisième édition.   New York : John Wiley and Sons, 1977.

  1. Grâce à Heisenberg, nous savons que ce n’est pas possible au niveau atomique, mais quand nous disions « petit », nous ne voulions pas dire aussi petit.

  2. Une unité astronomique est égale à la distance moyenne entre la Terre et le Soleil, soit environ 150 millions de kilomètres ou 93 millions de miles. C’est un grand nombre, mais je ne dirais pas qu’il est astronomique.

  3. En fait, c’est encore plus facile maintenant. Seuls les dinosaures comme les auteurs savent lire une horloge analogique à aiguilles.

  4. D’accord, celle-là n’est peut-être pas si familière. Elle l’est pour l’un des auteurs, mais c’est parce qu’il est britannique.

  5. Par « chercher », nous entendons utiliser Internet. Il n’y a pas de place dans ce livre pour des tables d’informations facilement trouvables en ligne. Nous avions besoin de l’espace pour toutes nos opinions, blagues et notes de bas de page inutiles.

  6. Ian Parberry est partagé à ce sujet. Bien qu’il soit britannique et pense qu’il devrait donc soutenir Newton, le directeur de thèse de son directeur de thèse... remontant 14 générations en arrière était Leibniz, et il se sent donc une certaine loyauté « familiale » envers lui.

  7. Pascal et Descartes sont des « cousins » dans la lignée de directeurs de thèse d’Ian Parberry, remontant à la 16e génération, mais il ne peut s’empêcher de penser à la Chanson du Philosophe des Monty Python chaque fois qu’il pense à Descartes.

  8. Surtout parce que c’est la chose la plus facile à créer pour des auteurs paresseux dans Photoshop.

  9. La vitesse choisie pour le lièvre ressemble un peu à la réalité, mais pour des raisons pédagogiques et pour que la Figure 11.1 tienne bien sur une page, la vitesse de la tortue est totalement inventée. Les puristes doivent se rappeler que c’est une histoire avec un lapin et une tortue qui parlent. Oh, et un furlong vaut 1/8 de mile.

  10. L’abréviation « fur » signifie « furlongs » et n’a rien à voir avec la fourrure du lapin.

  11. Bien qu’il puisse sembler soigneusement arrangé que le mouvement du lièvre soit décrit par une équation quadratique à coefficients entiers, nous verrons plus tard que ce n’est pas si arrangé qu’on pourrait le croire. La nature semble apprécier les équations quadratiques. Mais vous avez raison sur les coefficients entiers, qui ont été soigneusement choisis.

  12. Aucun jeu de mots n’était intentionnel. Malheureusement.

  13. Incluant des notions telles que la continuité, les limites à gauche, à droite, etc.

  14. Les auteurs définissent « avoir eu » comme : vous l’avez passé, vous l’avez compris et vous vous en souvenez.

  15. Spoiler : nous vous l’avons déjà donné dans cette section !

  16. Cela serait vrai si nous opérions dans une plage où la loi des gaz parfaits est une approximation valide, et que la température reste constante.

  17. Pour ainsi dire.

  18. Le lecteur attentif aura peut-être remarqué que, comme il n’y a qu’un nombre fini de grimpeurs à échantillonner, S(y)S(y) provient de données d’échantillons et présente donc des « marches » discontinues plutôt qu’une ligne lisse. À ces marches, la dérivée n’est pas vraiment définie. Pour les besoins de cette illustration, supposons que nous avons ajusté une courbe lisse aux données expérimentales afin d’obtenir une fonction continue.

  19. Nous venons de faire des guillemets dans l’air avec nos doigts en écrivant ce mot.

  20. Nous utilisons le mot « officiel » ici parce qu’il existe d’autres façons de définir la dérivée qui conduisent à des méthodes améliorées pour approximer numériquement les dérivées avec un ordinateur. Ces méthodes sont utiles lorsqu’une solution analytique est trop difficile ou trop lente à calculer.

  21. Pour tous les puristes mathématiques qui s’offusquent des lignes verticales aux discontinuités où la dérivée est mathématiquement indéfinie, la justification de l’ingénieur dans cette situation et d’autres similaires est que la formule mathématique n’est qu’un modèle de ce qui est en réalité une situation physique.

  22. Remarquez notre maladroit besoin de parenthèses avec (Δt)2(\Delta t)^{2} pour éviter la notation potentiellement confuse Δt2\Delta t^{2}.

  23. Supposez que les rails sont droits, de sorte que les axes du corps du train sont alignés avec les axes du monde, et qu’aucune rotation n’est nécessaire.

  24. Attention, t0t^{0} est indéfini quand t=0t = 0.

  25. Dans ce cas particulier, la série de Taylor reçoit le nom plus spécifique de série de Maclaurin.

  26. Nous soulignons que nous n’avons pas prouvé que ce sont les bonnes dérivées, parce que nous avons commencé avec le développement en série de Taylor, qui est en réalité défini en termes de dérivées.

  27. Le bâtiment que tout le monde appelle encore la Tour Sears.

  28. Eh bien, presque. Il y a certaines limitations que nous devons imposer à v(t)v(t). Par exemple, si elle explose et tend vers l’infini, il est probable, bien que pas certain, que le déplacement sera infini ou indéfini. Dans ce livre, nous nous concentrons sur les phénomènes physiques et nous évitons ces problèmes en supposant que nos fonctions seront bien définies.

  29. Cela revient aussi fréquemment lors des entretiens d’embauche.

  30. La racine négative nous donne l’autre point où la parabole infinie contenant la trajectoire du roulement à billes croise le trottoir.

  31. Nous utilisons 𝐩\mathbf{p}, abréviation de « position », plutôt que 𝐱\mathbf{x}, pour éviter l’hypothèse que la coordonnée xx est spéciale par rapport à yy ou zz.

  32. À une exception près — voir l’Exercice 8.

  33. Tout cela en supposant un mouvement de projectile idéal, qui ignore la résistance de l’air. Bien sûr, nous ne pouvons pas percuter un bâtiment adjacent, sinon le mouvement horizontal serait certainement pertinent.

  34. Un salut à tous les fans des Eagles qui ont compris cette blague. Les autres le découvriront à la longue. Prenez-le facilement, surmontez-le, et vous aurez un sentiment paisible et facile.

  35. En fait, « summierung » puisqu’il parlait allemand ; nous avons juste la chance que ça fonctionne aussi en anglais.

  36. Cette affirmation s’applique plus généralement qu’on ne le reconnaît habituellement. Par exemple, si vous avez eu quelques cours de calcul, remarquez que la technique d’intégration connue sous le nom d’« intégration par parties » n’est en réalité que la règle du produit de la différentiation à l’envers !

  37. Ni, comme nous l’avons appris de La Folle Journée de Ferris Bueller, facilement réinitialisé.

  38. Du moins pas à cette fin. Quand la chaîne de distribution casse à 5 miles après l’expiration de votre garantie, ces valeurs brutes sont très importantes.

  39. Oméga (ω\omega) est la lettre traditionnelle pour la fréquence angulaire. Pour voir d’où vient le calcul s/rs/r, considérez que la circonférence du cercle est 2πr2\pi r, et cette distance est parcourue à la vitesse ss. Par conséquent, la fréquence angulaire est 2πr/s2\pi r/s tours par seconde. Mais un tour est égal à 2π2\pi radians, donc le facteur 2π2\pi s’annule. C’est un exemple de la raison pour laquelle l’utilisation des radians est souvent si commode (à condition que nous travaillions de façon symbolique et que nous ne nous préoccupions pas des valeurs numériques des angles).

  40. Ne faites pas seulement confiance à la Figure 11.17, utilisez votre main gauche pour vérifier cela. Votre pouce est le premier argument, 𝐧̂\hat{\mathbf{n}} ; l’index est le deuxième argument, 𝐫\mathbf{r} ; et le majeur est le résultat, 𝐯\mathbf{v}. Regardez votre main avec le pouce pointé vers vous (c’est l’axe de rotation, souvenez-vous), et faites tourner votre main dans la direction de 𝐯\mathbf{v} (votre majeur), en effectuant la rotation autour de votre pouce. (La particule est à l’extrémité de votre index.) Votre main tournera dans le sens des aiguilles d’une montre de votre point de vue, ce qui est la définition d’une rotation positive selon la règle de la main gauche.

<< Sujets mathématiques issus des graphiques 3D

Retour en haut

Mécanique 2 : Dynamique linéaire et rotationnelle >>