Code covered by the BSD License
function Sbig=periodic_sg(S0) [Mx,My]=size(S0); Sbig=[S0(Mx+1-(1:end),My+1-(1:end)) S0(Mx+1-(1:end),:) S0(Mx+1-(1:end),My+1-(1:end)); S0(:,My+1-(1:end)) S0 S0(:,My+1-(1:end)); S0(Mx+1-(1:end),My+1-(1:end)) S0(Mx+1-(1:end),:) S0(Mx+1-(1:end),My+1-(1:end))];
Contact us