root/applications/dual/SIDP/SIDPS/matlab/vyst.m @ 1181

Revision 1181, 1.4 kB (checked in by zimamiro, 14 years ago)
Line 
1function los=vyst(system,apriori,N)
2los=zeros(3,3);
3a=2;
4
5%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
6apriori.P0=1;
7apriori.b0=0;
8los(:,1)=mc_study(system,apriori,N);
9
10apriori.b0=1;
11los(:,2)=mc_study(system,apriori,N);
12
13apriori.b0=10;
14los(:,3)=mc_study(system,apriori,N);
15los=los./(ones(3,1)*los(3,:));
16pom(:,:)=los(1:2,:);
17
18figure
19subplot(1,2,1)
20bar(pom,'grouped'); hold on
21colormap summer;
22set(gca,'ylim',[0 a]);
23plot(0:3,0*(1:4)+1,'r--');
24
25title('\fontsize{18} P_0 = 1')
26set(gca,'xtick',[]);
27text(0.9,-0.05,'\fontsize{18} CC')
28text(1.8,-0.05,'\fontsize{18} SIDPS')
29set( get(gcf, 'Children'), 'FontSize', 18);
30
31%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
32apriori.P0=10;
33apriori.b0=0;
34los(:,1)=mc_study(system,apriori,N);
35
36apriori.b0=1;
37los(:,2)=mc_study(system,apriori,N);
38
39apriori.b0=10;
40los(:,3)=mc_study(system,apriori,N);
41
42los=los./(ones(3,1)*los(3,:));
43pom(:,:)=los(1:2,:);
44
45subplot(1,2,2)
46bar(pom,'grouped'); hold on
47colormap summer ;
48set(gca,'ylim',[0 a]);
49
50title('\fontsize{18} P_0 = 10')
51legend('{}_{^{\fontsize{18}\theta}}^{\fontsize{16}\^}_{0} = 0','{}_{^{\fontsize{18}\theta}}^{\fontsize{16}\^}_{0} = 1','{}_{^{\fontsize{18}\theta}}^{\fontsize{16}\^}_{0} = 10')
52set(gca,'xtick',[]);
53plot(0:3,0*(1:4)+1,'r--');
54
55text(0.9,-0.05,'\fontsize{18} CC')
56text(1.8,-0.05,'\fontsize{18} SIDPS')
57set( get(gcf, 'Children'), 'FontSize', 18);
58
59
60end
Note: See TracBrowser for help on using the browser.