隨机信号处理实验报告一.docVIP

《随机信号分析与处理》 实验报告 指导教师 ：廖 红 华 班 级 ：0309411 学 号 ：030941103 姓 名 ：钱进红 2011-12-7 实验一 熟悉MATLAB的随机信号处理相关命令 一、实验目的 1、熟悉GUI格式的编程及使用。 2、掌握随机信号的简单分析方法 3、熟悉语音信号的播放、波形显示、均值等的分析方法及其编程 二、实验原理 1、语音的录入与打开 在MATLAB中，[y,fs,bits]=wavread(Blip,[N1 N2]);用于读取语音，采样值放在向量y中，fs表示采样频率(Hz)，bits表示采样位数。[N1 N2]表示读取从N1点到N2点的值。 2、时域信号的FFT分析 FFT即为快速傅变换，是离散傅变换的快速算法，它是根据离散傅变换的奇、偶、虚、实等特性，对离散傅叶变换的算法进行改进获得的。 其中X是序列，Y是序列的FFT。 3、均值 随机变量X 的均值也称为数学期望，它定义为 对于离散型随机变量，假定随机变量X 有N 个可能取值，各个取值的概率为则均值定义为 上式表明，离散型随机变量的均值等于随机变量的取值乘以取值的概率之和，如果取值是等概率的，那么均值就是取值的算术平均值，如果取值不是等概率的，那么均值就是概率加权和，所以，均值也称为统计平均值。 4、方差 定义 为随机过程的方差。方差通常也记为D【X（t）x(t) 的希尔伯特变换为x(t) 与1/πt 的卷积，即 因此，对x(t) 的希尔伯特变换可以看作为x(t) 通过一个冲击响应为1/πt 的线性滤波器。希尔伯特变换器在整个频域上具有恒为1 的幅频特性，为全通网络，在相位上则引入?π/2 和π/2的相移 6、自相关函数 设任意两个时刻，，定义 为随机过程X（t）的一阶导数存在，则定义 为随机过程 的一维概率密度。如果知道了随机过程的一维概率密度，那么也就知道了随机过程在所有时刻上随机变量的一维概率密度。 8、功率谱密度估计 定义 为随机信号的功率谱，它表示单位频带类信号的频率分量消耗在单位电阻上的平均功率的统计平均值， 9，椭圆滤波 椭圆滤波器（Elliptic filter）又称考尔滤波器（Cauer filter），是在通带和阻带等波纹的一种滤波器。椭圆滤波器相比其他类型的滤波器，在阶数相同的条件下有着最小的通带和阻带波动。它在通带和阻带的波动相同，这一点区别于在通带和阻带都平坦的巴特沃斯滤波器，以及通带平坦、阻带等波纹或是阻带平坦、通带等波纹的切比雪夫滤波器。 10.巴特沃斯滤波器 巴特沃斯滤波器的特点是通频带内的频率响应曲线最大限度平坦，没有起伏，而在阻频带则逐渐下降为零。 在振幅的对数对角频率的波特图上,从某一边界角频率开始,振幅随着角频率的增加而逐步减少,趋向负无穷大一阶巴特沃斯滤波器的衰减率为每倍频6分贝，每十倍频20分贝。二阶巴特沃斯滤波器的衰减率为每倍频12分贝、 三阶巴特沃斯滤波器的衰减率为每倍频18分贝、如此类推。巴特沃斯滤波器的振幅对角频率单调下降，并且也是唯一的无论阶数，振幅对角频率曲线都保持同样的形状的滤波器。只不过滤波器阶数越高，在阻频带振幅衰减速度越快。其他滤波器高阶的振幅对角频率图和低级数的振幅对角频率有不同的形状。 2.FFT变换 3.均值 上图所示为对语音信号求均值的图形，在求均值时需要对信号加窗处理，否则出现的就是一个点不便于分析，所以上图是对信号加汉宁窗分别对10点和20采样点求均值的图形，由原理知道不同的采样点的均值是不一样的，由图可以看出20点的均值比10点的均值的图形更平滑,更接近与平行横轴，当继续加大采样点数时，最后均值会平行于横轴，所以当采样点数达到足够大时均值是一个常数。 4.方差 如图所示是对两个不同的采样点数求方差的图形，同均值一样，求方差是也需要加床便于分析，由图可看书方差曲线总是在横轴的上方，即方差也总是大于0，图中64点和128点的方差对比可以看出不同的采样点数对应的方差不相同，64点的方差最大值为： ，128点的方差的最大值为：，所以当点数取的越大时对应得方差曲线越平滑，当点数达到足够大时也是一个常数。 5.功率谱密度： 上图所示为原始信号与加噪信号的功率谱密度的比较图形，上述两个图形中功率谱都大于或等于0，即功率谱为非负的函数，原始信号中在0.8处功率谱最大其他点的功率谱基本上都可看成是0即常数，加噪信号中，功率谱都在0到0.