索引矩阵 Matlab

Question

我正在尝试计算两个矩阵，其中一个是 400*2，另一个是 20*20。第一个包含 400 个点的 x 和 y 坐标，第二个包含另外 400 个点，其中 x 为 2nd(i,i)，y 为 (i,k)，i 和 k 都是 20 个数字计数器。我正在尝试计算这两个矩阵点之间的距离会给我 400*400 矩阵。我正在使用的代码是；

for i=1:20
    for j=1:400
        for k=1:20
        L(j,)=sqrt((C(j,1)-M(i,i))^2+(C(j,2)-M(k,i))^2);
        end
    end
end

C 是第一个矩阵，M 是第二个矩阵。现在我知道我让它听起来有点困难，但问题是我找不到一个计数器来给 L(j, ) 那部分。如果你们有任何想法请告诉...

编辑：首先认为有 2 个点向量。 X 为 -0.95:0.1:0.95 和 y 相同。 M 向量是所有 x 与 Y 配对，因此它将生成 20*20 矩阵或 400*2，但是我无法使其像 400*2 一样工作，因此使用 20*20。矩阵有点小我很难做到这一点，我对 M 的索引表明了这一点。如果您有更好的方法来配对它们，我也会很感激。

我的矩阵看起来像这样；

-0.95 -0.85 -0.75...
-0.95 -0.85 -0.75 ... 
.
.
.

Answer 1

上面的代码应该很慢......使用大量的for循环不是“matlab风格”
我要做的是：

Mdiag=repmat(diag(M)',20,1);
L=pdist2(C,[Mdiag(:) M(:)]);

第一行提取 M 的对角项并重复它们。

之后 [Mdiag(:) M(:)] 是第二组点的标准形式。

pdist2 是计算所有成对距离的函数。

Answer 2

它们的顺序重要吗？我假设不会，只要您知道该顺序是什么。我会用这个：

for i=1:20
    for j=1:400
        for k=1:20
            L(j,(i-1)*20+k) = sqrt((C(j,1)-M(i,i))^2+(C(j,2)-M(k,i))^2);
        end
    end
end