xa是31维矩阵 xa是决策变量 约束条件 sum(xa)-sum(xa,2)'== 0 ，matlab使用linprog代码

如果您需要使用 `linprog` 函数来求解该约束条件的最优解，可以使用以下 MATLAB 代码：

n = 31; % 矩阵维度
f = ones(n, 1); % 目标函数
Aeq = ones(1, n) - eye(n); % 约束条件矩阵
beq = zeros(1, 1); % 约束条件值
lb = zeros(n, 1); % 变量下界
ub = ones(n, 1); % 变量上界
options = optimoptions('linprog', 'Display', 'off'); % 设置优化选项
[x, fval, exitflag, output] = linprog(f, [], [], Aeq, beq, lb, ub, [], options); % 求解

其中，`n` 表示决策变量 `xa` 的维度，`f` 表示目标函数，`Aeq` 表示约束条件矩阵，`beq` 表示约束条件值，`lb` 和 `ub` 分别表示变量的下界和上界。

在本例中，约束条件矩阵 `Aeq` 的构造方式为 `ones(1, n) - eye(n)`，其中 `ones(1, n)` 表示一个长度为 `n` 的全 1 向量，`eye(n)` 表示一个 `n` 阶单位矩阵。所以 `Aeq` 的每一行都为全 1 向量减去单位矩阵的一行，即 `sum(xa) - xa(i) == 0`，其中 `i` 表示第 `i` 行。

最后，调用 `linprog` 函数求解得到最优解 `x`，`fval` 为最优值，`exitflag` 表示求解器的退出标志，`output` 包含求解器的详细信息。