基于Matlab的光学实验仿真Word格式.docx
- 文档编号:20602078
- 上传时间:2023-01-24
- 格式:DOCX
- 页数:69
- 大小:2.20MB
基于Matlab的光学实验仿真Word格式.docx
《基于Matlab的光学实验仿真Word格式.docx》由会员分享,可在线阅读,更多相关《基于Matlab的光学实验仿真Word格式.docx(69页珍藏版)》请在冰豆网上搜索。
应用金斯莱克(Kingslak)论述的方法对包含像差的典型双光束干涉图样进行了仿真分析。
并在此基础上讨论了各类像差对干涉条纹的影响。
(4)光学V系统仿真。
讨论了利用光学4,系统进行光学信息处理的基本原理,并据此提出了仿真光学够系统的算法。
孀仿真的铲系统实现了网格滤波以及图像相加、相减等光学信息处理实验。
(5)设计制作光学实验仿真软件界面。
针对Matlab应用程序无法脱离其运行环境的缺点.结合Mateom4.5将Matlab的M文件转化成相同功能的C++代码,然后在Visualc++环境下将其编译成外部可执行文件,使整个软件可脱离Matlab环境运行。
关键词:
光学实验仿真,干涉,衍射,像差,傅里叶变换,Matlab,Matcom4。
5,VisualC++。
・国家“十五”重点教材建设项目资助
OpticalExperimentSimulationsWithMatlab
Abstract
Inrecentdays,opticalexperimentsimulations,asanovelalternativemethoddifferentfromtheroyandexperimenttoexplorenewthings,areattractingmoreandmoreattentionfrombothscientistsandteachers.Basedonmanysimulativeexperiments,wepresenttheopticalsimulativealgorithmandresultsbyusingMatlabprogramandgiveanintroductiontoalltheresultsfromthesireulations.ThemainworkiSasfollows:
(1)First,simulationsofopticalinterferencearepresented.Twosphericalwavesandmultiplesphericalwavesinterference,andYang’Sinterferenceexperimentsaresimulatedbyusingdifferentmethods.Twospheficalwavesinterferenceisfocusedonthe0n-planeorout—planeofthespotlightsourceplaneandtheobservingplane.Then,weanalyzehowthetemporalcoherenceofopticalwavesinfluencethecontrastoftheinterferencepattem.Furthermore,equalintenseandunequalintensesphericalwavesinterferencesareinvestigated.Finally,Yang’Sdoubleslitsandcircleaperturessimulativeresultsaredetailedlypresented.
(2)Second,simulationsofopticaldiffractionarepresented.TheFresnelandFranhuferdiffractionpatterns.ofregularshapesscreensuchasslits,rectangularapertures,circleapertures,discs,triangleapertures,Ronchigratings,andsomeirregularshapesscreen,areobservedinthedifferentpropagationdistance.ThisclearlygivestherelationanddifferenceofthetWOdiffractionpatterns.
(3)Third,opticalaberrationsaresimulated.Basedonphase-intensitytransformprinciple,thesimulationresultsofSeidelaberrationandZernikepolynomialareprovided,AccordingtotheKingslak’Stheory,thetypicaldoublebeamsinterferencepatternsinfluencedbytheaberrationsarefurtheranalyzed.
(4)Forth,simulationsof4fopticalsystemarepresented.Theopticalinformationprocessingsuchasgriddlefilterandimageplus&substractarerealizedbyusingthesystem.
(5)Finally,asoftisprogrammedforthesimulationoftheopticalexperiment.Itconsistsofalltheabovesimulativeprograms,anddosenotdependontheenvironmentofMatlaborVisualC++,BasedontheMatcom4.5.nlem.filesofMatlabaretranslatedintoC++codes.thenasuitableapplicationprograminterfacearedesignedbyusingVisualC++.Keywords:
Opticalexperimentsimulations;
Interference;
Difhaction;
Aberrations;
Fouriertransform;
Matlab;
Matcom4.5;
VisualC++;
II
第一章绪论
1.1光学实验仿真的国内外研究现状
在计算机飞速发展的今天,光学实验仿真受到越来越多的科研工作者和教育工作者的广泛关注。
其应用主要有两个方面:
第一是在科学计算方面,利用仿真实验的结果指导实际实验,减少和避免贵重仪器的损伤:
第二是在光学教学方面,将抽象难懂的光学概念和规律,由仿真实验过程直观地描述,让学生饶有兴趣的掌握知识。
在科学计算方面,国外的光学实验仿真是在模拟设计和优化光学系统的过程中发展起来的。
在这方面,美国走在最前面,其中最具代表性的是劳伦斯利弗莫尔实验室光传输模拟计算软件Prop92【1,2】及大型总体优化设计软件CHAINOP和PROPSuITEn其中软件Prop92是采用Fortran语言编写,最大采样点数为512"512。
这些软件在大型激光器NOVA的设计制造中发挥了重要作用。
另外法国也开发完成其具有自身特点的光传输软件Miro[41。
该软件采用c++编写,可以运行于多个平台。
俄罗斯普通物理实验室在上世纪90年代初推出光传输软件Fresnel2.0t”,经过几年的发展,目前已经推出Fresnel4.0版本。
该软件是一套成熟的商业化光传输软件,在处理光传输的数值模拟算法上有独到之处,主要体现在其快速傅里叶变换的计算效率很高;
软件采用特殊方法能够处理小于计算分辨率的灰尘点的衍射过程以及截止频率小于计算网格分辨最小频率的滤波过程等。
另外,该软件图形显示界面友好,运行稳定可靠,如图1.1所示。
图I.1Fresnel软件界蕊
我国用于科学研究的光学实验计算机数值仿真软件虽开发较晚,但也已经取得了显著成绩。
特别是1999年,神光一III原型装置TIL分系统集成实验的启动为高功率固
体激光驱动器的计算机数值模拟的研究创造了条件。
目前已基本完成的SG99t6】光传输
模拟计算软件的开发,推出的标准版本基本能稳定运行,对SG99主要计算模块的验证结果表明SG99对能流放大、线性传输、非线性传输的计算是合理可靠的,其中线性传输的计算模块的计算精度与国外同类软件Fresnel十fl当:
目前该软件已经应用于神光一III主机可行性论证的工作中。
在光学教学方面,国外已有相关的配有光盘演示光学实验的教材【7l’该教材主要针对高年级学生和研究生使用。
其中不仅详尽的介绍了几何光学、物理光学、光学成像技术及图像处理技术,而且利用现在普遍使用的软件工具Matlab对它们进行了系统的仿真。
也有针对理科和工科低年级学生使用的光学教材【81,该教材使用MatchCad绘制各种逼真的光学仪器.创造出仿真的光学实验室。
学生可利用其进行探索和发现性学习,充分调动学生的积极性。
还有网络版光学教材【9J,该教材采用Mathematica进行光学仿真计算,结合LiveGraphic3DJaval.1的动画制作功能在网络上实时演示各种光学实验的结果图。
我国光学教材在利用计算机仿真方面相对落后,至今没有同类教材出现。
在2003年北京举行的网络教育软件展上,有关光学实验的网络教学软件都偏重于理论分析方面,对计算机应用于光学实验的仿真方面未给与充分重视。
结合国家十五教材建设计划,在光学实验仿真方面进行大量的研究,各项研究工作将在后续各章中~一介绍。
1.2光学实验仿真研究的意义
在工程设计领域中,人们通过对研究对象建立模型,用计算机程序实现系统的运行过程和得到运算结果,寻找出最优方案,然后再予以物理实现,此即为计算机仿真科学¨
0l。
在计算机同益普及的今天,计算机仿真技术作为虚拟实验手段已经成为计算机应用的一个重要分支。
它是继理论分析和物理实验之后,认识客观世界规律性的一种新型手段。
计算机仿真过程是以仿真程序的运行来实现的。
仿真程序运行时,首先要对描述系统特性的模型设置一定的参数值,并让模型中的某些变量在指定的范围内变化,通过计算可以求得这种变量在不断变化的过程中,系统运动的具体情况及结果。
仿真程序在运行过程中具有以下多种功能【111:
(1)计算机可以显示出系统运动时的整个过程和在这个过程中所产生的各种现象和状态。
具有观测方便,过程可控制等优点;
(2)可减少系统外界条件对实验本身的限制,方便地设置不同的系统参数,便于研究和发现系统运动的特性;
(3)借助计算机的高速运算能力,可以反复改变输入的实验条件、系统参数,大大提高实验效率。
因此.计算机仿真具有良好的可控制性(参数可根据需要调整)、无破坏性(不会因为设计上的不合理导致器件的损坏或事故的发生)、可复现性(排除多种随机因
素的影响,如温度、湿度等)、易观察性(能够观察某些在实际实验当中无法或者难以观察的现象和难以实现的测量,捕捉稍纵即逝的物理现象,可以记录物理过程的每一个细节)和经济性(不需要贵重的仪器设备)等特点。
在光学仪器设计和优化过程中,计算机的数值仿真已经成为不可缺少的手段。
通过仿真计算,可以大幅度节省实验所耗费的人力物力,特别是在一些重复实验工作强度较大且对实验器材、实验环境等要求较苛刻的情况下.如在大型激光仪器的建造过程中,结合基准实验的仿真计算结果可为大型激光器的设计和优化提供依据。
仿真光学实验也可应用于基础光学教学。
光学内容比较抽象,如不借助实验,学生很难理解,如光的干涉、菲涅耳衍射、夫琅禾费衍射等。
国外著名的光学教材I坨1配有大量的图片(包括计算和实验获得的图片),来形象地说明光学中抽象难懂的理论。
光学实验~般需要稳定的环境.高精密的仪器,因此在教室里能做的光学实验极为有限,而且也受到授课时间的限制。
为了克服光学实验对实验条件要求比较苛刻的缺点,可采用计算机仿真光学实验,特别是光学演示实验,配合理论课的进行.把光学课程涉及的大多数现象展示在学生面前,以加深对光学内容的理解。
如利用计算机仿真联合变换相关实验,可以得到清晰的相关峰,而在实验中液晶光阀的分辨率较低,很难得到清晰的相关峰:
又如光学菲涅耳衍射与夫琅禾费衍射,初学者不易理解,如果通过光学仿真实验,可以计算出它们之间的演化规律,清楚地说明二者之间的联系与区别。
学生们可以根据对光学原理和规律的理解,自己设置在仿真光学实验中的可控参数,探索和发现光学世界的奥秘,调动学习的积极性。
1_3MATLAB语言用于计算机仿真的优势
Matlab是MathWorks公司于1982年推出的一套高性能的数值计算和可视化软件。
它集数值分析、矩阵运算、信号处理和图形显示于一体,构成了一个方便、界面友好的用户环境。
它还包括了ToolBox(工具箱)的各类问题的求解工具,可用来求解特定学科的问题。
其特点旧”'J副是:
(1)可扩展性:
Matlab最重要的特点是易于扩展,它允许用户自行建立指定功能的M文件。
对于一个从事特定领域的工程师来说,不仅可利用Matlab所提供的函数及基本工具箱函数,还可方便地构造出专用的函数。
从而大大扩展了其应用范围。
当前支持Matlab的商用Toolbox(T具箱)有数百种之多。
而由个人开发的Toolbox则不可计数。
f21易学易用性:
Matlab不需要用户有高深的数学知识和程序设计能力,不需要用户深刻了解算法及编程技巧。
(3)高效性:
Matlab语句功能十分强大,一条语句可完成十分复杂的任务。
如m语句可完成对指定数据的快速傅里叶变换,这相当于上百条C语言语句的功能。
它大
第一聿绪论
大加快了工程技术人员从事软件开发的效率。
据MathWorks公司声称,Matlab软件中所包含的Matlab源代码相当于70万行C代码。
由于Matlab具有如此之多的特点,在欧美高等院校,Matlab已成为应用于线性代数、自动控制理论、数理统计、数字信号处理、时间序列分析、动态系统仿真等高级课程的基本教学工具:
在研究单位、工业部门,Matlab也被广泛用于研究和解决各种工程问题。
当前在全世界有超过40万工程师和科学家使用它来分析和解决问题。
1.4本文的主要工作
本文的主要目的是探讨利用Matlab进行光学实验仿真的方法及实现程序。
在大量仿真实验的基础上建立基于Matlab的光学实验仿真系统。
各章的具体内容如下:
第二章:
介绍了Matlab用于光学实验仿真的基本理论,分析利用Matlab实现光波波前、光束传播、规则图形绘图、图像显示等的可行性。
第三章:
讨论光的干涉实验仿真的理论基础和仿真方法,详细给出两列球面波的干涉、多光束干涉、杨氏实验的仿真方法及仿真结果。
第四章:
讨论光的衍射实验仿真的理论基础和仿真方法,具体给出菲涅耳衍射和夫琅禾费衍射的仿真方法及仿真结果。
第五章:
讨论光学波像差仿真的理论基础和仿真方法,具体给出赛德尔像差和泽尼克多项式以及组合像差的仿真结果,分析波像差对干涉条纹的影响。
第六章:
在前述几章工作的基础上详细讨论光学可系统的仿真,以及基于移系统的光学滤波和相干光学信息处理等的仿真。
第七章:
制作了初步应用程序界面,给出光学实验仿真系统。
在Matcom4.5的基础上,探讨如何将Matlab的M文件转化成相同功能的C++代码,然后应用Visual给出一个合理的应用程序界面。
最后为结论与展望,总结本文的一些主要工作与思路.并对后续工作方向提出一些个人的看法。
c十+
西北工业大学顶士学位论文
第二章MATLAB用于光学实验仿真的基本理论
在光学实验的仿真过程中,理解波前的概念是十分必要的。
“波前”一词,过去人们常指的是一个等相面(波面),或走在最前箍的波面,后一种含义只对冲击波一类非定态波有意义【l“。
研究定态托波时,“波前”~词泛指波场中任~曲面,更多地指一个平面,如记录介质、感光底片、接收屏幕、透明的黑自画面等所在的平面,或透镜前后的某个平面。
实际问题中,一列波携带着许多信息,如频率、波长和传播方向(二者包含在波矢k中)、振幅分布、位相分布、传播速度等信息全部包含在三维的复振幅分布函数中。
然而,通常光学系统中的一个元件只和波场中某个波前打交道.也就是说,与它有关的只是这个波前上的信息,所以通常只关心复振幅在波前上的二维分布问题。
在光学实验仿真过程中.需要解决的主要问题是如何准确地描述复振幅在某个空侧平面上的离散分布。
计算机处理的图像是离散化的,所谓图像离散化,就是把图像分割成如图2.1(b)所示的称为像素的小区域,每个像素的亮度或灰度值(graylevelvalue)用一个整数来表示‘”】。
一幅NxN个像素的数字图像.其灰度值可用Ⅳ行Ⅳ列的矩阵U来表示,如图2.1(c)所示。
这样,对数字图像的各种运算处理就变成对矩阵u的各种运算。
ZI
U=一:
五,厶,
厶.厶:
f曲Picturetb]Digltalpicture
(Ⅳ×
Ⅳ)pixdn1Ma”“(Ⅳ一Ⅳ)
图2.1图像的离散化和矩阵表示
2.】利用MATLAB语言实现光波波前
Matlab语言描述光学波前的优越性具体体现在三个方面:
是Maflab语言可以方便地确定光学平面,特别是同轴平面,并能方便地描述光学平面所满足的函数关系;
二是Matlab数学函数库提供了计算光束线性传播所需要的二维快速离散傅里叶变换;
三是Matlab语言具有强大的绘图功能,可以方便地显示各种光学场的计算结果。
2.1.1剥用Matl西语言实现光波波前
二维x—y直角坐标可唯一地确定一个平面(垂直于z轴的平面),用这个平面表示光波场的波前.可以在此平面上对光波场的复振幅作相应的各种运算。
用Mafl:
,b语言
可以方便确定这样一个平面。
利用Maflab庠函数meshgridOflal生成两个N×
N维矩阵
变量rn和n,其中矩阵m按行变换,矩阵n按列变换,按照这两个矩阵寻址可确定用矩阵U表示的平面上的任意像素,如图2.2所示。
可以根据需要确定矩阵的阶数Ⅳ,即平面的采样点数。
1
212l2l12….Ⅳ2…Ⅳ
NNL——————]r————J
工.
U=…ⅣJL12一:
^,^:
{q~{Hl…{Mq
由库函数meshgridO确定的各个平面的坐标原点位于平面的中心,使所有参与运算的平面之间保持很好的同轴性。
调用库函数cart2pol0[’91可进行直角坐标和极坐标之间的相互转换。
另外矩阵u中的元素可以为复数,且其中的所有元素都可进行加、减、乘、除等运算,由此可解决复振幅在波前上二维分布的计算问题。
2.1.2平面光波的波前实现
平面光波是光波的一种最简单形式(201。
在某一时刻t,具有相同位相的波动构成~组平面,每个平面都与波动的传播方向垂直。
点光源发出的光波经透镜准直,或者把点光源移到无穷远,可以近似获得平面波。
沿^方向传播的单色平面波.在光场中P(x,Y,z)点产生的复振幅可以表示为
u(x,Y,z)=aexp【ik(xcosct+ycosfl+zcosy)】
COS2Gr+COS2fl+COS2y=l
由此可将式(2.1.1)改写为:
(2.1.1)式中,a表示常量振幅,COSa'、cosfl、cosg'J口传播方向的方向余弦,它们满足:
(2.1.2)
吣∥,z’i=轰篆芹粤COS鬟器COS五荔矬。
。
口+心。
蒯cz^。
,=口expD—切√1—2口一2∥居xp【j七(xcos口+ycos∥)】
对于在确定方向上传播的平谣波,以及所选定的垂直z轴的x-y平面,上式中的第一个位相因子exp[jkz(1.COS2G-COS2历“2]是常数位相因子,不随x∥平面坐标变化。
因此,可以引入一个复常数4。
令
A=口e×
p【jk√1一cos2口一cos2pJ
以及平面波传播方向余弦COSLZ、cosp而变化。
(2.1.4)式中,位相因子表示平面波在Ⅵ,平面上产生均匀相移。
其大小随着x-y平面位置z,
西北工业大学硕士学位论文7
x-y平面上的复振幅分布可以表示为
£,G,y)=Aexp[jk(xeosa+yeosfl)](2.1.5)
通常称exp[jk(xcosa+ycos#)]为平面波的线性位相因子。
若平面上的复振幅分布的表达式中包含这一因子,可知它代表一个方向余弦COStz、cosfl的平面波经过该平面,见国2-3(a)。
等位相线的方程是
XCOS口+ycos口=C(2.16)
式中,C为某一常量。
不同C值所对应的等位相线是一些平行斜线,见图2.3(b)。
嫌夕
k。
~
(a)@fl。
’rh、
图2.3(a)传播过程中的平面波波前:
(b)平面波在x-y平面上的等位相线
要实现(2.1,5)式的离散化,可利用Matlab库函数meshgridO产生两个矩阵变量m和n。
其中矩阵m按行变换,代表离散化的x坐标变量;
矩阵玎按列变换,代表离散化的Y坐标变量。
按照这两个矩阵寻址来确平面光波波前的矩阵。
其中振幅A、方向余弦cost:
涌:
lcosfl为常量,与坐标变量无关。
设平面上x、Y方向的采样点数分别为肼和Ⅳ,则采样间距缸、缈分别为:
Ax=‘砌Ay2%r/I/
(2-1-7)
式中L和4分别为X和y方向的积分范围。
离散化的平面光波波前的表达式为:
(厂∞Ax,"缈)=Aexp[jkm缸cos口+尬ycosp)】
(2.1.s)
式中m、行及U均为M×
N维矩阵。
由相应程序实现的平面光波波前如图2.4所示。
(a】二维包裹位相(b)连续位相(c)十拶位福
图2.4与z轴和y轴夹角为450的平面光波波前位相图
2.1.3球面波的波前实现
从点光源发出的光,其波面表现为球面,因此球面是基本的波面形式f2D】。
单色的
发散球面波在光场中任意一点P所产生的复振幅
V(P)=aexp[jkr】
光源的距离。
口表示距点光源单位距离处的振幅。
对于会聚球面波,则有(2.1.9)式中,波数k=-2冗/A,表示单位长度上产生的位相变化。
,表示观察点P(x,Y,z)到点
or(e)=-%xp[-j鼢】
点光源或会聚点位于坐标原点,(2.1.10)
r:
仁2+y2+z2y”
点光源或会聚点位于空间任意一点S(xo,Yo,zo)时,(2.1.11)
振i再石葡可二习(z.1.12)
图2.5球面波在xy平NY-O向等位相线
考虑光波在某一特定平面上产生的复振幅分布的数学描述,见图2.5。
设点光源S(工o,Yo,O)位于xo-yo平面,考察与其相距:
(z>O)的x-y平面上光场分布。
,可以写为:
,:
推可二≯硼=z
—(x-xoY+—@-yoY<<12‘(2.1.13)若在x吵平面上只考虑一个对S点张角不大的范围,有
利用二项式展开,并略去高阶项,得到
,。
+垒=
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 基于 Matlab 光学 实验 仿真