|
Revision 1351, 1.2 kB
(checked in by zimamiro, 13 years ago)
|
|
|
| Line | |
|---|
| 1 | function sidps(regulator_parameters, compare_parameters,system,apriori)
|
|---|
| 2 | compare_parameters.student=tinv((1-compare_parameters.alpha)^(1/(regulator_parameters.num_of_candidates-1)),regulator_parameters.num_of_candidates-1);
|
|---|
| 3 | compare_parameters.rinott=4;
|
|---|
| 4 | regulator_parameters.dim=sum(regulator_parameters.H_tilde_on(:));
|
|---|
| 5 | u=zeros(2,1);
|
|---|
| 6 |
|
|---|
| 7 | [H_tilde C]=init_hyperstate(regulator_parameters,apriori); %vytvori H_tilde a mi, nastavi H0 a mi*
|
|---|
| 8 |
|
|---|
| 9 | for j=1:regulator_parameters.n_iter
|
|---|
| 10 | j
|
|---|
| 11 | regulator_parameters.search_region=regulator_parameters.gama^(j-1)*regulator_parameters.search_region_init;
|
|---|
| 12 |
|
|---|
| 13 | %prolez mrizku a iteruj rizeni
|
|---|
| 14 | for k=1:regulator_parameters.n_grid^regulator_parameters.dim
|
|---|
| 15 | H0=add_H_tilde0(regulator_parameters,H_tilde,k);
|
|---|
| 16 |
|
|---|
| 17 | u(:)=C(:,k);
|
|---|
| 18 | candidates=generate_candidates(u, regulator_parameters);
|
|---|
| 19 | C(:,k)=choose_best_control(H_tilde,C,H0,k,apriori,regulator_parameters,candidates,compare_parameters,system);
|
|---|
| 20 |
|
|---|
| 21 | end
|
|---|
| 22 |
|
|---|
| 23 | vypis=H_tilde;
|
|---|
| 24 | save 'H_tilde.txt' vypis -ASCII;
|
|---|
| 25 | save 'C.txt' C -ASCII;
|
|---|
| 26 | end
|
|---|
| 27 |
|
|---|
| 28 | end
|
|---|
| 29 |
|
|---|
