Sdric Admin
Messages : 60 Date d'inscription : 12/10/2010
| Sujet: Re: Historique Mar 30 Nov - 16:20 | |
| Picasso du pool (version 2.01) cf. code matlab : - Code:
-
clear all close all
date=['07-Oct-2010';'26-Oct-2010';'29-Oct-2010';'08-nov-2010';'15-nov-2010';'22-nov-2010';'29-nov-2010';'02-dec-2010'];
pooler=['alex ';'cedric';'flo ';'ju ';'phil ';'rudo ';'steti ';'tibi '];
score=[ 0 0 0 0 0 0 0 0 ; 309 366 343 391 380 520 385 451 ; 368 408 425 479 484 603 477 511 ; 636 563 594 642 701 831 669 715 ; 822 763 849 804 852 969 821 847 ; 1131 944 993 952 980 1173 1131 1195; 1208 1167 1131 1089 1126 1282 1194 1315; 1274 1212 1226 1146 1167 1315 1266 1420; ]; fin_pool=datenum('10-apr-2011','dd-mmm-yyyy');
%%% score-average score_ave=score-mtimes(mean(score')',ones(1,8));
Couleurs=[ 0 0 0; %alex 0.5 0 0; %cedric 1 0 0; %flo 1 0.5 0; %ju 1 1 0; %Phil 0 1 0; %rudo 0.5 1 1; %steti 0 0 1;]; %tibi
%%% conversion date en num time=datenum(date,'dd-mmm-yyyy');
%%% prediction fin du pool S=size(score); score_pred=score(S(1),:)+(score(S(1),:)-score(S(1)-4,:))*(fin_pool-time(S(1)))/(time(S(1))-time(S(1)-4)); score_ave_pred=score_pred-mtimes(mean(score_pred')',ones(1,8)); for k=1:8 rnk=find(sort(score_pred,'descend')==score_pred(k)); if rnk==1 rank(k,:)=' 1er '; else rank(k,:)=[num2str(rnk),'ieme']; end end
%%% legende et echelle M=max(double(int64(min(min(score_ave))/50-1))*-50,double(int64(max(max(score_ave))/50+1))*50); %M2=max(double(int64(min(score_ave_pred)/50-1))*-50,double(int64(max(score_ave_pred)/50+1))*50); %M=max(M,M2); MINI=-M; MAXI=M; fig=figure('position',[10 10 1000 800])
hold on for k=1:8 h2=plot(0,0); set(h2,'color',Couleurs(k,:),'linewidth',10) end h=legend([pooler(1,:),' pred : ',rank(1,:),' (',num2str(int64(score_pred(1))),')'],... [pooler(2,:),' pred : ',rank(2,:),' (',num2str(int64(score_pred(2))),')'],... [pooler(3,:),' pred : ',rank(3,:),' (',num2str(int64(score_pred(3))),')'],... [pooler(4,:),' pred : ',rank(4,:),' (',num2str(int64(score_pred(4))),')'],... [pooler(5,:),' pred : ',rank(5,:),' (',num2str(int64(score_pred(5))),')'],... [pooler(6,:),' pred : ',rank(6,:),' (',num2str(int64(score_pred(6))),')'],... [pooler(7,:),' pred : ',rank(7,:),' (',num2str(int64(score_pred(7))),')'],... [pooler(8,:),' pred : ',rank(8,:),' (',num2str(int64(score_pred(8))),')']); set(h,'location','NorthWest','fontsize',10) time_line=[ '15-Oct-2010'; '01-Nov-2010'; '15-Nov-2010'; '30-Nov-2010'; '15-Dec-2010'; '30-Dec-2010'; '15-Jan-2011'; '30-Jan-2011'; '15-Feb-2011'; '30-Feb-2011'; '15-Mar-2011'; '30-Mar-2011'; '15-Apr-2011'; '30-Apr-2011';]; time_stick=datenum(time_line, 'dd-mmm-yyyy'); time_line=datestr(time_stick, 'dd/mm'); set(gca,'Xtick',time_stick,'XTickLabel',time_line,'fontsize',14) xlabel('Date','fontsize',20) ylabel('Écart à la moyenne','fontsize',20) axis([min(time) max(time) MINI MAXI]) %axis([min(time) fin_pool MINI MAXI])
%%%% graphics %%%
% les aires % S=size(score); for k=1:8 %h1=fill([time(S(1)) time(1) time'],[MINI MINI score_ave(:,k)'],'k'); h1=fill([fin_pool time(1) time' fin_pool],[MINI MINI score_ave(:,k)' score_ave_pred(k)],'k'); set(h1,'FaceColor',Couleurs(k,:),'facealpha',0.1,'edgeColor',Couleurs(k,:),'edgealpha',0.1) end % les lignes % S=length(time); for k=1:8 plot(time,score_ave(:,k),'color',Couleurs(k,:),'linewidth',2) plot([time(S) fin_pool],[score_ave(S,k),score_ave_pred(k)] ,'--k','color',Couleurs(k,:),'linewidth',2) end S=size(time_line); for k=[-7:-1 1:7] plot([min(time) fin_pool],k*[100 100],'-k') end for k=1:S(1) plot([time_stick(k) time_stick(k)],[MINI MAXI],'-k') end plot([min(time) fin_pool],[0 0],'--k','linewidth',3) annotation(fig,'rectangle',[0.845 0.0375 0.072 0.04],... 'FaceColor',[0.8 0.8 0.8],... 'Color',[0.8 0.8 0.8]);
time_line=datestr(time_stick, 'dd/mm'); MAXI=double(int64(max(max(score))/50+1))*50; %%%% COMPARAISON FIGURE 1 %%% mat_comp=[ 1 5; 1 6; 1 7; 1 8; 2 5; 2 6; 2 7; 2 8; 3 5; 3 6; 3 7; 3 8; 4 5; 4 6; 4 7; 4 8;]; pos=[ 0.03 0.75+0.03 0.24-0.03 0.25-0.03; 0.25+0.03 0.75+0.03 0.24-0.03 0.25-0.03; 0.50+0.03 0.75+0.03 0.24-0.03 0.25-0.03; 0.75+0.03 0.75+0.03 0.24-0.03 0.25-0.03; 0.03 0.50+0.03 0.24-0.03 0.25-0.03; 0.25+0.03 0.50+0.03 0.24-0.03 0.25-0.03; 0.50+0.03 0.50+0.03 0.24-0.03 0.25-0.03; 0.75+0.03 0.50+0.03 0.24-0.03 0.25-0.03; 0.03 0.25+0.03 0.24-0.03 0.25-0.03; 0.25+0.03 0.25+0.03 0.24-0.03 0.25-0.03; 0.50+0.03 0.25+0.03 0.24-0.03 0.25-0.03; 0.75+0.03 0.25+0.03 0.24-0.03 0.25-0.03; 0.03 0.00+0.03 0.24-0.03 0.25-0.03; 0.25+0.03 0.00+0.03 0.24-0.03 0.25-0.03; 0.50+0.03 0.00+0.03 0.24-0.03 0.25-0.03; 0.75+0.03 0.00+0.03 0.24-0.03 0.25-0.03;]; fig=figure('position',[10 10 1000 1000]); for f=1:16 axes('Parent',fig,'Position',pos(f,:)); hold on for k=mat_comp(f,:) plot(time,score(:,k),'color',Couleurs(k,:),'linewidth',2) end h1=fill([time' sort(time,'descend')'],[min(score'),sort(max(score'),'descend')],'k'); set(h1,'FaceColor',[0.1 0.1 0.1],'facealpha',0.1,'EdgeColor',[0.1 0.1 0.1],'Edgealpha',0.1) h=legend(pooler(mat_comp(f,1),:),pooler(mat_comp(f,2),:)); set(h,'location','NorthWest') set(gca,'Xtick',time_stick,'XTickLabel',time_line,'fontsize',8) axis([min(time) max(time) 0 MAXI])
S=size(time_line);
for k=1:S(1) plot([time_stick(k) time_stick(k)],[0 MAXI],'-k') end end
%%%% COMPARAISON FIGURE 2 %%% mat_comp=[ 8 5; 8 6; 8 7; 4 3; 7 5; 7 6; 3 2; 4 2; 6 5; 2 1; 3 1; 4 1;]; pos=[ 0.03 0.66+0.03 0.24-0.03 0.33-0.03; 0.25+0.03 0.66+0.03 0.24-0.03 0.33-0.03; 0.50+0.03 0.66+0.03 0.24-0.03 0.33-0.03; 0.75+0.03 0.66+0.03 0.24-0.03 0.33-0.03; 0.03 0.33+0.03 0.24-0.03 0.33-0.03; 0.25+0.03 0.33+0.03 0.24-0.03 0.33-0.03; 0.50+0.03 0.33+0.03 0.24-0.03 0.33-0.03; 0.75+0.03 0.33+0.03 0.24-0.03 0.33-0.03; 0.03 0.00+0.03 0.24-0.03 0.33-0.03; 0.25+0.03 0.00+0.03 0.24-0.03 0.33-0.03; 0.50+0.03 0.00+0.03 0.24-0.03 0.33-0.03; 0.75+0.03 0.00+0.03 0.24-0.03 0.33-0.03;]; fig=figure('position',[10 10 1000 750]); for f=1:12 %subplot(3,4,f) axes('Parent',fig,'Position',pos(f,:)); hold on for k=mat_comp(f,:) plot(time,score(:,k),'color',Couleurs(k,:),'linewidth',2) end h1=fill([time' sort(time,'descend')'],[min(score'),sort(max(score'),'descend')],'k'); set(h1,'FaceColor',[0.1 0.1 0.1],'facealpha',0.1,'EdgeColor',[0.1 0.1 0.1],'Edgealpha',0.1) h=legend(pooler(mat_comp(f,1),:),pooler(mat_comp(f,2),:)); set(h,'location','NorthWest') set(gca,'Xtick',time_stick,'XTickLabel',time_line,'fontsize',8) axis([min(time) max(time) 0 MAXI])
S=size(time_line);
for k=1:S(1) plot([time_stick(k) time_stick(k)],[0 MAXI],'-k') end end
| |
|