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

1、基于Matlab实现近场和远场衍射1 预备知识1.1 Matlab简介 Matlab是美国MathWorks公司出品的商业数学软件，用于算法开发、数据可视化、数据分析以及数值计算的高级技术计算语言和交互式环境。 Matlab可以进行矩阵运算、绘制函数和数据、实现算法、创建用户界面、连接其他编程语言的程序等，主要应用于工程计算、控制设计、信号处理与通讯、图像处理、信号检测、金融建模设计与分析等领域。1.2 傅里叶变换函数Matlab提供了一套计算快速傅里叶变换的函数，它们包含一维、二维和N维离散傅里叶变换函数fft、fft2和fftN。这些离散傅里叶变换的逆变换函数为ifft、ifft2和iff

2、tN。FFT变换完成后，低频部分分布在矩阵的四个角上，需要通过移频将零频及低频量移到矩阵中心，这样才能与建立的坐标系一致，移频通过Matlab自带函数fftshift()完成。2 基本思路要对圆孔进行衍射模拟，首先要设计出圆孔的图形。选取单位振幅的平面波垂直照射，一般取值为1，方向垂直xy平面，平行与z轴。确定观察屏距离孔径的距离。将上述的所有的量代入到公式中，计算出衍射的光场分布。Matlab模拟。3 菲涅尔衍射实例 由已知菲涅尔衍射中心点的光强为：由图可以看出菲涅尔衍射图样的中心光强和距离z有关，其亮暗是交替变化的。程序：%圆孔菲涅尔衍射clear;z0=5;%范围2到10 x0=0;%范

3、围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;endmax1=max(fFZP);max2=max(max1);scale=1.0/max2;fFZP=fFZP.*scale;R=127*(1+imag(fFZP);figure(1)image(R);axis square onaxis off图1 菲

4、涅尔圆孔衍射4圆孔的夫琅禾费衍射仿真 利用衍射积分或者傅里叶分析可以求得在观察屏上任意径向坐标r处的复振幅分布为：其中，光强分布为：程序：clear all%圆孔的夫琅禾费衍射lamda=632.8e-9;a=input(请输入孔径大小)f=1;Np=2048;%调色板n=400;m=400;ym=4000*lamda*f;ys=linspace(-ym,ym,n);xs=ys;for i=1:n for j=1:m r=sqrt(xs(i).2+ys(j).2); x=2*pi*a*r./(lamda*f); I(i,j)=(2*BESSELJ(1,x).2./x.2; endendIr=I

5、.*Np;subplot(1,2,1)image(Ir)colormap(gray)subplot(1,2,2)mesh(Ir);colormap(gray)Case 1孔径大小为a=0.0005Case 2孔径大小为a=0.001 Case 3孔径大小为a=0.015矩形的夫琅禾费衍射仿真 利用衍射积分或者傅里叶分析可以求得单色光照明时在观察屏上的复振幅分布为：光强分布为：程序：f=1;a=input(a=);b=input(b=);%矩孔大小lmda=632.8e-9;%波长xmax=0.008;%观察屏大小ymax=xmax;N=1500;x=linspace(-xmax,xmax,N)

6、;y=linspace(-ymax,ymax,N);for n=1:Nfor m=1:Nalpha=pi*x(n)*a/(lmda*f);beta=pi*y(m)*a/(lmda*f);I(n,m)=(sin(alpha)/(alpha)2*(sin(beta)/(beta)2;endendX,Y=meshgrid(x,y);Np=512Ir=Np*I;figuremesh(X,Y,Ir);colormap(gray);axis square;zlabel(光强);hold onfigureimshow(Ir);Case 1 a=b=0.001时 Case 2 a=b=0.01时 由例二可以发

7、现当孔径增大到0.01mm尺度时，衍射现象已经不是很明显，当孔径继续增大时小孔不再衍射，光线沿直线传播。6 任意图形Np=2048;a=imread(C:Users帅Desktop1.bmp);grid onfigure(1)imshow(a,)afft=fft2(a);%傅里叶变换aabs=abs(afft);%取模aabss=fftshift(aabs);%取模后的矩阵进行频谱位移I=aabss.2;Ir=I.*Np;figureimshow(Ir,)colormap(gray);figuremesh(Ir);figureplot(Ir);Case 1 Case 2 Case 3 7 结论：1实现了利用FFT模拟菲涅尔衍射和夫琅禾费衍射。2利用模拟得到的图像可以直观地理解衍射。3利用模拟可以更加方便地变换参数以及计算其影响。4 对衍射现象的模拟可以是在教学中，将抽象难懂的概念、规律通过实验仿真生动、形象地表现出来，这样可以拓宽知识面，提高能力。内容总结（1）基于Matlab实现近场和远场衍射1 预备知识1.1 Matlab简介 Matlab是美国MathWorks公司出品的商业数学软