在R中动态构建线性规划约束矩阵的实现方法咨询

首先得指出你原代码里的两个小问题：一是if (sum(cars[[i]])) > 0的语法错误，括号位置不对，应该写成if (sum(cars[[i]]) > 0)；二是提前预定义矩阵大小的思路会导致索引错位——因为你跳过了零向量对应的i，矩阵里会留下空行或者NA，反而达不到你要的效果。

针对你的需求，我给你两种实用的方案，不用提前管矩阵大小，自动剔除零向量：

方案一：先筛选再绑定（高效推荐） #

既然你只想保留非零向量的元素，不如先把列表里符合条件的元素筛选出来，再直接合并成矩阵。这种方法比循环动态增长效率高得多，尤其适合大规模列表：

# 第一步：筛选出所有非零向量（若需严格判断全零，可改用!all(vec == 0)）
filtered_cars <- cars[sapply(cars, function(x) sum(x) > 0)]

# 第二步：把筛选后的列表转成矩阵（前提是所有向量长度一致，符合约束矩阵要求）
B <- do.call(rbind, filtered_cars)

sapply(cars, function(x) sum(x) > 0)会生成一个逻辑向量，标记每个列表元素是否是非零向量；用这个逻辑向量索引列表，就得到了只包含有效约束的子列表；最后do.call(rbind, ...)把列表里的向量按行合并成矩阵。

方案二：动态追加行（适合坚持用循环的场景） #

如果你一定要用循环的方式逐步添加行，可以初始化矩阵为NULL，每次遇到非零向量就用rbind追加。注意这种方法在列表规模很大时效率会低一些，因为每次rbind都会复制整个矩阵：

# 初始化空矩阵
B <- NULL

# 遍历列表里的每个向量
for (vec in cars) {
  # 判断是否为非零向量（可根据需求调整判断逻辑）
  if (sum(vec) > 0) {
    # 追加到矩阵末尾
    B <- rbind(B, vec)
  }
}

这里直接遍历列表元素，不用管索引，就避免了原代码里索引错位的问题，而且不需要提前定义矩阵大小。

另外补充一点：如果你的零向量判断是“所有元素都为0”，用all(vec == 0)来判断会更严谨（比如避免向量元素和为0但不是零向量的情况，比如c(1,-1)），可以把判断条件改成if (!all(vec == 0))，根据你的实际需求调整就行。

内容的提问来源于stack exchange，提问作者Jorden