语音信号处理过程.docVIP

《现代语音处理技术及应用期中考试》 姓名： 李红启 班级： A1011 学号： 36 一、实验原理 汉明窗加窗 给输入语音信号加窗 x(m)=s(m)w(n-m) 2、倒谱和复倒谱的实现 倒谱主要应用在语音信号同态处理中，同态处理方法是一种设法将非线性问题转化为线性问题来进行处理方法，它 能将两个信号通过乘法合成的信号，或通过卷积合成的信号分开。卷积同态系统的典范表示如图1 所示，它由三部分组成：第一部分为卷积同态系统的特征系统D*[ ]，其输入是若干信号的卷积组合，而输出为若干信号的加法组合。 利用z 变换表示，卷积组合可变为乘法组合，再利用对数特性，可将乘法组合变为加法组合，再进行z 逆变换，输出信号仍为加法组合，这就构成了卷积同态系统的特征系统D*[ ]，如图1 所示。 图 1 卷积同态系统的特征系统 这给出一种卷积同态系统构成方法。如果把输入序列局限于为绝对可加的情形，则输入信号的 z 变换具有包括单位圆在内的一个收敛域。也就是说，可以求得这个序列的傅里叶变换。在此情况下，图1中用傅里叶变换运算代替z 变换运算是恰当的。卷积同态系统可表示为图 2。 图 2 复倒谱的计算 由于复倒谱出现了多值性和arg[X(ejw )]的不连续性导致我们在分析时会出现很多问题，这里不一一讨论。于是我们提出了倒谱的概念。倒谱 c(n)由下式给出： 这里我们对 x(n)的 X(ejw )作了取模运算，这样算出的 c(n)是实数谱并且很好的解决了复倒谱出现的问题。为得到倒频谱的近似式， 可对输入序列进行离散傅里叶变换，取其模的对数，再计算离散傅里叶反变换。 图 3 倒谱的计算 图 2 是求复倒谱的原理图，图 3 是求倒谱的原理图。本软件设计中就是根据图 3 的原理来计算一段语音信号的倒谱的。 由于经过了同态处理后将声道冲击响应和激励信号变为加信号，并且，声道冲击响应处于低时部分，而激励信号处于高时部分，可以通过倒谱滤波器将声道冲击响应和激励信号分开 实验程序： 录音程序 x=wavrecord(16000,8000,1) wavwrite(x,a.wav) plot(x) 浊音的程序 [y,Fs,bits]=wavread(a.wav);% 读出信号，以及采样率和采样位数 y=y(:,1); % 单声道处理 y1=y(6000:5400); % 取浊音 y2=y1.*hamming(length(y1)); % 加汉明窗 sigLength=length(y2); % 求y2的长度 Y=zeros(Fs,1); Y=20*log10(abs(fft(y2,sigLength))); % 求浊音的频谱 y=y; %转置成向量，便于处理 a=fft(y1); %对y1做DFT变换 d=log(a); %取对数 c=ifft(d); %做IDFT变换 c=c(1:60); a1=fft(y1); %对y1做DFT变换 b1=abs(a1); %取模远算 d1=log(b1);%取对数 c1=ifft(d1); %做IDFT变换 c2=c1(1:26); k=c1(1:26); k1=fft(k); k2=exp(k1); k3=ifft(k2); h=c1(27：200); h1=fft(h); h2=exp(h1); h3=ifft(h2); sigLength=length(y2); Y1=20*log10(abs(fft(k3,sigLength))); sigLength=length(y2); Y2=20*log10(abs(fft(h3,sigLength))); figure(1); subplot(4,2,1); plot((1:length(y1))/8,y2); % 以时间为横轴 title(加汉明窗的信号幅度); xlabel(t/ms);ylabel(x(n)); subplot(4,2,2); plot(Y); % 一点对应1HZ title(对数幅度谱); xlabel(Frequence(Hz)); ylabel(20lg|Xn(ejw)|); subplot(4,2,3); plot(c); title(浊音段复倒谱); xlabel(t/ms);ylabel(c(n)); subplot(4,2,4); plot(c2); title(浊音段倒谱); xlabel(t/ms);ylab