control_wind_error.m

计算大气后向散射的Rayleigh光谱宽度

light_speed = 3E8;
laser_frequency = 8.458E14;
lambda = light_speed/laser_frequency;
air_refraction_index = 1.0029;
plate_space = 12.5E-3;
nh = air_refraction_index*plate_space;
m = round(2*laser_frequency*air_refraction_index*plate_space/light_speed)
nh0 = m*light_speed/(2*laser_frequency)
delta_nh = 1E-9;
delta_frequency = (m*light_speed)/(2*(nh0+delta_nh))-laser_frequency
%%
A=0.01;
R = 0.64;
F = 4*R/(1-R)^2;
summit = (1-A/(1-R))^2;
trans1 = summit/(1+F*(sin(2*pi*(laser_frequency - 2.5E9)*(nh0+delta_nh)/light_speed))^2)
trans2 = summit/(1+F*(sin(2*pi*(laser_frequency + 2.5E9)*(nh0+delta_nh)/light_speed))^2)
RR = (trans1-trans2)/(trans1+trans2)

%%

%wind_speed = delta_frequency*light_speed/(2*laser_frequency)

%sin(2*laser_frequency*nh0*pi/light_speed)
%windspeedaxis = linspace(-50,50,101)';
%windspeed_divide = length(windspeedaxis);

%frequencyshiftaxis = 2*laser_frequency*windspeedaxis/light_speed;
% frequencyshift_divide = windspeed_divide;
%         tic
%         etalon1_trans = linspace(0,0,frequencyshift_divide)';
%         for k=1:frequencyshift_divide
%             etalon1_trans(k) = quad(@(x)double_edge_integrand...
%                 (x,frequencyshiftaxis(k),etalon1_position,etalon_fwhm),...
%                 laser_frequency-1E12,laser_frequency+1E12,1E-15);
%         k
%         end
%         etalon2_trans = flipud(etalon1_trans);
%         R_sensitivity_method1 = (etalon1_trans-etalon2_trans)./(etalon1_trans+etalon2_trans);
%         plot(R_sensitivity_method1,windspeedaxis);
%         toc
% save('sensitivity_matrix.dat','sensitivity_matrix','-double','-ascii');