2.6 Die Wellengleichung auf einem endlichen Intervall

Es liegt das Problem

                                 2
                       Lu = utt- c uxx = f(x,t)        f¨ur  < x < t, t > 0
                    {  u(x,0) = f(x), ut(x,0) = y(x)  f¨ur  0 < x < l
(P) = P (f,f,y, B1,B2)  u(0,t) = B1(t), u(l,t) = B2(t)  f¨ur  t > 0
                       [ux(0,t) = B1(t), ux(l,t) = B2(t)] f¨ur  t > 0
                       [aux(0,t)+ a2u(0,t) = B1(t)]

vor. Würden Lösungen in C2((0,l),x(t > 0))  /~\ C1([0,l] × (t > 0)) gesucht, so hat man die Verträglichkeitsbedingungen zu beachten.

                                '           '
f(0) = B1(0), f(l) = B2(0), y(0) = B 1(0), y(l) = B 2(0)

Führen wir diese Schritte nun ausführlich aus. Der erste Schritt ist eine Transformation der Randbedingungen zu Null. Wir verwenden den Ansatz u(x,t) = U(x,t) + w(x,t), welcher das Problem für w liefert:

Lw(x,t) = f(x,t)- LU (x,t) = f~(x,t)

w(x,0) = f(x)- U (x,0) = ~f(x); wt(x,0) = y(x)- Ut(x,0) = ~y(x)

w(0,t) = B1(t)- U (0,t)=!0
                         }
w(l,t) = B2(t)- U(l,t)=!0   Bedingung an U

Wähle U(x,t) = B1(t) + x
l(B2(t)- B1(t)). Damit löst dann w folgendes Problem mit homogenen Randbedingungen:

Lw(x, t) = ~f(x,t)

w(x,0) = ~f(x), wt(x,0) = ~y(x)

w(0,t) = w(l,t) = 0

Nach Schritt 2 Wir zerlegen in P1(0,f,0,0,0), P2(0,0,y,0,0) und P3(f,0,0,0,0):

  {  Lu = 0                     {  Lu = 0                      { Lu = f
P1   u(x,0) = f(x), ut(x,0) = 0 , P2 u(x,0) = 0, ut(x,0) = y(x) , P3 u(x,0) = ut(x,0) = 0
     u(0,t) = u(l,t) = 0             u(0,t) = u(l,t) = 0             u(0,t) = u(l,t) = 0

Die Probleme P1 und P2 beschreiben die freie Schwingung einer eingespannten Saite (Separationsmethode). Bei P3 handelt es sich um eine erzwungene Schwingung (äußere Kräfte) (DUHAMEL-Methode). Wir verwenden nach Schritt 3 den Ansatz u(x,t) = a(x)b(t) für P1, P2. Durch Einsetzen folgt dann:

a(x)¨b(t)- c2a''(x)b(t) = 0

   ¨      ''
-12 b(t)-= a-(x)= m = const.
c  b(t)   a(x)

 ''
a (x)- ma(x) = 0 f¨ur 0 < x < l, a(0) = a(l) = 0

¨b(t)- c2mb(t) = 0 fu¨r t > 0, b(0) = 0

Die Lösung des Randwertproblems ergibt sich zu a(x) = exp(cx). Daraus erhalten wir c2 = m und somit c = ± V~ -
 m. Die allgemeine Lösung lautet damit:

              V~ --            V~ --
a(x) = C1exp ( mx)+ C2 exp(-  mx)

Aus a(0) = C1 + C1!
= 0 folgt C 1 = -C2 und weiterhin:

|------------- V~ ---------- V~ ----|
a(x)-=-C1[exp-(-mx)--exp-(----mx)]-

    !            V~ -          V~ -          V~ -          V~ -
a(l)= 0 = C1[exp( ml)- exp (-  ml)] ==> [exp( ml)- exp(- ml)] = 0

Aus exp(  V~ -)
 2  ml = 1 = exp(2pi.n) erhalten wir  V~ -
 ml = npi und daraus wiederum  V~ --
  m = nlpi.

|---------[---(np----)------(--np---)]-------(np--)-----------(np--)-----------------|
an(x) = C1 exp  --i.x  - exp - --i.x   = ~C sin  --x  , a(x) = sin--x  f¨ur n = 1, 2, 3, ...
-----------------l--------------l---------------l---------------l---------------------

Mit m = - 2 2
nl2p- lautet das Problem für b(t):

        2 2
¨b(t)+ c2n-p2--b(t) = 0 mitb˙(0) = 0
         l

Diese Differentialgleichung besitzt folgende Lösung:

         (     )
bn(t) = cos npct  f¨ur n = 1, 2, 3,...
           l

Wir erhalten schließlich das Gesamtergebnis:

|-----------(np--)----(np--)-----------------|
|un(x,t) = sin  -l x cos -l-ct  f¨ur n = 1, 2, 3, ...
---------------------------------------------|

Des weiteren ist die Bedingung u(x,0) = f(x) zu erfüllen. Wir machen dabei folgenden Ansatz: Gesucht sind Zahlen An (n = 1, 2, 3, ...) derart, daß für f(x) gilt:

       sum  oo         ||      sum  oo      (np  )   (np  )||
f(x) =   Anun(x, t)||   =    An sin  ---x cos ---ct||
      n=1         |t=0   n=1        l        l   |t=0

Falls die Reihe „genügend konvergent“ ist, gilt:

        oo 
f(x) =  sum  A sin(np-x)
      n=1  n     l

Hierbei handelt es sich gerade um die Fourierreihe von f(x). Das ist erfüllt, wenn die An die Fourierkoeffizienten der nach -l < x < x ungeraden fortgesetzten 2l-periodischen Funktion f(x) mit f|[0,l]f sind.

-----------------------------------------------
|       integral l                                    |
|    2          (np- )                        |
An =  l  f(q)sin   l q dq = ^fn f¨ur n = 1, 2, 3, ...
-------0---------------------------------------

Hierbei handelt es sich um die Fouriertransformierte. Die Lösung von P1 ist dann:

|---------------------------------|
|         oo  sum       (np  )    (np  ) |
u1(x,t) =   ^fn sin  l-x  cos  -l ct |
---------n=1-----------------------

u1xx, u1tt muß eine konvergente Reihe liefern, also muß  sum n=1 oo n2|f^ n|2 konvergieren. Als Übung kann durch viermaliges partielles Integrieren gezeigt werden:

     integral t      (np  )
^fn =   f(q) sin  --q  dq
    0          l

Beim Lösen des Problems P2 erhalten wir analog folgende Bedingungen:

      n2p2-      n2p2-2
m = -  l2 , b(t)+ l2 c b(t) = 0 mit b(0) = 0

Diese Gleichung wird gelöst durch:

         (np   )
bn(t) = sin-l-ct

                    |--(----)---(-----)-|
un(x,t) = an(x)bn(t) =|sin npx  sin  npct  |
                    ------l--------l----|

Dies erfüllt dann alle Bedingungen bis auf y(x) = ut(x,0). Gesucht sind also Zahlen Bn mit:

        oo  sum           (    )
y(x) =    Bn npcsin np-x
       n=1   l       l

Dies läßt sich wieder durch Fouriertransformation verwirklichen:

|--------------------------------------|
|  np               2 integral l      (np  )   |
|Bn---= f^n mit ^yn = -  y(q) sin  --q  dq|
|   l               l0          l      |
----------------------------------------

Damit gilt für das Problem P2:

|----------------------------------|
|         sum  oo  ^y    (np  )   (  p  )|
u2(x,t) =   --0-lsin  --x  sin  n-ct |
---------n=1npc------l---------l----

Durch Addition der Einzellösungen u1(x,t) und u2(x,t), also durch

                    sum  oo  [           ^        ]   (    )
(u1 + u2) = u0(x,t) =  f^n cos(wnt)+ yn-sin(wnt) sin np-x
                   n=1             wn              l

ergibt sich die Lösung von:

Lu0 = 0

u0(x,0) = f(x), ut(x,0) = y(x)

u0(0,t) = u0(l,t) = 0

durch Superposition unendlich vieler Lösungen (Reihen). Hierbei handelt es sich um eine Entwicklung nach stehenden Wellen. Wir setzen nun:

            ^
^fncos(wnt)+  ynsin(wnt) = ancoswn(t+ dn)
            wn

Wie berechnet man nun die an und dn? Dazu benötigen wir erst einmal die Additionstheoreme:

an coswn(t+ dn) = an cos(wnt)cos(wndn)- an sin(wnt)sin(wndn)

Durch Koeffizientenvergleich ergibt sich:

                      y^n
^fn = an cos(wndn) und - w-= an sin(wndn)
                       n

dn folgt aus der Division und anschließenden Anwendung des Arkustangens. Die an erhält man durch Quadrieren und Addition beider Gleichungen. Daraus resultiert dann:

          oo  sum                   (np- )
u0(x,t) =    an coswn(t+ dn)sin  l x
         n=1

                           (    )
u(0n)(x,t) = an coswn(t +dn)sin np-x
                             l

Dabei handelt es sich um eine stehende Welle:

Bei u0(n)(x,t) schwingt jeder Punkt mit derselben Frequenz wn = np
-lc. Des weiteren gilt ja bekanntlich wn = 2pnn = 2p
cnc. Daraus folgen dann die möglichen Wellenlängen cn und Frequenzen nn auf der Saite der Länge l:

|---------------------------|
cn = 2p-c = 2l, nn =-c-= n-c|
-----wn----n-------cn----2l--

n1 = c2l- wird Grundfrequenz genannt; die anderen Frequenzen ergeben sich durch nn = nn1.

                      V~ --            V~ -
                       T-         n-  T-
rutt- Tuxx = 0 mit c = r und nn = 2l  r

u0 löst P(0,f,y,0,0), also utt - c2uxx = 0. Die allgemeine Lösung der Gleichung ist u(x,t) = F(x - ct) + G(x + ct) (Entwicklung nach fortschreitenden Wellen).

u(n)(x,t) = ^f cos(w t)sin(np-x)+  ^yn-sin(w t)sin(np-x)=
 0         n     n      c      wn     n       l
          1   [  (np        )     (np       )]   1^yn [   (np       )     (np        )]
        = 2 ^fn sin  -l (x+ ct) + sin -l-(x - ct)  +  2w-- cos -l-(x - ct) - cos  l-(x+ ct)
                                                   n
(2.8)

Für die Lösung von P2(0,0,y,0,0) gilt:

            [          ]
          oo  sum   y^n            (np  )
u2(x,t) =    wn-sin(wnt) sin -l-x
         n=1

Wir betrachten nun wieder P3(f,0,0,0,0):

                 Lu = f(x,t)         f¨ur  0 < x < l, t > 0
              {  u(0,t) = u(l,t) = 0   f¨ur  t > 0
P3(f,0,0,0,0) =   u(x,0) = ut(x,0) = 0 f¨ur  0 < x < l

Hier bietet sich dann die Methode von DUHAMEL an. Löse für 0 < j < t das Problem P:

      Lu = 0                     f¨ur  0 < x < l, t > j
    { u(0,t) = u(l,t) = 0         f¨ur  t > j
P =   u(x,j) = 0, ut(x,j) = f (x,j) f¨ur  0 < c < l

Man erhält dadurch also das Problem P2 mit verschiedenen Anfangsbedingungen.

           oo  sum  ^fn(j)             (np--)
u(x,t,j) =     wn  sin wn(t- j)sin  l x
         n=1

         integral  l       (    )
^fn(j) = 2   f(q,j)sin np-q  dq
       l 0           l

Mit dem Prinzip von DUHAMEL erhält man nun:

          integral t
u3(x,t) =    u(x,t,j)dj
        j=0

|-----------------------------------------------------------------|
|         integral  t integral  l[o sum  o  2 1              (np  )   (np )]            |
u3(x,t) =           -.---sin wn(t- j)sin  --x  sin  ---q  f(q,j)dqdj |
|        j=0 q=0  n=1-l-wn----------- -----l--------l----           |
|                              G(x,q,t-j)                           |
------------------------------------------------------------------

G(x,q,t - j) ist die sogenannte GREENsche Funktion des hier betrachteten Anfangswert-Randwert-Problems.

 integral   integral 
    (LG)f  = f(x,t)

Wir machen die Probe:

u   - c2u    = f(x,t)
 3tt     3xx

          @  integral t  integral l[o o  sum  2            (np  )   (np  )]
u3tt(x,t) = --          -coswn(t- j)sin ---x sin ---q  f (q,j)dqdj =
          @t0 q=0  n=1 l                l        l
           integral l [                     ]            integral t integral  l[                                  ]
               oo  sum  2   (np--)   (np--)                  o sum  o  2                (np- )   (np- )
        =        l sin  l x  sin   l q  f(q,j)dq -           lwnsinwn(t- j)sin   l x sin   l q  f(q,j)
         q=0  n=1                               0 q=0  n=1
(2.9)

             integral t  integral l[                                       ]
                  o sum  o  2-1-n2p2-             (np- )   (np- )
u3xx(x,t) = -          lwn  l2 sinwn(t- j)sin   l x sin   l q  f(q,j)dqdj
            0 q=0

                integral t  integral l[o sum  o  2 1 n2p2             (np  )    (np )]
-c2u3xx(x,t) = +          ----c2 -2--sin wn(t- j)sin  --x  sin  --q   f(q,j)dqdj
               0 q=0     lwn    l                  l        l

Darüber hinaus gilt ja nun:

  n2p2  1        1
c2--2-.---= w2n .--= wn
   l   wn       wn

Damit erhalten wir schließlich:

               oo   |_   integral l                 _| 
u   -c2u   =  sum    2    f(q,t)sin (np-q) dq  sin(np-x)=  f(x,t)
 3tt     3xx  n=1  |_ l             l      _|      l
                   q=0

Es handelt sich gerade um die Fouriertransformierte von f(x,t). Die Probe geht folglich auf!

Problem P(0,f, Y,B1,B2):

Wir behandeln nun die zweite Methode zur Lösung von P(0,f,Y,B1,B2):

PIC PIC

2.6.1 Transmissionsprobleme

Angenommen, wir haben zwei Saiten mit verschiedenen physikalischen Eigenschaften (r, c), die aber an einer Stelle fest miteinander verbunden sind: Die Wellengleichungen für getrennten Probleme lauten dann:

     { utt- c21uxx = 0 f¨ur x < 0 und t > 0
Lu =   utt- c22uxx = 0 f¨ur x > 0 und t > 0 mit u(x,0) = ut(x,0) = 0 f¨ur x > 0

Die Lösung des Anfangswertproblems utt - c2uxx = 0 (für x  (- R, t > 0) ist:

            (   x)                            (  x)
        { P  t- c   fu¨r  x < ct      u(x,0) = P - c   }
u(x,t) =                        mit            '(  x)   x < 0
          0         fu¨r  x > ct      ut(x,0) = P - c

u(x,0) = ut(x,0) = 0 f¨ur x > 0

PIC

         [  (       )    ]      integral 0  (   )        (     )     (    )
u(x,t) = 1 P  - x--ct  + 0 + 1-    P' - s  ds = 1P t-  x-+ 1P  t- x-- 1P(0)
        2        c          2cx- ct      c      2       c   2      c   2- -
                                                                       !=0

Gesucht ist

        {  u(1)(x,t)  f¨ur  x < 0, t > 0
u(x,t) =   u(2)(x,t)  f¨ur  x > 0, t > 0

mit utt -c12uxx = 0 für t > 0, x < 0 und utt -c22uxx = 0 für t > 0, x > 0 mit den folgenden Bedingungen:

          (    )
u(x,0) = P - x-  f¨ur x < 0, u(x,0) = 0 f¨ur x > 0
             c1

           (  x )
ut(x,0) = P' ---  f¨ur x < 0, ut(x,0) = 0 f¨ur x > 0
              c1

Die Übergangsbedingungen (Matching Conditions) stellen nun einen Zusammenhang zwischen den beiden Problemen her:

u(1)(0,t) = u(2)(0,t)

u(x1)(0,t) = u(x2)(0,t)

Wir können die allgemeine Lösung der beiden Gleichungen hinschreiben:

            (      )     (      )
u(1)(x,t) = F1 t- x-  + G1  t+ x-
                 c1           c1

            (    x )     (    x )
u(2)(x,t) = F2 t- c-  + G2  t+ c-
                  2            2

Zu bestimmen sind nun noch die Funktionen F1 und G1. Mittels der Übergangsbedingungen folgt nun:

u(1)(0,t) = P(t)+ G (t) != F (t) = u(2)(0,t)
                 1      2

u(1)(0,t)- 1-P'(t)+ -1G' (t)=!- -1F '(t) = u(2)(0,t)
 x        c1      c1  1      c2 2      x

Man erhält dann folgende einfache Differentialgleichungen:

         2c2
F2'(t) = c-+-c-P'
        1   2

G'(t) = c2--c1P'
  1    c1 + c2

Durch Integration erhält man dann:

         2c2
F2(t) = c1-+c2P(t)

G1(t) = c2---c1P(t)
        c1 + c2

Also gilt:

Man definiert nun Reflexions- und Transmissionskoeffizient wie folgt:

|----------------------|
|T = --2c2--, R = c2--c1|
-----c1 +-c2----c1-+c2-|

Hierbei gilt nun darüber hinaus R + 1 = T.