0

我正在尝试为每个条目在 2 个标量和一个矩阵之间绘制一个 3D 图。在编译时，第 11 行出现“子矩阵定义错误”错误。代码：

i_max= 3;
u = zeros(4,5);
a1 = 1;
a2 = 1;
a3 = 1;
b1 = 1;
hx = linspace(1D-6,1D6,13);
ht = linspace(1D-6,1D6,13);
for i = 1:i_max
    for j = 2:4
        u(i+1,j)=u(i,j)+(ht*(a1*u(i,j))+b1+(((a2*u(i,j+1))-(2*a2*u(i,j))+(a2*u(i,j-1)))*(hx^-2))+(((a3*u(i,j+1))-(a3*u(i,j-1)))*(0.5*hx^-1)));
        plot(ht,hx,u(i+1,j));
    end
end

完整的错误信息：

-->exec('C:\Users\deba123\Documents\assignments and lecture notes\Seventh Semester\UGP\Scilab\Simulation1_Plot.sce', -1)
+(((a3*u(i,j+1))-(a3*u(i,j-1)))*(0.5*hx^-1)))
                                          !--error 15 
Submatrix incorrectly defined.
at line      11 of exec file called by :    
emester\UGP\Scilab\Simulation1_Plot.sce', -1

请帮忙。

4

1 回答 1

0

对于 3 维图形，您需要 2 个参数向量和一个用于函数值的矩阵。所以我将 u 扩展为张量。在代码中的每个操作中，我都添加了术语的当前维度。现在，对您的计算进行透明处理。对于绘图，您必须使用plot3d（单个值）或surf（表面）命令。在 3-dim 图中，您需要两个 map 2 向量(hx,ht)，其中 n 和 m 为 scalar z。因此，您的结果达到了 (nxm) 矩阵。这是，你想做什么？目前，每个条目都有 13 个值u(i,j,:)，但每个数字都需要 (13x13)。也许eval3d-function 可以帮助你。

i_max= 3;
u = zeros(4,5,13);
a1 = 1;
a2 = 1;
a3 = 1;
b1 = 1;
hx = linspace(1D-6,1D6,13); // 1 x 13
ht = linspace(1D-6,1D6,13); // 1 x 13

for i = 1:i_max
    for j = 2:4
        u(i+1,j,:)= u(i,j)... 
                  + ht*(a1*u(i,j))*b1... // 1 x 13 
                  +(((a2*u(i,j+1)) -(2*a2*u(i,j)) +(a2*u(i,j-1)))*(hx.^-2))... // 1 x 13 
                  +(((a3*u(i,j+1))-(a3*u(i,j-1)))*(0.5*hx.^-1)) ... // 1 x 13   
                  + hx*ones(13,1)*ht; // added to get non-zero values
                  z = squeeze( u(i+1,j, : ))'; // 1x13

                  // for a 3d-plot: (1x13, 1x13, 13x13) 
                  figure()
                  plot3d(ht,hx, z'* z ,'*' ); //

    end
end
于 2015-10-12T09:49:18.933 回答