我目前正在为我的语音处理课程设计一个项目，并且刚刚完成了西班牙语的一个口头单词的时间波形图以及宽带/窄带频谱图的制作。
该项目的下一部分如下：
根据时间，频率和功率谱密度绘制每个字信号的3D图。分析时间步长应为20ms，并应使用75％重叠的汉明窗和FFT计算功率密度。选择随时间和频率变化而最能突出显示信号特征的视角。
我希望有人可以为我提供一些有关如何开始执行此部分的指导。我从这里的“频谱图和瞬时频率”标题下开始，但是不确定如何将PSD添加到脚本中。
谢谢
我将举一个例子。
我将生成线性线性调频信号。
Fs = 1000;
t = 0:1/Fs:2;
y = chirp(t,100,2,300,'linear');
然后，我将定义fft和hamming窗口的数量。
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”是定义的汉明窗，“ nOvl”是重叠数，“ nfft”是fft数，“ Fs”是采样频率，“ psd”使结果pxx为幂光谱密度。
最后，我将使用瀑布图来绘制“ pxx”。
waterfall(f,t,pxx')
xlabel('frequency(Hz)')
ylabel('time(sec)')
zlabel('PSD')
FFT的长度（对应于20ms）取决于信号的采样频率。
编辑：在绘制瀑布图时，我转置了pxx来更改t和f轴。
本文收集自互联网，转载请注明来源。
如有侵权，请联系[email protected] 删除。
我来说两句