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

segnale_lidar_new=(segnale_lidar-fondo)*fs;
% è importante ricordare che i grafici devono essere in scala
% logaritmica, per cui una volta ottenuta la figura occorre
% modificare la scala delle y dal menù edit axis property
figure
plot(quote_binnate,segnale_lidar_new,'-g',quote_binnate,
molecolare,'-b')
xlabel('quote [m]')
%---------------------------------------------------------------
%-------------------RCS del segnale lidar-----------------------
% Una volta ottenuto il file mediato su 30 minuti devo
% moltipliacarlo per il quadrato della quota
RCS_segnale_lidar=segnale_lidar_new.*quote_binnate.^2;
% anche l'RCS deve essere in scala logaritmica
figure
plot(quote_binnate,RCS_segnale_lidar)
xlabel('quote [m]')
ylabel('Range corrected signal')
%---------------------------------------------------------------
%-------------------------Errore--------------------------------
% A questo punto occorre calcolare l'errore,visto come la
% deviazione standard dei files mediati ogni 30 minuti
errore=std(RCS_segnale_lidar);
%---------------------------------------------------------------
%----------------Calcolo del segnale molecolare-----------------
%------------------dovuto alle sole molecole--------------------
% Per poter calcolare il segnale molecolare dovuto alle sole
% molecole e non agli aerosol è necessario determinare dapprima
% i valori del number density alle varie quote.Questi valori ci
% vengono forniti richiamando il programma "ussa1976",che
% calcola inoltre i valori di pressione e temperatura alle varie
% quote.
number_density=import_density2('C:\Documents and Settings\
CORISTA\Desktop\luisa\tavole','densità.dat');
% sono 128 valori,si esculde quello relativo alla quota 0 m
quote_60_m=import_quote('C:\Documents and Settings\CORISTA\
Desktop\luisa\tavole da ussa1976','TABLE2_1.DAT');
% sono 128 valori,si esculde quello relativo alla quote 0 m
% Noto il number density è possibile calcolare il beta
% molecolare attraverso la seguente formula:
%
% beta_molecolare = Ng * dsigma/domega
101