matlab 短时傅里叶变换怎么得到相位谱

短时傅里叶变换（Short-Time Fourier Transform，STFT）是一种时频分析方法，将信号分成多个短时窗口，并对每个窗口进行傅里叶变换。STFT得到的是频率和时间信息，因此可以得到幅度谱和相位谱。

在MATLAB中，可以使用`spectrogram`函数进行STFT计算，并通过`angle`函数获得相位谱。示例代码如下：

% 生成示例信号
fs = 1000; % 采样率
t = 0:1/fs:1-1/fs; % 时间向量
f0 = 50; % 基频
x = sin(2*pi*f0*t); % 生成正弦信号

% 计算STFT
window = hann(256); % 窗口函数
noverlap = 128; % 重叠长度
nfft = 1024; % FFT长度
[S,F,T] = spectrogram(x, window, noverlap, nfft, fs);

% 获取相位谱
phase = angle(S);

% 显示结果
figure;
subplot(211);
plot(t, x);
xlabel('Time (s)');
ylabel('Amplitude');
title('Original Signal');
subplot(212);
surf(T, F, abs(S));
xlabel('Time (s)');
ylabel('Frequency (Hz)');
zlabel('Magnitude');
title('Spectrogram');

这段代码首先生成一个50Hz的正弦信号，然后使用`spectrogram`函数计算STFT，得到幅度谱和相位谱。最后使用`surf`函数将幅度谱绘制成三维图像。