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.
Evaluation of the wide angle data<br />
function d = diffuson (theta,l,R)<br />
z0 = 2/3 * l * (1 + R) / (1 - R) ;<br />
mu = cos(theta) ;<br />
d = ( mu .* (z0/l + mu./(mu+1)) ) ./ (z0/l + 1/2) ;<br />
end<br />
function c = cooperon (theta,k,l,R,D,tau)<br />
z0 = 2/3 * l * (1 + R) / (1 - R) ;<br />
mu = cos(theta) ;<br />
q = k * sin(theta) ;<br />
L_a = sqrt (D * tau) ;<br />
q_a = sqrt (q.^2 + 1/(D*tau)) ;<br />
c = ( l * (l + L_a)^2 * ( (1 - exp(-2*q_a*z0)) ./ ...<br />
(q_a*l) + 2*mu ./ (mu+1) ) ) ...<br />
./ ( L_a^2 * ( l + (1 - exp(-2*z0/L_a)) * L_a ) * ...<br />
( l*q_a + (mu+1) ./ (2*mu) ).^2 ) ;<br />
c_add = (8 * pi * a(k,l,R)) / (3 * k^2 * l * (l + 2 * z0)) * mu ./ (mu + 1) ;<br />
c = c - c_add ;<br />
end<br />
function a = a (k,l,R)<br />
theta = 0.00001*pi : 0.00001*pi : pi/2 ;<br />
z0 = 2/3 * l * (1 + R) / (1 - R) ;<br />
mu = cos(theta) ;<br />
q = k * abs(sin(theta)) ;<br />
E_A = trapz ( theta, ...<br />
3/(8*pi) * 1./(q*l + (mu+1)./(2*mu)).^2 ...<br />
.* ( (1 - exp(-2*q*z0))./(q*l) + (2*mu)./(mu+1) ) ...<br />
.* sin(theta) ) ;<br />
E_BC = trapz ( theta, ...<br />
1/(k*l)^2 * mu./(mu+1) ...<br />
.* sin(theta) ) ;<br />
a = E_A / E_BC ;<br />
end<br />
85