Hilfe bei Eingabe für Funktion (Levenberg Marquardt Alg.)

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

Partner:

Forum

Option

[Erweitert]

• Diese Seite per Mail weiterempfehlen

Gehe zu:

helpsearcher

Forum-Newbie


	Beiträge: 1

	Anmeldedatum: 08.05.13

	Wohnort: ---

	Version: ---

Verfasst am: 08.05.2013, 17:04 Titel: Hilfe bei Eingabe für Funktion (Levenberg Marquardt Alg.)

Hallo Matlab-Programmierer!
Ich arbeite zurzeit an einem Programm, das mithilfe des Levenberg-Marquardt-Algorithmus Werte bestimmt.
Aufgerufen wird es folgendermaßen:

Code:

[x,ssq,cnt] = LMFnlsq(res,x0);

Funktion ohne Link?

Es war auch ein kleines Beispiel dabei:
--------------------------------------------------------

Code:

a = 10;
b = -8;
c = -.2 % Model parameters
t = (0:.1:1)'; % should be a column vector
f = a + b*(1-exp(t/c)) + .2*randn(size(t));
% Column vector of residuals
res = @(x) x(1) + x(2)*(1-exp(t/x(3))) - f;
% solve it with rather bad estimate of parameters
x0 = [5,3,-.1];
[x,ssq,cnt] = LMFnlsq(res,x0);
figure;
plot(t,f, t,res(x)+f,'r'), grid

Funktion ohne Link?

------------------------------------------------------
Bei meinem Programm ist f (siehe Beispiel) eine gefittete Funktion, also vom Typ cfit:
------------------------------------------------------

Code:

myfit(t) = r1*(1-exp(-t/tau1))+r2*(1-exp(-t/tau2))+r3*(1-exp(-t/tau3))
Coefficients (with 95% confidence bounds):
r1 = 82.06 (72.67, 91.45)
r2 = 12.85 (10.52, 15.17)
r3 = 82.41 (72.3, 92.52)
tau1 = 14.38 (12.25, 16.51)
tau2 = 0.06762 (0.0008762, 0.1344)
tau3 = 3.31 (2.883, 3.737)

Funktion ohne Link?

------------------------------------------------------
Allerdings ist f (siehe Beispiel) vom Typ double und meine Funktion(myfit) vom Typ cfit. Umwandeln könnte ich diese
Funktion in einen string mit :
------------------------------------------------------

Code:

funktion_fit=formula(myfit);

Funktion ohne Link?

------------------------------------------------------
Auf die Koeffizienten gelange ich mit:
------------------------------------------------------

Code:

coeff_val=coeffvalues(myfit);
coeff_fit=coeffnames(myfit);

Funktion ohne Link?

------------------------------------------------------
Das nächste Problem ist die "function handle" bei:
------------------------------------------------------

Code:

res=@(x) x(1) + x(2)*(1-exp(t/x(3)))-f;

Funktion ohne Link?

------------------------------------------------------
Die Funktion: x(1) + x(2)*(1-exp(t/x(3))) liegt bei mir ebenfalls als string vor (wird dynamisch erstellt):
------------------------------------------------------

Code:

funktion_string =

'x(1)*exp(-t/(x(2)))+x(3)*exp(-t/(x(4)))+x(5)*exp(-t/(x(6)))'

Funktion ohne Link?

------------------------------------------------------
Wie kann ich das in die Form wie beim Beispiel bringen.
Ich würde mich sehr freuen wenn ihr mir helfen könntet.

Vielen Dank

[EDITED, Jan, Bitte Code-Umgebung verwenden - Danke!]

LMFnlsq.m

Beschreibung:

Hier die Funktion

Download

Dateiname:

LMFnlsq.m

Dateigröße:

17.8 KB

Heruntergeladen:

435 mal


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.