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

%
% dove Ng rappresenta il numero di molecole di gas x unità di
% volume e dsigma/domega la sezione d'urto
% esplicitando:
%
% beta_molecolare=number_density*5.45*(550/lambda)^4*10^-28
lambda=input('inserisci lunghezza d_onda:');
% Coefficiente di backscattering
% il number density ha le dimensioni di 1/m^3,mentre la sezione
% d'urto di cm^2/sr,per cui,passando tutto in metri occore
% moltiplicare per 10^-4
beta_molecolare=(5.45*(550/lambda)^4)*10^-32*number_density;
% coefficiente di estinzione dell'atmosfera
alfa_molecolare=beta_molecolare*8*pi/3;
% incremento relativo alla quota
dz=60;
spessore_ottico=sum(alfa_molecolare)*dz;
% E' possibile a questo punto calcolare il segnale molecolare
% applicando la seguente formula:
%
% segnale molecolare =
% = number_density * e^(-2*spessore ottico)/quota^2
segn_molecolare=number_density*exp(-2*spessore_ottico)./
quote_60_m.^2;
segn_molecolare=segn_molecolare';
% Una volta ottenuto il molecolare devo moltipliacarlo per il
% quadrato della quota,in modo da ottenere l'RCS del molecolare
RCS_molecolare=segn_molecolare(1:end-8).*quote_binnate.^2;
% (1x120)
figure
plot(quote_binnate,RCS_molecolare)
xlabel('quote (m)')
ylabel('RCS del molecolare')
%---------------------------------------------------------------
%---------------------------------------------------------------
%--------------------Algoritmo di Klett-------------------------
%---------------------------------------------------------------
% Occorre determinare,a questo punto,il valore della quota a
% partire dal quale il segnale laser e il molecolare coincidono.
quota_rif=input('inserire quota di riferimento');
102