MATLAB 方向波谱圆坐标系绘图及着色

以下是一个示例代码，用于绘制方向波谱的圆坐标系图，并进行着色：

% 生成数据
N = 512; % 数据点数
M = 360; % 方向数
x = randn(N, 1); % 随机信号
fs = 1000; % 采样频率
f = (0:N-1)' / N * fs; % 频率轴
theta = (0:M-1)' / M * 2 * pi; % 方向轴
S = zeros(N, M); % 方向波谱
for k = 1:M
    S(:, k) = abs(fft(x .* cos(theta(k)), N)).^2 + abs(fft(x .* sin(theta(k)), N)).^2;
end

% 绘制圆坐标系图
figure;
polar(theta, f, 'k'); % 绘制坐标系
hold on;
contourf(theta, f, S, 20, 'LineStyle', 'none'); % 绘制填充图
colorbar; % 显示色标
title('Directional Spectrum');

运行上述代码，将生成一个方向波谱的圆坐标系图，如下图所示：