Matlab在信息光学中应用

1、2021/3/10讲解：XX1 Matlab在信息光学中的应用 一些衍射现象的计算机模拟 李祥艳-05102133 2021/3/10讲解：XX2 一、单缝衍射模拟 2021/3/10讲解：XX3 2021/3/10讲解：XX4 单缝衍射计算机模拟代码 clc; clear; a=-2*pi:0.0001*pi:2*pi; p1=(1-sinc(a).2;%方便下面着色 p2=sinc(a).2; figure; plot(a,p2); xlabel(kasin); ylabel(光强I/I0); title(单缝衍射强度分布); lgray=zeros(256,3); for i=0:255

2、 lgray(i+1,:)=(255-i)/255; end figure; imagesc(p1) title(单缝衍射模拟图); colormap(lgray) 2021/3/10讲解：XX5 多缝衍射 2021/3/10讲解：XX6 2021/3/10讲解：XX7 2021/3/10讲解：XX8 2021/3/10讲解：XX9 2021/3/10讲解：XX10 2021/3/10讲解：XX11 2021/3/10讲解：XX12 2021/3/10讲解：XX13 2021/3/10讲解：XX14 2021/3/10讲解：XX15 多缝衍射代码 clc; clear; n=input(多缝衍

3、射缝数); a=-2*n*pi:0.0001*pi:2*pi*n; p1=1-(sin(n*a)./sin(a).2;%方便下面着色 p2=(sin(n*a)./sin(a).2; figure; plot(a,p2); xlabel(kasin); ylabel(光强I/I0); title(多缝缝衍射强度分布); lgray=zeros(256,3); for i=0:255 lgray(i+1,:)=(255-i)/255; end figure; imagesc(p1) title(多缝衍射模拟图); colormap(lgray); 2021/3/10讲解：XX16 平面光栅衍射 2

4、021/3/10讲解：XX17 2021/3/10讲解：XX18 2021/3/10讲解：XX19 2021/3/10讲解：XX20 2021/3/10讲解：XX21 2021/3/10讲解：XX22 2021/3/10讲解：XX23 2021/3/10讲解：XX24 2021/3/10讲解：XX25 平面衍射计算机模拟代码 clc; clear; n=input(平面衍射光栅数); d=-n*pi:0.0001*pi:pi*n; b=d/15; p1=1-(sinc(b).*sin(n*d)./sin(d).2;%方便下面着色 p2=(sinc(b).*sin(n*d)./sin(d).2;

5、 figure; plot(d,p2); xlabel(kasin); ylabel(光强I/I0); title(多缝缝衍射强度分布); lgray=zeros(100,3); for i=0:99 lgray(i+1,:)=(99-i)/99; end figure; imagesc(p1) title(多缝衍射模拟图); colormap(lgray); 2021/3/10讲解：XX26 圆孔衍射 2021/3/10讲解：XX27 2021/3/10讲解：XX28 2021/3/10讲解：XX29 圆孔衍射代码 clear; N=1; K=6; x,y=meshgrid(linspace

6、(0,N+1,800); z=x+i*y; u=0; for m=1:N; for n=1:N; zk=abs(z-m+n*i)*K; u=u+0.1*besselj(4,zk)./zk; r=1-u; A=1-abs(u).2; end ip=imshow(A,) 2021/3/10讲解：XX30 矩孔衍射 2021/3/10讲解：XX31 2021/3/10讲解：XX32 矩孔衍射代码 x=zeros(30); x(10:20,10:20)=1 imshow(x,-1 5) x=fft2(x,256,256); x=fftshift(x); imshow(log(abs(x),-1 5);

7、 2021/3/10讲解：XX33 2021/3/10讲解：XX34 2021/3/10讲解：XX35 2021/3/10讲解：XX36 2021/3/10讲解：XX37 巴比涅互补原理 2021/3/10讲解：XX38 2021/3/10讲解：XX39 2021/3/10讲解：XX40 2021/3/10讲解：XX41 2021/3/10讲解：XX42 2021/3/10讲解：XX43 圆孔衍射代码 x,map=imread(yb.bmp); figure,imshow(x) figure,subplot(1,3,1);imshow(x) x=fft2(x,256,256); subplot

8、(1,3,2);imshow(log(abs(x),-1 5);title(未移至中心频谱) x=fftshift(x); subplot(1,3,3);imshow(log(abs(x),-1 5);title(移至中心频谱) 2021/3/10讲解：XX44 多缝衍射计算机模拟 2021/3/10讲解：XX45 2021/3/10讲解：XX46 2021/3/10讲解：XX47 2021/3/10讲解：XX48 旋转对频谱的影响代码 x1,map=imread(j1.bmp);figure;subplot(2,3,1);i mshow(x1);x1=fft2(x1,256,256);x1=

9、fftshift(x1);su bplot(2,3,4);imshow(log(abs(x1),-1 5); x1,map=imread(j2.bmp);subplot(2,3,2);imshow (x1);x1=fft2(x1,256,256);x1=fftshift(x1);subplot(2, 3,5);imshow(log(abs(x1),-1 5); x1,map=imread(j3.bmp);subplot(2,3,3);imshow (x1);x1=fft2(x1,256,256);x1=fftshift(x1);subplot(2, 3,6);imshow(log(abs(x1),-1 5); 2021/3/10讲解：XX49 缝宽对夫瑯禾菲衍射的影响 2021/3/10讲解：XX50 缝宽对夫瑯禾菲衍射的影响代码 igure; for k=1:4; u=zeros(50,50); u(20:30,25:(24+k)=1; subplot(2,4,k);imshow(u); end; for k=5:8; u=zeros(50,50); u(20:30,25:(20+k)=1; u=fft2(u,256,