获取给定索引前后(窗口+/-1)的向量索引

您可以通过隐式单例扩展(R2016b或更新的MATLAB,固有于倍频程)来实现这一点

idx = [2, 6, 9]; % some vector of integers
% Use reshape with [] to tell MATLAB "however many rows it takes"
neighbours = reshape( idx + [-1;1], [], 1 );

>> neighbours = [1; 3; 6; 8; 8; 10];

如果你不知道 idx 是行或列,可以通过使用

neighbours = reshape( idx(:)' + [-1,1], [], 1)

如果不想使用隐式扩展(并再次处理行或列 idx公司 ),可以像这样使用重塑

neighbours = reshape( [idx(:)-1, idx(:)+1]', [], 1 )

注意:您可能还希望将整个过程包装在对的调用中 unique . 在我的示例中,您可以得到索引 8 有两次,我不确定这在你的情况下是可取的还是不可取的。

然而 唯一的 执行排序(除非使用 'stable' 标记,但这可能会使其变得更慢),因此如果要删除重复项,最好使用原始方法:

% Remove duplicates and sort the result using unique 
neighbours = unique( [idx-1, idx+1] );

嗯,我终于找到了这个倍频程矩阵操作:

vec([idx-1, idx+1]')
ans =

5297
5299
6275
6277
6832
6834
...
20187
20189

将Wolfie的解决方案改编为最短的纯倍频程代码:

[idx-1, idx+1]' (:)  
( idx(:)' + [-1; 1] )(:) 

idx = ( find(diff(data(:,3)) > 0 )' + [-1; 1] )(:) 用作一个衬里

... 和 [idx , data(idx,3)] 并排显示索引和数据