algorithm - 用于创建数据矩阵的迭代算法

Question

假设我想创建以下数据矩阵

为了创建这种类型的矩阵，我已经编写了代码

function [ x ]=create_matrix1(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

但是当我们有大量数据时，创建这样的矩阵就成了大问题，例如，如果我们像这样从 matlab 生成数据

load ampoutput2.mat
[m n]=size(y)

m =

      500000


n =

     1

如果 p=2000 那么它当然会崩溃，计算机将停止工作，是否有任何迭代版本可以创建这个矩阵？非常感谢

Answer 1

我想你想要的是 hankel 函数。但是您可能不想显式构造矩阵，这取决于您想要用它做什么。通过调用 flip 并使用 FFT，您可以在 O(N log N) 时间内将这种类型的矩阵与向量相乘。