Lösbarkeit fmincon vs. lsqlin

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

Partner:

Forum

Option

[Erweitert]

• Diese Seite per Mail weiterempfehlen

Gehe zu:

wima2016

Forum-Newbie


	Beiträge: 5

	Anmeldedatum: 14.04.16

	Wohnort: ---

	Version: ---

Verfasst am: 14.04.2016, 17:06 Titel: Lösbarkeit fmincon vs. lsqlin

Hallo zusammen,

ich habe eine Zielfunktion zu minimieren:

Code:

fun=@(x) sqrt(sum((x-d).^2)).

Funktion ohne Link?

Statt dieser Funktion können wir auch die folgende Funktion betrachten:

Code:

fun=@(x) sum((x-d).^2) ,

Funktion ohne Link?

sodass wir ein linear least squares Problem haben. Zusätzlich gibt es auch Nebenbedingungen, deshalb verwende ich lsqlin:

Code:

Eye = eye(size(d,1))
[c,fval,residual,exitflag,output,lambda] = lsqlin(Eye,d,A,b,[],[],lb,ub).

Funktion ohne Link?

Für manche Input-Argumente bekomme ich die Nachricht

Code:

'Exiting: the constraints are overly stringent;
no feasible starting point found.'

Funktion ohne Link?

Deshalb habe ich es für die selben Argumente mal mit fmincon probiert:

Code:

f=zeros(length,1);
x0=linprog(f,A,b,[],[],lb,ub);
[c,fval]=fmincon(fun,x0,A,b,[],[],lb,ub);

Funktion ohne Link?

Jetzt das (für mich) überraschende Ergebnis: fmincon liefert für die selben Input-Argumente die Nachricht

Code:

'Local minimum found that satisfies the constraints.'

Funktion ohne Link?

Daher meine Frage: Wie kann es sein, dass fmincon ein Minimum findet, während lsqlin für die selben Input-Argumente kein Ergebnis findet?

Danke für eure Hilfe!

Harald

Forum-Meister


	Beiträge: 24.501

	Anmeldedatum: 26.03.09

	Wohnort: Nähe München

	Version: ab 2017b

Verfasst am: 14.04.2016, 17:37 Titel:

Hallo,

stelle bitte ein Beispiel zur Verfügung, an dem man das Verhalten reproduzieren kann.

Code:

f=zeros(length,1);

Funktion ohne Link?

erscheint mir wenig sinnvoll und erzeugt auch eine Fehlermeldung.

Grüße,
Harald

wima2016

Themenstarter

Forum-Newbie


	Beiträge: 5

	Anmeldedatum: 14.04.16

	Wohnort: ---

	Version: ---

Verfasst am: 14.04.2016, 18:08 Titel: Beispiel

Hallo,

Für bspw. die folgenden Parameter ist dies der Fall:

A= [0 0 0 0 -1 0 0 0;
0,956717930986567 0,0432820690134333 0 0 0 -1 0 0;
0 0 0,739292771154371 0,260707228845629 0 0 -1 0;
0 0 0 0 0 0 0 -1;
-1,38875538768198 0,171163690970216 0 0 0 0 0 0;
9,46622844285199 -3,72756482557360 0 0 0 0 0 0;
0,217591696711763 -1,38875538768198 0,171163690970216 0 0 0 0 0;
-4,73866361727840 9,46622844285200 -3,72756482557360 0 0 0 0 0;
0 0,217591696711763 -1,38875538768198 0,171163690970217 0 0 0 0;
0 -4,73866361727840 9,46622844285200 -3,72756482557361 0 0 0 0;
0 0 0,217591696711763 -1,38875538768198 0 0 0 0;
0 0 -4,73866361727840 9,46622844285200 0 0 0 0;
0 0 0 0 -1,30510677804058 0,112573951480094 0 0;
0 0 0 0 7,64454761066144 -2,45161049889983 0 0;
0,956717930986567 0,0432820690134333 0 0 0,192532826560482 -1,30510677804058 0,112573951480094 0;
-0,956717930986567 -0,0432820690134333 0 0 -4,19293711176161 7,64454761066144 -2,45161049889983 0;
0 0 0,739292771154371 0,260707228845629 0 0,192532826560482 -1,30510677804058 0,112573951480094;
0 0 -0,739292771154371 -0,260707228845629 0 -4,19293711176161 7,64454761066144 -2,45161049889983;
0 0 0 0 0 0 0,192532826560482 -1,30510677804058;
0 0 0 0 0 0 -4,19293711176161 7,64454761066144;
73,8903467618132 -32,5329111301404 0 0 0 0 0 0;
-25,5912942965706 45,7221193910067 -20,1308250944360 0 0 0 0 0;
0 -15,8354678855411 28,2920881173324 -12,4566202317913 0 0 0 0;
0 0 -9,79872452905236 17,5066742465212 0 0 0 0;
0 0 0 0 41,4284673344520 -15,2856855608204 0 0;
0 0 0 0 -8,93754096589453 14,1633215302800 -5,22578056438547 0;
0 0 0 0 0 -3,05551410743947 4,84207332969023 -1,78655922225076;
0 0 0 0 0 0 -1,04460124953701 1,65537964240751]

b= [-0,595451833081332;
0;
0;
0,
-0,416678044586527;
2,50566447351887;
-0,130641764601194;
0,130641764601194;
0;
0;
0;
0;
-0,742443073927751;
3,79659441151445;
0;
0;
0;
0;
0;
0;
19,1094331804346;
0;
0;
0;
19,9589856999065;
0;
0;
0]

lb =[0,317959817501480;
0,157210572517622;
0,0414816884607183;
0,00431907643367183;
0,597674187954456;
0,350362382106325;
0,114752826581231;
0,0153456105926976]

ub=[0,393461321453417;
0,291523921804051;
0,198972403792105;
0,123620445067927;
0,707871165935499;
0,588403789955047;
0,455276427507082;
0,323442021523119]

d=[0,317959817501480;
0,157210572517622;
0,0414816884607183;
0,00431907643367183;
0,597674187954456;
0,350362382106325;
0,114752826581231;
0,0153456105926976]

Und hier dann

Code:

f=zeros(8,1);

Funktion ohne Link?

Sorry, kürzeres Beispiel habe ich nicht gefunden.

Viele Grüße
Diana

Harald

Forum-Meister


	Beiträge: 24.501

	Anmeldedatum: 26.03.09

	Wohnort: Nähe München

	Version: ab 2017b

Verfasst am: 14.04.2016, 19:16 Titel:

Hallo,

wenn du den Startwert auch an lsqlin zur Verfügung stellst, findet lsqlin eine Lösung:

Code:

[c,fval,residual,exitflag,output,lambda] = lsqlin(Eye,d,A,b,[],[],lb,ub,x0)

Funktion ohne Link?

Wenn du statt des Standard-Algorithmus interior-point verwendest, findet lsqlin auch ohne den Startpunkt eine Lösung:

Code:

opts = optimoptions('lsqlin', 'Algorithm', 'interior-point');
[c,fval,residual,exitflag,output,lambda] = lsqlin(Eye,d,A,b,[],[],lb,ub,[], opts)

Funktion ohne Link?

Zu Empfehlungen zur Solverwahl siehe hier:
http://www.mathworks.com/help/optim.....ing-a-solver.html#bsbwx1l

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