我在MATLAB代碼中有一個錯誤，但找不到

Question

我正在嘗試使用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個元素。 我看不到錯誤。 它在哪里？

謝謝。

Answer 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