01.01.2015 Views

Rapport Bilan Carbone de l'université Paris Diderot.

Rapport Bilan Carbone de l'université Paris Diderot.

Rapport Bilan Carbone de l'université Paris Diderot.

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

Annexe 5.5 : Différents scripts du calcul <strong>de</strong> l’influence <strong>de</strong>s mo<strong>de</strong>s <strong>de</strong><br />

transport ainsi que le calcul <strong>de</strong>s distances (MatLab)<br />

Calcul <strong>de</strong>s distances totales parcourues par<br />

mo<strong>de</strong> <strong>de</strong> transports et évolution <strong>de</strong> l’erreur<br />

associée<br />

clear all<br />

close all<br />

[dv dc dd dp db dt dm dg dverr dcerr d<strong>de</strong>rr<br />

dperr dberr dterr dmerr dgerr ndt dtot]...<br />

= routine_lecture('donnees_up<strong>de</strong>tu.xlsx');<br />

%profil=[dv,dc,dd,dp,db,dt,dm,dg]';<br />

profilerr=[dverr,dcerr,d<strong>de</strong>rr,dperr,dberr,dterr<br />

,dmerr,dgerr];<br />

dtot=dtot.*ndt;<br />

M=[dv,dc,dd,dp,db,dt,dm,dg]';<br />

dtot2=sum(M);<br />

dmi=1; %valeur minimum <strong>de</strong> trajet<br />

dma=850; %valeur maximum <strong>de</strong> trajet<br />

ik=find(dtot2>=dmi&dtot2=2<br />

for col=1:8<br />

tmp=Resum(1:k,col);<br />

ResumTmpTot(k-1,col)=mean(tmp,1);<br />

ResumerrTmpTot(k-<br />

1,col)=std(tmp,1);<br />

end<br />

end<br />

end<br />

comp<br />

txt=strcat('DistParcourueEtudiants',num2str(k)<br />

, 'Universités.xls')<br />

txt2=strcat('ErrTotParcourueEtudiants',num2str<br />

(k), 'Universités.xls')<br />

xlswrite(txt,Resum) % matrice <strong>de</strong> données<br />

moyennées pour toutes les Universités<br />

xlswrite(txt2,Resumerrco) % matrice <strong>de</strong>s<br />

erreurs moyennées pour toutes les Universités<br />

n=size(D);<br />

b=n(1,2);<br />

%concaténation <strong>de</strong>s données<br />

for j=1:b<br />

ResumeTot(1,j)=mean(Resum(:,j)); %<br />

Distance totale moyennées parcourue par mo<strong>de</strong><br />

<strong>de</strong> transport<br />

ResumeTot(2,j)=std(Resum(:,j)); % Erreur<br />

absolue sur les données<br />

ResumeTot(3,j)=mean(Resumerrco(:,j)); %<br />

Erreur sur les constantes<br />

end<br />

'Voici le résumé <strong>de</strong>s calculs'<br />

ResumeTot<br />

txt3=strcat('DistTotParcourueEtudiants',num2st<br />

r(k), 'Universités.xls')<br />

xlswrite(txt3,ResumeTot)<br />

TotalKm=sum(ResumeTot(1,:))<br />

% figure<br />

% explo<strong>de</strong>=[0 0 0 0 0 0 1 0];<br />

%<br />

pie(ResumeTot(1,:),explo<strong>de</strong>,{'Voi','Covoit','2<br />

Roues','VoP','Bus','TER','MRT','TGV'})<br />

%% Répartition <strong>de</strong>s distributions <strong>de</strong>s distances<br />

totales calculées par mo<strong>de</strong><br />

%% <strong>de</strong> transport<br />

col =<br />

['Voi';'Cov';'2Ro';'VoP';'Bus';'TER';'MRT';'TG<br />

V'];<br />

n=size(Resum);<br />

b=n(1,2);<br />

figure<br />

for i=1:b;<br />

subplot(4,2,i)<br />

hist(Resum(:,i));<br />

title(col(i,:))<br />

xlabel('Distance(km)','Fontsize',8)<br />

ylabel('Nombre Universités','Fontsize',8)<br />

end<br />

txt=strcat('Répartition<strong>de</strong>sdistmo<strong>de</strong><strong>de</strong>transports<br />

',num2str(k),'Universités.png');<br />

saveas(gcf,txt)<br />

%% Evolution <strong>de</strong>s erreurs liées au calcul en<br />

fonction <strong>de</strong>s mo<strong>de</strong>s <strong>de</strong> transport<br />

Resulttmp=ResumerrTmpTot./ResumTmpTot; %<br />

erreur relative<br />

Result=ResumTmpTot+ResumTmpTot.*Resulttmp;<br />

Result2=ResumTmpTot-ResumTmpTot.*(Resulttmp);<br />

140 Université <strong>Paris</strong> Di<strong>de</strong>rot – <strong>Paris</strong> 7 <strong>Paris</strong> Di<strong>de</strong>rot Développement Durable

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!