- •1 Grundlagen
- •1.1 Logik und Mengen
- •1.1.1 Aussagenlogik
- •1.1.2 Mengen
- •1.2 Zahlen
- •1.2.1 Natürliche Zahlen
- •1.2.2 Ganze Zahlen
- •1.2.3 Rationale Zahlen
- •1.2.4 Reelle Zahlen
- •1.2.5 Ordnung
- •1.2.6 Intervalle
- •1.2.7 Betrag und Signum
- •1.2.8 Summe und Produkt
- •1.3 Potenz und Wurzel
- •1.3.1 Potenzen
- •1.3.2 Potenzgesetze
- •1.3.3 Wurzeln
- •1.3.4 Binomischer Satz
- •1.4 Trigonometrie
- •1.4.1 Trigonometrie im rechtwinkligen Dreieck
- •1.5 Gleichungen und Ungleichungen
- •1.5.1 Lineare Gleichungen
- •1.5.2 Potenzgleichungen
- •1.5.3 Quadratische Gleichungen
- •1.5.4 Wurzelgleichungen
- •1.5.5 Ungleichungen
- •1.6 Beweise
- •1.6.1 Direkter Beweis
- •1.6.2 Indirekter Beweis
- •1.6.3 Konstruktiver Beweis
- •1.6.4 Vollständige Induktion
- •1.7 Aufgaben
- •2 Lineare Gleichungssysteme
- •2.1 Einführung
- •2.2 Gauß-Algorithmus
- •2.2.1 Äquivalenzumformungen
- •2.2.2 Vorwärtselimination
- •2.2.3 Rückwärtseinsetzen
- •2.2.4 Gaußsches Eliminationsverfahren
- •2.2.5 Rechenschema
- •2.3 Spezielle Typen linearer Gleichungssysteme
- •2.3.1 Lineare Gleichungssysteme ohne Lösung
- •2.3.2 Lineare Gleichungssysteme mit unendlich vielen Lösungen
- •2.3.3 Systeme mit redundanten Gleichungen
- •2.3.4 Unterbestimmte lineare Gleichungssysteme
- •2.3.5 Überbestimmte lineare Gleichungssysteme
- •2.3.6 Homogene lineare Gleichungssysteme
- •2.3.7 Lineare Gleichungssysteme mit Parametern
- •2.4 Numerische Verfahren
- •2.4.1 Jakobi-Iteration
- •2.4.2 Gauß-Seidel-Iteration
- •2.5 Anwendungen
- •2.5.1 Produktion
- •2.5.2 Netzwerkanalyse in der Elektrotechnik
- •2.6 Aufgaben
- •3 Vektoren
- •3.2 Vektorrechnung ohne Koordinaten
- •3.2.1 Addition und Subtraktion
- •3.2.2 Skalare Multiplikation
- •3.2.3 Skalarprodukt
- •3.2.4 Vektorprodukt
- •3.2.5 Spatprodukt
- •3.2.6 Lineare Abhängigkeit und Komponentenzerlegung
- •3.3 Vektoren in Koordinatendarstellung
- •3.3.1 Koordinatendarstellung
- •3.3.2 Addition und Subtraktion
- •3.3.3 Skalare Multiplikation
- •3.3.4 Skalarprodukt
- •3.3.5 Vektorprodukt
- •3.3.6 Spatprodukt
- •3.3.7 Lineare Abhängigkeit und Komponentenzerlegung
- •3.4 Punkte, Geraden und Ebenen
- •3.4.1 Kartesisches Koordinatensystem
- •3.4.2 Parameterdarstellung von Geraden und Ebenen
- •3.4.3 Parameterfreie Darstellung von Geraden und Ebenen
- •3.4.4 Schnitte von Geraden und Ebenen
- •3.4.5 Abstände
- •3.4.6 Winkel
- •3.5 Anwendungen
- •3.5.1 Kraft
- •3.5.2 Arbeit
- •3.5.3 Drehmoment
- •3.6 Aufgaben
- •4 Matrizen
- •4.2 Rechnen mit Matrizen
- •4.2.1 Addition, Subtraktion und skalare Multiplikation
- •4.2.2 Multiplikation von Matrizen
- •4.3 Determinanten
- •4.3.1 Determinante einer (2,2)-Matrix
- •4.3.2 Determinante einer (3,3)-Matrix
- •4.3.3 Determinante einer (n,n)-Matrix
- •4.4 Inverse Matrix
- •4.4.1 Invertierbare Matrizen
- •4.4.2 Inverse einer (2,2)-Matrix
- •4.4.3 Inverse Matrix und lineares Gleichungssystem
- •4.5 Lineare Abbildungen
- •4.5.1 Matrizen als Abbildungen
- •4.5.2 Kern, Bild und Rang
- •4.6 Eigenwerte und Eigenvektoren
- •4.7 Numerische Verfahren
- •4.7.1 Potenzmethode
- •4.8 Anwendungen
- •4.9 Aufgaben
- •5 Funktionen
- •5.1 Einführung
- •5.1.2 Wertetabelle
- •5.1.3 Schaubild
- •5.1.4 Explizite und implizite Darstellung
- •5.1.6 Funktionsschar
- •5.1.7 Verkettung von Funktionen
- •5.2 Polynome und rationale Funktionen
- •5.2.1 Potenzfunktionen mit ganzen Hochzahlen
- •5.2.2 Polynome
- •5.2.3 Gebrochenrationale Funktionen
- •5.3 Eigenschaften
- •5.3.1 Symmetrie
- •5.3.2 Periode
- •5.3.3 Monotonie
- •5.3.4 Beschränktheit
- •5.4 Sinus, Kosinus und Tangens
- •5.4.2 Eigenschaften
- •5.5 Grenzwert und Stetigkeit
- •5.5.1 Zahlenfolgen
- •5.5.2 Grenzwert einer Funktion
- •5.5.3 Stetigkeit
- •5.5.4 Asymptotisches Verhalten
- •5.6.1 Exponentialfunktionen
- •5.6.2 Die e-Funktion
- •5.6.3 Hyperbelfunktionen
- •5.7 Umkehrfunktionen
- •5.7.1 Das Prinzip der Umkehrfunktion
- •5.7.2 Wurzelfunktionen
- •5.7.3 Arkusfunktionen
- •5.7.4 Logarithmusfunktionen
- •5.7.5 Area-Funktionen
- •5.8 Numerische Verfahren
- •5.8.1 Berechnung von Funktionswerten
- •5.8.2 Bisektionsverfahren
- •5.9 Anwendungen
- •5.9.1 Messwerte
- •5.9.2 Industrieroboter
- •5.10 Aufgaben
- •6.1 Steigung und Ableitungsfunktion
- •6.1.3 Ableitungsfunktion
- •6.1.5 Höhere Ableitungen
- •6.2 Ableitungstechnik
- •6.2.1 Ableitungsregeln
- •6.2.2 Ableitung der Umkehrfunktion
- •6.2.5 Zusammenfassung
- •6.3 Regel von Bernoulli-de l’Hospital
- •6.4 Geometrische Bedeutung der Ableitungen
- •6.4.1 Neigungswinkel und Schnittwinkel
- •6.4.2 Monotonie
- •6.4.3 Krümmung
- •6.4.4 Lokale Extrema
- •6.4.5 Wendepunkte
- •6.4.6 Globale Extrema
- •6.5 Numerische Verfahren
- •6.5.2 Newton-Verfahren
- •6.5.3 Sekantenverfahren
- •6.6 Anwendungen
- •6.6.1 Fehlerrechnung
- •6.6.2 Extremwertaufgaben
- •6.7 Aufgaben
- •7 Integralrechnung
- •7.1 Flächenproblem
- •7.1.1 Integralsymbol
- •7.1.2 Integral als Grenzwert von Summen
- •7.1.3 Bestimmtes Integral
- •7.2 Zusammenhang von Ableitung und Integral
- •7.2.1 Integralfunktion
- •7.2.2 Stammfunktion
- •7.2.3 Bestimmtes Integral und Stammfunktion
- •7.2.4 Mittelwertsatz der Integralrechnung
- •7.3 Integrationstechnik
- •7.3.1 Integrationsregeln
- •7.3.2 Integration durch Substitution
- •7.3.3 Partielle Integration
- •7.3.4 Gebrochenrationale Funktionen
- •7.3.5 Uneigentliche Integrale
- •7.4 Länge, Flächeninhalt und Volumen
- •7.4.1 Flächeninhalte
- •7.4.2 Bogenlänge
- •7.4.3 Rotationskörper
- •7.5 Numerische Verfahren
- •7.5.1 Trapezregel
- •7.5.2 Romberg-Verfahren
- •7.6 Anwendungen
- •7.6.2 Schwerpunkte und statische Momente ebener Flächen
- •7.7 Aufgaben
- •8 Potenzreihen
- •8.1 Unendliche Reihen
- •8.2 Potenzreihen und Konvergenz
- •8.3 Taylor-Reihen
- •8.4 Eigenschaften
- •8.5 Numerische Verfahren
- •8.5.1 Berechnung von Funktionswerten
- •8.6 Anwendungen
- •8.6.1 Normalverteilung in der Statistik
- •8.7 Aufgaben
- •9 Kurven
- •9.1 Parameterdarstellung
- •9.2 Kegelschnitte
- •9.3 Tangente
- •9.4 Krümmung
- •9.5 Bogenlänge
- •9.6 Numerische Verfahren
- •9.6.1 Bézier-Kurve
- •9.7 Anwendungen
- •9.7.1 Mechanik
- •9.7.2 Straßenbau
- •9.8 Aufgaben
- •10 Funktionen mit mehreren Variablen
- •10.1.2 Schaubild einer Funktion mit mehreren Variablen
- •10.1.3 Schnittkurven mit Ebenen und Höhenlinien
- •10.2 Grenzwert und Stetigkeit
- •10.2.1 Grenzwert einer Funktion mit mehreren Variablen
- •10.2.2 Stetigkeit
- •10.3.3 Gradient und Richtungsableitung
- •10.3.5 Höhere partielle Ableitungen
- •10.3.6 Extremwerte
- •10.4 Ausgleichsrechnung
- •10.4.1 Methode der kleinsten Fehlerquadrate
- •10.4.2 Ausgleichsrechnung mit Polynomen
- •10.4.3 Lineare Ausgleichsrechnung
- •10.5 Vektorwertige Funktionen
- •10.6 Numerische Verfahren
- •10.6.1 Mehrdimensionales Newton-Verfahren
- •10.6.2 Gradientenverfahren
- •10.7 Anwendungen
- •10.7.1 Fehlerrechnung
- •10.8 Aufgaben
- •11 Komplexe Zahlen und Funktionen
- •11.1.1 Komplexe Zahlen
- •11.1.2 Gaußsche Zahlenebene
- •11.1.3 Polarkoordinaten
- •11.1.4 Exponentialform
- •11.2 Rechenregeln
- •11.2.1 Gleichheit
- •11.2.2 Addition und Subtraktion
- •11.2.3 Multiplikation und Division
- •11.2.4 Rechnen mit der konjugiert komplexen Zahl
- •11.2.5 Rechnen mit dem Betrag einer komplexen Zahl
- •11.3 Potenzen, Wurzeln und Polynome
- •11.3.1 Potenzen
- •11.3.2 Wurzeln
- •11.3.3 Fundamentalsatz der Algebra
- •11.4 Komplexe Funktionen
- •11.4.1 Ortskurven
- •11.4.2 Harmonische Schwingungen
- •11.4.3 Transformationen
- •11.5 Anwendungen
- •11.5.1 Komplexe Wechselstromrechnung
- •11.6 Aufgaben
- •12.1 Einführung
- •12.1.3 Richtungsfeld und Orthogonaltrajektorie
- •12.2.1 Separation der Variablen
- •12.2.2 Lineare Substitution
- •12.3.3 Allgemeine Eigenschaften
- •12.4.1 Allgemeine Form
- •12.4.2 Freie Schwingung
- •12.4.3 Harmonisch angeregte Schwingung
- •12.4.4 Frequenzgänge
- •12.5.1 Eliminationsverfahren
- •12.5.2 Zustandsvariablen
- •12.5.5 Stabilität
- •12.6 Numerische Verfahren
- •12.6.1 Polygonzugverfahren von Euler
- •12.7 Anwendungen
- •12.7.1 Temperaturverlauf
- •12.7.2 Radioaktiver Zerfall
- •12.7.3 Freier Fall mit Luftwiderstand
- •12.7.4 Feder-Masse-Schwinger
- •12.7.5 Pendel
- •12.7.6 Wechselstromkreise
- •12.8 Aufgaben
- •13 Fourier-Reihen
- •13.1 Fourier-Analyse
- •13.1.1 Periodische Funktionen
- •13.1.2 Trigonometrische Polynome
- •13.1.3 Fourier-Reihe
- •13.1.4 Satz von Fourier
- •13.1.5 Gibbssches Phänomen
- •13.2 Komplexe Darstellung
- •13.2.1 Komplexe Fourier-Reihe
- •13.2.3 Spektrum
- •13.2.4 Minimaleigenschaft
- •13.3 Eigenschaften
- •13.3.1 Symmetrie
- •13.3.2 Integrationsintervall
- •13.3.3 Mittelwert
- •13.3.4 Linearität
- •13.3.5 Ähnlichkeit und Zeitumkehr
- •13.3.6 Zeitverschiebung
- •13.4 Aufgaben
- •14 Verallgemeinerte Funktionen
- •14.1 Heaviside-Funktion
- •14.2 Dirac-Distribution
- •14.3 Verallgemeinerte Ableitung
- •14.4 Faltung
- •14.5 Aufgaben
- •15 Fourier-Transformation
- •15.1 Integraltransformation
- •15.1.4 Transformation gerader und ungerader Funktionen
- •15.1.5 Darstellung mit Amplitude und Phase
- •15.2 Eigenschaften
- •15.2.1 Linearität
- •15.2.2 Zeitverschiebung
- •15.2.3 Amplitudenmodulation
- •15.2.4 Ähnlichkeit und Zeitumkehr
- •15.3 Inverse Fourier-Transformation
- •15.3.2 Vertauschungssatz
- •15.3.3 Linearität
- •15.4.3 Multiplikationssatz
- •15.4.5 Faltung
- •15.5 Periodische Funktionen
- •15.5.1 Fourier-Transformation einer Fourier-Reihe
- •15.5.3 Grenzwertbetrachtung
- •15.6 Anwendungen
- •15.6.1 Lineare zeitinvariante Systeme
- •15.7 Aufgaben
- •16 Laplace-Transformation
- •16.1 Bildbereich
- •16.2 Eigenschaften
- •16.2.1 Linearität
- •16.2.2 Ähnlichkeit
- •16.2.3 Zeitverschiebung
- •16.2.4 Dämpfung
- •16.3.2 Integration
- •16.3.3 Faltung
- •16.3.4 Grenzwerte
- •16.4 Transformation periodischer Funktionen
- •16.5 Rücktransformation
- •16.7 Anwendungen
- •16.7.1 Regelungstechnik
- •16.8 Aufgaben
- •17 z-Transformation
- •17.1 Transformation diskreter Signale
- •17.1.2 z-Transformation und Laplace-Transformation
- •17.2 Eigenschaften
- •17.2.1 Linearität
- •17.2.2 Verschiebung
- •17.2.3 Dämpfung
- •17.4 Anwendungen
- •17.4.1 Zeitkomplexität von Quicksort
- •A Anhang
- •A.1 Ableitungsregeln
- •A.2 Ableitungen
- •A.3 Potenzreihen
- •A.4 Integralregeln
- •A.5 Integrale
- •A.6 Fourier-Reihen
- •A.7 Fourier-Transformationen
- •A.8 Laplace-Transformationen
- •A.9 Griechisches Alphabet
- •A.10 Bedeutende Mathematiker
- •Literaturverzeichnis
- •Sachwortverzeichnis
12.3 Lineare Di erenzialgleichungen |
475 |
12.3.3 Allgemeine Eigenschaften
Lineare Di erenzialgleichungen besitzen besonders schöne Eigenschaften. In diesem Abschnitt zeigen wir, dass bei homogenen Di erenzialgleichungen Linearkombinationen von Lösungen wieder Lösungen ergeben. Außerdem werden wir sehen, dass die allgemeine Lösung einer homogenen Di erenzialgleichung aus elementaren Lösungsbausteinen, den sogenannten Fundamentallösungen, zusammengesetzt wird. Schließlich stellen wir das sogenannte Superpositionsprinzip für inhomogene Di erenzialgleichungen vor.
Kennt man bereits zwei Lösungen y1 und y2 einer homogenen linearen Di erenzialgleichung, so kann man diese beiden Lösungen mit beliebigen Konstanten multiplizieren und addieren und erhält wieder eine Lösung.
Satz 12.7 (Linearkombination von Lösungen)
Sind y1 und y2 Lösungen einer homogenen linearen Di erenzialgleichung, dann ist auch
y(x) = C1 y1(x) + C2 y2(x)
mit beliebigen Konstanten C1 und C2 eine Lösung dieser Di erenzialgleichung.
Zum Nachweis dieser Eigenschaft setzt man y(x) = C1 y1(x) + C2 y2(x) in die Di erenzialgleichung aus Definition 12.11 ein:
an(x) ŠC1 y1(n) + C2 y2(n)•+. . .+a1(x) (C1 y1′ + C2 y2′ )+a0(x) (C1 y1 + C2 y2) = 0.
Fasst man die Terme mit y1 und mit y2 jeweils zusammen, so entsteht
C1 Šan(x)y1(n) + . . . + a0(x)y1• + C2 Šan(x)y2(n) + . . . + a0(x)y2• = 0.
´¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¸¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¶ |
´¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¸¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¶ |
0 |
0 |
Die beiden Ausdrücke in den Klammern sind jeweils null, denn y1 und y2 sind ja Lösungen der homogenen Gleichung. Also löst y als Linearkombination von y1 und y2 ebenfalls die homogene Di erenzialgleichung.
Beispiel 12.21 (Linearkombination von Lösungen)
Die homogene Di erenzialgleichung |
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
y′′ 9 y |
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
hat die beiden+ |
Lösungen= |
y1 |
( |
x |
) = |
sin |
( |
3 x |
) |
und y2 |
x |
cos |
( |
3 x . Also ist auch |
|||||||||||
y |
( |
x |
) = |
C |
sin 3 |
x |
|
|
|
( |
|
|
|
|
( ) = |
|
|
) |
|||||||
|
|
|
|
C cos |
3 x |
|
|
|
|
|
|
|
|
||||||||||||
|
|
1 |
|
( ) + |
|
|
2 |
|
|
|
) |
|
|
C |
1 und |
C |
2. |
Ì |
|||||||
eine Lösung. Dies gilt für beliebige Konstanten |
|
|
Die allgemeine Lösung einer homogenen linearen Di erenzialgleichung ist ähnlich strukturiert wie die Lösungsmenge eines homogenen linearen Gleichungssystems. An die Stelle der Basisvektoren treten bei Di erenzialgleichungen die sogenannten Fundamentallösungen. Ein Fundamentalsystem setzt sich dann aus Fundamentallösungen zusammen.
476 |
12 Gewöhnliche Di erenzialgleichungen |
Definition 12.15 (Fundamentalsystem)
Die n Lösungen y1, y2, . . ., yn einer linearen Di erenzialgleichung n-ter Ordnung bilden ein Fundamentalsystem, wenn aus
C1 y1(x) + C2 y2(x) + . . . + Cn yn(x) = 0 für alle x R
folgt, dass C1 = 0, C2 = 0, . . ., Cn = 0. Man nennt dann y1, y2, . . ., yn Fundamentallösungen der Di erenzialgleichung.
Bei linearen Di erenzialgleichungen besteht der zentrale Sachverhalt darin, die allgemeine Lösung mithilfe von Fundamentallösungen zu ermitteln. Auf die theoretischen Details gehen wir jedoch nicht ein.
Satz 12.8 (Allgemeine Lösung und Fundamentallösung)
Die allgemeine Lösung einer homogenen linearen Di erenzialgleichung n-ter Ordnung erhält man als Linearkombination von n Fundamentallösungen
y(x) = C1 y1(x) + C2 y2(x) + . . . + Cn yn(x).
Dabei bilden die n Fundamentallösungen y1, y2, . . ., yn ein Fundamentalsystem und C1, C2, . . ., Cn sind beliebige Konstanten.
Es stellt sich nun die Frage, wann n Funktionen ein Fundamentalsystem bilden. Diese Fragestellung ist ähnlich zu der Frage nach n linear unabhängigen Vektoren und wird deshalb auch mithilfe einer Determinanten beantwortet. Diese Determinante ist nach dem polnischen Mathematiker und Philosophen Joseph Marie Wronski benannt.
Satz 12.9 (Wronski-Kriterium)
Die n Funktionen y1, y2, . . ., yn bilden genau dann ein Fundamentalsystem, wenn die Wronski-Determinante für alle reellen Zahlen x nicht null ist, also
R |
|
y1 |
( |
x |
) |
|
y2 |
( |
x |
) |
. . . |
yn |
( |
x |
) |
|
R |
0 |
für alle x R. |
|||||||
R |
|
y1 |
x |
|
y2 |
x |
. . . |
yn |
x |
|
R |
|||||||||||||||
R |
|
|
′ |
|
|
|
|
|
|
|
′ |
|
|
|
|
|
′ |
|
|
|
|
|
|
R |
|
|
R |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
R |
|
|
R |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
R |
|
|
R |
|
|
|
( |
|
|
) |
|
|
|
( |
|
|
) |
|
|
( |
|
|
) |
|
R |
|
|
||
R |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
R |
|
|
||||||||
R |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
R |
|
|
R |
y |
1 |
|
|
|
|
|
x |
y |
2 |
|
|
|
|
x |
. . . yn |
|
|
|
|
x |
|
R |
|
|
|
R |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
R ≠ |
|
|||
R |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
R |
|
|
R |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
R |
|
|
|||
R |
|
( − ) |
|
|
|
|
( − ) |
|
|
|
( − ) |
|
|
|
|
R |
|
|
||||||||
R |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
R |
|
|
R |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
R |
|
|
R |
|
n 1 |
|
|
|
|
|
n 1 |
|
|
|
|
n 1 |
|
|
|
|
|
R |
|
|
|||||
R |
|
|
|
|
|
( |
|
) |
|
|
|
|
( ) |
|
|
|
|
( |
|
) |
R |
|
|
|||
R |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
R |
|
|
R |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
R |
|
|
Beispiel 12.22 (Wronski-Kriterium)
a)Wir untersuchen, ob die beiden Funktionen y1(x) = cos (ω x) und y2(x) = sin (ω x) für ω ≠ 0 ein Fundamentalsystem bilden. Dazu berechnen wir die Wronski-Determinante
|
|
cos |
ω x |
sin |
ω x |
|
ω cos2 |
|
ωx |
|
ω sin2 |
|
ωx |
|
ω. |
W |
− |
ω sin |
(ω x) |
ω cos |
(ω x) |
W = |
( |
) + |
( |
) = |
|||||
|
( ) |
|
( ) |
|
|
|
|
|
Die Wronski-Determinante ist für alle reellen Zahlen x ungleich null. Somit bilden die beiden Funktionen ein Fundamentalsystem.
12.3 Lineare Di erenzialgleichungen |
477 |
b)Die drei Polynome y1(x) = 1, y2(x) = x und y3(x) = x2 bilden ein Fundamentalsystem. Denn die Wronski-Determinante
R |
0 |
1 |
2 x R |
|
2 |
|
R |
1 |
x |
x2 |
R |
|
|
R |
0 |
0 |
2 |
R |
|
|
R |
R |
|
|
|||
R |
|
|
|
R |
|
|
R |
|
|
|
R |
= |
|
R |
|
|
|
R |
|
|
R |
|
|
|
R |
|
|
R |
|
|
|
R |
|
|
R |
|
|
|
R |
|
|
R |
|
|
|
R |
|
|
R |
|
|
|
R |
|
|
R |
|
|
|
R |
|
|
R |
|
|
|
R |
|
|
ist für alle reellen Zahlen x ungleich null.
c)Die beiden Funktionen y1(x) = ea x und y2(x) = eb x bilden für a ≠ b ein Fundamentalsystem. Die Wronski-Determinante
W |
ea x |
eb x |
W = (b − a)e( |
a |
+ |
b |
) |
x |
|
aea x |
beb x |
|
|
|
|
||||
ist für alle reellen Zahlen x ungleich null, falls a und b verschieden sind. |
Ì |
Eine weitere Eigenschaft linearer Di erenzialgleichungen betri t die Störfunktionen. Ist eine Störfunktion die Summe von Teilstörfunktionen, so kann zunächst für jede Teilstörfunktion separat eine partikuläre Lösung bestimmt werden. Diese einzelnen partikulären Teillösungen überlagern sich bei linearen Di erenzialgleichungen zur gesuchten partikulären Lösung.
Satz 12.10 (Superposition)
Wenn bei einer linearen Di erenzialgleichung die Störfunktion r durch Addition oder Subtraktion einzelner Störfunktionen zusammengesetzt ist
r(x) = r1(x) ± r2(x) ± r3(x) ± . . . ,
kann man die partikuläre Lösung für jede einzelne Störfunktion ri getrennt ermitteln.
Das Prinzip der Superposition basiert auf der Linearität der hier betrachteten Di erenzialgleichungen. Ist y1 eine partikuläre Lösung für die Störfunktion r1 und y2 eine partikuläre Lösung für r2, so gilt für die Summe y = y1 + y2 durch Einsetzen in die Di erenzialgleichung aus Definition 12.11
an(x) Šy1(n) + y2(n)• + . . . + a1(x) (y1′ + y2′ ) + a0(x) (y1 + y2) = r(x).
Durch entsprechendes Zusammenfassen
Šan(x) y1(n) + . . . + a0(x) y1• +Šan(x) y2(n) + . . . + a0(x) y2• = r(x)
´¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¸¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¶ ´¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¸¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¶ |
|
||||||||||
|
r1 |
x |
|
|
|
|
|
r2 |
x |
|
|
erkennt man, dass |
y |
( |
y |
|
y |
|
eine partikuläre |
Lösung für die Störfunktion r |
r ist. |
||
|
|
) |
|
2 |
|
( ) |
1 + 2 |
||||
|
|
= |
|
|
1 + |
|
|
|
|
478 |
12 Gewöhnliche Di erenzialgleichungen |
Beispiel 12.23 (Superposition)
Wir betrachten eine Di erenzialgleichung, bei der die Störfunktion durch Addition dreier Störfunktionen zusammengesetzt ist:
y′ |
+ |
2y |
= |
5 e3 x |
4 x |
+ |
4 cos |
|
2 x |
) |
. |
||||
|
|
|
1 |
|
|
+r2 |
x |
r3 ( |
|
|
|||||
|
|
|
|
² |
° |
´¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¸¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¶ |
|
||||||||
|
|
|
|
r |
|
( |
x |
|
( |
) |
( |
x |
|
|
|
|
|
|
|
|
|
) |
|
|
) |
|
|
Die Di erenzialgleichung hat die homogene Lösung yh(x) = C e−2 x y1(x) = e−2 x. Mit der Formel aus Satz 12.6 können wir zu jeder partikuläre Lösung berechnen:
und die Fundamentallösung einzelnen Störfunktion eine
|
( |
|
) = |
e−2 x |
5 e3 x |
|
= |
e−2 x |
5 e5 x dx |
|
= |
e3 x |
|
|
||||||
yp1 |
x |
|
|
dx |
|
|
|
|||||||||||||
e 2 x |
|
|
|
|||||||||||||||||
|
|
|
2 x |
4−x |
|
|
2 x |
|
2 x |
|
|
|
|
|
||||||
yp2 |
( |
x |
) = |
e− |
S |
|
|
dx |
= |
e− |
S |
4 x e |
|
dx |
|
= |
2 x |
− |
1 |
|
e 2 x |
|
|
||||||||||||||||||
|
|
|
2 x S |
4−cos |
2 x |
|
2 x S |
|
|
|
2 x |
|
|
|||||||
yp3 |
(x) = |
e− |
S |
e−2(x |
) |
dx = e− |
S |
4 cos (2 x) e |
|
dx = |
cos (2 x) + sin (2 x) |
Dabei haben wir zur Berechnung der Stammfunktionen Formeln aus Anhang A.5 verwendet. Alles in allem ergibt sich die allgemeine Lösung
y(x) = yh(x) + yp1(x) + yp2(x) + yp3(x) = C e−2 x + e3 x + 2 x − 1 + cos (2 x) + sin (2 x)
für die gesamte Ausgangsdi erenzialgleichung. |
Ì |
Das Superpositionsprinzip funktioniert nicht nur für Störfunktionen, die aus endlich vielen Teilfunktionen zusammengesetzt sind, sondern auch für die Überlagerung von unendlich vielen Teilfunktionen. Man kann Superposition deshalb auch auf Funktionen, die in Form von Potenzreihen oder Fourier-Reihen dargestellt werden, anwenden.
In Beispiel 12.23 erkennt man, dass sich Störfunktionen und partikuläre Lösungen jeweils entsprechen. Wenn die Störfunktion ein Polynom vom Grad n ist, dann ist auch die zugeordnete partikuläre Lösung ein Polynom vom Grad n. Entsprechendes gilt für die Exponentialfunktion, wobei Störfunktion und partikuläre Lösung denselben Exponenten besitzen. Auch eine Störfunktion in Form eines Kosinus erzeugt eine partikuläre Lösung in Form einer harmonischen Schwingung, wobei die Kreisfrequenz erhalten bleibt. Eine Störfunktion erzeugt also eine partikuläre Lösung vom selben Typ. Dieses Prinzip werden wir in Abschnitt 12.3.4 verwenden, um partikuläre Lösungen mithilfe geeigneter Störansätze zu ermitteln.
12.3.4 Di erenzialgleichungen mit konstanten Koe zienten
Lineare Di erenzialgleichungen, bei denen die Koe zienten keine echte Funktionen ai(x), sondern nur konstante Zahlen ai sind, lassen sich systematisch lösen.
12.3 Lineare Di erenzialgleichungen |
479 |
Definition 12.16 (Lineare Di erenzialgleichung mit konstanten Koe zienten)
Eine lineare Di erenzialgleichung, bei der die Faktoren vor den Ableitungen keine echten Funktionen, sondern nur Konstanten a0, a1, . . ., an sind
an y(n) + an−1 y(n−1) + . . . + a2 y′′ + a1 y′ + a0 y = r(x),
nennt man eine lineare Di erenzialgleichung mit konstanten Koe zienten.
Entsprechend unserer Lösungsstrategie für lineare Di erenzialgleichungen betrachten wir zunächst nur homogene Di erenzialgleichungen. Die allgemeine Lösung einer homogenen linearen Di erenzialgleichung erster Ordnung mit konstanten Koe zienten besteht aus Exponentialfunktionen. Dies legt die Vermutung nahe, dass die allgemeine Lösung auch bei linearen Di erenzialgleichungen mit konstanten Koe zienten höherer Ordnung aus Exponentialfunktionen bestehen.
Exponentialansatz
Bei linearen homogenen Di erenzialgleichungen mit konstanten Koe zienten versucht man, ein Fundamentalsystem aus Exponentialfunktionen zu bestimmen:
y(x) = eλx
Wenn wir eine homogene lineare Di erenzialgleichung der Ordnung n mit konstanten Koe zienten mit dem Exponentialansatz y(x) = eλ x versuchen zu lösen, dann erhalten wir aus der homogenen Gleichung in Definition 12.16
an λneλx |
+ |
an |
− |
1 λn−1eλx |
+ |
. . . |
+ |
a2 λ2eλx |
+ |
a1 λeλx |
+ |
a0 eλx |
= |
0. |
||||||
y |
n |
) |
|
y |
n |
1 |
) |
|
y |
′′ |
y |
y |
|
|||||||
( |
|
|
|
|
( |
|
− |
|
|
|
|
|
±′ |
|
° |
|
|
|||
´¹¹¹¹¸¹¹¹¹¶ |
|
|
|
´¹¹¹¹¹¹¹¹¹¹¹¹¸¹¹¹¹¹¹¹¹¹¹¹¹¶ |
|
|
|
² |
|
|
|
|
|
|
Diese Gleichung können wir durch eλx teilen. Damit haben wir x aus der Gleichung eliminiert:
anλn + an−1λn−1 + . . . + a2λ2 + a1λ + a0 = 0.
Definition 12.17 (Charakteristische Gleichung)
Zur homogenen linearen Di erenzialgleichung
an y(n) + an−1 y(n−1) + . . . + a2 y′′ + a1 y′ + a0 y = 0,
gehört die charakteristische Gleichung
anλn + an−1λn−1 + . . . + a2λ2 + a1λ + a0 = 0.
Sie entsteht aus der Di erenzialgleichung durch den Exponentialansatz y(x) = eλx.
480 |
12 Gewöhnliche Di erenzialgleichungen |
Die charakteristische Gleichung entsteht also aus der Di erenzialgleichung dadurch, dass man die Ableitungen y(k) der gesuchten Funktion y durch Potenzen λk ersetzt. Dabei ist zu beachten, dass insbesondere y′ der Variablen λ1 = λ und y dem Wert λ0 = 1 entspricht. Hat die Di erenzialgleichung die Ordnung n, so besteht die charakteristische Gleichung auf der linken Seite aus einem Polynom vom Grad n mit der Variablen λ. Dieses Polynom nennt man auch charakteristisches Polynom.
Die Lösung einer linearen Di erenzialgleichung mithilfe der charakteristischen Gleichung stellt eine gewaltige Vereinfachung dar. Das Problem wird durch den Exponentialansatz soweit reduziert, dass man nur noch die Nullstellen eines Polynoms berechnen muss.
Beispiel 12.24 (Charakteristische Gleichung)
Zu der homogenen linearen Di erenzialgleichung betrachten wir die charakteristische Gleichung:
y(5) − 3 y′′′ + 52 y′′ = 0 Ô λ5 − 3 λ3 + 52 λ2 = 0 .
Den Faktor λ2 kann man ausklammern. Dadurch sind λ1,2 = 0 Lösungen der charakteristischen Gleichung. Außerdem ist λ3 = −4 eine Lösung. Eine Polynomdivision mit λ + 4 ergibt
λ5 − 3 λ3 + 52 λ2 = λ2 (λ + 4) ‰λ2 − 4 λ + 13Ž .
Die beiden letzten Lösungen erhält man aus der quadratischen Gleichung: |
|
||||||||||
λ4,5 |
|
4 |
|
√ |
|
|
|
2 |
|
3 i . |
|
|
|
16 |
52 |
|
|
|
|||||
|
= |
|
± |
|
2 − |
|
= |
|
± |
|
Ì |
Definition 12.18 (Eigenwerte und Eigenfunktionen)
Die Lösungen der charakteristischen Gleichung
anλn + an−1λn−1 + . . . + a2λ2 + a1λ + a0 = 0
bezeichnet man als Eigenwerte λ1, λ2, . . ., λn. Die Lösungsfunktionen, die sich daraus durch den Exponentialansatz ergeben, nennt man Eigenfunktionen
y1(x) = eλ1x, y2(x) = eλ2x, . . . , yn(x) = eλnx.
Eigenwerte und Eigenfunktionen können reell oder komplex sein.
Die Bezeichnung Eigenwert wird klar, wenn wir den Zusammenhang zwischen einer linearen Di erenzialgleichung und einem System von Di erenzialgleichungen erster Ordnung betrachten, siehe Abschnitt 12.5.3. Lineare Systeme können in Matrixschreibweise formuliert werden. Die Eigenwerte der Systemmatrix, siehe Definition 4.19, entsprechen dann den Nullstellen des charakteristischen Polynoms.
Wenn alle Eigenwerte reell und paarweise verschieden sind, dann sind die Eigenfunktionen genau unsere gesuchten Fundamentallösungen. In Beispiel 12.24 haben wir jedoch bereits gesehen, dass die charakteristische Gleichung unter Umständen auch Paare konjugiert komplexer Zahlen als Lösungen haben kann. Außerdem können Lösungen, also Nullstellen des charakteristischen Polynoms, mehrfach vorkommen.
12.3 Lineare Di erenzialgleichungen |
481 |
Fälle bei Eigenwerten
Bei den Eigenwerten einer linearen Di erenzialgleichung mit konstanten Koe zienten sind folgende Fälle zu beachten:
(A) |
Einfache reelle Eigenwerte |
(B) |
Mehrfache reelle Eigenwerte |
(C) |
Einfache komplexe Eigenwerte |
(D) |
Mehrfache komplexe Eigenwerte |
Komplexe Eigenwerte und Eigenfunktionen sind natürlich nicht unser Ziel, sondern nur Mittel zum Zweck. Wir werden Methoden aufzeigen, mit denen wir aus komplexen Eigenwerten und Eigenfunktionen reelle Lösungen bestimmen.
(A) Einfache reelle Eigenwerte
Jeder einfache reelle Eigenwert λ einer linearen Di erenzialgleichung mit konstanten Koe zienten erzeugt genau eine reelle Fundamentallösung
y(x) = eλ x.
Beispiel 12.25 (Einfache reelle Eigenwerte)
a) Die Di erenzialgleichung mit der charakteristischen Gleichung
y′′ − 7y′ + 12y = 0 Ô λ2 − 7λ + 12 = 0
besitzt die Eigenwerte λ1 = 3 und λ2 = 4. Daraus ergeben sich die Eigenfunktionen
y1(x) = e3 x, y2(x) = e4 x.
In Beispiel 12.22 haben wir bereits gesehen, dass zwei Exponentialfunktionen mit unterschiedlichem Exponenten ein Fundamentalsystem bilden. Die allgemeine Lösung der Di e- renzialgleichung ist
y(x) = C1e3 x + C2e4 x .
b) Die charakteristische Gleichung der Di erenzialgleichung
y′′ + 3 y′ = 0 Ô λ2 + 3 λ = 0
liefert die Eigenwerte λ1 = 0 und λ2 = −3. Die entsprechenden Eigenfunktionen sind
y1(x) = e0 x = 1, y2(x) = e−3 x.
Exponentialfunktionen mit unterschiedlichen Exponenten bilden ein Fundamentalsystem, siehe Beispiel 12.22. Somit ist
y(x) = C1 + C2e−3 x
die allgemeine Lösung der Di erenzialgleichung. |
Ì |
482 |
12 Gewöhnliche Di erenzialgleichungen |
Bei linearen Di erenzialgleichungen mit konstanten Koe zienten stellen die Eigenwerte eine direkte Verbindung zwischen der allgemeinen Lösung der homogenen Gleichung und der Di erenzialgleichung selbst her. Man kann dadurch Di erenzialgleichungen konstruieren, die vorgegebene Lösungen haben. Dazu geht man in umgekehrter Reihenfolge wie bei der Lösung einer Di erenzialgleichung vor. Aus der allgemeinen Lösung bestimmt man die Fundamentallösungen und daraus die Eigenwerte. Mit den Eigenwerten bildet man die charakteristische Gleichung, woraus sich die Di erenzialgleichung ergibt.
Beispiel 12.26 (Konstruktion einer Di erenzialgleichung aus Eigenwerten)
Wir suchen eine Di erenzialgleichung mit der allgemeinen Lösung
y(x) = C1e−2 x + C2e5 x.
Die Fundamentallösungen sind y1(x) = e−2 x und y2(x) = e5 x. Die Eigenwerte müssen λ1 = −2 und λ2 = 5. sein. Deshalb können wir die charakteristische Gleichung in der Form
(λ + 2) (λ − 5) = λ2 − 3 λ − 10
angeben. Daraus kann man die gesuchte Di erenzialgleichung ablesen:
y′′ |
− |
3 y′ |
− |
10 y |
= |
0. |
Ì |
|
|
|
|
|
Bei mehrfachen reellen Eigenwerten behilft man sich damit, dass man die erste Fundamentallösung entsprechend oft mit x multipliziert. Dieser Trick funktioniert natürlich nur unter den Voraussetzungen, dass das Wronski-Kriterium erfüllt ist und dass durch die Multiplikation mit x wieder eine Lösung der homogenen Di erenzialgleichung entsteht. Im Fall eines doppelten reellen Eigenwerts werden wir diese beiden Bedingungen nachrechnen. Bei einem n-fachen reellen Eigenwert verläuft der Nachweis analog.
(B) Mehrfache reelle Eigenwerte
Jeder doppelte reelle Eigenwert λ einer linearen Di erenzialgleichung mit konstanten Koe zienten erzeugt genau zwei reelle Fundamentallösungen
y1(x) = eλ x, y2(x) = x eλ x.
Falls die Vielfachheit des reellen Eigenwerts größer als zwei ist, entstehen jeweils durch Multiplikation mit x neue Fundamentallösungen
y3(x) = x2 eλ x, y4(x) = x3 eλ x, . . .
Das Wronski-Kriterium für die beiden Funktionen y |
1( |
x |
) = |
eλ x |
und y |
2 |
( |
x |
) = |
x eλ x |
ist er- |
|||
füllt, denn die Wronski-Determinante |
|
|
|
|
|
|
|
|
||||||
eλ x |
x eλ x |
1 |
x |
W = e2 λ x |
|
|
|
|
|
|
|
|||
W λ eλ x |
(1 + λ x) eλ x W = e2 λ x W |
λ |
1 + λ x |
|
|
|
|
|
|
|
ist für alle reellen Zahlen x ungleich null.
12.3 Lineare Di erenzialgleichungen |
483 |
Die Funktion y2 ist tatsächlich auch eine Lösung. Zum Nachweis betrachten wir die homogene lineare Di erenzialgleichung zweiter Ordnung und ihre charakteristische Gleichung:
a2 y′′ + a1 y′ + a0 y = 0, Ô a2λ2 + a1λ + a0 = 0.
Die charakteristische Gleichung hat genau dann eine doppelte reelle Lösung, wenn auch die Ableitung nach λ null ergibt, also 2 a2 λ + a1 = 0. Wir berechnen die erste und die zweite Ableitung von y2 und setzen alles in die Di erenzialgleichung ein:
a |
λ λ2 x |
Ž e |
λ x |
+ |
a 1 λ x |
) |
eλ x |
+ |
a |
x eλ x |
= |
0. |
|||||||||
2 |
‰2 y |
+ x |
|
|
1 |
(y+ x |
|
|
|
0 y2 x |
|
||||||||||
|
´¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¸¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¶ |
|
|
|
|
|
´¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¸¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¶ |
|
|
|
|
² |
|
|
|||||||
|
|
′′ |
|
|
|
|
|
|
|
|
′ |
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
Diese Gleichung( können) |
|
|
|
|
|
x |
|
|
|
|
|
|
zusammenfassen: |
||||||||
wir durch e(λ ) teilen und anders( ) |
|||||||||||||||||||||
x ‰a2λ |
2 |
a λ |
+ |
a |
0Ž + 2 |
a |
λ |
a |
= 0 |
. |
|
|
|
|
|
||||||
|
+ 01 |
|
|
|
|
2 0 + |
1 |
|
|
|
|
|
|
||||||||
´¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¸¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¶ |
´¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¸¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¶ |
|
|
|
|
|
|
|
|
Somit ist auch y2 eine Lösung der linearen homogenen Di erenzialgleichung.
Beispiel 12.27 (Mehrfache reelle Eigenwerte)
a) Die Lösung der Di erenzialgleichung bestimmen wir aus der charakteristischen Gleichung:
y′′ + 2y′ + y = 0 Ô λ2 + 2λ + 1 = 0
Der doppelte reelle Eigenwert λ1,2 = −1 ergibt die erste Fundamentallösung y1(x) = e−x. Diese multiplizieren wir mit x und erhalten die zweite Fundamentallösung y2(x) = x e−x. Die Linearkombination der Fundamentallösungen liefert die allgemeine Lösung
y(x) = C1 e−x + C2 x e−x.
b) Aus der charakteristischen Gleichung der Di erenzialgleichung |
|
|||||
y′′′ |
3y′′ |
3y′ |
y 0 |
λ3 3λ2 3λ 1 |
λ 1 3 |
0 |
erhalten+ wir |
+einen |
+dreifachen= Ôreellen Eigenwert+ + λ1+,2,3= ( 1+. Somit) = |
ist |
|||
y(x) = C1 e−x + C2 x e−x + C3 x2 e−x |
= − |
|
||||
die allgemeine Lösung der Di erenzialgleichung. |
|
Ì |
(C) Einfache komplexe Eigenwerte
Jedes einfache konjugiert komplexe Paar Eigenwerte λ1,2 = a ± i b einer linearen Di e- renzialgleichung mit konstanten Koe zienten erzeugt genau zwei reelle Fundamentallösungen
y1(x) = ea x cos (b x), y2(x) = ea x sin (b x).
Der Realteil erzeugt Lösungsanteile in Form von Exponentialfunktionen und die Imaginärteile erzeugen Sinusund Kosinusanteile.
484 |
12 Gewöhnliche Di erenzialgleichungen |
Die Lösungen der charakteristischen Gleichung sind die Nullstellen eines reellen Polynoms. Aus Satz 11.14 wissen wir, dass komplexe Nullstellen immer als konjugiert komplexes Paar λ = a ± i b auftreten. Die beiden Funktionen
z1,2(x) = e(a ±i b)x = ea x‰ cos (b x) ± i sin (b x)Ž
lassen sich nach dem Satz von Euler in Realund Imaginärteil zerlegen. Sie sind zwar Lösungen der Di erenzialgleichung, allerdings sind diese beiden Funktionen komplex und nicht rein reell. Durch einen Trick kann man aus diesen beiden komplexen Lösungen zwei reelle Lösungen erzeugen: Aufgrund der Linearität wissen wir, dass die Linearkombination aus zwei Lösungen wieder eine Lösung ergibt. Deshalb sind die beiden Lösungen
y1 |
|
x |
1z1 |
|
x |
|
1 z2 |
|
x |
|
|||||
|
( |
) = |
|
|
|
( |
|
) + |
|
|
|
|
( |
|
) |
|
1 |
|
|
1 |
|
|
|||||||||
|
|
|
2 |
|
|
|
|
|
2 |
|
|
|
|
|
|
y2(x) = |
|
z1(x) − |
|
z2(x) |
|||||||||||
2 i |
2 i |
=ea x cos (b x)
=ea x sin (b x)
zwei reelle Lösungen.
Die Berechnung der Wronski-Determinante
R |
a ea x cos b x |
|
b ea x sin |
|
b x |
|
R |
|
|
ea x cos |
( |
b x |
) |
R |
|
|
|
|
||
R |
|
|
|
|
|
|
R |
|
|
|
|
|
|
R |
|
|
|
|
|
|
R |
|
|
|
|
|
|
R |
( |
) − |
|
( |
|
) |
R |
|
|
||||
R |
|
|
|
|
|
|
R |
|
|
|
|
|
|
a ea x sin |
b x |
b ea x cos |
|
b x |
|
R |
|
|
ea x sin |
( |
b x |
) |
R |
|
|
|
|
R |
||
|
|
|
|
|
|
R |
|
|
|
|
|
|
R |
|
|
|
|
|
|
R |
|
|
|
|
|
|
R |
( |
) + |
|
( |
|
) |
R |
|
|
R |
||||
|
|
|
|
|
|
R |
|
|
|
|
|
|
R |
verläuft ähnlich zu den Berechnungen in Beispiel 12.22. Mit ein paar Umformungen erkennt man, dass die Wronski-Determinante den Wert e2 a x b hat. Der Imaginärteil b ist bei einer echt komplexen Zahl immer ungleich null. Somit sind y1 und y2 tatsächlich zwei reelle Fundamentallösungen.
Beispiel 12.28 (Einfache komplexe Eigenwerte)
a) Bei der Di erenzialgleichung |
|
|
|
|
|
|
|
|||||||||||||
y′′ |
|
9y |
|
|
0 |
|
|
|
|
λ2 |
+ |
9 0 |
|
|
|
|
||||
sind die+ Eigenwerte= Ô1 2 |
|
|
|
rein= |
imaginär. Die allgemeine Lösung der Di erenzialgleichung |
|||||||||||||||
|
|
|
|
|
|
|
λ |
|
, |
|
|
|
3 i |
|
|
|
|
|
|
|
besteht aus |
harmonischen Schwingungen: |
|
||||||||||||||||||
|
|
|
|
|
= ± |
|
|
|
|
|
|
|
|
|||||||
y(x) = C1 cos (3x) + C2 sin (3x) . |
|
|
||||||||||||||||||
b) Die charakteristische Gleichung der Di erenzialgleichung |
|
|||||||||||||||||||
y′′ |
|
4y′ 20y 0 |
|
|
Ô |
|
λ2 4λ 20 0 |
|
||||||||||||
hat die− Nullstellen+ = 1 |
, |
2 |
4 i |
. Der− |
Realteil+ =der Eigenwerte steht bei der Lösung |
|
||||||||||||||
|
|
|
e |
|
|
|
λ |
|
= |
2 |
± |
|
|
|
|
|
||||
y x |
|
2 x |
|
C1 cos 4 |
|
2 sin 4 |
|
|
|
|||||||||||
( ) = |
|
( |
|
|
|
x |
+ |
C |
|
|
|
|
x |
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
) |
Ì |
|||||
in der |
e |
-Funktion und die Imaginärteile erzeugen Sinusund Kosinusschwingungen. |
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
12.3 Lineare Di erenzialgleichungen |
485 |
(D) Mehrfache komplexe Eigenwerte
Jedes doppelte konjugiert komplexe Paar Eigenwerte λ = a ± i b einer linearen Di e- renzialgleichung mit konstanten Koe zienten erzeugt genau vier reelle Fundamentallösungen
y1 |
x |
ea x cos |
b x , |
y2 |
x |
ea x sin |
b x , |
y3 |
(x) = xea x cos |
(b x), |
y4 |
(x) = xea x sin |
(b x). |
||
|
( |
) = |
( ) |
|
( |
) = |
( ) |
Falls die Vielfachheit des komplexen Eigenwerts größer als zwei ist, entstehen jeweils durch Multiplikation mit x neue Fundamentallösungen
y5(x) = x2ea x cos (b x), y6(x) = x2ea x sin (b x), . . .
Beispiel 12.29 (Doppelte komplexe Eigenwerte)
Die Di erenzialgleichung
y(4) + 4 y′′′ + 8 y′′ + 8 y′ + 4 y = 0
hat die charakteristische Gleichung
λ4 + 4 λ3 + 8 λ2 + 8 λ + 4 = 0.
Durch Probieren gelangt man auf die Faktorisierung (λ2 + 2λ + 2)2 = 0 und damit auf die doppelten Eigenwerte λ1,2 = 1 + i und λ3,4 = 1 − i. Somit ist
y(x) = ex (C1 cos x + C2 sin x + C3x cos x + C4x sin x)
die allgemeine Lösung der Di erenzialgleichung. |
Ì |
Lösung einer homogenen linearen Di erenzialgleichung
Die allgemeine Lösung einer homogenen linearen Di erenzialgleichung der Ordnung n mit konstanten Koe zienten kann man durch folgende Schritte bestimmen:
(1)Berechne alle Eigenwerte aus der charakteristischen Gleichung anλn + an−1λn−1 + . . . + a1λ + a0 = 0.
(2)Bestimme zu jedem Eigenwert λi die Eigenfunktionen yi. Dabei sind Spezialfälle bei mehrfachen und komplexen Eigenwerten zu beachten.
(3)Die allgemeine Lösung besteht aus einer Linearkombination der Eigenfunktionen:
yh(x) = C1y1(x) + C2y2(x) + . . . + Cnyn(x).
486 |
12 Gewöhnliche Di erenzialgleichungen |
Beispiel 12.30 (Eigenwerte)
Die Eigenwerte der Di erenzialgleichung
y(5) − 3 y′′′ + 52 y′′ = 0
sind λ1,2 = 0, λ3 = −4 und λ4,5 = 2 ± 3 i, siehe Beispiel 12.24. Es kommen also sowohl eine doppelte reelle Lösung als auch ein komplex konjugiertes Paar als Eigenwerte vor. Die allgemeine Lösung dieser homogenen Di erenzialgleichung ist somit
y |
x |
) = |
C1 |
+ |
C2 x |
+ |
C3 e−4 x |
+ |
e−2 x |
‰ |
C1 cos |
3 x |
) + |
C2 sin |
( |
3 x |
)Ž |
. |
Ì |
( |
|
|
|
|
|
|
( |
|
|
|
|
Aufgrund der einfachen Struktur von linearen Di erenzialgleichungen erzeugen Störfunktionen in Form eines Polynoms in der Regel partikuläre Lösungen, die auch Polynome sind. Dasselbe Prinzip gilt auch für Exponentialfunktionen und gedämpfte und ungedämpfte harmonische Schwingungen, siehe Beispiel 12.23. Deshalb wählt man als Ansatz für eine partikuläre Lösung einfach die Störfunktion in einer etwas allgemeineren Form mit zunächst freien Koe zienten. Die folgende Übersicht zeigt die Ansätze für einige Typen von Störfunktionen.
Störansatztabelle
Störfunktion |
Ansatz für partikuläre Lösung |
|
|
Polynom vom Grad n |
Polynom vom Grad n |
r(x) = a0 + a1x + . . . + anxn |
yp(x) = A0 + A1x + . . . + Anxn |
Exponentialfunktion |
Exponentialfunktion |
r(x) = a ek x |
yp(x) = A ek x |
Harmonische Schwingung |
Harmonische Schwingung |
r(x) = a1 cos (ω x) + a2 sin (ω x) |
yp(x) = A1 cos (ω x) + A2 sin (ω x) |
Gedämpfte harmonische Schwingung |
Gedämpfte harmonische Schwingung |
r(x) = ek x (a1 cos (ω x) + a2 sin (ωx)) |
yp(x) = ek x (A1 cos (ω x) + A2 sin (ωx)) |
Falls die Störfunktion ein Polynom vom Grad n ist, dann verwendet man für die partikuläre Lösung ein Polynom vom selben Grad. Auch der Exponent k der e-Funktion und die Kreisfrequenz ω von Sinus und Kosinus sind bei der Störfunktion und dem entsprechenden Ansatz für eine partikuläre Lösung identisch.
Beim Störansatz muss stets der vollständige Term angesetzt werden. Besteht die Störfunktion nur aus an xn, so benötigt man für die partikuläre Lösung trotzdem einen Ansatz mit einem kompletten Polynom vom Grad n. Auch wenn die Störfunktion nur aus einem Sinus besteht, so braucht man für die partikuläre Lösung trotzdem einen Ansatz mit einer Sinusund einer Kosinusfunktion. Entsprechendes gilt auch für gedämpfte Schwingungen.
12.3 Lineare Di erenzialgleichungen |
487 |
Für inhomogene Di erenzialgleichungen gilt das Superpositionsprinzip. Besteht also die Störfunktion aus einer Summe von mehreren Termen, so kann zunächst für jeden Term separat eine partikuläre Lösung mittels Störansatz bestimmt werden. Anschließend werden die einzelnen partikulären Lösungen zu einer gesamten partikulären Lösung für die vollständige Störfunktion addiert.
Die Koe zienten im Störansatz versucht man, durch Einsetzen des Ansatzes in die Differenzialgleichung zu bestimmen. Dadurch erhält man ein lineares Gleichungssystem für die Koe zienten. Lassen sich die Koe zienten bestimmen, so war der Ansatz zielführend. Falls sich die Koe zienten nicht bestimmen lassen, so liegt in der Regel Resonanz vor. Bevor wir uns mit dem Erkennen und Behandeln von Resonanz beschäftigen, betrachten wir zunächst ein paar Beispiele ohne Resonanz.
Beispiel 12.31 (Störansatz)
a) Die Störfunktion der Di erenzialgleichung
y′′ − 3 y′ + 2 y = 4 − 8 x + 4 x2
besteht aus einem Polynom vom Grad 2. Laut Tabelle ist der Störansatz für eine partikuläre Lösung yp(x) = A0 + A1 x + A2 x2. Das Einsetzen der Funktion mit ihren Ableitungen in die Di erenzialgleichung kann mithilfe des Tableaus
yp |
x |
|
= |
A0 |
+ |
A1 x |
+ |
A2 x2 |
S ( |
2 |
|
yp |
(x) |
|
A1 |
2A2 x |
|
|
3) |
||||
′ |
(x) |
|
= |
2A2 |
+ |
|
|
|
|
S (−1) |
|
yp |
|
|
|
|
|
||||||
′′( ) |
2 |
= |
|
+ (2A1 − 6A2) x |
+ 2A2 x |
2 |
S ( |
) |
|||
4 − 8 x + 4 x |
|
= 2A0 − 3A1 + 2A2 |
|
|
|
erfolgen. In der letzten Gleichung steht sowohl auf der linken als auch auf der rechten Seite ein Polynom. Da die Gleichung für alle x erfüllt sein soll, müssen die Koe zienten der Polynome links und rechts übereinstimmen. Dadurch erhalten wir ein lineares Gleichungssystem mit drei Gleichungen für unsere gesuchten Größen A0, A1 und A2
2A2 |
+ |
|
|
|
= |
4 |
6A2 |
2A1 |
|
2A0 |
8 |
||
−2A2 |
3A1 |
+ |
= −4 |
|||
|
− |
|
|
= |
|
das die eindeutige Lösung A2 = 2, A1 = 2 und A0 = 3 besitzt. Eine partikuläre Lösung lautet
yp(x) = 3 + 2 x + 2 x2.
b) Bei der Di erenzialgleichung
y′′ − 3 y′ + 2 y = e3 x
lautet der Störansatz für eine partikuläre Lösung
yp |
( |
x |
Ae3 x, y |
p′ |
x |
|
′′ |
( |
x |
) = |
9Ae3 x. |
|
|
|
|
||||
|
3Ae3 x, yp |
|
|
|
|
||||||||||||||
|
) = |
|
|
|
|
( ) = |
|
|
|
|
|
|
|
|
|
|
|||
Diesen Ansatz setzt man in die Di erenzialgleichung ein |
|
|
|
||||||||||||||||
9A e3 x 3 3A e3 x 2A e3 x e3 x |
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
A |
|
1 |
|
partikuläre Lösung lautet yp |
|
x |
|
1 e3 x |
. |
|
||||||
und erhält− |
2 . Eine+ |
= |
|
|
|
|
|
( |
|
) = |
2 |
Ì |
|||||||
|
|
|
|
= |
|
|
|
|
|
|
|
|
|
|
|
|
488 |
12 Gewöhnliche Di erenzialgleichungen |
Bei inhomogenen Di erenzialgleichungen kann folgendes Phänomen auftreten: Ist die Störfunktion eine spezielle Lösung der homogenen Di erenzialgleichung, so verstärkt die Anregung durch diese Störfunktion die spezielle Lösung. Dieses Phänomen nennt man Resonanz.
Resonanz
Bei einer linearen Di erenzialgleichung mit konstanten Koe zienten spricht man von Resonanz, falls die Störfunktion in der allgemeinen Lösung der homogenen Di erenzialgleichung enthalten ist.
Beispiel 12.32 (Störansatz bei Resonanz)
Wir betrachten die Di erenzialgleichung aus Beispiel 12.31 mit einer anderen Störfunktion:
y′′ − 3 y′ + 2 y = e2 x .
Der Ansatz für eine partikuläre Lösung laut Störansatztabelle
|
( ) = |
A e2 x, y |
p′ ( ) = |
2 A e2 x |
′′ |
( |
|
) = |
4 A e2 x |
yp |
x |
x |
, yp |
|
x |
|
ergibt beim Einsetzen in die Di erenzialgleichung
4 A e2 x − 3 2 A e2 x + 2 A e2 x = e2 x.
´¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¸¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¹¶
0
Diese Gleichung ist für kein A erfüllbar. Somit führt dieser Ansatz nicht auf eine partikuläre Lösung. Es stellt sich die Frage, ob dieses Verhalten vorhersehbar ist. Bei der Störfunktion
r |
x |
) = |
e3 x |
in Beispiel 12.31 ist dieses Problem nicht aufgetreten. Die Beantwortung der Frage |
||||||||||||
|
( |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ist mithilfe der homogenen Lösung möglich. Die charakteristische Gleichung |
||||||||||||||||
|
|
|
λ2 |
|
3 λ |
2 |
|
0 |
|
|
|
|
|
|
|
|
liefert die Eigenwerte− + = |
λ1 |
= |
1 und λ2 |
= |
2. Die Störfunktion r x |
e2 x ist in der homogenen Lösung |
||||||||||
|
|
|
yh |
|
x C e |
x |
|
|
e |
2 x |
( ) = |
|
||||
|
|
|
|
( |
) = |
1 |
+ |
C |
2 |
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
enthalten. Es liegt Resonanz vor. In diesem Fall multipliziert man den Ansatz mit x:
|
|
p′ |
|
|
|
|
|
|
|
|
′′ |
|
|
|
+ |
|
|
) |
e2 x. |
|
|
|
yp x Ax e2 x, y |
x A 1 2 x e2 x, yp |
x A 4 |
4 x |
|
|
|||||||||||||||||
Eingesetzt( in) die= |
Di erenzialgleichung( ) = ( |
erhält+ ) man |
( ) = ( |
|
|
|
|
|
|
|||||||||||||
A 4 4 x e2 x |
3 A 1 2 x e2 x |
|
|
2Axe2 x e2 x. |
|
|
|
|
|
|
|
|
|
|
||||||||
Nach der (Division+ )durch− e2 x(vereinfacht+ ) |
sich+ |
die Gleichung= |
zu 4A |
|
|
3A |
|
|
1. Bemerkenswert hierbei |
|||||||||||||
ist, dass auf der linken Seite die Terme in |
x verschwinden. Die Gleichung ist also für A 1 lösbar |
|||||||||||||||||||||
|
x e |
2 x |
. |
|
|
− |
|
|
= |
|
= |
|
||||||||||
und eine partikuläre Lösung lautet yp |
( |
x |
) = |
|
|
|
|
|
|
|
|
|
|
Ì |
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
12.3 Lineare Di erenzialgleichungen |
489 |
Resonanztest
Bevor man eine partikuläre Lösung einer linearen Di erenzialgleichung mit konstanten Koe zienten berechnet, muss man unbedingt überprüfen, ob Resonanz vorliegt. Eine lineare Di erenzialgleichung mit konstanten Koe zienten kann man auf Resonanz überprüfen, indem man der Störfunktion einen Eigenwert zuordnet und diesen Eigenwert mit den Eigenwerten der homogenen Di erenzialgleichung vergleicht. Ist der zugeordnete Eigenwert der Störfunktion auch ein Eigenwert der homogenen Di erenzialgleichung, so liegt Resonanz vor.
Beispiel 12.33 (Resonanztest)
Die Di erenzialgleichung
y′′ − 3 y′ + 2 y = e2 x
besitzt die charakteristische Gleichung
λ2 − 3 λ + 2 = 0.
Daraus bestimmen wir die Eigenwerte λ1 |
2 und λ2 |
= |
1. Der Störfunktion können wir den Ei- |
|||||
genwert |
λ |
= |
2 |
zuordnen. Somit liegt |
Resonanz vor. |
Ì |
||
|
|
= |
|
Zuordnung von Eigenwerten bei Störfunktionen
Störfunktion |
Zugeordneter Eigenwert |
|
|
Polynom vom Grad n |
λ = 0 |
r(x) = a0 + a1x + . . . + anxn |
|
Exponentialfunktion |
λ = k |
r(x) = a ek x |
|
Harmonische Schwingung |
λ = ±ω i |
r(x) = a1 cos (ω x) + a2 sin (ω x) |
|
Gedämpfte harmonische Schwingung |
λ = k ± ω i |
r(x) = ek x (a1 cos (ω x) + a2 sin (ωx)) |
Auch im Resonanzfall lässt sich eine partikuläre Lösung finden. Man multipliziert den Ansatz laut Störansatztabelle mit dem Faktor xs. Dabei ist s die Vielfachheit des resonanten Eigenwerts.
Resonanzansatz
Wenn bei einer linearen Di erenzialgleichung mit konstanten Koe zienten einfache Resonanz vorliegt, multipliziert man den Störansatz mit x. Allgemein wird bei Resonanz mit einem Eigenwert der Vielfachheit s der Störansatz mit der entsprechenden Potenz xs multipliziert.
490 |
12 Gewöhnliche Di erenzialgleichungen |
Beispiel 12.34 (Ansätze bei Störfunktionen)
Die folgende Tabelle gibt eine Übersicht über die bisherigen Beispiele zu linearen Di erenzialgleichungen mit konstanten Koe zienten. Bei jeder homogenen Di erenzialgleichung werden zu unterschiedlichen Störfunktionen jeweils der entsprechende Ansatz für eine partikuläre Lösung angegeben.
Di erenzialgleichung |
Störfunktion |
Ansatz für partikuläre Lösung |
|
|
|||||||||||||||||||||||||||||||||||||||||||||||
Charakteristische Gleichung |
|
|
|
|
|
|
|
|
|
(Resonanzfälle sind rot markiert) |
|||||||||||||||||||||||||||||||||||||||||
Eigenwerte |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
y′ |
+ |
2 y |
= |
0 |
|
|
|
|
|
|
|
e−2 x |
|
|
|
|
x A e−2 x |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||
λ |
|
2 |
|
0 |
|
|
|
|
|
|
|
|
|
e2 x |
|
|
|
|
|
|
|
A e2 x |
|
|
|
|
|
A2 sin x |
|
e |
|
|
2 x |
|
|
|
|||||||||||||||
λ1+ |
|
|
=2 |
|
|
|
|
|
|
|
|
|
|
|
|
e |
|
2 x cos x |
|
|
( |
A1 cos x |
+ |
) |
− |
|
|
|
|||||||||||||||||||||||
|
= − |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
− |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
y′′ |
+ |
9 y |
|
|
= |
0 |
|
|
|
|
|
|
|
cos |
( |
3 x |
|
|
x |
( |
A1 cos |
( |
3 x |
) + |
A2 sin |
( |
3 x |
)) |
|
|
|||||||||||||||||||||
λ2 |
9 |
|
|
|
|
|
|
|
|
|
|
|
|
sin |
|
) |
|
x |
A |
1 |
cos |
3 x |
2 |
|
|
|
3 x |
|
|
||||||||||||||||||||||
+ |
= |
|
0 |
|
|
|
|
|
|
|
|
( |
3 x |
) |
|
|
( |
|
|
( |
) + |
A sin |
( |
)) |
|
|
|||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
λ1,2 |
= ± |
|
3 i |
|
|
|
|
|
|
|
|
e3 x cos |
|
3 x |
|
|
|
A1 cos |
|
3 x |
|
|
A2 sin |
|
|
3 x |
|
e3 x |
|||||||||||||||||||||||
y′′ |
− |
|
+ |
12 y |
= |
0 |
e3 |
x |
|
|
|
( |
) |
|
( |
|
|
|
|
x |
|
( |
|
|
|
) + |
|
|
|
|
( |
|
|
)) |
|
|
|||||||||||||||
7 y′ |
|
|
|
|
|
|
|
x A e3 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
λ2 |
− |
7 λ 12 |
4= |
0 |
|
e4 x |
|
|
|
|
|
|
x A e4 x |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||
λ1 |
|
3, λ+2 |
|
|
|
|
|
e |
|
3 x |
|
|
|
|
|
A e |
3 x |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||
y′′ |
= |
3 y′ |
|
|
= |
|
|
|
|
|
0 |
|
|
−x |
|
|
|
|
|
|
|
|
|
|
|
−x |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
− |
+ |
2 y |
|
= |
|
e− |
|
|
|
|
|
|
x A e− |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||
λ2 |
|
|
|
|
|
|
|
|
|
|
|
e |
− |
|
|
|
|
|
|
x A e |
− |
2 x |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
− |
3 λ 2 |
|
0 |
|
|
|
|
2 x |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||
λ1 |
|
|
1,+λ2 |
= |
|
|
|
2 |
|
ex |
|
|
|
|
|
|
|
A ex |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
y′′ |
= |
|
− |
|
|
|
|
|
|
= − |
|
|
e−3 |
x |
|
|
|
|
x A e−3 |
x |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
+ |
3 y′ |
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||
λ2 |
3 λ |
|
=0 |
|
|
|
|
|
|
|
e3 x |
|
|
|
|
|
|
|
A e3 x |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
λ1 |
+ |
|
0 |
, |
|
|
= |
|
|
|
|
3 |
|
|
x2 |
|
|
|
|
|
|
x |
‰ |
A0 |
|
A1 x |
|
|
A2 x2 |
Ž |
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
λ2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||
y′′ |
= |
2 y′ |
|
|
= − |
|
|
|
|
|
|
x |
|
|
|
|
|
|
|
|
|
|
|
+x |
|
|
|
+ |
|
|
|
|
|
|
|
|
|
|
|||||||||||||
+ |
|
|
y |
0 |
|
e− |
|
|
|
|
|
|
x2 A e− |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||
λ2 |
2 λ |
|
+1 |
|
=0 |
|
|
|
ex |
|
|
|
|
|
|
|
|
A ex |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
+ |
= − |
|
+ |
|
= |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
+ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
λ1,2 |
1 |
|
|
|
|
|
|
|
|
|
x |
|
|
|
|
|
|
|
|
|
A0 |
A1 x |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
y′′ |
− |
|
|
|
20 y |
= |
0 |
e2 |
x |
cos |
|
4 x |
) |
x |
( |
|
|
|
|
|
( |
4 x |
) + |
A2 sin |
( |
4 x |
)) |
e2 |
x |
||||||||||||||||||||||
4 y′ |
|
|
|
|
A1 cos |
|
|||||||||||||||||||||||||||||||||||||||||||||
λ2 |
4 |
|
|
|
+20 |
|
|
|
|
e |
|
|
|
|
|
( |
|
|
|
|
e |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||
− |
|
|
λ |
|
+ |
|
|
= |
0 |
|
|
2 x |
|
|
|
|
|
|
|
|
A |
|
|
2 x |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
1 |
|
|
2 |
|
|
|
|
|
|
sin |
|
4 x |
|
|
|
|
|
A1 cos |
|
4 x |
|
|
A2 sin |
|
|
4 x |
|
|
|
||||||||||||||||||||||
λ ,2 |
= |
|
± |
4 i |
|
|
|
|
|
( |
) |
|
|
|
|
( |
) + |
( |
) |
|
Ì |
||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|