Matlab在信息光学中应用.ppt

1、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(单缝衍射模拟图

2、); colormap(lgray),多缝衍射,多缝衍射代码,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) titl

3、e(多缝衍射模拟图); 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,:)=(9

4、9-i)/99; end figure; imagesc(p1) title(多缝衍射模拟图); colormap(lgray);,圆孔衍射,圆孔衍射代码,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,),矩孔衍射,矩孔衍射代码,x=zeros(30); x(10:20,10:20)=1 imshow(x

5、,-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(移至中心频谱),多缝衍射计算机

6、模拟,旋转对频谱的影响代码,x1,map=imread(j1.bmp);figure;subplot(2,3,1);imshow(x1);x1=fft2(x1,256,256);x1=fftshift(x1);subplot(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);,缝宽对夫瑯禾菲衍射的影响,缝宽对夫瑯禾菲衍射的影响代码,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;