如何在MATLAB中高效查找数组指定索引区间的最大值（无循环）

首先得指出你原代码的问题：pv = max(w(s_on(1:end):s_off(1:end)))之所以失效，是因为Matlab会把这个索引表达式解析成从s_on(1)到s_off(end)的单一连续区间，而不是逐个对应s_on(i)到s_off(i)的独立区间。比如你的示例里，这会直接取w中第5到152位的全局最大值，完全不是你想要的每个区间的最大值。

下面给你几种无需显式循环的解决方案：

方法1：使用accumarray（高效推荐） #

accumarray是Matlab处理分组计算的神器，适合大规模数据或大量区间的场景：

% 示例数据
s_on = [5 19 78 101];
s_off = [10 28 97 152];
w = rand(1,200);

% 生成每个元素对应的区间标签（无循环）
starts = zeros(1, max(s_off) + 1);
starts(s_on) = 1;
starts(s_off + 1) = -1;
labels = cumsum(starts(1:end-1)); % 每个位置的区间编号，0表示不属于任何区间

% 过滤有效区间元素并分组求最大值
valid_mask = labels > 0;
pv = accumarray(labels(valid_mask)', w(valid_mask)', [], @max);

原理：通过累加起始/结束标记生成区间编号，再用accumarray按编号分组取最大值，全程没有显式循环，效率拉满。

方法2：使用arrayfun（简洁直观） #

如果你的区间数量不多，arrayfun的写法会更简洁（本质是隐式循环，但代码可读性高）：

pv = cell2mat(arrayfun(@(a,b) max(w(a:b)), s_on, s_off, 'UniformOutput', false));

原理：逐个对s_on和s_off的对应元素调用匿名函数，计算单个区间的最大值，最后把返回的cell数组转成普通数组。

方法3：矩阵索引法（适合区间长度相近的场景） #

通过构造索引矩阵，把每个区间的元素放到矩阵的一行，再逐行取最大值：

max_interval_len = max(s_off - s_on + 1);
% 构造每个区间的索引矩阵，不足长度的位置用NaN填充（max会忽略NaN）
idx_matrix = bsxfun(@plus, s_on, 0:max_interval_len-1);
% 把超出区间的索引对应的w值设为NaN
w_padded = [w, NaN(1, max_interval_len)];
pv = max(w_padded(idx_matrix), [], 2);

原理：利用矩阵的广播特性生成索引，用NaN填充无效位置后，逐行取最大值即可得到每个区间的结果。

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