>> A = [1 2 3; 4 5 6];
>> B = repmat(A,2,3)


function B = repmat(A, M, N)
% 第一步：獲取矩陣的行數和列數
[m,n] = size(A); 
% 第二步：生成索引向量
mind = (1:m)'; 
nind = (1:n)';
% 第三步：通過索引向量創建索引矩陣
mind = mind(:,ones(1, M));
nind = nind(:,ones(1, N));
% 第四步：創建輸出數組
B = A(mind,nind);