基于Matlab的大气气溶胶粒子散射特性研究.pdfVIP

高 新 技 术 SCIENCE＆TE0HNOLo。Y 基于Matlab的大气气溶胶粒子散射特性研究 贺雅奇 (1．西安电子科技大学研究生院 西安 710071；2．陕西理工学院物理系 陕西汉中 723001) r 摘 要：本文主要利用Matlab壕程软件模拟大气气溶胶柱子的模型蛄构和电磁散射的传播过程。通过数值计算得到了气溶胶粒子散射 特性的一般规律。 关键词：气溶胶柱子 Matlab 电磁散射 中图分类号：O4 文献标识码：A 文章编号：1672-379t(2oo8)1l(b)一0009—0I 大气气溶胶是大气辐射遥感过程 中十分重要的修正量 ，也是 基于散射原理的能见度器测理论中的主要影响因子 ，因此研究气 q(．，，，e)=i7 I【(，m，e)l：一l(N~m,0)】 溶胶粒子的散射特性具有重要意义。长期 以来 ，对气溶胶粒子散 射特性的理论研究作为一种实际问题的理论抽象 ，该模型的散射 ： 业盟 (7) 特性成为许 多相关 问题研究的核心和基础 。 ∑(2n+1)【()f+Ib,，(，m)l：】 1气溶胶粒子模型及散射特性的Matlab实现 1．1气溶胶粒子的结构模型 大气气溶胶是由大气介质和混合于其中的固体或液体颗粒物 舢)=喜 c0s。( c。so)】8() 组成的体系。大气气溶胶是 由固体粒子、液体粒子或固体及液体 混合粒子悬浮于气体介质中的一种悬胶体 ，其降落速度极小。溶 e)=喜 os。)州 c。se)】9() 胶通常也指悬浮于大气中的微小粒子 ，粒径范围可以从0．001pm到 几十微米 。利用Mandelbrot曾提出的整体与局部具有 自相似特征 其 (c。s9)=d p se)；7c(cosO)= 1 (c。so)(1。) 的集合分形理论 ，将大气气溶胶粒子的分形维数进行确定 。可以 把粒子体系的数 目N与包围溶胶粒子体系的最小球半径 R之间的 对于 多粒子 体系有 关系表述成更普遍的形式 ： ( 帅 川 …) N：p(／心)。， (1) 式 中D，就是所涉及到的烟尘 溶胶粒子的分形维数 ，R也称 为 回旋半径R ，其计算公式为 ： g(Lm．O)= ，r，帆e r，州 (12) = 2 (11)式的物理意义是，橐一方向的散射辐射能量和各向同性散 射时该方向的散射辐射能量的比值 ；(12)式的意义是气溶胶粒子在 式中的r是单体 f的球心到体系中心 的距离。 茸十值 糯 时 的 角 备 布 函 特 fill 值 #；相 函 特 单就气溶胶散射特性 的数值模拟而言 ，气溶胶粒子的复杂多 变的结构对其散射特性影响更大 ，故对其空间结构的模拟显得更 2Matlab编程计算结果 为重要。 目前 比较成熟和与广泛应用的理论模型主要有 DLA和 由式(6)～(12)利用 Matlab语言编程计算 (源程序略)可以得到 CCA两种。通过以上所建立的 CCA模