bnbguicb.m

工具箱 (分枝定界法工具箱) matlab 中使用

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

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

switch action
  
case('init')
   mlock;
   data.id='BNBGUI20 BNB20';
   data.fun='';
   data.xtag={[]};
   data.x0=[];
   data.xstatus=[];
   data.xlb=[];
   data.xub=[];
   data.A=[];
   data.B=[];
   data.Aeq=[];
   data.Beq=[];
   data.nonlcon='';
   data.settings=[0];
   data.options=optimset('fmincon');
   data.options=optimset(data.options,'Display','iter');
   data.options.MaxSQPIter=1000;
   data.partag={[]};
   data.par=[];
   data.numx=0;
   data.numpar=0;
   data.selx=0;
   data.selpar=0;
   data.selopt=2;
   data.optliststring={'DerivativeCheck' 'Diagnostics' 'DiffMaxChange' 'DiffMinChange' ...
         'Display' 'GradConstr' 'GradObj' 'Hessian' 'LargeScale' ...
         'MaxFunEvals' 'MaxIter' 'MaxPCGIter' 'MaxSQPIter' 'PrecondBandWidth' 'TolCon' ...
         'TolFun' 'TolPCG' 'TolX'};
   popuponoff={'popupmenu' {'on' 'off'}}; 
   editposscalar={'edit' {'scalar' 'positive scalar'}};
   editposint={'edit' {'scalar' 'positive integer'}};
   edit0posintinf={'edit' {'scalar' '0 positive integer inf'}};
   data.optlistcells={popuponoff popuponoff editposscalar editposscalar ...
         {'popupmenu' {'off' 'iter' 'final'}} popuponoff popuponoff popuponoff popuponoff ...
         editposint editposint editposint editposint edit0posintinf editposscalar ...
         editposscalar editposscalar editposscalar};      
   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 2.0 for BNB 20';
   handles.main=[];
   handles.fun=[];
   handles.set=[];
   handles.err=[];
case('main')
   if ~ih(handles.main)
      handles.main=guimain;
   end;
   if data.numx>0
      onoffx(handles,'on');
      handle=findobj(handles.main,'Tag','xlist');
      set(handle,'String',num2str([1:data.numx]'));
      set(handle,'Value',data.selx);
      updatex(handles,data);
   else onoffx(handles,'off'); end;
   if data.numpar>0
      onoffpar(handles,'on');
      handle=findobj(handles.main,'Tag','parlist');
      set(handle,'String',num2str([1:data.numpar]'));
      set(handle,'Value',data.selpar);
      updatepar(handles,data);
   else onoffpar(handles,'off'); end;
   handle1=findobj(handles.main,'Tag','results');
   handle2=findobj(handles.main,'Tag','resultsslider');
   initwindowslider(handle1,handle2,data.results.msg);
case('save')
   [file,path]=uiputfile('bnb.mat','Select mat-file to save data.');
   if file~=0
      BNBdata=data;
      eval(['save ',fullfile(path,file),' BNBdata;'],'handles.err=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;'],'handles.err=uierror(''File not found.''); check=0;');
      if check==1
         if exist('BNBdata') & isstruct(BNBdata) ...
               & isfield(BNBdata,'id') & strcmp(BNBdata.id,'BNBGUI20 BNB20')
            mlock;
            data=BNBdata;
            bnbguicb('main');
            if ih(handles.fun), bnbguicb('function'); end;
            if ih(handles.set), bnbguicb('settings'); end;
         else handles.err=uierror('Data has not been updated.'); end;
      end;
   end;
case('X -> x0')
   if all(size(data.results.X)==size(data.x0)),
      data.x0=data.results.X;
      if data.numx>0, updatex(handles,data); end;
   end;
case('Z X t c fail -> base')
   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 base workspace.');
case('help')
   edit bnbhelp.txt;
case('copyright')
   msg={'E.C. Kuipers' 'e-mail:' '[e-mail]' 'FI-Lab' ...
         'Applied Physics' 'Rijksuniversiteit Groningen'};
   handles.err=uierror(msg);
   handle=findobj(handles.err,'Tag','error BNB GUI');
   set(handle,'Name','copyright BNB GUI');
   handle=findobj(handles.err,'Tag','error');
   set(handle,'String','copyright');
case('resultsslider')
   handle1=findobj(handles.main,'Tag','results');
   handle2=findobj(handles.main,'Tag','resultsslider');
   updatewindowslider(handle1,handle2);
case('xlist')
   data.selx=selectx(handles);
   updatex(handles,data);
case('xtag')
   handle=gcbo;
   [xtag,succes]=callbackmultiedit(handle,{'edit' {'string'}});
   data.xtag{data.selx}=xtag;
case('x0')
   handle=gcbo;
   [x0,succes]=callbackmultiedit(handle,{'edit' {'scalar' ''}});
   if succes, data.x0(data.selx)=x0; else set(gcbo,'String',data.x0(data.selx)); end; 
case('xlb')
   handle=gcbo;
   [xlb,succes]=callbackmultiedit(handle,{'edit' {'scalar' ''}});
   if succes, data.xlb(data.selx)=xlb; else set(gcbo,'String',data.xlb(data.selx)); end;
case('xub')
   handle=gcbo;
   [xub,succes]=callbackmultiedit(handle,{'edit' {'scalar' ''}});
   if succes, data.xub(data.selx)=xub; else set(gcbo,'String',data.xub(data.selx)); end;
case('continuous')
   data.xstatus(data.selx)=0;
   updatex(handles,data);
case('integer')
   data.xstatus(data.selx)=1;
   updatex(handles,data);
case('fixed')
   data.xstatus(data.selx)=2;
   updatex(handles,data);
case('parlist')
   data.selpar=selectpar(handles);
   updatepar(handles,data);
case('partag')
   handle=gcbo;
   [partag,succes]=callbackmultiedit(handle,{'edit' {'string'}});
   data.partag{data.selpar}=partag;
case('par')
   handle=gcbo;
   [par,succes]=callbackmultiedit(handle,{'edit' {'scalar' ''}});
   if succes, data.par(data.selpar)=par; else set(gcbo,'String',data.par(data.selpar)); end;
case('function')
   if ~ih(handles.fun)
      handles.fun=guifun;
   end;
   set(findobj(handles.fun,'Tag','fun'),'String',data.fun(1:min(12,size(data.fun,2))));
   set(findobj(handles.fun,'Tag','nonlcon'),'String',data.nonlcon(1:min(12,size(data.nonlcon,2))));
   set(findobj(handles.fun,'Tag','numx'),'String',num2str(data.numx));
   set(findobj(handles.fun,'Tag','numpar'),'String',num2str(data.numpar));
   updatefunclear(handles,data);
case('settings')
   if ~ih(handles.set)
      handles.set=guiset;
   end;
   set(findobj(handles.set,'Tag','set1'),'Value',data.settings(1));
   set(findobj(handles.set,'Tag','optlist'),'String',data.optliststring);
   set(findobj(handles.set,'Tag','optlist'),'Value',data.selopt);
   updateoptedit(handles,data);
case('optimize')
   par=num2cell(data.par);
   data.settings(2)=handles.main;
   [data.results.errmsg,data.results.Z,data.results.X,...
      data.results.t,data.results.c,data.results.fail]=...
      bnb20(data.fun,data.x0,data.xstatus,data.xlb,data.xub,...
      data.A,data.B,data.Aeq,data.Beq,data.nonlcon,...
      data.settings,data.options,par{:});
   if ~isempty(data.results.X)
      disp('X is:');
      disp(data.results.X);
      data.results.X(find(data.xstatus==1))=round(data.results.X(find(data.xstatus==1)));
   end;
   data.results.msg=updateresults(data);
   handle1=findobj(handles.main,'Tag','results');
   handle2=findobj(handles.main,'Tag','resultsslider');
   initwindowslider(handle1,handle2,data.results.msg);
case('hide main')
   if ih(handles.main), set(handles.main,'Visible','off'); end;
   if ih(handles.fun), set(handles.fun,'Visible','off'); end;
   if ih(handles.set), set(handles.set,'Visible','off'); end;
case('show main')
   if ih(handles.main), set(handles.main,'Visible','on'); end;
   if ih(handles.fun), set(handles.fun,'Visible','on'); end;
   if ih(handles.set), set(handles.set,'Visible','on'); end;
case('quit main')
   if ih(handles.fun), close(handles.fun); end;
   if ih(handles.set), close(handles.set); end;
   munlock;
   handles.main=[];
case('fun')
   file=uigetfile('*.m','Select m-file for fun(x).');
   if file~=0
      if checkextension(file,'.m')
         data.fun=file(1:(size(file,2)-2));
         set(findobj(handles.fun,'Tag','fun'),'String',data.fun(1:min(12,size(data.fun,2))));
         set(findobj(handles.fun,'Tag','clear fun'),'BackgroundColor',[1 0 0]);
      else handles.err=uierror({'No m-file.' 'fun has not been updated.'}); end;
   end; 
case('edit fun')
   if ~isempty(data.fun), eval(['edit ',data.fun]); else edit; end;
case('clear fun')
   data.fun='';
   set(findobj(handles.fun,'Tag','fun'),'String','');
   set(findobj(handles.fun,'Tag','clear fun'),'BackgroundColor',[0.7529 0.7529 0.7529]);
case('A')
   [succes,A]=loadmatrix('A');
   if succes==1
      data.A=A;
      set(findobj(handles.fun,'Tag','clear A'),'BackgroundColor',[1 0 0]);
   end;
case('clear A')
   data.A=[];
   set(findobj(handles.fun,'Tag','clear A'),'BackgroundColor',[0.7529 0.7529 0.7529]);
case('B')
   [succes,B]=loadmatrix('B');
   if succes==1
      data.B=B;
      set(findobj(handles.fun,'Tag','clear B'),'BackgroundColor',[1 0 0]);
   end;
case('clear B')
   data.B=[];
   set(findobj(handles.fun,'Tag','clear B'),'BackgroundColor',[0.7529 0.7529 0.7529]);   
case('Aeq')
   [succes,Aeq]=loadmatrix('Aeq');
   if succes==1
      data.Aeq=Aeq;
      set(findobj(handles.fun,'Tag','clear Aeq'),'BackgroundColor',[1 0 0]);
   end;
case('clear Aeq')
   data.Aeq=[];
   set(findobj(handles.fun,'Tag','clear Aeq'),'BackgroundColor',[0.7529 0.7529 0.7529]);   
case('Beq')
   [succes,Beq]=loadmatrix('Beq');
   if succes==1
      data.Beq=Beq;
      set(findobj(handles.fun,'Tag','clear Beq'),'BackgroundColor',[1 0 0]);
   end;
case('clear Beq')
   data.Beq=[];
   set(findobj(handles.fun,'Tag','clear Beq'),'BackgroundColor',[0.7529 0.7529 0.7529]);   
case('lincon')
   file=uigetfile('*.mat','Select mat-file for matrixes A,B, Aeq and Beq.');
   if file~=0
      if checkextension(file,'.mat')
         foundfile=1;
         eval(['load ',file,' A B Aeq Beq;'],'foundfile=0; handles.err=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;
            updatefunclear(handles,data);
            if size(errmsg,2)>0, handles.err=uierror(errmsg); end;
         end;
      else handles.err=uierror({'No mat-file.' 'lincon has not been updated.'}); end;
   end;
case('nonlcon')
   file=uigetfile('*.m','Select m-file for nonlcon(x).');
   if file~=0
      if checkextension(file,'.m')
         data.nonlcon=file(1:(size(file,2)-2));
         set(findobj(handles.fun,'Tag','nonlcon'),'String',data.nonlcon(1:min(12,size(data.nonlcon,2))));
         set(findobj(handles.fun,'Tag','clear nonlcon'),'BackgroundColor',[1 0 0]);
      else handles.err=uierror({'No m-file.' 'nonlincon has not been updated.'}); end;
   end;
case('edit nonlcon')
   if ~isempty(data.nonlcon), eval(['edit ',data.nonlcon]); else edit; end;
case('clear nonlcon')
   data.nonlcon='';
   set(findobj(handles.fun,'Tag','nonlcon'),'String','');
   set(findobj(handles.fun,'Tag','clear nonlcon'),'BackgroundColor',[0.7529 0.7529 0.7529]);   
case('numx')
   handle=gcbo;
   [numx,succes]=callbackmultiedit(handle,{'edit' {'scalar' '0 positive integer'}});
   if succes
      data.numx=numx;
      if numx==0
         data.x0=[];
         data.xstatus=[];
         data.xlb=[];
         data.xub=[];
         data.xtag={[]};
         onoffx(handles,'off');
      else
         lx=size(data.x0,1);
         if numx>lx
            data.x0=[data.x0;zeros(numx-lx,1)];
            data.xstatus=[data.xstatus;zeros(numx-lx,1)];
            data.xlb=[data.xlb;zeros(numx-lx,1)];
            data.xub=[data.xub;ones(numx-lx,1)];
            data.xtag{numx,1}=[]; 
         else
            data.x0=data.x0(1:numx);
            data.xstatus=data.xstatus(1:numx);
            data.xlb=data.xlb(1:numx);