我说的是海浪谱中的PM谱

非常抱歉，我理解错了你的问题。以下是绘制海浪谱中的PM谱的MATLAB代码：

% 设置参数
fs = 1; % 采样频率
N = 256; % FFT长度
M = 120; % 每帧采样点数
P = 60; % 重叠点数

% 读取海浪高度数据
data = load('waves.mat');
waves = data.waves;

% 分帧
frames = buffer(waves, M, P, 'nodelay');

% 加窗
win = hamming(M);
frames_win = frames .* repmat(win, 1, size(frames, 2));

% 计算短时功率谱
frames_power = abs(fft(frames_win, N)).^2 / N;

% 计算短时相位谱
frames_phase = angle(fft(frames_win, N));

% 计算PM谱
pm_spec = zeros(N/2+1, size(frames_power, 2));
for i = 1:size(frames_power, 2)
    pm_spec(:, i) = diff(frames_phase(:, i)) / (2*pi/fs);
end

% 绘图
f = linspace(0, fs/2, N/2+1);
t = linspace(0, size(waves, 1)/fs, size(pm_spec, 2));
imagesc(t, f, pm_spec);
axis xy;
xlabel('Time (s)');
ylabel('Frequency (Hz)');
title('PM Spectrum of Wave Height');
colorbar;

该代码也使用了STFT来计算短时功率谱和短时相位谱，然后通过相邻帧相位差的差分来计算PM谱。最后使用MATLAB的`imagesc`函数来绘制PM谱。需要将`waves.mat`替换为实际的海浪高度数据文件名。