Here's how I understood the question. If you have A = magic(3) and B = eye(3)--that is,
A = [8 1 6; 3 5 7; 4 9 2];
B = [1 0 0; 0 1 0; 0 0 1];
then you can subtract the rows of B from each row of A with
m = size(A,1);
C = repelem(A,m,1) - repmat(B,m,1);
This approach gives
C = [7 1 6; 8 0 6; 8 1 5; 2 5 7; 3 4 7; 3 5 6; 3 9 2; 4 8 2; 4 9 1];
But in the example given in the question, this approach would produce a 676x5100 matrix, not a 2601x5100 matrix.
