查看: 307|回复: 0|关注: 0

[未答复] 线性规划intlinprog函数问题

[复制链接]

新手

10 麦片

财富积分


050


2

主题

2

帖子

0

最佳答案
发表于 2018-7-20 14:49:48 | 显示全部楼层 |阅读模式
本帖最后由 胡焕_zWbdw 于 2018-7-20 14:50 编辑

各位大佬, 下面这个函数解出来的x现在是3组进行配比,  怎么修改这个函数固定计算出两组物资配比结果

  1. f = [630.07;24.99;14.0;630.07;660.0;630.07;465.68;631.45;750.07;875.5;630.07;630.07;630.07;630.07;630.07;0.0;411.23;635.71;648.78;696.33;987.65;854.66;987.65;1561.74;1561.74;772.19;860.92;609.7;609.7;639.33;679.96;157.52;0.0;621.44];

  2. Aeq = [1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0];

  3. beq =10;

  4. A= [-21.23 -0.842 -0.472 -21.23 -22.0 -21.23 -22.0 -22.0 -16.0 -19.68 -21.23 -21.23 -21.23 -21.23 -21.23 -0.0 -8.856 -21.8649 -21.8779 -15.197 -21.35 -19.35 -21.35 -18.88 -18.88 -16.3 -22.1 -21.071 -21.071 -21.059 -20.309 -5.308 -0.0 -20.939; 21.23 0.842 0.472 21.23 22.0 21.23 22.0 22.0 16.0 19.68 21.23 21.23 21.23 21.23 21.23 0.0 8.856 21.8649 21.8779 15.197 21.35 19.35 21.35 18.88 18.88 16.3 22.1 21.071 21.071 21.059 20.309 5.308 0.0 20.939;-0.95 -0.04 -0.0 -0.95 -1.0 -0.95 -0.7 -0.95 -2.0 -1.2 -0.95 -0.95 -0.95 -0.0 -0.0 -0.0 -0.21 -0.9162 -0.9246 -0.79 -0.0 -0.95 -0.0 -1.06 -1.06 -0.8 -1.2 -0.95 -0.95 -0.95 -0.95 -0.24 -0.0 -0.94;0.95 0.04 0.0 0.95 1.0 0.95 0.7 0.95 2.0 1.2 0.95 0.95 0.95 0.0 0.0 0.0 0.21 0.9162 0.9246 0.79 0.0 0.95 0.0 1.06 1.06 0.8 1.2 0.95 0.95 0.95 0.95 0.24 0.0 0.94;-21.65 -0.86 -0.48 -21.65 -15.0 -21.65 -15.0 -21.65 -17.0 -21.8 -21.65 -21.65 -21.65 -21.65 -21.65 -0.0 -3.28 -15.8262 -16.1546 -13.46 -16.98 -21.65 -16.98 -20.2 -20.2 -14.2 -21.4 -21.61 -21.61 -21.65 -21.65 -5.41 -0.0 -21.35;21.65 0.86 0.48 21.65 15.0 21.65 15.0 21.65 17.0 21.8 21.65 21.65 21.65 21.65 21.65 0.0 3.28 15.8262 16.1546 13.46 16.98 21.65 16.98 20.2 20.2 14.2 21.4 21.61 21.61 21.65 21.65 5.41 0.0 21.35];

  5. b= [-100.0;200.0;-5.0;15.0;-50.0;300.0];

  6. lb = [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;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;0.0;0.0;0.0;0.0;0.0;0.0;0.0;0.0];

  7. ub = [10.0;10.0;10.0;10.0;10.0;10.0;10.0;10.0;10.0;10.0;10.0;10.0;10.0;10.0;10.0;10.0;10.0;10.0;10.0;10.0;10.0;10.0;10.0;10.0;10.0;10.0;10.0;10.0;10.0;10.0;10.0;10.0;10.0;10.0];

  8. intcon = [1.0, 2.0, 3.0, 4.0, 5.0, 6.0, 7.0, 8.0, 9.0, 10.0, 11.0, 12.0, 13.0, 14.0, 15.0, 16.0, 17.0, 18.0, 19.0, 20.0, 21.0, 22.0, 23.0, 24.0, 25.0, 26.0, 27.0, 28.0, 29.0, 30.0, 31.0, 32.0, 33.0, 34.0];

  9. [x,fval,exitflag] = intlinprog(f,intcon, A,b,Aeq,beq,lb,ub);

  10. LP:                Optimal objective value is 2618.426829.                                          



  11. Cut Generation:    Applied 1 mir cut.                                                               

  12.                    Lower bound is 2625.468807.                                                      



  13. Heuristics:        Found 1 solution using rounding.                                                

  14.                    Upper bound is 3387.850000.                                                      

  15.                    Relative gap is 22.49%.                                                         



  16. Branch and Bound:



  17.    nodes     total   num int        integer       relative                                          

  18. explored  time (s)  solution           fval        gap (%)                                         

  19.        8      0.01         2   2.737740e+03   1.809572e+00                                          

  20.       93      0.01         2   2.737740e+03   0.000000e+00                                          



  21. Optimal solution found.



  22. Intlinprog stopped because the objective value is within a gap tolerance of the optimal value, options.AbsoluteGapTolerance = 0 (the default value). The intcon variables are integer within tolerance, options.IntegerTolerance = 1e-05 (the default value).



  23. >> x



  24. x =



  25.          0

  26.        5.0000

  27.          0

  28.          0

  29.          0

  30.          0

  31.         4.0000

  32.          0

  33.      1.0000

  34.          0

  35.          0

  36.          0

  37.          0

  38.          0

  39.          0

  40.          0

  41.          0

  42.          0

  43.          0

  44.          0

  45.          0

  46.          0

  47.          0

  48.          0

  49.          0

  50.          0

  51.          0

  52.          0

  53.          0

  54.          0

  55.          0

  56.          0

  57.          0

  58.          0

复制代码


您需要登录后才可以回帖 登录 | 注册

本版积分规则

关闭

站长推荐上一条 /3 下一条

快速回复 返回顶部 返回列表