实验三语音信号采集实验
学号：0900220501 姓名：陈娟
一实验目的
（1）了解模拟信号（语音信号）到数字信号的转换原理；
（2）熟悉DSP对数据的保存方法；
（3）熟悉CCS对数据的分析功能；
二实验内容
本实验要求使用AD转换芯片将模拟信号转换成为数字信号，使用DSP对转换后的数字信号读取并保存，并利用CCS对这些采集到的数据进行分析。
三实验原理
DSP的应用系统一般包括模数转换器（AD）、信号处理芯片（DSP）、数模转换器（DA）等主要器件，DSP系统首先将模拟信号经过一个或者多个硬件滤波器，或者其它的信号预处理，到达AD转换成为数字信号，传输到DSP，DSP对子这个信号进行采集、处理、分析，如果有必要再经过DA，转换成为模拟信号输出，实验中可以由示波器查看输出的信号波形。
DSP系统中，首要的就是数据采集。实验中，我们选择了语音信号进行采样，然后将采样后的数据从DA输出，从DA输出的通过耳机还原成为话音，也可以采集确定的信号，通过示波器观看DA转换后的波形，或者使用CCS的画图功能查看采集到的数据。
在DSP芯片中，数据的存储格式都是二进制的格式，实际中从AD读取的数据或者输出到DA的数据也必须市二进制的数据格式，而DSP仿真软件可以显示多种的数据格式，分别是16进制数、16位无符号数、二进制数、32位无符号数、32位有符号数、浮点数等多种类型。不同的数据类型有不同的显示结果，所以在查看采集到的数据的时候，一定需要明白采集到的数据格式，必须保持一致。一般从AD读取的数据都是二进制格式，但是位数不一定，从8位到32位不等，根据具体选用的AD所决定。
四实验步骤
本实验从AD采集数据和硬件系统联系较大，不同的硬件系统有不同的设置方法，这些设置方法都会在硬件系统的说明书中有详细的说明。例如，在TI的DSK实验系统上采用的AD芯片型号是TLC320AD50，有关AD50的设置可以参考该器件的PDF文档资料。
AD50采用过采样（over sample）的∑-Δ技术提供从数字到模拟和模拟到数字的高分辨率低速信号转换，该器件包括两个串行的同步转换通道，在DAC之前有一个插入滤波器和DAC之后有一个抽取滤波器，其它高级功能有片内时序和控制。∑-Δ结构在低系统速度和低价格下产生高分辨率的模数和数模转换。
对参考程序的说明
本程序从AD50采集数据然后从AD50将数据输出（AD50将AD和DA集成在同一芯片内）。对数据的采集和发送通过缓冲串口0实现。在程序运行上，需要将语音信号（或者确定信号）送到DSK板的音频信号（麦克风）输入端，然后将耳机（或者示波器）接到DSK 的信号输出端。在输入固定频率信号后运行程序，可以看到信号经过DSP复制，从示波器看到复制后的信号的波形，与信号发生器的波形一致。其中对DSP,串口0，AD50芯片的初始化都是在头文件里面具体定义的。由中断服务程序完成数据的采集和转换发送的任务.
五．任务
依照上例做语音采集，然后对采集到的数据进行高通滤波。
滤波器输入输出关系为
y[n]=-6.76195(x[n]+x[n-2])+13.456335x[n-1]
输入为（音频文件dualfre.wav）
x[n]={cos(0.1n)+cos(0.4n)}
输出为
y[n]=cos(0.4n)
主程序：
#include "Initial.h"
interrupt void BSP0_int() // Interrupt function
{
static int Data,data1[3];
int i;
for(i=2;i>0;i--)
data1[i]=data1[i-1];
data1[0]=DRR10;
Data=-6.76195*data1[0]+13.456335*data1[1]-6.76195*data1[2];
DXR10=Data&0xFFFF
}
void main()
{
asm("\t ssbx intm"); // Interrupt disable
InitDSP();
InitMcBSP0();
InitAD50();
asm("\t rsbx intm"); // Interrupt enable
while(1)
{
asm("\t nop"); // Waiting for interrupt.
asm("\t nop");
}
}
实验总结：
通过本次实验，我了解模拟信号（语音信号）到数字信号的转换原理，熟悉DSP对数据的保存方法以及对CCS对数据的分析功能。
语音信号处理实验 班级： 学号： 姓名：
实验一 基于MATLAB 的语音信号时域特征分析（2学时） 1） 短时能量 （1）加矩形窗 a=wavread('mike.wav'); a=a(:,1); subplot(6,1,1),plot(a); N=32; for i=2:6 h=linspace(1,1,2.^(i-2)*N);%形成一个矩形窗，长度为2.^(i-2)*N En=conv(h,a.*a);% 求短时能量函数En subplot(6,1,i),plot(En); if (i==2) ,legend('N=32'); elseif (i==3), legend('N=64'); elseif (i==4) ,legend('N=128'); elseif (i==5) ,legend('N=256'); elseif (i==6) ,legend('N=512'); end end 00.51 1.52 2.5 3 x 10 4 -1 100.5 1 1.5 2 2.5 3x 10 4 024 N=3200.5 1 1.5 2 2.5 3x 10 4 05 N=6400.5 1 1.5 2 2.5 3x 10 4 0510 N=12800.5 1 1.5 2 2.5 3x 10 4 01020 N=2560 0.5 1 1.5 2 2.5 3x 10 4 02040 N=512 （2）加汉明窗 a=wavread('mike.wav'); a=a(:,1); subplot(6,1,1),plot(a); N=32;
for i=2:6 h=hanning(2.^(i-2)*N);%形成一个汉明窗，长度为2.^(i-2)*N En=conv(h,a.*a);% 求短时能量函数En subplot(6,1,i),plot(En); if (i==2), legend('N=32'); elseif (i==3), legend('N=64'); elseif (i==4) ,legend('N=128'); elseif (i==5) ,legend('N=256'); elseif (i==6) ,legend('N=512'); end end 00.51 1.52 2.5 3 x 10 4 -1 100.5 1 1.5 2 2.5 3x 10 4 012 N=3200.5 1 1.5 2 2.5 3x 10 4 024 N=6400.5 1 1.5 2 2.5 3x 10 4 024 N=12800.5 1 1.5 2 2.5 3x 10 4 0510 N=2560 0.5 1 1.5 2 2.5 3x 10 4 01020 N=512 2） 短时平均过零率 a=wavread('mike.wav'); a=a(:,1); n=length(a); N=320; subplot(3,1,1),plot(a); h=linspace(1,1,N); En=conv(h,a.*a); %求卷积得其短时能量函数En subplot(3,1,2),plot(En); for i=1:n-1 if a(i)>=0 b(i)= 1;
语音识别系统实验报告 专业班级：信息安全 学号： 姓名：
目录 一、设计任务及要求 (1) 二、语音识别的简单介绍 2.1语者识别的概念 (2) 2.2特征参数的提取 (3) 2.3用矢量量化聚类法生成码本 (3) 2.4VQ的说话人识别 (4) 三、算法程序分析 3.1函数关系 (4) 3.2代码说明 (5) 3.2.1函数mfcc (5) 3.2.2函数disteu (5) 3.2.3函数vqlbg (6)
3.2.4函数test (6) 3.2.5函数testDB (7) 3.2.6 函数train (8) 3.2.7函数melfb (8) 四、演示分析 (9) 五、心得体会 (11) 附：GUI程序代码 (12) 一、设计任务及要求 实现语音识别功能。 二、语音识别的简单介绍
基于VQ的说话人识别系统，矢量量化起着双重作用。在训练阶段，把每一个说话者所提取的特征参数进行分类，产生不同码字所组成的码本。在识别(匹配)阶段，我们用VQ方法计算平均失真测度(本系统在计算距离d时，采用欧氏距离测度)，从而判断说话人是谁。 语音识别系统结构框图如图1所示。 图1 语音识别系统结构框图 2.1语者识别的概念 语者识别就是根据说话人的语音信号来判别说话人的身份。语音是人的自然属性之一，由于说话人发音器官的生理差异以及后天形成的行为差异，每个人的语音都带有强烈的个人色彩，这就使得通过分析语音信号来识别说话人成为可能。用语音来鉴别说话人的身份有着许多独特的优点，如语音是人的固有的特征，不会丢失或遗忘；语音信号的采集方便，系统设备成本低；利用电话网络还可实现远程客户服务等。因此，近几年来，说话人识别越来越多的受到人们的重视。与其他生物识别技术如指纹识别、手形识别等相比较，说话人识别不仅使用方便，而且属于非接触性，容易被用户接受，并且在已有的各种生物特征识别技术中，
电子工程系 信号与系统课程实验报告 2011-----2012学年第一学期 专业: 电子信息工程技术班级: 学号 : 姓名: 指导教师: 实常用连续时间信号的实现
一、实验目的 （1）了解连续时间信号的特点； （2）掌握连续时间信号表示的向量法和符号法； （3）熟悉MATLAB Plot函数等的应用。 二、实验原理 1、信号的定义 信号是随时间变化的物理量。信号的本质是时间的函数。 2、信号的描述 1）时域法 时域法是将信号表示成时间的函数f(t)来对信号进行描述的方法。信号的时间特性指的是信号的波形出现的先后，持续时间的长短，随时间变化的快慢和大小，周期的长短等。 2）频域（变换域）法 频域法是通过正交变换，将信号表示成其他变量的函数来对信号进行描述的方法。一般常用的是傅立叶变换。信号的频域特性包括频带的宽窄、频谱的分布等。 信号的频域特性与时域特性之间有着密切的关系。 3、信号的分类 按照特性的不同，信号有着不同的分类方法。 （1）确定性信号：可以用一个确定的时间函数来表示的信号。 随机信号：不可以用一个确定的时间函数来表示，只能用统计特性加以描述的信号。 （2）连续信号：除若干不连续的时间点外，每个时间点在t上都有对应的数值信号。离散信号：只在某些不连续的点上有数值，其他时间点上信号没有定义的信号。 （3）周期信号：存在T，使得等式f(t+T)=f(t)对于任意时间t都成立的信号。非周期信号：不存在使得等式f(t+T)=f(t)对于任意时间t都成立的信号。 绝对的周期信号是不存在的，一般只要在很长时间内慢走周期性就可以了。 （4）能量信号：总能量有限的信号。 功率信号:平均功率有限切非零的信号。 （5）奇信号：满足等式f(t)=--f(--t)的信号。偶信号：满足等式f(t)=f(--t)的信号。 三、涉及的MATLAB函数 1、plot函数 功能：在X轴和Y轴方向都按线性比例绘制二维图形。 调用格式： Plot（x,y）:绘出相x对y的函数线性图。 Plot（x1,y1,x2,y2,…..）:会出多组x对y的线性曲线图。 2、ezplot函数 功能：绘制符号函数在一定范围内的二维图形。简易绘制函数曲线。 调用格式： Ezplot (fun):在[-2π，2π]区间内绘制函数。 Ezplot (fun，［min，max］)：在［min,max］区间内绘函数。 Ezplot (funx,funy):定义同一曲面的函数，默认的区间是[0, 2π]。】 3、sym函数 功能：定义信号为符号的变量。 调用格式：sym(fun):fun为所要定义的表达式。 4、subplot函数
实验一：信号的时域分析 一、实验目的 1.观察常用信号的波形特点及产生方法 2.学会使用示波器对常用波形参数的测量 二、实验仪器 1.信号与系统试验箱一台（型号ZH5004） 2.40MHz双踪示波器一台 3.DDS信号源一台 三、实验原理 对于一个系统特性的研究，其中重要的一个方面是研究它的输入输出关系，即在一特定的输入信号下，系统对应的输出响应信号。因而对信号的研究是对系统研究的出发点，是对系统特性观察的基本手段与方法。在本实验中，将对常用信号和特性进行分析、研究。 信号可以表示为一个或多个变量的函数，在这里仅对一维信号进行研究，自变量为时间。常用信号有：指数信号、正弦信号、指数衰减正弦信号、复指数信号、Sa(t)信号、钟形信号、脉冲信号等。 1、信号：指数信号可表示为f(t)=Ke at。对于不同的a取值，其波形表现为不同的形式，如下图所示： 图1―1 指数信号 2、信号：其表达式为f（t）=Ksin（ωt＋θ），其信号的参数：振幅K、角频率ω、与初始相位θ。其波形如下图所示：
图1－2 正弦信号 3、指数衰减正弦信号：其表达式为其波形如下图： 图1－3 指数衰减正弦信号 4、Sa(t)信号：其表达式为：。Sa(t)是一个偶函数，t= ±π,±2π,…,±nπ时，函数值为零。该函数在很多应用场合具有独特的运用。其信号如下图所示：
图1－4 Sa（t）信号 5、钟形信号（高斯函数）：其表达式为：其信号如下图所示： 图1－5 钟形信号 6、脉冲信号：其表达式为f(t)=u(t)-u(t-T)，其中u(t)为单位阶跃函数。其信号如下图所示： 7、方波信号：信号为周期为T，前T/2期间信号为正电平信号，后T/2期间信号为负电平信号，其信号如下图所示 U(t)
电 子 科 技 大 学 信号与系统实验报告 学生姓名： 杜杰 学 号： 2014030103007 指导教师：张鹰 一、实验室名称：信号与系统实验室 二、实验项目名称：连续系统的幅频特性测量 三、实验原理： 正弦波信号)cos()(0t A t x ω=输入连续LTI 系统，输出)(t y 仍为正弦波信号。 图信号输入连续LTI 系统 图中， )(cos()()(000ωωωj H t j H A t y ∠+=) 通过测量输入)(t x 、输出)(t y 的正弦波信号幅度，计算输入、输出的正弦波信号幅度比值，可以得到系统的幅频特性在0ω处的测量值)(0ωj H 。改变0ω可以测出不同频率处的系统幅频特性。 四、实验目的： 使学生对系统的频率特性有深入了解。 五、实验内容： 实验内容（一）、低通滤波器的幅频特性测量 实验内容（二）、带通滤波器的幅频特性测量 六、实验器材（设备、元器件）： 数字信号处理实验箱、信号与系统实验板的低通滤波器模块U11、高通滤波器模块U21、PC 机端信号与系统实验软件、＋5V 电源 七、实验步骤：打开PC 机端软件SSP.EXE ，在下拉菜单“实验选择”中选择 “实验三”；使用串口电缆连接计算机串口和实验箱串口，打开实验箱电源。 实验内容（一）、低通滤波器的幅频特性测量 (x ) (t y
实验步骤： 1、信号选择：按实验箱键盘“3”选择“正弦波”，再按“＋”或“－”依次 选择一个频率。 2、连接接口区的“输入信号1”和“输出信号”，如下图所示。点击SSP软 件界面上的按钮，观察输入正弦波。将正弦波频率值和幅度值(Vpp/2, Vpp为峰-峰值)记录于表1。 观察输入正弦波的连线示意图 3、按下图的模块连线示意图连接各模块。 实验三实验内容（一）模块连线示意图 4、点击SSP软件界面上的按钮，观察输入正弦波通过连续系统的 响应波形；适当调整X、Y轴的分辨率可得到如下图所示的实验结果。将输出正弦波的幅度值(Vpp/2, Vpp为峰-峰值)记录于表1中。 输入正弦波和响应波形
开放实验项目报告 项目名称：语音识别机器人 专业 学生姓名 班级学号 指导教师 指导单位 2012/2013学年第一学期 一．设计背景
在科学日新月异的今天，电子设备的便捷化，人性化，智能化已成为不可逆转的潮流，而语音控制智能，更是其中研究发展的热点。凌阳SPCE061以其便捷的操作，可靠的性能，成为了各位电子爱好者的首选。本实验采用凌阳61板和运动小车（迷你型）模组设计的语音控制小车。凌阳板嵌入小车模型顶部。语音处理技术不仅包括语音的录制和播放，还涉及语音的压缩编码和解码、语音的识别等各种处理技术。本设计的语音控制小车，借助于SPCE061A在语音处理方面的特色，不仅具有前进、后退、左转、右转、停止等基本程序控制功能，而且还具备语音控制功能。 二．总流程图
三．主要模块 1、凌阳SPCE061是继μ’nSP?系列产品SPCE500A等之后凌阳科技推出的又一款16 位结构的微控制器。与SPCE500A不同的是，在存储器资源方面考虑到用户的较少资源的需求以及便于程序调试等功能，SPCE061A里只内嵌32K字的闪存（FLASH ）。较高的处理速度使μ’nSP?能够非常容易地、快速地处理复杂的数字信号。因此，与SPCE500A相比，以μ’nSP?为核心的SPCE061A 微控制器是适用于数字语音识别应用领域产品的一种最经济的选择。 其性能如下： A、16 位μ’nSP?微处理器； B、工作电压(CPU) VDD 为2.4~3.6V (I/O) VDDH 为2.4~5.5V C、CPU 时钟：0.32MHz~49.152MHz ； D、内置2K 字SRAM； E、内置32K FLASH； F、可编程音频处理； G、晶体振荡器; H、系统处于备用状态下(时钟处于停止状态)，耗电仅为2μA@3.6V ； I、2 个16 位可编程定时器/计数器(可自动预置初始计数值)； J、2 个10 位DAC(数-模转换)输出通道； K、32 位通用可编程输入/输出端口； L、14 个中断源可来自定时器A / B ，时基，2 个外部时钟源输入，键唤醒；
通信与信息工程学院 信息处理综合实验报告 班级：电子信息工程1502班 指导教师： 设计时间：2018/10/22-2018/11/23 评语： 通信与信息工程学院 二〇一八年 实验题目：语音信号分析与处理 一、实验内容 1. 设计内容 利用MATLAB对采集的原始语音信号及加入人为干扰后的信号进行频谱分析，使用窗函数法设计滤波器滤除噪声、并恢复信号。 2．设计任务与要求 1. 基本部分
（1）录制语音信号并对其进行采样；画出采样后语音信号的时域波形和频谱图。 （2）对所录制的语音信号加入干扰噪声，并对加入噪声的信号进行频谱分析；画出加噪后信号的时域波形和频谱图。 （3）分别利用矩形窗、三角形窗、Hanning窗、Hamming窗及Blackman 窗几种函数设计数字滤波器滤除噪声，并画出各种函数所设计的滤波器的频率响应。 （4）画出使用几种滤波器滤波后信号时域波形和频谱，对滤波前后的信号、几种滤波器滤波后的信号进行对比，分析信号处理前后及使用不同滤波器的变化；回放语音信号。 2. 提高部分 （5）录制一段音乐信号并对其进行采样；画出采样后语音信号的时域波形和频谱图。 （6）利用MATLAB产生一个不同于以上频段的信号；画出信号频谱图。 （7）将上述两段信号叠加，并加入干扰噪声，尝试多次逐渐加大噪声功率，对加入噪声的信号进行频谱分析；画出加噪后信号的时域波形和频谱图。 （8）选用一种合适的窗函数设计数字滤波器，画出滤波后音乐信号时域波形和频谱，对滤波前后的信号进行对比，回放音乐信号。 二、实验原理 1.设计原理分析 本设计主要是对语音信号的时频进行分析，并对语音信号加噪后设计滤波器对其进行滤波处理，对语音信号加噪声前后的频谱进行比较分析，对合成语音信号滤波前后进行频谱的分析比较。 首先用PC机WINDOWS下的录音机录制一段语音信号，并保存入MATLAB软件的根目录下，再运行MATLAB仿真软件把录制好的语音信号用audioread函数加载入MATLAB仿真软件的工作环境中，输入命令对语音信号进行时域，频谱变换。 对该段合成的语音信号，分别用矩形窗、三角形窗、Hanning窗、Hamming窗及Blackman窗几种函数在MATLAB中设计滤波器对其进行滤波处理，滤波后用命令可以绘制出其频谱图，回放语音信号。对原始语音信号、合成的语音信号和经过滤波器处理的语音信号进行频谱的比较分析。 2.语音信号的时域频域分析 在Matlab软件平台下可以利用函数audioread对语音信号进行采样,得到了声音数据变量y,同时把y的采样频率Fs=44100Hz放进了MATALB的工作空间。
语音信号的数字滤波 一、实验目的： 1、掌握使用FFT进行信号谱分析的方法 2、设计数字滤波器对指定的语音信号进行滤波处理 二、实验内容 设计数字滤波器滤除语音信号中的干扰（4 学时） 1、使用Matlab的fft函数对语音信号进行频谱分析，找出干扰信号的频谱； 2、设计数字滤波器滤除语音信号中的干扰分量，并进行播放对比。 三、实验原理 通过观察原语音信号的频谱，幅值特别大的地方即为噪声频谱分量，根据对称性，发现有四个频率的正弦波干扰，将它们分别滤掉即可。采用梳状滤波器，经过计算可知，梳状滤波器h[n]={1，A，1}的频响|H(w)|=|A+2cos(w)|，由需要滤掉的频率分量的频响w,即可得到A，进而得到滤波器的系统函数h[n]。而由于是在离散频域内进行滤波，所以令w=(2k*pi/N)即可。 对原信号和四次滤波后的信号分别进行FFT变换，可以得到它们的幅度相应。最后，将四次滤波后的声音信号输出。 四、matlab代码 clc;clear;close all; [audio_data,fs]=wavread('SunshineSquare.wav'); %读取未处理声音 sound(audio_data,fs); N = length(audio_data); K = 0:2/N:2*(N-1)/N; %K为频率采样点
%sound(audio_data,fs); %进行一次FFT变换 FFT_audio_data=fft(audio_data); mag_FFT_audio_data = abs(FFT_audio_data); %画图 figure(1) %原信号时域 subplot(2,1,1);plot(audio_data);grid; title('未滤波时原信号时域');xlabel('以1/fs为单位的时间');ylabel('采样值'); %FFT幅度相位 subplot(2,1,2);plot(K,mag_FFT_audio_data);grid; title('原信号幅度');xlabel('以pi为单位的频率');ylabel('幅度'); %构造h[n]={1，A，1}的梳状滤波器，计算A=2cosW，妻子W为要滤掉的频率%由原信号频谱可知要分四次滤波，滤掉频响中幅度大的频率分量 %第一次滤波 a = [1,0,0,0];%y[n]的系数 [temp,k]=max(FFT_audio_data); A1=-2*cos(2*pi*k/N); h1=[1,A1,1]; audio_data_h1 = filter(h1,a,audio_data); FFT_audio_data_h1=fft(audio_data_h1);
实验一 语音信号的时域分析 一、 实验目的、要求 （1）掌握语音信号采集的方法 （2）掌握一种语音信号基音周期提取方法 （3）掌握语音信号短时能量和短时过零率计算方法 （4）了解Matlab 的编程方法 二、 实验原理 语音是一时变的、非平稳的随机过程，但由于一段时间内(10-30ms)人的声带和声道形状的相对稳定性，可认为其特征是不变的，因而语音的短时谱具有相对稳定性。在语音分析中可以利用短时谱的这种平稳性，将语音信号分帧。 10～30ms 相对平稳，分析帧长一般为20ms 。 语音信号的分帧是通过可移动的有限长度窗口进行加权的方法来实现的。几种典型的窗函数有：矩形窗、汉明窗、哈宁窗、布莱克曼窗。 语音信号的能量分析是基于语音信号能量随时间有相当大的变化，特别是清音段的能量一般比浊音段的小得多。定义短时平均能量 [][]∑∑+-=∞-∞=-=-= n N n m m n m n w m x m n w m x E 122)()()()( 下图说明了短时能量序列的计算方法，其中窗口采用的是直角窗。 过零就是信号通过零值。对于连续语音信号，可以考察其时域波形通过时间轴的情况。而对于离散时间信号，如果相邻的取样值改变符号则称为过零。由此可以计算过零数，过零数就是样本改变符号的次数。单位时间内的过零数称为平
均过零数。 语音信号x (n )的短时平均过零数定义为 ()[]()[]()()[]()[]() n w n x n x m n w m x m x Z m n *--=---= ∑∞ -∞=1sgn sgn 1sgn sgn 式中，[]?sgn 是符号函数，即 ()[]()()()()???<-≥=01 01sgn n x n x n x 短时平均过零数可应用于语音信号分析中。发浊音时，尽管声道有若干个共振峰，但由于声门波引起了谱的高频跌落，所以其语音能量约集中干3kHz 以下。而发清音时．多数能量出现在较高频率上。既然高频率意味着高的平均过零数，低频率意味着低的平均过零数，那么可以认为浊音时具有较低的平均过零数，而清音时具有较高的平均过零数。然而这种高低仅是相对而言，没有精确的数值关系。 短时平均过零的作用 1.区分清/浊音： 浊音平均过零率低，集中在低频端； 清音平均过零率高，集中在高频端。 2.从背景噪声中找出是否有语音，以及语音的起点。 基音是发浊音时声带震动所引起的周期性，而基音周期是指声带震动频率的倒数。基音周期是语音信号的重要的参数之一，它描述语音激励源的一个重要特征，基音周期信息在多个领域有着广泛的应用，如语音识别、说话人识别、语音分析与综合以及低码率语音编码，发音系统疾病诊断、听觉残障者的语音指导等。因为汉语是一种有调语言，基音的变化模式称为声调，它携带着非常重要的具有辨意作用的信息，有区别意义的功能，所以，基音的提取和估计对汉语更是一个十分重要的问题。 由于人的声道的易变性及其声道持征的因人而异，而基音周期的范围又很宽，而同—个人在不同情态下发音的基音周期也不同，加之基音周期还受到单词发音音调的影响，因而基音周期的精确检测实际上是一件比较困难的事情。基音提取的主要困难反映在：①声门激励信号并不是一个完全周期的序列，在语音的
实用 中南大学 信息科学与工程学院 语音信号处理 实验报告 指导老师：覃爱娜 学生班级：信息0704 学生名称：阮光武 学生学好：0903070430 提交日期：2010年6月18日
实验一 语音波形文件的分析和读取 一、实验的任务、性质与目的 本实验是选修《语音信号处理》课的电子信息类专业学生的基础实验。通过实验： （1）掌握语音信号的基本特性理论：随机性，时变特性，短时平稳性，相关性等； （2）掌握语音信号的录入方式和*.WAV音波文件的存储结构； （3）使学生初步掌握语音信号处理的一般实验方法。 二、实验原理和步骤： WAV文件格式简介 WAV文件是多媒体中使用了声波文件的格式之一，它是以RIFF格式为标准。每个WAV文件的头四个字节就是“RIFF”。WAV文件由文件头和数据体两大部分组成，其中文件头又分为RIFF/WAV文件标识段和声音数据格式说明段两部分。常见的WAV声音文件有两种，分别对应于单声道（11.025KHz采样率、8Bit的采样值）和双声道（44.1KHz采样率、16Bit的采样值）。采样率是指声音信号在“模拟→数字”转换过程中，单位时间内采样的次数；采样值是指每一次采样周期内声音模拟信号的积分值。对于单声道声音文件，采样数据为8位的短整数（short int 00H-FFH）；而对于双声道立体声声音文件，每次采样数据为一个16位的整数（int），高八位和低八位分别代表左右两个声道。WAV文件数据块包含以脉冲编码调制（PCM）格式表示的样本。在单声道WAV文件中，道0代表左声道，声道1代表右声道；在多声道WAV文件中，样本是交替出现的。WAV文件的格式见表1。
中南大学 信号与系统试验报告 姓名: 学号: 专业班级:自动化 实验一 基本信号的生成 1．实验目的 ● 学会使用MATLAB 产生各种常见的连续时间信号与离散时间信号； ● 通过MATLAB 中的绘图工具对产生的信号进行观察，加深对常用信号的 理解； ● 熟悉MATLAB 的基本操作，以及一些基本函数的使用，为以后的实验奠 定基础。 2．实验内容 ⑴ 运行以上九个例子程序，掌握一些常用基本信号的特点及其MATLAB 实现方法；改变有关参数，进一步观察信号波形的变化。 ⑵ 在 k [10:10]=- 范围内产生并画出以下信号： a) 1f [k][k]δ=； b) 2f [k][k+2]δ=； c) 3f [k][k-4]δ=； d) 4f [k]2[k+2][k-4]δδ=-。
源程序： k=-10:10; f1k=[zeros(1,10),1,zeros(1,10)]; subplot(2,2,1) stem(k,f1k) title('f1[k]') f2k=[zeros(1,8),1,zeros(1,12)]; subplot(2,2,2) stem(k,f2k) title('f2[k]') f3k=[zeros(1,14),1,zeros(1,6)]; subplot(2,2,3) stem(k,f3k) title('f3[k]') f4k=2*f2k-f3k; subplot(2,2,4) stem(k,f4k) title('f4[k]') ⑶ 在 k [0:31]=范围内产生并画出以下信号： a) ()()k k 144f [k]sin cos π π=； b) ()2k 24f [k]cos π =； c) ()()k k 348f [k]sin cos π π=。 请问这三个信号的基波周期分别是多少？ 源程序： k=0:31; f1k=sin(pi/4*k).*cos(pi/4*k); subplot(3,1,1) stem(k,f1k) title('f1[k]') f2k=(cos(pi/4*k)).^2; subplot(3,1,2) stem(k,f2k) title('f2[k]') f3k=sin(pi/4*k).*cos(pi/8*k); subplot(3,1,3) stem(k,f3k) title('f3[k]') 其中f1[k]的基波周期是4, f2[k]的基波周期是4, f3[k]的基波周期是16。
通信工程学院12级1班罗恒2012101032 实验一语音信号的低通滤波和短时分析综合实验 一、实验要求 1、根据已有语音信号，设计一个低通滤波器，带宽为采样频率的四分之一，求输出信号； 2、辨别原始语音信号与滤波器输出信号有何区别，说明原因； 3、改变滤波器带宽，重复滤波实验，辨别语音信号的变化，说明原因； 4、利用矩形窗和汉明窗对语音信号进行短时傅立叶分析，绘制语谱图并估计基音周期，分析两种窗函数对基音估计的影响； 5、改变窗口长度，重复上一步，说明窗口长度对基音估计的影响。 二、实验目的 1.在理论学习的基础上，进一步地理解和掌握语音信号低通滤波的意义，低通滤波分析的基本方法。 2.进一步理解和掌握语音信号不同的窗函数傅里叶变化对基音估计的影响。 三、实验设备 1.PC机； 2.MATLAB软件环境； 四、实验内容 1.上机前用Matlab语言完成程序编写工作。 2.程序应具有加窗（分帧）、绘制曲线等功能。 3.上机实验时先调试程序，通过后进行信号处理。 4.对录入的语音数据进行处理，并显示运行结果。 5. 改变滤波带宽，辨别与原始信号的区别。 6.依据曲线对该语音段进行所需要的分析，并且作出结论。 7.改变窗的宽度（帧长），重复上面的分析内容。 五、实验原理及方法 利用双线性变换设计IIR滤波器（巴特沃斯数字低通滤波器的设计），首先要设计出满足指标要求的模拟滤波器的传递函数Ha(s),然后由Ha(s)通过双线性变换可得所要设计的IIR滤波器的系统函数H(z)。如果给定的指标为数字滤波器的指标，则首先要转换成模拟滤波器的技术指标，这里主要是边界频率Wp和Ws 的转换，对ap和as指标不作变化。边界频率的转换关系为∩=2/T tan(w/2)。接着，按照模拟低通滤波器的技术指标根据相应设计公式求出滤波器的阶数N和3dB截止频率∩c ；根据阶数N查巴特沃斯归一化低通滤波器参数表，得到归一化传输函数Ha(p)；最后，将p=s/ ∩c 代入Ha(p)去归一，得到实际的模拟滤波器传输函数Ha(s)。之后，通过双线性变换法转换公式s=2/T((1-1/z)/(1+1/z))得到所要设计的IIR滤波器的系统函数H(z)。
通信工程学院12级1班 罗恒 2012101032 实验二 基于MATLAB 的语音信号频域特征分析 一、 实验要求 要求根据已有语音信号，自己设计程序，给出其倒谱、语谱图的分析结果，并根据频域分析方法检测所分析语音信号的基音周期或共振峰。 二、 实验目的 信号的傅立叶表示在信号的分析与处理中起着重要的作用。因为对于线性系统来说，可以很方便地确定其对正弦或复指数和的响应，所以傅立叶分析方法能完善地解决许多信号分析和处理问题。另外，傅立叶表示使信号的某些特性变得更明显，因此，它能更深入地说明信号的各项红物理现象。 由于语音信号是随着时间变化的，通常认为，语音是一个受准周期脉冲或随机噪声源激励的线性系统的输出。输出频谱是声道系统频率响应与激励源频谱的乘积。声道系统的频率响应及激励源都是随时间变化的，因此一般标准的傅立叶表示虽然适用于周期及平稳随机信号的表示，但不能直接用于语音信号。由于语音信号可以认为在短时间内，近似不变，因而可以采用短时分析法。 三、 实验设备 1.PC 机； 2.MATLAB 软件环境； 四、 实验内容 1.上机前用Matlab 语言完成程序编写工作。 2.程序应具有加窗（分帧）、绘制曲线等功能。 3.上机实验时先调试程序，通过后进行信号处理。 4.对录入的语音数据进行处理，并显示运行结果。 5.依次给出其倒谱、语谱图的分析结果。 6. 根据频域分析方法检测所分析语音信号的基音周期或共振峰。 五、 实验原理及方法 1、短时傅立叶变换 由于语音信号是短时平稳的随机信号，某一语音信号帧的短时傅立叶变换的定义为： 其中w(n -m)是实窗口函数序列，n 表示某一语音信号帧。令n -m=k'，则得到 ()()()jw jwm n m X e x m w n m e ∞-=-∞= -∑
语音控制小车实验报告 专业： 学号： 姓名： 2014年01月12日
一、实验目的 语音控制小车以SPCE061A单片机为核心，采用语音识别技术，可通过语音命令对其行驶状态进行控制。本次实验的主要目的： 1．通过简单的I/O 操作实现小车的前进、后退、左转、右转功能； 2．配合SPCE061A 的语音特色，利用系统的语音播放和语音识别资源，实现语音控制的功能； 3．在行走过程中声控改变小车运动状态； 4．在超出语音控制范围时使小车停车。 二、实验内容 1、SPCE061A简介 SPCE061A是一款性价比很高的十六位单片机，使用它可以非常方便灵活的实现语音的录放，该芯片拥有8路10位精度的ADC，其中一路为音频转换通道，并且内置有自动增益电路。这为实现语音录入提供了方便的硬件条件。两路10位精度的DAC，只需要外接功放（SPY0030A）即可完成语音的播放。该单片机具有一套易学易用的指令系统和集成开发环境，在此环境中，它支持标准 C 语言编程，也支持 C 语言与汇编语言的互相调用。另外还提供了语音录放的库函数，只要了解库函数的使用，就可以很容易的完成语音的录放、识别等功能，这些都为软件开发提供了方便的条件。 SPCE061A特性： 16位μ’nSP微处理器； 工作电压：内核工作电压VDD为 3.0V~3.6V(CPU)，I/O口工作电压VDDH为VDD~5.5V(I/O)； CPU时钟：0.32MHz~49.152MHz； 内置2K 字 SRAM； 内置32K 闪存 ROM； 可编程音频处理； 晶体振荡器； 系统处于备用状态下(时钟处于停止状态)，耗电小于 2μA@3.6V； 2 个 16 位可编程定时器/计数器(可自动预置初始计数值)； 2 个 10 位 DAC(数-模转换)输出通道； 32 位通用可编程输入/输出端口；
语音信号处理实验报告 专业班级电子信息1203 学生姓名钟英爽 指导教师覃爱娜 完成日期2015年4月28日 电子信息工程系 信息科学与工程学院
实验一语音波形文件的分析和读取 一、实验学时：2 学时 二、实验的任务、性质与目的： 本实验是选修《语音信号处理》课的电子信息类专业学生的基础实验。通过实验 （1）掌握语音信号的基本特性理论：随机性，时变特性，短时平稳性，相关性等； （2）掌握语音信号的录入方式和*.WAV音波文件的存储结构； （3）使学生初步掌握语音信号处理的一般实验方法。 三、实验原理和步骤： WAV 文件格式简介 WAV 文件是多媒体中使用了声波文件的格式之一，它是以RIFF格式为标准。每个WAV 文件的头四个字节就是“RIFF”。WAV 文件由文件头和数据体两大部分组成，其中文件头又分为RIFF/WAV 文件标识段和声音数据格式说明段两部分。常见的WAV 声音文件有两种，分别对应于单声道（11.025KHz 采样率、8Bit 的采样值）和双声道（44.1KHz 采样率、16Bit 的采样值）。采样率是指声音信号在“模拟→数字”转换过程中，单位时间内采样的次数；采样值是指每一次采样周期内声音模拟信号的积分值。对于单声道声音文件，采样数据为8 位的短整数（short int 00H-FFH）；而对于双声道立体声声音文件，每次采样数据为一个16 位的整数（int），高八位和低八位分别代表左右两个声道。WAV 文件数据块包含以脉冲编码调制（PCM）格式表示的样本。在单声道WAV 文件中，道0 代表左声道，声道1 代表右声道；在多声道WAV 文件中，样本是交替出现的。WAV 文件的格式 表1 wav文件格式说明表
Harbin Institute of Technology 匹配滤波器实验报告 课程名称：信号检测理论 院系：电子与信息工程学院 姓名：高亚豪 学号：14SD05003 授课教师：郑薇 哈尔滨工业大学
1. 实验目的 通过Matlab 编程实现对白噪声条件下的匹配滤波器的仿真，从而加深对匹配滤波器及其实现过程的理解。通过观察输入输出信号波形及频谱图，对匹配处理有一个更加直观的理解，同时验证匹配滤波器具有时间上的适应性。 2. 实验原理 对于一个观测信号()r t ，已知它或是干扰与噪声之和，或是单纯的干扰， 即 0()()()()a u t n t r t n t +?=?? 这里()r t ，()u t ，()n t 都是复包络，其中0a 是信号的复幅度，()u t 是确知的归一化信号的复包络，它们满足如下条件。 2|()|d 1u t t +∞ -∞=? 201||2 a E = 其中E 为信号的能量。()n t 是干扰的均值为0，方差为0N 的白噪声干扰。 使该信号通过一个线性滤波系统，有效地滤除干扰，使输出信号的信噪比在某一时刻0t 达到最大，以便判断信号的有无。该线性系统即为匹配滤波器。 以()h t 代表系统的脉冲响应，则在信号存在的条件下，滤波器的输出为 0000()()()d ()()d ()()d y t r t h a u t h n t h τττττττττ+∞+∞+∞ =-=-+-???
右边的第一项和第二项分别为滤波器输出的信号成分和噪声成分，即 00()()()d x t a u t h τττ+∞ =-? 0 ()()()d t n t h ?τττ+∞ =-? 则输出噪声成分的平均功率（统计平均）为 2 20E[|()|]=E[|()()d |]t n t h ?τττ+∞ -? **00*000200 =E[()(')]()(')d d '=2()(')(')d d ' 2|()|d n t n t h h N h h N h ττττττδττττττττ+∞+∞+∞+∞+∞ ---=?? ?? ? 而信号成分在0t 时刻的峰值功率为 22 20000|()||||()()d |x t a u t h τττ+∞ =-? 输出信号在0t 时刻的总功率为 22000E[|()|]E[|()()|]y t x t t ?=+ 22**0000002200E[|()||()|()()()()] |()|E[|()|] x t t x t t t x t x t t ????=+++=+ 上式中输出噪声成分的期望值为0，即0E[()]0t ?=，因此输出信号的功率 成分中只包含信号功率和噪声功率。 则该滤波器的输出信噪比为 222000022000|||()()d ||()|E[|()|]2|()|d a u t h x t t N h τττρ?ττ+∞ +∞-==?? 根据Schwartz 不等式有
模拟电子技术课程设计 语 音 放 大 器 姓名：伍慧兰 学号：2015550828 班级：15通信工程1班 指导老师：罗光明 目录 一、设计目的 (2) 二、知识点和设计内容 (2) 三、设计方案 (3) 四、实验原理与参考电路 (4) （一）实验原理图如图1-2 (4) （二）实验原理 (5) 1) 前置放大器 (5) 2) 有源带通滤波器 (5) 3) 功率放大器 (6)
五、实验的主要元器件 (7) （一）元器件清单 (7) （二）部分器件的使用介绍 (8) 1) LM324芯片 (8) 2) TDA2030引脚图与应用电路参数 (12) 六、实验步骤 (13) （一）电路仿真实验 (13) （二）硬件实物实验 (19) 1) 前置放大器的焊接与调试 (19) 2) 有源带通滤波器 (20) 七、实验中的问题提出与解决方法 (24) 八、注意事项 (25) 九、实验感想 (26) 参考资料 (27) 语音放大器设计 一、设计目的 1、了解语音识别知识； 2、掌握集成运算放大器的工作原理及其应用； 3、掌握低频小信号放大电路、带通滤波器和功放电路的设计方法； 4、培养应用现代工具对模拟电子系统进行仿真测试、制作调试、故障检查及分析的能力； 5、培养市场素质、工艺素质、自主学习能力、分析问题解决问题的能力以及团队精神； 6、培养文献查阅与综述和撰写课程设计报告的能力。 二、知识点和设计内容 本实验的知识点为分立元件放大器或集成运放、有源滤波器、集
成功率放大器；涉及电子电路各个模块之间的联合调试技术。 三、设计方案 语音放大器设计的基本设计思路 分析可知本语音放大器应包括输入电路、前置放大器、有源带通滤波器、功率放大器、扬声器等几部分组成，如图1-1所示。 前置放大器可采用集成运算放大器，有源带通滤波器可采用LPF 和HPF 串联构成，功率放大电路选用集成功放。 设计的性能指标 通常语音信号非常微弱，需要经过放大、滤波、功率放大后才能驱动扬声器发声。假设语音信号为峰峰值不大于10mV 频率范围100Hz~3kHz 的正弦波，要求驱动8Ω1W 的扬声器。具体性能指标如下： 1、前置放大器：输入信号Uid ≤10mV ；输入阻抗Ri ≥100k Ω 2、有源带通滤波器：通带100Hz~3kHz ；增益Au=1 3、功放：最大不失真输出功率Pomax ≥1W ；负载阻抗R L =8Ω 4、输出功率连续可调；直流输出电压≤50mV ；静态电源电流≤100mA 输入 电路 前置 放大 带通 滤波 功率 放大 图1-1 语音放大电路原理框图
脉搏、语音及图像信号的傅里叶分析 一、实验简介 任何波形的周期信号均可用傅里叶级数来表示。傅里叶级数的各项代表了不同频率的正弦或余弦信号，即任何波形的周期信号都可以看作是这些信号（谐波）的叠加。利用不同的方法，可以从周期信号中分解出它的各次谐波的幅值和相位。也可依据信号的傅里叶级数表达式，将各次谐波按表达式的要求叠加得到所期望的信号。 二、实验目的 1、了解常用周期信号的傅里叶级数表示。 2、了解周期脉搏信号、语音信号及图像信号的傅里叶分析过程 3、理解体会傅里叶分析的理论及现实意义 三、实验仪器 脉搏语音实验仪器，数字信号发生器，示波器 四、实验原理 1、周期信号傅里叶分析的数学基础 任意一个周期为T 的函数f(t)都可以表示为傅里叶级数： 00010000 1 ()(cos sin ) 21()() 1 ()cos()()1 ()sin()()n n n n n f t a a n t b n t a f t d t a f t n t d t b f t n t d t π π π ππ πωωωωπ ωωωπ ωωωπ ∞ =-- - =++=== ∑??? 其中0ω为角频率，称为基频，0a 为常数，n a 和n b 称为第n 次谐波的幅值。任何
周期性非简谐交变信号均可用上述傅里叶级数进行展开，即分解为一系列不同次谐波的叠加。 对于如图1所示的方波，一个周期内的函数表达式为: (0t<)2() (-t 0) 2 h f t h ππ? ≤??=? ?-≤
通信原理多种信号的产生实验实验目的： 1.了解多种时钟信号的产生方法 2.掌握用数字电路产生伪随机序列码的实现方法 3.了解PCM编码中的收、发帧同步信号的产生过程 实验设备及仪器名称： 1.双踪示波器
实验原理： 帧同步信号识别、提取与分析： 由于数字通信系统传输的是一个接一个按节拍传送的数字信号单元,即码元,因而在接收端必须按与发送端相同的节拍进行接收,否则,会因收发节拍不一致而导致接收性能变差。此外,为了表述消息的内容,基带信号都是按消息内容进行编组的,因此,编组的规律在收发之间也必须一致。在数字通信中,称节拍一致为"位同步",称编组一致为"帧同步".在时分复用通信体统中,为了正确地传输信息,必须在信息码流中插入一定数量的帧同步码,它可以是一组特定的码组,也可以是特定宽度的脉冲,可以集中插入,也可以分散插入。若信息数据中含有与帧同步码完全相同的码元序列,则系统将进入错误的同步维持状态,由于在这里是连续传输以24位为周期的周期信号,所以此状态将维持下去。但在实际的信息传输中不会连续传输这种周期信号,因此连续几帧都输出假识别信号的概率极小,所以这种错误的同步维持状态存在的时间是短暂的。 伪随机码： 多种数字信号产生及形成实验：TP101、TP102、TP103 帧同步信号识别、提取与分析实验：TP103、TP104 伪随机码、特殊码观察测量分析实验：TP107、TP108 简易正弦波信号产生及用途实验：TP105、TP106
实验方案： 接好电源，打开开关K2、K3、K100，对应的发光管D6、D7、D1亮，电路加电工作，测量各点波形： TP101、TP301、TP102、TP103、TP104：K102的1、2脚送入信号，分别测试观 察波形并记录。 TP105: K102的1、2脚送入2KHz方波，观察并记录；2、3脚送入1KHz方形，观 察并记录。 TP106: K102的1、2脚送入2KHz简易正弦波形，观察并记录；2、3脚送入1KHz 简易正弦波形，观察并记录。 TP107、TP108：常态下测试观察并记录。 TP108：按下“复位”后，再按下“开始”与“FSK”，观察波形并记录；再按下“复位”后， 再按下“开始”与“PSK”，观察波形并记录。 实验数据：