DGL mit Randwertproblem in Simulink lösbar?

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

Partner:

Forum

Option

[Erweitert]

• Diese Seite per Mail weiterempfehlen

Gehe zu:

Bernd1207

Forum-Newbie


	Beiträge: 1

	Anmeldedatum: 24.10.12

	Wohnort: ---

	Version: ---

Verfasst am: 24.10.2012, 17:30 Titel: DGL mit Randwertproblem in Simulink lösbar?

Ich nutze unten aufgeführte Matlab Funktion(en) zum Lösen einer DGL mit Randwertptoblem (Solver bvp4). Auf Grund von Kopplungsgründen mit einer anderen Software müsste diese Funktion in Simulink implementiert werden (S-Function?) mit dem Input der 3 Größen T_g, T_w und alpha. Leider übersteigt dies meine Fähigkeiten, falls dies überhaupt möglich ist... Kann jemand helfen?

Code:

function t_kopf=thermo(T_g, T_w, alpha)
% Solver initialisieren: Raster von 0-0.02m in 4 Schritten
solinit = bvpinit(linspace(0,0.022,4),@Thermo_init);
% Randwertproblem-Solver
sol = bvp4c(@Thermo_differentialEq,@Thermo_bc,solinit);

erg = deval(sol,0.022);
t_kopf = erg(1);

% Initiale Lösungsfunktion
function yinit = Thermo_init(x)
% yinit = [T_t(x) ; T_t'(x)]
yinit = [ T_w+x ; 1 ];
end

% Randbedingungen
function res = Thermo_bc(T_tA,T_tB)
% T_t(0) = T_w
% T_t'(L) = 0 K/m (= adiabater Kopf des Thermoelements)

res = [ T_tA(1)-T_w ; T_tB(2) ];
end

% Differentialgleichung
function dT = Thermo_differentialEq(x,T_t)
emiss = 0.74;
U = 0.00942;
A_q = 0.00000707;
sig = 5.67*10^(-8);
conduct_lambda = 18.2;

dT = [ T_t(2) ; ((alpha*U)/(conduct_lambda*A_q))*(T_t(1) - T_g) + ((emiss*U*sig)/(conduct_lambda*A_q))*(T_t(1)^4 - T_w^4) ];
end

end

Funktion ohne Link?

Es wäre mir eine sehr große Hilfe! Vielen Dank!

Harald

Forum-Meister


	Beiträge: 24.492

	Anmeldedatum: 26.03.09

	Wohnort: Nähe München

	Version: ab 2017b

Verfasst am: 24.10.2012, 18:22 Titel:

Hallo,

generell ist Simulink meines Wissens nicht zum Lösen von Randwertproblemen gedacht.
Ich sehe zunächst zwei Möglichkeiten - entweder das Problem umformulieren oder die Simulation mit verschiendenen Startwerten laufen lassen und iterativ herausfinden, für welche Startwerte die richtigen Randwerte herauskommen.

Die Frage ist allerdings, was genau diese "Kopplungsgründe" sind. Man kann aus MATLAB heraus so ziemlich alles aufrufen, was z.B. ein C- oder Active X-Interface hat.

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.