Changeset 1185 for applications/dual/SIDP/SIDPS/matlab/sidp3.m
- Timestamp:
- 09/14/10 13:38:21 (14 years ago)
- Files:
-
- 1 modified
Legend:
- Unmodified
- Added
- Removed
-
applications/dual/SIDP/SIDPS/matlab/sidp3.m
r1181 r1185 8 8 compare_parameters.rinott=4; 9 9 10 %apriori.eta0=abs(apriori.y0)/system.sigma;11 apriori.eta0=0;10 apriori.eta0=abs(apriori.y0)/system.sigma; 11 %apriori.eta0=0; 12 12 13 13 apriori.beta0=abs(apriori.b0)/sqrt(apriori.P0); … … 33 33 for k=1:sidp_parameters.n_grid 34 34 candidates=generate_candidates(C(k), search_region, sidp_parameters.num_of_candidates); 35 C(k)=choose_best_control(H,C,k,apriori .eta0,sidp_parameters.horizont,candidates,compare_parameters);35 C(k)=choose_best_control(H,C,k,apriori,sidp_parameters.horizont,candidates,compare_parameters); 36 36 end 37 37 … … 49 49 end 50 50 51 %subplot(2,1,1);52 %plot(H,C,['',col(j)]); hold on; a=max(a,max(C(:)));51 % subplot(2,1,1); 52 % plot(H,C,['',col(j)]); hold on; a=max(a,max(C(:))); 53 53 end 54 54 55 %xlabel('\fontsize{18} \beta'); ylabel('\fontsize{18} \nu^{ (2)}_0'); set(gca,'ylim',[0 a]); set(gca,'xlim',[0 H(end)]); legend('\pi_1','\pi_2','\pi_3','\pi_4');56 %subplot(2,1,2);57 %plot(0:sidp_parameters.n_iter,[w' x' ones(size(x,2),1)]); xlabel('\fontsize{18} iterace'); ylabel('\fontsize{18} relativn�tr�'); legend('mean','median');set(gca,'ylim',[0.99 max(w(1),x(1))]); set(gca,'xtick',0:sidp_parameters.n_iter); set( get(gcf, 'Children'), 'FontSize', 18);55 %xlabel('\fontsize{18} \beta'); ylabel('\fontsize{18} \nu^{ (2)}_0'); set(gca,'ylim',[0 a]); set(gca,'xlim',[0 H(end)]); legend('\pi_1','\pi_2','\pi_3','\pi_4'); 56 %subplot(2,1,2); 57 %plot(0:sidp_parameters.n_iter,[w' x' ones(size(x,2),1)]); xlabel('\fontsize{18} iterace'); ylabel('\fontsize{18} relativn�tr�'); legend('mean','median');set(gca,'ylim',[0.99 max(w(1),x(1))]); set(gca,'xtick',0:sidp_parameters.n_iter); set( get(gcf, 'Children'), 'FontSize', 18); 58 58 %[w;x] 59 59 end