《Matlab在信息光学中应用.ppt》由会员分享，可在线阅读，更多相关《Matlab在信息光学中应用.ppt（52页珍藏版）》请在三一办公上搜索。
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;endfigure;imagesc(p1)title(单缝衍射模拟图);colormap(lgray
2、),多缝衍射,多缝衍射代码,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;endfigure;imagesc(p1)title(多缝衍射模拟图);colormap(lgray);,平面光栅衍
3、射,平面衍射计算机模拟代码,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;endfigure;imagesc(p1)title(多缝衍射模拟图);col
4、ormap(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;endip=imshow(A,),矩孔衍射,矩孔衍射代码,x=zeros(30);x(10:20,10:20)=1imshow(x,-1 5)x=fft2(x,256,256);x=fftshift(x);imshow(log(abs(x),-1 5);,巴比涅互
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(移至中心频谱),多缝衍射计算机模拟,旋转对频谱的影响代码,x1,map=imread(j1.bmp);figure;subplot(2,3,1);imshow(x1);x1=fft2
6、(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;u=fft2(u,256,256);u=fftshift(u);subplot(2,4,k);imshow(log(abs(u),-1 1);title(夫瑯禾费衍射);end;,空域位移对夫瑯禾菲衍射的影响,
