%yalmip设置

 options = sdpsettings('solver' ,'cplex');

 options.showprogress = 1;%1为设置显示yalmip现在在做什么

 options.verbose = 2;%设置显示信息程度，1为适度显示，2为完全显示。

 [model,recoverymodel,diagnostic,internalmodel] = export(Constraints,Objective,options);%输出Yalmip模型

 %定义一个Cplex类实体

 cplex_milp = Cplex('Milp for HTC');

 cplex_milp.Model.sense = 'minimize';

 cplex_milp.Model.obj = model.f;

 cplex_milp.Model.lb = model.lb;

 cplex_milp.Model.ub = model.ub;

 cplex_milp.Model.A = [model.Aineq;model.Aeq];

 cplex_milp.Model.lhs = [-Inf.*ones(size(model.bineq,1),1);model.beq];

 cplex_milp.Model.rhs = [model.bineq;model.beq];

 cplex_milp.Model.ctype = model.ctype';
 cplex_milp.Param.mip.tolerances.mipgap.Cur = zk.solver_error;%相对MIP间隔容差，（精度）

 cplex_milp.Param.threads.Cur = 1;%设置将由任何 CPLEX 并行优化器调用的并行线程的缺省数量。1为单线程

 cplex_milp.Param.mip.display = 1;%决定在混合整数优化期间cplex报告到屏幕的内容。默认为2.

 cplex_milp.Param.timelimit = 15;%设置对优化器的调用的最长时间（以秒为单位）

 cplex_milp.writeModel('hmilp.lp');%输出数学模型

 cplex_milp.solve();