arrays - 在 matlab 中创建 3D 矩阵

Question

我想绘制安全系数(FS，对特定区域发生山体滑坡风险的量化)随时间的演变情况。

计算公式如下:

effcohesion=0;
rootcohesion=0;
gammat=12.9E3;
gammaw=9810;
c=0;
deltac=0;
m=0.5;
z=2.5;
phi=16;
slope=rand(20,20)*30       % slope of a gridpoint in area

Strength = c + deltac + (gammat - gammaw.*m).*z.*(cosd(slope).^2);
Stress = gammat.*z.*(sind(slope)).*(cosd(slope));
Part = tand(phi);
FS2 = (Strength./Stress).*(Part)

现在。 m 的值(= 地下水位高度，决定 FS)全年变化，因此不是恒定的。我有一个包含降水量、蒸发量等数据的文件，但为了使其不太复杂，我在这里假设 m 只是一年中某一天的函数:

mnew=zeros(365,1);
for t=1:365
mnew(t)=(m+t)/150;
end

我现在有一个 FS 为 20x20 点的数据集，其中 m =0.5 (=FS2) 和一个包含一年中 m 演变的文件 (= mnew)。

我现在如何创建一个 3D 矩阵，其中 (1) 存储 FS 的空间变化(即 20x20 矩阵上的 FS 值)和 (2) 全年中 FS 随 m 函数变化的时间演变。最终，我想要一个矩阵，其中包含 FS 的空间和时间演化。

第 1 层 = 第 1 天所有 20x20 点的 FS

第 2 层 = 第 2 天 20x20 点的 FS

等等

有人可以帮助我吗？

提前致谢!

Answer 1

“3D 矩阵”更恰本地称为 3 阶数组。为此，只需将 FS2 计算结果粘贴到时间循环内即可。使用适当的 mnew 代替 m 来计算 FS2。然后将 FS3 层(等级 3 数组)设置为 FS2。

然后，第 1 层(第 1 天)由 FS3(:,:,1) 给出，第 2 层由 FS3(:,:,2) 给出，依此类推.

m0=0.5;

% Sizes of array
n1 = 20;
n2 = 20;
n3 = 365;

FS3 = zeros(n1, n2, n3);

mnew=zeros(n3,1);
for t=1:n3

    mnew(t)=(m0+t)/150;

    effcohesion=0;
    rootcohesion=0;
    gammat=12.9E3;
    gammaw=9810;
    c=0;
    deltac=0;

    m = mnew(t);

    z=2.5;
    phi=16;
    slope=rand(n1,n2)*30;       % slope of a gridpoint in area

    Strength = c + deltac + (gammat - gammaw.*m).*z.*(cosd(slope).^2);
    Stress = gammat.*z.*(sind(slope)).*(cosd(slope));
    Part = tand(phi);
    FS2 = (Strength./Stress).*(Part);

    FS3(:,:,t) = FS2;

end