我在MATLAB代码中有一个错误,我找不到它

时间:2014-03-08 09:34:57

标签: matlab mathematical-optimization

我正在尝试使用MATLAB R2012b进行简单的线性编程优化。我想出了这段代码:

clear;
y=1663;
f=[100; 100; 100; 100; 100; 100; 100; 100];
A=[1/867.2 1/704.7 1/619.2 1/679 1/770.3 1/658.5 1/749 1/783.1
  -105.2 -74.4 -94.3 -88.2 -108.9 999 -112.7 -91.4 
  -94.1 -71.2 -87.6 -87.8 -99.4 999 -100 -81.2 
  0.5 0.6 4.5 3 28 15 4.5 12.5
  0.78 2.11 0.01 0.01 0.56 0.09 0.02 0.9
  93.63 16.3 0.01 0.01 1.44 0.09 0.03 37.64
  0.82 3.23 0.03 0.36 3.13 29.97 0.04 16.16];
b=[y/0.775; -95*y; -85*y; 10*y; 1*y; 35*y; 18*y];
Aeq=[1; 1; 1; 1; 1; 1; 1; 1];
beq=y;
lb=zeros(8,1);
hb=[550; 250; 50; 200; 50; 350; 400];
linprog(f, A, b, Aeq, beq, lb, hb);

当我尝试运行此源时,它表示 A中的行数必须与b的元素数相同。

A有7行,b有7个元素,可以看出。我看不出错误。它在哪里?

感谢。

1 个答案:

答案 0 :(得分:1)

你必须重新定义Aeq。由于您只有一个等式约束,因此它应该是1xn。我还将x作为linprog的输出。

   Aeq=[1, 1, 1, 1, 1, 1, 1, 1];
   x=linprog(f, A, b, Aeq, beq, lb, hb);

结果是

x =

  200.7466
   59.3486
   30.1325
   62.8837
   35.7733
    1.6669
  341.6570
  930.7914