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.

Evaluation of the wide angle data<br />

% remove NANs and data with huge errors <strong>from</strong> damaged photodiodes<br />

m = 1 ;<br />

for n = 1 : 256<br />

if isnan (cone(n))<br />

elseif cone_delta(n) > 0.1*P_samp<br />

else<br />

Cone(m) = cone(n) ;<br />

Delta(m) = cone_delta(n) ;<br />

Theta(m) = theta(n) ;<br />

m = m + 1 ;<br />

end<br />

end<br />

cone = Cone ;<br />

cone_delta = Delta ;<br />

theta = Theta ;<br />

clear Theta Cone Delta<br />

% fold data (i.e. -pi/2 < theta < pi/2 -> 0 < theta < pi/2)<br />

m = 1 ;<br />

n = 1 ;<br />

while n < size(theta,2)<br />

if theta(n+1) > theta(n)<br />

Theta(m) = theta(n) ;<br />

Theta(m+1) = theta(n+1) ;<br />

Cone(m) = cone(n) ;<br />

Cone(m+1) = cone(n+1) ;<br />

Delta(m) = cone_delta(n) ;<br />

Delta(m+1) = cone_delta(n+1) ;<br />

n = n + 1 ;<br />

else % theta(n) == theta(n+1)<br />

Theta(m) = theta(n) ;<br />

Cone(m) = (cone(n) + cone(n+1)) / 2 ;<br />

Delta(m) = 0.5 * sqrt ( (cone_delta(n)).^2 + (cone_delta(n+1)).^2 ) ;<br />

n = n + 2 ;<br />

end<br />

m = m + 1 ;<br />

end<br />

cone = Cone ;<br />

cone_delta = Delta ;<br />

theta = Theta ;<br />

%% Derive experimental and theoretical cones<br />

% Correct cone with albedo mismatch<br />

corr_cone = aleph * cone ;<br />

corr_cone_delta = aleph * cone_delta ;<br />

83

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

Saved successfully!

Ooh no, something went wrong!