Spirale logarithmique |
|||
|
|||
Pour charger le fichier Excel
|
|||
Si on construit une succession de triangles semblables par similitude d'un triangle initial on obtient une courbe appelée spirale logarithmique. Nous allons donc nous poser quelques questions autour de cet objet qui fascina tout d'abord Descartes puis plus particulièrement Johann Bernoulli qui en fit même graver une sur sa tombe accompagée des mots « Eadem mutata resurgo », citation que vous pouvez facilement traduire. Choisissons donc une similitude directe dans le plan complexe, de centre O :
qui n'est jamais que
la traduction complexe des suites géométriques. On peut changer le centre
évidemment, mais ça ne change rien au problème, la situation sera alors
déplacée au point fixe de la similitude et les résultats seront identiques. On
peut également compliquer en prenant des similitudes indirectes, mais on risque
de perdre la rotation autour de O dans le cas de la symétrie glissée. |
|||
, |
|||
soit |
|||
|
|||
Sur la fig. 1 nous avons |
|||
, |
|||
soit un rapport de similitude de 0,99 et un angle de rotation de 0,0055. On remarque immédiatement que les arguments des affixes zn des points Mn forment une suite arithmétique alors que leurs modules forment une suite géométrique. |
|||
Figure 1 : spirale définie par itération de similitudes. |
|||
Equation polaire |
|||
Prenons par exemple un morceau de spirale commençant au point A d'affixe 1 et finissant à un point B d'affixe |
|||
; |
|||
les deux paramètres k et caractérisent entièrement la spirale par la similitude z'=az. Divisons l'angle en n parties identiques, d'angle ; pour arriver à B avec OB=k, il faut effectuer n similitudes de rapport et d'angle correspondant à autant de triangles dont les côtés ont pour longueurs et .Le côté sera alors porteur de la tangente à la spirale au point ; nous allons donc chercher la direction de ce vecteur, puis en faisant tendre n vers l'infini obtenir une relation importante caractérisant la spirale. |
|||
|
|
|
|
Le vecteur en question a pour affixe |
|||
; |
|||
sa direction par rapport à est donnée par |
|||
|
|||
dont nous allons chercher la cotangente. On a |
|||
. |
|||
Le premier membre est indéterminé, mais on trouve les limites avec un D.L. ou en utilisant les limites suivantes en 0 : |
|||
et ; |
|||
le deuxième membre tend vers 0 et devient nul. |
|||
|
|||
donc tend vers une valeur fixe b, ce qui signifie que tend vers un angle fixe égal à |
|||
. |
|||
Passons en coordonnées polaires et cherchons l'équation polaire de la spirale : nous laissons au lecteur la joie de démontrer que la tangente à une courbe est telle que |
|||
|
|||
où est l'angle entre le rayon vecteur et la tangente. ; pour notre spirale le résultat est donc que , simple équation différentielle du premier ordre qui nous donne |
|||
|
|||
ou plutôt |
|||
|
|||
avec les notations habituelles. Dans notre cas où f(0)=1 (point A) nous avons a=1. |
|||
En coordonnées paramétriques la courbe est alors |
|||
|
|||
ou encore avec les complexes : |
|||
. |
|||
Comme on a |
|||
|
|||
et |
|||
|
|||
et on retrouve bien nos similitudes ; ces relations permettent d'obtenir la plupart des caractéristiques de la spirale comme présentées dans |
|||
http://www.mathcurve.com/courbes2d/logarithmic/logarithmic.shtml |
|||
par exemple.
|
|||
Calcul de la longueur de la spirale |
|||
Toujours avec les mêmes paramètres, ici encore avec |
|||
|
|||
nous cherchons la longueur de spirale comprise entre A et B : on peut évidemment repartir sur les suites géométriques en calculant la somme des longueurs , ce qui aboutit au même genre de calculs que précédemment, mais on peut également utiliser les formules : pour une courbe en paramétriques la longueur de courbe est donnée par |
|||
, |
|||
ce qui donne en polaires |
|||
; |
|||
nous calculons donc ici |
|||
. |
|||
Remarquons qu'avec k=1 et la spirale est un cercle de rayon a auquel cas b=0, |
|||
, et lorsque b tend vers 0 (ouf.).
|
|||
Aire de la spirale |
|||
Cherchons l'aire de notre spirale préférée pour t compris entre 0 et : là encore on peut réutiliser les triangles puisqu'on connaît tous les angles de chaque triangle ainsi que les côtés. L'aire de sera |
|||
, |
|||
il reste à intégrer, soit |
|||
; |
|||
comme dans le cas précédent si k=1 et on a un cercle de rayon a et on retrouve bien l'aire du cercle .
|
|||
Une application |
|||
Johann Bernoulli pour montrer les possibilités du « nouveau calcul » change la loi de Newton en |
|||
|
|||
pour une force en |
|||
|
|||
et résout le problème de l'attraction des corps pour une telle loi. |
|||
Reprenons la technique déjà vue pour résoudre Newton à l'aide du produit vectoriel et du produit scalaire : nous avons dit que la trajectoire se situait dans un plan que nous assimilons au plan complexe en prenant un repère tel que et ; l'affixe de M est alors . Dérivons z deux fois : |
|||
|
|||
puis |
|||
|
|||
qui doit être dans la même direction que z ; il faut donc que |
|||
; |
|||
si on multiplie cette égalité par r, on fait apparaître la dérivée de : |
|||
|
|||
d'où (1) |
|||
, |
|||
c constante. En fait nous avons réutilisé la même technique que dans le livre mais en passant par les complexes. |
|||
Avec notre force centrale en |
|||
|
|||
nous obtenons donc |
|||
, |
|||
k=constante, soit en introduisant (1) |
|||
|
|||
nous procédons comme dans le livre en posant |
|||
, |
|||
soit |
|||
|
|||
et |
|||
|
|||
et |
|||
. |
|||
Il est quand même amusant de retrouver l'équation de l'oscillateur harmonique dans cette situation ! Le plus gros du travail est maintenant fait ; supposons que |
|||
|
|||
(négatif), nous obtenons |
|||
, |
|||
soit |
|||
|
|||
d'où le type de courbes suivantes : |
|||
fig. 2 : spirale |
|||
si A ou B est nul on retrouve nos spirales logarithmiques simples. |
|||
Maintenant si |
|||
|
|||
(positif), on a |
|||
, |
|||
soit en réels d'où |
|||
|
|||
dont voici un exemple : |
|||
fig. 3 : spirale hyperbolique |
|||
On peut se dire que ce type de calcul ne présente qu'un intérêt limité en tout cas en termes physiques, mais on peut soulever quelques questions de fond à ce propos : dans la loi de Newton la gravitation impose la force en 1/r2 obtenue de manière empirique (on peut parler de potentiel en 1/r, mais ça ne fait que déplacer le problème) ; de la même manière en électromagnétisme on obtient également une force en 1/r2 pour des particules chargées, force dont le vecteur est le photon ; l'explication couramment admise est que une particule émet des photons dans toutes les directions et qu'à une distance r le flux de photons reçus est inversement proportionnel à la dispersion dans l'espace de ces photons. Pour une quantité K de photons émis, à la distance r on retrouvera |
|||
|
|||
photons par unité de surface ce qui « explique » la force en question. |
|||
Dans le cas de la gravitation nous devons faire intervenir une variable supplémentaire : le temps (qui n'intervient pas pour l'électromagnétisme car le temps propre d'un photon est nul et donc fixé une fois pour toutes). Supposons qu'une particule, le graviton, serve de vecteur à la force gravitationnelle et qu'une masse quelconque émette K gravitons dans toutes les directions de l'espace-temps à 4 dimensions ; à la distance r on retrouvera |
|||
|
|||
gravitons (inversement proportionnel au volume d'une sphère), soit une force en 1/r3. |
|||
Si cette idée est correcte pour quelle raison la loi de Newton serait-elle en 1/r2 ? A priori le temps n'est pas trop variable dans notre système solaire et un peu comme dans le cas du photon ne doit pas intervenir réellement ; on fixe donc cette variable ce qui nous ramène dans un espace classique à 3 dimensions et on récupère notre loi en 1/r2, les coniques et tout le tremblement. La stabilité du temps à l'intérieur du système solaire est due certainement à la masse du soleil et on doit voir apparaître des variations temporelles lorsqu'on s'en éloigne : les sondes Voyager qui sont en train de quitter le système solaire subissent actuellement une accélération imprévue qui pourrait s'expliquer par ce biais. Une autre situation non-expliquée est celle de la formation des galaxies spirales (comme la nôtre.) : si nous envisageons une force en 1/r3 par rapport au centre de la galaxie le problème se résout de lui-même avec un temps « instable » ; chaque système solaire, chaque étoile posséde alors un temps propre, la galaxie dans son ensemble fonctionne dans un intervalle de temps suffisamment large pour autoriser ce type de configuration. Là encore ce ne sont que des idées personnelles mais il est toujours intéressant d'en émettre. Un travail qui ne doit pas être très amusant serait de refaire le travail précédent dans les quaternions.
|
|||
Pour d'autres exemples voir |
|||
http://irem.u-strasbg.fr/irem/ouvert/Ouvert97/Les spirales.pdf |
|||
|
|||
|
|||
|
|||
|