用MATLAB语言模拟光衍射实验

圆盘衍射matlab
在MATLAB中，我们可以使用以下步骤来模拟圆盘衍射：
1. 首先，我们需要定义一个圆形的光源。

这可以通过circle函数完成，例如：source circle = [0 0; 1 0];
2. 接下来，我们需要定义一个光栅，可以通过meshgrid 函数来实现。

例如：mesh grid = meshgrid(linspace(-1,1), linspace(-1,1), linspace(-1,1));
3. 然后，我们需要定义一个光源到光源的距离，例如：source to circle = 1;
4. 接着，我们需要定义一个光源到光栅的距离，例如：circle to grid = 1;
5. 最后，我们需要定义一个光源到光栅的距离，例如：circle to grid = 1;
6. 然后，我们可以通过以下公式来计算光栅的衍射图像：
diff image = (mesh grid.^2).*(source circle.^2).*(1 - (circle to grid.^2).*(source to circle.^2));
7. 最后，我们可以将这个公式放入MATLAB的代码中来实现圆孔衍射的模拟。

以上步骤只是一个基本的示例，具体的实现可能会根据你的具体需求有所不同。

phi = 0:pi/10:2*pi; % x的范围从0到2πz = exp(i * phi);u = real(z);%取实部v = image(z);%取虚部figure(1);%画图plot(phi,u);%画实部的变化曲线figure(2);plot(phi,v)%画虚部的变化曲线clear;%清除内存close all;%关闭所有窗口x=linspace(0,2*pi,100);%设置横坐标范围及网格y1=sqrt(x);%开根号figure(1);plot(x,y1,'r')%做出y1的曲线,'r'代表图线为红色clear;%清除内存close all;%关闭所有窗口x=linspace(-2*pi,2*pi,100);%100个数横坐标从-2pi到2pi y1=besselj(1,x);%第一类贝塞尔函数,一阶y2=besselj(0,x);%第一类贝塞尔函数,0阶figure(1);plot(x,y1,'r');figure(2);plot(x,y2,'b');clear;close all;[x,y]=meshgrid(-2:.05:2,-2:.05:2);%设置二维网格z=sinc(x).*sinc(y);figure(1);imshow(z);%二维灰度图figure(2);imagesc(z);%二维彩图figure(3);mesh(z);%三维图像clear;close all;a=5;x=linspace(-5,5,100);y1=exp(-pi.*(x/a).^2);%第一类贝塞尔函数,一阶b=2;y2=exp(-pi.*(x/b).^2);%第一类贝塞尔函数,一阶figure(1);plot(x,y1,'r');figure(2);plot(x,y2,'b');clear;close all;a=1;b=1;[x,y]=meshgrid(-2:.05:2,-2:.05:2);%设置二维网格z=exp(-pi.*((x/a).^2+(y/b).^2));%高斯函数figure(1);imagesc(z);figure(2);mesh(z);clear;close all;M=100;%x方向抽样点数N=100;%y方向抽样点数w=10;%缝的宽度d=10;%缝偏离中心的位置obj1=zeros(M,N);%赋值为0的一个二维矩阵,大小为M*N obj1(:,(N/2-w/2):N/2+w/2)=1.0;%设置缝在y向中心,且缝的宽度为wfigure(1);imshow(obj1);obj2=zeros(M,N);obj2(:,(N/2-d)-w/2:(N/2-d)+w/2)=1.0;%缝偏离中心位置d个像素figure(2);imshow(obj2);clear;close all;M=100;%x方向抽样点数N=100;%y方向抽样点数w=10;%缝的宽度d=10;%缝偏离中心的位置obj1=zeros(M,N);%赋值为0的一个二维矩阵,大小为M*N obj1(:,(N/2+d)-w/2:(N/2+d)+w/2)=1.0;%设置缝在y向中心,且缝的宽度为wobj2=zeros(M,N);obj2(:,(N/2-d)-w/2:(N/2-d)+w/2)=1.0;%缝偏离中心位置d个像素obj3=obj1+obj2;figure(1);imshow(obj3);clear;close all;M=100;N=100;w=10;%缝宽obj=zeros(M,N);obj((M/2-w/2):(M/2+w/2),(N/2-w/2):(N/2+w/2))=1.0;%缝处于中心位置figure(1);imshow(obj);figure(2);imagesc(obj);figure(3);mesh(obj);。

成绩：《工程光学》综合性练习二题目：基于matlab的衍射系统仿真学院精密仪器与光电子工程学院专业测控技术与仪器年级20**级班级**班姓名20**年**月综合练习大作业二一、要求3-4人组成小组，对下面给出的各题目利用Matlab等工具进行仿真。

练习结束时每组提交一份报告及仿真程序。

在报告中应注明各仿真结果所对应的参数，如屏与衍射屏间距、孔径形状尺寸等。

二、仿真题目1.改变观察屏与衍射屏间距，观察观察屏上发生的衍射逐渐由菲涅耳衍射转为夫琅和费衍射1)原理图：S点光源发出的波长lam=500纳米S点发出光线经过单缝，缝宽a；单缝到衍射屏的距离L'2)Matlab代码clear;clcl=10;%l=input('单缝到衍射屏的距离L=');a=0.2;%a=input('单缝的宽度(mm)a=');lam=500e-6;%lam=input('波长(nm)');x=-1:0.001:1;%接收屏边界y=x./sqrt(x.^2+l^2);z=a.*y/lam;I=1000*(sinc(z)).^2;%计算接受屏某点光强subplot(2,1,1)%绘制仿真图样及强度曲线image(2,x,I)colormap(gray(3))title('单缝衍射条纹')subplot(2,1,2)plot(x,I)title(光强分布)3)初始仿真图样（d=10）4)改变d之后的图样(d=1000)5)变化规律根据衍射屏以及接受屏的相对位置不同，由此产生菲涅尔衍射和夫琅禾费衍射的区别，根据我们模拟的情况得到菲涅尔衍射和夫琅禾费衍射的明显不同是夫琅禾费衍射条件下：中央有一条特别明亮的亮条纹，其宽度是其他亮条纹的两倍；其他亮条纹的宽度相等，亮度逐渐下降。

2.改变孔径形状、尺寸，观察图样变化1）原理图矩孔衍射：透镜焦距：1000mm;照射光波长：500nm;孔高：a（mm）；孔宽：b（mm）；圆孔衍射：圆孔直径：r（mm）；照射光波长：500nm;照射光波长：500nm;2）matlab代码矩孔衍射：focallength=1000;lambda=500;a=2.0;b=2.0;resolution=64;center=(resolution)/2;A=zeros(resolution,resolution);for i=1:1:resolutionfor j=1:1:resolutionif abs(i-center)<a*10/2&abs(j-center)<b*10/2 A(j,i)=255;endendendE=ones(resolution,resolution);k=2*pi*10000/focallength/lambda;imag=sqrt(-1);for m=1:1:resolutionx=m-center;for n=1:1:resolutiony=n-center;C=ones(resolution,resolution);for i=1:1:resolutionp=i-center;for j=1:1:resolutionq=j-center;C(j,i)=A(j,i)*exp(-imag*k*(x*p+y*q)); endendE(n,m)=sum(C(:));endendE=abs(E);I=E.^2;I=I.^(1/3);I=I.*255/max(max(I));L=I;I=I+256;CM=[pink(255).^(2/3);gray(255)];Colormap(CM);edge=(resolution-1)/20;[X,Y]=meshgrid([-edge:0.1:edge]);x=linspace(-edge,edge,resolution);y=linspace(-edge,edge,resolution);subplot(1,2,1);surf(x,y,L);axis([-edge,edge,-edge,edge,0,255]);caxis([0,511]);subplot(1,2,2);image(x,y,I);axis([-edge,edge,-edge,edge,0,511]);view(2);axis square;圆孔衍射：clearlmda=500e-9;%波长r=1.2e-3;%f=1;%焦距N=19;K=linspace(-0.1,0.1,N);lmda1=lmda*(1+K);xm=2000*lmda*f;xs=linspace(-xm,xm,2000);ys=xs;z0=zeros(2000);[x,y]=meshgrid(xs);for i=1:19s=2*pi*r*sqrt(x.^2+y.^2)./(lmda1(i));z=4*(besselj(1,s)./(s+eps)).^2;%光强公式z0=z0+z;endz1=z0/19;subplot(1,2,1)imshow(z1*255);%平面图xlabel('x')ylabel('y')subplot(1,2,2)mesh(x,y,z1)%三维图colormap(gray)xlabel('x')ylabel('y')zlabel('光强')3）仿真图样：矩孔衍射：a=1，b=2a=2，b=2可知：矩孔在一个维度上展宽一定倍数将导致衍射图样在相同维度上缩短相同倍数，同时能量会更向中心亮斑集中。

B(i,:)=(sin(alpha)./alpha).^2.*(sin(N
*beta)./sin(beta)).^2;
B1=B/max(B);
end
NC=255;
Br=(B/max(B))*NC;
subplot(1,2,1)
image(xs,ys,Br);
colormap(gray(NC));
subplot(1,2,2)
( t 2 r )
dE k( ) cos
ds
r
( t 2 r )
E(P) k( ) cos S
Hale Waihona Puke dsr这便是计算衍射场的常用公式，即惠更斯-菲涅尔原理的表达式。
单色平行光通过单缝,可看作N条光束的叠加来计算, 单缝衍射的几何关系如图所示.把单缝看作一排等间 隔的光源,N个光源分布,P处的光强为这N个光源照射 结果的合成。
波程差 相位差
=

2
光的衍射现象是光的波动性的一种表现，通过对光的各种衍射现象的研究可
以深入具体地了解光的波动性。惠更斯原理认为：传播中的波面上任何一点
都可以认为是一个新的次波源，由这些次波源发出的次波是球面波，这些次
波的公共包络面就是下一个时刻的波面。菲涅尔将其进一步具体化，认为可
以将波面S分成许多面积元ds，每个ds在P点的振动可写为：
N=256;
%确定灰度的等级
Br=(B/max(B))*N; %定标
subplot(1,2,1)
image(xm,xs,Br);
colormap(gray(N)); %色调处理
subplot(1,2,2)
plot(B,xs,'k');
杨氏双缝干涉