bnbguicb.m

大部分都是数学建模的资料,大家要多多支持啊!

function BNBGUICB(action,file);
% BNBGUICB Callback function for BNBGUI.
% Do not run this file from the Matlab prompt.
persistent data;

if nargin<1, help bnbguicb; action=''; end;

switch action
  
case('init')
   data.id='BNBGUI10 BNB18';
   data.fun='';
   data.x0=[];
   data.xstatus=[];
   data.xlb=[];
   data.xub=[];
   data.A=[];
   data.B=[];
   data.Aeq=[];
   data.Beq=[];
   data.nonlcon='';
   data.settings=[0 0 0];
   data.maxSQPiter=1000;
   data.options1=optimset('fmincon');
   data.options2=optimset('fmincon');
   data.par=[];
   data.numx=0;
   data.numpar=0;
   data.selx=0;
   data.selpar=0;
   data.selopt=2;
   data.seloptpar=1;
   data.optliststring={'DerivativeCheck' 'Diagnostics' 'DiffMaxChange' ...
         'DiffMinChange' 'Display' 'GradConstr' 'GradObj' ...
         'LargeScale' 'MaxFunEvals' 'MaxIter' 'MaxPCGIter' 'PrecondBandWidth' ...
         'TolCon' 'TolFun' 'TolPCG' 'TolX'};
   data.results.errmsg='';
   data.results.Z=inf;
   data.results.X=[];
   data.results.t=0;
   data.results.c=0;
   data.results.fail=0;
   data.results.msg='BNB GUI for BNB18';
case('main')
   handle=findobj('Tag','main BNB GUI');
   if isempty(handle)
      mlock;
      guimain;
      handle=findobj('Tag','main BNB GUI');
      CloseRequestFcn='BNBGUICB(''quit main''); closereq;';
      set(handle,'CloseRequestFcn',CloseRequestFcn);
   end;
   if data.numx>0
      onoffx('on');
      handle=findobj('Tag','xlist');
      set(handle,'String',num2str([1:data.numx]'));
      set(handle,'Value',data.selx);
      updatex(data);
   else onoffx('off'); end;
   if data.numpar>0
      onoffpar('on');
      handle=findobj('Tag','parlist');
      set(handle,'String',num2str([1:data.numpar]'));
      set(handle,'Value',data.selpar);
      updatepar(data);
   else onoffpar('off'); end;
   showresults(data.results.msg);
case('resultsslider')
   showresults(data.results.msg)
case('X -> x0')
   if all(size(data.results.X)==size(data.x0)),
      data.x0=data.results.X;
      if data.numx>0, updatex(data); end;
   end;
case('Z X t c fail -> workspace')
   assignin('base','Z',data.results.Z);
   assignin('base','X',data.results.X);
   assignin('base','t',data.results.t);
   assignin('base','c',data.results.c);
   assignin('base','fail',data.results.fail);
   disp('The variables Z, X, t, c and fail are loaded in the workspace.');
case('clc workspace')
   clc;
case('copyright')
   errmsg={'E.C. Kuipers' 'e-mail:' '[e-mail]' 'FI-Lab' ...
         'Applied Physics' 'Rijksuniversiteit Groningen'};
   handle1=guierr;
   set(handle1,'UserData',errmsg);
   handle2=findobj('Tag','error');
   set(handle2,'String','Copyright:');
   bnbguicb('errorslider');
case('xlist')
   data.selx=selectx;
   updatex(data);
case('x0')
   data.selx=selectx;
   data.x0(data.selx)=strtonum(get(gcbo,'String'));
case('xlb')
   data.selx=selectx;
   data.xlb(data.selx)=strtonum(get(gcbo,'String'));
case('xub')
   data.selx=selectx;
   data.xub(data.selx)=strtonum(get(gcbo,'String'));
case('continuous')
   data.selx=selectx;
   data.xstatus(data.selx)=0;
   updatex(data);
case('integer')
   data.selx=selectx;
   data.xstatus(data.selx)=1;
   updatex(data);
case('fixed')
   data.selx=selectx;
   data.xstatus(data.selx)=2;
   updatex(data);
case('parlist')
   data.selpar=selectpar;
   updatepar(data);
case('par')
   data.selpar=selectpar;
   data.par(data.selpar)=strtonum(get(gcbo,'String'));
case('save')
   [file,path]=uiputfile('bnb.mat','Select mat-file to save data.');
   if file~=0
      BNBdata=data;
      eval(['save ',fullfile(path,file),' BNBdata;'],'uierror(''Data has not been saved.'')');
   end;
case('load')
   if ~exist('file'), file=uigetfile('bnb.mat','Select mat-file to load settings.'); end;
   if file~=0
      check=1;
      eval(['load ',file,' BNBdata;'],'uierror(''File not found.''); check=0;');
      if check==1
         if exist('BNBdata') & isstruct(BNBdata) ...
               & isfield(BNBdata,'id') & strcmp(BNBdata.id,'BNBGUI10 BNB18')
            data=BNBdata;
            bnbguicb('main');
            handle=findobj('Tag','function BNB GUI');
            if ~isempty(handle), bnbguicb('function'); end;
            handle=findobj('Tag','settings BNB GUI');
            if ~isempty(handle), bnbguicb('settings'); end;
         else uierror('Data has not been updated.'); end;
      end;
   end;
case('help')
   edit bnbhelp.txt;
case('function')
   handle=findobj('Tag','function BNB GUI');
   if isempty(handle), guifun; end;
   set(findobj('Tag','fun'),'String',data.fun);
   set(findobj('Tag','nonlcon'),'String',data.nonlcon);
   set(findobj('Tag','numx'),'String',num2str(data.numx));
   set(findobj('Tag','numpar'),'String',num2str(data.numpar));
case('settings')
   handle=findobj('Tag','settings BNB GUI');
   if isempty(handle), guiset; end;
   set(findobj('Tag','set1'),'Value',data.settings(1));
   set(findobj('Tag','set2'),'Value',data.settings(2));
   set(findobj('Tag','set3'),'Value',data.settings(3));
   set(findobj('Tag','maxSQPiter'),'String',data.maxSQPiter);
   set(findobj('Tag','opt1'),'Value',data.selopt==1);
   set(findobj('Tag','opt2'),'Value',data.selopt==2);
   set(findobj('Tag','optlist'),'String',data.optliststring);
   set(findobj('Tag','optlist'),'Value',data.seloptpar);
   updateopt(data);
case('optimize')
   par=num2cell(data.par);
   [data.results.errmsg,data.results.Z,data.results.X,...
      data.results.t,data.results.c,data.results.fail]=...
      bnb18(data.fun,data.x0,data.xstatus,data.xlb,data.xub,...
      data.A,data.B,data.Aeq,data.Beq,data.nonlcon,...
      data.settings,data.options1,data.options2,data.maxSQPiter,par{:});
   data.results.msg=updateresults(data.results);
   showresults(data.results.msg);
case('quit main')
   handle=findobj('Tag','function BNB GUI');
   if ~isempty(handle), close(handle); end;
   handle=findobj('Tag','settings BNB GUI');
   if ~isempty(handle), close(handle); end;
   munlock;
case('fun')
   fun=uigetfile('*.m','Select m-file for fun(x).');
   if fun~=0
      if all(fun(size(fun,2)-1:size(fun,2))=='.m')
         data.fun=fun(1:size(fun,2)-2);
         set(findobj('Tag','fun'),'String',data.fun(1:min(12,size(data.fun,2))));
      else uierror({'No m-file.' 'fun has not been updated.'}); end;
   end; 
case('funedit')
   if ~isempty(data.fun), eval(['edit ',data.fun]); else edit; end;
case('A')
   [succes,A]=loadmatrix('A');
   if succes==1, data.A=A; end;
case('B')
   [succes,B]=loadmatrix('B');
   if succes==1, data.B=B; end;
case('Aeq')
   [succes,Aeq]=loadmatrix('Aeq');
   if succes==1, data.Aeq=Aeq; end;
case('Beq')
   [succes,Beq]=loadmatrix('Beq');
   if succes==1, data.A=A; end;
case('lincon')
   file=uigetfile('*.mat','Select mat-file for matrixes A,B, Aeq and Beq.');
   if file~=0
      if all(file(size(file,2)-3:size(file,2))=='.mat')
         foundfile=1;
         eval(['load ',file,' A B Aeq Beq;'],'foundfile=0; uierror(''File not found.'');');
         if foundfile==1
            errmsg={};
            if exist('A'), data.A=A; 
            else errmsg{size(errmsg,2)+1}='A has not been updated'; end; 
            if exist('B'), data.B=B;
            else errmsg{size(errmsg,2)+1}='B has not been updated'; end;
            if exist('Aeq'), data.Aeq=Aeq;
            else errmsg{size(errmsg,2)+1}='Aeq has not been updated'; end;
            if exist('Beq'), data.Beq=Beq;
            else errmsg{size(errmsg,2)+1}='Beq has not been updated'; end;
            if size(errmsg,2)>0, uierror(errmsg); end;
         end;
      else uierror({'No mat-file.' 'lincon has not been updated.'}); end;
   end;
case('nonlcon')
   nonlcon=uigetfile('*.m','Select m-file for nonlcon(x).');
   if nonlcon~=0
      if all(nonlcon(size(nonlcon,2)-1:size(nonlcon,2))=='.m')
         data.nonlcon=nonlcon(1:size(nonlcon,2)-2);
         set(findobj('Tag','nonlcon'),'String',data.nonlcon(1:min(12,size(data.nonlcon,2))));
      else uierror({'No m-file.' 'nonlincon has not been updated.'}); end;
   end;
case('nonlconedit')