|
Revision 706, 0.7 kB
(checked in by smidl, 16 years ago)
|
|
eol-native
|
-
Property svn:eol-style set to
native
|
| Line | |
|---|
| 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; |
|---|
| 15 | C2.b0 = 0; |
|---|
| 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 |
|---|
| 23 | c.b = randn(); |
|---|
| 24 | c.seed = seeds(i); |
|---|
| 25 | c.controller = C1; |
|---|
| 26 | Mmc=iterativemc(c); |
|---|
| 27 | losses(i) = sum((Mmc.y-yr*ones(size(Mmc.y))).^2); |
|---|
| 28 | |
|---|
| 29 | c.controller = C2; |
|---|
| 30 | Mmc=iterativemc(c); |
|---|
| 31 | losses2(i) = sum((Mmc.y-yr*ones(size(Mmc.y))).^2); |
|---|
| 32 | end |
|---|
| 33 | [min(losses) median(losses) max(losses)] |
|---|
| 34 | [min(losses2) median(losses2) max(losses2)] |
|---|
| 35 | |
|---|
