平行平面腔自再现模FoxLi数值迭代解法及MATLAB实现

axes(handles.axes3>

LOGO=imread('hit.jpg'>。 imshow(LOGO>

% --- Outputs from this function are returned to the command line.v4bdyGious function varargout = Laser_GUI_OutputFcn(hObject, eventdata, handles> J0bm4qMpJ9 % varargout cell array for returning output args (see VARARGOUT>。XVauA9grYP % hObject handle to figure

% eventdata reserved - to be defined in a future version of MATLABbR9C6TJscw % handles structure with handles and user data (see GUIDATA>pN9LBDdtrd

% Get default command line output from handles structureDJ8T7nHuGT varargout{1} = handles.output。

function wavelength_Callback(hObject, eventdata, handles>QF81D7bvUA % hObject handle to wavelength (see GCBO>

% eventdata reserved - to be defined in a future version of MATLAB4B7a9QFw9h % handles structure with handles and user data (see GUIDATA>ix6iFA8xoX

% Hints: get(hObject,'String'> returns contents of wavelength as textwt6qbkCyDE % str2double(get(hObject,'String'>> returns contents of wavelength as a doubleKp5zH46zRk

% --- Executes during object creation, after setting all properties.Yl4HdOAA61 function wavelength_CreateFcn(hObject, eventdata, handles>ch4PJx4BlI % hObject handle to wavelength (see GCBO>

% eventdata reserved - to be defined in a future version of MATLABqd3YfhxCzo % handles empty - handles not created until after all CreateFcns calledE836L11DO5

% Hint: edit controls usually have a white background on Windows.S42ehLvE3M % See ISPC and COMPUTER.

if ispc && isequal(get(hObject,'BackgroundColor'>, get(0,'defaultUicontrolBackgroundColor'>>501nNvZFis set(hObject,'BackgroundColor','white'>。 end

function Length_Callback(hObject, eventdata, handles>jW1viftGw9 % hObject handle to Length (see GCBO>

% eventdata reserved - to be defined in a future version of MATLABxS0DOYWHLP

% handles structure with handles and user data (see GUIDATA>LOZMkIqI0w

% Hints: get(hObject,'String'> returns contents of Length as textZKZUQsUJed % str2double(get(hObject,'String'>> returns contents of Length as a doubledGY2mcoKtT

% --- Executes during object creation, after setting all properties.rCYbSWRLIA function Length_CreateFcn(hObject, eventdata, handles>FyXjoFlMWh % hObject handle to Length (see GCBO>

% eventdata reserved - to be defined in a future version of MATLABTuWrUpPObX % handles empty - handles not created until after all CreateFcns called7qWAq9jPqE

% Hint: edit controls usually have a white background on Windows.llVIWTNQFk % See ISPC and COMPUTER. if ispc && isequal(get(hObject,'BackgroundColor'>, get(0,'defaultUicontrolBackgroundColor'>>yhUQsDgRT1 set(hObject,'BackgroundColor','white'>。 end

function area_Callback(hObject, eventdata, handles>MdUZYnKS8I % hObject handle to area (see GCBO>

% eventdata reserved - to be defined in a future version of MATLAB09T7t6eTno % handles structure with handles and user data (see GUIDATA>e5TfZQIUB5

% Hints: get(hObject,'String'> returns contents of area as texts1SovAcVQM % str2double(get(hObject,'String'>> returns contents of area as a doubleGXRw1kFW5s

% --- Executes during object creation, after setting all properties.UTREx49Xj9 function area_CreateFcn(hObject, eventdata, handles>8PQN3NDYyP % hObject handle to area (see GCBO>

% eventdata reserved - to be defined in a future version of MATLABmLPVzx7ZNw % handles empty - handles not created until after all CreateFcns calledAHP35hB02d

% Hint: edit controls usually have a white background on Windows.NDOcB141gT % See ISPC and COMPUTER.

if ispc && isequal(get(hObject,'BackgroundColor'>, get(0,'defaultUicontrolBackgroundColor'>>1zOk7Ly2vA set(hObject,'BackgroundColor','white'>。 end

function times_Callback(hObject, eventdata, handles>fuNsDv23Kh % hObject handle to times (see GCBO>

% eventdata reserved - to be defined in a future version of MATLABtqMB9ew4YX % handles structure with handles and user data (see GUIDATA>HmMJFY05dE

% Hints: get(hObject,'String'> returns contents of times as textViLRaIt6sk % str2double(get(hObject,'String'>> returns contents of times as a double9eK0GsX7H1

% --- Executes during object creation, after setting all properties.naK8ccr8VI function times_CreateFcn(hObject, eventdata, handles>B6JgIVV9ao % hObject handle to times (see GCBO>

% eventdata reserved - to be defined in a future version of MATLABP2IpeFpap5 % handles empty - handles not created until after all CreateFcns called3YIxKpScDM

% Hint: edit controls usually have a white background on Windows.gUHFg9mdSs % See ISPC and COMPUTER. if ispc && isequal(get(hObject,'BackgroundColor'>, get(0,'defaultUicontrolBackgroundColor'>>uQHOMTQe79 set(hObject,'BackgroundColor','white'>。 end

% --- Executes on button press in pushbutton2.

function pushbutton2_Callback(hObject, eventdata, handles>IMGWiDkflP % hObject handle to pushbutton2 (see GCBO>

% eventdata reserved - to be defined in a future version of MATLABWHF4OmOgAw % handles structure with handles and user data (see GUIDATA>aDFdk6hhPd

global steps L k a

lamda=str2num(get(handles.wavelength,'String'>>。 L=str2num(get(handles.Length,'String'>>。 a=str2num(get(handles.area,'String'>>。 N=str2num(get(handles.times,'String'>>。 k=2*pi/lamda。 steps=500。

x=linspace(-a,a,steps>。 u_=ones(1,steps>。 for m=1:N

for mm=1:steps

u0(mm>=QU(x(mm>,u_>。 end

axes(handles.axes1>

plot(x,abs(u0>/abs(u0(steps/2>>>

axis([-a a 0 1.3]>

angle_u0=angle(u0>/pi*180。

angle_u0=angle_u0-angle_u0(steps/2>。 axes(handles.axes2> plot(x,angle_u0>

u_=u0/max(abs(u0>>。 end

% --- Executes on button press in pushbutton3.

function pushbutton3_Callback(hObject, eventdata, handles>ozElQQLi4T % hObject handle to pushbutton3 (see GCBO>

% eventdata reserved - to be defined in a future version of MATLABCvDtmAfjiA % handles structure with handles and user data (see GUIDATA>QrDCRkJkxh global steps L k a

lamda=str2num(get(handles.wavelength,'String'>>。 L=str2num(get(handles.Length,'String'>>。 a=str2num(get(handles.area,'String'>>。 N=str2num(get(handles.times,'String'>>。 k=2*pi/lamda。 steps=500。

x=linspace(-a,a,steps>。

[X,Y]=meshgrid(linspace(-a,a,50>,linspace(-a,a,50>>。4nCKn3dlMX u_=ones(1,steps>。 for m=1:N

for mm=1:steps

u0(mm>=QU(x(mm>,u_>。 end

abs_u0=abs(u0>/abs(u0(steps/2>>。

Z_A=sqrt((interp1(x,abs_u0,X>>.^2+(interp1(x,abs_u0,Y>>.^2>。ijCSTNGm0E axes(handles.axes1> surf(X,Y,Z_A>

axis([-a a -a a 0 1.5]>

angle_u0=angle(u0>/pi*180。

angle_u0=angle_u0-angle_u0(steps/2>。

Z_P=sqrt((interp1(x,angle_u0,X>>.^2+(interp1(x,angle_u0,Y>>.^2>。vfB1pxanfk axes(handles.axes2> surf(X,Y,Z_P>

axis([-a a -a a -10 90]> u_=u0/max(abs(u0>>。 end

% --------------------------------------------------------------------JbA9VhEou1