function MatOut = rep(MatIn,REPN)
% 輸入矩陣大小
   [N_D,N_L] = size(MatIn);
% 計算
   Ind_D = rem(0:REPN(1)*N_D-1,N_D) + 1;
   Ind_L = rem(0:REPN(2)*N_L-1,N_L) + 1;
% 生成輸出矩陣
   MatOut = MatIn(Ind_D,Ind_L);