毕业设计(论文)
便、是适用于现场应用的理想仪器，其性能和功能可适应有严格要求的研发应用需要。内置的信号源及可选的特性更优化了仪器用于分析和查找噪音、振动与声学问题，评测控制系统的功能，以及评估和解决了旋转机器问题，并定性与评估控制系统参数。近几年，随着数字信号处理器(DSP)和PC机的发展,先进的动态信号分析仪都采用数据采集处理系统加PC机显示管理程序来实现。国内外对动态信号分析仪的研究尤其是国内有了较快的发展。
国外多通道动态信号分析仪的技术已经基本成熟。产品的体积小，速度快，功能强大，性能优异，操作简单但是价格比较昂贵。主要研制的公司有:美国惠普公司(HP)，美国安捷伦公司(Agilent惠普公司重组后的测量公司)，美国尼高力仪器技术公司(Nioclet)，美国Dactorn公司和美国国家仪器有限公司(NI)，
惠普公司和安捷伦公司的产品比较早，功能上比较单一，外形类似一般的示波器。如：HP35660A、HP35665A、 HP35670A、Agilent 35670A。
NI公司的产品硬件是一块动态信号采集板，而上层则是使用它的Labview软件开发的虚拟仪器。
而Dactron和 Nicolet的动态信号分析仪都采用了最新的技术，一般都是采集处理系统再配合PC机的上层管理软件PC机一般工作Windows操作系统，操作简单，功能强大。如 Dactron 公司：Photon便携式动态信号分析仪、Spectra动态信号分析仪、Book Focus实时动态信号分析仪Nicolet 公司：Horizon、 Compass动态信号分析仪[8]。
国内自主研发了一批动态信号分析仪器，以北京万博振通公司的产品为代表。北京万拨通公司生产的BVM-102高档触摸型双通道信号分析仪，采用5.6寸液晶显示屏，在Windows xp平台下运行的控制管理软件，信号频率范围高达200kHz，体积小(190*150*64mm)。国内的动态信号分析仪大体功能与国外接近，但是国外产品由于技术比较成熟，所以仪器各方面指标都比较高，另外国外厂商在注重仪器本身的研发的同时注重外围维护和管理软件环境的研发，产品附加值较高[7]。
1.3主要研究内容与创新
1.3.1研究内容与意义
数字滤波是常见的信号处理分析方法。本课题主要讨论应用于冲击测试数据的一些常用数字滤波算法如中位值法、算术平均法、中位值平均法、限幅法、限速法、一阶滞后法、低通法等。基于MATLAB环境实现了各常用数字滤波算法并分析、讨论各算法的特点及选用原则。针对冲击测试采样数据的噪声特点，挑选出适合的算法对冲击测试采样数据进行滤波，分析各滤波算法的滤波效果并且进行相关优化以实现更优的滤波效果。
1.3.2课题的创新点
课题的创新点主要有三点：第一、算法应用于实际环境的冲击测试采样数据，其中所受的噪声干扰如脉冲干扰、工频干扰具有代表性。对它进行的滤波算法研究有一定实际价值。第二、算法基于MATLAB实现，充分利用了MATLAB
5
毕业设计(论文)
的强大计算仿真功能。MATLAB作为一款强大的数学、金融、通信、信号处理的交互式软件，我们可以很方便的在其中编写程序、仿真模型。第三、综合实现了各常用数字滤波算法，总结了各个算法的滤波特点与选用原则。
1.3.3 研究目的与技术指标 研究目的：
毕业设计是一个全面运用各方面知识与能力的综合实践。希望能够通过此次毕业设计能使得我们加深对常用数字滤波知识的理解并尽可能完成滤波算法的优化，同时也使自己掌握分析解决问题的方法，提高自己解决实际问题的能力。
技术指标：
1、冲击测试滤波算法的设计方法，即数字滤波算法的研究与比较，从中选出适合冲击测试系统的滤波算法。
2、冲击测试滤波算法的效果分析，即对某算法应用于具体采样数据的滤波效果分析，进行相关调试或优化以达到更加理想的效果。
6
毕业设计(论文)
第二章 数字滤波基础
2.1数字滤波算法概念
信号分析与处理中，常常会遇到有用信号中叠加了无用噪声的问题。这些噪声有的是与信号同时产生的，有时是传输过程中混入的，有时噪声会大于有用的信号，从而淹没掉信号。因此，从接收到的信号中消除或减弱干扰噪声，就成为信号传输与处理中十分重要的问题。根据有用信号与噪声的不同特性，消除或减弱噪声，提取有用的信号的过程称为滤波[9]。滤波技术是信号分析处理的重要分支。无论是信号的获取、传输还是信号的处理和交换都离不开滤波技术，它对信号安全可靠和灵活有效地传递是至关重要的。
数字信号处理 (Digital Signal Processing)是一门涉及多学科并广泛于各个科学和工程领域的新兴学科[5]。它是将信号以数字方式表示并处理的理论和技术。数字信号处理与模拟信号处理是信号处理的子集。现代社会已经进入数字化、信息化时代，DSP技术因此得到广泛应用和不断发展。而数字滤波作为其重要分支也飞速发展。
数字滤波是用数字处理的方法按预定的要求滤除干扰信号，获得有用信号。其最初的设想是代替模拟滤波[10]。不同于模拟滤波通过无源或有源器件搭建模拟电路的滤波方式，数字滤波是由数字乘法器、加法器和延时单元组成的一种计算方法。其功能是对输入离散信号的数字代码进行运算处理，以达到改变信号频谱的目的。由于电子计算机技术和大规模集成电路的发展，数字滤波已可用计算机软件实现，也可用大规模集成数字硬件实时实现。数字滤波是一个离散时间系统，应用数字滤波处理模拟信号时，首先须对输入模拟信号进行限带、抽样和模数转换。数字滤波输入信号的抽样率应大于被处理信号带宽的两倍，其频率响应具有以抽样频率为间隔的周期重复特性，且以折叠频率即1/2抽样频率点呈镜像对称。为得到模拟信号，数字滤波处理的输出数字信号须经数模转换、平滑。
数字滤波方式具有精度高、可靠性高、灵活易用(可程控改变特性)、便于集成等特点。
精度高： 数字滤波方式具有比模拟滤波器更高的精度，甚至能够实现后者在理论上也无法达到的性能。例如，对于数字滤波器来说很容易就能够做到一个 1000Hz 的低通滤波器允许 999Hz 信号通过并且完全阻止 1001Hz 的信号，模拟滤波器无法区分如此接近的信号。
可靠性高： 数字滤波方式还具有模拟滤波方式不能比拟的可靠性。传统的模拟滤波器由电子元件组成，其电路特性会随着时间、温度、电压的变化而漂移，而数字电路则没有这种问题，只要在数字电路的工作环境下，数字滤波器就能够稳定、可靠的持续工作。
灵活易用：数字滤波方式有低通、高通、带通、带阻和全通等类型。它可以是时不变的或时变的、因果的或非因果的、线性的或非线性的。应用最广的是线性、时不变数字滤波器，以及FIR数字滤波器，而这些设计都不需改变硬件，只要修改参数和程序就可实现。
7
毕业设计(论文)
便于集成： 数字滤波器一般由寄存器、延时器、加法器和乘法器等基本数字电路实现。随着大规模集成电路技术的成熟，数字滤波器可以很方面的集成与实现，其应用领域也因此越来越广。
数字滤波是语音处理、图像处理、模式识别、频谱分析等应用的基本处理算法。离开了数字滤波几乎是寸步难行[5]。
2.2数据采样与频谱分析原理
2.2.1 时域抽样定理[11]
时域抽样定理给出了连续信号抽样过程中信号不失真的约束条件：对于基带信号，信号抽样频率fsam大于等于2倍的信号最高频率fm，即fsam?2fm。时域抽样是把连续信号X(t)变成适于数字系统处理的离散信号X[k]。对连续信号
X(t)以间隔T抽样，则可得到的离散序列为X[k]?X(kT)t?kT。
x(t)0T2Tx[k]t012k 图2-1 连续信号抽样的离散序列
若X[k]?X(kT)t?kT，则信号X(t)与X[k]的频谱之间存在：
X(ej?1)?Tn?????X?j(??n?sam)? (2-1)
其中,X(ej?)的频谱为X(j?),X[k]的频谱为X(ej?)。
π/T为抽样角频率，可见，信号时域抽样导致信号频谱的周期化。?sam?2fsam?1T为抽样频率。数字角频率Ω与模拟角频率ω的关系为：Ω=ωT。
2.2.2 离散傅立叶变换(DFT)[12]
8
毕业设计(论文)
有限长序列x(n)的离散傅立叶变换(DFT)为
?knX(k)?DFT[x(n)]??x(n)WN,0?n?N?1 (2-2)
n?0N?1
逆变换为
1x(n)?IDFT[X(k)]?N?knX(k)W,0?n?N?1?Nn?0N?1 (2-3)
2.2.3 快速傅立叶变换(FFT)[13]
在各种信号序列中，有限长序列占重要地位。对有限长序列可以利用离散
傅立叶变换(DFT)进行分析。DFT不但可以很好的反映序列的频谱特性，而且易于用快速算法(FFT)在计算机上进行分析。
有限长序列的DFT是其z变换在单位圆上的等距离采样，或者说是序列傅立叶的等距离采样，因此可以用于序列的谱分析。FFT是DFT的一种快速算法，它是对变换式进行一次次分解，使其成为若干小数据点的组合，从而减少运算量。
MATLAB为计算数据的离散快速傅立叶变换，提供了一系列丰富的数学函数，主要有fft、ifft、fft2 、ifft2, fftn、ifftn和fftshift、ifftshift等。当所处理的数据的长度为2的幂次时，采用基-2算法进行计算，计算速度会显著增加。所以，要尽可能使所要处理的数据长度为2的幂次或者用添零的方式来添补数据使之成为2的幂次。 fft函数调用方式： ○1Y=fft(X)； ○2Y＝fft(X,N)； ○3Y＝fft(X,[],dim)或Y＝fft(X,N,dim)。 函数ifft的参数应用与函数fft完全相同。
2.2.4 频谱分析原理[14]
时域分析只能反映信号的幅值随时间的变化情况，除单频率分量的简单波形外，很难明确提示信号的频率组成和各频率分量大小，而频谱分析能很好的解决此问题。由于从频域能获得的主要是频率信息，所以本节主要介绍频率(周期)的估计与频谱图的生成。 1、频率、周期的估计
对于Y(k?f)，如果当k?f = f时，Y(k?f)取最大值，则f为频率的估计值，由于采样间隔的误差，f也存在误差，其误差最大为?f / 2。 周期T=1/f。
从原理上可以看出，如果在标准信号中混有噪声，用上述方法仍能够精确地估计出原标准信号的频率和周期，这个将在下一章做出验证 2、频谱图
9