Question

我有通过管道的流体速度的3D流量数据。我知道管的直径，并查看了速度场，并找到了管两端xy平面的场中心。所以我基本上有一条穿过管子中心轴的线。我想要NaN直径以外的所有数据点。为此，我使用的方程式给出了距3D中一条点的距离，我在这里找到了mathworld.wolfram.com/Point-LineDistance3-Dimensional.html。然后我创建了一个if语句，表明小于直径的点将是NaN。

我是matlab的新手所以我不知道我现在将如何绘制这个。

%%

diff_axis = end_axis-start_axis;
diff_axis_mag = (diff_axis(1)^2 + diff_axis(2)^2 + diff_axis(3)^2)^0.5;

[rw col pl] = size(X);

for j = 1:col

    for i = 1:rw

        for k = 1:pl

             x_curr = X(i,j,k);
             y_curr = Y(i,j,k);
             z_curr= Z(i,j,k);

             x0 = [x_curr y_curr z_curr]

             t = - dot((start_axis-x0),(diff_axis))./(diff_axis_mag)^2;

             d = sqrt(((start_axis(1) - x0(1)) + (end_axis(1) - start_end(1))*t)^2 + ((start_axis(2)-x0(2))+(end_axis(2)-start_end(2))*t)^2+((start_axis(3)-x0(3))+(end_axis(3)-start_end(3))*t)^2);

             if (d > D)

                 x_curr=NaN
                 y_curr=NaN
                 z_curr=NaN

             end




        end

    end


end

Answer 1

很高兴为您的X，Y和Z提供解释性名称。我猜他们是流量组件，diff_axis是轴坐标？这是一个非常麻烦的符号。

你在循环中做的是你取点值（X,Y,Z），将它们复制到临时常量，然后如果它们掉出来将它们设置为NaN。但问题是通常你不在MATLAB中逐点绘制。所以像x_curr这样的临时家伙将会失败。此外，在MATLAB中执行操作的最佳方法是尽可能避免循环。

你可以做的是先创建一个面具

%// remember to put a dot like in `.^` for entrywise array operations
diff_axis_mag = sqrt(diff_axis(1).^2 + diff_axis(2).^2 + diff_axis(3).^2);
%// are you sure you need to include the third axis? 
%// then it is a ball, not a tube

%// create a binary mask
mask = diff_axis_mag < tube_radius

X(~mask) = NaN;
Y(~mask) = NaN;
Z(~mask) = NaN;

然后，您可以使用quiver3或 stream3

我如何从下面的代码中绘制for循环？

1 个答案: