MATLAB Functions for Mie Scattering and Absorption
MATLAB Functions for Mie Scattering and Absorption
MATLAB Functions for Mie Scattering and Absorption
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
15<br />
<strong>Mie</strong> angular scattering: m=5+0.4i, x=1<br />
90<br />
0.4<br />
120<br />
60<br />
0.3<br />
150<br />
0.2<br />
30<br />
0.1<br />
180 0<br />
210<br />
330<br />
Figure 2: Angular <strong>Mie</strong>-scattering diagram<br />
of S<br />
2<br />
1<br />
(upper Half circle) <strong>and</strong> of<br />
2<br />
2<br />
S (lower half circle). Here, scattering<br />
in the backward hemisphere is<br />
slightly larger than in the <strong>for</strong>ward<br />
hemisphere ( >1, a large number of spherical harmonics have to be computed. Here we<br />
consider the example of m=2+0.01i. Figure 4 shows the <strong>Mie</strong> Efficiencies over the x<br />
range from 0 to 25. As a result of the Extinction Paradox (van de Hulst, 1957), Q ext<br />
approaches the value 2 <strong>for</strong> very large x. The computed values are slightly higher in<br />
Figure 4. For increasing x, Qext decreases to 2.09 at x=100, to 2.01 at x=2800 <strong>and</strong><br />
to Qext=2.0014 <strong>for</strong> x=40’000. The program still works at x=70'000 (Qext =2.0012),<br />
but <strong>for</strong> x>80'000 NaN values are returned.