0 Daumen
902 Aufrufe

Aufgabe:


Konvergenz gegen die Exponentialfunktion.
Gegeben sei das Anfangswertproblem
y=y+1,y(0)=0 y^{\prime}=y+1, \quad y(0)=0
mit exakter Lösung
y(x)=ex1 y(x)=\mathrm{e}^{x}-1
Sei xend >0 x_{\text {end }}>0 beliebig. Wir verwenden eine konstante Schrittweite h=xend N h=\frac{x_{\text {end }}}{N} mit einem NN N \in \mathbb{N} . Die (endliche) Folge (yk)k=0,1,,N \left(y_{k}\right)_{k=0,1, \ldots, N} werde mit dem expliziten Euler-Verfahren berechnet. Zeigen Sie durch direktes Ausrechnen
limNyN=y(xend ) \lim \limits_{N \rightarrow \infty} y_{N}=y\left(x_{\text {end }}\right)


Hallo zusammen, könnte mir jemand bitte dabei helfen?

Danke im Voraus! :)

Avatar von

1 Antwort

0 Daumen
 
Beste Antwort

Das Eulerverfahren ist definiert durch

yk+1=yk+h(yk+1) y_{k+1} = y_k + h (y_k + 1) mit y0=0 y_0 = 0 .

Also y1=h=(1+h)1 y_1 = h = (1+h)-1

y2=(1+h)21 y_2 = (1+h)^2 - 1 und allgemein

yk=(1+h)k1=(1+xendN)k1 y_k = (1+h)^k - 1 = \left(1+\frac{x_{end}}{N} \right)^k - 1

Damit gilt dann yN=(1+xendN)N1 y_N = \left( 1+\frac{x_{end}}{N} \right)^N - 1

Weiter gilt limN(1+xendN)N=exend \lim_{N \to \infty} \left( 1+\frac{x_{end}}{N} \right)^N = e^{x_{end}}

In Summe also limNyN=exend1 \lim_{N \to \infty} y_N = e^{x_{end}}-1

Avatar von 39 k


Danke erstmal für deine Antwort!

Könntest du mir bitte erklären Wie du darauf kommst?

Berechne y2 y_2 und auch y3 y_3   dann siehst Du es.

y2 = h + h(h+1) = h+h+h = 2h + h2

Binomischer Lehrsatz????

https://www.mathelounge.de/877828/lokaler-diskretisierungsfehler-im-…


Könntest du mir bitte nur bei b) helfen?

Danke im Voraus!

Ein anderes Problem?

Stell deine Frage