版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
1、高级数字信号处理大作业 学 院 专 业 学 号 学生姓名 2016年 月 日1.仿真题目For this computer experiment involving the LMS algorithm, use a first-order, autoregressive (AR) process to study the effects of ensemble averaging on the transient characteristics of the LMS algorithm for real data.Consider an AR process of order one, des
2、cribed by the difference equationwhere is the (one and only) parameter of the process and is a zero-mean white-noise process of variance . To estimate the parameter , use an adaptive predictor of order one, as depicted in Fig. P4.Given different sets of AR parameters, fixed the step-size parameter,
3、and the initial condition . ·Please plot the transient behavior of weight including the single realization and ensemble-averaged result. ·The transient behavior of the squared prediction error. ·Draw the experimental learning curves, what is the result when the step-size parameter is
4、reduced.2.题目分析已知满足一阶AR过程,其差分方程为 (2-1)这里a是这个过程的参数,是零均值,方差为的白噪声。为了估计参数a,我们使用一阶自适应预测器,初始值,自适应算法采用LMS算法。预测器抽头权值的LMS自适应算法形式为 (2-2)其中 (2-3)是预测误差。题目要求我们用计算机对一阶AR过程随机信号的LMS滤波进行仿真,通过绘制不同过程参数情况下,单一实现和集平均时一阶自适应预测器权值和平方预测误差的瞬态特性图,来研究实时数据集平均对LMS算法瞬态特性的影响。最后,绘出不同步长因子情况的学习曲线图后,分析学习曲线在步长因子减小情况下的变化情况。仿真的目的是通过计算机仿真来研
5、究输入到LMS滤波器的信号参数、步长因子对LMS算法性能的影响。因此在仿真实验中需要改变题目中的AR参数a和步长因子,并对实验得到的权系数和均方误差的瞬态特性进行分析,来获得这些因素对LMS算法性能的影响情况。3.仿真原理(1). LMS算法LMS(Least-Mean-Square)又称最小均方算法是自适应滤波器中常用的一种算法,与维纳算法不同的是,其系统的系数是随输入序列而改变的。维纳算法中截取输入序列自相关函数的一段构造系统的最佳系数。而LMS算法则是对初始化的滤波器系数依据最小均方误差准则进行不断修正来实现的。因此,理论上讲LMS算法的性能在同等条件下要优于维纳算法,但是LMS算法是在
6、一个初始化值的基础上进行逐步跟踪调整得到的,这就使得系统在进入稳定之前有了一个调整的时间,这个时间受到算法步长因子的控制,在一定数值范围内,增大会减小调整时间,但超过这个值范围时系统收效性能会降低或者不再收敛,的最大取值为的迹。LMS算法表达式:设计参数:时刻的输入数据矢量时刻的期望响应时刻的滤波器系数矢量系数的数目步长参数初始化: (3-1)计算:对于有 (3-2) (3-3) (3-4)(2). 实验数据的产生的方差未知,步长因子未知,这两个需要我们自己给定。输入数据和噪声也没有给出,需要我们自己产生。 的产生:由知,给定,对于的全部唯一确定。 的产生:由Yule-Walker方程 (3-
7、5)可以得到 (3-6)4.仿真过程(1).实验涉及到的各个参数条件设定如下:AR参数;步长因子0.01、0.04、0.08;高斯白噪声的功率为;AR序列的长度为1000;实验次数为100。(2).仿真步骤:实验开始先定义以上参变量。对于不同的参数、步长因子和试验次数M分别进行N次迭代来获得每个点的权值系数。对于不同的参数和步长因子分别计算它们的M次集平均权值特性和均方误差特性。绘出三幅图,分别是该一阶自适应预测器的权值和平方误差瞬时特性图以及不同步长因子下的学习曲线图。5.仿真结果通过运行附录程序可以得出该一阶自适应预测器权值的瞬态特性图(=0.04)如下所示:图1 一阶自适应预测器权值的瞬
8、态特性(=0.04)图一中,显示了平均的学习曲线及一个实现,其中实线所示为某一单独实验中的得到的权值瞬态特性曲线,虚线表示的是在100次实验后得到的一个平均结果。观察两条曲线发现,虚线较实线平滑。这是因为集平均实现采用的是100次试验的平均处理,平滑了单一处理中梯度噪声的影响。图2一阶自适应预测器的平方预测误差瞬时特性(=0.04)图二中,是相应的一阶自适应预测器的平方预测误差与迭代次数的关系图,可见,LMS算法单一实现的学习曲线呈现出严重的噪声的形式,但经总体平均后得到了一条较稳定的曲线,即=0.04时的一阶自适应预测器的学习曲线。图3 一阶自适应预测器的学习曲线(变步长) 在图三中,我们可
9、以看到当步长参数分别取0.01、0.04和0.08时均方误差的学习曲线。图中显示,影响收敛速度,也影响稳态值。当时算法大约在100次迭代收敛,而当时大约需要250次迭代才收敛,而前者的稳态值要比后者的高。当时需要700次的迭代。综上所述随着步长因子的减小,LMS算法就需要进行更多的迭代才能收敛,即到达最优点的时间越长。6.结论由以上实验结果分析可归纳为:(1)单一实现的轨迹和学习曲线明显是随机的或是“有噪声的”,而随着实验次数的增多,其平均结果趋于平滑。即总体平均有一个平滑的作用。(2)LMS算法的收敛速度依赖于步长。在允许的范围内值取得越大,收敛速度越快,反之步长越小,收敛速度就越慢。LMS
10、算法是一种相对简单而性能又十分优越的自适应算法。在对LMS算法进行应用或设计的时候,我们可以用变步长方法来缩短其自适应收敛过程,为了达到快速收敛的目的,必须选择合适的步长因子的值,一个可能的策略是尽可能多的减少瞬时平方误差,即用瞬时平方误差作为均方误差的简单估计。LMS算法的上述特点使得LMS滤波器越来越多的应用于更广泛的领域。附 录:clcclearN=1000; M=100; w=zeros(M,N,3,2); f=zeros(M,N,3,2);for l=1:2, if l=1 a=0.95; else a=-0.197; end; for d=1:3, if d=1 u=0.01; e
11、lse u=0.04*(d-1); end; for k=1:M, v=0.169*randn(1,N); x(1)=1; for n=2:N, x(n)=-a*x(n-1)+v(n); f(k,n,d,l)=x(n)-w(k,n,d,l)*x(n-1); w(k,n+1,d,l)=w(k,n,d,l)+u*x(n-1)*f(k,n,d,l); end end endendfor l=1:2 for d=1:3 for n=1:N, fea(n,d,l)=0; wea(n,d,l)=0; for m=1:M fea(n,d,l)=fea(n,d,l)+f(m,n,d,l)2; wea(n,d,
12、l)=wea(n,d,l)+w(m,n,d,l); end fea(n,d,l)=fea(n,d,l)/M; wea(n,d,l)=wea(n,d,l)/M; end endendn=1:N;figure(1)plot(n,w(1,n,2,1),'r-',n,wea(n,2,1),'b-',n,w(1,n,2,2),'r-',n,wea(n,2,2),'b-');title('一阶自适应预测器权值的瞬态特性(u=0.04)')xlabel('迭代次数'),ylabel('抽头权值')
13、legend('it单一实现','it100次实验集平均实现')text(400,-0.01,'a=-0.197')text(400,-0.8,'a=+0.95')figure(2)plot(n,f(1,n,2,1),'r-',n,fea(n,2,1),'b-');title('一阶自适应预测器的平方预测误差瞬时特性(u=0.04)')axis(0 1000 0.001 1);xlabel('迭代次数'),ylabel('平方误差')legend('it单次误差','it100次重复实验平均误差')figure(3)semilogy(n,fea(n,1,1),'b-',n,fea(n,
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 日本订购合同范例
- 大型装修工程合同范例
- 生鲜市场租赁合同范例
- 传媒公司会务合同范例
- 商业计划合同范例
- 车展模特合同范例
- 物业工程改造合同范例
- 油居间协议合同范例
- 堆垛货架租赁合同范例
- 电商线下商家合同范例
- 膳食营养课件教学课件
- 国开(内蒙古)2024年《创新创业教育基础》形考任务1-3终考任务答案
- 民族团结进步教育-中华民族共同体概论学习通超星期末考试答案章节答案2024年
- 2024入团知识题库(含答案)
- 2023年秋季国家开放大学-01467-土木工程力学(本)期末考试题带答案
- JCT908-2013 人造石的标准
- HCCDP 云迁移认证理论题库
- 工程结算单(样本)
- 日常物业管理服务流程图
- 洁净室内洁净度测试记录填写范例
- 铁建设110铁路建设项目施工图投资检算管理暂行
评论
0/150
提交评论