多目標非線性規(guī)劃程序-資料下載頁
2025-07-22 05:46本頁面
【正文】 mbent = x。 separation = 0。 end。 % STEP 6 SELECTION if separation == 1 amp。 ~isempty(K) dzsep=1。 for i=1:size(K,1) dxsepc = abs(round(x(K(i)))x(K(i)))。 if dxsepc=errx | convflag==0 xsepc = x。 xsepc(K(i))=round(x(K(i)))。 dzsepc = abs(feval(fun,xsepc,varargin{:})z)。 if dzsepcdzsep dzsep=dzsepc。 ixsep=K(i)。 end。 end。 end。 % STEP 7 SEPARATION if xchoice(ixsep)==0 % XCHOICE==0 branch=1。 domain=[xlb(ixsep) xub(ixsep)]。 while branch==1 xboundary=(domain(1)+domain(2))/2。 if x(ixsep)xboundary domainA=[domain(1) floor(xboundary)]。 domainB=[floor(xboundary+1) domain(2)]。 else domainA=[floor(xboundary+1) domain(2)]。 domainB=[domain(1) floor(xboundary)]。 end。 stacksize=stacksize+1。 stackx0(:,stacksize)=x。 stackxlb(:,stacksize)=xlb。 stackxlb(ixsep,stacksize)=domainB(1)。 stackxub(:,stacksize)=xub。 stackxub(ixsep,stacksize)=domainB(2)。 if domainA(1)==domainA(2) stacksize=stacksize+1。 stackx0(:,stacksize)=x。 stackxlb(:,stacksize)=xlb。 stackxlb(ixsep,stacksize)=domainA(1)。 stackxub(:,stacksize)=xub。 stackxub(ixsep,stacksize)=domainA(2)。 branch=0。 else domain=domainA。 branch=1。 end。 end。 else % XCHOICE~=0 L=find(xchoice==xchoice(ixsep))。 M=intersect(K,L)。 [dummy,N]=sort(x(M))。 part1=M(N(1:floor(size(N)/2)))。 part2=M(N(floor(size(N)/2)+1:size(N)))。 stacksize=stacksize+1。 stackx0(:,stacksize)=x。 O = (1sum(stackx0(part1,stacksize)))/size(part1,1)。 stackx0(part1,stacksize)=stackx0(part1,stacksize)+O。 stackxlb(:,stacksize)=xlb。 stackxub(:,stacksize)=xub。 stackxub(part2,stacksize)=0。 stacksize=stacksize+1。 stackx0(:,stacksize)=x。 O = (1sum(stackx0(part2,stacksize)))/size(part2,1)。 stackx0(part2,stacksize)=stackx0(part2,stacksize)+O。 stackxlb(:,stacksize)=xlb。 stackxub(:,stacksize)=xub。 stackxub(part1,stacksize)=0。 if size(part2,1)==1, stackxlb(part2,stacksize)=1。 end。 end。 elseif separation==1 amp。 isempty(K) fail=fail+1。 end。end。 % STEP 8 OUTPUT t=toc。Z = z_incumbent。X = x_incumbent。errmsg=39。39。eval([39。warning 39。,currentwarningstate])。 function CON=BNBCON(x,A,B,Aeq,Beq,xlb,xub,nonlcon,varargin)。if isempty(A), CON1=[]。 else CON1 = max(A*xB,0)。 end。if isempty(Aeq), CON2=[]。 else CON2 = abs(Aeq*xBeq)。 end。CON3 = max(xlbx,0)。CON4 = max(xxub,0)。if isempty(nonlcon) CON5=[]。 CON6=[]。else [C Ceq]=feval(nonlcon,x,varargin{:})。 CON5 = max(C,0)。 CON6 = abs(Ceq)。end。CON = max([CON1。 CON2。 CON3。 CON4。 CON5。 CON6])。8 / 8
黨政相關相關推薦
鄂ICP備17016276號-1