Changeset 1185 for applications/dual

Timestamp:

09/14/10 13:38:21 (14 years ago)

Author:

zimamiro

Message:

 

Location:

applications/dual/SIDP/SIDPS

Files:

• matlab/Beta.txt (modified) (1 diff)
• matlab/C.txt (modified) (1 diff)
• matlab/choose_best_control.m (modified) (3 diffs)
• matlab/generate_realization.m (modified) (1 diff)
• matlab/interpol.m (modified) (1 diff)
• matlab/nastav.m (modified) (2 diffs)
• matlab/rid_y0.m (added)
• matlab/sidp3.m (modified) (3 diffs)
• text/SIDPS.pdf (modified) (previous)

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

@@ -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.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

@@ -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
+  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

@@ -1 @@
-function best_control=choose_best_control(H,C,index,eta0,horizont,candidates,compare_parameters)
+function best_control=choose_best_control(H,C,index,apriori,horizont,candidates,compare_parameters)
 num_of_candidates=size(candidates,1);
 n_0=compare_parameters.n0;
@@ -8 +8 @@
     C(index)=candidates(i);
     for j=1:n_0
-        realization1(i,j)=generate_realization(eta0,index,H,C,horizont);
+        realization1(i,j)=generate_realization(apriori,index,H,C,horizont);
     end
 end
@@ -40 +40 @@
                 C(index)=candidates(i);
                 for j=1:n_1
-                    realization2(j)=generate_realization(eta0,index,H,C,horizont);
+                    realization2(j)=generate_realization(apriori,index,H,C,horizont);
                 end
                 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

@@ -1 @@
-function realization=generate_realization(eta,index,H,C,horizont)
+function realization=generate_realization(apriori,index,H,C,horizont)
+eta=apriori.eta0;
+switch(apriori.Hc)
+    case(0)
+        eta=0;
+    case(1)
+        eta=eta*rand;
+    case(2)
+        %eta=eta*abs((1+randn)/5^(index-1));
+        %eta=eta/10^(-floor(-H(index)/H(end)*size(H,2))-1);
+        eta=eta/10^H(index);
+    case(3)
+        %eta=eta/2*abs((1+randn)/5^(index-1));
+        eta=eta/10^H(index)*rand;
+end
+
 realization=0;
 beta=H(index);

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

@@ -5 +5 @@
 else
     if (pom(end)<0 || pom(end)==0)
-        rizeni=C(end);
+        rizeni=C(end)*2^pom(end);
     else
         for i=2:size(H,2)

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

@@ -32 +32 @@
     n_grid_max=30;
  
-    y0=1;
-    b0=1;  
+    y0=100;
+    b0=0;  
     P0=10;
+    Hc=0;
     
     n0=20;
@@ -60 +61 @@
 apriori.b0=b0;
 apriori.P0=P0;
+apriori.Hc=Hc;
 end

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

@@ -8 +8 @@
 compare_parameters.rinott=4;
 
-%apriori.eta0=abs(apriori.y0)/system.sigma;
-apriori.eta0=0;
+apriori.eta0=abs(apriori.y0)/system.sigma;
+%apriori.eta0=0;
 
 apriori.beta0=abs(apriori.b0)/sqrt(apriori.P0);
@@ -33 +33 @@
         for k=1:sidp_parameters.n_grid
             candidates=generate_candidates(C(k), search_region, sidp_parameters.num_of_candidates);
-            C(k)=choose_best_control(H,C,k,apriori.eta0,sidp_parameters.horizont,candidates,compare_parameters);
+            C(k)=choose_best_control(H,C,k,apriori,sidp_parameters.horizont,candidates,compare_parameters);
         end
         
@@ -49 +49 @@
         end
         
-        %subplot(2,1,1);
-        %plot(H,C,['',col(j)]); hold on; a=max(a,max(C(:)));
+    %    subplot(2,1,1);
+    %    plot(H,C,['',col(j)]); hold on; a=max(a,max(C(:)));
     end
     
-   % 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');
-   % subplot(2,1,2);
-   % 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); 
+    %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');
+    %subplot(2,1,2);
+    %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); 
 %[w;x]
 end