M1 Cour Ylm
M1 Cour Ylm
M1 Cour Ylm
gradient
divergence
rotationnel
1
Opérateurs classiques en coordonnées sphériques
Laplacien
2
Résolution de l’équation de Laplace (1/5)
Un champ quelconque sur une sphère doit satisfaire l'équation de Laplace loin des
sources (∆P = 0). donc si on cherche une base sur laquelle exprimer ce champ P,
les "fonctions-vecteurs" de cette base doivent aussi satisfaire cette équation. On
doit donc chercher les fonctions V(r,θ,ϕ) qui satisfont ∆V = 0.
Soit :
l'équation devient :
Solution du terme en r
4
Résolution de l’équation de Laplace (3/5)
Solution du terme en ϕ
5
Résolution de l’équation de Laplace (4/5)
Solution du terme en θ
On obtient :
6
Résolution de l’équation de Laplace (5/5)
Solution du terme en θ
Les fonctions solutions de cette équation sont les polynômes de Legendre, qui
n'existent que pour des valeurs entières et positive de l et des valeurs entières
de m entre -l et +l. Ces polynômes s'écrivent :
En fait ces fonctions sont moins compliquées qu'il n'y parait puisque ce sont
tout simplement des produits de sinus et cosinus de puissances l et m.
Il est assez facile de se rendre compte que chaque Polynôme de Legendre passe
un certain nombre de fois par zero quand θ varie entre 0 et π.
Le nombre de passage par zero d'un polynôme de degré l et d'ordre m est
exactement l-m.
7
Polynômes de Legendre
P 00 (cos(θ)) = 1
P 10 (cos(θ)) = cosθ
P 11 (cos(θ)) = -sinθ
P 20 (cos(θ)) = ½ (3cos2θ−1)
P 2 (cos(θ)) = 3 sin2θ
2
8
Harmoniques Sphériques
Pour un champ que l'on peut définir sur une surface sphérique (r=constante), par
exemple la topographie terrestre, on supprime la dépendance en r pour obtenir les
harmoniques sphériques surfaciques :
9
Harmoniques Sphériques
m=0 m=1
L=1
10
Décomposition d’un champ sur la base des harmoniques sphériques
+ 3
+
4
11
Décomposition d’un champ sur la base des harmoniques sphériques
12
Décomposition d’un champ sur la base des harmoniques sphériques
Tout champ qui a un Laplacien nul peut donc "s'écrire" sur la base des
harmoniques sphériques :
En fait, Clm est un nombre complexe tel que : Clm = (alm +i blm)
et l'exponentielle complexe eimϕ s'exprime aussi : eimϕ = cos(mϕ) + i sin(mϕ)
On obtient donc :
13
Décomposition d’un champ sur la base des harmoniques sphériques
P étant un champ réel, il est clair que le terme imaginaire de la recomposition doit
toujours être nul. Cela impose que la somme sur l et m des termes imaginaires soit
nulle. La somme pour m allant de -l à l peut être réécrite comme une somme pour
m de -l à 1, 0, et 1 à l. On obtient alors :
Il est facile de voire que pour les valeurs positive et négatives de m, les termes
imaginaires en i(almsin(m ϕ) + blm cos(m ϕ)) vont s'annuler 2 à 2 à condition que :
al-m = alm et bl-m = -blm ; et que pour les terme m = 0 (al0+ i.bl0) il faut tout
simplement que bl0 = 0. On obtient alors :
Alm et Blm sont les coefficients de ce champ, tout comme une courbe à une
dimension peut être représentées par la suite infinie de ses coefficients dans une
décomposition en série de Fourier de Fourier. l est le degré, m est l'ordre.
14
Correspondance degré - longueur d’onde
λ = 2 π RT / l
15
Spectre d’un champ
On a deux champs (A et B) définis en tous points sur la sphère (θ et ϕ). La décomposition sur la base des
harmoniques sphériques de ces deux champs est constituées des deux suites de nombres complexes Alm et Blm.
le coefficient de corrélation entre les deux champs au degré l est :
Si C(l)=1 alors les coefficients des champs A et B au degré l sont proportionnels. Pour que les deux champs
A et B soient complètement identiques, il faut que tous les coefficients de corrélation valent 1 et que les
coefficients de proportionnalité de tous les degrés valent 1 aussi.
C(l)=1 (ou = -1) implique qu'il y a 100% de chance que les deux champs soient parfaitement corrélés (ou
anti-corrélés) au degré l considéré.
Plus précisement : au degré l il y a N=2l+1 termes (m varie de -l à l). Le coefficient de corrélation est r =
C(l). Or, dans un problème à N degrés de liberté, la probabilité pour que la variable t = (N/1-r2) soit plus
petite qu'une certaine valeur donnée t0, est la student's t-distribution: Q(t,N). La valeur 1-Q(t,N) est donc le
niveau de confiance auquel l'hypothèse d'une corrélation due au hasard est infirmée.
Par exemple : r=1 => t = (qqsoit N) => Q(t,N) = 0 => 1 - Q(t,N) = 1 = 100%
C(l)=1 => il y a 100% de chance que la corrélation entre les champs A et B au degré l ne soit pas due au
17 hasard
Courbe de corrélation
On pourrait en déduire que la source à l'origine des deux champs est distincte à grande
longueur d'onde (ce n'est pas le même phénomène qui est à l'origine des deux champs), et
peut être commune à courte longueur d'onde...
18
Application à la topographie terrestre
19
Application à la topographie terrestre
20
Application à la topographie terrestre
21
Application à la topographie terrestre
22
D’autres planètes….
23