资源描述：
《语音信号处理过程》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、《现代语音处理技术及应用期中考试》姓名：李红启班级：A1011学号：36一、实验原理1、汉明窗加窗给输入语音信号加窗x(m)=s(m)w(n-m)2、倒谱和复倒谱的实现倒谱主要应用在语音信号同态处理中，同态处理方法是一种设法将非线性问题转化为线性问题来进行处理方法，它能将两个信号通过乘法合成的信号，或通过卷积合成的信号分开。卷积同态系统的典范表示如图1所示，它由三部分组成：第一部分为卷积同态系统的特征系统D*[]，其输入是若干信号的卷积组合，而输出为若干信号的加法组合。利用z变换表示，卷积组合可变为乘法组合，再利用对数特性，可将乘法组合变为加法组
2、合，再进行z逆变换，输出信号仍为加法组合，这就构成了卷积同态系统的特征系统D*[]，如图1所示。图1卷积同态系统的特征系统这给出一种卷积同态系统构成方法。如果把输入序列局限于为绝对可加的情形，则输入信号的z变换具有包括单位圆在内的一个收敛域。也就是说，可以求得这个序列的傅里叶变换。在此情况下，图1中用傅里叶变换运算代替z变换运算是恰当的。卷积同态系统可表示为图2。图2复倒谱的计算由于复倒谱出现了多值性和arg[X(ejw)]的不连续性导致我们在分析时会出现很多问题，这里不一一讨论。于是我们提出了倒谱的概念。倒谱c(n)由下式给出：这里我们对x(n
3、)的X(ejw)作了取模运算，这样算出的c(n)是实数谱并且很好的解决了复倒谱出现的问题。为得到倒频谱的近似式，可对输入序列进行离散傅里叶变换，取其模的对数，再计算离散傅里叶反变换。图3倒谱的计算图2是求复倒谱的原理图，图3是求倒谱的原理图。本软件设计中就是根据图3的原理来计算一段语音信号的倒谱的。由于经过了同态处理后将声道冲击响应和激励信号变为加信号，并且，声道冲击响应处于低时部分，而激励信号处于高时部分，可以通过倒谱滤波器将声道冲击响应和激励信号分开二、实验程序：1、录音程序x=wavrecord(16000,8000,1)wavwrite(
4、x,'a.wav')plot(x)1、浊音的程序[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(
5、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=
6、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);%一点对应1HZtitle('对数幅度谱');xlabel('Frequence(Hz)');ylabel('20lg
7、Xn(ejw)
8、');subplot(4,2,3);plot(c);title('浊音段复倒谱');xla
9、bel('t/ms');ylabel('c(n)');subplot(4,2,4);plot(c2);title('浊音段倒谱');xlabel('t/ms');ylabel('c(n)');subplot(4,2,5);plot((1:length(k3))/8,k3);title('所得的冲激响应');xlabel('t/ms');ylabel('c(n)');subplot(4,2,6);plot((1:length(h3))/8,h3);title('所得的激励信号');xlabel('t/ms');ylabel('h(n)');subp
10、lot(4,2,7);plot(Y1);title('平滑的对数幅度谱');xlabel('Frequence(Hz)');subplo
