1、信号与系统课题Matlab语音信号处理姓名:陈冠英学号：12350004专业：自动化一、Matlab简介 MATLAB是功能强大的科学及计算软件，它不但具有以矩阵计算为根底的强大数学计算和分析功能，而且还具有丰富的可视化图形表现功能和方便的程序设计能力。MATLAB的应用领域极为广泛，除了数学计算和分析外，还被广泛地应用于自动控制、系统仿真、数字信号处理、图形图像分析、数理统计、人工智能、虚拟现实技术、通信工程、金融系统等领域，因此，MATLAB是面向21世纪的计算机程序设计及科学计算语言。 随着MATLAB/SIMULINK通信、信号处理专业函数库和专业工具箱的成熟，它们在通信理论研究、算法
2、设计、系统设计、建模仿真和性能分析验证等方面的应用也更加广泛。MATLAB软件集数值分析、信号处图形显示于一体，且界面友好、具有强大的专业函数库和工具箱，在数字信号处理的科学研究中越来越重要的计算和仿真验证工具。二、语音信号的录制使用的声卡设备采集一段语音信号，并将其保存在中。 图1 PC机的语音信号采集过程由于本机的系统为WIN8，其自带的录音机所保存的录音格式不是WAV，而MATLAB只识别WAV的语音信号，于是我通过一个录音软件Moo0 Voice Recorder 1.40录制一段语音“我爱信号与系统由于接下来的分析处理。图2 录音界面三、 信号的时域分析与频域分析3.1 语音信号的采
3、样与回放我们使用Matlab里面的wavread(filename)函数对其进行读取。其函数具体如下 y,fs,nbits= wavread(filename); 其中y就是音频信号，fs是采样频率，nbit是采样精度。由于录制得到的语音信号为双声道，我们只取其中一声道进行分析。程序片段如下：运行后即可听到“我爱信号与系统的语音。3.2时域分析 我将其转化为了时域模式。 程序片段： 图3 语音时域图这样我们就得到了语音的时域图像，可以看到信号的波动与“我爱信号与系统的口音相同。3.3频域分析 我们已经学了连续信号的傅里叶变换，在这里我们需要处理数字信号。我们使用FFT函数来分析信号的频域特性。
4、 FFT是离散傅立叶变换的快速算法，可以将一个信号变换到频域。有些信号在时域上是很难看出什么特征的，但是如果变换到频域之后，就很容易看出特征了。这就是很多信号分析采用FFT变换的原因。另外，FFT可以将一个信号的频谱提取出来，这在频谱分析方面也是经常用的。在使用FFT(x,N)时需要注意:假设采样频率为Fs，采样点数为N，做FFT之后，某一点nn从1开始表示的频率为：Fn=(n-1)xFs/N；该点的模值除以N/2就是对应该频率下的信号的幅度对于直流信号是除以N；该点的相位即是对应该频率下的信号的相位。 由于FFT结果的对称性，通常我们只使用前半局部的结果，即小于采样频率一半的结果。 下面是程
5、序片段：图4 频谱图上 FFT模值下将FFT模值得很坐标缩小一下，我们得到以下列图： 图5 FFT模值缩小横坐标可以看到信号的频率主要分布在100Hz-1KHZ的频段上，即语音信号为低频信号。四、 信号加噪与分析MATLAB中产生高斯白噪声非常方便，可以直接应用两个函数，一个是WGN，另一个是AWGN。WGN用于产生高斯白噪声，AWGN那么用于在某一信号中参加高斯白噪声。在这里我们直接使用randn(),在MATLAB中可以用randn产生均值为0方差为1的正态分布白噪声。程序片段： 所得图像如下：图6 噪音时域图上信号加噪音时域图下图7 噪音频域图上信号加噪音频域图下运行程序后，首先会听到沙
6、沙的噪声，接着听到混有沙沙声的“我爱信号与系统。我们可以看到噪声的频率范围遍布低频与高频，使得信号在低频与高频均有值，我们再下面需要滤去高频成分。五、进行幅度调制与分析为了减少在传输时的耗损，人们一般是先对传输信号进行特殊处理，然后再传递。把原始的待传信号托附到高频振荡的过程称为调制。我们知道音频信号的频率在我们的听觉范围内，音频信号在无线传输的过程中会受到各种声音的干扰而产生能量消耗！因此限制了传输的距离！所以在现实的传送过程中要现对语音信号进行调制再发送出去！ 设正弦载波为 式中，A为载波幅度；cw为载波角频率；0j为载波初始相位(通常假设0j=0).调制信号基带信号为m(t)。根据调制的定义，振幅调制信号已调信号一般可以表示为设调制信号m(t)的频谱为M(w)，那么已调信号Sm(t)的频谱Sm(w)：在这里我们使用载波函数modulate=cos(2023xxt),载波频率为10KHz, 录音的频率主要分布在100Hz-1Khz。程序片段如下：运行程序后，由于在程序中回放了录音，会听到声音，但是声音与之前差异很大，伴有吱吱响。我取了调制前后时域与频域波形进行比拟，如以下列图：图8 调制前时域与频域左上右上 调制后时域与频域左下右下
