dblquadgui.m

数值分析程序中的一个例子

function varargout = dblquadGUI(varargin)
%DBLQUADGUI M-file for dblquadGUI.fig
%      DBLQUADGUI, by itself, creates a new DBLQUADGUI or raises the existing
%      singleton*.
%
%      H = DBLQUADGUI returns the handle to a new DBLQUADGUI or the handle to
%      the existing singleton*.
%
%      DBLQUADGUI('Property','Value',...) creates a new DBLQUADGUI using the
%      given property value pairs. Unrecognized properties are passed via
%      varargin to dblquadGUI_OpeningFcn.  This calling syntax produces a
%      warning when there is an existing singleton*.
%
%      DBLQUADGUI('CALLBACK') and DBLQUADGUI('CALLBACK',hObject,...) call the
%      local function named CALLBACK in DBLQUADGUI.M with the given input
%      arguments.
%
%      *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 dblquadGUI

% Last Modified by GUIDE v2.5 06-Nov-2005 16:34:55

% Begin initialization code - DO NOT EDIT
gui_Singleton = 1;
gui_State = struct('gui_Name',       mfilename, ...
                   'gui_Singleton',  gui_Singleton, ...
                   'gui_OpeningFcn', @dblquadGUI_OpeningFcn, ...
                   'gui_OutputFcn',  @dblquadGUI_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 dblquadGUI is made visible.
function dblquadGUI_OpeningFcn(hObject, eventdata, handles, varargin)
set(handles.fcn,'string','tan(x.^2+y.^2)');
set(handles.a,'string','0');
set(handles.b,'string','pi/3');
set(handles.c,'string','0');
set(handles.d,'string','pi/6');
set(handles.eps,'string','0.5*1e-9');
% 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   unrecognized PropertyName/PropertyValue pairs from the
%            command line (see VARARGIN)

% Choose default command line output for dblquadGUI
handles.output = hObject;

% Update handles structure
guidata(hObject, handles);

% UIWAIT makes dblquadGUI wait for user response (see UIRESUME)
% uiwait(handles.figure1);


% --- Outputs from this function are returned to the command line.
function varargout = dblquadGUI_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 fcn_Callback(hObject, eventdata, handles)
fcn=get(handles.fcn,'string');
set(handles.fcn,'string',fcn);
% hObject    handle to fcn (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 fcn as text
%        str2double(get(hObject,'String')) returns contents of fcn as a double


% --- Executes during object creation, after setting all properties.
function fcn_CreateFcn(hObject, eventdata, handles)
% hObject    handle to fcn (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 INV_fcn_Callback(hObject, eventdata, handles)
% hObject    handle to INV_fcn (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 INV_fcn as text
%        str2double(get(hObject,'String')) returns contents of INV_fcn as a double


% --- Executes during object creation, after setting all properties.
function INV_fcn_CreateFcn(hObject, eventdata, handles)
% hObject    handle to INV_fcn (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 Cal.
function Cal_Callback(hObject, eventdata, handles)
format long
fcn=get(handles.fcn,'string');
a=str2num(get(handles.a,'string'));
b=str2num(get(handles.b,'string'));
c=str2num(get(handles.c,'string'));
d=str2num(get(handles.d,'string'));
eps=str2num(get(handles.eps,'string'));
tic;
I=dblquad(fcn,a,b,c,d,'quad8');
toc;
handles.I=I;
set(handles.INV_fcn,'string',num2str(I,10));
set(handles.CT,'string',num2str(toc));
guidata(hObject,handles);

% hObject    handle to Cal (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)



% --- 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



function c_Callback(hObject, eventdata, handles)
 c= get(gcbo,'String');
 set(handles.c,'string',c); 
% hObject    handle to c (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 c as text
%        str2double(get(hObject,'String')) returns contents of c as a double


% --- Executes during object creation, after setting all properties.
function c_CreateFcn(hObject, eventdata, handles)
% hObject    handle to c (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 during object creation, after setting all properties.
function b_CreateFcn(hObject, eventdata, handles)
% hObject    handle to b (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 during object creation, after setting all properties.
function d_CreateFcn(hObject, eventdata, handles)
% hObject    handle to d (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 during object creation, after setting all properties.
function CT_CreateFcn(hObject, eventdata, handles)
% hObject    handle to CT (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 Closef.
function Closef_Callback(hObject, eventdata, handles)
close;
% hObject    handle to Closef (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)

function [I2 i]=MultipleIntegral(fcn,a,b,c,d,m0,n0,eps);
i=0;
I1=0;
I2=0;
fxy=inline(fcn,'x','y');
while(i==0|abs(I2-I1)>eps)
    m=m0*2^i;
    n=n0*2^i;
    h=(b-a)/m;
    k=(d-c)/n;
    x=a:h:b;
    y=c:k:d;
    X=zeros(m,1);
    Y=zeros(n,1);
    for p=1:m
        X(p,1)=(x(p)+x(p+1))/2;
    end
    for q=1:n
        Y(q,1)=(y(q)+y(q+1))/2;
    end 
    S=zeros(m,n);
    for p=1:m
        for q=1:n
             S(p,q)=h*k/4*(fxy(x(p),y(q))+fxy(x(p),y(q+1))...
             +fxy(x(p+1),y(q))+fxy(x(p+1),y(q+1)));
        end
    end
    TMN=sum(sum(S));
    S=zeros(m,1);
    for p=1:m
        S(p,1)=fxy(X(p,1),y(1))+fxy(X(p,1),y(n+1));
    end
    T2MN=1/4*TMN+h*k/8*sum(sum(S));
    S=zeros(n,1);
    for q=1:n
         S(q,1)=fxy(x(1),Y(q,1))+fxy(x(m+1),Y(q,1));
    end
    T2MN=T2MN++h*k/8*sum(sum(S));
    S=zeros(m,n-1);
    for p=1:m
        for q=2:n
            S(p,q-1)=fxy(X(p),y(q));
        end
    end
    T2MN=T2MN+h*k/4*sum(sum(S));
     S=zeros(m-1,n);
    for p=2:m
        for q=1:n
            S(p-1,q)=fxy(x(p),Y(q));
        end
    end
    T2MN=T2MN+h*k/4*sum(sum(S));
    S=zeros(m,n);
    for p=1:m
        for q=1:n
            S(p,q)=fxy(X(p),Y(q));
        end
    end
    T2MN=T2MN+h*k/4*sum(sum(S));
    i=i+1;
    I1=I2;
    I2=4/3*T2MN-TMN/3;
    clear x y X Y S;
end
 
function a_Callback(hObject, eventdata, handles)
 a= get(gcbo,'String');
 set(handles.a,'string',a); 
% 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



function b_Callback(hObject, eventdata, handles)
 b= get(gcbo,'String');
 set(handles.b,'string',b); 
% hObject    handle to b (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 b as text
%        str2double(get(hObject,'String')) returns contents of b as a double



function d_Callback(hObject, eventdata, handles)
 d= get(gcbo,'String');
 set(handles.d,'string',d); 
% hObject    handle to d (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 d as text
%        str2double(get(hObject,'String')) returns contents of d as a double


function eps_Callback(hObject, eventdata, handles)
 eps= get(gcbo,'String');
 set(handles.eps,'string',eps); 
% hObject    handle to eps (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 eps as text
%        str2double(get(hObject,'String')) returns contents of eps as a double


% --- Executes during object creation, after setting all properties.
function eps_CreateFcn(hObject, eventdata, handles)
% hObject    handle to eps (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 CT_Callback(hObject, eventdata, handles)
% hObject    handle to CT (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 CT as text
%        str2double(get(hObject,'String')) returns contents of CT as a double