显示要串联的数组纬度不一致，如何解决

Question

原代码显示要串联的数组纬度不一致，修改后有显示the number of rows in A must be the same as the numberof elememts of b

Answer 1

引自免费微信小程序：皆我百晓生

在您提供的MATLAB代码中，存在两个主要问题导致了错误：

1. 错误地使用vertcat（垂直连接）函数。
2. 在定义约束条件矩阵Aineq时，行数与b的元素数量不匹配。

首先，让我们修复第一个问题。在原始代码中，您尝试对f进行重塑操作，但是使用了vertcat而不是正确的 reshape 函数。将这一行更改为正确使用 reshape 的语句：

f = -reshape(900 * ones(1, 12) - 450, [], 1);

接下来，我们处理第二个问题。根据错误信息，“The number of rows in A must be the same as the number of elements of b”，这表明您的约束矩阵Aineq的行数应与向量b的长度相同。由于b是一个4x1的向量，因此Aineq应该有4行。在当前的Aineq定义中，您创建了一个9x4的矩阵，而实际上只需要一个4x4的矩阵。我们可以重新构建Aineq来满足这个要求：

Aineq = [ones(4, 4) zeros(4, 4)];

这里，我们将所有水库的供水上限约束合并到同一个矩阵中。由于额外用水量已经包含在基本生活用水量的限制里，所以不再需要额外的约束。

修正后的完整代码如下：

c = [160 130 220 170; 140 130 190 150; 190 200 230 10000]; % 引水管理费
b = [30; 70; 10; 10]; % 基本生活用水量
a = [50; 60; 50]; % 水库供水上限
c_extra = [50; 70; 20; 40]; % 额外用水量

% 定义目标函数系数矩阵
f = -reshape(900 * ones(1, 12) - 450, [], 1);

% 定义约束条件矩阵
Aeq = [eye(4) eye(4) eye(4)]; % 基本生活用水量约束
Aineq = [ones(4, 4) zeros(4, 4)]; % 水库供水上限约束和额外用水量约束
beq = b; % 等式约束的常数项
bineq = a + c_extra; % 不等式约束的常数项

% 求解线性规划问题
[x, fval] = linprog(f, Aineq, bineq, Aeq, beq, zeros(12, 1));

% 输出结果
disp('最优供水方案：');
for i = 1:3
    for j = 1:4
        fprintf('从水库%d向区域%d供水：%0.2f\n', i, j, x((i-1)*4+j));
    end
end
disp(['最大利润：', num2str(-fval)]);

这段代码应该能成功运行并提供最优供水方案以及最大利润。如果您遇到其他问题，随时告诉我。