matlab语音采样系统设计,对采集的语音信号进行matlab设计与分析课程设计报告 30...

对采集的语音信号进行matlab设计与分析课程设计报告 30

语 音 信 号 采 集 与 处 理 主 题： Matlab与信息系统仿真课程设计 组 长： XXX 组 员： XXXXXXX 班 级： 电子信息工程1201 目 录 第一章 课程设计情况概述 1.1 课程设计目的…………………………………………………3 1.2 课程设计内容…………………………………………………3 1.3 小组分工………………………………………………………3 第二章 具体实现方案及小结 2.1 第(1)、(2)版块：录制语以及频谱分析………………………………4 2.1.1 板块要求…………………………………………………………………4 2.1.2 版块分析及具体实现代码………………………………………………4 2.1.3 Matlab仿真结果图……………………………………………………5 2.2 第(3)版块：采样…………………………………………………………7 2.2.1 版块要求…………………………………………………………………7 2.2.2 板块分析及具体实现代码 ……………………………………………7 2.3 第(4)版块：对语音实现几种常见变化……………………………………13 2.3.1 版块要求………………………………………………………………13 2.3.2 板块分析及具体实现代码……………………………………………13 2.4 第(5)版块：调制、解调………………………………………………25 2.4.1 板块要求………………………………………………………………25 2.4.2 版块分析及具体实现代码……………………………………………25 2.4.3 Matlab仿真结果图……………………………………………………28 2.5 第(6)版块：加噪、滤波………………………………………………32 2.5.1 板块要求………………………………………………………………32 2.5.2 版块分析及具体实现代码……………………………………………32 2.6 小结………………………………………………………………………41 第一章 课程设计情况概述 1.1课程设计目的 本课程设计要求掌握Matlab语言特性、数学运算和图形表示；掌握 Matlab程序设计方法在信息处理方面的一些应用；掌握在 Windows 环境下，利用Matlab进行语音信号 采集与处理的基本方法。 综合设计实验方案设计5天分析题目，查找相关资料，利用已掌握知识、原理，提出设计具体方案，拟订设计步骤，确定方案，编写程序代码并上机调试，修改完善设计，记录数据，分析结果。 1.2 课程设计内容 这次的课程设计大致分为两个阶段，一周为一个阶段。 1.关于Matlab的基础训练 Matlab在信号与系统、随机信号处理中的基本应用 2.综合设计：语音信号的采集与处理 (1)利用 windows 自带的录音机或者其它录音软件，录制几段语音信号(要有几种不同的声音，要有男声、女声)。 (2)对录制的语音信号进行频谱分析，确定该段语音的主要频率范围，由此频率范围判断该段语音信号的特点(低沉 or 尖锐？) (3)利用采样定理，对该段语音信号进行采样，观察不同采样频率(过采样、欠采样、临界采样)对信号的影响。 (4)实现语音信号的快放、慢放、倒放、回声、男女变声。 (5)对采集到的语音信号进行调制与解调，观测调制与解调前后信号的变化。 (6)对语音信号加噪，然后进行滤波，分析不同的滤波方式对信号的影响。 1.3 小组分工 第(1)、(2)版块XXX，第(3)版块XXX，第(4)版块XXX，第(5)版块XXX，第(6)版块XXX，报告最终由组长统一整理。 第二章 具体实现方案 2.1 第(1)、(2)版块:录制语音以及频谱分析 2.1.1 版块要求： ①利用windows自带的录音机或者其它录音软件，录制几段语音信号(要有几种不同的声音，要有男声、女声)。 ②对录制的语音信号进行频谱分析，确定该段语音的主要频率范围，由此频率范围判断该段语音信号的特点(低沉还是尖锐？) 2.1.2 版块分析及具体实现代码： 先录制几段音频，再利用Matlab编程对其进行频谱分析。 [x1,fs,nbits]=wavread( D:\matlab\语音样本.wav ); %读取音频文件，识别的是音频存放地址 sound(x1,fs); %利用matlab中发声函数sound语句 N=length(x1); M=log(N)/log(2); Q=round(M); %round函数是对所得数取整 S=2^Q; %求得与N最接近的2的n次幂 n=[0:N-1]; %从0取到N-1 X=fft(x1); %进行傅里叶变换 Fs=1*fs; f=n/N*Fs; %图2的横坐标 figure(1); subplot(2,1,1);plot(n,x1); %绘制原音频信号图 title( 原始音频信号 ); Amp=abs(X)/(S/2); Amp(1,:)=Amp(1,:)/2; subplot(2,1,2);plot(f,Amp); axis([0,3000,0,0.01]) %绘制频谱图 title( 音频的幅度谱 ); 在设计算法的过程中，本段程序的难点在于如何找到与length(x1)最接近的2的n次幂，以及傅立叶变换之后的频谱值(即图中的纵坐标值)。通过先对length(x1)求2的对数，再利用round函数进行四舍五入求整，最终得到想要的最接近length(x1)的2的n次幂。 2.1.3 Matlab仿真结果图： 从图中可以读出这段音频信号的频率大概在100～800赫兹，比较低沉，一般来说应该是男声。 从图中可以读出这段音频信号的频率大概在700～1400赫兹，比较尖锐，一般来说应该是女声。 2.2 第(3)版块：采样 2.2.1 版块要求 利用采样定理，对该段语音信号进行采样，观察不同采样频率(过采样、欠采样、临界采样)对信号的影响。 2.2.2 板块分析及具体实现代码 这里涉及到采样下重采样的问题。所谓重采样就是对一系列离散的数字信号进行插值(升频)或去值(降频