processamento di dati lidar per l'analisi dell'evoluzione ... - CO.RI.STA

%ylabel('logaritmo RCS')
esp=exp(S-Sm);
% figure
% plot(quote_binnate,esp)
esp_corr=esp.*K;
% coefficiente di backscattering totale
beta=esp./(1/betam+2*integrale(esp_corr,dz,posiz));
figure
plot(quote_binnate,beta)
xlabel('quote [m]')
ylabel('beta')
% Noto il beta totale è possibile infine calcolare la nostra
% incognita,ovvero il coefficiente di backscattering aerosolico,
% come differenza del beta appena calcolato e del beta
% molecolare
beta_aer=beta-beta_molecolare(1:end-8)';
figure
plot(quote_binnate,beta_aer)
xlabel('quote [m]')
ylabel('beta aerosolico')
% Per ottenere l'alfa aerosolico occorre moltiplicare il beta
% aerosolico appena calcolato per il lidar ratio
alfa_aer=beta_aer.*K;
figure
plot(quote_binnate,alfa_aer)
xlabel('quote [m]')
ylabel('alfa aerosolico')
%---------------------------------------------------------------
%-----------------Salvataggio del beta aerosolico---------------
name=lista_files(i,1:11);
file_name=[path_write,'\', name,'.mat'];
% salvataggio dati
save(file_name,'beta_aer');
%---------------------------------------------------------------
% Dopo aver effettuato la procedura per i 30 files analizzati,
% bisogna considerare i successivi 30 e ripetere il procedimento
% fino all'ultimo blocco di 30 files.
% Per passare ai successivi 30 è necessario incrementare il
105