function B=shift(A)
s=size(A);
Nmid=floor(s(2)/2);
Mmid=floor(s(1)/2);
h1=zeros(s);
Mmid=Mmid;
Nmid=Nmid;
if Mmid<1
	Mmid=1;
   end;
if Nmid<1
   Nmid=1;
end;

h1(Mmid,Nmid)=1;
if (s(1)>1) & (s(2)>1)
   B=ifft2(fft2(A).*fft2(h1));
else
   B=ifft(fft(A).*fft(h1));
end;