我正在努力解决以下问题。
我有2个矩阵,让我们说A和B哪个维度可以及时改变,我想创建一个新的矩阵C,如下所示:
我想从矩阵A中的所有列中减去矩阵B中的所有列。因此,如果例如矩阵A是[1 2; 5 6; 8 5]和矩阵B [2 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
答案 0 :(得分:3)
我没有查看您的代码,但您可以使用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
为not necessary:
C = reshape(permute(A, [1 3 2]) - B, size(A,1), []);
答案 1 :(得分:0)
使用逐列相同大小的矩阵的示例非常容易处理:
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