MATLAB 中的 3-D 图包含：时间、频率和功率谱密度

Question

我目前正在为我的语音处理课程做一个项目，并且刚刚完成了一个时间波形图以及两个 wide/narrow 西班牙语口语词 (aire) 的频带频谱图。

接下来的部分项目如下：

Make a 3-D plot of each word signal, as a function of time, frequency and power spectral density. The analysis time step should be 20ms, and power density should be computed using a 75%-overlapped Hamming window and the FFT. Choose a viewing angle that best highlights the signal features as they change in time and frequency.

我希望有人能为我提供一些关于如何开始这一部分的指导。我开始查看 频谱图和瞬时频率 标题下的 here，但不确定如何将 PSD 添加到脚本中。

谢谢

Answer 1

我给大家举个例子

我要生成线性调频信号。

Fs = 1000;
t = 0:1/Fs:2;
y = chirp(t,100,2,300,'linear');

然后，我将定义 fft 和 hamming 的数量 window。

nfft=128;
win=hamming(nfft);

然后我将定义重叠长度，nfft 的 75%。

nOvl=nfft*0.75;

然后，我正在使用频谱图函数执行 STFT。

[s,f,t,pxx] = spectrogram(y,win,nOvl,nfft,Fs,'psd');

'y'为时间信号，'win'为汉明定义window，'nOvl'为重叠数，'nfft'为fft数，'Fs'为采样频率，'psd'使结果pxx为功率谱密度。

最后，我将使用瀑布图绘制 'pxx'。

waterfall(f,t,pxx')
xlabel('frequency(Hz)')
ylabel('time(sec)')
zlabel('PSD')

FFT 的长度，对应于 20ms，取决于信号的采样频率。

编辑：在绘制瀑布图时，我转置 pxx 以更改 t 和 f 轴。

MATLAB 中的 3-D 图包含：时间、频率和功率谱密度

3-D Plot in MATLAB Containing: Time, Frequency and Power Spectral Density

matlab

speech-recognition