| 1 | function los=vyst2(regulator_parameters,system,apriori,N) | 
|---|
| 2 | los=zeros(4,3); | 
|---|
| 3 |  | 
|---|
| 4 | %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% | 
|---|
| 5 | regulator_parameters.horizont=5; | 
|---|
| 6 | apriori.P0=1; | 
|---|
| 7 | apriori.b0=0; | 
|---|
| 8 | [m los]=mc_study(regulator_parameters,system,apriori,N,[regulator_parameters.horizont apriori.P0 apriori.b0]); | 
|---|
| 9 |  | 
|---|
| 10 | figure | 
|---|
| 11 | hist(min(los(1,:),10),100);    AXIS([0 10 0 150]);   text(8,120,'\fontsize{18} CE'); | 
|---|
| 12 | title('P_0 = 1, {}_{^{\fontsize{18}\theta}}^{\fontsize{16}\^}_{0} = 0'); set( get(gcf, 'Children'), 'FontSize', 18);set(get(get(gcf, 'Children'),'Title'), 'FontSize', 18); | 
|---|
| 13 | text(9.4,-8,'\fontsize{18}<') | 
|---|
| 14 |  | 
|---|
| 15 | figure | 
|---|
| 16 | hist(min(los(2,:),10),100);    AXIS([0 10 0 150]);   text(8,120,'\fontsize{18} CC') | 
|---|
| 17 | title('P_0 = 1, {}_{^{\fontsize{18}\theta}}^{\fontsize{16}\^}_{0} = 0');set( get(gcf, 'Children'), 'FontSize', 18);set(get(get(gcf, 'Children'),'Title'), 'FontSize', 18); | 
|---|
| 18 | text(9.4,-8,'\fontsize{18}<') | 
|---|
| 19 |  | 
|---|
| 20 | figure | 
|---|
| 21 | hist(min(los(3,:),10),100);    AXIS([0 10 0 150]);   text(8,120,'\fontsize{18} SIDP') | 
|---|
| 22 | title('P_0 = 1, {}_{^{\fontsize{18}\theta}}^{\fontsize{16}\^}_{0} = 0');set( get(gcf, 'Children'), 'FontSize', 18);set(get(get(gcf, 'Children'),'Title'), 'FontSize', 18); | 
|---|
| 23 | text(9.4,-8,'\fontsize{18}<') | 
|---|
| 24 |  | 
|---|
| 25 | figure | 
|---|
| 26 | hist(min(los(4,:),10),100);    AXIS([0 10 0 150]);   text(8,120,'\fontsize{18} DP') | 
|---|
| 27 | title('P_0 = 1, {}_{^{\fontsize{18}\theta}}^{\fontsize{16}\^}_{0} = 0');set( get(gcf, 'Children'), 'FontSize', 18);set(get(get(gcf, 'Children'),'Title'), 'FontSize', 18); | 
|---|
| 28 | text(9.4,-8,'\fontsize{18}<') | 
|---|
| 29 |  | 
|---|
| 30 | apriori.b0=1; | 
|---|
| 31 | [m los]=mc_study(regulator_parameters,system,apriori,N,[regulator_parameters.horizont apriori.P0 apriori.b0]); | 
|---|
| 32 | text(9.4,-8,'\fontsize{18}<') | 
|---|
| 33 |  | 
|---|
| 34 | figure | 
|---|
| 35 | hist(min(los(1,:),10),100);    AXIS([0 10 0 150]);   text(8,120,'\fontsize{18} CE'); | 
|---|
| 36 | title('P_0 = 1, {}_{^{\fontsize{18}\theta}}^{\fontsize{16}\^}_{0} = 1');set( get(gcf, 'Children'), 'FontSize', 18);set(get(get(gcf, 'Children'),'Title'), 'FontSize', 18); | 
|---|
| 37 | text(9.4,-8,'\fontsize{18}<') | 
|---|
| 38 |  | 
|---|
| 39 | figure | 
|---|
| 40 | hist(min(los(2,:),10),100);    AXIS([0 10 0 150]);   text(8,120,'\fontsize{18} CC') | 
|---|
| 41 | title('P_0 = 1, {}_{^{\fontsize{18}\theta}}^{\fontsize{16}\^}_{0} = 1');set( get(gcf, 'Children'), 'FontSize', 18);set(get(get(gcf, 'Children'),'Title'), 'FontSize', 18); | 
|---|
| 42 | text(9.4,-8,'\fontsize{18}<') | 
|---|
| 43 |  | 
|---|
| 44 | figure | 
|---|
| 45 | hist(min(los(3,:),10),100);    AXIS([0 10 0 150]);   text(8,120,'\fontsize{18} SIDP') | 
|---|
| 46 | title('P_0 = 1, {}_{^{\fontsize{18}\theta}}^{\fontsize{16}\^}_{0} = 1');set( get(gcf, 'Children'), 'FontSize', 18);set(get(get(gcf, 'Children'),'Title'), 'FontSize', 18); | 
|---|
| 47 | text(9.4,-8,'\fontsize{18}<') | 
|---|
| 48 |  | 
|---|
| 49 | figure | 
|---|
| 50 | hist(min(los(4,:),10),100);    AXIS([0 10 0 150]);   text(8,120,'\fontsize{18} DP') | 
|---|
| 51 | title('P_0 = 1, {}_{^{\fontsize{18}\theta}}^{\fontsize{16}\^}_{0} = 1');set( get(gcf, 'Children'), 'FontSize', 18);set(get(get(gcf, 'Children'),'Title'), 'FontSize', 18); | 
|---|
| 52 | text(9.4,-8,'\fontsize{18}<') | 
|---|
| 53 |  | 
|---|
| 54 | apriori.b0=10; | 
|---|
| 55 | [m los]=mc_study(regulator_parameters,system,apriori,N,[regulator_parameters.horizont apriori.P0 apriori.b0]); | 
|---|
| 56 | text(9.4,-8,'\fontsize{18}<') | 
|---|
| 57 |  | 
|---|
| 58 | figure | 
|---|
| 59 | hist(min(los(1,:),10),100);    AXIS([0 10 0 150]);    text(8,120,'\fontsize{18} CE'); | 
|---|
| 60 | title('P_0 = 1, {}_{^{\fontsize{18}\theta}}^{\fontsize{16}\^}_{0} = 10');set( get(gcf, 'Children'), 'FontSize', 18);set(get(get(gcf, 'Children'),'Title'), 'FontSize', 18); | 
|---|
| 61 | text(9.4,-8,'\fontsize{18}<') | 
|---|
| 62 |  | 
|---|
| 63 | figure | 
|---|
| 64 | hist(min(los(2,:),10),100);    AXIS([0 10 0 150]);   text(8,120,'\fontsize{18} CC') | 
|---|
| 65 | title('P_0 = 1, {}_{^{\fontsize{18}\theta}}^{\fontsize{16}\^}_{0} = 10');set( get(gcf, 'Children'), 'FontSize', 18);set(get(get(gcf, 'Children'),'Title'), 'FontSize', 18); | 
|---|
| 66 | text(9.4,-8,'\fontsize{18}<') | 
|---|
| 67 |  | 
|---|
| 68 | figure | 
|---|
| 69 | hist(min(los(3,:),10),100);    AXIS([0 10 0 150]);   text(8,120,'\fontsize{18} SIDP') | 
|---|
| 70 | title('P_0 = 1, {}_{^{\fontsize{18}\theta}}^{\fontsize{16}\^}_{0} = 10');set( get(gcf, 'Children'), 'FontSize', 18);set(get(get(gcf, 'Children'),'Title'), 'FontSize', 18); | 
|---|
| 71 | text(9.4,-8,'\fontsize{18}<') | 
|---|
| 72 |  | 
|---|
| 73 | figure | 
|---|
| 74 | hist(min(los(4,:),10),100);    AXIS([0 10 0 150]);   text(8,120,'\fontsize{18} DP') | 
|---|
| 75 | title('P_0 = 1, {}_{^{\fontsize{18}\theta}}^{\fontsize{16}\^}_{0} = 10');set( get(gcf, 'Children'), 'FontSize', 18);set(get(get(gcf, 'Children'),'Title'), 'FontSize', 18); | 
|---|
| 76 | text(9.4,-8,'\fontsize{18}<') | 
|---|
| 77 |  | 
|---|
| 78 | %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% | 
|---|
| 79 | apriori.P0=10; | 
|---|
| 80 | apriori.b0=0; | 
|---|
| 81 | [m los]=mc_study(regulator_parameters,system,apriori,N,[regulator_parameters.horizont apriori.P0 apriori.b0]); | 
|---|
| 82 | text(9.4,-8,'\fontsize{18}<') | 
|---|
| 83 |  | 
|---|
| 84 | figure | 
|---|
| 85 | hist(min(los(1,:),10),100);    AXIS([0 10 0 150]);   text(8,120,'\fontsize{18} CE'); | 
|---|
| 86 | title('P_0 = 10, {}_{^{\fontsize{18}\theta}}^{\fontsize{16}\^}_{0} = 0');set( get(gcf, 'Children'), 'FontSize', 18);set(get(get(gcf, 'Children'),'Title'), 'FontSize', 18); | 
|---|
| 87 | text(9.4,-8,'\fontsize{18}<') | 
|---|
| 88 |  | 
|---|
| 89 | figure | 
|---|
| 90 | hist(min(los(2,:),10),100);    AXIS([0 10 0 150]);   text(8,120,'\fontsize{18} CC') | 
|---|
| 91 | title('P_0 = 10, {}_{^{\fontsize{18}\theta}}^{\fontsize{16}\^}_{0} = 0');set( get(gcf, 'Children'), 'FontSize', 18);set(get(get(gcf, 'Children'),'Title'), 'FontSize', 18); | 
|---|
| 92 | text(9.4,-8,'\fontsize{18}<') | 
|---|
| 93 |  | 
|---|
| 94 | figure | 
|---|
| 95 | hist(min(los(3,:),10),100);    AXIS([0 10 0 150]);   text(8,120,'\fontsize{18} SIDP') | 
|---|
| 96 | title('P_0 = 10, {}_{^{\fontsize{18}\theta}}^{\fontsize{16}\^}_{0} = 0');set( get(gcf, 'Children'), 'FontSize', 18);set(get(get(gcf, 'Children'),'Title'), 'FontSize', 18); | 
|---|
| 97 | text(9.4,-8,'\fontsize{18}<') | 
|---|
| 98 |  | 
|---|
| 99 | figure | 
|---|
| 100 | hist(min(los(4,:),10),100);    AXIS([0 10 0 150]);   text(8,120,'\fontsize{18} DP') | 
|---|
| 101 | title('P_0 = 10, {}_{^{\fontsize{18}\theta}}^{\fontsize{16}\^}_{0} = 0');set( get(gcf, 'Children'), 'FontSize', 18);set(get(get(gcf, 'Children'),'Title'), 'FontSize', 18); | 
|---|
| 102 | text(9.4,-8,'\fontsize{18}<') | 
|---|
| 103 |  | 
|---|
| 104 | apriori.b0=1; | 
|---|
| 105 | [m los]=mc_study(regulator_parameters,system,apriori,N,[regulator_parameters.horizont apriori.P0 apriori.b0]); | 
|---|
| 106 |  | 
|---|
| 107 | figure | 
|---|
| 108 | hist(min(los(1,:),10),100);    AXIS([0 10 0 150]);   text(8,120,'\fontsize{18} CE'); | 
|---|
| 109 | title('P_0 = 10, {}_{^{\fontsize{18}\theta}}^{\fontsize{16}\^}_{0} = 1');set( get(gcf, 'Children'), 'FontSize', 18);set(get(get(gcf, 'Children'),'Title'), 'FontSize', 18); | 
|---|
| 110 | text(9.4,-8,'\fontsize{18}<') | 
|---|
| 111 |  | 
|---|
| 112 | figure | 
|---|
| 113 | hist(min(los(2,:),10),100);    AXIS([0 10 0 150]);   text(8,120,'\fontsize{18} CC') | 
|---|
| 114 | title('P_0 = 10, {}_{^{\fontsize{18}\theta}}^{\fontsize{16}\^}_{0} = 1');set( get(gcf, 'Children'), 'FontSize', 18);set(get(get(gcf, 'Children'),'Title'), 'FontSize', 18); | 
|---|
| 115 | text(9.4,-8,'\fontsize{18}<') | 
|---|
| 116 |  | 
|---|
| 117 | figure | 
|---|
| 118 | hist(min(los(3,:),10),100);    AXIS([0 10 0 150]);   text(8,120,'\fontsize{18} SIDP') | 
|---|
| 119 | title('P_0 = 10, {}_{^{\fontsize{18}\theta}}^{\fontsize{16}\^}_{0} = 1');set( get(gcf, 'Children'), 'FontSize', 18);set(get(get(gcf, 'Children'),'Title'), 'FontSize', 18); | 
|---|
| 120 | text(9.4,-8,'\fontsize{18}<') | 
|---|
| 121 |  | 
|---|
| 122 | figure | 
|---|
| 123 | hist(min(los(4,:),10),100);    AXIS([0 10 0 150]);   text(8,120,'\fontsize{18} DP') | 
|---|
| 124 | title('P_0 = 10, {}_{^{\fontsize{18}\theta}}^{\fontsize{16}\^}_{0} = 1');set( get(gcf, 'Children'), 'FontSize', 18);set(get(get(gcf, 'Children'),'Title'), 'FontSize', 18); | 
|---|
| 125 | text(9.4,-8,'\fontsize{18}<') | 
|---|
| 126 |  | 
|---|
| 127 | apriori.b0=10; | 
|---|
| 128 | [m los]=mc_study(regulator_parameters,system,apriori,N,[regulator_parameters.horizont apriori.P0 apriori.b0]); | 
|---|
| 129 |  | 
|---|
| 130 | figure | 
|---|
| 131 | hist(min(los(1,:),10),100);    AXIS([0 10 0 150]);   text(8,120,'\fontsize{18} CE'); | 
|---|
| 132 | title('P_0 = 10, {}_{^{\fontsize{18}\theta}}^{\fontsize{16}\^}_{0} = 10');set( get(gcf, 'Children'), 'FontSize', 18);set(get(get(gcf, 'Children'),'Title'), 'FontSize', 18); | 
|---|
| 133 | text(9.4,-8,'\fontsize{18}<') | 
|---|
| 134 |  | 
|---|
| 135 | figure | 
|---|
| 136 | hist(min(los(2,:),10),100);    AXIS([0 10 0 150]);   text(8,120,'\fontsize{18} CC') | 
|---|
| 137 | title('P_0 = 10, {}_{^{\fontsize{18}\theta}}^{\fontsize{16}\^}_{0} = 10');set( get(gcf, 'Children'), 'FontSize', 18);set(get(get(gcf, 'Children'),'Title'), 'FontSize', 18); | 
|---|
| 138 | text(9.4,-8,'\fontsize{18}<') | 
|---|
| 139 |  | 
|---|
| 140 | figure | 
|---|
| 141 | hist(min(los(3,:),10),100);    AXIS([0 10 0 150]);   text(8,120,'\fontsize{18} SIDP') | 
|---|
| 142 | title('P_0 = 10, {}_{^{\fontsize{18}\theta}}^{\fontsize{16}\^}_{0} = 10');set( get(gcf, 'Children'), 'FontSize', 18);set(get(get(gcf, 'Children'),'Title'), 'FontSize', 18); | 
|---|
| 143 | text(9.4,-8,'\fontsize{18}<') | 
|---|
| 144 |  | 
|---|
| 145 | figure | 
|---|
| 146 | hist(min(los(4,:),10),100);    AXIS([0 10 0 150]);   text(8,120,'\fontsize{18} DP') | 
|---|
| 147 | title('P_0 = 10, {}_{^{\fontsize{18}\theta}}^{\fontsize{16}\^}_{0} = 10');set( get(gcf, 'Children'), 'FontSize', 18);set(get(get(gcf, 'Children'),'Title'), 'FontSize', 18); | 
|---|
| 148 | text(9.4,-8,'\fontsize{18}<') | 
|---|
| 149 |  | 
|---|
| 150 | end | 
|---|