Coherent Backscattering from Multiple Scattering Systems - KOPS ...

MATLAB codes
Angular intensity distribution of single scattering
The following MATLAB script calculates the angular distribution of the singly scattered intensity
after the sequence ‘circular polarizer – spherical scatterer – circular polarizer’ in both
backscattering setups, as derived in secs. 2.2 and 5.2.3.
function mie_setup
refractive_particle = 1.52 ;
refractive_surround = 1.33 ;
wavelength = 575e-9 ;
diameter = 150e-6 ;
% refractive index of the scatterer
% refractive index of the surrounding medium
% wavelength of the scattered light
% particle diameter
a = diameter / 2 ;
m = refractive_particle / refractive_surround ;
k = 2 * pi * refractive_surround / wavelength ;
x = k * a ;
angle = 180 : 0.01 : 181 ;
theta = angle / 180 * pi ;
mu = cos(theta) ;
% or larger angular range for wide angle setup
S1 = 0 ;
S2 = 0 ;
N = 50 + 1.03 * x ;
N = fix (N) ;
% avoids numerical overflow in Bessel functions
for n = 1 : N
A = J(n,m*x) ;
B = J(n,x) ;
C = dJ(n,x) ;
D = dJ(n,m*x) ;
E = H(n,x) ;
F = dH(n,x) ;
an = ( m * A * C - B * D ) / ( m * A * F - E * D ) ;
bn = ( A * C - m * B * D ) / ( A * F - m * E * D ) ;