• p1=1-(sinc(b).*sin(n*d)./sin(d)).^2;%方便下面着色
• p2=(sinc(b).*sin(n*d)./sin(d)).^2;
• figure;
• plot(d,p2);
• xlabel('kasinθ');
• ylabel('光强I/I0');
• title('多缝缝衍射强度分布');
讲解：XX
26
2021/3/10
讲解：XX
27
2021/3/10
讲解：XX
28
圆孔衍射代码
• clear; • N=1; • K=6; • [x,y]=meshgrid(linspace(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,[])
讲解：XX
20
2021/3/10
讲解：XX
21
2021/3/10
讲解：XX
22
2021/3/10
讲解：XX
23
2021/3/10
讲解：XX
24
平面衍射计算机模拟代码
• clc;
• clear;
• n=input('平面衍射光栅数');
• d=-n*pi:0.0001*pi:pi*n;
• b=d/15;
2021/3/10
讲解：XX
13
2021/3/10
讲解：XX
14
多缝衍射代码
• clc;
• clear;
• n=input('多缝衍射缝数');
• 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;
2021/3/10
讲解：XX
4
多缝衍射
2021/3/10
讲解：XX
5
2021/3/10
讲解：XX
6
2021/3/10
讲解：XX
7
2021/3/10
讲解：XX
8
2021/3/10
讲解：XX
9
2021/3/10
讲解：XX
10
2021/3/10
讲解：XX
11
2021/3/10
讲解：XX
12
45
2021/3/10
讲解：XX
46
2021/3/10
讲解：XX
47
旋转对频谱的影响代码
• [x1,map]=imread('j1.bmp');figure;subplot(2,3,1);i mshow(x1);x1=fft2(x1,256,256);x1=fftshift(x1);su bplot(2,3,4);imshow(log(abs(x1)),[-1 5]);
2021/3/10
讲解：XX
32
2021/3/10
讲解：XX
33
2021/3/10
讲解：XX
34
2021/3/10
讲解：XX
35
2021/3/10
讲解：XX
36
巴比涅互补原理
2021/3/10
讲解：XX
37
2021/3/10
讲解：XX
38
2021/3/10
讲解：XX
39
2021/3/10
Matlab在信息光学中的应用
——一些衍射现象的计算机模拟 李祥艳-05102133
2021/3/10
讲解：XX
1
一、单缝衍射模拟
2021/3/10
讲解：XX
2
2021/3/10
讲解：XX
3
单缝衍射计算机模拟代码
• 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 • lgray(i+1,:)=(255-i)/255; • end • figure; • imagesc(p1) • title('单缝衍射模拟图'); • colormap(lgray)
5]);title('未移至中心频谱') • x=fftshift(x); • subplot(1,3,3);imshow(log(abs(x)),[-1
5]);title('移至中心频谱')
2021/3/10
讲解：XX
43
多缝衍射计算机模拟
2021/3/10
讲解：XX
44
2021/3/10
讲解：XX
讲解：XX
40
2021/3/10
讲解：XX
41
2021/3/10
讲解：XX
ቤተ መጻሕፍቲ ባይዱ42
圆孔衍射代码
• [x,map]=imread('yb.bmp'); • figure,imshow(x) • figure,subplot(1,3,1);imshow(x) • x=fft2(x,256,256); • subplot(1,3,2);imshow(log(abs(x)),[-1
• 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
讲解：XX
15
平面光栅衍射
2021/3/10
讲解：XX
16
2021/3/10
讲解：XX
17
2021/3/10
讲解：XX
18
2021/3/10
讲解：XX
19
2021/3/10
2021/3/10
讲解：XX
29
矩孔衍射
2021/3/10
讲解：XX
30
2021/3/10
讲解：XX
31
矩孔衍射代码
• 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]);
• lgray=zeros(100,3);
• for i=0:99
•
lgray(i+1,:)=(99-i)/99;
• end
• figure;
• imagesc(p1)
• title('多缝衍射模拟图');
• colormap(lgray);
2021/3/10
讲解：XX
25
圆孔衍射
2021/3/10