球对称梯度折射率介质中光线传输轨迹的计算_石市委.pdfVIP

第14卷第12期 宜宾 学院 学 报 2014年12月 Vol.14,No.12 JournalofYibinUniversity December,2014 球对称梯度折射率介质中光线传输轨迹的计算 石市委，王佩红，杨 群，张 苗，刘 泉，孙兆奇 （安徽大学物理与材料科学学院，安徽合肥230039） 摘 要：通过将二阶的光线微分方程降阶为一阶微分方程组，并用四阶的Runge-Kutta方法对微分方程组进行迭代求解，从而 求出球梯度折射率介质中的光线传输轨迹. 该方法精确度高，适用于任意梯度分布的球梯度介质中的光线传输轨迹的计算. 关键词：梯度折射率；光线方程；光线轨迹 中图分类号：O439 文献标志码：A 文章编号：1671-5365(2014)12-0030-03 —————————————————————— TheRayTracinginMediumwithSphericalSymmetryGradientRefractiveIndex SHIShiwei,WANGPeihong,YANGQun,ZHANGMiao,LIUQuan,SUNZhaoqi (SchoolofPhysicsMaterialScience,AnhuiUniversity,Hefei,Anhui230039,China) Abstract: Byreducingtheraydifferentialequationtolowerorderequationsandusingthefourth-orderRunge-Kutta method,araytracingmethodwasobtainedforcalculatingthelighttraceinmediumwithsphericalsymmetrygradientre⁃ fractiveindex.Thismethodhashighprecisionandcanbeusedforraytracinginanysphericalsymmetrygradientrefrac⁃ tiveindexdistribution. Keywords:gradientrefractiveindex;raydifferentialequation;raytracing 梯度折射率光学介质因其具有均匀折射率光学 梯度介质进行光线传输轨迹的计算，对本科生和研 介质不能比拟的优异性能，近几十年来受到广泛关 究生在光学课程学习或梯度折射率器件的研究时进 注[1]. 在折射率梯度分布的介质中，光线的传输轨迹 行简单的梯度折射率器件性能的计算也有帮助. 不再是直线，而是按照光线微分方程的约束沿曲线 1 光线追迹迭代公式的推导 轨迹传输[2]. 光线微分方程是二阶的矢量微分方程， 光线在梯度折射率介质中传输必须依照光线微 通常只有极少数的几种梯度折射率分布曲线才有此 分方程进行，光线方程为[1]： 方程的解析解，如Maxwell鱼眼透镜模型、Lunerburg d dr [1-3] æ ön  =∇n (1) ds ds 透镜模型等 . 对于任意梯度的折射率分布介质， è ø 很难通过求解解析解的方法获得，通常是通过将微 这是一个二阶微分方程，式中r为描述光线轨 分方程降阶为两个一阶的微分方程组，然后通过严 迹的位置矢量，s为光线路程， 为光线所到达位置n 格的数值计算来求解[4]. 此过程在光学设计中被称 的折射率值. 对球对称梯度折射率介质，n可表示 为光线追迹，该过程的具体实施是较为复杂的. 在 为r 的函数，即n(r). 目前已知的球对称梯度折射率 通常的光学教材中，往往只是给出非常简单的过程 介质模型中，仅有Maxwell鱼眼模型和Luneburg介质 框架，学生很难通过这些