计算机实现近场和远场衍射.pdfVIP

基于 Matlab 实现近场和远场衍射 1 预备知识 1.1 Matlab 简介 Matlab 是美国 MathWorks 公司出品的商业数学软件，用于算法 开发、数据可视化、数据分析以及数值计算的高级技术计算语言和交 互式环境。 Matlab 可以进行矩阵运算、绘制函数和数据、实现算法、创建用 户界面、连接其他编程语言的程序等，主要应用于工程计算、控制设 计、信号处理与通讯、图像处理、信号检测、金融建模设计与分析等 领域。 1.2 傅里叶变换函数 Matlab 提供了一套计算快速傅里叶变换的函数，它们包含一维、 二维和 N 维离散傅里叶变换函数 fft 、fft2 和 fftN 。这些离散傅里叶变 换的逆变换函数为 ifft、ifft2 和 ifftN。 FFT 变换完成后，低频部分分布在矩阵的四个角上，需要通过移 频将零频及低频量移到矩阵中心，这样才能与建立的坐标系一致，移 频通过 Matlab 自带函数fftshift()完成。 2 基本思路 （1）要对圆孔进行衍射模拟，首先要设计出圆孔的图形。 （2 ）选取单位振幅的平面波垂直照射，一般取值为 1，方向垂直xy 平面，平行与 z 轴。 （3 ）确定观察屏距离孔径的距离。 （4 ）将上述的所有的量代入到公式中，计算出衍射的光场分布。 （5 ）Matlab 模拟。 3 菲涅尔衍射实例 由已知菲涅尔衍射中心点的光强为： 2 2  k  2    I 0, 0  U 0, 0  4 sin   4 sin    4z   2z  由图可以看出菲涅尔衍射图样的中心光强和距离z 有关，其亮暗 是交替变化的。 程序： %圆孔菲涅尔衍射 clear; z0=5;%范围 2 到 10 x0=0;%范围 ROWS=256; COLS=256; colormap(gray(255)) sigma=1/z0; y0=-x0; y=-12.8; for r=1:COLS, x=-12.8; for c=1:ROWS, fFZP(r,c)=exp(j*sigma*(x-x0).^2+j*sigma*(y-y0).^2); x=x+0.1; end y=y+0.1; end max1=max(fFZP); max2=max(max1); scale=1.0/max2; fFZP=fFZP.*scale; R=127*(1+imag(fFZP)); figure(1) image(R); axis square on axis off 图 1 菲涅尔圆孔衍射 4 圆孔的夫琅禾费衍射仿真 利用衍射积分或者傅里叶分析可以求得在观察屏上任意径向坐标 r 处的复振幅分布为： ka2  kr2   2J 1 kar / z   U r   exp j kz  exp j     j 2f  2f   kar / z  其中 2 2 r  x  y 光强分布为： 2 2  ka2   2J 1 kar / z   I r      