1、随机信号分析与处理实验报告完整版GUI内附完整函数代码随机信号分析与处理实验报告指导教师 :班 级 : 学 号 : 姓 名 :实验一熟悉MATLAB勺随机信号处理相关命令一、 实验目的1、 熟悉GUI格式的编程及使用。2、 掌握随机信号的简单分析方法3、 熟悉语音信号的播放、波形显示、均值等的分析方法及其编程二、 实验原理1、 语音的录入与打开在MATLAB, y,fs,bits=wavread(Blip,N1 N2); 用于读取语音,采样值放在向量 y中,fs表示采样频率(Hz),bits表示采样位数。N1 N2表示读取从N1点到N2点的值。2、 均匀分布白噪声在matlab中,有x=ran
2、d (a,b)产生均匀白噪声序列的函数,通过与语言信号的叠加来分析其特性。3、 均值随机变量X的均值也称为数学期望,它定义为e+oc对于离散型随机变量,假定随机变量 X有N个可能取值,各个取值的概率为- p y -1则均值定义为(X) = .rfPf/=1上式表明,离散型随机变量的均值等于随机变量的取值乘以取值的概率之和,如果取值是等概率的,那么 均值就是取值的算术平均值,如果取值不是等概率的,那么均值就是概率加权和,所以,均值也称为统计 平均值。4、方差定义为随机过程r 的方差。方差通常也记为 d【x(t)】,随机过程的方差也是时间 t的函数,由方差的 定义可以看岀,方差是非负函数。5、自相
3、关函数设任意两个时刻t1,t2,定义 :Rx (叩2)= EX(tJX(t2) = Jq JX1X2 f (X1, X2,t1,t2)dX1dX2为随机过程X(t)的自相关函数,简称为相关函数。自相关函数可正,可负,其绝对值越大表示相关性越 强。6. 哈明(hamming)窗0.54+0.46(10.100)- (10.101)B = 1.3 f , A = -43dB , D= -6dB/oct.哈明窗本质上和汉宁窗是一样的,只是系数不同。哈明窗比汉宁窗消除旁瓣的效果好一些而且主瓣稍窄,但是旁瓣衰减较慢是不利的方面。适当地改变系数,可得到不同特性的窗函数。三、实验结果分析1,语言信号 (a)
4、 y的原始信号图为原始语音信号的原始图及取不同样值时的效果,取样值越多变化越快。2.加入均匀分布的白噪声后的信号3.求均值b)均值1E0.40.20-32的平均64的平均10 20 30 40 50 60 70 90 90 100上图所示为对语音信号求均值的图形,在求均值时需要对信号加窗处理,否则岀现的就是一个点不便于分析,所以上图是对信号加汉宁窗分别对 32点和64采样点求均值的图形,由原理知道不同的采样点的均值是不一样的,由图可以看岀64点的均值比32点的均值的图形更平滑,更接近与平行横轴,当继续加大采样 点数时,最后均值会平行于横轴,所以当采样点数达到足够大时均值是一个常数。4.求方差)
5、.432的方差64的方差C)方差由图可知:自相关函数可正,可负,其绝对值越大表示相关性越强。00 10 20 30 40 50 60 70 80 90 10如图所示是对两个不同的采样点数求方差的图形,同均值一样,求方差是也需要加床便于分析,由图可看书方差曲线总是在横轴的上方,即方差也总是大于 0,图中32点和64点的方差对比可以看出不同的采样点数对应的方差不相同,5 自相关函数四实验心得通过这次实验使我对随机信号处理这门课有了新的认识与理解,对 matlab这款软件有了更新的认识与了解,初次认识了什么事 GUI,什么是工程与设计及懂得了查找资料的主要性。同时,通过实验也发现了自 己很多不足与知
6、识的短缺,特别是有关专业方面的知识点掌握的不是很好以后应该抓紧补回来,在今后的 学习中也知道了什么是学习的重要方向。实验二随机信号处理的工程编程实现一、 实验目的1、 熟悉各种随机信号分析及处理方法。2、 掌握运用MATLAB的统计工具包和信号处理工具包绘制概率密度的方法二、 实验原理1. 平滑滤波平滑滤波可以与中值滤波结合使用, 对应的线性平滑器可以仅仅用低阶的低通滤波器 (如果采用高阶的系统,则将抹掉信号中应该保存的不连续性) 。2. IIR数字滤波器设计原理利用双线性变换设计IIR滤波器(只介绍巴特沃斯数字滤波器的设计) ,首先要设计岀满足指标要求的模拟滤波器的传递函数 H a(s),然
7、后由HS)通过双线性变换可得所要设计的 |R滤波器的系统函数H (z)。3. 协方差设两个随机变量 X和 丫,定义: covX,Y)二 EX_E(X)Y_E(Y) = E(XY)_E(X)E(Y)为X和丫的协方差。cov( X ,Y )其相关函数为:rxy : : 由此可见协方差的相关性与x和丫是Jd(X )D(Y)密切相关的,表征两个函数变化的相似性4. 互相关互相关函数定义为:Rxy魚,t2)= EX(tJY(t2) = ” xyfXY(x, y,t“t2)dxdy如果X(t)与Y(t)是相互独立的,则一定是不相关的。反之则不一定成立。它是两个随机过程联合统计 特性中重要的数字特征。5.最
8、大似然法 功率谱估值的表达式 给定信号x(n),其最似然法功率谱估值为式中;t为采样时间间隔;Rx为信号x(n)的自相关矩阵;R为Rx的逆矩阵;T为转置运算;*为取共轭 值。满足上述要求的滤波器系数 a的表达式为式中由上式可以看岀,滤波器系数与信号的自相关函数和 E有关。可以看为,滤波器将根据输入的信号及所要求的频率而调整其系数,使所关心的频率分量能完全通过, 而使其他频率分量的输岀功率最小。 因此,它能得到比使用固定的窗口函数的周期图法更高的分辨率。6时域信号的FFT分析FFT即为快速傅里叶变换,是离散傅里叶变换的快速算法,它是根据离散傅里叶变换的奇、偶、虚、实等特性,对离散傅里叶变换的算法
9、进行改进获得的。在 MATLAB的信号处理工具箱中函数 FFT的一种调用格式为y - fft(x)其中X是序列,Y是序列的FFT。7白噪声的检测与分析白噪声信号是一个均值为零的随机过程,任一时刻是均值为零的随机变量。而服从高斯分布的白噪声即称为高斯白噪声。8希尔伯特变换及性质x(t)的希尔伯特变换为x(t)与1/ n t的卷积,即因此,对x(t)的希尔伯特变换可以看作为 x(t)通过一个冲击响应为1/ n t的线性滤波器。希尔伯特变换 器在整个频域上具有恒为1的幅频特性,为全通网络,在相位上则引入 -n /2和n /2的相移9功率谱密度估计定义 1 I 2Gx)=*2齐冷创jTXt(.存宀水为
10、随机信号的功率谱,它表示单位频带类信号的频率分量消耗在单位电阻上的平均功率的统计平均值, 10瑞利分布概率密度函数raylpdf()使用方法:Y=raylpdf(x,b)参数为b的瑞利分布在x处的概率密度函数。11 .瑞利分布概率分布函数raylcdf()使用方法:Y=raylcdf(x,b)参数为b的瑞利分布在x处的概率分布函数值12. 巴特沃斯滤波器巴特沃斯滤波器的特点是通频带内的频率响应曲线最大限度平坦, 没有起伏,而在阻频带则逐渐下降为零。在振幅的对数对角频率的波特图上_ ,从某一边界角频率开始,振幅随着角频率的增加而逐步减少 ,趋向负无穷大。一阶巴特沃斯滤波器的衰减率为每倍频 6分贝
11、,每十倍频20分贝。二阶巴特沃斯滤波器的衰减率为每倍频12分贝、三阶巴特沃斯滤波器的衰减率为每倍频 18分贝、如此类推。巴特沃斯滤波器的振幅对角频率单调下降,并且也是唯一的无论阶数,振幅对角频率曲线都保持同样的形状的滤波器。只不过滤波器 阶数越高,在阻频带振幅衰减速度越快。其他滤波器高阶的振幅对角频率图和低级数的振幅对角频率有不 同的形状。13. 直方图法估计随机序列概率密度的估计使用方法:hist(y,x), 画出用矢量y表示的随机序列的直方图,参数 x表示计算直方图划分的单元,也是用矢量表示。例:x=-2.9:0.1:2.9;y=normrnd(0,1,1000,1);hist(y,x);
12、14. 概率密度如果FY(x,t)的一阶导数存在,则定义dxa x生成n*n的包含标准正态分布的随机矩阵。生成的m*n随机矩阵16.韦伯分布随机变量分布之一。又称韦伯分布、韦氏分布或威布尔分布,由瑞典物理学家 Waloddi Weibull 于1939年引进,是可靠性分析及寿命检验的理论基础。威布尔分布(山型极值分布)记为W(k,a,b)。威布尔分布:在可靠性工程中被广泛应用,尤其适用于机电类产品的磨损累计失效的分布形式。 由于它可以利用概率值很容易地推断岀它的分布参数,被广泛应用于各种寿命试验的数据处理。四、实验结果及分析1.平滑滤波由图知,经过平滑滤波后,原始信号的峰值变化减小了,信号的频
13、谱变得平滑了很多。说明平滑滤波对信 号具有很好的平滑效果。2. IIR低通滤波(a) y的原始信号10 20 30 40 50 60 70 80 90 100取值范围II(b)经过IIR低通谑波后111 H 1 1 1J-110 911L 1 1 J11()102030 40 50 60 708090 100取值范围由图分析,可知经过低通滤波后信号的波形变化很大, 说明原始信号主要分布在高频不分, 低频占据很少。3.高通滤波20, r r , L| IIR高通泯液后信亏频谓40 1 1 1 J 1 1 1 1 J 0 10 20 30 40 50 60 70 80 90 100取值范围由图可知
14、:经过滤波后基本不改变原始信号的波形幅值及频谱, 说明该语言信号高频部分占据很多。4.最大似然估计由图可知:语音信号的最大似然估计满足由 0到0.5变化,呈直线变化,随着选取点数的增加越缓慢。5.互相关(a)正态分布JS始信号50-51l I I 1 1 11110 10 20 30 40 50 60 70 30 90 100(b) y的原始信号10 -1 LI T 1 1 1 K 010 20 30 40 60 60 70 60 90 1010010韦伯伯号/卜) /vo 10 2030 4050 60 70 80 901001 1 1 12 1,3 14 1.5 16 1.7 1.6 1,9 2由图可知:协方差反映了两个之间的相关程度, 又图可知两个的差别变化很大。 在大约1.7的时