Matlab數(shù)字衍射光學實驗(一)

1、phi = 0:pi/10:2*pi; % x的范圍從0到2z = exp(i * phi);u = real(z);%取實部v = image(z);%取虛部figure(1);%畫圖plot(phi,u);%畫實部的變化曲線figure(2);plot(phi,v)%畫虛部的變化曲線clear;%清除內(nèi)存close all;%關閉所有窗口x=linspace(0,2*pi,100);%設置橫坐標范圍及網(wǎng)格y1=sqrt(x);%開根號figure(1);plot(x,y1,'r')%做出y1的曲線,'r'代表圖線為紅色clear;%清除內(nèi)存close all

2、;%關閉所有窗口x=linspace(-2*pi,2*pi,100);%100個數(shù)橫坐標從-2pi到2piy1=besselj(1,x);%第一類貝塞爾函數(shù),一階y2=besselj(0,x);%第一類貝塞爾函數(shù),0階figure(1);plot(x,y1,'r');figure(2);plot(x,y2,'b');clear;close all;x,y=meshgrid(-2:.05:2,-2:.05:2);%設置二維網(wǎng)格z=sinc(x).*sinc(y);figure(1);imshow(z);%二維灰度圖figure(2);imagesc(z);%二維彩

3、圖figure(3);mesh(z);%三維圖像clear;close all;a=5;x=linspace(-5,5,100);y1=exp(-pi.*(x/a).2);%第一類貝塞爾函數(shù),一階b=2;y2=exp(-pi.*(x/b).2);%第一類貝塞爾函數(shù),一階figure(1);plot(x,y1,'r');figure(2);plot(x,y2,'b');clear;close all;a=1;b=1;x,y=meshgrid(-2:.05:2,-2:.05:2);%設置二維網(wǎng)格z=exp(-pi.*(x/a).2+(y/b).2);%高斯函數(shù)fig

4、ure(1);imagesc(z);figure(2);mesh(z);clear;close all;M=100;%x方向抽樣點數(shù)N=100;%y方向抽樣點數(shù)w=10;%縫的寬度d=10;%縫偏離中心的位置obj1=zeros(M,N);%賦值為0的一個二維矩陣,大小為M*Nobj1(:,(N/2-w/2):N/2+w/2)=1.0;%設置縫在y向中心,且縫的寬度為wfigure(1);imshow(obj1);obj2=zeros(M,N);obj2(:,(N/2-d)-w/2:(N/2-d)+w/2)=1.0;%縫偏離中心位置d個像素figure(2);imshow(obj2); cle

5、ar;close all;M=100;%x方向抽樣點數(shù)N=100;%y方向抽樣點數(shù)w=10;%縫的寬度d=10;%縫偏離中心的位置obj1=zeros(M,N);%賦值為0的一個二維矩陣,大小為M*Nobj1(:,(N/2+d)-w/2:(N/2+d)+w/2)=1.0;%設置縫在y向中心,且縫的寬度為wobj2=zeros(M,N);obj2(:,(N/2-d)-w/2:(N/2-d)+w/2)=1.0;%縫偏離中心位置d個像素obj3=obj1+obj2;figure(1);imshow(obj3);clear;close all;M=100;N=100;w=10;%縫寬obj=zeros(M,N);obj(M/2-w/2):(M/2+w/2),(N/2-w/2):(N/2+w