mie散射物理参量的一种改进数值算法.docxVIP

mie散射物理参量的一种改进数值算法 1 考虑复折射率对计算结果的影响主要表现 首先，可以严格地分解mie和debie运动，但mie运动的计算并不复杂，因为mie分散计算非常复杂。另一方面，直接编程mie理论可能会导致严重的计算错误，尤其是在计算复变量bessel函数时。到1968年Dave最先发表了完整的算法,以后Lentz和Wiscombe又针对部分计算提出了改进算法,但Dave的算法不能应用到大尺度参数的情况下,Wiscombe的方案详细讨论了计算时间,循环次数,得到了准确的结果,但是他没有考虑复折射率较大的粒子,安庆师范学院的张杰等人对大气中具有复折射率气溶胶粒子的散射光学特性,余其铮等人利用连分式进行计算,通过对计算次数进行拟合,提高了计算的速度,但是对于尺度参数在100以外的拟合公式没有讨论。王式民等利用后向递推办法大大提高了计算速度,但是文中提到的递推表达式是和粒子的复折射率相关的表达式。朱震等人在以往Mie散射计算的基础上提出一种新的算法,此算虽然区别于传统的递推算法,精度高,但仍然使用了递推算法,不可避免的会出现积累误差。前面所说的改进都是在同一种算法上改进,而姬丰等人利用Matlab对散射参量进行数值模拟,得到了正确的结果,但是它没有考虑折射率虚部较大和粒子尺度参数较大对计算结果的影响。而我们的数值计算完全不依赖连分式或者后向递推算法,直接利用Matlab中内