《数字信号处理》课程设计lpc.ppt

《数字信号处理》课程设计———语音的线性预测编码 提 纲 课程设计目的 课程设计内容 课程设计要求 课程设计设备 课程设计原理 课程设计过程 课程设计目的 了解和熟悉语音信号线性预测分析的基本原理 掌握自相关互相关的计算方法 掌握基音周期的求法 熟悉MATLAB的编程 课程设计内容 结合语音信号处理中的线性预测编码理论、求解线性预测方程的matlab函数以及求基音周期的方法，实现对一段语音的分析（用MATLAB编写函数实现线性预测编码，提取特征）、合成（利用提取的特征，合成原来的语音信号）。 课程设计要求 使用MATLAB语言实现基于线性预测编码（LPC）的语音信号分析与合成 将合成的语音信号与原语音信号作比较，评估合成效果 课程设计设备 计算机 MATLAB编程环境 Cool Edit软件 课程设计原理 语音的发声过程可以用周期冲击脉冲序列（浊音）或随机噪声序列（清音）通过IIR滤波器来模拟。IIR滤波器的系数就是用线性预测编码器编码出来的结果。滤波器的阶数取10就可以了（可能保证语音信号的清晰）。 语音信号在20-30ms内可以近似看成稳定的，在这一段时间内，K系数可以看成是稳定不变的。 课程设计原理 课程设计原理 用于语音分析与合成的线性预测编码（LPC）方法是基于将声道按一线性全极点（IIR）滤波器建模得到的，这个滤波器具有的系统函数为 课程设计过程 y=wavread('lpc'); %读取语音信号 N=round(length(y)/space-0.5); %将语音信号分为段长为160样本的若干段 [a,g]=lpc(s,p); %利用lpc函数求解预测系数ap % 预测系数的自相关序列 for k=1:p ra(k)=a(1:p-k)*a(1+k:p)'; end %预测系数与相关系数的互相关 for k=11:space re(k)=ra(1:p)*r(k:-1:k-p+1)'; end 课程设计过程 %周期冲激串发生器 for k=0:round(space/period(i)-0.5) x(k*period(i)+1)=1*1600*G(i)*(-1)^k; end %白噪声发生器 x=randn(1,space)*G(i); %画出语音波形 subplot(211);plot(y);title('输入语音波形'); subplot(212);plot(result);title('合成后的语音波形'); %播放声音文件 wavplay(y,8000); wavplay(result,8000) 课程设计过程 课程设计过程 演示 原始语音 合成语音 原始语音编码速率为64kbit/ s 合成语音编码速率为3.6kbit/ s，压缩为原来的9/160 更详细的内容可参考：[美]维纳.K.恩格尔，约翰.G.普罗科斯著，刘树棠译.数字信号处理——使用MATLAB[M].西安：西安交通大学出版社，2002：p381-384. 谢 谢！ * * 指导教师 李红莲 白噪声 发生器 周期冲激 串发生器 全极点 滤波器 浊音与清音开关 合成语音信号 语音信号产生的方框图模型 通过使观察到的样本s(n)和预测值之间的平方误差之和最小，来确定滤波器的系数，再通过这些系数以及相关系数确定语音桢是浊音还是清音，以及浊音的基音周期 语音信号分析与合成方框图 读语音信号 分帧 求LPC系数 求音调周期 是浊音？ 浊音合成 清音合成 合成后的语音信号与原语音信号作比较 课 程 设 计 过 程 shuru 原始语音波形 合成语音波形