| 1 | function losses=mc_study(regulator_parameters,system,range,n)
|
|---|
| 2 | [H u]=load_H;
|
|---|
| 3 |
|
|---|
| 4 | for i=1:n
|
|---|
| 5 | apriori(3)=range.P0;
|
|---|
| 6 | system.b=range.b0 + apriori(3)*randn;
|
|---|
| 7 | apriori(2)=range.b0 + apriori(3)*randn;
|
|---|
| 8 | apriori(1)=range.y0 + range.y0_range*randn;
|
|---|
| 9 | while (system.b==0)
|
|---|
| 10 | system.b=range.b0 + apriori(3)*randn;
|
|---|
| 11 | end
|
|---|
| 12 | while (apriori(2)==0)
|
|---|
| 13 | apriori(2)=range.b0 + apriori(3)*randn;
|
|---|
| 14 | end
|
|---|
| 15 | seed=randn(1,size(H,1)+2);
|
|---|
| 16 | j=1;
|
|---|
| 17 | losses(j,i)=ce(regulator_parameters,system,apriori,seed); j=j+1;
|
|---|
| 18 | losses(j,i)=cc(regulator_parameters,system,apriori,seed);j=j+1;
|
|---|
| 19 | %losses(j,i)=ps(regulator_parameters,system,apriori,seed);j=j+1;
|
|---|
| 20 | losses(j,i)=rizeni(H,u,system,apriori,seed);j=j+1;
|
|---|
| 21 | %losses(j,i)=triv(regulator_parameters,system,apriori,seed);j=j+1;
|
|---|
| 22 | losses(j,i)=alstr(regulator_parameters,system,apriori,seed);j=j+1;
|
|---|
| 23 |
|
|---|
| 24 | end
|
|---|
| 25 | end |
|---|
