第三时域分析精要.ppt

女声汉语拼音a的一帧信号（在采样频率为22050Hz的情况下，取20ms作为一帧），短时过零率为46。 男声汉语拼音s的一帧信号（在采样频率为22050Hz的情况下，取20ms作为一帧），短时过零率为183。 在实际应用中，短时平均过零率容易受到A/D转换是的直流偏移、50Hz交流电源的干扰以及噪声的影响。 减少这些干扰可以有两种方法：一种是采用带通滤波器消除信号中的直流和50Hz低频分量； Bandpass filte Hb w Hb(ejw) x(n) y(n) 3.3 短时平均过零率 另一种是用过门限率来修改过零率，减少随机噪声的影响。 过门限率反映了穿过正负门限的次数，如果存在随机噪声， 只要信号没有超过[-T,T]的范围，就没有有过零率产生。 T -T 3.3 短时平均过零率 短时能量、平均幅度、过零率 M 文件 function [zcr,e2,am]=STAV(vin,Fs,N,L) if nargin 3, N= 100; end if nargin 4, L= fix(N/2); end //分帧窗长N，重叠L x=buffer(vin,N,L); //短时能量 ie2=sum(x.^2); //短时平均幅度 iam=sum(abs(x)); //短时过零率 ncz=buffer(diff(sign(vin+1e-9)),N,L); izcr=sum(abs(ncz))/2/N; if nargout == 0 subplot(221);plot((1:length(ie2))*N/2/Fs,ie2); subplot(222);plot((1:length(iam))*N/2/Fs,iam); subplot(223);plot((1:length(izcr))*N/2/Fs,izcr); plen=length(izcr); subplot(224);plot((1:plen)*N/2/Fs,[ie2(1:plen)/max(ie2)+2;iam(1:plen)/max(iam)+1;izcr]); else zcr = izcr; e2 = ie2; am=iam; end 3.4 短时自相关函数 相关分析的依据： 用于确定两个信号在时域内的相似性。常用的物理量为自相关函数和互相关函数。当两个信号的互相关函数大时，则说明一个信号可能是另一个信号的时间滞后或提前；当互相关函数为0时，则两个信号完全不同。 自相关函数用于研究信号本身，如波形的同步性和周期性。 1.???? 定义： 时域离散确定信号的自相关函数定义为： 时域离散随机信号的自相关函数定义为：  若信号为一周期信号，周期为P，则 3.4 短时自相关函数 上式说明，周期信号的自相关函数也是一个同样周期的周期信号。 自相关函数具有下述性质： (1)??对称性 R（k）= R（-k） (2) 在k = 0处为最大值，即对于所有k来说，|R(k)|≤R(0) (3) 对于确定信号，值R（0）对应于能量，而对于随机信号，R（0）对应于平均功率 在上述的第（2）个性质中，如果是一个周期为NP的信号，则在取样0，± NP ，± 2NP ，…处，其自相关函数也是最大值，因此可以根据自相关函数的最大值的位置来估计周期信号的周期值。 3.4 短时自相关函数 对于语音来说，可以定义短时自相关函数为 2．实现 因为 所以有 3.4 短时自相关函数 定义 则 可以写成： 上式表明，序列经过一个冲激响应为 的数字滤波器滤波即得到短