光学设计 MATLAB 南邮

function varargout = A4(varargin) % A4 M-file for A4.fig % A4, by itself, creates a new A4 or raises the existing % singleton*. % % H = A4 returns the handle to a new A4 or the handle to % the existing singleton*. % % A4('CALLBACK',hObject,eventData,handles,...) calls the local % function named CALLBACK in A4.M with the given input arguments. % % A4('Property','Value',...) creates a new A4 or raises the % existing singleton*. Starting from the left, property value pairs are % applied to the GUI before A4_OpeningFunction gets called. An % unrecognized property name or invalid value makes property application % stop. All inputs are passed to A4_OpeningFcn via varargin. % % *See GUI Options on GUIDE's Tools menu. Choose "GUI allows only one % instance to run (singleton)". % % See also: GUIDE, GUIDATA, GUIHANDLES % Edit the above text to modify the response to help A4 % Last Modified by GUIDE v2.5 12-Oct-2010 12:32:29 % Begin initialization code - DO NOT EDIT gui_Singleton = 1; gui_State = struct('gui_Name', mfilename, ... 'gui_Singleton', gui_Singleton, ... 'gui_OpeningFcn', @A4_OpeningFcn, ... 'gui_OutputFcn', @A4_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 A4 is made visible. function A4_OpeningFcn(hObject, eventdata, handles, varargin) % This function has no output args, see OutputFcn. % hObject handle to figure % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) % varargin command line arguments to A4 (see VARARGIN) % Choose default command line output for A4 handles.output = hObject; % Update handles structure guidata(hObject, handles); % UIWAIT makes A4 wait for user response (see UIRESUME) % uiwait(handles.figure1); % --- Outputs from this function are returned to the command line. function varargout = A4_OutputFcn(hObject, eventdata, handles) % varargout cell array for returning output args (see VARARGOUT); % hObject handle to figure % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) % Get default command line output from handles structure varargout{1} = handles.output; function N1_Callback(hObject, eventdata, handles) % hObject handle to N1 (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) % Hints: get(hObject,'String') returns contents of N1 as text % str2double(get(hObject,'String')) returns contents of N1 as a double % --- Executes during object creation, after setting all properties. function N1_CreateFcn(hObject, eventdata, handles) % hObject handle to N1 (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles empty - handles not created until after all CreateFcns called % Hint: edit controls usually have a white background on Windows. % See ISPC and COMPUTER. if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor')) set(hObject,'BackgroundColor','white'); end function N2_Callback(hObject, eventdata, handles) % hObject handle to N2 (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) % Hints: get(hObject,'String') returns contents of N2 as text % str2double(get(hObject,'String')) returns contents of N2 as a double % --- Executes during object creation, after setting all properties. function N2_CreateFcn(hObject, eventdata, handles) % hObject handle to N2 (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles empty - handles not created until after all CreateFcns called % Hint: edit controls usually have a white background on Windows. % See ISPC and COMPUTER. if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor')) set(hObject,'BackgroundColor','white'); end function a_Callback(hObject, eventdata, handles) % hObject handle to a (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) % Hints: get(hObject,'String') returns contents of a as text % str2double(get(hObject,'String')) returns contents of a as a double % --- Executes during object creation, after setting all properties. function a_CreateFcn(hObject, eventdata, handles) % hObject handle to a (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles empty - handles not created until after all CreateFcns called % Hint: edit controls usually have a white background on Windows. % See ISPC and COMPUTER. if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor')) set(hObject,'BackgroundColor','white'); end % --- Executes on button press in pushbutton1. function pushbutton1_Callback(hObject, eventdata, handles) close; % --- Executes on button press in pushbutton2. function pushbutton2_Callback(hObject, eventdata, handles) n1=str2double(get(handles.N1,'String')); n2=str2double(get(handles.N2,'String')); a1=str2double(get(handles.a,'String'))*pi/180; b1=asin(n1*sin(a1)/n2); if(n1<=n2) rs1=(n1*cos(a1)-n2*cos(b1))/(n1*cos(a1)+n2*cos(b1)); set(handles.rs,'String',rs1); rp1=(n2*cos(a1)-n1*cos(b1))/(n2*cos(a1)+n1*cos(b1)); set(handles.rp,'String',rp1); ts1=2*n1*cos(a1)/(n1*cos(a1)+n2*cos(b1)); set(handles.ts,'String',ts1); tp1=2*n1*cos(a1)/(n2*cos(a1)+n1*cos(b1)); set(handles.tp,'String',tp1); R1=0.5*(rs1^2+rp1^2); set(handles.R,'String',R1) a=0:0.01:pi/2; b=asin(n1*sin(a)/n2); rs=(n1*cos(a)-n2*cos(b))./(n1*cos(a)+n2*cos(b)); rp=(n2*cos(a)-n1*cos(b))./(n2*cos(a)+n1*cos(b)); ts=2*n1*cos(a)./(n1*cos(a)+n2*cos(b)); tp=2*n1*cos(a)./(n2*cos(a)+n1*cos(b)); R=0.5*(rs.^2+rp.^2); axes(handles.axes1) hold off plot(a,rs,'-');grid on; hold on plot(a,rp,'--'); plot(a,ts,'.'); plot(a,tp,'-.'); legend rs rp ts tp; xlabel('角度（单位：弧度）','Color', 'm'); ylabel('', 'Color', 'm'); title('rs rp ts tp 随角度变化关系', 'Color','b'); axes(handles.axes5) plot(a,R);grid on; xlabel('角度（单位：弧度）','Color', 'm'); ylabel('反射率R', 'Color', 'm'); title('反射率R随角度变化关系', 'Color','b'); legend R end if(n1>n2) if(a1<asin(n2/n1)) rs1=(n1*cos(a1)-n2*cos(b1))/(n1*cos(a1)+n2*cos(b1)); rp1=(n2*cos(a1)-n1*cos(b1))/(n2*cos(a1)+n1*cos(b1)); else rs1=1; rp1=1; end set(handles.rs,'String',rs1); set(handles.rp,'String',rp1); set(handles.ts,'String',''); set(handles.tp,'String',''); R1=0.5*(rs1^2+rp1^2); set(handles.R,'String',R1) a=0:0.01:(asin(n2/n1)+0.001); b=asin(n1*sin(a)/n2); rs=(n1*cos(a)-n2*cos(b))./(n1*cos(a)+n2*cos(b)); rp=(n2*cos(a)-n1*cos(b))./(n2*cos(a)+n1*cos(b)); R=0.5*(rs.^2+rp.^2); axes(handles.axes1) hold off plot(a,rs,'-');grid on; hold on plot(a,rp,'--'); axis([0,pi/2,-1,1]) legend rs rp ; xlabel('角度（单位：弧度）','Color', 'm'); ylabel('', 'Color', 'm'); title('光密到光疏 rs rp 随角度变化关系', 'Color','b'); axes(handles.axes5) plot(a,R);grid on; axis([0,pi/2,-0.01,1]) xlabel('角度（单位：弧度）','Color', 'm'); ylabel('反射率R', 'Color', 'm'); title('反射率R随角度变化关系', 'Color','b'); legend R end funct