ODE Problem

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

Partner:

Forum

Option

[Erweitert]

• Diese Seite per Mail weiterempfehlen

Gehe zu:

mike793

Forum-Century


	Beiträge: 119

	Anmeldedatum: 11.08.08

	Wohnort: ---

	Version: ---

Verfasst am: 24.09.2009, 12:23 Titel: ODE Problem

Hallo zusammen

Ich habe ein großes Problem mit dem ODE Befehl

Ich will ein Datenfile Einlesen und mit diesen Daten einen Verlauf den der ODE Befehl verwerten kann erzeugen. das geht ja sehr gut schon, aber mein Problem ist das mit dem ODE 45 Befehl!!!! ich bekomme hier immer die gleiche Fehlermedlung die lautet:

Code:

??? Attempted to access u(0); index must be a positive integer or
logical.

Error in ==> DIFF_Step at 42
u(llll)=u(end);

Error in ==> odearguments at 110
f0 = feval(ode,t0,y0,args{:}); % ODE15I sets args{1} to yp0.

Error in ==> ode45 at 173
[neq, tspan, ntspan, next, t0, tfinal, tdir, y0, f0, odeArgs,
odeFcn, ...

Error in ==> xpunkt_stoch at 19
[t,z] = ode45 (@DIFF_Step ,0:0.01:20,z0,[],m1,m2,c1,c2,d1,d2);

Funktion ohne Link?

Ich stelle auc mal den code rein um einen besseren überblick zu bekommen!!!!!

Code:

function dzdt = DIFF_Step (t,z,m1,m2,c1,c2,d1,d2)

fid=fopen('C:\1Ordner_DA\Sch.ASC','r');
list=dlmread('C:\1Ordner_DA\Sch.ASC');
fclose(fid);

%lll=0:0.01:22;
%t=0:0.01:20;
tt=t(1:end-1);

tt=tt*100;
t_n=floor(tt);
t1=t_n+1;
t2=t_n+2;
t_1=tt(1:end-1);
t_2=tt(2:end);
delta_t=t_2-t_1;
delta_L=list(t2)-list(t1);
delta_L=delta_L(1:end-1);
k=delta_L'./delta_t;
hh=tt-t_n;
hh=hh(1:end-1);
list=list(1:end-2);
rr=k.*(hh);
t1=t1(1:end-1);
u=list(t1)'+rr;
tt=tt(1:end-1);
tt=tt/100;

for llll=1:length(t)
if llll<= length(u)
u(llll)=u(llll);
else
u(llll)=u(end);
end

end

plot(llll,u)

dzdt = [z(2)
(-d1-d2)/m1*z(2)+(-c1-c2)/m1*z(1)+c2/m1*z(3)+d2/m1*z(4)+c1/m1*u+d1/m1*u'
z(4)
-d2/m2*z(4)-c2/m2*z(3)+d2/m2*z(2)+c2/m2*z(1)];

%%%%und

clear;
% ------------------------------------------
t0 = 0;
tf = 20;
z0 = [0 0 0 0];
tol = 1.e-3;
trace = 0;
%Parameter
%--------------------------------------------
m1=50;
m2=500;
d1=40;
d2=2800;
c1=210000;
c2=25000;

% --------------------------------------------
[t,z] = ode45 (@DIFF_Step ,0:0.01:20,z0,[],m1,m2,c1,c2,d1,d2);
%

Funktion ohne Link?

DANKE!!![/code]

Harald

Forum-Meister


	Beiträge: 24.448

	Anmeldedatum: 26.03.09

	Wohnort: Nähe München

	Version: ab 2017b

Verfasst am: 24.09.2009, 12:43 Titel:

Ein Problem, das ich sehe: in DIFF_STEP wird t als Skalar übergeben, d.h. tt ist wahrscheinlich leer, und von da an wirds zumindest nicht so laufen wie geplant.

Grüße,
Harald

mike793

Themenstarter

Forum-Century


	Beiträge: 119

	Anmeldedatum: 11.08.08

	Wohnort: ---

	Version: ---

Verfasst am: 24.09.2009, 12:49 Titel:

Danke für deine rasche antwort!!!!

Bist du dir da sicher das im ode Befehl die Zeit also t als skalar übergeben wird?

Denn interessant ist wenn ich im Abschnitt:

Code:

for llll=1:length(t)
if llll<= length(u)
u(llll)=u(llll);
else
u(llll)=u(end); %%%-70
end

end

Funktion ohne Link?

statt u(end) jetzt -70 hineinschreibe integriert er mir also nimm er nur den letzten wert!!!!

Soll ich nun das Mfile so auslegen das er mit jeder neuen Zeit den richtigen Wert eraussucht oder weiter mit vektoren?

Danke

mike793

Themenstarter

Forum-Century


	Beiträge: 119

	Anmeldedatum: 11.08.08

	Wohnort: ---

	Version: ---

Verfasst am: 24.09.2009, 13:12 Titel:

hallo zusammen nochmals

ich habe nun das Problem mit sklaren gelöst nun kommt aber als ergebnis heraus NaN im Vektor das ist auch nicht die lösung!!!!!

was meint ihr wie kann an das lösen!!!!

Code:

t00=t*100;
tn=floor(t00);
t1=tn+1;
t2=tn+2;
delta_list=list(t2)-list(t1);
delta_t=t-tn;
k=delta_list/delta_t;
tt=t-tn;
u=list(t1)+k*tt;%% der code für den sklara!!!

Funktion ohne Link?

Harald

Forum-Meister


	Beiträge: 24.448

	Anmeldedatum: 26.03.09

	Wohnort: Nähe München

	Version: ab 2017b

Verfasst am: 24.09.2009, 14:47 Titel:

Ich kann der Logik nicht ganz folgen, sehe aber ein Problem:
t = 0
-> t00 = 0
-> tn = 0
-> delta_t = 0
-> (möglicherweise) delta_list = 0
-> k = 0 / 0 = NaN

aber was soll das t00 = t*100 eigentlich? Und was soll t-tn sein (das eine ist mit 100 skaliert, das andere nicht...)?

Hoffe, das hilft.

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