该楼层疑似违规已被系统折叠 隐藏此楼查看此楼
这是代码：
clc,clear;n1=input('Ⅰ型车的数量n1=');
n2=input('Ⅱ型车的数量n2=');
n3=input('Ⅲ型车的数量n3=');
a=[-4,-3,-3,-2,-2,-2,-1,-1,-1,-1,0,0,0,0,0,-5,-4,-4,-3,-3,-3,-2,-2,-2,-2,-1,-1,-1,-1,-1,0,0,0,0,0,0,-4,-3,-2,-1,0,-4,-3,-2,-1,0,-5,-4,-3,-2,-1,0;
0,-1,0,-2,-1,0,-3,-2,-1,0,-5,-3,-2,-1,0,0,-1,0,-2,-1,0,-4,-2,-1,0,-5,-4,-2,-1,0,-6,-5,-4,-2,-1,0,0,-1,-2,-3,-5,0,-1,-2,-3,-5,0,-1,-2,-4,-5,-6;
0,0,-1,0,-1,-2,0,-1,-2,-3,0,-1,-2,-3,-4,0,0,-1,0,-1,-2,0,-1,-2,-3,0,-1,-2,-3,-4,0,-1,-2,-3,-4,-5,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0;
-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0];
aeq=[1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,-1,-1,-1,-1,-1,0,0,0,0,0,0;
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,0,0,0,0,0,0,0,0,0,0,-1,-1,-1,-1,-1,-1]
f=[1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1;1]
intcon=1:52;
b=[-n1 -n2 -n3 0]';
beq=[0 0]';
lb=zeros(52,1);
[x,fval]=intlinprog(f,intcon,a,b,aeq,beq,lb)
现在需要增加一个约束条件：x37+x38+...+x41 为4的整数倍。求大神指点。