suponha que eu tenha uma matriz como esta:

a = 
    1    2
    3    4

Eu quero dobrar o tamanho da matriz e criar algo assim:

aa = 
     1    1    2    2
     1    1    2    2
     3    3    4    4
     3    3    4    4

Dessa maneira, cada elemento da primeira matriz se propaga para quatro elementos na matriz maior.

a(i,j) == aa(2*i-1, 2*j-1)
       == aa(2*i  , 2*j-1)
       == aa(2*i-1, 2*j)
       == aa(2*i  , 2*j)

Existe alguma função pré-definida para fazer isso?

definitivamente eu posso fazer isso por dois loops, mas eu quero o jeito mais fácil e mais limpo!