matlab，FFT三维频谱图及其三维中心频谱图代码

以下是 Matlab 代码示例，用于生成三维频谱图及其三维中心频谱图：

% 生成三维频谱图及其三维中心频谱图

% 生成测试数据
x = rand(64,64,64);

% 计算三维FFT
y = fftn(x);

% 计算三维频谱
z = abs(y);

% 计算三维中心频谱
shifted_z = fftshift(z);

% 绘制三维频谱图
figure
isosurface(z,0.5)
axis equal
xlabel('X')
ylabel('Y')
zlabel('Z')
title('三维频谱图')

% 绘制三维中心频谱图
figure
isosurface(shifted_z,0.5)
axis equal
xlabel('X')
ylabel('Y')
zlabel('Z')
title('三维中心频谱图')

以上代码中，我们首先生成了一个 64x64x64 的测试数据，然后对其进行三维FFT计算，得到三维频谱，再计算三维中心频谱。最后，我们使用 `isosurface` 函数绘制了三维频谱图和三维中心频谱图。