本问题已经有最佳答案，请猛点这里访问。

我有一个向量，例如

vector = [1 2 3]

我想在自身内部重复n次，即如果n = 3，则最终结果为：

vector = [1 2 3 1 2 3 1 2 3]

对于任何n值，如何实现呢？ 我知道我可以执行以下操作：

newvector = vector;

for i = 1 : n-1

newvector = [newvector vector];

end

不过，这似乎有点麻烦。 还有更有效的方法吗？

尝试

repmat([1 2 3],1,3)

我将让您检查repmat的文档。

这是比repmat或reshape快一个数量级的方法

进行此类操作的最佳方法之一是使用Tony's Trick。通常发现Repmat和Reshape比Tony的技巧要慢，因为它直接使用Matlabs固有的索引。为了回答您的问题，

假设您要像r=[1 2 3 1 2 3 1 2 3...]那样平铺行向量r=[1 2 3] N的时间，然后，

c=r'

cc=c(:,ones(N,1));

r_tiled = cc(:)';

对于大型N，此方法与reshape或repmat相比可节省大量时间。

编辑：回复@ Li-aung Yip的疑问

我进行了一个小型的Matlab测试，以检查repmat和tony's trick之间的速度差异。使用下面提到的代码，我计算了从基本向量A=[1:N]构造相同平铺向量的时间。结果表明，是的，托尼的绝技更快，尤其适用于较大的N。欢迎人们自己尝试一下。如果必须循环执行这种操作，那么这么多的时间差就很关键。这是我使用的小脚本；

N= 10 ;% ASLO Try for values N= 10, 100, 1000, 10000

% time for tony_trick

tic;

A=(1:N)';

B=A(:,ones(N,1));

C=B(:)';

t_tony=toc;

clearvars -except t_tony N

% time for repmat

tic;

A=(1:N);

B=repmat(A,1,N);

t_repmat=toc;

clearvars -except t_tony t_repmat N

两种方法的时间(以秒为单位)如下：

N = 10，time_repmat = 8e-5，time_tony = 3e-5

N = 100，time_repmat = 2.9e-4，time_tony = 6e-5

N = 1000，time_repmat = 0.0302，time_tony = 0.0058

N = 10000，time_repmat = 2.9199，time_tony = 0.5292

我的内存不允许我超过N = 10000。我相信，对于N = 100000，两种方法之间的时间差会更大。我知道，这些时间对于不同的机器可能会有所不同，但是时间顺序的相对差异将会存在。而且，我知道，平均时间本来可以是更好的指标，但是我只想显示两种方法之间在时间消耗上的数量级差异。我的机器/操作系统的详细信息如下：

相关机器/操作系统/ Matlab详细信息：Athlon i686 Arch，Ubuntu 11.04 32位，3GB内存，Matlab 2011b

首先;什么？！其次，人们想知道为什么repmat()比手动执行要慢。第三，任何使用此功能的人最好在其旁边添加评论。

在对Tonys Trick进行了一些研究之后，似乎在14年前撰写本文时速度更快。自那时以来，MATLAB已有了很大改进，Tonys Trick可能不再比repmat快。 (您应该编写一个基准并对此进行测试。))

会的，谢谢你的链接。

@ Li-aungYip是的，数学论坛上的发布日期是2010年12月13日。此外，回复您的早期评论：C是40年前编写的，这并不意味着它比Matlab要慢。不要依靠观点/当代传统智慧做出客观决定。

很好的答案：关于程序性能的陈述有证据支持。我在SO上看到很多关于程序性能的断言，很少有证据支持。我看不到自己经常使用Tonys Trick，但后来我并没有太多使用repmat，现在我知道如果发现repmat太慢，该去哪里找。

好吧，我该死。辛苦了:)

2018年更新我复制了上面的实验(一定要将您的代码计时到函数M文件中！)，并使用repmat始终获得更快的结果。就是说，即使repmat已过时，我们已有bsxfun十多年了，现在却隐含了单例扩展。如今，repmat的合法用法很少而且相去甚远。

基于Abhinav的答案和一些测试，我编写了一个始终比repmat()快的函数！

它使用相同的参数，除了第一个参数必须是向量而不是矩阵。

function vec = repvec( vec, rows, cols )

%REPVEC Replicates a vector.

%   Replicates a vector rows times in dim1 and cols times in dim2.

%   Auto optimization included.

%   Faster than repmat()!!!

%

%   Copyright 2012 by Marcel Schnirring

if ~isscalar(rows) || ~isscalar(cols)

error('Rows and cols must be scaler')

end

if rows == 1 && cols == 1

return  % no modification needed

end

% check parameters

if size(vec,1) ~= 1 && size(vec,2) ~= 1

error('First parameter must be a vector but is a matrix or array')

end

% check type of vector (row/column vector)

if size(vec,1) == 1

% set flag

isrowvec = 1;

% swap rows and cols

tmp = rows;

rows = cols;

cols = tmp;

else

% set flag

isrowvec = 0;

end

% optimize code -> choose version

if rows == 1

version = 2;

else

version = 1;

end

% run replication

if version == 1

if isrowvec

% transform vector

vec = vec';

end

% replicate rows

if rows > 1

cc = vec(:,ones(1,rows));

vec = cc(:);

%indices = 1:length(vec);

%c = indices';

%cc = c(:,ones(rows,1));

%indices = cc(:);

%vec = vec(indices);

end

% replicate columns

if cols > 1

%vec = vec(:,ones(1,cols));

indices = (1:length(vec))';

indices = indices(:,ones(1,cols));

vec = vec(indices);

end

if isrowvec

% transform vector back

vec = vec';

end

elseif version == 2

% calculate indices

indices = (1:length(vec))';

% replicate rows

if rows > 1

c = indices(:,ones(rows,1));

indices = c(:);

end

% replicate columns

if cols > 1

indices = indices(:,ones(1,cols));

end

% transform index when row vector

if isrowvec

indices = indices';

end

% get vector based on indices

vec = vec(indices);

end

end

随时使用所有数据测试该功能，并给我反馈。当您发现可以改善的地方时，请告诉我。

总是更快？我们可以看到一些数字，跨平台/操作系统/矢量大小/等等吗？