数字信号处理课程设计报告基于MATLAB的语音信号的特技处理.docx

1、数字信号处理课程设计报告基于MATLAB的语音信号的特技处理数字信号处理课程设计报告-基于MATLAB的语音信号的特技处理xxxx 数字信号处理课程设计报告 题目:基于MATLAB 的语音信号的特技处理 系 (院): 计算机工程学院 专 业: 通信工程 班 级: 通信xx班 学 号: xxxxxxxx 姓 名: xxx 指导教师: xxx 学年学期: 2009 2010 学年 第 1 学期 2009年12月 18 日 设计任务书 课题 基于MATLAB 的语音信号的特技处理 名称 1. 巩固所学的数字信号处理理论知识， 理解信号的采集、处理、传输、显示和存储过程; 设计 2. 综合运用专业及基

2、础知识，解决实际工程技术问题的能力; 目的 3. 学习资料的收集与整理，学会撰写课程设计报告。 1. 微型电子计算机(PC); 实验 环境 2. 安装Windows 2000以上操作系统，MATLAB等开发工具。 1. 选择一个语音信号作为分析的对象，或录制一段各人自己的语音信号，并对其进行频谱分析; 然后在时域用数字信号处理方法将信号加入延时和混响， 再分析其频谱，并与原始信号频谱进行比较;最后设计一个信号处理系统界面。 2. 利用课余时间去图书馆或上网查阅课题相关资料，深入理解课题含义及设计要求，任务 注意材料收集与整理; 要求 3. 在第15周末之前完成预设计，并请指导教师审查，通过后方

3、可进行下一步工作; 4. 结束后，及时提交设计报告(含纸质稿、电子稿)，要求格式规范、内容完整、结论正确，正文字数不少于3000字(不含代码)。 工作进度计划 序号 起止日期 工 作 内 容 2009.12.142009.12.14 在预设计的基础上，进一步查阅资料，完善设计方案。 1 2009.12.142009.12.17 设计总体方案，构建、绘制流程框图，编写代码，上机调试。 2 2009.12.172009.12.18 测试程序，完善功能，撰写设计报告。 3 2009.12.18 参加答辩，根据教师反馈意见，修改、完善设计报告。 4 指导教师(签字): 年 月 日 摘要 语音是人们交流

4、思想和进行社会活动的最基本手段,我们要对语音信号进行测定并将其转变为另一种形式，以提高我们的通信能力。 语音信号数字特技与时基压扩的概念近年来,随着计算机及大规模数字集成电路的迅速发展,语音数字信号处理得到了相应的发展。语音信号分析模拟、语音合成、语音识别等的研究已较成熟。各种声码器、声控器、语声识别系统、语声合成器等已逐渐有商品出现。语音数字特技处理是从语音信号数字处理中发展出来的一个新的分支,近年来开始逐渐为人们所重视。语音数字特技处理是一种用数字信号处理方法对语音信号进行某种变换、组合、压缩、扩展、合成、仿真的处理技术,以求得特殊的语音效果,如:语音变速重放、音调变换、语声伪造、仿真合成

5、等等。语音数字特技处理作为一种新技术在国外广泛地被应用于广播、教育、文艺、公安等各个领域。 关键词 Matlab 语音信号 特技处理 频谱分析 1 引言 . 1 2 课题综述 . 1 2.1 课题来源 . 1 2.2 预期目标 . 1 2.3 面对的问题及需解决的关键技术 . 2 3 系统分析及设计 . 2 3.1设计的总体方案图 . 2 3.2涉及的基础知识 . 3 3.2.1混响与延时 . 3 3.2.2离散傅立叶变换 . 4 3.2.3滤波器设计 . 5 3.3解决问题的基本思路 . 5 3.4 详细流程图 . 5 4 代码编写 . 6 5 程序调试 . 9 5.1 调试过程与步骤 .

6、9 5.2 发现的问题 . 9 5.3 解决的办法 . 9 6运行与测试 . 10 6.1运行程序 . 10 结 论 . 13 致 谢 . 14 参 考 文 献 . 15 数字信号处理课程设计报告 1 引言 数字信号处理是随着计算机技术的发展而迅速发展起来的一门新兴而古老的学科，它在新的领域如生物医学工程、声学、雷达、地震不、语音通信、数据通信、核科学等学科发挥着重要的作用，而它所采用的各种方法及众多应用已有悠久的历史;同时也是一门具有很强的理论性与实践性，且理论和技术发展都十分迅速的前沿性学科。随着数字化时代的来临，科学技术的进步而生产发展需求的与日俱增，促进了数字信号处理学科的发展，产生了

7、各种巧妙的信号处理算法;特别是计算机技术的飞速发展，为数字信号处理增添了巨大的生命力。 数字信号处理主要是研究用数字或符号序列表示和处理信号。处理的目的可以是削弱信号中的多余内容，滤除混杂的噪声和干扰，或者是将信号变换为容易分析和识别的形式，便于估计和选择它的特征参数。例如通过分析和运算，可以估计脑电图或心电图中的某种特征参数，帮助医生查找病因和分析病情，确定合理的治疗方案;又如，信号在传输时，要受到各种干扰，包括失真、衰落和混入的背景噪声，信号处理要排除这些干扰。 声音信号是一维连续信号,而计算机只能处理离散信号。为了从离散信号还原连续信号,根据采样定理,可以确定采样频率的最小值。wav文件

8、是一种数字声音文件格式,本课程设计基于Matlab分析了wav声音文件频谱与声音的关系。通过采集个人的一段声音进行频谱分析等处理，然后设计数字滤波器处理这个原始声音的wav文件,并比较滤波以后输出声音信号与原声音信号的异同。 2 课题综述 2.1 课题来源 近年来,随着计算机及大规模数字集成电路的迅速发展,语音数字信号处理得到了相应的发展。语音信号分析模拟、语音合成、语音识别等的研究已较成熟。、各种声码器、声控器、语声识别系统、语声合成器等已逐渐有商品出现。 2.2 预期目标 选择一个语音信号作为分析的对象，或录制一段各人自己的语音信号，并对其进行频谱分析; 然后在时域用数字信号处理方法将信号

9、加入延时和混响， 再分析其频谱，并与原始信号频谱进行比较;最后设计一个信号处理系统界面。 1 数字信号处理课程设计报告 2.3 面对的问题及需解决的关键技术 1(研究语音信号的特点; 2(探讨语音分析与识别、处理的基本理论基础、关键技术; 3(写出各步骤的MATLAB的程序代码; 4(分析采集的语音信号的时域波形与频谱;分析滤波前后语音信号的波形与频谱; 5(滤波器的设计可采用图形化设计工具FDATool和图形化信号处理设计工具SPTool; 3 系统分析及设计 3.1设计的总体方案图 利用Windows下的录音机或其他软件，录制一段自己的语音信号，时间控制在1s左右，并对录制的信号进行采样

10、语音信号的频谱分析，画出采样后语音信号的时域波形和频谱图 将信号加入延时和混响，再分析其频谱，并与原始信号频谱进行比较 几种特殊类型的滤波器:无限个回声滤波器，全通结构的混响器，并画出滤波器的频域响应。用自己设计的滤波器对采集的语音信号进行滤波 分析得到信号的频谱，画出滤波后信号的时域波形和频谱，并对滤波前后的信号进行对比，分析信号的变化 回放语音信 图3-1 语音信号的特技处理设计方案框图 2 数字信号处理课程设计报告 3.2涉及的基础知识 3.2.1混响与延时 (1)混响效果主要是用于增加音源的融合感。自然音源的延时声阵列非常密集、复杂，所以模拟混响效果的程序也复杂多变。常见参数有以下几种

11、: 混响时间:能逼真的模拟自然混响的数码混响器上都有一套复杂的程序，其上虽然有很多技术参数可调，然而对这些技术参数的调整都不会比原有的效果更为自然，尤其是混响时间。 高频滚降:此项参数用于模拟自然混响当中，空气对高频的吸收效应，以产生较为自然的混响效果。一般高频混降的可调范围为0.11.0。此值较高时，混响效果也较接近自然混响;此值较低时，混响效果则较清澈。 扩散度:此项参数可调整混响声阵密度的增长速度，其可调范围为010，其值较高时，混响效果比较丰厚、温暖;其值较低时，混响效果则较空旷、冷僻。 预延时:自然混响声阵的建立都会延迟一段时间，预延时即为模拟次效应而设置。 声阵密度:此项参数可调整

12、声阵的密度，其值较高时，混响效果较为温暖，但有明显的声染色;其值较低时，混响效果较深邃，切声染色也较弱。 频率调制:这是一项技术性的参数，因为电子混响的声阵密度比自然混响稀疏，为了使混响的声音比较平滑、连贯，需要对混响声阵列的延时时间进行调制。此项技术可以有效的消除延时声阵列的段裂声，可以增加混响声的柔和感。 调治深度:指上述调频电路的调治深度。 (2)延时就是将音源延迟一段时间后，再欲播放的效果处理。依其延迟时间的不同，可分别产生合唱、镶边、回音等效果。 当延迟时间在335ms之间时人耳感觉不到滞后音的存在，并且他与原音源叠加后，会因其相位干涉而产生梳状滤波效应，这就是镶边效果。如果延迟时间

13、在50ms以上时，其延迟音就清晰可辨，此时的处理效果才是回音。回音处理一般都是用于产生简单的混响效果。 延时、合唱、镶边、回音等效果的可调参数都差不多，具体有以下几项: *延时时间(Dly)，即主延时电路的延时时间调整。 *反馈增益(FB Gain)，即延时反馈的增益控制。 *反馈高频比(Hi Ratio)，即反馈回路上的高频衰减控制。 3 数字信号处理课程设计报告 *调制频率(Freq)，指主延时的调频周期。 *调制深度(Depth)，指上述调频电路的调制深度。 *高频增益(HF)，指高频均衡控制。 *预延时(Ini Dly)，指主延时电路预延时时间调整。 *均衡频率(EQ F)，这里的频率

14、均衡用于音色调整，此为均衡的中点频率选择。 由于延时产生的效果都比较复杂多变，如果不是效果处理专家，建议使用设备提供的预置参数，因为这些预置参数给出的处理效果一般都比较好。 3.2.2离散傅立叶变换 MATLAB的信号处理工具箱中函数FFT和IFFT用于快速傅立叶变换和逆变换。在下面介绍这些函数。 函数FFT用于序列快速傅立叶变换。 函数的一种调用格式为 y=fft(x) 其中，x是序列，y是序列的FFT，x可以为一向量或矩阵，若x为一向量，y是x的FFT。且和x相同长度。若x为一矩阵，则y是对矩阵的每一列向量进行FFT。 如果x长度是2的幂次方，函数fft执行高速基,2FFT算法;否则fft

15、执行一种混合基的离散傅立叶变换算法，计算速度较慢。 函数FFT的另一种调用格式为 y=fft(x,N) 式中，x，y意义同前，N为正整数。 函数执行N点的FFT。若x为向量且长度小于N，则函数将x补零至长度N。若向量x的长度大于N，则函数截短x使之长度为N。若x 为矩阵，按相同方法对x进行处理。 经函数fft求得的序列y一般是复序列，通常要求其幅值和相位。MATLAB提供求复数的幅值和相位函数:abs，angle，这些函数一般和FFT同时使用。 函数abs(x)用于计算复向量x的幅值，函数angle(x)用于计算复向量的相角，介于 和 之间，以弧度表示。 函数unwrap(p)用于展开弧度相位

16、角p ,当相位角绝对变化超过 时，函数把它扩展至 。 用MATLAB工具箱函数fft进行频谱分析时需注意: (1) 函数fft返回值y的数据结构对称性。 (2) 频率计算。 (3) 作FFT分析时，幅值大小与FFT选择点数有关，但不影响分析结果。 4 数字信号处理课程设计报告 3.2.3滤波器设计 单回声滤波器的系统函数: H(z)= 1 (3-1) 1，a，z,(,R)a无限个回声滤波器的系统函数: H(z)= 1 (3-2) z,(,R)1,a，z,(,R)a全通结构的混响器的系统函数: H(z)= 1 (3-3) a，z,(,R)1，a，z*(,R)a函数filter 函数filter的

17、调用格式为 y=filter(b,a,x) 该格式采用数字滤波器对数据进行滤波，既可以用于IIR滤波器，也可以用于FIR滤波器。其中向量b和a分别表示系统函数的分子、分母多项式的系数，若a,1，此时表示FIR滤波器，否则就是IIR滤波器。该函数是利用给出的向量b和a，对x中的数据进行滤波，结果放入向量y。 3.3解决问题的基本思路 (1)利用Windows下的录音机或其他软件，录制一段自己的语音信号，时间控制在1s左右，并对录制的信号进行采样。 (2)语音信号的频谱分析，画出采样后语音信号的时域波形和频谱图。 (3)将信号加入延时和混响，再分析其频谱，并与原始信号频谱进行比较。 (4)设计几种

18、特殊类型的滤波器:单回声滤波器，多重回声滤波器，无限个回声滤波器，全通结构的混响器，并画出滤波器的频域响应。 (5)用自己设计的滤波器对采集的语音信号进行滤波。 (6)分析得到信号的频谱，画出滤波后信号的时域波形和频谱，并对滤波前后的信号进行对比，分析信号的变化。 (7)回放语音信号。 3.4 详细流程图 5 数字信号处理课程设计报告 用wavread从自己的电脑导入Windows下1s语音信号 画出采样后语音信号的时域波形与频谱图。N=2001,和调用fft函数快速傅里叶变换。 对采样后语音信号进行延时30，在调整采样信号的长度，两者相加，就得到混响。 无限个回声滤波器 全通结构的混响器 单

19、回声滤波器 Bz=0,0,0,0,0,0,0,0,0,0,Bz1=a,0,0,0,0,0,0,0,0,0,1a=0.5; 1; ; y2=x+z*0.5; Y2=fft(y2,2001); Az=1,0,0,0,0,0,0,0,0,0,Az1=1,0,0,0,0,0,0,0,0,0,a-a; ; yy1=filter(Bz,Az,x); yy2=filter(Bz1,Az1,x); YY1=fft(yy1,2001); YY2=fft(yy2,2001); sound(yy1,fs,bits); sound(yy2,fs,bits); sound(y2,fs,bits); 图3-2 语音信号的

20、特技处理设计程序流程图 4 代码编写 %原始信号 x1=wavread(wlx.wav); %读取信号 figure;plot(x1);title(原始信号); sound(5*x1,40000); %对原始声音的回放 6 数字信号处理课程设计报告 %采样信号 x,fs,bits=wavread(wlx.wav,10000 40000); x=x(:,1); %只取单声道 sound(5*x,fs); %对声音的回放 n1=0:2000; N=size(x,1); figure;subplot(2,1,1);plot(x);title(采样后语音信号的时域波形); Y=fft(x,2001);

21、 subplot(2,1,2); %对信号做2001点FFT变换 plot(n1(1:1000),Y(1:1000); title(采样后语音信号的频谱图); %信号幅度 %延时的信号 x,fs,bits=wavread(wlx.wav,10000 40000); x=x(:,1); %只取单声道 n1=0:2000; z=zeros(200,1);x; figure(2); %信号的延时 subplot(2,1,1);plot(z); title(延时后的时域图); %画出延时后的信号时域图 Z=fft(z,2001); subplot(2,1,2);plot(n1(1:1000),Z(1:

22、1000); %延时后的信号频谱图 title(延时后的频谱图 ); sound(5*z,fs); %混响后信号 x,fs,bits=wavread(wlx.wav,10000 40000); %读取语音信号 x=x(:,1); %只取单声道 n1=0:2000; z=zeros(200,1);x; %对语音信号进行延时 x=x;zeros(200,1); %使语音信号与延时后信号同等长度 y1=x+z; %信号的混响 figure(3);subplot(2,1,1);plot(y1); title(混响的时域图); %混响时域图 Y1=fft(y1,2001); %对混响信号2001点的FF

23、T变换 subplot(2,1,2);plot(n1(1:1000),Y1(1:1000); %混响频谱图 title(混响的频谱图); 7 数字信号处理课程设计报告 sound(5*y1,fs); %回放混响后的信号 %单回声滤波器 x,fs,bits=wavread(wlx.wav,10000 40000); %读取语音信号 x=x(:,1); %只取单声道 n1=0:2000; a=0.5; %a取小于等于1 z=zeros(200,1);x; %对语音信号进行延时 x=x;zeros(200,1); %使语音信号与延时后信号同等长度 y2=x+z*0.5; %信号经单回声滤波 figure;subplot(2,1,1);plot(y2);title(单回声滤波器时域图); Y2=fft(y2,2001); %对单回声信号做2001点的FFT变换 subplot(2,1,2);plot(n1(1:1000),Y2(1:1000); %单回声信号频谱图 title(单回声滤波器频谱图); sound(5*y2,fs,bits); %回放经单回声滤波