MATLAB矩阵列减法

Question

我正在努力解决以下问题。

我有两个矩阵，假设A和B的维数可以随时间变化，我想创建一个新的矩阵C，如下所示：

我想从矩阵A的所有列中减去矩阵B中的所有列，因此，如果矩阵A为1 2;5 6;8 5和矩阵B2 2;2 1;3 5，则矩阵C为1-2 1- 2 -2 -2- 2-2;5-2 5-1 6-2 6-1;8-3 8-5 5-3 5-5。

我试过这段代码，但是在某个地方有个错误，我找不到它。请您更正代码，使其工作，或提供一个更好的解决方案，以获得正确的结果？

  matrixA=[4 5; 4 7;9 6];
  matrixB=[3 6 4;5 8 9;4 1 2];

[m,n]=size(matrixA);
[o,p]=size(matrixB);

C = zeros(m,n*p);
for p = 1:n
  for q = 1:p
    C(:,3*p+q-3) = matrixA(:,p)-matrixB(:,q);
  end
end

Answer 1

我还没有研究过您的代码，但是您可以使用permute、bsxfun、reshape实现更简单和矢量化的代码

A = [1 2;5 6;8 5];
B = [2 2;2 1;3 5];
C = reshape(bsxfun(@minus, permute(A, [1 3 2]), B), size(A,1), []);

从R2016b开始，bsxfun是不必了

C = reshape(permute(A, [1 3 2]) - B, size(A,1), []);

Answer 2

列大小相等的矩阵的示例非常容易处理：

matrixA = [1 2; 5 6; 8 5];
matrixB = [2 2; 2 1; 3 5];

A = repelem(matrixA,1,2);
B = repmat(matrixB,1,2);
C = A - B

C =
     -1    -1     0     0
      3     4     4     5
      5     3     2     0

但是对于处理不匹配的列大小，事情会变得更复杂一些：

matrixA = [1 2; 5 6; 8 5];
matrixB = [2 2; 2 1; 3 5];

A = repelem(matrixA,1,size(matrixB,2));
B = repmat(matrixB,1,size(matrixA,2));
C = (B - A) .* -1

C =
    -1    -1     0     0
     3     4     4     5
     5     3     2     0

matrixA = [4 5; 4 7; 9 6];
matrixB = [3 6 4; 5 8 9; 4 1 2];

A = repelem(matrixA,1,size(matrixB,2));
B = repmat(matrixB,1,size(matrixA,2));
C = (B - A) .* -1

C =
     1    -2     0     2    -1     1
    -1    -4    -5     2    -1    -2
     5     8     7     2     5     4