完成图像的傅立叶变换，在同一幅图像中以子图的形式分别显示原始图像、实部、虚部、幅度谱、 相位谱、幅度谱重构图像、相位谱重构图像并进行相应标注，给出matlab代码

好的，以下是Matlab实现图像傅立叶变换，并以子图的形式显示的代码。

首先，读取原始图像并进行傅立叶变换。

img = imread('lena.jpg');
img_gray = rgb2gray(img);
f = fft2(img_gray);
fshift = fftshift(f);

接下来，计算出实部、虚部、幅度谱、相位谱。

real_part = real(fshift);
imag_part = imag(fshift);
magnitude_spectrum = log(abs(fshift));
phase_spectrum = angle(fshift);

然后，我们可以对幅度谱和相位谱进行反变换，得到幅度谱重构图像和相位谱重构图像。

magnitude_spectrum_recon = ifft2(exp(magnitude_spectrum));
phase_spectrum_recon = ifft2(exp(1i*phase_spectrum));

最后，我们可以将所有的图像以子图的形式显示出来。

subplot(2, 4, 1), imshow(img_gray), title('Input Image');
subplot(2, 4, 2), imshow(real_part, []), title('Real Part');
subplot(2, 4, 3), imshow(imag_part, []), title('Imaginary Part');
subplot(2, 4, 4), imshow(magnitude_spectrum, []), title('Magnitude Spectrum');
subplot(2, 4, 5), imshow(phase_spectrum, []), title('Phase Spectrum');
subplot(2, 4, 6), imshow(abs(magnitude_spectrum_recon), []), title('Magnitude Spectrum Reconstruction');
subplot(2, 4, 7), imshow(abs(phase_spectrum_recon), []), title('Phase Spectrum Reconstruction');

代码如上所示，运行后会生成一个包含所有子图的画布。

需要注意的是，Matlab中傅立叶变换的结果也是一个复数，因此需要分别计算出实部和虚部。幅度谱和相位谱是复数的模和辐角，需要使用log、abs和angle函数分别计算出来。反变换时，需要使用exp函数分别计算出幅度和相位。