2

はじめまして、ここに小さな問題があります。

ベクトル v (サイズ(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ループを使用してこれを行う方法は知っていますが、ベクトル化によっても可能だと思います-ここで問題は-私はベクトル化にあまり熟練していません。

ヒントなどの Thx :))

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 に答える