Guten Tag alle zusammen!
Ich bin Matlab-Neuling und soll folgende Aufgabe analytisch - im Idealfall nicht numerisch und ohne Schleifen - lösen:
Ich bin Wissenschaftler und meine Inputvektoren sind 3x1 - Vektoren (Messdaten).
Nun möchte ich ein lineares Gleichungssystem nach a und b auflösen, unter Bedingungen:
Für einen Wert funktioniert es: Dabei ist a in der Erstberechnung ohne Bedingungen negativ und b positiv. Dann bekomme ich brauchbare Werte.
Allerdings für den umgekehrten Fall (b negativ) funktioniert es nicht! ('epmty sym: 0-by-1').
D.h. ich konnte den solver bisher nur erfolgreich für einen Parameter anwenden, aber nicht für die gesamte Matrix. Deswegen jeweils eine Zeile für p1, p2, p3.
Ich stecke hier ziemlich fest und bin dankbar für jeden Hinweis!
Muss ich evtl einen anderen Solver nehmen, um eine simples lineares Gleichungssystem zu lösen?
Ja, das war zuerst mein Plan, numerisch und mit Matrizen.
Ich sollte das analytisch lösen und der Solver gibt nur leere structs zurück, wenn ich mit der kompletten Matrix löse. :-/
Aber ich werds nochmal mit lsqlin probieren!
_________________
Isqlin eignet sich von Grund auf her für dieses Regressionsproblem.
Allerdings habe ich constraints eingeben (wh nicht richtig).
Für ein Parameter kommt 0 heraus, obwohl beide Ergebniswerte größer als 0 sein sollten.
isqlin
% numerical solver: lsqlin
% one parameter per line/row
C = [source_10_measurements_ini source_11_measurements_ini ];
d = drain_12_measurements_ini;
A = [0,1; 0,1];
b = [1,1];
Aeq = [0,0; 0,0];
beq = [0,0];
lb = [0;0];
ub = [1;1];
x = lsqlin(C,d,A,b,Aeq,beq,lb,ub);
_________________
Mit strengen Ungleichungen sind gemeint, dass die Werte z.B. nicht null annehmen können. Dies ist jedoch im Solver nicht vorgesehen.
Der 'interior-point' - Algorithmus ist per default eingestellt. D.h. wenn options in ' x = lsqlin(C,d,A,b,Aeq,beq,lb,ub,x0,options)' nicht definiert ist, ist er default-Algorithmus aktiv.
Genau!
EndMembers und d sollten unendlich viele Elemente bzw Spalten haben:
EndMembers = [A; B; C] (drei Zeilen)
und d wie EndMembers in dem Fall neun Spalten,
d.h. length = 9
Wenn ich Isqlin ausführe, erscheint die Meldung:
'Error using lsqlin (line 216)
The number of rows in C must be equal to the length of d.'
Warum bloß?!
Die Bedingungen sind, dass alle 3 Elemente von x positiv sind und zusammen 1 ergeben (100%).
Bei mir erscheint zudem folgende Meldung:
"Error using vertcat
Dimensions of matrices being concatenated are not consistent." wegen [A; B; C]
Hier nochmal verständlicher:
Code:
A = [21.87390.09861.38841.19220.68160.27480.19080.12910.0887];
B = [20.96760.09751.46421.15940.65850.25950.17560.11700.0775];
C = [18.85160.08961.73511.08430.60420.21470.13320.07600.0416];
%field value
D = [21.45070.09751.39521.16340.66100.26470.18070.11890.0799];
%% Isqlin % numerical solver: lsqlin % one parameter per line/row % L = length(D); % L = 9
EndMembers = [A; B; C];
d = D;
A = [];
b = [];
Aeq = [1;1;1];
beq = 1;
lb = [0;0;0];
ub = [1;1;1];
x = lsqlin(EndMembers,d,A,b,Aeq,beq,lb,ub);
% The number of rows in 'EndMembers' must be equal to the length of d. disp(x)
Ok, habe das wie folgt geschrieben, nach mehreren Fehlermeldungen:
Code:
A = [21.87390.09861.38841.19220.68160.27480.19080.12910.0887];
B = [20.96760.09751.46421.15940.65850.25950.17560.11700.0775];
C = [18.85160.08961.73511.08430.60420.21470.13320.07600.0416];
%field value
D = [21.45070.09751.39521.16340.66100.26470.18070.11890.0799];
EndMembers = [A;B;C]';
d = D;
A = [];
b = [];
Aeq = [111];
beq = 1;
lb = [0;0;0];
ub = [1;1;1];
x = lsqlin(EndMembers,d,A,b,Aeq,beq,lb,ub);
% The number of rows in EndMembers must be equal to the length of d. disp(x)
Optimization completed because the objective function is non-decreasing in
feasible directions, to within the default value of the optimality tolerance,
and constraints are satisfied to within the default value of the constraint tolerance.
<stopping criteria details>
0.5346
0.4654
0.0000
"
Die Summe ist zwar 1, jedoch ist C (= Inputvektor) auf gar keinen Fall 0.
Edit: C ist nicht 0, sondern 1.4507e-07
_________________
hmm, mir ist nicht wirklich eine bessere berechnete Lösung bekannt und wenn, kann ich es noch nicht beurteilen. Da mir bspw ein tieferes Verständnis der Unterschiede der Funktionen Isqlin, Regress und Regression fehlt.
Es soll mein Mischungsmodell widerspiegeln (geht hier nicht um Konzentrationen oder logaritmische Größen!), deswegen sollten alle Komponenten messbar sein, auch wenn verschwindend gering.
Das führt mich auch grad zum nächsten Punkt: Alle wissenschaftlichen Daten sind nicht validiert bzw. zu veröffentlichen, wenn die Ungenauigkeit bzw. der Fehler nicht definiert wurde: In dem Fall 3 % von jedem Meßwert - jedem Wert in den Inputvektoren A, B und C.
Wie kann ich diese berücksichtigen` bzw. als extra Matrix miteinbinden?
Wird das durch A und b definiert?!
--> optimality tolerance?
--> constraint tolerance?
wenn alle Messgrößen gleichermaßen fehlerbehaftet sind, würde ich erwarten, dass sich diese Fehler herausmitteln, was die zu erwartenden Ergebnisse angeht. Klar ist aber auch, dass dies das Ergebnis verfälschen könnte, wenn die Messfehler beispielsweise die Beiträge einer Komponente systematisch über- oder unterschätzen.
Man könnte verschiedene (im Zweifelsfall zufällige) Störungen auf die Messungen legen und dann analysieren, wie unterschiedliche Messungen die Ergebnisse beeinflussen. Man bekäme quasi eine Verteilung der Ergebnisse.
Nein, tatsächlich sind nicht alle Messgrößen gleichermaßen fehlbehaftet. Habe das Beispiel nur aus Vereinfachungsgründen genannt:
Tatsächlich ist es eine Matrix, die jeder Messgröße also jedem Wert von den Inputvektoren einen Fehler zuordnet: Standardabweichungen der im Labor ermittelten Werte.
Wenn man das aber in A schreibt, würde der Vektor so groß sein wie C [https://de.mathworks.com/help/optim/ug/lsqlin.html?requestedDomain=www.mathworks.com]
Klar, dass das Ergebnis dann anders ist.
Viele Grüße
Patlab
_________________
Socratrees: Ich weiss, dass ich nichts weiss.
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