1

假设我们有以下一维数组

x(1),x(2),......x(n)

其中n是样本的长度,假设我们要使用以下规则创建由这些元素组成的矩阵,使用一些参数L,我们的二维数组由以下元素组成

(1 row) x(1),x(2),.....x(l)
(2 row)  x(2),x(3),.....x(l+1)
(3 row) x(3),x(4),.......x(l+2)
.
.
.
.(m row)  x(m),x(m+1).... x(n)

我试过以下代码

function [ x ]=create_matrix(b,n,L)
for i=1:n/2 
    x(i,end)=b(i:L);
end;

其中 b 是输入一维向量,n 是长度(b),L 是矩阵行长度,但是当我运行此代码时,出现以下错误

??? Undefined function or variable "x".

Error in ==> create_matrix at 4
    x(i,end)=b(i:L);

也许在matlab中声明二维数组有问题,或者其他什么?实际上我知道实现起来很简单,但我自己做不到,请帮忙

4

2 回答 2

2

代码:

function [ x ]=create_matrix(b,l)
    n = length(b);
    m = n-l+1;
    x = zeros(m,l);
    for i=1:m
        x(i,:)=b(i:i+l-1);
    end;
end

例子:

EDU>> create_matrix(1:10,4)

ans =

     1     2     3     4
     2     3     4     5
     3     4     5     6
     4     5     6     7
     5     6     7     8
     6     7     8     9
     7     8     9    10

此外,这个矩阵是一个翻转的托普利兹矩阵。所以同样的事情可以通过以下方式实现:

fliplr(toeplitz(4:10,4:-1:1))

输出:

ans =

     1     2     3     4
     2     3     4     5
     3     4     5     6
     4     5     6     7
     5     6     7     8
     6     7     8     9
     7     8     9    10
于 2013-03-22T17:29:34.347 回答
2

另一种选择(可能在性能上更快)是使用bsxfun,例如

m=7;
l=4;
bsxfun(@plus,[1:m]',0:l-1)

ans =
 1     2     3     4
 2     3     4     5
 3     4     5     6
 4     5     6     7
 5     6     7     8
 6     7     8     9
 7     8     9    10

或用于通用矢量x使用

x(bsxfun(@plus,x(1:m)',0:x(l-1)))
于 2013-03-22T18:36:30.323 回答