Euler-Vorwärts-Verfahren in Matlab und Anfangswert

WICHTIG: Der Betrieb von goMatlab.de wird privat finanziert fortgesetzt. - Mehr Infos...

Partner:

Forum

Option

[Erweitert]

• Diese Seite per Mail weiterempfehlen

Gehe zu:

ianC

Gast


	Beiträge: ---

	Anmeldedatum: ---

	Wohnort: ---

	Version: ---

Verfasst am: 09.02.2009, 15:22 Titel: Euler-Vorwärts-Verfahren in Matlab und Anfangswert

Hallo Fachleute,
ich bin relativ oberflächlich in Matlab eingearbeitet und habe für euch eine sicherlich banale Frage.
Ich muss eine DGL mit Hilfe des eulerschen Streckungsverfahren in Matlab lösen, dabei verwende ich folgendes Skript:

Code:

% Beispiel zur numerischen Integration einer gewöhnlichen
% Differentialgleichung nach dem Euler'schen Streckenzugverfahren.

clear; clf % Löschen alter Variablen und alte Abbildung

% Festlegung von Integrationsintervall und Schrittweite

tl = 0.; % untere Grenze des Integrationsintervalls
tr = 5.; % obere Grenze des Integrationsintervalls
dt = 0.1; % Schrittweite des numerischen Schemas
fh = @(t,x) ((t-2).^2)+t./x; % Definition der Funktion über ein Handle

t=[tl:dt:tr];x(1)=2; % Randbedingung vorgeben

for i=2:length(t);
x(i)=x(i-1) + dt*feval(fh,t(i-1),x(i-1));
end

plot(t,x,':ob');
% %%%%%%%%%%%%%%Ende numerische Lösung%%%%%%%%%%%%%%%%%%%%%%%%%%%

% Darstellung von numerischer und analytischer Lösung
xlabel('Zeit t','Fontsize',13);
ylabel('x','Fontsize',13);
title('Euler-Verfahren Vorwaerts','Fontsize',14);

Funktion ohne Link?

Wenn ich mir die Lösung darstellen lasse fängt die y-achse mit x=2 halt erst bei 2 an, dabei möchte ich die y-achse von 0 bis 15 laufen lassen.
Wenn ich für x(0)=1 einsetze erhalte ich eine Fehlermeldung, dass x positiv oder ein logical sein muss.
Wisst ihr wie ich die Darstellung von 0 bis 15 auf der y-achse erzielen kann?

Verzeihung für die banale Frage und danke im voraus für hilfreiche Antworten.

Edit by Bijick: Code-Formatierung ergänzt. Bitte beim nächsten Mal selbst machen (Code-Button mittig über dem Eingabefenster). Danke

Epfi

Forum-Meister


	Beiträge: 1.134

	Anmeldedatum: 08.01.09

	Wohnort: ---

	Version: ---

Verfasst am: 09.02.2009, 15:32 Titel:

Die Indizierung von Vektoren fängt bei Matlab immer bei 1 an. Das "nullte" Element gibt es im Gegensatz zu C und fast allen anderen Sprachen nicht. Das erste Element hat den Index 1. x(0) ist also grundsätzlich etwas, womit Matlab nichts anfangen kann.

Außerdem empfiehlt es sich, weder i noch j als Zähler zu benutzen, weil die in Matlab für die imaginäre Einheit reserviert sind ( i * i = -1).

Zum eigentlich Problem:

Code:

axis = [x_min x_max y_min y_max];

Funktion ohne Link?

könnte Dir weiterhelfen.


Einstellungen und Berechtigungen
Beiträge der letzten Zeit anzeigen:	Du kannst Beiträge in dieses Forum schreiben. Du kannst auf Beiträge in diesem Forum antworten. Du kannst deine Beiträge in diesem Forum nicht bearbeiten. Du kannst deine Beiträge in diesem Forum nicht löschen. Du kannst an Umfragen in diesem Forum nicht mitmachen. Du kannst Dateien in diesem Forum posten Du kannst Dateien in diesem Forum herunterladen

Impressum | Nutzungsbedingungen | Datenschutz | FAQ |

RSS

Hosted by:

Copyright © 2007 - 2024 goMatlab.de | Dies ist keine offizielle Website der Firma The Mathworks

MATLAB, Simulink, Stateflow, Handle Graphics, Real-Time Workshop, SimBiology, SimHydraulics, SimEvents, and xPC TargetBox are registered trademarks and The MathWorks, the L-shaped membrane logo, and Embedded MATLAB are trademarks of The MathWorks, Inc.