Vektor als input einer ODE-Funktion, Schleife benötigt

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

Partner:

Forum

Option

[Erweitert]

• Diese Seite per Mail weiterempfehlen

Gehe zu:

Paul_ME

Forum-Newbie


	Beiträge: 6

	Anmeldedatum: 11.02.20

	Wohnort: ---

	Version: ---

Verfasst am: 17.02.2020, 12:36 Titel: Vektor als input einer ODE-Funktion, Schleife benötigt

Hallo zusammen,
ich benutze zurzeit eine ODE-Funktion, bei der ich als input einen Anfangsvektor gebe. Dieser Anfangsvektor habe ich IC genannt (siehe Code unten). Die Funktion berechnet neue Werte und gibt sie mir als Vektor zurück. Dieser neue Vektor habe ich Y-Vektor genannt (wird in einem anderen Skript zusammengestellt) und hat genau so viele Elemente wie der Anfangsvektor. Ich möchte nun diesen neuen Y-Vektor als input verwenden.

Die Fuktion soll immer den aktuellen Y-Vektor mit den neuesten Ergebnissen als input nehmen. Und dies so lange machen, wie meine Zeitvariable einen bestimmten Wert (t_process) erreicht hat (es werden also unteschiedliche Male sein, da ich t_process variieren werde). Meine Frage ist nun: Welche Art von Schleife soll ich am besten verwenden? und könnt ihr mir einen Tipp geben, wie ich das schaffe, dass der input immer der Y-Vektor mit den aktuellen Ergebnissen ist?

Lieben Dank im Voraus!

Code:

t = 0:1:t_process; % Zeitvariable
t_step = 0*3600 + 30*60 + 0;
N_time = 1;
m = 0;

while t <= t_process

% time integration with Matlab-ODE-Solver ode15s (ordinary differential equations system solver):
[t, Y] = ode15s(@(t,Y)PDESys_Paula(t, Y, L, W, epsilon, d_hydr, h_htx, a,b,f,...
d_part,tau, N_z, N_x, z, dz, x, dx, p_in_Value, T_in_htf,...
u_in_htf, N_z_plt, z_plt, dz_plt, T_plate, L_plt,...
p_th_set_Value, T_htf_out_set_Value, geoCase, htxCase,...
pinCase, procCase, contCase, pltCase, mtlCase, t_wait, R_vap, m_zeo, H_F, T_F, N_htx),...
[(N_time-1)*t_step N_time*t_step], IC);

N_t = length(t);

% Saving results:
% -------------------------------------------------------------------------
for k = 1:1:N_t

% save solution time vector
m = m + 1;
t_sol(m,1) = t(k);
t_sol = vertcat(t_sol);

% save solution for water uptake X
for i = 1:1:N_x
X_sol(m,1:N_z,i) = Y(k,((3*N_z)*(i-1)+1):((3*N_z)*(i-1)+N_z));
X_sol = vertcat(X_sol);
end

% save solution for temperature T
for i = 1:1:N_x
T_sol(m,1:N_z,i) = Y(k,((3*N_z)*(i-1)+(N_z+1)):((3*N_z)*(i-1)+2*N_z));
T_sol = vertcat(T_sol);
end

% save solution for pressure p
for i = 1:1:N_x
p_sol(m,1:N_z,i) = Y(k,((3*N_z)*(i-1)+(2*N_z+1)):((3*N_z)*(i-1)+3*N_z));
p_sol = vertcat(p_sol);
end

end

N_time = N_time + 1;
t = N_time*t_step;

end

Funktion ohne Link?


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.