Nullstelle einer Symbolische Gleichung mit fminsearch finden

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

Partner:

Forum

Option

[Erweitert]

• Diese Seite per Mail weiterempfehlen

Gehe zu:

Hunni

Forum-Newbie


	Beiträge: 4

	Anmeldedatum: 13.08.17

	Wohnort: ---

	Version: ---

Verfasst am: 26.08.2017, 12:08 Titel: Nullstelle einer Symbolische Gleichung mit fminsearch finden

Hallo

ich habe eine symbolische Gleichung, die ich gerne mit fminsearch lösen würde.
Die Gleichung möchte ich mit diesem Code lösen. Im Moment verwende ich solve, jedoch ist das sehr langsam. In dem Struct M ist meine Gleichung gespeichert, mit den ganzen Parametern.
Die Gleichungen sind in der Form
$0 = c_{ein} + R - c_{aus}$

$R = f(c_{aus}) = k*c_{aus}$

$c_{ein} ist bekannt.$

Code:

function [ sim ] = contisimData( Para,IC,Temp,Tref,V,Para_isoth,M )

T = [Tref, Temp];

M.g = subs(M.g,[M.kiso,M.k, M.T],[Para_isoth, Para, T]);
expr = (M.x' + (M.vr./V).*M.S*M.g' - M.y' == zeros(max(size(M.y)),1));
expr = subs(expr, M.x, IC);
%Das function handle ist mein Problem
Eqsystem = matlabFunction(expr);

sim_out = fminsearch(Eqsystem,0);

sim = zeros(1,max(size(M.y)));
for i = 1:max(size(M.y))

name = string(['y', num2str(i)]);
sim(1,i) = double(sim_out(1).(name{1}));

end

end

Funktion ohne Link?

Mein Function Handle sieht dann wie folgt aus

Eqsystem =

function_handle with value:

@(y1,y2,y3,y4,y5)[-conj(y1)-(exp(-3.44803341351754).*(conj(y3.^(3.0./2.0)).*conj(y1)-(1.0e1.^9.180946579194002e-1.*conj(y5).*1.0e1)./conj(sqrt(y3))).*4.97.....

Mein Problem ist nun, dass man mit fminsearch kein Function Handle lösen kann, bei dem mehr input variablen vorhanden sind. Kann man das Function Handle irgendwie konvertieren, dass ich aus y1 y(1) machen kann?

Danke

Harald

Forum-Meister


	Beiträge: 24.492

	Anmeldedatum: 26.03.09

	Wohnort: Nähe München

	Version: ab 2017b

Verfasst am: 26.08.2017, 20:52 Titel:

Hallo,

Code:

eqsystem2 = @(y) eqsystem(y(1), y(2), y(3), y(4), y(5))

Funktion ohne Link?

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.