我需要取向量子集的最大值，然后在向量中移动。例如列向量
a=[1;2;3;4;5;6;7]
如何获取max(a(1:3)), max(a(2:4)), ...max(a(5:end))并将所有输出放在另一个向量中？我可以使用 for 循环轻松地做到这一点，但是我正在寻找一种优雅的方式来使用 MATLAB 中的矩阵运算，最好是在一行代码中（尽管我意识到 MATLAB 中的矩阵运算可能会这样做正在使用 for 循环来实现）。
谢谢！
问问题
1207 次
5 回答
2
假设您的移位和窗口长度是恒定的（在您的情况下分别为 1 和 3），您可以使用nlfilter在向量上定义一般的滑动窗口操作：
a = [1:7];
fun = @(x) max(x(:));
b = nlfilter(a, [1 4], fun);
于 2012-10-12T22:04:00.403 回答
1
对于您提供的示例：
b = max([a(1:end-2) a(2:end-1) a(3:end)], [], 2)
我觉得还有一个更优雅、更通用的accumarray解决方案，虽然我现在没有时间想出一个:)
于 2012-10-13T10:10:15.050 回答
1
这是我的尝试。首先，定义“窗口”长度（在您的情况下为 3）。接下来，用于为每个窗口bsxfun的向量创建索引。a最后一步是一个简单的max函数调用：
flen = 3;
idx = bsxfun(@plus, [0:flen-1]', 1:numel(a)-flen+1);
max(a(idx))
它是通用的，即您可以更改窗口大小。
它是基本的 matlab，应该比 快nlfilter，虽然我现在无法验证它。
编辑这里是上述bsxfun解决方案和建议的简要性能比较im2col：
a = rand(10^7,1);
tic;
idx = bsxfun(@plus, [0:2]', 1:numel(a)-2);
m1 = max(a(idx));
toc
tic;
m2 = max(im2col(a,[3 1],'sliding'));
toc;
isequal(m1, m2)
Elapsed time is 0.839869 seconds.
Elapsed time is 1.797665 seconds.
ans =
1
bsxfun当然可以，而且速度快两倍多。
于 2012-10-13T11:33:36.907 回答
0
如果您有信号处理工具箱，则可以使用该buffer功能。它接受一个向量并生成一个矩阵，其中每一列都包含滑动窗口内的向量值。
在你的情况下：
>> a = (1:7)';
>> buffer(a,3,2,'nodelay')
ans =
1 2 3 4 5
2 3 4 5 6
3 4 5 6 7
然后，只需应用 MATLAB 的max函数来计算每列的最大值。总之，只需一行代码，
>> max(buffer(a,3,2,'nodelay'))
ans =
3 4 5 6 7
于 2015-01-06T15:00:10.640 回答
-1
感谢所有的帮助。这给了我一些想法，尤其是使用 nlfitler 的建议。事实证明，nlfitler 调用 im2col 并取其最大值对我有用。特别是因为提供的样本向量是包含随机变量的实向量的简化，因此存在的 bsxfun 选项将不起作用。
最大值（im2col（a，[3 1]，'滑动'））
于 2012-10-13T17:08:30.567 回答