wobei R ein konstanter (3 x 1)-Vektor und VC eine konstante (3 x 3)-Matrix ist (beide mit Werten zwischen 0 und 1). EW steht für einen Erwartungswert, den ich über folgende Schleife abzubilden versuche:
Code:
n = 100000;
z = zeros(3,1);
for i=1:n
J = ( lognrnd(M, N) - 1) * (1 + x' * ( lognrnd(M, N) - 1))^(-3);
z = [ z(1,1) + J(1,1);
z(2,1) + J(2,1);
z(3,1) + J(3,1)];
end
ew = [ z(1,1) / n;
z(2,1) / n;
z(3,1) / n];
M und N sind wieder (3 x 1)-Vektoren mit Werten zwischen 0 und 1.
Wenn ich nun versuche, diese Konstruktion mit fsolve nach x aufzulösen (und x0 = [ 1; 1; 1 ] vorgebe), dann wird zwar die Optimierung durchlaufen, jedoch ohne x zu verändern!!!
Beim Versuch das Problem einzukreisen, glaube ich schon festgestellt zu haben, dass das Problem die Kombination aus Schleife, x-Werten und Zufallszahlen sein dürfte. Aber ohne die kann ich den Erwartungswert nicht bestimmen.
Ich bin echt mit meinem Latein am Ende.
Hat jemand von Euch vielleicht eine Idee???
P.S.: Keine Ahnung, ob ich das richtige Forum für mein Problem erwischt habe, ansonsten bitte den Thread einfach verschieben.
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