Annette EickerAPMG 1 1 19.01.2014 Annette Eicker Kugelfunktionen.
-
Upload
eleonore-schneidman -
Category
Documents
-
view
134 -
download
9
Transcript of Annette EickerAPMG 1 1 19.01.2014 Annette Eicker Kugelfunktionen.
Annette Eicker APMG 1
1
11.04.23
Annette Eicker
Kugelfunktionen
Annette Eicker APMG 1
2
11.04.23
Wiederholung: Gravitationspotential
PotentialPotential
)()(1
QP
P dGV rrrr
r
z
y
x
PrAufpunkt:
w
v
u
QrQuellpunkt:
FeldstärkeFeldstärke
dwdvduwvuwzvyux
wzG
dwdvduwvuwzvyux
vyG
dwdvduwvuwzvyux
uxG
zzyxV
yzyxV
xzyxV
zyx
),,()()()(
),,()()()(
),,()()()(
,,
,,
,,
,,
3222
3222
3222
g
dwdvduwvuwzvyux
GzyxV ),,()()()(
1),,(
222
Wir haben diese Integrale gelöst
Kugel Hohlkugel Kugeloberfläche
(für einfache Körper)
Annette Eicker APMG 1
3
11.04.23
Kugel mit homogener Dichteverteilung
z
z
z
z
Vzg z
2
3
3
4
z
GRzg z
zG3
4
2
3
3
4
z
GRzg z
z
GRzV
3
4 3
z
GRzV
3
4 3
3
22
2 zRG
zV
Annette Eicker APMG 1
4
11.04.23
Hohlkugel mit homogener Dichteverteilung
z
z
z
PV r
Pzg r
Annette Eicker APMG 1
5
11.04.23
)(3
4 31
32 RR
r
GV
r
Potential im Außenraum
Potential im Außenraum )( 2Rr
GesamtmasseGesamtmasse
31
323
4RRM
( )GM
Vr
r
r
GRV
3
4 3 r
GesamtmasseGesamtmasse
3
3
4RM
( )GM
Vr
r
Direktes Problem: Das Gravitationspotential ist durch eine gegebene Massenverteilung eindeutig bestimmt
Direktes Problem: Das Gravitationspotential ist durch eine gegebene Massenverteilung eindeutig bestimmt
Inverses Problem: Aus einem Gravitationspotential lässt sich nicht eindeutig auf die Massenverteilung schließen.
Inverses Problem: Aus einem Gravitationspotential lässt sich nicht eindeutig auf die Massenverteilung schließen.
( )r R
Annette Eicker APMG 1
6
11.04.23
Massenverteilung und Potential
Inverses Problem:Aus einem Gravitationspotential lässt sich nicht eindeutig auf die
Massenverteilung schließen.
Inverses Problem:Aus einem Gravitationspotential lässt sich nicht eindeutig auf die
Massenverteilung schließen.
Annahme:Massenveränderungen erfolgen nur auf der Erdoberfläche,
Annahme einer dünnen Schicht
Annahme:Massenveränderungen erfolgen nur auf der Erdoberfläche,
Annahme einer dünnen Schicht
In diesem Fall ist das inverse Problem eindeutig.In diesem Fall ist das inverse Problem eindeutig.
FeldstärkeFeldstärke
)()(3 QQ
QP
QPP dG rr
rr
rrrg
Potential
mit der Flächendichte
Potential
mit der Flächendichte
)()(1
QP
P dGV rrrr
r
2/)( mkgQr
Annette Eicker APMG 1
7
11.04.23
Kugeloberfläche mit konstanter Flächendichte
z
z
z
z
Vzg z
zV
z
GRzV
24
z
GRzV
24
4V z RG
0zg z 2
24
z
GRzg z
2
24
z
GRzg z
Annette Eicker APMG 1
8
11.04.23
Divergenz der GravitationsfeldstärkeDivergenz der Gravitationsfeldstärke
Divergenz und Laplaceoperator
Die Gravitationsfeldstärke lässt sich durch Gradientenbildung aus dem Potential bestimmen
Die Gravitationsfeldstärke lässt sich durch Gradientenbildung aus dem Potential bestimmen
2
2
2
2
2
2
divz
V
y
V
x
V
zV
yV
xV
z
y
x
V
g
Vg
LaplaceoperatorLaplaceoperator
2
2
2
2
2
2
:zyx
Mit Hilfe des Laplaceoperators lassen sich die Quellen des Gravitationsfeldes bestimmen
Mit Hilfe des Laplaceoperators lassen sich die Quellen des Gravitationsfeldes bestimmen
2
2
2
2
2
2
z
V
y
V
x
VV
Annette Eicker APMG 1
9
11.04.23
Laplace- und Poissongleichung
Für beliebige Massenanordnungen gilt:Für beliebige Massenanordnungen gilt:
Außerhalb der Massen:
Laplacegleichung
Außerhalb der Massen:
Laplacegleichung
02
2
2
2
2
2
z
V
y
V
x
VV
Innerhalb der Massen:
Poissongleichung
Innerhalb der Massen:
Poissongleichung
Gz
V
y
V
x
VV 4
2
2
2
2
2
2
Funktionen,die die Laplacegleichung erfüllen,
nennt manharmonische Funktionen.
Funktionen,die die Laplacegleichung erfüllen,
nennt manharmonische Funktionen.
Annette Eicker APMG 1
10
11.04.23
Gaußscher Integralsatz
Gaußscher IntegralsatzGaußscher Integralsatz
dd ngg
n
Integration aller Quellen eines Volumens =Integraler Fluss durch die Oberfläche
Anwendung auf das GravitationsfeldAnwendung auf das Gravitationsfeld
dg
dV
dV
dG4
dG 4
GM4
Gaußsche FormelGaußsche Formel
GMd 4
ng
Annette Eicker APMG 1
11
11.04.23
Repräsentation des Gravitationspotentials
(Kugelfunktionen)
Repräsentation des Gravitationspotentials
(Kugelfunktionen)
Annette Eicker APMG 1
12
11.04.23
Bisher: Für einfache Körper lassen sich auch „einfache“ Formeln angeben.
Kugel Hohlkugel Kugeloberfläche
Jetzt: Für die reale Erde lässt sich das Schwerefeld nicht so leicht beschreiben…
Annette Eicker APMG 1
13
11.04.23
Gravitationspotential
product_type gravity_fieldmodelname ITG-Grace03comment static field from 2002-08 to 2007-04 of GRACE dataearth_gravity_constant 3.986004415e+14radius 6378136.6max_degree 180key L M C S sigma C sigma Send_of_head =========================================================================gfc 0 0 1.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00gfc 1 0 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00gfc 1 1 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00gfc 2 0 -4.841692718699e-04 0.000000000000e+00 6.469883774458e-13 0.000000000000e+00
gfc 2 1 -2.654790999243e-10 1.475393314283e-09 6.108979511966e-13 6.355307212507e-13
gfc 2 2 2.439383367978e-06 -1.400273635220e-06 6.254221806143e-13 6.423410956098e-13
gfc 3 0 9.571610348416e-07 0.000000000000e+00 4.908157850872e-13 0.000000000000e+00
gfc 3 1 2.030461736678e-06 2.482003394707e-07 4.904543816334e-13 5.118675157415e-13
gfc 3 2 9.047877724984e-07 -6.190053685183e-07 5.459595906001e-13 5.482674767117e-13
gfc 3 3 7.213217237276e-07 1.414349090196e-06 5.163836126113e-13 5.163483433061e-13
gfc 4 0 5.399657665980e-07 0.000000000000e+00 3.758481731782e-13 0.000000000000e+00
gfc 4 1 -5.361573220519e-07 -4.735672404588e-07 3.874699557956e-13 3.973550735355e-13
gfc 4 2 3.505015650151e-07 6.624798955603e-07 4.501829959710e-13 4.398486563277e-13
gfc 4 3 9.908565738322e-07 -2.009566568843e-07 4.776067657084e-13 4.766590461153e-13
gfc 4 4 -1.885196275153e-07 3.088038091544e-07 4.556511148108e-13 4.565287154679e-13
gfc 5 0 6.867029195170e-08 0.000000000000e+00 2.444626602968e-13 0.000000000000e+00
gfc 5 1 -6.292117216708e-08 -9.436975416042e-08 2.510438328896e-13 2.640391465060e-13
...
product_type gravity_fieldmodelname ITG-Grace03comment static field from 2002-08 to 2007-04 of GRACE dataearth_gravity_constant 3.986004415e+14radius 6378136.6max_degree 180key L M C S sigma C sigma Send_of_head =========================================================================gfc 0 0 1.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00gfc 1 0 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00gfc 1 1 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00 0.000000000000e+00gfc 2 0 -4.841692718699e-04 0.000000000000e+00 6.469883774458e-13 0.000000000000e+00
gfc 2 1 -2.654790999243e-10 1.475393314283e-09 6.108979511966e-13 6.355307212507e-13
gfc 2 2 2.439383367978e-06 -1.400273635220e-06 6.254221806143e-13 6.423410956098e-13
gfc 3 0 9.571610348416e-07 0.000000000000e+00 4.908157850872e-13 0.000000000000e+00
gfc 3 1 2.030461736678e-06 2.482003394707e-07 4.904543816334e-13 5.118675157415e-13
gfc 3 2 9.047877724984e-07 -6.190053685183e-07 5.459595906001e-13 5.482674767117e-13
gfc 3 3 7.213217237276e-07 1.414349090196e-06 5.163836126113e-13 5.163483433061e-13
gfc 4 0 5.399657665980e-07 0.000000000000e+00 3.758481731782e-13 0.000000000000e+00
gfc 4 1 -5.361573220519e-07 -4.735672404588e-07 3.874699557956e-13 3.973550735355e-13
gfc 4 2 3.505015650151e-07 6.624798955603e-07 4.501829959710e-13 4.398486563277e-13
gfc 4 3 9.908565738322e-07 -2.009566568843e-07 4.776067657084e-13 4.766590461153e-13
gfc 4 4 -1.885196275153e-07 3.088038091544e-07 4.556511148108e-13 4.565287154679e-13
gfc 5 0 6.867029195170e-08 0.000000000000e+00 2.444626602968e-13 0.000000000000e+00
gfc 5 1 -6.292117216708e-08 -9.436975416042e-08 2.510438328896e-13 2.640391465060e-13
...
Wie kann man das Gravitationsfeld der Erde abspeichern?
- Benötigt wird das gesamte Feld im Außenraum (Werte auf der Erdoberfläche, Bahnberechung von Satelliten o. Mond)
Wie kann man das Gravitationsfeld der Erde abspeichern?
- Benötigt wird das gesamte Feld im Außenraum (Werte auf der Erdoberfläche, Bahnberechung von Satelliten o. Mond)
Annette Eicker APMG 1
14
11.04.23
Approximation
Beispiel: Bestimmung der Temperaturabhängigkeit eines EDMsBeispiel: Bestimmung der Temperaturabhängigkeit eines EDMs
x
y
Approximation durch ein Polynom
Vorteile:
- Es müssen nur wenige Polynomkoeffizienten abgespeichert werden.
- Es kann auch zwischen den Messwerten Funktionswerte berechnet werden.
Approximation durch ein Polynom
Vorteile:
- Es müssen nur wenige Polynomkoeffizienten abgespeichert werden.
- Es kann auch zwischen den Messwerten Funktionswerte berechnet werden.
Annette Eicker APMG 1
15
11.04.23
Approximation
Approximation des Potentials durch räumliche PolynomeApproximation des Potentials durch räumliche Polynome
2310
239
338
237
336
352
342
332
323
31
22625
2242322
221
1312110),,(
yzaxzazazyaya
xyzaxyazxayxaxa
zayzayaxzaxyaxa
zayaxaazyxV
Gruppierung in Polynome mit homogenem Grad nGruppierung in Polynome mit homogenem Grad n
0
),,(),,(n
n zyxhzyxf
Homogenes Polynom vom Grad nHomogenes Polynom vom Grad n
( , , ) i k ln nm
i k l n
h x y z a x y z
Beispiel: Homogenes Polynom vom Grad 5Beispiel: Homogenes Polynom vom Grad 5
54225 2315),,( zxyzyxzyxh
Annette Eicker APMG 1
16
11.04.23
Homogene Polynome
Homogenes Polynom vom Grad nHomogenes Polynom vom Grad n
( , , ) i k ln nm
i k l n
h x y z a x y z
Es gilt:Es gilt:
),,(),,( zyxhrrzryrxh nn
n
Beweis:Beweis:
),,( rzryrxhn
n i k lnm
i k l n
r a x y z
( ) ( ) ( )i k lnm
i k l n
a rx ry rz
( )i k l i k l
nmi k l n
a r x y z
),,( zyxhr nn
Es gibtlinear unabhängige homogene Polynome von Grad n
Es gibtlinear unabhängige homogene Polynome von Grad n
)2)(1(21 nn
Beispiel:
n=2
Beispiel:
n=2
12 ( 1)( 2) 6n n
Polynome:
2 2 2, , , , ,x y z xy xz yz
Annette Eicker APMG 1
17
11.04.23
Homogene harmonische Polynome
Das Gravitationspotential ist im Außenraum harmonisch
=> Die Approximation sollte auch harmonisch sein
=> Alle Basisfunktionen sollten harmonisch sein
Das Gravitationspotential ist im Außenraum harmonisch
=> Die Approximation sollte auch harmonisch sein
=> Alle Basisfunktionen sollten harmonisch sein
02
2
2
2
2
2
z
V
y
V
x
VV
Approximation des Potentials durch Summe von homogenen harmonischen Polynomen
Approximation des Potentials durch Summe von homogenen harmonischen Polynomen
0
),,(),,(n
n zyxHzyxV
Homogenes harmonisches Polynom vom Grad n dargestellt als Linearkombination von linear unabhängigen Basispolynomen
Homogenes harmonisches Polynom vom Grad n dargestellt als Linearkombination von linear unabhängigen Basispolynomen
m
nmnmn zyxHazyxH ),,(),,(
Beispiel: Grad n=2Beispiel: Grad n=2
22225
2224
23
22
21
2),,(
),,(
),,(
),,(
),,(
xyzzyxH
yxzyxH
yzzyxH
xzzyxH
xyzyxH
Wie viele Basisfunktionen gibt es?Wie viele Basisfunktionen gibt es?
Annette Eicker APMG 1
18
11.04.23
Homogene harmonische Polynome
Satz: Der Laplace Operator angewendet auf ein homogenes Polynom hn liefert ein homogenes Polynom vom Grad n-2
Satz: Der Laplace Operator angewendet auf ein homogenes Polynom hn liefert ein homogenes Polynom vom Grad n-2
Es gibt 2n+1 linear unabhängige harmonische homogene Polynome Hnm von Grad nEs gibt 2n+1 linear unabhängige harmonische homogene Polynome Hnm von Grad n
BeweisBeweis
( , , ) i k ln nm
i k l n
h x y z a x y z
2 2 2( , , ) ( 1) ( 1) ( 1)i k l i k l i k l
n nm nm nmi k l n
h x y z a i i x y z a k k x y z a l l x y z
Es gibt linear unabhängige homogene Polynome von Grad n
Das Ergebnis der Laplace Operation ist durch linear unabhängige homogene Polynome von Grad n-2 eindeutig bestimmt.
Es verbleiben linear unabhängige homogene harmonische Polynome.
Es gibt linear unabhängige homogene Polynome von Grad n
Das Ergebnis der Laplace Operation ist durch linear unabhängige homogene Polynome von Grad n-2 eindeutig bestimmt.
Es verbleiben linear unabhängige homogene harmonische Polynome.
)2)(1(21 nn
)2)2)((1)2((21 nn
12)2)2)((1)2(()2)(1( 21
21 nnnnn
Annette Eicker APMG 1
19
11.04.23
Homogene harmonische Polynome
Das Gravitationspotential ist im Außenraum harmonisch
=> Die Approximation sollte auch harmonisch sein
=> Alle Basisfunktionen sollten harmonisch sein
Das Gravitationspotential ist im Außenraum harmonisch
=> Die Approximation sollte auch harmonisch sein
=> Alle Basisfunktionen sollten harmonisch sein
02
2
2
2
2
2
z
V
y
V
x
VV
Approximation des Potentials durch Summe von homogenen harmonischen Polynomen
Approximation des Potentials durch Summe von homogenen harmonischen Polynomen
0
),,(),,(n
n zyxHzyxV
Homogenes harmonisches Polynom vom Grad n dargestellt als Linearkombination von linear unabhängigen Basispolynomen
Homogenes harmonisches Polynom vom Grad n dargestellt als Linearkombination von linear unabhängigen Basispolynomen
12
1
),,(),,(n
mnmnmn zyxHazyxH
Beispiel: Grad n=2Beispiel: Grad n=2
22225
2224
23
22
21
2),,(
),,(
),,(
),,(
),,(
xyzzyxH
yxzyxH
yzzyxH
xzzyxH
xyzyxH
n
nmnmnmn zyxHazyxH ),,(),,(
oder wenn man negative Indizies einführt:
Grad nGrad n Ordnung mOrdnung m
Annette Eicker APMG 1
20
11.04.23
Kugelflächenfunktionen
Sphärische PolarkoordinatenSphärische Polarkoordinaten
cos
sinsin
sincos
r
r
r
z
y
x
r
homogene harmonische Polynomehomogene harmonische Polynome
),(
)cos,sinsin,sin(cos
)cos,sinsin,sincos(),,(
nmn
nmn
nmnm
Yr
Hr
rrrHzyxH
Darstellung des PotentialsDarstellung des Potentials
0
),(),,(n
n
nmnmnm
n YarrV
Approximation von Funktionen auf der KugelApproximation von Funktionen auf der Kugel
0
),(),(n
n
nmnmnmYaf
Kugelflächenfunktion:homogenes harmonisches
Polynom beschränkt auf die Kugeloberfläche
Kugelflächenfunktion:homogenes harmonisches
Polynom beschränkt auf die Kugeloberfläche
Annette Eicker APMG 1
21
11.04.23
Approximation durch Kugelflächenfunktionen
Approximation durch Kugelflächenfunktionen
Annette Eicker APMG 1
22
11.04.23
Approx. durch Kugelflächenfunktionen
Grad nAnzahl der
Koeffizienten
4 25
8 81
16 289
30 961
60 3721
120 14641
240 58081
0
),(),(n
n
nmnmnmYaf
Annette Eicker APMG 1
23
11.04.23
Approx. durch Kugelflächenfunktionen
Grad nAnzahl der
Koeffizienten
4 25
8 81
16 289
30 961
60 3721
120 14641
240 58081
0
),(),(n
n
nmnmnmYaf
Annette Eicker APMG 1
24
11.04.23
Approx. durch Kugelflächenfunktionen
Grad nAnzahl der
Koeffizienten
4 25
8 81
16 289
30 961
60 3721
120 14641
240 58081
0
),(),(n
n
nmnmnmYaf
Annette Eicker APMG 1
25
11.04.23
Approx. durch Kugelflächenfunktionen
Grad nAnzahl der
Koeffizienten
4 25
8 81
16 289
30 961
60 3721
120 14641
240 58081
0
),(),(n
n
nmnmnmYaf
Annette Eicker APMG 1
26
11.04.23
Approx. durch Kugelflächenfunktionen
Grad nAnzahl der
Koeffizienten
4 25
8 81
16 289
30 961
60 3721
120 14641
240 58081
0
),(),(n
n
nmnmnmYaf
Annette Eicker APMG 1
27
11.04.23
Approx. durch Kugelflächenfunktionen
Grad nAnzahl der
Koeffizienten
4 25
8 81
16 289
30 961
60 3721
120 14641
240 58081
0
),(),(n
n
nmnmnmYaf
Annette Eicker APMG 1
28
11.04.23
Approx. durch Kugelflächenfunktionen
Grad nAnzahl der
Koeffizienten
4 25
8 81
16 289
30 961
60 3721
120 14641
240 58081
0
),(),(n
n
nmnmnmYaf
Annette Eicker APMG 1
29
11.04.23
Kugelfunktionen
Approximation des PotentialsApproximation des Potentials
0
),(),,(n
nnYrrV
Laplacesche KugelflächenfunktionenLaplacesche Kugelflächenfunktionen
n
nmnmnmn YaY ),(),(
Die Reihe konvergiert nur für r<1Die Reihe konvergiert nur für r<1
Approximation des PotentialsApproximation des Potentials
01
),(1
),,(n
nnY
rrV
Diese Reihe ist ebenfalls harmonischund konvergiert für r>1
Beweis: nächste Folien
Diese Reihe ist ebenfalls harmonischund konvergiert für r>1
Beweis: nächste Folien
Annette Eicker APMG 1
30
11.04.23
Laplace OperatorLaplace Operator
Laplace und Beltrami Operator
Laplace Operator in sphärischen KoordinatenLaplace Operator in sphärischen Koordinaten
2
2
2
2
2
2
z
V
y
V
x
VV
2
2
2222
2
22
2
sin
1
tan
112
V
r
V
r
V
rr
V
rr
VV
2
2
22
2
22
2
sin
1
tan
112
VVV
rr
V
rr
V
Vrr
V
rr
V *22
2 12
Beltrami Operator (Laplace Operator beschränkt auf die Kugeloberfläche)Beltrami Operator (Laplace Operator beschränkt auf die Kugeloberfläche)
2
2
22
2*
sin
1
tan
1
VVV
V
Annette Eicker APMG 1
31
11.04.23
Kugelfunktionen
Approximation des PotentialsApproximation des Potentials
0
),(),,(n
nnYrrV
Die Reihe konvergiert nur für r<1Die Reihe konvergiert nur für r<1
Approximation des PotentialsApproximation des Potentials
01
),(1
),,(n
nnY
rrV
Diese Reihe ist ebenfalls harmonischund konvergiert für r>1
Beweis: nächste Folien
Diese Reihe ist ebenfalls harmonischund konvergiert für r>1
Beweis: nächste Folien
Wir wissen, dass die Laplacegleichung
hierfür gilt
Wir wissen, dass die Laplacegleichung
hierfür gilt
2*
2 2
2 10
V VV V
r r r r
Ziel: zeigen, dass dann die
Laplacegleichung auch
dafür gilt!
Annette Eicker APMG 1
32
11.04.23
Laplace OperatorLaplace Operator
Laplace und Beltrami Operator
Vrr
V
rr
VV *
22
2 12
0),( nnYr
),(11
nnY
r*
3 3 3
( 1)( 2) 2( 1) 1( , ) ( , ) ( , )n n nn n n
n n nY Y Y
r r r
*3 3
( 1) 1( , ) ( , )n nn n
n nY Y
r r
2 2 2 *( 1) ( , ) 2 ( , ) ( , )n n nn n nn n r Y nr Y r Y
2 2 *( 1) ( , ) ( , ) 0n nn nn n r Y r Y
( 2) 3: :n nr r
*3 3
( 1) 1( , ) ( , ) 0n nn n
n nY Y
r r
0
* ( , ) ( 1) ( , )n nY n n Y
Annette Eicker APMG 1
33
11.04.23
Kugelfunktionen
Approximation des Potentials für r<1Approximation des Potentials für r<1
0
),(),,(n
nnYrrV
Laplacesche KugelflächenfunktionenLaplacesche Kugelflächenfunktionen
n
nmnmnmn YaY ),(),(
Approximation des Potentials für r>1Approximation des Potentials für r>1
01
),(1
),,(n
nnY
rrV
Kugelflächenfunktionen
Annette Eicker APMG 1
34
11.04.23
Herleitung der Kugelfunktionen über die Lösung der Laplacegleichung
Herleitung der Kugelfunktionen über die Lösung der Laplacegleichung
Annette Eicker APMG 1
35
11.04.23
Gesucht sind alle Lösungen der LaplacegleichungGesucht sind alle Lösungen der Laplacegleichung
Lösung der Laplace Gleichung
012 *22
2
Vrr
V
rr
VV
Annahme: Lösung entspricht folgendem SeparationsansatzAnnahme: Lösung entspricht folgendem Separationsansatz
),()(),,( YrfrV
In die Laplacegleichung eingesetzt:In die Laplacegleichung eingesetzt:
0),()(1),()(2),()( *22
2
Yrf
rr
Yrf
rr
Yrf
0),()()(),(2)(
),( *22
2
Y
r
rf
r
rf
r
Y
r
rfY
),(),(
1)(
)(
2)(
)(*
2
22
YYr
rf
rf
r
r
rf
rf
r
Beltrami OperatorBeltrami Operator
2
2
22
2*
sin
1
tan
1
YYY
Y
2
( ): : ( , )f r
Yr
Annette Eicker APMG 1
36
11.04.23
Lösung der Laplace Gleichung
),(),(
1)(
)(
2)(
)(*
2
22
YYr
rf
rf
r
r
rf
rf
r
Die linke Seite hängt nur von r ab, die rechte nur von
=> Damit beide Seiten immer gleich sind, müssen beide Seiten konstant sein.
=> Wahl der Konstanten
Die linke Seite hängt nur von r ab, die rechte nur von
=> Damit beide Seiten immer gleich sind, müssen beide Seiten konstant sein.
=> Wahl der Konstanten
),(
)1( nn
Linke Seite:Linke Seite:
)1()(
)(
2)(
)( 2
22
nnr
rf
rf
r
r
rf
rf
r
)()1()(
2)(
2
22 rfnn
r
rfr
r
rfr
DGL hat zwei linear unabhängige LösungenDGL hat zwei linear unabhängige Lösungen
nrrf )(1)1(
2 )( nrrf
Rechte Seite:Rechte Seite:
),(),(
1)1( *
Y
Ynn
),()1(),(* YnnY
Annette Eicker APMG 1
37
11.04.23
Lösung der Laplace Gleichung
Rechte Seite:Rechte Seite:
),()1(),(* YnnY
Annahme: Lösung entspricht folgendem SeparationsansatzAnnahme: Lösung entspricht folgendem Separationsansatz
)()(),( pgY
Beltrami OperatorBeltrami Operator
2
2
22
2*
sin
1
tan
1
YYY
Y
YnnYYY
)1(sin
1
tan
12
2
22
2
eingesetzt:eingesetzt:2 2
2 2 2
( ) ( ) ( ) ( ) ( )( ) ( 1) ( ) ( )
tan sin
p g p p gg n n g p
2
2
2
22 )(
)(
1)()1(
)(
sin
cos)(
)(
sin
g
gpnn
pp
p
2
: ( )
sin
( )
g
p
Annette Eicker APMG 1
38
11.04.23
Linke Seite:Linke Seite:
22
22
)()1()(
sin
cos)(
)(
sinmpnn
pp
p
Lösung der Laplace Gleichung
Die linke Seite hängt nur von ab, die rechte nur von
=> Damit beide Seiten immer gleich sind, müssen beide Seiten konstant sein.
=> Wahl der Konstanten
Die linke Seite hängt nur von ab, die rechte nur von
=> Damit beide Seiten immer gleich sind, müssen beide Seiten konstant sein.
=> Wahl der Konstanten
2m
Rechte Seite:Rechte Seite:
22
2 )(
)(
1m
g
g
)()( 2
2
2
gm
g
DGL hat zwei linear unabhängige Lösungen
DGL hat zwei linear unabhängige Lösungen
mg cos)(1 mg sin)(2
0)(sin
sin)1()(
cos)(
sin2
2
2
p
mnn
pp
2
2
2
22 )(
)(
1)()1(
)(
sin
cos)(
)(
sin
g
gpnn
pp
p
Lösung der DGL sind die sogenanntenzugeordneten legendreschen Polynome
Lösung der DGL sind die sogenanntenzugeordneten legendreschen Polynome
)(cos)( mnPp
Annette Eicker APMG 1
39
11.04.23
DGL hat zwei linear unabhängige Lösungen
DGL hat zwei linear unabhängige Lösungen
mg cos)(1 mg sin)(2
Linke Seite:Linke Seite:
Lösung der Laplace Gleichung
Die linke Seite hängt nur von ab, die rechte nur von
=> Damit beide Seiten immer gleich sind, müssen beide Seiten konstant sein.
=> Wahl der Konstanten
Die linke Seite hängt nur von ab, die rechte nur von
=> Damit beide Seiten immer gleich sind, müssen beide Seiten konstant sein.
=> Wahl der Konstanten
2m
Rechte Seite:Rechte Seite:
22
2 )(
)(
1m
g
g
)()( 2
2
2
gm
g
0)(sin
sin)1()(
cos)(
sin2
2
2
p
mnn
pp
2
2
2
22 )(
)(
1)()1(
)(
sin
cos)(
)(
sin
g
gpnn
pp
p
22
22
)()1()(
sin
cos)(
)(
sinmpnn
pp
p
Lösung der DGL sind die sogenanntenzugeordneten legendreschen Polynome
Lösung der DGL sind die sogenanntenzugeordneten legendreschen Polynome
)(cos)( mnPp
00 (cos ) 1P
01 (cos ) cosP 11 (cos ) sinP
0 22 (cos ) 1/ 2(3cos 1)P 12 (cos ) 3sin cosP 2 22 (cos ) 3sinP
03 (cos )P
Annette Eicker APMG 1
40
11.04.23
Lösung der Laplace Gleichung
Gesucht sind alle Lösungen der Laplacegleichung
Gesucht sind alle Lösungen der Laplacegleichung
0V
Spezielle Lösung
mit
Spezielle Lösung
mit
)()()(),,( pgrfrV
nrrf )(1)1(
2 )( nrrf
mg cos)(1
mg sin)(2
)(cos)( mnPp
Die allgemeine Lösung ist eine Linearkombination aller LösungenDie allgemeine Lösung ist eine Linearkombination aller Lösungen
mPsmPcr
rVn
n
m
mnnm
mnnmn
sin)(coscos)(cos1
),,(0 0
1
mPsmPcrrVn
n
m
mnnm
mnnm
n sin)(coscos)(cos),,(0 0
Vergleich: Fourier-ReiheVergleich: Fourier-Reihe
1
( ) cos( ) sin( )n nn
f t a n t b n t
Kugelfunktionsentwicklung kann man interpretieren als 2D Fourier-Reihe auf der Kugel