MATLAB中的三维图包含：时间，频率和功率谱密度

Question

我目前正在为我的语音处理课程开发一个项目，并且刚刚完成了一个时间波形图以及一个宽/窄频谱图西班牙语口语（aire）。

项目的下一部分是如下：

使每个字信号的3-d情节，如时间，频率和功率谱密度的函数。分析时间步长应为20ms，功率密度应使用75％重叠的Hamming窗口和FFT进行计算。选择最适合突出显示信号特征随时间和频率变化的视角。

我希望有人能够为我提供一些指导，告诉我如何开始做这部分。我已经开始在频谱图和瞬时频率标题下寻找here，但不知道如何将PSD添加到脚本。

谢谢

Answer 1

我打算给你举个例子。

我要生成一个线性啁啾信号。

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

然后，我要定义fft和hamming窗口的数量。

nfft=128; 
win=hamming(nfft); 

然后我要定义重叠长度75％的nfft。

nOvl=nfft*0.75; 

然后，我使用光谱图功能执行STFT。

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

“y”为时间信号，“赢得”被定义汉明窗，“NOVL”是重叠的数目，“N FFT个”是FFT的数，“FS”是采样频率，和“PSD”使结果pxx作为功率谱密度。

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

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

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

编辑：在绘制瀑布图时，我调整了pxx来改变t和f轴。