Elliptische Partielle Differentialgleichung

-Δ u(x,y)=f(x,y)

u(0,y)=u(L,y)=u(x,0)=u(x,L)=0

x,y ∈[0,L]×[0,L]

Lösen des Eigenwertproblems

"ContinuousFourierTransform2D_1.gif"

weil wenn man eine solche Funktion finden würde, dann wäre die Lösung von -Δ u(x,y)=f(x,y) ja einfach u(x,y)=f(x,y)/λ.

Als Ansatz u(x,y)=X(x) Y(y), dann wird aus der Gleichung

"ContinuousFourierTransform2D_2.gif"

X''(x)Y(x)+X(x)Y''(y)+λ X(x)Y(y)=0

oder

X''(x)/X(x)+Y''(y)/Y(y)+λ=0

-X''(x)/X(x)=Y''(y)/Y(y)+λ=const=Λ

X''(x)+      Λ          X(x)=0   mit X(0)=X(L)=0
Y''(x)+(Λ-λ) Y(x)=0   mit Y(0)=Y(L)=0

Man muss also nur einmal das Eigenwertproblem
u''(x)+λ  u(x)=0
mit den Randbedingungen u(0)=u(L)=0 lösen. Für die Differentialgleichung gibr es als Lösungen

"ContinuousFourierTransform2D_3.gif"

Für ein allgemeines λ lassen sich die Randbedingungen nicht erfüllen, ("ContinuousFourierTransform2D_4.gif" muss 0 sein sein wegen der Bedingung u(0)=0), die zweite Randbedingung lässt sich nur erfüllen wenn λ spezielle Werte an nimmt, nämlich

"ContinuousFourierTransform2D_5.gif"

Es gibt abzählbar unendlich viele Lösungen "ContinuousFourierTransform2D_6.gif"  die alle nur bis auf eine multiplikative Konstante bestimmt sind, n=1,2,3..

die Eigenfunktionen des zweidimensionalen Problems sind dann "ContinuousFourierTransform2D_7.gif" und die Eigenwerte "ContinuousFourierTransform2D_8.gif"

für die Randbedingungen u'(0)=u'(L)=0  wären die Eingenfunktionen "ContinuousFourierTransform2D_9.gif" mit n=1,2,3... und "ContinuousFourierTransform2D_10.gif" gefunden worden

für die periodischen Randbedingungen u'(0)=u'(L)  ergibt sich "ContinuousFourierTransform2D_11.gif" , da man nicht immer die zwei Konstanten "ContinuousFourierTransform2D_12.gif" und "ContinuousFourierTransform2D_13.gif" schreiben will, wird diese Eigenfunktion zu "ContinuousFourierTransform2D_14.gif"

die Lösung des elliptischen Problems ergibt sich als Fourier-Reihe "ContinuousFourierTransform2D_15.gif"

Hilbert-Raum

Jeder Satz Eigenfunktionen eines elliptischen Eigenwertproblems bildet eine vollständige Basis im Hilbertraum der Funktionen

Das Skalarprodukt wird durch ein Integral definiert

"ContinuousFourierTransform2D_16.gif"

also im eindimensionalen über dem Interval x∈[0,L]

"ContinuousFourierTransform2D_17.gif"

und im zweidimensionalen x,y∈[0,L]×[0,L]

"ContinuousFourierTransform2D_18.gif"

Jede Funktion, für die das Skalarprodukt berechnet werden kann kann als Reihe

"ContinuousFourierTransform2D_19.gif" mit "ContinuousFourierTransform2D_20.gif" dargestellt werden. Analog gilt im zweidimensionalen
"ContinuousFourierTransform2D_21.gif" mit "ContinuousFourierTransform2D_22.gif".

"ContinuousFourierTransform2D_23.gif"

"ContinuousFourierTransform2D_24.gif"

"ContinuousFourierTransform2D_25.gif"

"ContinuousFourierTransform2D_26.gif"

"ContinuousFourierTransform2D_27.gif"

"ContinuousFourierTransform2D_28.gif"

f(x,y)

"ContinuousFourierTransform2D_29.gif"

Graphics:ρ(x, y)

"ContinuousFourierTransform2D_31.gif"

"ContinuousFourierTransform2D_32.gif"

"ContinuousFourierTransform2D_33.gif"

"ContinuousFourierTransform2D_34.gif"

"ContinuousFourierTransform2D_35.gif"

"ContinuousFourierTransform2D_36.gif"

u(x,y)

"ContinuousFourierTransform2D_37.gif"

"ContinuousFourierTransform2D_38.gif"

Graphics:u(x, y)

Spikey Created with Wolfram Mathematica 7.0