hk = ifftshift(ifft(HK));