如何将matlab中的一个向量plot以另一个向量作为参数？

Question

我正在尝试优化我正在编写的 function 的速度，并尝试尽可能多地使用向量。 我是 Matlab 的新手，矢量化有时对我来说是可以理解的，但我需要一些额外的帮助。 这是我当前的代码：

请注意，oracle() function 代表一个随机形状的 object，如果您输入一个 1x2 矩阵，它将返回矩阵（或在我们的例子中，x-y 坐标）是否在随机 object 内。

图片中的代码

function area = MitchellLitvinov_areaCalc(n)
    % create random coordinate vectors, with bounds from (3, 14)
    x = rand(n, 1) * 11 + 3;
    y = rand(n, 1) * 11 + 3;

    % find every point that is inside of oracle
    inOracle = oracle([x y]);

    % calculate the proportion, and multiply total area by proportion to find area
    % of oracle
    numPointsInOracle = nnz(inOracle);
    area = numPointsInOracle/n * (11*11);

    % create variable to store number of points in the area, and create a
    % matrix with size [numPoints, 2] to hold x and y values
    oracleCoordinates = zeros(numPointsInOracle, 2);

    % HERE IS WHERE I NEED ASSISTANCE!!!!!!!!!
    % find the points that are in the oracle shape
    index = 0; % start index at 0 % is the index of our oracleCoordinates matrix
    for i = 1:n % have to go through every point again to get their index

        % if point is inside oracle, increase index and record
        % coordinates
        if (inOracle(i) == 1) % see if point is in oracle
            index = index + 1; 
            oracleCoordinates(index, 1) = x(i, 1);
            oracleCoordinates(index, 2) = y(i, 1);
        end
    end

    % plot all points inside the oracle
    scatter(oracleCoordinates(:,1), oracleCoordinates(:,2))
    title("Oracle Shape")
    xlim([3, 14]);
    ylim([3, 14]);
end

是的，即使接近最大使用量 memory，代码也会运行得相当快。 但出于速度原因，我希望它完全矢量化，如果我需要重新使用此代码进行成像。 目前，要计算我使用矢量的面积，但要实际再现图像，我需要创建一个单独的存储矩阵，并手动使用索引/追加来传输 oracle function 内的点。我想知道是否有任何直接的“快捷方式”都可以让我的绘图更快一些。

Answer 1

您可以使用一个数组作为另一个数组中 select 某些项目的索引。 例如，使用您的变量名称：

oracleCoordinates(:,1) = x(inOracle == 1);
oracleCoordinates(:,2) = y(inOracle == 1);

这应该给出与您问题中的代码相同的结果，而无需使用循环。