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 />
function R = reflectivity (n1,n2)<br />
if n1 > n2<br />
theta_totalreflection = asin (n2/n1) ;<br />
else<br />
theta_totalreflection = pi/2 ;<br />
end<br />
theta1 = 0 : pi/100 : theta_totalreflection ;<br />
theta2 = theta_totalreflection+pi/100 : pi/100 : pi/2 ;<br />
theta = [ theta1, theta2 ] ;<br />
R_parallel = ( n1 * cos(theta1) ...<br />
- n2 * sqrt( 1 - (n1 / n2 * sin(theta1)).^2 ) ).^2 ./ ...<br />
( n1 * cos(theta1) ...<br />
+ n2 * sqrt( 1 - (n1 / n2 * sin(theta1)).^2 ) ).^2 ;<br />
R_parallel = [ R_parallel , ones(size(theta2)) ] ;<br />
R_perp = ( n1 * sqrt( 1 - (n1 / n2 * sin(theta1)).^2 ) ...<br />
- n2 * cos(theta1) ).^2 ./ ...<br />
( n1 * sqrt( 1 - (n1 / n2 * sin(theta1)).^2 ) ...<br />
+ n2 * cos(theta1) ).^2 ;<br />
R_perp = [ R_perp , ones(size(theta2)) ] ;<br />
C1 = trapz ( theta , ...<br />
(R_parallel + R_perp) / 2 .* sin(theta) .* cos(theta) ) ;<br />
C2 = trapz ( -theta , ...<br />
(R_parallel + R_perp) / 2 .* sin(-theta) .* (cos(-theta)).^2 ) ;<br />
R = ( 3 * C2 + 2 * C1 ) / ( 3 * C2 - 2 * C1 + 2 ) ;<br />
end<br />
86