2

Geetings,这是一个小问题。

我有一个向量 v (size(v) = T),一个正数 P,P < T 和正数 N,N < T。并且想要制作大小为 PxN 的矩阵 M,这样:

M = [v(T),     v(T-1), ....., v(T-N+2),   v(T-N+1);
     v(T-1),   v(T-2), ....., v(T-N+1),   v(T-N)  ;
     v(T-2),   v(T-3), ....., v(T-N),     v(T-N-1);
        .        .       .       .           .    ;
        .        .       .       .           .    ;
     v(T-P+1), v(T-P), ....., v(T-P-N+3), v(T-P-N+2)]

它认为 T > P + N。

举个例子吧。

v = [1, 2, 3, 4, 5]';
P = 3;
N = 3;
M = [5, 4, 3;
     4, 3, 2;
     3, 2, 1]

我知道如何使用 for 循环来做到这一点,但我也认为可以通过矢量化来做到这一点 - 这里的问题是 - 我在矢量化方面不太熟练。

谢谢提示等:))

4

2 回答 2

7

使用汉克尔矩阵:

v = [1 2 3 4 5];
T = length(v);
P = 3;
N = 3;
out = hankel(v(T:-1:T-P+1),v(T-P+1:-1:T-P-N+2));
于 2013-10-24T08:28:55.840 回答
3

你可以利用bsxfun

v = [1, 2, 3, 4, 5]';
P = 3;
N = 3;
idx = bsxfun(@minus, length(v):-1:length(v)-N+1, (0:P-1)');

result  = v(idx);

这将导致

result =

     5     4     3
     4     3     2
     3     2     1
于 2013-10-24T08:28:14.387 回答