《随机信号分析与处理》实验报告完整版GUI内附完整函数代码.docx
- 文档编号:28410840
- 上传时间:2023-07-13
- 格式:DOCX
- 页数:13
- 大小:303.41KB
《随机信号分析与处理》实验报告完整版GUI内附完整函数代码.docx
《《随机信号分析与处理》实验报告完整版GUI内附完整函数代码.docx》由会员分享,可在线阅读,更多相关《《随机信号分析与处理》实验报告完整版GUI内附完整函数代码.docx(13页珍藏版)》请在冰豆网上搜索。
《随机信号分析与处理》实验报告完整版GUI内附完整函数代码
随机信号分析与处理》
实验报告
指导教师:
班级:
学号:
姓名:
实验一熟悉MATLAB勺随机信号处理相关命令
一、实验目的
1、熟悉GUI格式的编程及使用。
2、掌握随机信号的简单分析方法
3、熟悉语音信号的播放、波形显示、均值等的分析方法及其编程
二、实验原理
1、语音的录入与打开
在MATLAB^,[y,fs,bits]=wavread('Blip',[N1N2]);用于读取语音,采样值放在向量y中,fs表示采
样频率(Hz),bits表示采样位数。
[N1N2]表示读取从N1点到N2点的值。
2、均匀分布白噪声
在matlab中,有x=rand(a,b)产生均匀白噪声序列的函数,通过与语言信号的叠加来分析其特性。
3、均值
随机变量X的均值也称为数学期望,它定义为
e+oc
对于离散型随机变量,假定随机变量X有N个可能取值,各个取值的概率为-py--1则均值定义为
£(X)=£.rfPf
/=1
上式表明,离散型随机变量的均值等于随机变量的取值乘以取值的概率之和,如果取值是等概率的,那么均值就是取值的算术平均值,如果取值不是等概率的,那么均值就是概率加权和,所以,均值也称为统计平均值。
4、方差
定义
为随机过程<r>的方差。
方差通常也记为d【x(t)】,随机过程的方差也是时间t的函数,由方差的定义可以看岀,方差是非负函数。
5、自相关函数
设任意两个时刻t1,t2,定义:
:
:
:
Rx(叩2)=E[X(tJX(t2)]=JqJX1X2f(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)y的原始信号
图为原始语音信号的原始图及取不同样值时的效果,取样值越多变化越快。
2.加入均匀分布的白噪声后的信号
3.求均值
1E
—
0.4
0.2
0
-32的平均
64的平均
102030405060709090100
上图所示为对语音信号求均值的图形,在求均值时需要对信号加窗处理,否则岀现的就是一个点不便于分
析,所以上图是对信号加汉宁窗分别对32点和64采样点求均值的图形,由原理知道不同的采样点的均值
是不一样的,由图可以看岀64点的均值比32点的均值的图形更平滑,更接近与平行横轴,当继续加大采样点数时,最后均值会平行于横轴,所以当采样点数达到足够大时均值是一个常数。
4.求方差
).4
32的方差
64的方差
由图可知: 自相关函数可正,可负,其绝对值越大表示相关性越强。 0 010203040506070809010 如图所示是对两个不同的采样点数求方差的图形,同均值一样,求方差是也需要加床便于分析,由图可看 书方差曲线总是在横轴的上方,即方差也总是大于0,图中32点和64点的方差对比可以看出不同的采样 点数对应的方差不相同, 5•自相关函数 四•实验心得 通过这次实验使我对随机信号处理这门课有了新的认识与理解,对matlab这款软件有了更新的认识与 了解,初次认识了什么事GUI,什么是工程与设计及懂得了查找资料的主要性。 同时,通过实验也发现了自己很多不足与知识的短缺,特别是有关专业方面的知识点掌握的不是很好以后应该抓紧补回来,在今后的学习中也知道了什么是学习的重要方向。 实验二随机信号处理的工程编程实现 一、实验目的 1、熟悉各种随机信号分析及处理方法。 2、掌握运用MATLAB^的统计工具包和信号处理工具包绘制概率密度的方法 二、实验原理 1.平滑滤波 平滑滤波可以与中值滤波结合使用,对应的线性平滑器可以仅仅用低阶的低通滤波器(如果采用高阶的 系统,则将抹掉信号中应该保存的不连续性)。 2.IIR数字滤波器设计原理 利用双线性变换设计IIR滤波器(只介绍巴特沃斯数字滤波器的设计),首先要设计岀满足指标要求的 模拟滤波器的传递函数Ha(s),然后由H'S)通过双线性变换可得所要设计的||R滤波器的系统函数 H(z)。 3.协方差 设两个随机变量X和丫,定义: covX,Y)二E{[X_E(X)][Y_E(Y)]}=E(XY)_E(X)E(Y) 为X和丫的协方差。 cov(X,Y) 其相关函数为: rxy: —: 由此可见协方差的相关性与x和丫是 Jd(X)D(Y) 密切相关的,表征两个函数变化的相似性 4.互相关 互相关函数定义为: Rxy魚,t2)=E{X(tJY(t2)}=”xyfXY(x,y,t“t2)dxdy 如果X(t)与Y(t)是相互独立的,则一定是不相关的。 反之则不一定成立。 它是两个随机过程联合统计特性中重要的数字特征。 5.最大似然法功率谱估值的表达式给定信号x(n),其最似然法功率谱估值为 式中;t为采样时间间隔;Rx为信号x(n)的自相关矩阵;R为Rx的逆矩阵;T为转置运算;*为取共轭值。 满足上述要求的滤波器系数a的表达式为 式中 由上式可以看岀,滤波器系数与信号的自相关函数和E有关。 可以看为,滤波器将根据输入的信号及 所要求的频率而调整其系数,使所关心的频率分量能完全通过,而使其他频率分量的输岀功率最小。 因此, 它能得到比使用固定的窗口函数的周期图法更高的分辨率。 6•时域信号的FFT分析 FFT即为快速傅里叶变换,是离散傅里叶变换的快速算法,它是根据离散傅里叶变换的奇、偶、 虚、实等特性,对离散傅里叶变换的算法进行改进获得的。 在MATLAB的信号处理工具箱中函数FFT 的一种调用格式为 y-fft(x) 其中X是序列,Y是序列的FFT。 7•白噪声的检测与分析 白噪声信号是一个均值为零的随机过程,任一时刻是均值为零的随机变量。 而服从高斯分布的白噪声即 称为高斯白噪声。 8希尔伯特变换及性质 x(t)的希尔伯特变换为x(t)与1/nt的卷积,即 因此,对x(t)的希尔伯特变换可以看作为x(t)通过一个冲击响应为1/nt的线性滤波器。 希尔伯特变换器在整个频域上具有恒为1的幅频特性,为全通网络,在相位上则引入-n/2和n/2的相移 9•功率谱密度估计 定义1I2 Gx④)=£*2齐冷创j T Xt(".存⑴宀水 为随机信号的功率谱,它表示单位频带类信号的频率分量消耗在单位电阻上的平均功率的统计平均值,10•瑞利分布概率密度函数 raylpdf() 使用方法: Y=raylpdf(x,b) 参数为b的瑞利分布在x处的概率密度函数。 11.瑞利分布概率分布函数 raylcdf() 使用方法: Y=raylcdf(x,b) 参数为b的瑞利分布在x处的概率分布函数值 12.巴特沃斯滤波器 巴特沃斯滤波器的特点是通频带内的频率响应曲线最大限度平坦,没有起伏,而在阻频带则逐渐下降为零。 在振幅的对数对角频率的波特图上_,从某一边界角频率开始,振幅随着角频率的增加而逐步减少,趋向负无 穷大。 一阶巴特沃斯滤波器的衰减率为每倍频6分贝,每十倍频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); 14.概率密度 如果FY(x,t)的一阶导数存在,则定义 dx ax 生成n*n的包含标准正态分布的随机矩阵。 生成的m*n随机矩阵 16.韦伯分布 随机变量分布之一。 又称韦伯分布、韦氏分布或威布尔分布,由瑞典物理学家WaloddiWeibull于 1939年引进,是可靠性分析及寿命检验的理论基础。 威布尔分布(山型极值分布)记为W(k,a,b)。 威布尔分布: 在可靠性工程中被广泛应用,尤其适用于机电类产品的磨损累计失效的分布形式。 由于它可以利用概率值很容易地推断岀它的分布参数,被广泛应用于各种寿命试验的数据处理。 四、实验结果及分析 1.平滑滤波 由图知,经过平滑滤波后,原始信号的峰值变化减小了,信号的频谱变得平滑了很多。 说明平滑滤波对信号具有很好的平滑效果。 2.IIR低通滤波 (a)y的原始信号 102030405060708090100 取值范围 II (b)经过IIR低通谑波后 1 1 1H111 J- 1 1 0 9 1 1 L11J 1 1 ( ) 10 20 3040506070 80 90100 取值范围 由图分析,可知经过低通滤波后信号的波形变化很大,说明原始信号主要分布在高频不分,低频 占据很少。 3.高通滤波 20,rr,L~|IIR高通泯液后信亏频谓 40111J1111J 0102030405060708090100 取值范围 由图可知: 经过滤波后基本不改变原始信号的波形幅值及频谱,说明该语言信号高频部分占据很 多。 4.最大似然估计 由图可知: 语音信号的最大似然估计满足由0到0.5变化,呈直线变化,随着选取点数的增加越 缓慢。 5.互相关 (a)正态分布JS始信号 5 0 -5 ■1 ™l~"I~"I~"1~"1~ 1111 0102030405060703090100 (b)y的原始信号 1 0- -1L IT「111K\\ 0 10203040606070609010[ 10 0 10 ~T~1IIn1H1~ 11111J1J 0204060SD10&120140160160200 由图可知: 相关性可正课负,在以0为基准上下改变。 大约在100的时候达到最大值,此后向两 边多趋于平稳,接近于0. 6.协方差 4 2 n Ca>韦伯伯号 /卜) \/ vo1020 3040 5060708090 100 111121,3141.5161.71.61,92 由图可知: 协方差反映了两个之间的相关程度,又图可知两个的差别变化很大。 在大约1.7的时
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 随机信号分析与处理 随机 信号 分析 处理 实验 报告 完整版 GUI 完整 函数 代码
![提示](https://static.bdocx.com/images/bang_tan.gif)