Question

我需要一个三维矩阵，但第一个维度不一样。所以我说NxT1（N乘T1），NxT2 NxT3和NxT4。我想将它们放在一个矩阵中，这样我就可以遍历每个维度。这是我的代码：

y2(:,:,1) = zeros(N,T2(:,1));
y2(:,:,2) = zeros(N,T2(:,2));
y2(:,:,3) = zeros(N,T2(:,3));
y2(:,:,4) = zeros(N,T2(:,4));
y2(:,1,:) = c/(1-rho);

for z=1:size(T2,2)
    for i=2:T2(:,z)
        for j=1:N
            y2(j,i,z) = y2(j,i-1)+randn;
        end
    end
end

我想基本上随机游走不同的时间范围。 T2 = [50,100,150,200]所以我的3维矩阵将包含4种不同时间规格的N次模拟。

Answer 1

我认为你想要的是一个数组，而不是一个矩阵。

c = 1.0;
rho = 0.5;
N = 100;
T2 = [50, 100, 150, 200];
for i = [1:length(T2)];
    y2{i} = zeros(N, T2(i));
    y2{i}(1,:) = c/(1-rho);
end;

for i = [1:length(T2)];
    for j = [2:N];
        for k = [1:T2(i)];
            y2{i}(j,k) = y2{i}(j-1,k) + randn()
        end;
    end;
end;

MATLAB，具有不同二维的三维矩阵

1 个答案: