课件：语音信号处理解析.ppt

*/81 8.5 短时傅里叶变换 频率分辨率Δf、取样周期T、加窗宽度N三者关系： 窗形状对短时傅立叶变换的影响 矩形窗——主瓣窄，衰减慢； 汉明窗——主瓣宽，衰减快； 窗宽对短时频谱的影响 窗宽长——频率分辨率高，能看到频谱快变化； 窗宽短——频率分辨率低，看不到频谱的快变化； */81 8.5 短时傅里叶变换 二、移动窗形状对 的影响 元音[a]的短时幅度谱（n=512） rectangular window hamming window */81 8.5 短时傅里叶变换 从短时频谱图中可知： 快速变化，由激励信号引起的； 慢速变化，声道滤波器的共振峰特性引起的； 采用汉明窗得到的短时频谱较矩形窗平滑，因而在语音分析中汉明窗用得较普遍。 hamming window */81 三、窗宽对 的影响 窗宽为64点（10kHz）元音[a]的短时频谱图 短窗的频率分辨率下降，频谱图中丢失了关于语音基音周期的信息，只保留慢变化。 矩形窗 汉明窗 8.5 短时傅里叶变换 */81 8.5 短时傅里叶变换 四、STFT的线性滤波实现 第1种： */81 8.5 短时傅里叶变换 第2种： */81 8.5 短时傅里叶变换 语音的时域分析和频域分析是语音分析的两种重要的方法，但是这两种方法均有局限性： 时域分析对语音信号的频率特性没有直观的了解； 频域特性中又没有语音信号随时间的变化关系。 因此人们致力于研究将时域分析和频域相结合，综合了频谱图和时域波形的优点，明显得展示了语音频谱随时间的变化情况。 五、语谱图 */81 8.5 短时傅里叶变换 */81 8.5 短时傅里叶变换 语谱图结构： 横轴为时间； 纵轴为频率； 强度用灰度来表示，指示特定频带的能量大小，即声音的强弱。 根据滤波器的带宽可将语谱图分为两类。带宽为300 Hz的宽带语谱图可以显示细致的时间结构，但谐波结构就不太清楚。对于一个有声带振动的声音来说，语谱图上有垂直的条纹，每一个条纹代表了声门的一次开放。而宽的横带则为共振峰。带宽为45 Hz的窄带语谱图使时间的结构模糊，但是频率的信息则较好。在较宽的共振峰带中可以看到各别的谐波频率。 */81 8.5 短时傅里叶变换 时域波形与语谱图 */81 8.5 短时傅里叶变换 共振峰与语谱图 在频域上，能量集中处就是共振峰之所在，在语谱图上就是颜色较深的位置。 在发元音时，音强较大，声带振动而呈现出基频及其谐振频率，也可以明显看到共振峰，能量集中在低频。 如果是发辅音，而且声带不振动，就看不到谐振频率。通常辅音的音强小，颜色看来就比较淡，而且能量较集中在高频。 若是在没有语音的空档，则语谱图上呈现的，就是有一段空白。 */81 使用相关函数 1. 使用tfrstft函数实现STFT 由于一般情况下MATLAB未提供该函数，需要自行安装。 下载最新的时频分析工具箱，假设为tftb； 将tftb复制到MATLAB安装目录下的toolbox目录下面 使用下列语句将该工具箱的路径添加到搜索路径下： addpath(genpath(C:\Program Files\MATLAB\R2009a\toolbox\tftb)) 在命令窗口中输入下面的命令： which tfrstft 若显示路径正确，则表示安装成功。 */81 使用相关函数 函数： [TFR,T,F]=TFRSTFT(X,T,N,H,TRACE) X：待分析信号. T：时间? ?? ? (默认值: 1:length(X))； N：频率点数 (默认值: length(X))； H：选择的平滑窗函数? ? (默认值: Hamming(N/4))； TRACE：如果不为零，跟踪计算过程? ?(default : 0)； TFR：得到的时频谱值（复数）； F：归一化频率。 */81 使用相关函数 2. 使用specgram函数： [b,f,t]=specgram(data,nfft,Fs,window,numoverlap); 目前，该函数已被spectrogram函数取代。 */81 使用相关函数 3. spectrogram函数 S?=?spectrogram(x) S?=?spectrogram(x,window) S?=?spectrogram(x,window,noverlap) S?=?spectrogram(x,window,noverlap,nfft) S?=?spectrogram(x,window,noverlap,nfft,fs) [S,F,T]?=?spectrogram(x,window,noverlap,F) [S,F,T]?=?spectrogram(x,window,nover