Revision 661, 0.7 kB
(checked in by smidl, 15 years ago)
|
doc
|
Rev | Line | |
---|
[650] | 1 | % configuracni struktura se bude jmenovat "c" |
---|
| 2 | c.b =1; |
---|
| 3 | c.sigma = 0.1; |
---|
| 4 | c.ndat = 50; |
---|
| 5 | yr = 1; |
---|
| 6 | |
---|
| 7 | % nejlepsi mozny regulator - zna b |
---|
| 8 | C1.class='exact_ctrl'; |
---|
| 9 | C1.yr = yr; |
---|
| 10 | C1.b = c.b; |
---|
| 11 | |
---|
| 12 | % CE regulator - b se odhaduje |
---|
| 13 | C2.class='ce_ctrl'; |
---|
| 14 | C2.yr = yr; |
---|
[661] | 15 | C2.b0 = 0; |
---|
[650] | 16 | C2.P0 = 1; |
---|
| 17 | |
---|
| 18 | % monte carlo study |
---|
| 19 | n=100; |
---|
| 20 | losses=zeros(10,1); |
---|
| 21 | seeds=32000*rand(1,n); |
---|
| 22 | for i=1:n |
---|
[661] | 23 | c.b = randn(); |
---|
[650] | 24 | c.seed = seeds(i); |
---|
[661] | 25 | c.controller = C1; |
---|
[650] | 26 | Mmc=iterativemc(c); |
---|
| 27 | losses(i) = sum((Mmc.y-yr*ones(size(Mmc.y))).^2); |
---|
[661] | 28 | |
---|
| 29 | c.controller = C2; |
---|
| 30 | Mmc=iterativemc(c); |
---|
| 31 | losses2(i) = sum((Mmc.y-yr*ones(size(Mmc.y))).^2); |
---|
[650] | 32 | end |
---|
[661] | 33 | [min(losses) median(losses) max(losses)] |
---|
| 34 | [min(losses2) median(losses2) max(losses2)] |
---|
| 35 | |
---|