Matlab在信息光学中应用ppt课件

• imagesc(p1)
• title('多缝衍射模拟图');
• colormap(lgray);
.
平面光栅衍射
.
.
.
.
.
.
.
.
.
平面衍射计算机模拟代码
• clc;
• clear;
• n=input('平面衍射光栅数');
• d=-n*pi:0.0001*pi:pi*n;
Байду номын сангаас
• b=d/15;
.
空域位移对夫瑯禾菲衍射的影响
.
.
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)
.
.
.
.
.
.
圆孔衍射代码
• [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
.
矩孔衍射
.
.
矩孔衍射代码
• 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);
.
圆孔衍射
.
.
.
圆孔衍射代码
• 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,[])
• [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]);
.
多缝衍射
.
.
.
.
.
.
.
.
.
.
多缝衍射代码
• 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;
• 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('多缝缝衍射强度分布');
.
缝宽对夫瑯禾菲衍射的影响
.
缝宽对夫瑯禾菲衍射的影响代码
• 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;
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);i mshow(x1);x1=fft2(x1,256,256);x1=fftshift(x1);su bplot(2,3,4);imshow(log(abs(x1)),[-1 5]);