bmwkraftur.is

Er með eitt dæmi sem ég er búinn að berjast við í c.a. 15klst og er að spá hvort einhver geti aðstoðað mig (Ef það heppnast þá er bjórkippa í boði(eða andvirði hennar)).

Þetta snýst semsagt um að skrifa kóða til að leysa kerfi af quadratic Riccati diffurjöfnum á forminu:

dX/dt = A'X + 1/2*B'*diag(C'X)*C'X -g

dx/dt = a'X +1/2*b'*diag(C'X)*C'X -g0

Endilega sendið mér PM ef þið getið bjargað mér,, er að verða gráhærður

Árni Björn ætti að geta reddað þér á no time fyrst það er bjór í boði

_________________
Fullt af CF sem rúllar

Haha ég las yfir þetta og ákvað svo að þetta væri ekki þess virði. Fékk hausverk bara við að lesa þetta

Ekki það að ætti breik í að hjálpa þér með þetta Jónas

_________________
Enginn bíll!

Hef ekki notad Matlab mikid ...

Datt nidur a thetta. Getur thu eitthvad notad kodan a bak vid thetta sem grunn?

http://ideas.repec.org/c/dge/qmrbcd/21.html

_________________
BMW E46 ///M3 04.2003
Land Rover Discovery 3 04.2007
Ski Doo Summit X 800 151" - 2006
5xE30,E32,E34,4xE36,E39,E53 - Seldir

Þetta á að vera vel þekkt jafna sem þarf að leysa í sambandi við "affine diffusion model" til að verðleggja hlutabréf og eignir, en ég hef ekkert fundið um þessa Riccati jöfnu, eina sem ég hef fundið hingað til er hvernig á að stjórna halla á gervihnöttum... en málið er það að ég er í hagfræði, ekki verkfræði

Það eru verða 3ár síðan ég notaði Matlab síðast og ég er ekki með Matlab uppsett á neinni nothæfri vél í dag svo mín hjálp er nánast einskis virði.

En ef ég ætti að leysa þetta myndi ég byrja á að lesa um ode föllin til að leysa diffurjöfnur í Matlab og svo reyna að útfæra það á Riccatti í framhaldi.

_________________
E36 M3 US | Jeep Grand Cherokee Limited V8

http://www.mathworks.com/support/soluti ... on=1-1C6YS

Notaðist við þetta og fékk:

function dX2dt = mRiccati2(t, X, A, B, C,g)
X = reshape(X, size(A));
dX2dt = A.'*X + 1/2*B.'*((C.'*X)*(C.'*X)) - g
dX2dt = dX2dt(:);

Smelli svo inn tölunum sem gefnar eru:

A = [-17.4 17.4 -9.309; 0 -0.226 0.879; 0 0 -0.362];
B = [0 0 1; 0 0 0; 0 0 0.00782];
C = [1 -3.42 4.27; -0.0943 1 0; 0 0 1];
g = [1;0;0];
X0 = zeros(size(A));

Og fæ:
>> [T X] = ode45(@mRiccati2,[0 30], X0, [], A,B,C,g)
??? Error using ==> minus
Matrix dimensions must agree.

en ef ég breyti aðeins og set g = 0 (scalar) þá fæ ég þetta til að virka en allar afleiður eru = 0 og X er alltaf 0...

Búinn að skoða þetta?: http://www.mathworks.com/support/soluti ... on=1-1C6YS

_________________
Þórður
'99 M5 SC // '89 M3 S85 // '04 Patrol GR

Benti einmitt á þetta í fyrri póstinum

Úps... var ekki kominn þegar ég replyaði....

_________________
Þórður
'99 M5 SC // '89 M3 S85 // '04 Patrol GR

Kein prob, fæ þessa villu að ég geti ekki dregið g vektorinn fá vegna dimensions, en kóðinn virðist vera réttur hrmpft

Mér nú veit ég hvernig Alpina líður þegar hann les tölvu troubleshooting þræði...

_________________
Enginn BMW

Svarið er X

_________________
BMW E39 530D '03 ///M-Tech - Stations
BMW E39 540i '00 - Sedans

E39 k1ng
Sævar P.

Vá þetta er hebreska fyrir mér.

_________________
Axel Jóhann Helgason. S. 695-7205


E34 Touring
Musso 2.9TDi

Nei það er klárlega 43 Svarið er ALLTAF 43