matlab fft 相位,从相位和幅度获取傅立叶变换 – Matlab

傅立叶变换F的幅度和相位定义为：

Mag = sqrt(Real(F)^2 + Imaginary(F)^2)

和

Phase = arctan(Imaginary(F)/Real(F))

我试图编写一个采用灰度图像矩阵的matlab代码，在矩阵上执行fft2()，然后计算变换的幅度和相位。然后我想计算傅里叶变换的虚部和实部。这是通过将前两个方程重新排列成：

Real = Mag/sqrt(1 + tan(Phase)^2)

和

Imaginary = Real*tan(Phase)

最后结合和逆fft2：

F = Real + i*Imaginary

image = ifft2(F)

我期望看到与输入相同的图像，但我得到垃圾。我的数学错了吗我的matlab mfile代码如下：

function y = forwardBackwardFFT(image)

F = fft2(image);

mag = sqrt(real(F).^2 + imag(F).^2);

phase = atan(imag(F)./real(F));

re = sqrt((mag.^2)./(1 + tan(phase).^2));

im = re.*tan(phase);

F = re + i*im;

f = ifft2(F);

subplot(1,2,1);

imshow(image);

Title('Original Image');

subplot(1,2,2);

imshow(f);

Title('Image after forward and backward FFT');

y = f;

非常感谢 ：)