Matlab数字衍射光学实验.

1、1)复指数函数cS=cos(0)+ Zsin(0),实部 cos(p),虚部 sin(p), phi = 0:pi/10:2*pi; % x 的范围从 0 到 2 兀z = exp(i * phi);u = real(z);% 取实部v = image;取虚部figure(l);% 画图plot(phi,u);%画实部的变化曲线figure(2);plot(phi,v)%画虚部的变化曲线2)开根号：sqrt (x),clear;%清除内存close all;%关闭所有窗口x=linspace(0,2*pi,100);%设置横坐标范围及网格 yl=sqrt(x);% 开根号figure(l);p

2、lotfxl/r1)%做出yl的曲线斤代表图线为红色3)函数：一维sine (x)函数,(sin(/rx/G)sme = a ) 7rx/ aclear;%清除内存 close all;%关闭所有窗口 x=linspace(-2*pi,2*pi,100);%100 个数横坐标从2pi 到 2pi yl=besselj(l,x);%第一类贝塞尔函数，一阶 y2=besselj(0,x);%第一类贝塞尔函数,0阶 figure(l);plot(x,y门;figure(2);plot(x,y2,b)；0.60.40.20 - 02 -0.4 - -0.6 -88450246二维sine函数:/Xy、

3、/ Vsine25=sine sineJja, b为正数。clear;close all;二维网格x,y=meshgrid(-2:.05:2,-2:.05:2);% 设:z=si nc(x).*si nc(y); figure(l);imshow;二维灰度图 figure(2);imagesc;二维彩图 figure(3);mesh(z);%三维图像102030405060708010203040506070803）高斯函数:维高斯函数Gauss/ 、( X= exp-TVX2丿J J,a, b为正数。Gauss 二工U b /GaussGauss bexp-7Tclear;close all

4、;a=l;b=l;x,y=meshgrid(-2:.05:2,-2:.05:2);% 设置二维网格 z=exp(-pi.*(x/a).A2+(y/b).A2);% 高斯函数 figure(l);imagesc(z);figure(2);mesh(z);203040506070801020304050607080100常用物体的Mat lab模拟1)单缝clear;close all;M=100;%x方向抽样点数N=100;%y方向抽样点数w=10;%缝的宽度d=10;%缝偏离中心的位置objl=zeros(M,N);%赋值为0的一个二维矩阵,大小为M*N objl(:z(N/2-w/2):N/

5、2+w/2)=1.0;%设置缝在y向中心，且缝的宽 度为wfigure(l);imshow(objl);obj2=zeros(M,N);obj2(:/(N/2-d)-w/2:(N/2-d)+w/2)=1.0;%缝偏离中心位置 d 个像素figure(2);imshow(obj2);2)双缝在单缝的基础上，编程两个不同位置的单缝，然后相加，即可得出双缝clear;close all;M=100;%x方向抽样点数N=100;%y方向抽样点数w=10;%缝的宽度d=10;%缝偏离中心的位置objl=zeros(M,N);%赋值为0的一个二维矩阵,大小为M*N objl(:z(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=objl+obj2;figure(l);imshow(obj3);3)矩形孔clear;close all;M=100;N=100;w=10;%缝宽obj=zeros(M,N);obj(M/2-w/2):(MZ2+w/2)/(NZ2-w/2):(NZ2+w/2)=1.0;% 缝处于中心位置figure(l);ims