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

%
%
% beta(z)=exp(S(z)-Sm)/[1/betam + 2*integ(K*exp(S(z)-Sm)dz)
%
% dove K è il LIDAR RATIO e
%
% S(z)=log(RCS)-2*integ((1-3K/8pi)*alfa_mol dz)
% lidar ratio,variabile a seconda della quota
quota=input('inserire quota per il lidar ratio');
% indice corrispondente alla quota indicata
val_1=trova_posiz(quota,quote_binnate);
% valore del lidar ratio
K1(1:val_1)=input('inserire lidar ratio');
quota=input('inserire quota per il lidar ratio');
val_2=trova_posiz(quota,quote_binnate);
val=val_2-val_1;
K2(1:val)=input('inserire lidar ratio');
quota=input('inserire quota per il lidar ratio');
val_3=trova_posiz(quota,quote_binnate);
val_n=val_3-val_2;
K3(1:val_n)=input('inserire lidar ratio');
%quota=input('inserire penultima quota per il lidar ratio');
% l'ultima quota che si considera è sempre pari a 7140 m
quota=7140;
val_4=trova_posiz(quota,quote_binnate);
val_m=val_4-val_3;
K4(1:val_m+1)=input('inserire lidar ratio');
% si crea a questo punto il vettore contenente tutti i valori
% del lidar ratio da inserire nell'algoritmo di klett
K=[K1 K2 K3 K4];
% beta molecolare alla quota di riferimento
betam=beta_molecolare(posiz);
alfa_molecolare=alfa_molecolare(1:end-8)';
alfa_corr=(1-3*K/(8*pi)).*alfa_molecolare; %(1x120)
% logaritmo dell'RCS
S=log(RCS_segnale_lidar)-2*
integrale(alfa_corr,dz,posiz);
%(1x120)
% logaritmo dell'RCS alla quota di riferimento
Sm=S(posiz);
%figure
%plot(quote_binnate,S)
%xlabel('quote [m]')
104