• Keine Ergebnisse gefunden

Numerik II (Numerik gewöhnlicher Differentialgleichungen)

N/A
N/A
Protected

Academic year: 2021

Aktie "Numerik II (Numerik gewöhnlicher Differentialgleichungen)"

Copied!
2
0
0

Wird geladen.... (Jetzt Volltext ansehen)

Volltext

(1)

BTU COTTBUS

LEHRSTUHLNUMERISCHE MATHEMATIK UND

WISSENSCHAFTLICHESRECHNEN Prof. Dr. G. Bader, Dipl.-Math. Friedemann Kemm

Numerik II (Numerik gewöhnlicher Differentialgleichungen)

Serie 3

www.math.tu-cottbus.de/˜kemm/lehre/ode/index.html

Aufgaben zur Übung am 30.04.04

Liegengebliebene Aufgaben der beiden letzten Übungsblätter Hausaufgaben zur Abgabe am 30.04.04

1. Man bestimme alle Lösungen der AWA y0=p

1−y2, y(0) =−1

2. Zur Konstruktion von Differenzenformeln für die skalare Differentialgleichung u0(t) = f(t,u(t)) über die Taylor-Entwicklung benötigt man die k-ten Ableitungen der Lösung u(t). Dazu wird der von der rechten Seite f abhängige Differentialoperator D definiert:

D(ϕ(t,u)):=ϕt(t,u) +f(t,u)ϕu(t,u)

Die k-te Ableitung von u läßt sich durch (k-1)-malige Anwendung von D auf f berechnen:

u0 = D0(f) = f(t,u) u(k) = Dk−1(f) Dk(ϕ) := D(Dk−1(ϕ)) Desweiteren sei der folgende Operator definiert:

Dn(ϕ):=

n k=0

n k

nϕ

∂tn−k∂uk·fk

Beweisen Sie die Relation D(Dn(ϕ)) =Dn+1(ϕ) +n Dn−1u)D(f) und stellen Sie die ersten vier Ableitungen von u mit Hilfe der Operatoren Dkdar.

Konstruieren Sie die Ableitungen für das AWP y0=y2,y(0) =1 ! 3. Gegeben sei die Differentialgleichung zweiter Ordnung

u00(t) + (u0(t))2=0. (1)

Ihre allgemeine Lösung lautet

u(t) =ln(c1t+c2) mit reellen Konstanten c1und c2.

(2)

a) Stellen Sie für diese Differentialgleichung entsprechend Satz 1.3 im Skriptum die Matrix- Differentialgleichung für die Abhängigkeit von den Anfangsdaten auf.

b) Lösen Sie das Anfangswertproblem für (1) mit den Anfangsdaten u(0) =1; u0(0) =0

und lösen Sie die Matrix-Differentialgleichung aus dem ersten Teil für diese Nominallö- sung.

Referenzen

ÄHNLICHE DOKUMENTE

Test the stability of the solution to (2) with respect to changes in the initial conditi- ons by changing x 0 = 80, y 0 = 30 by a unit amount in each direction (four different

Hint: The midpoint rule is only stable if |hλ n−1 (A)| < 2, where λ n−1 (A) is the lar- gest eigenvalue of A.. Implement the trapezoidal method on the linear system of

for solving ordinary differential equations and systems of ordinary differential equations using step size control, where fun is a function handle, y0 an initial condition, t =