Coherent Backscattering from Multiple Scattering Systems - KOPS ...
Coherent Backscattering from Multiple Scattering Systems - KOPS ...
Coherent Backscattering from Multiple Scattering Systems - KOPS ...
You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
MATLAB codes<br />
end<br />
if n == 1<br />
p1 = 0 ;<br />
p = 1 ;<br />
t = mu ; % because t = n * mu * p - (n+1) * p1 = 1 * mu * 1 - 2 * 0<br />
else<br />
p0 = p1 ;<br />
p1 = p ;<br />
p = (2*n-1) / (n-1) * mu .* p1 - n / (n-1) * p0 ;<br />
t = n * mu .* p - (n+1) * p1 ;<br />
end<br />
S1 = S1 + (2*n+1) / (n*(n+1)) * (an*p + bn*t) ;<br />
S2 = S2 + (2*n+1) / (n*(n+1)) * (an*t + bn*p) ;<br />
S11 = 1/2 * ( S2 .* conj(S2) + S1 .* conj(S1) ) ;<br />
S33 = 1/2 * ( S1 .* conj(S2) + S2 .* conj(S1) ) ;<br />
I1 = S11 ;<br />
I2 = 1/2 * ( S11 + S33 ) ;<br />
% intensity after the scattering particle<br />
% after second transmission of circular polarizer<br />
plot (angle,I1,angle,I2)<br />
xlabel (’scattering angle [deg]’)<br />
ylabel (’intensity [a.u.]’)<br />
legend (’after scatterer’,’after circular polarizer’)<br />
end<br />
function J = J (n,x)<br />
J = sqrt(x*pi/2) * besselj(n+1/2,x) ;<br />
end<br />
function H = H (n,x)<br />
H = sqrt(x*pi/2) * ( besselj(n+1/2,x) + i * bessely(n+1/2,x) ) ;<br />
end<br />
function dJ = dJ (n,x)<br />
dJ = sqrt(pi/(2*x)) * ( (1+n) * besselj(n+1/2,x) - x * besselj(n+3/2,x) ) ;<br />
end<br />
function dH = dH (n,x)<br />
dH = sqrt(pi/(2*x)) * ( (1+n) * ( besselj(n+1/2,x) + i * bessely(n+1/2,x) ) ...<br />
- x * ( besselj(n+3/2,x) + i * bessely(n+3/2,x) ) ) ;<br />
end<br />
80