function r=checkint(x)%判断x(i)是不是整数了。是的话r(i)返回1，不是的话，返回0%输入参数：x   X向量%输出参数：r   R向量for i=1:length(x)    if(min(abs(x(i)-floor(x(i))),abs(x(i)-ceil(x(i))))<1e-3)        r(i)=1;    else        r(i)=0;    endend

function val=isrowinmat(arow,mat)%用来判断mat中是否包含与arow一样的向量%输入变量：arow    向量%         mat     矩阵%输出变量：val     1表示有，0表示没有val=0;rows=size(mat,1);for i=1:rows    temp=(mat(i,:)==arow);    if length(find(temp==0))==0        val=1;        return;    else        val=0;    end; end

function [x,fval,exitflag,output,lambda]=linprogdis(ifint,f,A,b,Aeq,beq,lb,ub,x0,options)% 用法%    [x,fval,exitflag,output,lambda]=lpint(ifint.f,A,b,Aeq,beq)%    [x,fval,exitflag,output,lambda]=lpint(ifint,f,A,b,Aeq,beq,lb)%    [x,fval,exitflag,output,lambda]=lpint(ifint,f,A,b,Aeq,beq,lb,ub)%    [x,fval,exitflag,output,lambda]=lpint(ifint,f,A,b,Aeq,beq,lb,ub,x0)%    [x,fval,exitflag,output,lambda]=lpint(ifint,f,A,b,Aeq,beq,lb,ub,x0,options)if nargin<10, options=[];  endif nargin<9,  x0=[];       endif nargin<8,  ub=inf*ones(size(f));      endif nargin<7,  lb=zeros(size(f));      end[x,fval,exitflag,output,lambda]=linprog(f,A,b,Aeq,beq,lb,ub,x0,options);if exitflag<=0        %表示线性规划没有最优解    return endv1=find(ifint==1);  %找到需要整数规划的变量的下标temp=x(v1);%如果不是要求整数规划的就可以返回了。if isempty(temp)    returnendv2=find(checkint(temp)==0);if isempty(v2)   %都是整数，得到最众解    returnendk=v1(v2(1));temp1=zeros(1,length(f));temp1(k)=1;low=floor(x(k));if isrowinmat([temp1,low],[A,b])==1    thisA=A;    thisb=b;else    thisA=[A;temp1];    thisb=b;    thisb(end+1)=low;end[x1,fval1,exitflag1,output1,lambda1]=linprogdis(ifint,f,thisA,thisb,Aeq,beq,lb,ub,x0,options);temp2=zeros(1,length(f));temp2(k)=-1;high=-ceil(x(k));if isrowinmat([temp2,high],[A,b])==1    thisA=A;    thisb=b;else    thisA=[A;temp2];    thisb=b;    thisb(end+1)=high;end[x2,fval2,exitflag2,output2,lambda2]=linprogdis(ifint,f,thisA,thisb,Aeq,beq,lb,ub,x0,options);if (isempty(v2) && ((exitflag1>0 && exitflag2<=0 && fval<=fval)||(exitflag2>0 && exitflag1<=0 && fval<=fval2)||(exitflag1>0 && exitflag2>0 && fval<=fval1 && fval<=fval2)))    disp(‘error call‘);    return ; %表示都是整数endif exitflag1>0&&exitflag2<=0     x=x1;     fval=fval1;     exitflag=exitflag1;     output=output1;     lambda=lambda1;elseif exitflag1<=0&&exitflag2>0     x=x2;     fval=fval2;     exitflag=exitflag2;     output=output2;     lambda=lambda2;elseif exitflag1>0 && exitflag2>0    if fval1<fval2        x=x1;        fval=fval1;        exitflag=exitflag1;        output=output1;        lambda=lambda1;    else         x=x2;         fval=fval2;         exitflag=exitflag2;         output=output2;         lambda=lambda2;    endend