请写一段matlab代码，要求对图像I进行二维快速傅里叶变换

在MATLAB中，你可以使用`fft2`函数来进行二维快速傅里叶变换（FFT）。假设`I`是你的灰度图像，这里是一段简单的示例代码：

% 加载或读取图像数据
I = imread('image.jpg'); % 替换为你的图像文件名

% 确保图像为双精度浮点数
I = single(I);

% 对图像进行二维快速傅里叶变换
F = fft2(I);

% 计算实部和虚部
F_mag = sqrt(real(F).^2 + imag(F).^2);
F_phase = atan2(imag(F), real(F));

% 如果需要将频率轴转换为空间轴，可以使用ifftshift和ifft2
F_shifted = ifftshift(F);
I_back = abs(ifft2(F_shifted)); % 取绝对值得到重构后的图像近似于原像

% 显示原始图像和反变换后的图像
subplot(1, 2, 1), imshow(uint8(I)), title('Original Image');
subplot(1, 2, 2), imshow(uint8(I_back)), title('Reconstructed Image (Magnitude)');

在这个例子中，`F`存储了频域信息，`F_mag`和`F_phase`分别包含了幅度和相位信息。注意`imshow`函数显示的是图像的强度值，不是复数。