Euler Implizit zur Lösung einer DGL 1.Ordnung in Matrixform

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

Partner:

Forum

Option

[Erweitert]

• Diese Seite per Mail weiterempfehlen

Gehe zu:

Monika

Forum-Anfänger


	Beiträge: 17

	Anmeldedatum: 15.01.16

	Wohnort: ---

	Version: ---

Verfasst am: 01.04.2016, 13:34 Titel: Euler Implizit zur Lösung einer DGL 1.Ordnung in Matrixform

Hallo zusammen!

Ich habe ein Problem zur Berechnung einer DGL 1. Ordnung in Matrixschreibweise. Ich möchte das mit dem Euler-Rückwärts (implizit) Verfahren lösen, allerdings bekomme ich einige Fehlermeldungen.

Im Hauptskript rufe ich die Funktion "EulerImplizit" wie folgt auf:

Code:

dt = 1e-5; %Zeitschrittweite
tstart = 0; %Startzeitpunkt
tend = 1e-1; %Endzeitpunkt
tr = 1e-3; %Lade-Zeit für Dehnung
te =[tstart:dt:tend];
nt = length(te); %Anzahl Zeitschritte

A=rand(4,6);
D=rand(4,4);
L=rand(4,4);
e=4;

E=zeros(nt,6);
E(2:end,1)=1;
for i=1:(tr/dt)
E(1+i,1)=i*1e-2*0.01;
end
E(2+(tr/dt):end,1)=0.01;

xi=EulerImplizit(E,L,D,A,e);

Funktion ohne Link?

Der Euler-Rückwärts-Algorithmus lautet in der m-Datei "EulerImplizit"wie folgt:

Code:

function xi=EulerImplizit(L,D,A,E,e)
dt = 1e-5; %Zeitschrittweite
tstart = 0; %Startzeitpunkt
tend = 1e-1; %Endzeitpunkt
te =[tstart:dt:tend];
nt = length(te); %Anzahl Zeitschritte
xi = zeros(nt,e);
xip = zeros(nt,e);
for i=2:nt
dudu = xi(i-1,:)' + dt*(L_ab\D_ab*xi(i-1,:)' + L_ab\A_ai*E(i-1,:)');
xip(i,:) = dudu';

dudu = xi(i-1,:)' + dt*(L_ab\D_ab*xip(i,:)' + L_ab\A_ai*E(i,:)');
xi(i,:) = dudu';
end
end

Funktion ohne Link?

Folgende Fehlermeldungen bekomme ich dazu:

Error using \
Matrix dimensions must agree.

Error in EulerImplizit (line 12)
dudu = xi(i-1,: )' + dt*(L_ab\D_ab*xi(i-1,: )' + L_ab\A_ai*E(i-1,: )');

Error in Ordnungsreduktion_mit_Euler (line 154)
xi=EulerImplizit(E,L_ab,D_ab,A_ai,e);

Sieht jemand, was ich falsch gemacht habe? Rufe ich die Funktion richtig auf und ist die Parameterübergabe korrekt (damit kenne ich mich nicht so gut aus)? Und ich verstehe nicht, warum der Teil mit der Inversen (L_ab\D_ab bzw. L_ab\A_ai) nicht funktioniert.

Ich bin für jede Hilfe sehr dankbar!

Grüße,
Monika

Harald

Forum-Meister


	Beiträge: 24.501

	Anmeldedatum: 26.03.09

	Wohnort: Nähe München

	Version: ab 2017b

Verfasst am: 01.04.2016, 14:44 Titel:

Hallo,

setze doch einen Haltepunkt und überprüfe die Dimensionen beim Aufruf.

Mich verwundert auch, dass du keine Fehlermeldung bzgl. nicht definierter Variablen (z.B. L_ab) bekommst.

Grüße,
Harald


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 - 2025 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.