use kron
:
K = kron(X,Y)
returns the Kronecker tensor product of X and Y. The result is a large array formed by taking all possible products between the elements of X and those of Y. If X is m-by-n and Y is p-by-q, then kron(X,Y) is m*p-by-n*q.
In your case:
kron(a,[1 1])
will give you what you wanted
Some alternatives answers for your question:
reshape([a ; a],1,[])
reshape([a'*[1 1]]',1,[])