Schreiben Sie Eulers Antwort bei x = 8 in geschlossener Form auf, ohne einen einzigen Schritt auszuführen, und erklären Sie dann, warum die Halbierung von h den Fehler nicht halbiert. Dies ist y′ = x − y mit Startwert (0, 1), h = 0,2 und 40 Schritten.
-
Beginnen Sie mit der Antwort, die Euler anstrebt. Die Gleichung ist linear, x − 1 ist eine partikuläre Lösung, und die Anfangsbedingung legt den Koeffizienten des abklingenden Teils auf 2 fest.
y′=x−y, y(0)=1 ⟹ y=x−1+2e−x,y(8)=7+2e−8=7.000671
-
Eulers Regel für diese Gleichung ist kein allgemeines nichtlineares Durcheinander — sie ist eine affine Rekursion. Und dasselbe x − 1, das die Differentialgleichung gelöst hat, löst die Rekursion exakt, Schritt für Schritt. Euler trifft den linearen Teil der Antwort vollkommen exakt.
yn+1=yn+h(xn−yn)=(1−h)yn+hxn,yn=xn−1 ⟹ yn+1=xn−1+h=xn+1−1
-
Der gesamte Fehler steckt also in dem, was übrig bleibt, und was übrig bleibt, ist eine geometrische Folge mit dem Quotienten 1 − h. Vierzig Faktoren von 0,8 landen auf dem Endpunkt des Bedienfelds, ohne dass überhaupt ein Schritt ausgeführt wurde.
yn=xn−1+2(1−h)n⟹y40=7+2(0.8)40=7+2.658456×10−4=7.000266
-
Stellen Sie nun die beiden abklingenden Teile nebeneinander: Der wahre ist e−x, Eulers ist (1 − h)x/h. Beide sind Exponentialfunktionen in x, und δ ist die gesamte Differenz zwischen ihren Raten. Seine Reihe beginnt bei xh/2 — dieses führende h ist der gesamte Inhalt des Begriffs “erster Ordnung”.
(1−h)x/h=e−x−δ,δ=x(h−ln(1−h)−1)=2xh+3xh2+⋯=0.925742
-
Der Fehler ist das, was übrig bleibt, wenn die beiden Exponentialfunktionen subtrahiert werden. Erste Ordnung sagt voraus, dass die Halbierung von h ihn halbiert, sodass das Verhältnis 2 betragen sollte.
Eh=2e−x(1−e−δ): E0.2=0.000405, E0.1=0.000234,E0.1E0.2=1−e−0.4288411−e−0.925742=0.3487370.603763=1.731285
Antwort
Euler endet bei 7,000266 gegenüber einem exakten Wert von 7,000671, was einem globalen Fehler von 0,000405 entspricht; halbiert man die Schrittweite, fällt er nur auf 0,000234, ein Verhältnis von 1,731285, wo erste Ordnung 2 versprochen hatte. Das Versprechen ist nicht gebrochen — es hat noch gar nicht begonnen. Es setzt voraus, dass δ ≈ xh/2 klein ist, und bei x = 8 mit h = 0,2 beträgt diese Zahl 0,8, sodass der lineare Term des Fehlers noch immer mit jedem dahinterstehenden Term konkurriert; um das beobachtete Verhältnis auf 1,9933 anzuheben, benötigt man h = 0,005, was 1600 Schritten über dieselben acht Einheiten entspricht. Die geschlossene Form gibt die präzisere Interpretation preis: Euler approximiert diese Lösung nicht, sondern löst eine andere Gleichung exakt, wobei ein geometrisches Abklingen von 0,8n an die Stelle von e−x tritt — und 0,8 ist der Multiplikator 1 + hλ, den das Bedienfeld bereits anzeigt. Allein an diesem abklingenden Teil gemessen liegt Euler um 60% falsch. Der globale Fehler sieht nur deshalb winzig aus, weil der Teil, den er falsch wiedergibt, nahezu abgeklungen ist, bis x 8 erreicht.
Eh/2Eh=1−e−δh/21−e−δh⟶2(δ≪1),δh≈2xh=0.8