Coherent Backscattering from Multiple Scattering Systems - KOPS ...
Coherent Backscattering from Multiple Scattering Systems - KOPS ...
Coherent Backscattering from Multiple Scattering Systems - KOPS ...
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
MATLAB codes<br />
% Derive experimental diffuson and cooperon<br />
exp_diffuson = diffuson (theta,l,R) * P_samp ;<br />
exp_diffuson_delta = diffuson (theta,l,R) * Perror_samp ;<br />
exp_cooperon = (corr_cone - exp_diffuson) / P_samp ;<br />
exp_cooperon_delta = ...<br />
sqrt ( (corr_cone_delta / P_samp).^2 ...<br />
+ (exp_diffuson_delta / P_samp).^2 ...<br />
+ ((corr_cone - exp_diffuson) / P_samp^2 * Perror_samp).^2 ) ;<br />
% Calculate theory curve<br />
theory_theta = 0 : 0.001*pi : pi/2 ;<br />
theory_cooperon = enhancement * cooperon (theory_theta,k,l,R,D,tau) ;<br />
%% Calculate energy<br />
e_theta = [ 0 , theta , pi/2 ] ;<br />
e_cooperon = [ 1 , exp_cooperon , 0 ] ;<br />
e_cooperon_delta = [ 0 , exp_cooperon_delta , 0 ] ;<br />
E = 2*pi ...<br />
* sum ( ( e_theta(:,2:end) - e_theta(:,1:end-1) ) ...<br />
* 1/2 .* ( e_cooperon(:,1:end-1) .* sin(e_theta(:,1:end-1)) ...<br />
+ e_cooperon(:,2:end) .* sin(e_theta(:,2:end)) ) ) ;<br />
E_delta = 2*pi ...<br />
* sqrt ( sum ( ( 1/2 * ( e_theta(1,2:end) - e_theta(1,1:end-1) ) ...<br />
.* e_cooperon_delta(1:end-1) ).^2 ...<br />
.* sin (1/2 * ( e_theta(1,2:end) + e_theta(1,1:end-1) )) ...<br />
+ ( 1/2 * ( e_theta(1,2:end) - e_theta(1,1:end-1) ) ...<br />
.* e_cooperon_delta(2:end) ).^2 ...<br />
.* sin (1/2 * ( e_theta(1,2:end) + e_theta(1,1:end-1) )) ) ) ;<br />
disp ([’E = ’,num2str(E),’+-’,num2str(E_delta)])<br />
%% Plot data<br />
angles = theta / pi * 180 ;<br />
theory_angles = theory_theta / pi * 180 ;<br />
plot ( angles, exp_cooperon, ’b’, ...<br />
theory_angles, theory_cooperon, ’r’, ...<br />
angles, exp_cooperon+exp_cooperon_delta, ’:b’, ...<br />
angles, exp_cooperon-exp_cooperon_delta, ’:b’ )<br />
xlabel (’scattering angle [deg]’)<br />
ylabel (’cooperon’)<br />
legend (’backscattering data’,’theory’)<br />
end<br />
84