Show
Ignore:
Timestamp:
09/14/10 13:38:21 (14 years ago)
Author:
zimamiro
Message:
 
Location:
applications/dual/SIDP/SIDPS/matlab
Files:
1 added
7 modified

Legend:

Unmodified
Added
Removed

  • applications/dual/SIDP/SIDPS/matlab/Beta.txt

    r1181 r1185  
    1   1.0000000e-004  7.3794257e-001  1.4757851e+000  2.2136277e+000  2.9514703e+000  3.6893128e+000  4.4271554e+000  5.1649980e+000  5.9028405e+000  6.6406831e+000 
     1  1.0000000e-004  7.0280615e-001  1.4055123e+000  2.1082184e+000  2.8109246e+000  3.5136307e+000  4.2163369e+000  4.9190430e+000  5.6217492e+000  6.3244553e+000 

  • applications/dual/SIDP/SIDPS/matlab/C.txt

    r1181 r1185  
    1   9.6000000e-001  4.1904762e-002  1.3714286e-001  3.2380952e-002  9.1428571e-002  6.6666667e-002  1.3333333e-002  1.3333333e-002  0.0000000e+000  7.6190476e-003 
     1  3.0533333e+000  1.5523810e+000  2.9333333e-001  1.9809524e-001  1.2000000e-001  1.0476190e-001  3.0476190e-002  7.8095238e-002  0.0000000e+000  2.4761905e-002 

  • applications/dual/SIDP/SIDPS/matlab/choose_best_control.m

    r1181 r1185  
    1 function best_control=choose_best_control(H,C,index,eta0,horizont,candidates,compare_parameters) 
     1function best_control=choose_best_control(H,C,index,apriori,horizont,candidates,compare_parameters) 
    22num_of_candidates=size(candidates,1); 
    33n_0=compare_parameters.n0; 
    … …  
    88    C(index)=candidates(i); 
    99    for j=1:n_0 
    10         realization1(i,j)=generate_realization(eta0,index,H,C,horizont); 
     10        realization1(i,j)=generate_realization(apriori,index,H,C,horizont); 
    1111    end 
    1212end 
    … …  
    4040                C(index)=candidates(i); 
    4141                for j=1:n_1 
    42                     realization2(j)=generate_realization(eta0,index,H,C,horizont); 
     42                    realization2(j)=generate_realization(apriori,index,H,C,horizont); 
    4343                end 
    4444                mean_values(i)=mean_values(i)/(1+n_1/n_0)+mean(realization2)/(1+n_0/n_1); 

  • applications/dual/SIDP/SIDPS/matlab/generate_realization.m

    r1181 r1185  
    1 function realization=generate_realization(eta,index,H,C,horizont) 
     1function realization=generate_realization(apriori,index,H,C,horizont) 
     2eta=apriori.eta0; 
     3switch(apriori.Hc) 
     4    case(0) 
     5        eta=0; 
     6    case(1) 
     7        eta=eta*rand; 
     8    case(2) 
     9        %eta=eta*abs((1+randn)/5^(index-1)); 
     10        %eta=eta/10^(-floor(-H(index)/H(end)*size(H,2))-1); 
     11        eta=eta/10^H(index); 
     12    case(3) 
     13        %eta=eta/2*abs((1+randn)/5^(index-1)); 
     14        eta=eta/10^H(index)*rand; 
     15end 
     16 
    217realization=0; 
    318beta=H(index); 

  • applications/dual/SIDP/SIDPS/matlab/interpol.m

    r1181 r1185  
    55else 
    66    if (pom(end)<0 || pom(end)==0) 
    7         rizeni=C(end); 
     7        rizeni=C(end)*2^pom(end); 
    88    else 
    99        for i=2:size(H,2) 

  • applications/dual/SIDP/SIDPS/matlab/nastav.m

    r1181 r1185  
    3232    n_grid_max=30; 
    3333  
    34     y0=1; 
    35     b0=1;   
     34    y0=100; 
     35    b0=0;   
    3636    P0=10; 
     37    Hc=0; 
    3738     
    3839    n0=20; 
    … …  
    6061apriori.b0=b0; 
    6162apriori.P0=P0; 
     63apriori.Hc=Hc; 
    6264end 

  • applications/dual/SIDP/SIDPS/matlab/sidp3.m

    r1181 r1185  
    88compare_parameters.rinott=4; 
    99 
    10 %apriori.eta0=abs(apriori.y0)/system.sigma; 
    11 apriori.eta0=0; 
     10apriori.eta0=abs(apriori.y0)/system.sigma; 
     11%apriori.eta0=0; 
    1212 
    1313apriori.beta0=abs(apriori.b0)/sqrt(apriori.P0); 
    … …  
    3333        for k=1:sidp_parameters.n_grid 
    3434            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); 
    3636        end 
    3737         
    … …  
    4949        end 
    5050         
    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(:))); 
    5353    end 
    5454     
    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);  
    5858%[w;x] 
    5959end