matlab的pwelch函数画图设置0-2pi或者-pi-pi

在 MATLAB 中，可以使用 `pwelch` 函数绘制功率谱密度图。默认情况下，`pwelch` 函数返回的频率向量是基于采样率和 FFT 长度计算得到的，单位为 Hz。如果需要将频率轴转换为弧度单位，则可以将频率向量乘以 $2\pi$。如果需要将频率轴限制在 $[-\pi, \pi]$ 或 $[0, 2\pi]$ 范围内，则可以使用 `unwrap` 函数对相位进行展开。

以下是一个示例代码，演示了如何使用 `pwelch` 函数绘制功率谱密度图，并将频率轴转换为弧度单位，以及将频率轴限制在 $[-\pi, \pi]$ 范围内：

% 生成信号
fs = 1000; % 采样率
t = 0:1/fs:1; % 时间向量
x = sin(2*pi*50*t) + sin(2*pi*150*t); % 信号

% 对信号进行功率谱密度估计
nfft = 256; % FFT 长度
[Pxx,f] = pwelch(x, hamming(nfft), nfft/2, nfft, fs);

% 将频率向量转换为弧度单位
w = f * 2*pi;

% 绘制功率谱密度图
plot(w, 10*log10(Pxx));
xlabel('Frequency (rad)');
ylabel('Power/Frequency (dB/rad)');

% 将频率轴限制在 [-pi, pi] 范围内
xlim([-pi, pi]);
xticks([-pi, -pi/2, 0, pi/2, pi]);
xticklabels({'-\pi', '-\pi/2', '0', '\pi/2', '\pi'});

在上述示例代码中，使用 `xticks` 和 `xticklabels` 函数将 x 轴上的刻度值设置为 $[-\pi, -\pi/2, 0, \pi/2, \pi]$，并使用 `xlim` 函数将 x 轴范围限制在 $[-\pi, \pi]$ 范围内。