Matlab在信息光学中应用

矩孔衍射

矩孔衍射代码
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]);

巴比涅互补原理

圆孔衍射代码
[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 5]);title('未移至中心频谱') x=fftshift(x); subplot(1,3,3);imshow(log(abs(x)),[-1 5]);title('移至中心频谱')

缝宽对夫瑯禾菲衍射的影响

缝宽对夫瑯禾菲衍射的影响代码
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,256);u=fftshift(u); subplot(2,4,k);imshow(log(abs(u)),[-1 1]);title('夫瑯禾费衍射'); end;
ห้องสมุดไป่ตู้
多缝衍射

多缝衍射代码
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; 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);
Matlab在信息光学中的应用
——一些衍射现象的计算机模拟 李祥艳-05102133

一、单缝衍射模拟

单缝衍射计算机模拟代码
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)

空域位移对夫瑯禾菲衍射的影响



平面光栅衍射

平面衍射计算机模拟代码
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; 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);

多缝衍射计算机模拟

旋转对频谱的影响代码
[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]); [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]);

圆孔衍射

圆孔衍射代码
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,[])