光学基于matlab GUI干涉衍射仿真【含Matlab源码 1723期】.zip

%双缝干涉可改D或f和间距 %等倾干涉可改厚度和角度 %等厚干涉可改波长 %衍射的焦距是公共的 %距孔衍射是以长为基的正方孔 % HUAS 201611030117 function varargout = GS_YS(varargin) gui_Singleton = 1; gui_State = struct('gui_Name', mfilename, ... 'gui_Singleton', gui_Singleton, ... 'gui_OpeningFcn', @GS_YS_OpeningFcn, ... 'gui_OutputFcn', @GS_YS_OutputFcn, ... 'gui_LayoutFcn', [] , ... 'gui_Callback', []); if nargin && ischar(varargin{1}) gui_State.gui_Callback = str2func(varargin{1}); end if nargout [varargout{1:nargout}] = gui_mainfcn(gui_State, varargin{:}); else gui_mainfcn(gui_State, varargin{:}); end % End initialization code - DO NOT EDIT % --- Executes just before GS_YS is made visible. function GS_YS_OpeningFcn(hObject, eventdata, handles, varargin) handles.output = hObject; % Update handles structure guidata(hObject, handles); % --- Outputs from this function are returned to the command line. function varargout = GS_YS_OutputFcn(hObject, eventdata, handles) varargout{1} = handles.output; % --- Executes on button press in pushbutton1_ShF. function pushbutton1_ShF_Callback(hObject, eventdata, handles) lambda=632e-9; % d=0.002; %双缝间距d d=str2num(get(handles.ShF_jianju,'string')); D=str2num(get(handles.D_or_f,'string')); % ymax=0.035; xs=ymax; ny=101; ys=linspace(-ymax,ymax,ny); for i=1:ny L1=sqrt((ys(i)-d/2).^2+D/2); L2=sqrt((ys(i)+d/2).^2+D/2); phi=2*pi*(L1-L2)/lambda; % b(i,:)=4*cos(phi/2).^2; %每个元素的平方， b(i,:)=4*cos(phi/2)^2; %整体的平方 end figure(gcf); nclevels=255; br=(b/4.014)*nclevels; % subplot(1,2,1) % title('双缝干涉'); figure(1); image(xs,ys,br); colormap(gray(nclevels)); % subplot(1,2,2) % plot(b(:),ys) % --- Executes on button press in pushbutton2_DQ. function pushbutton2_DQ_Callback(hObject, eventdata, handles) f=0.2; lambda=500*10^(-9); % d=input('薄膜厚度='); d=0.00004;%薄膜厚度e-5量级 % theta=input('倾斜角度='); theta=0.8;%倾斜角度0.1-1 d=str2num(get(handles.DQ_hd,'string'));%获取 theta=str2num(get(handles.DQ_jd,'string'));%倾斜角度0.1-1 rMax=f*tan(theta/2); N=501; for i=1:N x(i)=(i-1)*2*rMax/(N-1)-rMax; for j=1:N y(i)=(j-1)*2*rMax/(N-1)-rMax; r(i,j)=sqrt(x(i)^2+y(i)^2); delta(i,j)=2*d/sqrt(1+r(i,j)^2/f^2); Phi(i,j)=2*pi*delta(i,j)/lambda; B(i,j)=4*cos(Phi(i,j)/2)^2; end end NCLevels=255; Br=B/4*NCLevels; % title('等倾干涉'); figure(2); image(x,y,Br); colormap(gray(NCLevels)); axis square; % --- Executes on button press in pushbutton3_DH. function pushbutton3_DH_Callback(hObject, eventdata, handles) % lamda = 0.000632; lamda = str2num(get(handles.lamada,'string')); [X,Y] = meshgrid(-2:0.01:2); for R=3200:-500:200 I=cos(pi/lamda*(X.^2+Y.^2)/(R)+pi/2); % title('等厚干涉'); figure(3); imshow(I) %将光强显示为灰度值 % pause(0.8) end function DQ_hd_Callback(hObject, eventdata, handles) % --- Executes during object creation, after setting all properties. function DQ_hd_CreateFcn(hObject, eventdata, handles) if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor')) set(hObject,'BackgroundColor','white'); end function DQ_jd_Callback(hObject, eventdata, handles) % --- Executes during object creation, after setting all properties. function DQ_jd_CreateFcn(hObject, eventdata, handles) if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor')) set(hObject,'BackgroundColor','white'); end function ShF_jianju_Callback(hObject, eventdata, handles) % --- Executes during object creation, after setting all properties. function ShF_jianju_CreateFcn(hObject, eventdata, handles) if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor')) set(hObject,'BackgroundColor','white'); end function D_or_f_Callback(hObject, eventdata, handles) % --- Executes during object creation, after setting all properties. function D_or_f_CreateFcn(hObject, eventdata, handles) if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor')) set(hObject,'BackgroundColor','white'); end function lamada_Callback(hObject, eventdata, handles) % --- Executes during object creation, after setting all properties. function lamada_CreateFcn(hObject, eventdata, handles) if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor')) set(hObject,'BackgroundColor','white'); end function slit_wid_Callback(hObject, eventdata, handles) % --- Executes during object creation, after setting all properties. function slit_wid_CreateFcn(hObject, eventdata, handles) if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor')) set(hObject,'BackgroundColor','white'); end function rec_len_Callback(hObject, eventdata, handles) % --- Executes during object creation, after setting all properties. function rec_len_CreateFcn(hObject, eventdata, handles) if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor')) set(hObject,'BackgroundColor','white'); end function rec_wid_Callback(hObject, eventdata, handles) % --- Executes during object creation, after setting all properties. function rec_wid_CreateFcn(hObject, eventdata, handles) if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor')) set(hObject,'BackgroundColor','white'); end function focus_Callback(hObject, eventdata, handles) % --- Executes during object creation, after setting all properties. function focus_CreateFcn(hObject, eventdata, handles) if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor')) set(hObject,'BackgroundColor','white'); end % --- Executes on button press in pushbutton7_yuankong. function pushbutton7_yuankong_Callback(hObject, eventdata, handles) lambda=632e-9; f=str2num(get(handles.focus,'string')); a=2e-3; m=200; ym=2000*lambda*f; ys=linspace(-ym,ym,m); xs=ys; n=255; for i=1:m r=xs(i)^2+ys.^2; sinth=sqrt(r./(r+f^2)); x=2*pi*a*sinth./lambda; hh=(2*besselj(1,x)).^2./x.^2; %b(:,i)=(hh).*5000; b(:,i)=(hh).*10000; end figure(4); % imshow(b)%失去了颜色 % figure(2); % subplot(1,2,1); image(xs,ys,b) % subplot(1,2,2); % plot(hh,x) % colormap(gray(n)); % --- Executes on button press in pushbutton8_jukong. function pushbutton8_jukong_Callback(hObject, eventdata, handles) len= str2num(get(handles.rec_len,'string'));%获取距孔长 wid= str2num(get(handles.rec_wid,'string')); foc= str2num(get(handles.focus,'string')); lamr=650e-9; %lamg=520e-9; %lamb=450e-9; f=foc; %定义波长以及透镜长度 xmax=len; ymax=wid; % xmax=0.025e-4; % ymax=0.025e-4; sinix=0; siniy=0; %定义矩形孔尺寸以及x，y轴偏移量 cx=0.3; cy=0.3; %观察屏尺寸 N=10e2; x=linspace(-xmax,xmax,N); y=linspace(-ymax,ymax,N); %如果采用积分法需要使用的光源坐标坐标 x1=linspace(-cx,cx,N); y1=linspace(-cy,cy,N); %观察屏坐标 numx=length(x); numy=length(y); I=zeros(numx,numy,3); %定义三色矩阵 for t=1:numx for tt=1:numy alphar=pi*2*xmax/lamr*(x1(t)/f-sinix); % alphag=pi*2*xmax/lamg*(x1(t)/f-sinix); % alphab=pi*2*xmax/lamb*(x1(t)/f-sinix); betar=pi*2*xmax/lamr*(x1(tt)/f-siniy); % betag=pi*2*xmax/lamg*(x1(tt)/f-siniy); % betab=pi*2*xmax/lamb*(x1(tt)/f-siniy); I(t,tt,1)=(sin(alphar)/alphar)^2*(sin(betar)/betar)^2; Ix(t)=(sin(alphar)/alphar)^2; Iy(tt)=(sin(betar)/betar)^2; % I(t,tt,2)=(sin(alphag)/alphag)^2*(sin(betag)/betag)^2; % I(t,tt,3)=(sin(alphab)/alphab)^2*(sin(betab)/betab)^2; end end Ix=Ix/max(Ix); Iy=Iy/max(Iy); I=I/max(max(max(I))); % subplot(2,2,1); % imshow(255*I) % subplot(2,2,2); % plot(Ix,x) % subplot(2,2,3); % plot(y,Iy) figure(5); subplot(2,2,1); imshow(255*I) subplot(2,2,2); plo