Repeat element of a vector n times without loop.

106 次查看(过去 30 天)
Say I have a column vector x=[a;b;c]. I want to repeat each element n times to make a long length(x)*n vector. For example, for n=3, the answer would be:
ans=
a
a
a
b
b
b
c
c
c
Can anyone think of an elegant way to do this without looping?
Thanks,
Justin
  1 个评论
John
John 2015-12-9
U can use repmat it not exactly elegant but it will do the job
x=[a;b;c]; n=3;
newx = [repmat(x(1),n,1);repmat(x(2),n,1);repmat(x(3),n,1)]

请先登录,再进行评论。

采纳的回答

Azzi Abdelmalek
Azzi Abdelmalek 2012-8-28
编辑:Azzi Abdelmalek 2012-8-28
n=3 ; x=(1:3)' % example
r=repmat(x,1,n)';
r=r(:)'
  3 个评论
Jan
Jan 2012-8-29
I guess, you are right. repmat(1:3, 1, 2) = [1,2,3,1,2,3] but the OP wants [1,1,2,2,3,3]. Then r = repmat(1:3, 2, 1); r = r(:) avoid the expensive transposition of the matrix. Well, I admit that even reading this message will waste more time then millions of matrix transpositions will cost...

请先登录,再进行评论。

更多回答(6 个)

jack
jack 2015-11-23
I would use
repelem(X,3,1)
  3 个评论
Arif Billah
Arif Billah 2023-8-1
This should be chosen as the best 'correct' answer, thanks!
DGM
DGM 2023-8-2
编辑:DGM 2023-8-2
This is probably the more accepted answer today (hence the upvotes), but repelem() was not available until after the question was originally answered (R2015a).

请先登录,再进行评论。


Walter Roberson
Walter Roberson 2012-8-28
kron(x, ones(n,1))
  4 个评论
Abdelrahman Abdeltawab
编辑:Abdelrahman Abdeltawab 2018-12-13
Dear Walter Roberson,
why you did not use outer product and you chosen kronecker ( just curious ) because the guy's question was having vectors ?
Walter Roberson
Walter Roberson 2018-12-14
The * matrix multiplication operator cannot by itself repeat elements. You would need something like
(x.' * repmat(eye(length(x)), 1, n)).'
if you wanted to use the * operator to duplicate elements -- forcing you to call upon repmat() to duplicate elements.
Using the kronecker is a known idiom for duplicating data. It can be used for non-vectors too.
>> kron([1 2;3 4], ones(3,1))
ans =
1 2
1 2
1 2
3 4
3 4
3 4

请先登录,再进行评论。


Kevin Moerman
Kevin Moerman 2012-8-29
There is several others ways of doing it which in some cases are more efficient. Have a look at what the size of your vector is and compare the methods. Below I compare speeds and it appears that on my computer the third and fourth methods are mostly faster for large arrays.
n=100000; x=1:3;
a=zeros(n,numel(x)); b=a; c=a; d=a; %memory allocation
tic; a=repmat(x, n, 1); t1=toc; %Repmat method
tic; b=kron(x, ones(n,1)); t2=toc; %kron method
tic; c=x(ones(1,n),:); t3=toc; %indexing method
tic; d=ones(n,1)*x; t4=toc; %multiplication method
Kevin
  2 个评论
Walter Roberson
Walter Roberson 2021-9-13
format long g
n=100000; x=1:3;
a=zeros(n,numel(x)); b=a; c=a; d=a; %memory allocation
tic; a=repmat(x, n, 1); t1=toc %Repmat method
t1 =
0.000543
tic; b=kron(x, ones(n,1)); t2=toc %kron method
t2 =
0.006106
tic; c=x(ones(1,n),:); t3=toc %indexing method
t3 =
0.002276
tic; d=ones(n,1)*x; t4=toc %multiplication method
t4 =
0.001798

请先登录,再进行评论。


Justin Solomon
Justin Solomon 2012-8-28
Thanks guys, these all work perfectly!

Jianshe Feng
Jianshe Feng 2016-10-3
y = repmat(x,1,3); y = transpose(y); y = y(:);

Jianshe Feng
Jianshe Feng 2016-10-3
ind = [1;1;1;2;2;2;3;3;3]; x(ind)

类别

Help CenterFile Exchange 中查找有关 Loops and Conditional Statements 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by