matlab 矩阵列减,matlab – 如何从矩阵的每一行中减去一个向量？

Possible Duplicate:

07000

我正在寻找一种优雅的方式从矩阵的每一行中减去相同的向量。这是一种非常优雅的做法。

a = [1 2 3];

b = rand(7,3);

c(:,1) = b(:,1) - a(1);

c(:,2) = b(:,2) - a(2);

c(:,3) = b(:,3) - a(3);

此外，优雅的方式不能比这种方法慢。

我试过了

c = b-repmat(a,size(b,1),1);

似乎较慢。

编辑：赢家是这种方法。

c(:,1) = b(:,1) - a(1);

c(:,2) = b(:,2) - a(2);

c(:,3) = b(:,3) - a(3);

编辑：更多的方法，和tic toc结果：

n = 1e6;

m = 3;

iter = 100;

a = rand(1,m);

b = rand(n,m);

tic

c = zeros(size(b));

for i = 1:iter

c(:,1) = b(:,1) - a(1);

c(:,2) = b(:,2) - a(2);

c(:,3) = b(:,3) - a(3);

end

toc

tic

c = zeros(size(b));

for i = 1:iter

c(:,1) = b(:,1) - a(1);

c(:,2) = b(:,2) - a(2);

c(:,3) = b(:,3) - a(3);

end

toc

tic

c = zeros(size(b));

for i = 1:iter

for j = 1:3

c(:,j) = b(:,j) - a(j);

end

end

toc

tic

for i = 1:iter

c = b-repmat(a,size(b,1),1);

end

toc

tic

for i = 1:iter

c = bsxfun(@minus,b,a);

end

toc

tic

c = zeros(size(b));

for i = 1:iter

for j = 1:size(b,1)

c(j,:) = b(j,:) - a;

end

end

toc

结果

Elapsed time is 0.622730 seconds.

Elapsed time is 0.627321 seconds.

Elapsed time is 0.713384 seconds.

Elapsed time is 2.621642 seconds.

Elapsed time is 1.323490 seconds.

Elapsed time is 17.269901 seconds.