Verfasst am: 11.08.2010, 10:55
Titel: Problem mit lsqcurvefit:first order optimality less than opt
Hallo zusammen!
Ich habe eine Verteilung und würde da gerne eine Kurve drüber legen.
Genauer: Ich habe Produktionsdaten (Produzierte Mengen) und versuche diese mit einer Todzeit und einer Verzögerung dritter Ordnung anzunähern und suche die optimalen Parameter für die Todzeit und die Verzögerungszeit.
Jetzt sagt aber Matlab:
Code:
Optimization terminated: first-order optimality less than options. TolFun, and no negative/zero curvature detected in trust region model
Ich vermute, dass der Gradient zu klein ist damit Matlab optimieren kann.
In den Options habe ich "eps" als toleranz angegeben.
Wenn ich nun verschiedene Startwerte für die beiden Parameter angebe ändert sich die "RERNORM" entscheidend, also der "Fehler" ändert sich in der Größenordnung 10 (ein Startwert 0.00034, ein anderer Startwert 0.00034). Man kann also deutlioch den Unterschied erkenne und es ist bei weiterm nicht in der Nähe von "eps".
Allerdings ist der Fehler der gleiche.
Kann mir jemand sagen was falsch sein könnte?
Falls ihr mein Modell haben möchtet, kann ich das gerne hier posten.
da ist vermutlich gar nichts falsch. Optimization terminated ist neutral zu sehen, d.h. weder positiv noch negativ. Aller Wahrscheinlichkeit nach befindet man sich in einem lokalen (!) Minimum in der "Nähe" der Startwerte.
Gib auch mal in Google MATLAB "Optimization terminated" ein, da findest du mehr dazu, wie das zu interpretieren ist.
Die Toleranz bezieht sich auf die Änderung der Zielfunktion, nicht auf deren Wert.
Um zu sehen, wie gut die Optimierungsergebnisse sind, hilft oft auch schon ein Plot.
Zitat:
(ein Startwert 0.00034, ein anderer Startwert 0.00034)
Sorry, ich meinte ein Startwert gibt als RESNORM 0.0034, ein anderer 0.00034, also ein Faktor 10.
Von daher könnte man doch eigentlich erwarten, dass der mit der höhren RESNORM zumindest ein bisschen in die Gegend von dem anderen konvertiert (oder zumindest mal woanders hin konvertiert).
Das habe ich gemeint.
Oder ist meine Schrittweite zu klein oder zu groß?
Könnte das ein Problem sein?
% Input Data ist ein normierter Vektor, der die Daten der "zu fittenden" % Kurve enthält. % Um also das Verhalten von meinem Delay-Glied (Todzeit mit Verzögerung % dritter Ordnung vergleichen zu können, muss ich mein System mit einem % Puls der "Fläche 1" anregen, d.h. Pulshöhe = 1/timestep und Breite % timestep.
input_data = [0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;0;7.11044043571270e-05;0.000142209881996382;0.000213313481391913;0.000284418422390104;0.000355523363388296;0.000426628304386487;0.000497733245384678;0.000568838186382869;0.000639943127381060;0.000711048068379251;0.000759345764151608;0.000807643459923964;0.000855941155696320;0.000904238851468677;0.000952536547241033;0.00100083424301339;0.00104913193878575;0.00109742963455810;0.00114572733033046;0.00119402502610281;0.00110950808330917;0.00102498711570755;0.000940466148105925;0.000855945180504301;0.000771424212902678;0.000686903245301054;0.000602382277699431;0.000517861310097808;0.000433340342496184;0.000348819374894560;0.000402479456500309;0.000456143562914038;0.000509807669327767;0.000563471775741496;0.000617134540552565;0.000670798646966294;0.000724462753380023;0.000778126859793753;0.000831790966207482;0.000885455072621211;0.000889482563807561;0.000893507371788591;0.000897532179769621;0.000901556987750651;0.000905581795731680;0.000909606603712710;0.000913631411693740;0.000917656219674769;0.000921681027655799;0.000925705835636829;0.000947170136599660;0.000968635779165152;0.000990101421730643;0.00101156706429614;0.00103303270686163;0.00105449834942712;0.00107596533359527;0.00109743097616076;0.00111889661872625;0.00114036226129175;0.00110950540010385;0.00107864853891596;0.00104779167772806;0.00101693481654017;0.000986077955352274;0.000955221094164380;0.000924364232976485;0.000893507371788591;0.000862650510600697;0.000831794991015463;0.000885459097429192;0.000939123203842921;0.000992787310256650;0.00104645141667038;0.00110011686468677;0.00115378097110050;0.00120744507751423;0.00126110918392796;0.00131477329034169;0.00136844812957670;0.00149858358762999;0.00162871904568328;0.00175885450373658;0.00188898996178987;0.00201912541984316;0.00214926087789646;0.00227939633594975;0.00240953179400304;0.00253966725205634;0.00266981612613623;0.00297972634067552;0.00328963655521481;0.00359954676975409;0.00390947040031998;0.00421938061485927;0.00452929082939855;0.00483920104393784;0.00514911125847713;0.00545902147301641;0.00576887802344929;0.00562800974411325;0.00548714146477721;0.00534625976941456;0.00520539149007853;0.00506452321074249;0.00492365493140645;0.00478278665207041;0.00464191837273437;0.00450105009339833;0.00436020864611550;0.00432801018226726;0.00429581171841902;0.00426361325457079;0.00423141479072255;0.00419921632687431;0.00416701786302607;0.00413481939917783;0.00410260751930299;0.00407042247148136;0.00403821059160652;0.00397515526657039;0.00391209994153426;0.00384904461649812;0.00378598929146199;0.00372293396642586;0.00365987864138973;0.00359682331635360;0.00353376799131746;0.00347071266628133;0.00340764392521860;0.00332043975229629;0.00323323557937398;0.00314603140645167;0.00305882723352936;0.00297162306060705;0.00288441888768474;0.00279721471476243;0.00271001054184012;0.00262280636891781;0. 00253569610818172;0.00274901093117630;0.00296232575417087;0.00317564057716544;0.00338895540016002;0.00360227022315459;0.00381558504614916;0.00402889986914374;0.00424221469213831;0.00445552951513288;0.00466875042594123;0.00459630388228270;0.00452385733862416;0.00445141079496563;0.00437896425130710;0.00430651770764856;0.00423407116399003;0.00416162462033149;0.00408917807667296;0.00401673153301442;0.00394439231756872;0.00416709835918569;0.00438980440080267;0.00461251044241964;0.00483521648403662;0.00505792252565360;0.00528062856727057;0.00550333460888755;0.00572605406653113;0.00594876010814811;0.00617117099717981;0.00567209480753212;0.00517301861788444;0.00467394242823676;0.00417486623858908;0.00367577663291479;0.00317670044326711;0.00267762425361943;0.00217854806397175;0.00167947187432407;0.00752629701233930;0.00731566539466541;0.00710503377699152;0.00689440215931764;0.00668377054164375;0.00647313892396986;0.00626250730629597;0.00605187568862209;0.00584124407094820;0.00563061245327431;0.00541992717149401;0.00511940817557713;0.00481888917966024;0.00451837018374336;0.00421785118782648;0.00391733219190959;0.00361681319599271;0.00331629420007583;0.00301577520415894;0.00271525620824206;0.00241487137259121;0.00239206412736538;0.00236925688213954;0.00234644963691371;0.00232364239168787;0.00230083514646204;0.00227802790123620;0.00225522065601037;0.00223241341078453;0.00220960616555870;0.00218690624854569;0.00235728978640928;0.00252767332427287;0.00269805686213646;0.00286844040000005;0.00303882393786364;0.00320920747572723;0.00337959101359082;0.00354997455145441;0.00372035808931800;0.00389075504320819;0.00406516338905281;0.00423957173489743;0.00441398008074205;0.00458838842658667;0.00476279677243129;0.00493720511827591;0.00511161346412053;0.00528602180996515;0.00546043015580977;0.00563459701317553;0.00541994058752061;0.00520528416186570;0.00499062773621078;0.00477597131055586;0.00456131488490095;0.00434665845924603;0.00413200203359111;0.00391733219190959;0.00370267576625468;0.00348818033291900;0.00350562116750346;0.00352306200208793;0.00354050283667239;0.00355794367125685;0.00357538450584131;0.00359282534042577;0.00361026617501024;0.00362770700959470;0.00364514784417916;0.00366250818260400;0.00356725439371963;0.00347200060483526;0.00337674681595089;0.00328149302706652;0.00318623923818215;0.00309098544929779;0.00299573166041342;0.00290047787152905;0.00280522408264468;0.00270998370978691;0.00263887876878872;0.00256777382779053;0.00249666888679234;0.00242556394579415;0.00235445900479596;0.00228335406379776;0.00221224912279957;0.00214114418180138;0.00207003924080319;0.00199896113185821;0.00196676266800997;0.00193456420416173;0.00190236574031349;0.00187016727646526;0.00183796881261702;0.00180577034876878;0.00177357188492054;0.00174137342107231;0.00170917495722407;0.00167697649337583;0.00164477802952759;0.00161257956567936;0.00158038110183112;0.00154818263798288;0.00151598417413464;0.00148378571028641;0.00145158724643817;0.00141938878258993;0.00138719031874169;0.00135508576707968;0.00143692352936062;0.00151876129164155;0.00160059905392249;0.00168243681620343;0.00176427457848436;0.00184611234076530;0.00192793668701964;0.00200977444930057;0.00209161221158151;0.00217334264564962;0.00209418808868937;0.00201503353172912;0.00193587897476887;0.00185672441780862;0.00177756986084837;0.00169841530388812;0.00161927416295447;0.00154011960599422;0.00146096504903397;0.00138183732412692;0.00135768847624075;0.00133354096995723;0.00130939212207105;0.00128524327418487;0.00126109576790135;0.00123694692001518;0.00121279807212900;0.00118865056584548;0.00116450171795930;0.00114036226129175;0.00114036226129175;0.00114036226129175;0.00114036226129175;0.00114036226129175;0.00114036226129175;0.00114036226129175;0.00114036226129175;0.00114036226129175;0.00114036226129175;0.00114036628609973;0.00115512257376084;0.00116988020302462;0.00118463783228839;0.00119939546155217;0.00121415309081594;0.00122891072007972;0.00124366700774083;0.00125842463700461;0.00127318226626839;0.00128793318751886;0.00122487920408539;0.00116182522065192;0.00109877123721844;0.00103571591218231;0.000972661928748842;0.000909607945315370;0.000846552620279238;0.000783498636845767;0.000720444653412296;0.000657390669978824;0.000678849604531014;0.000700313905493845;0.000721779548059337;0.000743245190624829;0.000764710833190320;0.000786175134153152;0.000807640776718643;0.000829106419284135;0.000850572061849627;0.000872036362812458;0.000886795333678894;0.000901552962942670;0.000916310592206445;0.000931068221470221;0.000945825850733996;0.000960582138395111;0.000975339767658887;0.000990097396922662;0.00100485502618644;0.00101961265545021;0.00102095962452120;0.00102230122718154;0.00102364282984188;0.00102498443250223;0.00102632603516257;0.00102766763782291;0.00102900924048326;0.00103034950154094;0.00103169110420128;0.00103303270686163;0.00100085302545063;0.000968655903205057;0.000936457439356819;0.000904258975508582;0.000872061853263005;0.000839863389414767;0.000807664925566530;0.000775467803320952;0.000743269339472715;0.000711070875624477;0.000719093659533330;0.000727143275495389;0.000735192891457449;0.000743242507419508;0.000751290781778907;0.000759340397740966;0.000767390013703026;0.000775439629665085;0.000783489245627145;0.000791538861589204;0.000756023955964598;0.000720471485465502;0.000684919014966407;0.000649367886069972;0.000613815415570876;0.000578264286674441;0.000542711816175345;0.000507159345676250;0.000471608216779814;0.000436055746280719;0.000400504617384283;0.000364952146885188;0.000329399676386092;0.000293848547489657;0.000258296076990561;0.000222744948094126;0.000187192477595031;0.000151641348698595;0.000116088744039234;8.05368101812022e-05;8.45162719923123e-05;8.85410799733420e-05;9.25657537941056e-05;9.65905617751353e-05;0.000100615235595899;0.000104640043576929;0.000108664717397692;0.000112689525378722;0.000116714333359752;0.000120739007180515;0.000143521969398198;0.000166327873021373;0.000189135118247208;0.000211942363473042;0.000234749608698877;0.000257555512322052;0.000280362757547887;0.000303170002773722;0.000325977247999557;0.000348783151622731;0.000322023544959525;0.000295192833355321;0.000268360780148456;0.000241528726941592;0.000214698015337387;0.000187865962130523;0.000161033908923658;0.000134201855716794;0.000107370741631791;8.05390909057248e-05;8.71933059407612e-05;9.39011850822113e-05;0.000100609198383927;0.000107317077525378;0.000114024956666828;0.000120732835808278;0.000127440849109994;0.000134148728251444;0.000140856204912096;0.000147564218213812;0.000152933312060506;0.000158299722701879;0.000163666133343252;0.000169032543984625;0.000174398954625997;0.000179765365267370;0.000185131775908743;0.000190496844947456;0.000195863255588829;0.000201229666230202;0.000187849862898599;0.000174433836295166;0.000161019151294394;0.000147603124690962;0.000134187098087530;0.000120771071484098;0.000107355179040931;9.39394207580311e-05;8.05235283148648e-05;6.71077700319646e-05;7.64518984009891e-05;8.58429828631257e-05;9.52340673252623e-05;0.000104625151787399;0.000114016102089269;0.000123407186551406;0.000132798271013543;0.000142189757956477;0.000151580976578880;0.000160972195201282;0.000170362072221025;0.000179753290843427;0.000189144509465830;0.000198535728088233;0.000207926946710635;0.000217318165333038;0.000226709383955440;0.000236100602577843;0.000245491821200246;0.000254881698219988;0.000254012339696085;0.000253117490721637;0.000252223983349848;0.000251329134375399;0.000250434285400950;0.000249540778029161;0.000248645929054713;0.000247751080080264;0.000246857572708475;0.000245962723734026;0.000245067874759577;0.000244174367387789;0.000243279518413340;0.000242384669438891;0.000241491162067102;0.000240596313092653;0.000239701464118204;0.000238807956746416;0.000237913107771967;0.000237018258797518;0.000236124751425729;0.000235229902451280;0.000234335053476831;0.000233441546105043;0.000232546697130594;0.000231653189758805;0.000230758340784356;0.000229863491809907;0.000228969984438119;0.000228075135463670;0.000227180286489221;0.000226286779117432;0.000225391930142983;0.000224497081168534;0.000223603573796746;0.000222708724822297;0.000221813875847848;0.000220920368476059;0.000220025519501610;0.000219130670527161;0.000218237163155373;0.000217342314180924;0.000216447465206475;0.000215553957834686;0.000214659108860237;0.000213764259885789;0.000212870752514000;0.000211975903539551;0.000211081054565102;0.000210187547193313;0.000209292698218865;0.000208397849244416;0.000207504341872627;0.000206609492898178;0.000205714643923729;0.000204821136551941;0.000203926287577492;0.000203031438603043;0.000202137931231254;0.000201243082256805;0.000200349574885017;0.000199454725910568;0.000198559876936119;0.000197666369564330;0.000196771520589881;0.000195876671615432;0.000194983164243644;0.000194088315269195;0.000193193466294746;0.000192299958922957;0.000191405109948508;0.000190510260974059;0.000189616753602271;0.000188721904627822;0.000187827055653373;0.000186933548281584;0.000186038699307135;0.000185143850332686;0.000184250342960898;0.000183355493986449;0.000182460645012000;0.000181567137640211;0.000180672288665762;0.000179777439691314;0.000178883932319525;0.000177989083345076;0.000177094234370627;0.000176200726998838;0.000175305878024390;0.000174411029049941;];
%--------------
final_time = 62;
time_step = 0.1;
dead_time_start = 7; %Startwert für Todzeit
D_start = 16.25; % Startwert für Verzögerung
t=0:0.1:final_time;
t=t';
input_pulse=sum(input_data)/time_step;
simin=[t zeros(length(t),1)];
simin(2,2)=input_pulse; %Ich will ja nur einen Puls am Anfang haben
parameter_start=[dead_time_start D_start];
plot(t,simin(:,2))
output_des = input_data; % Das gegebene Signal soll "gefittet" werden.
hat irgendjemand von euch eine Idee, wieso MATLAB so reagiert bei dem Code?
Würde mich über Antworten freuen.
Liebe Grüße
Mattes
Einstellungen und Berechtigungen
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
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.
)
RSS