约束: 可以用MATLAB自带的整数规划函数intlinprog (找了半天用群体智能算的,发现没有合适的模板套,而且求出的也不一定是最优解,这个方法求出来的是最优解而且由于问题是线性的,所以求解也很快)