语音信号的时域处理语音信号的声道参数分析.doc

PAGE PAGE 2 Homework2 语音信号的时域处理语音信号的声道参数分析 刘超 211117班 小班号13 学号：020501 实验内容简介： 采集一段语音信号，或从文件获得一段语音信号，采样率为8KHZ，量化精度为16比特线性码； 分析帧长30ms (240样点)； 取出一帧语音信号，注意不要取无声（静音，幅度很小）帧。 利用哈明窗对输入帧语音加窗； 求短时自相关函数。 利用Durbin算法求输入语音帧的声道参数，A参数。 利用求出的A参数构造分析滤波器A（z），利用A(z)对输入语音信号进行滤波，得到残差信号。 编写程序完成上述计算 实验步骤： 用cooledit录了一段音，本次试验用的语音文件是刘江的。采样率为8KHZ，量化精度为16比特线性码。还是上次实验用的那个语音文件 利用哈明窗对输入帧语音加窗； 求短时自相关函数。 利用Durbin算法求输入语音帧的声道参数，A参数。 利用求出的A参数构造分析滤波器A（z），利用A(z)得到h(n),对输入语音信号进行滤波，得到残差信号。 画出输入语音和残差信号的波形图。 实验结论分析： 从图中明显可以看出，arpha是递减的，和理论分析一致，随着阶数的增加，预测误差能量减小； 从递减的情况看，p取6到8就可以了。 残差信号幅度明显比原始信号的幅度小很多，而且比较平滑。方便编码，可以减小存储空间。 浊音信号的arpha比清音的下降得快。但是起点要高一点。原来讲义中间的土有点问题。 画出的波形如下图所示： 实验心得和体会：通过这次实验，我对理论的理解更加深入了，对matlab的应用也更加熟练了。这次实验的内容看似很简单，因为所有的公式都已经给出，分析的方法都有了，连流程图都已经给好了，自己要做的仅仅只是把它们转化成计算机语言。但是实验过程中间还是遇到了不少麻烦，因为matlab中的数组不是从零开始的，所以给编程和调试都带来了很大的麻烦。整个实验，是和同寝室的刘江同学一组的，通过相互检查，找错，大大提高了程序跳时的速度。 因为调试很麻烦，不知道错误出在什么地方，所以我先自己计算出m=0，1，2，3，4的所有数值，然后带入循环体中间验算。 设想一，因为c++中间的数组是从零开始的，所以用c++编写源程序，画图的时候用matlab。 设想二，自己把Durbin算法重新算一次，推导一次，把下标弄得和matlab相吻合。 附录： 调试程序用的数据：自己算的。 %%%%%%%%%%% r =0.6871 0.5487 0.2126 -0.1592 -0.3948 -0.4093 -0.2419 -0.0030 0.1801 0.2534 o = 0.6871 0.2489 0.0445 -0.4270 0.1278 -9.4294 -8.6233 -8.6094 -7.9725 -6.2655 m=0: a0(0)=1 o(0)= 0.6871 m=1: k=-0.7986 a(0)=1 a(1)=-0.7986 o(1)= 0.2489 m=2: k= 0.9064 a(0)=1 a(1)= -1.5225 a(2)= o(2)= 0.0445 m=3: k= -1/0.0445*(-0.1592*1+ 0.2126 *(-1.5225)+0.5487* 0.9064)= -0.3249 a(0)=1 a(1)= -1.8170 a(2)= 1.4011 o(3)= 0.0398 m=4: k= -1/0.0445*(-0.1592*1+ 0.2126 *(-1.5225)+0.5487* 0.9064)= -0.3249 o(4)= -0.3560 a = 1 1 1 1 1 1 1 1 1 1 1 -0.7986 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1 -1.5222 0.9061 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1 -1.8129 1.3944 -0.3208 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1 -1.6979 0.8947 0.3289 -0.3584 1.0000 1.0000 1.0000 1.0000 1.0000