1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36
|
clc
clear
M='02' % mois de fevrier ;
cd(strcat('C:\Users\majda\Desktop\matrice_raport\Rap_lidar\',M,'\'))
x=1;
for J=[04 06 08 07] % les jours qui existe dans le mois de fevrier
if J<10
j=strcat('0',num2str(J));
else
j=num2str(J);
end
load(strcat('Rap_',M,'_',j,'.mat'))
li1200(1:2522,x)=Rap_li1200;
lio3s(1:1262,x)=Rap_lio3s;
x=x+1;
end
clear J M Rap_li1200 Rap_lio3s j x
for g=1:2522
moy_Map(g,1)=nanmean(Map(g,:));
moy_Arl(g,1)=nanmean(Arl(g,:));
end
clear g Map Arl
figure(1)
clf(1)
plot(moy_Map(856:2522),alt(856:2522),'r')
hold on
plot(moy_Arl(856:2522),alt(856:2522),'b') |
Partager