root/applications/dual/SIDP/SIDPS/matlab/generate_realization.m
@
1183
Revision 1181, 414 bytes (checked in by zimamiro, 14 years ago) |
---|
Line | |
---|---|
1 | function realization=generate_realization(eta,index,H,C,horizont) |
2 | realization=0; |
3 | beta=H(index); |
4 | for o=1:horizont |
5 | best_control=-(eta*beta/(1+beta^2)+interpol(beta,H,C)); |
6 | %best_control=-(eta*beta/(1+beta^2)+C(index)*abs(rand)); |
7 | s=randn; |
8 | pom=sqrt(1+best_control^2); |
9 | eta=abs(eta+beta*best_control+pom*s); |
10 | beta=abs(pom*beta+best_control*s); |
11 | realization=realization+eta^2; |
12 | end |
13 | end |
Note: See TracBrowser
for help on using the browser.