gl=abs(Y).^2; %功率谱 magif=ifft2(abs(f2));%幅度重构 pha=angle(f2);%取相位 phaif=ifft2(exp(j*pha));%相位重构 数据是x(i)，共N个点，采样频率是fsample