《基于Matlab的常用滤波算法研究(含代码)_毕业设计(论文)(34页).doc》由会员分享，可在线阅读，更多相关《基于Matlab的常用滤波算法研究(含代码)_毕业设计(论文)(34页).doc（34页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、-基于Matlab的常用滤波算法研究(含代码)_毕业设计(论文)-第 32 页 毕业设计(论文) UNDERGRADUATE PROJECT (THESIS)题 目: 冲击测试常用滤波算法研究学 院 专 业 学 号 学生姓名 指导教师 起讫日期 目录摘要2ABSTRACT3第一章 绪论41.1课题背景41.2国内外相关领域的研究41.3主要研究内容与创新51.3.1研究内容与意义51.3.2课题的创新点51.3.3 研究目的与技术指标6第二章 数字滤波基础72.1数字滤波算法概念72.2数据采样与频谱分析原理82.2.1 时域抽样定理82.2.2 离散傅立叶变换(DFT)82.2.3 快速傅立
2、叶变换(FFT)92.2.4 频谱分析原理92.3常用数字滤波算法基础102.3.1常用数字滤波算法分类102.3.2常用数字滤波算法特点112.3.3常用滤波算法相关原理132.4 冲击测试采样数据162.4.1噪声的特点与分类162.4.2冲击测试采样数据特点172.5 MATLAB简介172.5.1 MATLAB功能简介182.5.2 MATLAB的发展18第三章、冲击测试滤波算法设计及滤波效果分析203.1 冲击测试采样数据的分析203.2 滤波算法设计及效果分析213.2.1 中位值平均法的设计213.2.2限幅法和限速法的设计233.2.3一阶滞后法的设计253.2.4低通法的设计
3、26第四章 结论与展望344.1冲击测试的滤波算法总结344.2冲击测试的滤波算法展望34致 谢36参考文献37附录：程序代码清单38冲击测试常用滤波算法研究摘要动态信号分析仪是一种主要应用在噪声、振动分析、模型分析、电子设计和声学测试的工具，冲击测试和冲击谱分析是确定设备在经受外力冲撞或作用时的安全性、可靠性和有效性的实验方法，也是动态信号分析仪的一项重要功能。冲击测试采样数据往往会受到来自环境中的各种噪声干扰，有必要对其进行滤波以更好地分析其真实的冲击谱特性。本文中主要研究了一些应用于冲击测试数据滤波的常用滤波算法如中位值法、算术平均法、中位值平均法、限幅法、限速法、一阶滞后法、低通法等。
4、MATLAB是用于算法开发、数据可视化、数据分析以及数值计算的高级技术计算语言和交互式环境。本文在MATLAB环境下实现各常用数字滤波算法并讨论了各滤波算法的特点及其选用原则。针对冲击测试仪采样数据的噪声干扰特点，挑选出合适的算法对冲击测试采样数据进行滤波，分析比较各滤波算法的实际滤波效果并且进行相关优化以实现更优的滤波效果。论文的主要内容：第一章、绪论主要介绍了有关动态信号分析仪的课题背景、国内外发展情况，课题的研究意义及创新点等。第二章、介绍了数字滤波的基本原理、常用滤波算法分类与特点、频谱分析基本原理和环境噪声的相关知识。第三章、主要介绍各滤波算法的参数选择和具体设计流程，并对各算法的滤
5、波效果进行分析比较。第四章、对冲击测试常用滤波算法实现和滤波效果的分析总结及前景展望。关键词：数字信号处理、数字滤波、滤波算法、MATLABABSTRACTDynamic signal analyzer is a tool mainly applied in noise, vibration analysis, model analysis, electronic design and acoustic test. Impact test and impact spectrum analysis are experimental methods that examine the safety,
6、 reliability and validity of the equipment under external force impact. Impact test sample data often suffers from all kinds of noises in the environment, it is necessary to filter the data for better and real impact spectrum characteristics.In this paper, we mainly study commonly used digital filte
7、ring algorithm such as median value method, the arithmetic average method, weighted average method, clipping method, limited speed method, and low pass method and so on. MATLAB is a computing language and interactive environment for senior technical computing algorithm development, data visualizatio
8、n, data analysis and numerical analysis.Using MATLAB, we simulate the commonly used digital filtering algorithms and discuss the characteristics and selection principles of the algorithms. According to the characteristics of noise of the data sampling impact tester, a suitable algorithm it is select
9、ed, analyzed, justified to achieve better filtering effect.The reminder of the thesis is as follows:Chapter one introduces the project background, the development of the dynamic signal analyzer domestic and abroad, the research significance and innovations etc.Chapter two includes the basic principl
10、e of digital filtering, classification and characteristics of filtering algorithm commonly used, basic principle of spectrum analysis and the relevant knowledge of environmental noise.Chapter three is the design and effect analysis of digital filtering algorithm in details.The last chapter is the co
11、nclusion of the design and effect analysis of impact test filtering algorithms. Keywords: Digital signal processing; Digital filtering; Filtering algorithm; MATLAB第一章 绪论1.1课题背景动态信号分析仪是从表示物理量的电信号中分析其特性参数的仪器，由硬件和分析软件构成。动态信号分析仪可从时域、频域和幅值域分析被测信号，具有功能全、分析速度快、测量参数多、频率分辨力和幅值精度高等特点，广泛应用于如计算机制造、航空航天、电子、军事、生物
12、医学、通讯等高科技产品的质量检测和诊断分析1。 动态测试、信号处理、模态试验、环境试验、状态监测、故障诊断的核心手段是动态信号分析仪2。 而随着科技的进步，特别是微电子技术的迅猛发展，嵌入式微控制器的出现与使用，诞生了集动态信号采集、分析、存储、显示为一体的动态信号分析仪。该种动态信号分析仪体积小、重量轻、速度快、功能强大、方便携带于工程现场的使用，需要时可将所存的数据传送到计算机进行更详细的分析。冲击测试一般是确定军民用设备在经受外力冲撞或作用时产品的安全性、可靠性和有效性的一种试验方法。而冲击响应谱通常简称“冲击谱”，它是将受到机械冲击作用的一系列单自由度系统的最大响应(如位移、速度或加速
13、度)响应值随系统的固有频率而变化的频谱3。国家标准化组织(ISO)所属的技术委员会以及我国的国家标准，都已经把冲击谱作为规定冲击环境的方法之一。冲击谱是对设备实施抗冲击设计的分析基础，也是控制产品冲击环境模拟实验的基本参数。因此冲击响应谱分析是动态信号分析仪应具备的一项重要功能。由于在冲击测试中，采样数据不可避免地受到各种环境噪声的干扰。所以在进行冲击响应谱分析前，对冲击测试采样数据进行数字滤波处理是很有必要的。本文主要研究冲击测试的常用数字滤波算法，数字滤波根据有用信号与噪声的不同特性，消除或减弱噪声。它对信号安全可靠和有效灵活地传递是至关重要的4。数字滤波方式具有精度高、可靠性高、灵活易用
14、(可程控改变特性)、便于集成等特点。数字滤波是语音处理、图像处理、模式识别、频谱分析等应用的基本处理算法5。语音处理是最早应用数字滤波的领域之一，也是最早推动数字信号处理理论发展的领域之一。语音的去噪与增强技术已取得许多成果。目前，数字信号滤波在图像处理、数据压缩等方面取得了巨大的进展和成就。小波理论由于其局部分析性能的优异在图像处理中的应用研究得到迅速发展，尤其是在图像压缩、图像去噪等方面的应用研究6。而在数字通信、网络通信、图像通信、多媒体通信等应用中，离开了数字滤波几乎是寸步难行。1.2国内外相关领域的研究动态信号分析仪在电子测量领域中被称为频域中的“射频万用表”，具有较高的实用性，并得
15、到了广泛的应用7。它同时具备几种仪器的功能，坚固、轻便、是适用于现场应用的理想仪器，其性能和功能可适应有严格要求的研发应用需要。内置的信号源及可选的特性更优化了仪器用于分析和查找噪音、振动与声学问题，评测控制系统的功能，以及评估和解决了旋转机器问题，并定性与评估控制系统参数。近几年，随着数字信号处理器(DSP)和PC机的发展,先进的动态信号分析仪都采用数据采集处理系统加PC机显示管理程序来实现。国内外对动态信号分析仪的研究尤其是国内有了较快的发展。国外多通道动态信号分析仪的技术已经基本成熟。产品的体积小，速度快，功能强大，性能优异，操作简单但是价格比较昂贵。主要研制的公司有:美国惠普公司(HP
16、)，美国安捷伦公司(Agilent惠普公司重组后的测量公司)，美国尼高力仪器技术公司(Nioclet)，美国Dactorn公司和美国国家仪器有限公司(NI)，惠普公司和安捷伦公司的产品比较早，功能上比较单一，外形类似一般的示波器。如：HP35660A、HP35665A、 HP35670A、Agilent 35670A。NI公司的产品硬件是一块动态信号采集板，而上层则是使用它的Labview软件开发的虚拟仪器。而Dactron和 Nicolet的动态信号分析仪都采用了最新的技术，一般都是采集处理系统再配合PC机的上层管理软件PC机一般工作Windows操作系统，操作简单，功能强大。如 Dactr
17、on 公司：Photon便携式动态信号分析仪、Spectra动态信号分析仪、Book Focus实时动态信号分析仪Nicolet 公司：Horizon、 Compass动态信号分析仪8。国内自主研发了一批动态信号分析仪器，以北京万博振通公司的产品为代表。北京万拨通公司生产的BVM-102高档触摸型双通道信号分析仪，采用5.6寸液晶显示屏，在Windows xp平台下运行的控制管理软件，信号频率范围高达200kHz，体积小(190*150*64mm)。国内的动态信号分析仪大体功能与国外接近，但是国外产品由于技术比较成熟，所以仪器各方面指标都比较高，另外国外厂商在注重仪器本身的研发的同时注重外围维
18、护和管理软件环境的研发，产品附加值较高7。1.3主要研究内容与创新1.3.1研究内容与意义数字滤波是常见的信号处理分析方法。本课题主要讨论应用于冲击测试数据的一些常用数字滤波算法如中位值法、算术平均法、中位值平均法、限幅法、限速法、一阶滞后法、低通法等。基于MATLAB环境实现了各常用数字滤波算法并分析、讨论各算法的特点及选用原则。针对冲击测试采样数据的噪声特点，挑选出适合的算法对冲击测试采样数据进行滤波，分析各滤波算法的滤波效果并且进行相关优化以实现更优的滤波效果。1.3.2课题的创新点课题的创新点主要有三点：第一、算法应用于实际环境的冲击测试采样数据，其中所受的噪声干扰如脉冲干扰、工频干扰
19、具有代表性。对它进行的滤波算法研究有一定实际价值。第二、算法基于MATLAB实现，充分利用了MATLAB的强大计算仿真功能。MATLAB作为一款强大的数学、金融、通信、信号处理的交互式软件，我们可以很方便的在其中编写程序、仿真模型。第三、综合实现了各常用数字滤波算法，总结了各个算法的滤波特点与选用原则。 1.3.3 研究目的与技术指标研究目的：毕业设计是一个全面运用各方面知识与能力的综合实践。希望能够通过此次毕业设计能使得我们加深对常用数字滤波知识的理解并尽可能完成滤波算法的优化，同时也使自己掌握分析解决问题的方法，提高自己解决实际问题的能力。技术指标：1、冲击测试滤波算法的设计方法，即数字滤
20、波算法的研究与比较，从中选出适合冲击测试系统的滤波算法。2、冲击测试滤波算法的效果分析，即对某算法应用于具体采样数据的滤波效果分析，进行相关调试或优化以达到更加理想的效果。第二章 数字滤波基础2.1数字滤波算法概念信号分析与处理中，常常会遇到有用信号中叠加了无用噪声的问题。这些噪声有的是与信号同时产生的，有时是传输过程中混入的，有时噪声会大于有用的信号，从而淹没掉信号。因此，从接收到的信号中消除或减弱干扰噪声，就成为信号传输与处理中十分重要的问题。根据有用信号与噪声的不同特性，消除或减弱噪声，提取有用的信号的过程称为滤波9。滤波技术是信号分析处理的重要分支。无论是信号的获取、传输还是信号的处理
21、和交换都离不开滤波技术，它对信号安全可靠和灵活有效地传递是至关重要的。数字信号处理 (Digital Signal Processing)是一门涉及多学科并广泛于各个科学和工程领域的新兴学科5。它是将信号以数字方式表示并处理的理论和技术。数字信号处理与模拟信号处理是信号处理的子集。现代社会已经进入数字化、信息化时代，DSP技术因此得到广泛应用和不断发展。而数字滤波作为其重要分支也飞速发展。数字滤波是用数字处理的方法按预定的要求滤除干扰信号，获得有用信号。其最初的设想是代替模拟滤波10。不同于模拟滤波通过无源或有源器件搭建模拟电路的滤波方式，数字滤波是由数字乘法器、加法器和延时单元组成的一种计算
22、方法。其功能是对输入离散信号的数字代码进行运算处理，以达到改变信号频谱的目的。由于电子计算机技术和大规模集成电路的发展，数字滤波已可用计算机软件实现，也可用大规模集成数字硬件实时实现。数字滤波是一个离散时间系统，应用数字滤波处理模拟信号时，首先须对输入模拟信号进行限带、抽样和模数转换。数字滤波输入信号的抽样率应大于被处理信号带宽的两倍，其频率响应具有以抽样频率为间隔的周期重复特性，且以折叠频率即1/2抽样频率点呈镜像对称。为得到模拟信号，数字滤波处理的输出数字信号须经数模转换、平滑。数字滤波方式具有精度高、可靠性高、灵活易用(可程控改变特性)、便于集成等特点。精度高： 数字滤波方式具有比模拟滤
23、波器更高的精度，甚至能够实现后者在理论上也无法达到的性能。例如，对于数字滤波器来说很容易就能够做到一个 1000Hz 的低通滤波器允许 999Hz 信号通过并且完全阻止 1001Hz 的信号，模拟滤波器无法区分如此接近的信号。可靠性高： 数字滤波方式还具有模拟滤波方式不能比拟的可靠性。传统的模拟滤波器由电子元件组成，其电路特性会随着时间、温度、电压的变化而漂移，而数字电路则没有这种问题，只要在数字电路的工作环境下，数字滤波器就能够稳定、可靠的持续工作。灵活易用：数字滤波方式有低通、高通、带通、带阻和全通等类型。它可以是时不变的或时变的、因果的或非因果的、线性的或非线性的。应用最广的是线性、时不
24、变数字滤波器，以及FIR数字滤波器，而这些设计都不需改变硬件，只要修改参数和程序就可实现。便于集成： 数字滤波器一般由寄存器、延时器、加法器和乘法器等基本数字电路实现。随着大规模集成电路技术的成熟，数字滤波器可以很方面的集成与实现，其应用领域也因此越来越广。数字滤波是语音处理、图像处理、模式识别、频谱分析等应用的基本处理算法。离开了数字滤波几乎是寸步难行5。2.2数据采样与频谱分析原理2.2.1 时域抽样定理11时域抽样定理给出了连续信号抽样过程中信号不失真的约束条件：对于基带信号，信号抽样频率大于等于2倍的信号最高频率，即。时域抽样是把连续信号变成适于数字系统处理的离散信号。对连续信号以间隔
25、抽样，则可得到的离散序列为。图2-1 连续信号抽样的离散序列若，则信号与的频谱之间存在： (2-1)其中,的频谱为,的频谱为。可见，信号时域抽样导致信号频谱的周期化。为抽样角频率，为抽样频率。数字角频率与模拟角频率的关系为：=T。2.2.2 离散傅立叶变换(DFT)12有限长序列的离散傅立叶变换(DFT)为 (2-2)逆变换为 (2-3) 2.2.3 快速傅立叶变换(FFT)13在各种信号序列中，有限长序列占重要地位。对有限长序列可以利用离散傅立叶变换(DFT)进行分析。DFT不但可以很好的反映序列的频谱特性，而且易于用快速算法(FFT)在计算机上进行分析。有限长序列的DFT是其z变换在单位圆
26、上的等距离采样，或者说是序列傅立叶的等距离采样，因此可以用于序列的谱分析。FFT是DFT的一种快速算法，它是对变换式进行一次次分解，使其成为若干小数据点的组合，从而减少运算量。MATLAB为计算数据的离散快速傅立叶变换，提供了一系列丰富的数学函数，主要有fft、ifft、fft2 、ifft2, fftn、ifftn和fftshift、ifftshift等。当所处理的数据的长度为2的幂次时，采用基-2算法进行计算，计算速度会显著增加。所以，要尽可能使所要处理的数据长度为2的幂次或者用添零的方式来添补数据使之成为2的幂次。fft函数调用方式：Y=fft(X)；Yfft(X,N)；Yfft(X,d
27、im)或Yfft(X,N,dim)。函数ifft的参数应用与函数fft完全相同。2.2.4 频谱分析原理14时域分析只能反映信号的幅值随时间的变化情况，除单频率分量的简单波形外，很难明确提示信号的频率组成和各频率分量大小，而频谱分析能很好的解决此问题。由于从频域能获得的主要是频率信息，所以本节主要介绍频率(周期)的估计与频谱图的生成。1、频率、周期的估计对于Y(kf)，如果当kf=f时，Y(kf)取最大值，则f为频率的估计值，由于采样间隔的误差，f也存在误差，其误差最大为f/ 2。周期T=1/f。从原理上可以看出，如果在标准信号中混有噪声，用上述方法仍能够精确地估计出原标准信号的频率和周期，这
28、个将在下一章做出验证2、频谱图为了直观地表示信号的频率特性，工程上常常将Fourier变换的结果用图形的方式表示，即频谱图。以频率f为横坐标，|Y(f)|为纵坐标，可以得到幅值谱；以频率f为横坐标，argY(f)为纵坐标，可以得到相位谱；以频率f为横坐标，ReY(f)为纵坐标，可以得到实频谱；以频率f为横坐标，ImY(f)为纵坐标，可以得到虚频谱。根据采样定理，只有频率不超过Fs/2的信号才能被正确采集，即Fourier变换的结果中频率大于Fs/2的部分是不正确的部分，故不在频谱图中显示。即横坐标f0,Fs/2。2.3常用数字滤波算法基础2.3.1常用数字滤波算法分类图2-2 数字滤波算法分类
29、 图中为数字滤波算法的分类，包括经典滤波算法和现代滤波算法。本文研究的冲击测试数据滤波对于滤波要求相对简单，所以文中主要讨论几种经典数字滤波算法。 2.3.2常用数字滤波算法特点中位值法 A、方法： 连续进行N次数据采样(N一般取奇数) 把采样的N个数据按大小排列 取数列的中间值为本次有效值 B、优点： 能有效滤除由于偶尔因素或采样器本身不稳定而引起的脉冲干扰 对温度、液位等变化缓慢的被测参数比较有效 C、缺点： 对流量、速度等快速变化的参数则不宜采用 对周期性的随机干扰效果一般算术平均法 A、方法： 连续进行N次采样并对采样得到的这N个数据进行算术平均运算 把计算得到的算术平均值作为本次有效
30、值 当N值较大时：信号平滑度较高，但灵敏度较低 当N值较小时：信号平滑度较低，但灵敏度较高 B、优点： 对具有幅值变化不大的周期随机干扰滤波效果好，对毛刺有平滑作用 受随机干扰的信号的特点是有一个平均值，信号在某一数值附近上下波动 C、缺点： 对于测量速度较慢或要求数据计算速度较快的实时控制不适用 对脉冲干扰的滤波效果一般中位值平均滤波法(又称防脉冲干扰平均滤波法) A、方法： 相当于“中位值滤波法”+“算术平均滤波法” 连续采样N个数据，去掉一个最大值和一个最小值 然后计算N-2个数据的算术平均值作为本次的有效值 N值的选取：314 B、优点： 同时具有两种滤波法的优点 对于脉冲噪声干扰与幅
31、值变化不大的随机干扰都有良好的抑制作用，可消除干扰引起的采样值偏差 C、缺点： 测量速度比较慢，和算术平均滤波法一样 实际滤波效果可能不及“中位值滤波法”+“算术平均滤波法”的方式。限幅法 A、方法： 根据经验常识判断，确定两次采样允许的最大偏差值(设为y) 每次检测到新值时判断： 如果|y(n)-y(n-1) |y, 则取y(n)=y(n-1) B、优点： 能有效克服因偶然因素引起的脉冲干扰 C、缺点 无法抑制周期性的随机干扰限幅差值难确定，平滑度差限速法 A、方法： 根据经验判断, 确定最大偏差值y，最多取两次采样值并进行判断以得到本次采样的有效值 每次检测到新值时判断： 如果|y(n)-
32、y(n-1) | y, 则y(n)无效,但仍保留y(n),继续采样取得y(n+1) 如果|y(n+1)-y(n) | y, 则y(n) = B、优点： 能克服因偶然因素引起的脉冲干扰 既照顾了采样的实时性，又顾及了采样值变化的连续性 C、缺点 增加了较多的计算量，处理速度变慢 同样无法抑制周期性的随机干扰 差值难确定，平滑度差一阶滞后滤波法 A、方法： 一阶滞后滤波也叫一阶惯性滤波 取滞后系数a=01 本次滤波结果=(1-a)本次采样值+a上次滤波结果 B、优点： 对周期性干扰具有良好的抑制作用 效果相当硬件电路一阶惯性环节适用于波动频率较高的场合 C、缺点： 相位滞后，灵敏度低 滞后程度取决
33、于a值大小低通法 A、方法： 根据噪声的频率设定一个合适的截止频率wc 按wc及其他指标设计相应的低通数字滤波器 按频率滤波，使有用信号所在频率通过、噪声所在频率被滤除 B、优点： 能有效滤波截止频率范围外的噪声，且平滑度较高 参数易调整，可以方便地进行修改以达到最佳效果 C、缺点： 无法滤除与有用信号混频的噪声干扰无法滤除高于采样频率一半的干扰2.3.3常用滤波算法相关原理中位值平均滤波法：中位值滤波：对连续采样的N个数据进行排序，取排序在中间的值。算术平均滤波：连续进行N次采样并对采样得到的N个数据进行算术平均运算。中位值平均滤波法：综合了两者优势的算法，算法虽然简单，但是它是基于数学概率
34、的基本原理，其滤波效果是不错的。限幅法、限速法：实际信号中的原始信号因为信号的连续性，往往存在相邻两个采样值的一个最大偏差差值，超过这个偏差值基本就可以确定是噪声干扰。限幅法就是根据经验常识判断，确定两次采样允许的最大偏差值(设为y)，进行简单的程序判断，而限速法是在权衡采样实时性和采样变化的连续性后对限幅法的有限改进。一阶滞后算法：15常用的滤波器的传递函数是： (2-4) 其中，滤波器的滤波效果取决于滤波时间常数。可以模仿上式做成一阶惯性滤波器亦称低通滤波器，即将上式写成差分方程： (2-5)整理可得： (2-6)FIR滤波器相关知识： 16FIR滤波器的主要设计方法有窗函数法、最优化设计
35、法及约束最小二乘逼近法。在滤波器传统设计中，要得到其幅频和相频响应特性，需要根据这些方法进行大量的计算，这使得滤波器的设计缓慢，周期变长，不利于设计的优化。MATLAB信号处理工具箱中提供了基于滤波器设计方法的工具函数，编程中可根据设计要求直接调用相应的函数，方便快捷，FIR数字滤波器的设计方法主要是建立在对理想滤波器频率特性作某种近似的基础上，这些近似方法较多。本文选用窗函数法来设计FIR滤波器。窗函数设计滤波器的基本思想,是把给定的频率响应通过IDTFT(Inverse Discrete Time Fourier Transform)，求得脉冲响应，然后利用加窗函数对它进行截断和平滑，以实
36、现一个物理可实现且具有线性相位的FIR滤波器的设计目的。其核心是从给定的频率特性，通过加窗确定有限长单位脉冲响应序列。系统单位脉冲序列，它是一个有限长序列。假设是所要求的理想响应序列，则 (2-7)式(2-7)中是对应的单位脉冲响应序列，而滤波器的频率响应和单位脉冲响应序列是傅里叶变换对，则 (2-8)求得序列后，可得到 (2-9) 注意到，这里为无限长序列，因此是物理不可实现的。为了使系统变为物理可实现的，且使FIR滤波器实际频率响应尽可能逼近理想的频率响应，采用窗函数将无限脉冲响应截取一段来近似表示可得 (2-10) 由此可得 (2-11) 式(2-11)中，为窗口宽度，是物理可实现系统。
37、而窗函数所必须满足的特性：(1)窗谱主瓣尽可能地窄，以获得较陡的过渡带；(2)尽量减少窗谱的最大旁瓣的相对幅度，也就是能量尽量集中于主瓣，使肩峰和纹波减小，就可增多阻带的衰减。(3)窗函数的频率特性的旁瓣在当趋近于的过程中，其能量迅速趋于零。这就给窗函数序列的形状和长度选择提出了严格的要求。常用窗函数有如下几种8：矩形窗(Rectangle Window)、三角窗(Bartlett Window)、汉宁窗(Hanning Window)、海明窗(Hamming Window)、布莱克曼窗(Blackman Window)、凯泽窗(Kaiser Window)( =7865)，各种窗函数基本参数
38、比较如下表:表1 常用窗函数基本参数窗函数频谱性能指标加窗后滤波器性能指标旁瓣峰值主瓣宽度过渡带宽阻带最小衰减/dB/(2/N)/(2/N)/ dB矩形窗-1320.9-21三角形窗-2542.1-25汉宁窗-3143.1-44海明窗-4143.3-53布莱克曼窗-5765.5-74凯泽窗(=7.865)-575-80IIR滤波器相关知识 17IIR数字滤波器是一种离散时间系统，其系统函数为 (2-12) 假设MN，当MN时,系统函数可以看作一个IIR的子系统和一个(M-N)的FIR子系统的级联。IIR数字滤波器的设计实际上是求解滤波器的系数和，它是数学上的一种逼近问题，即在规定意义上(通常采
39、用最小均方误差准则)去逼近系统的特性。如果在S平面上去逼近，就得到模拟滤波器；如果在z平面上去逼近，就得到数字滤波器。模拟滤波器技术指标变换数字滤波器技术指标模拟滤波器设计方法变换数字滤波器H(Z)模拟滤波器H(S)S=f(Z)图2-3 IIR数字滤波器的设计步骤流程图IIR数字滤波器的设计可以使用模拟滤波器原型，借鉴成熟的模拟滤波器设计方法进行设计再转化成数字滤波器。常用的设计方法有脉冲响应不变法和双线性法，本文主要使用脉冲响应不变法进行设计。而常用的模拟滤波器原型有：巴特沃兹滤波器、切比雪夫滤波器、椭圆滤波器等，其幅频特性可逼近于理想滤波器的幅频响应特性。下图分别为巴特沃兹滤波器、切比雪夫
40、滤波器 I 型和II型、椭圆滤波器的幅频响应特性图。图2-4 常用滤波器幅频响应特性图2.4 冲击测试采样数据2.4.1噪声的特点与分类(1) 根据干扰的来源划分内部干扰指系统内部电子电路的各种干扰。外部干扰指由外界窜入到系统的各种干扰。(2) 根据干扰出现的规律划分固定干扰指系统附近固定的电气设备运行时发出的干扰。半固定干扰指某些偶然使用的电气设备(如行车、电钻)引起的干扰。随机干扰属于偶发性的干扰。如闪电、供电系统继电保护的动作等干扰。(3)根据干扰产生和传播的方式划分静电干扰磁场耦合干扰电磁辐射干扰 电导通路耦合干扰漏电耦合干扰(4) 根据干扰输入信号的关系划分串模(差模)干扰共模干扰其
41、中，根据干扰产生和传播的方式划分： 静电干扰静电干扰实际是电场通过电容耦合的干扰。从电路理论可知，电流流经一导体时，导体产生电场，这个电场可交连到附近的导体中，使它们感生出电位，这个电位就是干扰电压。从交流电路传输来看，干扰起因于导体与导线之间、元件之间的寄生电容。外部噪声源与导体之间的寄生电容耦合到电路，造成对电路的干扰。 磁场耦合干扰磁场耦合干扰是一种感应干扰。在连接信号源的传输线经过的空间总存在着交变电磁场。在诸如动力线、变压器、电动机、继电器、电风扇等附近，都会有这种磁场。这些交变的磁场穿过传输线形成的回路，将在传输线上或闭合导线上感应出交流干扰电压。 电磁辐射干扰在工厂内，各种大功率
42、高频、中频发生装置以及各种电火花机床，都将产生高频电磁波向周围空间辐射，形成电磁辐射干扰源。辐射能量是以与通信接收机接收无线电频率能量相同的方法耦合到电路中而产生干扰。 电导通路耦合干扰电导通路是指构成电回路的通路。电导通路耦合干扰是由各单元回路之间的公共阻抗干扰。由于接地电位不同而造成的干扰为这类干扰的主要表现形式。在数据采集系统中，“地”有两种含义：一是指大地，它是系统中各个设备的自然参考电位，二是指一个设备内部电源的参考电位。如果一个仪器的地线不与大地连接，则称为“浮地”，否则称为接地。 漏电耦合干扰漏电耦合干扰是由于仪器内部的电路绝缘不良而出现的漏电流引起的电阻耦合产生的干扰，如图所示
43、。在高输入阻抗器件组成的系统中，其阻抗与电路板绝缘电阻可以比拟，通过电路板产生漏电流，将形成干扰。2.4.2冲击测试采样数据特点作为冲击测试采样的电压数据，其在采样和传输过程中不可避免地会受到环境和系统电路本身的噪声干扰，噪声的干扰对冲击测试结果会产生许多不利的影响，我们有必要其进行数字滤波处理。根据冲击测试结果可知，冲击测试采样数据受到的主要干扰为电工噪声干扰和随机噪声干扰.。其中的随机噪声干扰主要是由数量少幅值变化大的脉冲干扰和数量多幅值变化小的随机干扰组成，而这些随机干扰普遍存在于环境中。电工干扰则属于波动频率较高的周期性干扰，是市电以电磁波的辐射形式对周围的电气设备和电子设备造成干扰，
44、可以导致设备运行异常，严重时甚至损坏设备。2.5 MATLAB简介 2.5.1 MATLAB功能简介MATLAB是一种科学计算软件。MATLAB 是矩阵实验室(matrix laboratory)的缩写。这是一种以矩阵为基础的交互式程序计算语言。早期的MATLAB主要用于解决科学和工程的复杂数学计算问题。由于它使用方便、输入便捷、运算高效、适应科技人员的思维方式，并且有绘图功能，有用户自行扩展的空间，特别受到用户的欢迎，使它成为在科技界广为使用的软件，也是国内外高校教学和科学研究的常用软件。MATLAB比较易学，它只有一种数据类型(即64位双精度二进制)，一种标准的输入输出语句，它用解释方式工作，不需要编译，一般入门后经过自学就可以掌握。如果有不清楚的地方，可