Parameterermittlung einer impliziten Funktion

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

Partner:

Forum

Option

[Erweitert]

• Diese Seite per Mail weiterempfehlen

Gehe zu:

thunder

Forum-Anfänger


	Beiträge: 11

	Anmeldedatum: 27.08.08

	Wohnort: ---

	Version: R2010a Unix (Ubuntu)

Verfasst am: 14.09.2010, 07:10 Titel: Parameterermittlung einer impliziten Funktion

Hallo zusammen,

ich möchte aus Messwerten (U,I) 7 Parameter einer impliziten nicht linearen Funktion ermitteln. Über die Matlab Hilfe bin ich auf lsqcurvefit gekommen.
Nur leider funktioniert bzw. versteh ich nicht so ganz wie das funktioniert.

zunächst habe ich meine Funktion definiert:

Code:

function F=myfun(x,U)

M=6;
N=105;

e=1.60217733e-19; % Elementarladung in As
k=1.380658e-23; % Boltzmann-Konstante in J/K
Tt=273.15+31.2;%T(1); % absolute Temperatur in K
Ut=((k*Tt)/e);

%x=[Iph Rs Ids Irs Rp m1 m2]
% 1 2 3 4 5 6 7

%F=I
F=(M*x(1)-M*x(3)*(exp((M*U+N*x(2)*I)/(M*N*x(6)*Ut))-1) ...
-M*x(4)*(exp((M*U+N*x(2)*I)/(M*N*x(7)*Ut))-1) ...
-(M*U+N*x(2)*I)/(M*N*x(5)));

Funktion ohne Link?

Der Aufmerksame Leser wird bemerken, dass es sich um das 2-Dioden-modell eines Solarmoduls handelt.

Wenn ich im nächsten Schritt meine die Funktion

Code:

function KL_Fit_einfach(U,I)
%x= [ Iph Rs Ids Irs Rp m1 m2]
x0 = [(12.2e-3); 13.2263/ 105; 2.5028e-08; 2.4575e-12; (1/0.0019)/6; 1.05; 2]; % Startwerte
[x,resnorm] = lsqcurvefit(@myfun_einfach,x0,U,I)

Funktion ohne Link?

fitten möchte habe ich das Problem, das in myfunc die Variable I nicht definiert ist. Wenn ich aber I in myfunc die Messwerte übergebe kommt nichts sinnvolles heraus.

muss ich als yWerte einen Vektor mit 0 übergeben und I in myfunc die Messwerte?

Das ist mir nicht ganz klar. Wäre super wenn mir hier jemand helfen könnte.

eupho

Forum-Meister


	Beiträge: 777

	Anmeldedatum: 07.01.09

	Wohnort: Marburg

	Version: R2009b

Verfasst am: 14.09.2010, 08:54 Titel:

Probier's mal mit:

Code:

fun = @(x)myfun(x, U, I);
[x, resnorm] = lsqcurvefit(fun, x0);

Funktion ohne Link?

thunder

Themenstarter

Forum-Anfänger


	Beiträge: 11

	Anmeldedatum: 27.08.08

	Wohnort: ---

	Version: R2010a Unix (Ubuntu)

Verfasst am: 14.09.2010, 09:59 Titel:

Hallo eupho,

dann bekomme ich immer die Fehler Meldung

Code:

??? Error using ==> lsqcurvefit at 175
LSQCURVEFIT requires four input arguments.

Error in ==> KL_Fit_einfach at 13
[x,resnorm] = lsqcurvefit(fun,x0)

Funktion ohne Link?

mein Code dazu

Code:

fun=@(x)myfun_einfach(x,U,I)

x0 = [(12.2e-3); 13.2263/ 105; 2.5028e-08; 2.4575e-12; (1/0.0019)/6; 1.05; 2]; % Startwerte
[x,resnorm] = lsqcurvefit(fun,x0)

Funktion ohne Link?

die Werte U und I liegen dabei im workspace.

eupho

Forum-Meister


	Beiträge: 777

	Anmeldedatum: 07.01.09

	Wohnort: Marburg

	Version: R2009b

Verfasst am: 14.09.2010, 10:59 Titel:

Ja, völlig korrekt. LSQCURVEFIT will die x/ydata sehen. Dann zurück zu deinem alten Code:

Code:

function F=myfun(x,U)

Funktion ohne Link?

ist das ein Schreibfehler? Weil dort fehlt ja das I als Input!

thunder

Themenstarter

Forum-Anfänger


	Beiträge: 11

	Anmeldedatum: 27.08.08

	Wohnort: ---

	Version: R2010a Unix (Ubuntu)

Verfasst am: 14.09.2010, 16:51 Titel:

das ist kein schreibfehler.

Wenn ich die Matlab Hilfe richtig verstehe werden nur Startwerte für x (also die Parameter die ermittelt werden sollen) und die x-Werte (Messwerte) übergeben.

Allerdings basiert das ja auf der Beispielfunktion ydata=x1*exp(x2*xdata).
Ist für mich explizit. Im Fall meiner Funktion heißt das ja das die y-Werte sowohl rechts als auch links des Gleichheitszeichens vorkommen. Etwa der Art:

ydata=x1*exp(x2*xdata+x3*ydata)
(nur noch komplizierter s.o.)

muss ich sowohl beim Aufruf von lsqcurvefit die y-Werte übergeben?

Code:

[x,resnorm] = lsqcurvefit(@myfun,x0,x-Werte,y-Werte)

Funktion ohne Link?

???
und in myfun auch noch mal? Oder bringe ich den y-Wert von links nach rechts das die Funktion Null wird und mit zeros einen Nullvektor übergeben, etwa so:

Code:

function F=myfun_einfach(x,U,I)
M=6;
N=105;
% U=x-Werte
% I=y-Werte
e=1.60217733e-19; % Elementarladung in As
k=1.380658e-23; % Boltzmann-Konstante in J/K
Tt=273.15+31.2;%T(1); % absolute Temperatur in K
Ut=((k*Tt)/e);
load I.mat %sonst ist I undefiniert

I=I./6;
%x=[Iph Rs Ids Irs Rp m1 m2]
% 1 2 3 4 5 6 7
F=(M*x(1)-M*x(3)*(exp((M*U+N*x(2)*I)/(M*N*x(6)*Ut))-1) ...
-M*x(4)*(exp((M*U+N*x(2)*I)/(M*N*x(7)*Ut))-1) ...
-(M*U+N*x(2)*I)/(M*N*x(5)))-I;

Funktion ohne Link?

In diesem Beipiel wäre F=0. muss ich in diesem Fall dann in lsqcurvefit für die y-Werte einen Vektor mit Nullen übergeben?


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