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 />

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

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

Saved successfully!

Ooh no, something went wrong!