毕业论文维纳滤波器维纳霍夫方程课程设计.doc

《毕业论文维纳滤波器维纳霍夫方程课程设计.doc》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。

1、电气工程学院信号处理课程设计报告设计题目：维纳滤波器专业：生物医学工程指导教师:叶立夏学生姓名:叶立夏学号:20104320144起迄日期:2013年12月20日—2014年1月15日如果有做课程设计的同学不懂的，可以联系我！Qq：745581242目录前言………………………………………………31设计任务及指标…………………………………31.1课程设计的内容和要求………………41.2对课程设计成果的要求………………42设计思想…………………………………………42.1概述……………………………………42

2、.2主要仪器设备及耗材…………………53课程设计的组成部分……………………………53.1具体操作…………………………………54实验分析…………………………………………94.1原始图像显示…………………………94.2噪声的强度对维纳滤波器的影响……134.3阶数对滤波效果的影响………………174.4数据长度对维纳滤波的影响…………215设计总结………………………………………265.1思考题…………………………………265.2实验心得………………………………266主要参考文献…………………………………27

3、附录……………………………………………27前言去除信号中的噪声影响是信号处理中的一个重要内容，而滤波则是实现这一功能的重要手段之一。滤波器可以分为两类，及经典滤波器和现代滤波器经典滤波器是假定输入信号中嘚瑟有用成分和希望去除的成分各自占有不同的频带。当输入信号通过一个滤波器是可将欲去除的成分有效的去除，如果信号和噪声的频谱相互重叠，那么经典滤波器将无能为力。现代滤波器理论研究的主要内容是从含有噪声的数据记录中估计出信号的某些特征灬信号本身。一旦信号被估计出，那么估计出的信号的信噪比将比原信号的高。现代

4、滤波器把信号和噪声都视为随机信号，利用它们的统计特征导出一套最佳的估值算法，然后用硬件或软件予以实现。现代滤波器理论源于维纳在20世纪40年代及其以后的工作，因此维纳滤波器便是这一类滤波器的典型代表。维纳滤波器，也是最小平方滤波器，其基本思路为：设计一个滤波器，使其与输入信号滤波后的输出与期望输出在最小平方意义下的最佳逼近。寻求最小均方误差的实质就是解维纳-霍夫方程。1.设计任务及指标1.1课程设计的内容和要求（包括原始数据、技术要求、工作要求等）：本设计的目的是产生用于信号滤波的维纳霍夫方程。并且要

5、求调节该滤波器的参数使该滤波器能够最好的还原原始波形，以适应不同原始信号都能够被提取出来。设计要求：1.已知信号的自相关函数和噪声的能量，编写程序求解维纳-霍夫方程，寻找最优滤波器。2.编写程序仿真信号，噪声和观察波形，然后把观察信号通过滤波器得到的信号估计与原始信号比较，观察是否达到了去噪的目的。3.选择不同信号（仿真信号，实际采集的心电，脑电信号），人工添加噪声，调整噪声的相对强度，观察滤波效果。1.2.对课程设计成果的要求：利用MATLAB软件编写程序。根据维纳滤波的原理，通过调节参数获得效果最

6、佳的滤波器，用正弦衰减信号作为测试信号观察滤波效果，再将心电信号和脑电信号作为输入，观察滤波结果；调节滤波器的参数（噪声强度，阈值，数据长度）观察参数对滤波器的影响。2.设计思想2.1概述根据正交原理可以推导出维纳-霍夫方程，满足该方程的滤波器输出信号的估计值与信号在最小均方误差意义下最接近。（2-1）根据滤波器的形式，维纳滤波器可以分为三种情况：非因果IIR型，因果IIR型，FIR型，对于实时性有要求的情况下用后两种形式。图2.1维纳滤波器对于FIR型维纳滤波器，维纳-霍夫方程的形式为（2-2）或者

7、写成矩阵形式（2-3）其中（2-4）（2-5）（2-6）这样，如果信号和噪声的二阶统计特性已知，则易求解（2-7）维纳滤波的均方误差是（2-8）2.2主要仪器设备及耗材MATLAB软件，计算机3.课程设计的组成部分3.1具体操作3.1.1.建立一个以WH命名的函数（源程序上半部分）。function[h,e]=WH(Rss,Rww,M)%求解维纳-霍夫方程的函数，其中M为信号的长度e1=10;%给e1赋初值e0=0;%给e0赋初值N=0;%给N赋初值%以下循环的目的是找出FIR滤波器合适的阶数%判据是

8、当阶数增加而均方误差没有明显下降时，则认为阶数足够whileabs(e0-e1)>1e-6%e1和e0不够接近则循环，确定阈值N=N+1;%使N递增，使其逐步靠近需要满足的条件e0=e1;Rxs=Rss(M:(M+N-1));Rxx=Rww(M:(M+N-1))+Rss(M:(M+N-1));%输入信号总的自相关函数R_xx=zeros(N);%产生一个N*N的全零矩阵作为仿真信号的自相关函数forj=1:Nforn=1:NR_xx(j,n)=Rxx(a