File:ODE Exp Imp Euler method.svg

Aus Wikimedia Commons, dem freien Medienarchiv
Zur Navigation springen Zur Suche springen

Originaldatei ‎(SVG-Datei, Basisgröße: 1.200 × 900 Pixel, Dateigröße: 10 KB)

Bildtexte

Kurzbeschreibungen

Ergänze eine einzeilige Erklärung, was diese Datei darstellt.

Beschreibung

[Bearbeiten]
Beschreibung
English: Compares the forward/explicit and backward/implicit Euler methods for the numerical solution of the ODE y'(x)=-200*(y-cos(x))-sin(x), y(0)=1. With the substituton z=y-cos(x), the ODE z'=-200*z, z(0)=0, results. Numerical integration of the second ODE should always result in the constant 0. By the nonlinear substitution in the first ODE, numerical errors are introduced in every step that are magnified by the factor 200. For not too small stepsizes, this results in explosive behavior of the explicit method whereas the implicit method remains close to the exact solution.
Datum
Quelle Eigenes Werk
Urheber LutzL


Lizenz

[Bearbeiten]
Ich, der Urheberrechtsinhaber dieses Werkes, veröffentliche es hiermit unter der folgenden Lizenz:
w:de:Creative Commons
Namensnennung Weitergabe unter gleichen Bedingungen
Diese Datei ist unter der Creative-Commons-Lizenz „Namensnennung – Weitergabe unter gleichen Bedingungen 3.0 nicht portiert“ lizenziert.
Dieses Werk darf von dir
  • verbreitet werden – vervielfältigt, verbreitet und öffentlich zugänglich gemacht werden
  • neu zusammengestellt werden – abgewandelt und bearbeitet werden
Zu den folgenden Bedingungen:
  • Namensnennung – Du musst angemessene Urheber- und Rechteangaben machen, einen Link zur Lizenz beifügen und angeben, ob Änderungen vorgenommen wurden. Diese Angaben dürfen in jeder angemessenen Art und Weise gemacht werden, allerdings nicht so, dass der Eindruck entsteht, der Lizenzgeber unterstütze gerade dich oder deine Nutzung besonders.
  • Weitergabe unter gleichen Bedingungen – Wenn du das Material wiedermischst, transformierst oder darauf aufbaust, musst du deine Beiträge unter der gleichen oder einer kompatiblen Lizenz wie das Original verbreiten.
GNU head Es ist erlaubt, die Datei unter den Bedingungen der GNU-Lizenz für freie Dokumentation, Version 1.2 oder einer späteren Version, veröffentlicht von der Free Software Foundation, zu kopieren, zu verbreiten und/oder zu modifizieren; es gibt keine unveränderlichen Abschnitte, keinen vorderen und keinen hinteren Umschlagtext.

Der vollständige Text der Lizenz ist im Kapitel GNU-Lizenz für freie Dokumentation verfügbar.

Du darfst es unter einer der obigen Lizenzen deiner Wahl verwenden.

Source

[Bearbeiten]

Using a small C program and gnuplot.

exp_imp_euler.c

[Bearbeiten]
/*
implizites und explizites Eulerverfahren für y'(x)=-2000*(y-cos(x))-sin(x)

*/

#include<stdio.h>
#include<math.h>

#define ALPHA 200

double f(double x,double y) {
	return -ALPHA*(y-cos(x))-sin(x);
	}
	
double explicit(double x,double y,double h) {
	return y+h*f(x,y);
}
	
double implicit(double x,double y, double h) {
	/* solve yp=y+h*f(xp,yp) for yp
		explicit solution is
		y=yp-h*f(xp,yp)
		 =(1+2000*h)*yp-2000*h*cos(xp)+h*sin(xp)
	*/
	double xp=x+h;
	
	return (y+h*(ALPHA*cos(xp)-sin(xp)))/(1+ALPHA*h);
}


int main(){
	double x,ye,yi,h;
	
	FILE *fp=fopen("imp_euler.dat","w");
	
	x=0; ye=yi=1; h=1.2e-2;
	
	for(;x<1; x+=h){
		fprintf(fp,"%8.5f\t%8.5f\t%8.5f\t%8.5f\n",x,cos(x),ye,yi);
		ye=explicit(x,ye,h);
		yi=implicit(x,yi,h);
	}
	return 0;
}

gnuplot commands

[Bearbeiten]
set term push
set term svg size 1200,900 fsize 28 linewidth 2

set output "exp_imp_euler.svg"

pl [:] [0:2] 'imp_euler.dat' u 1:3 ti 'explizit' w po pt 3, 'imp_euler.dat' u 1:4 w lin lt 3 lw 3 ti 'implizit'

unset outp
set term pop

Dateiversionen

Klicke auf einen Zeitpunkt, um diese Version zu laden.

Version vomVorschaubildMaßeBenutzerKommentar
aktuell16:30, 4. Jul. 2010Vorschaubild der Version vom 16:30, 4. Jul. 20101.200 × 900 (10 KB)LutzL (Diskussion | Beiträge){{Information |Description={{en|1=Compares the forward/explicit and backward/implicit Euler methods for the numerical solution of the ODE y'(x)=-200*(y-cos(x))-sin(x), y(0)=1. With the substituton z=y-cos(x), the ODE z'=-200*z, z(0)=0, results. Numerica

Du kannst diese Datei nicht überschreiben.

Keine Seiten verwenden diese Datei.

Globale Dateiverwendung

Die nachfolgenden anderen Wikis verwenden diese Datei:

Metadaten

Diese Datei enthält weitere Informationen, die in der Regel von der Digitalkamera oder dem verwendeten Scanner stammen. Durch nachträgliche Bearbeitung der Originaldatei können einige Details verändert worden sein.

Abgerufen von „https://commons.wikimedia.org/w/index.php?title=File:ODE_Exp_Imp_Euler_method.svg&oldid=431947054“