matlab已知频谱画相位代码

以下是一个示例代码，可以根据已知的频谱画出对应的相位图：

% 假设已知频谱 S，采样频率为 fs，采样点数为 N
% 频率向量
f = (0:N-1) / N * fs;
% 计算相位
phase = angle(S);
% 绘制相位图
plot(f, phase);
xlabel('Frequency (Hz)');
ylabel('Phase (rad)');
title('Phase Spectrum');

在这个示例中，我们首先定义了频率向量 `f`，然后通过 `angle` 函数计算出频谱 `S` 对应的相位，最后使用 `plot` 函数绘制相位图。