26.12.2013 Views

Coherent Backscattering from Multiple Scattering Systems - KOPS ...

Coherent Backscattering from Multiple Scattering Systems - KOPS ...

Coherent Backscattering from Multiple Scattering Systems - KOPS ...

SHOW MORE
SHOW LESS

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

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!