Kurvenanpassung mit LSQNONLIN

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

Partner:

Forum

Option

[Erweitert]

• Diese Seite per Mail weiterempfehlen

Gehe zu:

MiWo

Forum-Newbie


	Beiträge: 4

	Anmeldedatum: 15.02.11

	Wohnort: ---

	Version: ---

Verfasst am: 15.02.2011, 11:51 Titel: Kurvenanpassung mit LSQNONLIN

Hi,
ich habe Messdaten (Impedanzwerte in Abhängigkeit von f) und will die entsprechenden Kurven (Z über f) mittels eines Ersatzschaltbildes modellieren.
Allerdings tritt das Problem auf, dass meine Startwerte kaum verändert werden, obwohl die modellierte Kurve noch stark von der echten Messkurve abweicht.
Hat Jemand ne Idee, woran es liegen könnte?

Hier mein Code:
clc;

global A;
global res;

%Einlesen der Messwerte
f=M16Olivenoel3V60Grad.freq;

w=2*pi*f;

Zmess(:,1)=M16Olivenoel3V60Grad.Zpr;
Zmess(:,2)=M16Olivenoel3V60Grad.Zpi;
Zmess(:,3)=M16Olivenoel3V60Grad.Z;

A=[w,Zmess];

%Startvektor [ Rd W Cd Re Cpr]
p = [1e11 1e11 1e-11 1e11 1e-10 ];
plb= [1e0 1e0 0 1e0 1e-13 ];
pub= [1e30 1e30 1 1e30 5e-9 ];

optioens = optimset('Diagnostics','on','Display','final','FunValCheck','on','MaxFunEvals',50000,'MaxIter',10000,'ScaleProblem','Jacobian','TolFun',1e-38,'TolX', 1e-3 Cool

;

L=lsqnonlin('wfun1', p, plb, pub, optioens, Zmess(:,1), Zmess(:,2), Zmess(:,3));

%Berechnung der modellierten Kurve
H1=(L(1)+L(2)./(sqrt(i*w)));

Zmod(:,1)=real((H1./(1i.*w.*L(3)))./(H1+1./(1i.*w.*L(3)))+(L(4)./(1i.*w.*L(5)))./(L(4)+1./(1*1i.*w.*L(5))));
Zmod(:,2)=imag((H1./(1i.*w.*L(3)))./(H1+1./(1i.*w.*L(3)))+(L(4)./(1i.*w.*L(5)))./(L(4)+1./(1*1i.*w.*L(5))));
Zmod(:,3)=abs((H1./(1i.*w.*L(3)))./(H1+1./(1i.*w.*L(3)))+(L(4)./(1i.*w.*L(5)))./(L(4)+1./(1*1i.*w.*L(5))));

%Plot der modellierten Kurve
loglog(M16Olivenoel3V60Grad.freq,M16Olivenoel3V60Grad.Zpr,M16Olivenoel3V60Grad.freq,Zmod(:,1))
figure (2)
loglog(M16Olivenoel3V60Grad.freq,M16Olivenoel3V60Grad.Zpi,M16Olivenoel3V60Grad.freq,Zmod(:,2))
figure (3)
loglog(M16Olivenoel3V60Grad.freq,M16Olivenoel3V60Grad.Z,M16Olivenoel3V60Grad.freq,Zmod(:,3))

und meine Funktion wfun1:

function [res] =wfun1(parameter,Re,Im,Fe)

global A;
global res;

w=A(:,1);

H1=(parameter(1)+parameter(2)./(sqrt(i*w)));

res(:,1)=Re-real((H1./(1i.*w.*parameter(3)))./(H1+1./(1i.*w.*parameter(3)))+(parameter(4)./(1i.*w.*parameter(5)))./(parameter(4)+1./(1*1i.*w.*parameter(5))));
res(:,2)=Im-imag((H1./(1i.*w.*parameter(3)))./(H1+1./(1i.*w.*parameter(3)))+(parameter(4)./(1i.*w.*parameter(5)))./(parameter(4)+1./(1*1i.*w.*parameter(5))));
res(:,3)=Fe-abs((H1./(1i.*w.*parameter(3)))./(H1+1./(1i.*w.*parameter(3)))+(parameter(4)./(1i.*w.*parameter(5)))./(parameter(4)+1./(1*1i.*w.*parameter(5))));
end

Vielen Dank schon mal vorab.

MiWo

Themenstarter

Forum-Newbie


	Beiträge: 4

	Anmeldedatum: 15.02.11

	Wohnort: ---

	Version: ---

Verfasst am: 16.02.2011, 15:41 Titel:

Hat keiner ne Idee? Liegt es vielleicht nur an den Einstellungen? Startwerte, Begrenzungen, Iterationsschritte, etc.
Falls ja wie könnte ich bessere Einstellungen finden?

Thomas84

Forum-Meister


	Beiträge: 546

	Anmeldedatum: 10.02.10

	Wohnort: ---

	Version: ---

Verfasst am: 16.02.2011, 17:11 Titel:

wie lautet denn das Abrruchkriterium?

Code:

L=lsqnonlin('wfun1', p, plb, pub, optioens, Zmess(:,1), Zmess(:,2), Zmess(:,3));

Funktion ohne Link?

laut http://www.gomatlab.de/posting.php?.....8cf646310f1e2d75186f8cc71 kann man lsqnonlin mit maximal 5 Argumenten aufrufen. Bist du sicher das Zmess überhaupt verwendet wird?

viele Grüße
Thomas

MiWo

Themenstarter

Forum-Newbie


	Beiträge: 4

	Anmeldedatum: 15.02.11

	Wohnort: ---

	Version: ---

Verfasst am: 16.02.2011, 20:51 Titel:

Danke für deinen Rat.
Ich habe es mal ausprobiert (indem ich Zmess anders übergeben habe), aber es hat trotzdem nicht funktioniert.

MiWo

Themenstarter

Forum-Newbie


	Beiträge: 4

	Anmeldedatum: 15.02.11

	Wohnort: ---

	Version: ---

Verfasst am: 28.02.2011, 19:43 Titel:

Hallo nochmal,
ich kann mein Problem jetzt besser eingrenzen.
Ich Kann mein ESB in zwei ESBs aufteilen und komme dadurch zu einer etwas unkomplizierteren Methodik.

Nun will ich 1/(R+W/sqrt(i*w))+i*w*C an meine Messkurve anpassen. (w=2*pi*f).

Es funktioniert auch mit LSQNONLIN, d.h. die Parameter R,W,C werden angepasst. Leider wird aber W immer auf den unteren Grenzwert gesetzt.
Weiß jemand warum bzw. was man dagegen tun kann.


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.