现代信号处理第2章

2015-11-2机械工程及自动化研究所现代信号处理技术及应用第二章信号的时域分析西安交通大学机械工程学院研究生学位课程第二章信号的时域分析2.1信号的预处理2.2信号的采样2.3时域统计分析2.4相关分析及应用引言以时间为自变量描述物理量的变化是信号最基本、最直观的表达形式。在时域内对信号进行滤波、放大、统计特征计算、相关性分析等处理，统称为信号的时域分析。通过时域分析方法，可以有效提高信噪比，求取信号波形在不同时刻的相似性和关联性，获得反映机械设备运行状态的特征参数，为机械系统动态分析和故障诊断提供有效的信息。第二章信号的时域分析2.1信号的预处理2.2信号的采样2.3时域统计分析2.4相关分析及应用2.1信号的预处理传感器获取的信号往往比较微弱，并伴随着各种噪声。不同类型的传感器，其输出信号的形式也不尽相同。为了抑制信号中的噪声，提高检测信号的信噪比，便于信息提取，须对传感器检测到的信号进行预处理。所谓信号预处理，是指在对信号进行变换、提取、识别或评估之前，对检测信号进行的转换、滤波、放大等处理。2.1信号的预处理常用的信号预处理方法信号类型转换应变测力传感器、热电阻传感器输出的信号均为电阻信号，为了便于后续处理常用电桥将电阻信号转变为电压信号信号放大常用的信号放大器包括：测量放大器、隔离放大器、可编程增益放大器等信号滤波（本节重点介绍）去除均值在计算信号的标准差等统计量时，需要去除信号均值去除趋势项常用的趋势项消除方法有滤波法、多项式拟合法2.1.1信号的滤波处理信号滤波处理是消除或减弱干扰噪声，保留有用信号的过程。把实现滤波功能的系统称之为滤波器。滤波器可分为两大类，即经典滤波器和现代滤波器。1.经典滤波器定义：当噪声和有用信号处于不同的频带时，噪声通过滤波器将被衰减或消除，而有用信号得以保留分类根据幅频特性的不同，滤波器分为低通滤波器、高通滤波器、带通滤波器、带阻滤波器等类型。根据处理信号类型的不同，滤波器可分为模拟滤波器和数字滤波器。对于数字滤波器来说，根据滤波器的单位脉冲响应序列长度的无限和有限，数字滤波器可进一步分为无限冲击响应滤波器(IIR)和有限冲击响应滤波器(FIR)两类1)经典滤波器原理经典滤波概念和方法建立在频域分析基础上)()()(tntstx滤波后的信号为)(ty滤波器的传递函数或滤波器的频率响应函数)(/)()(XYH在噪声频带和有用成分频带分离的情况下，通过设计如下的滤波器函数时当时当0)(,0)(0)(,1)(SHSH)(*)()(thtxty)}({)(1HFth滤波器的单位脉冲响应函数()H)()()(NSX)()(SY(2.1.1)2)理想模拟滤波器理想模拟滤波器是一个理想化的模型，对其讨论有助于进一步了解和改进实际滤波器的性能理想模拟滤波器的幅频特性曲线|H(j)|(a)低通c|H(j)|(b)高通c|H(j)|(c)带通c1c2|H(j)|(d)带阻c1c2图2.1.1理想模拟滤波器由于理想低通滤波器具有矩形幅频特性和线性相位特性。同时，理想高通、带通和带阻滤波器均可以由理想低通滤波器串联得到。因此，以后均以理想低通滤波器为例来说明理想低通滤波器的矩形幅频、相频特性ccjjH)(1)(0|H(j)|-cc01幅频特性texitezin(j)0-0相频特性(2.1.4)理想模拟滤波器理想低通滤波器的单位脉冲响应函数)]([csin)()(sin)(000tttthccc00CCth(t)h(t)=2=3=4tccc(2.1.5)3)实际滤波器及其基本参数实际的滤波器为了物理上可实现，通常在通带和阻带之间设置过渡带H(j)01pss1-pc0.707cps通带边缘频率阻带边缘频率截止频率(滤波器的半功率点)实际滤波器及其基本参数实际滤波器的幅频特性幅值在通带和阻带内一般不严格为1和0。它们分别允许的波动量分别为ps波动的大小分别用通带和阻带内的衰减110|()|20lg20lg|()|20lg(1)|()|ppjjppjHeHeHe110|()|20lg20lg|()|20lg(1)|()|ssjjssjHeHeHe实际滤波器的参数还有：波纹幅度、带宽、品质因数和倍频程选择性等4)数字滤波器的设计数字滤波器有无限冲击响应IIR型滤波器和有限冲击响应FIR型滤波器之分IIR型数字滤波器的传递函数是NkkkMrrrzazbzH101)(jtzeFIR型数字滤波器的传递函数是10)()(NnnznhzH)(nh滤波器的单位脉冲响应函数数字滤波器的设计方法简述FIR型数字滤波器的设计方法主要是建立在对理想滤波器频率特性作某种近似的基础上的。这些近似方法有窗函数法、频率抽样法等；IIR型数字滤波器的设计属于间接设计法。IIR型数字滤波器目前最通用的设计方法是利用已经很成熟的模拟滤波器的设计方法来进行设计。而模拟滤波器的设计方法又有巴特沃斯(Butterworth)滤波器、切比雪夫(Chebyshev)和椭圆滤波器等不同的设计方法。数字滤波器的设计方法简述不论是FIR型数字滤波器还是IIR型数字滤波器的设计都包括三个步骤Matlab数字滤波器设计演示(FDATool.fda,SPTool.fda)2．现代滤波器当噪声频带和有用信号频带相互重叠时，经典滤波器就无法实现滤波功能现代滤波器也称统计滤波器，从统计的概念出发对信号在时域进行估计，在统计指标最优的意义下，用估计值去逼近有用信号，相应的噪声也在统计最优的意义下得以减弱或消除常用的统计滤波器有维纳滤波器和卡尔曼滤波器两类1）维纳滤波器20世纪40年代第二次世界大战期间，由于军事上的需要，Wiener提出并解决了平稳过程的最佳线性滤波问题采用线性最小均方误差估计准则，设计的最佳滤波器称为维纳(Wiener)滤波器)()()(tntstxdthxtsftt0)()()(ˆ22ˆ{()}min[()()]EetEststh(t):脉冲响应函数，按最小均方误差准则确定(2.1.8)(2.1.9)维纳滤波器维纳滤波器可根据t时刻，及t以前时刻的观测值x(t)，实现以下三个方面的应用2）卡尔曼滤波器维纳滤波器由于计算量大，难以作实时处理，故不能广泛应用，同时它对非平稳信号的滤波也无能为力。60年代初由于航天事业发展的需要，卡尔曼(Kalman)和布西(Bucy)在解决非平稳、多输入输出随机序列的估计问题中引入了状态变量，在克服维纳滤波某些局限的基础上，提出了被后人称为卡尔曼滤波的新滤波方法。该方法在雷达、通信、控制、生物和勘探等领域得到了广泛的应用卡尔曼滤波器卡尔曼滤波是线性最小均方误差滤波器的另一种处理方法。卡尔曼滤波建立在已知随机信号模型的基础上，包括模型的阶次、模型的参数和激励白噪声的统计特性。其原理是用信号前一时刻的估计值与测量值的误差项的加权平均作为当前时刻的估计值。ˆˆˆ()(1)[()(1)]nsnasnGxnacsn(2.1.10)信号前一时刻估计值测量值的误差项常数a和c分别表示参数模型和测量模型的参数，用Gn表示时刻的加权系数。维纳滤波器与卡尔曼滤波器对比卡尔曼滤波的特点是把信号的先验知识用信号的模型表达出来；在时域上引入状态变量法进行处理；采用递推型的线性最小均方误差算法。卡尔曼滤波和维纳滤波都是在应用随机信号和观测噪声的前二阶矩的统计特性，以线性最小均方估计解决随机信号的滤波问题。维纳滤波需要给出随机信号和噪声的有理谱形式；卡尔曼滤波则要求把随机信号规定为白噪声驱动的线性系统的输出。维纳滤波理论适应于平稳随机过程；卡尔曼滤波适用于有限初始时间的非平稳随机过程第二章信号的时域分析2.1信号的预处理2.2信号的采样2.3时域统计分析2.4相关分析及应用将连续信号转换成离散的数字序列过程就是信号的采样，它包含了离散和量化两个主要步骤。本节主要介绍采样过程中采样与混频、量化与误差、采样长度与分辨率及窗函数与泄露四方面的内容。2.2.1采样与混频设模拟信号为x(t)，间距为t的采样脉冲函数为p(t)()()+n=-pt=δt-nt采样过程是x(t)和p(t)相乘，得到离散信号x(nt)+n=-tt-nδtnx=tnx)()()(设x(t)的傅里叶变换为x()。采样脉冲函数p(t)的傅里叶变换为+m=-tmδt=P)2(2)((2.2.1)(2.2.2)模拟信号、采样脉冲函数及其频谱0tx(t)(a)原函数0X()(b)原函数频谱max-maxc0tp(t)1(c)采样冲击函数p(t)t0P()(d)采样冲击函数的频谱s=2/t离散信号及其频谱根据频域卷积定理可知，则式(2.2.1)所示离散信号x(nt)的傅里叶变换为离散信号及其频谱+m=-tm-Xt=X)2(2)(0nx(n)=x(t)p(t)(e)离散时间信号t0X()*P()(f)采样序列的频谱s=2/tmax-max频谱混叠，不能由频谱准确地恢复原信号(2.2.3)采样定理采样定理：为避免混叠，采样频率s必须不小于信号中最高频率max的两倍max2s即或max/tmax1/(2)tf实际中采样频率的选取往往留有余地，一般选取采样频率s为处理信号中最高频率的2.5～4倍。另外，由于测量信号中的高频部分往往是由干扰引起的噪声或我们不感兴趣的频谱，因此采样前须先对信号进行低通滤波(又称抗混滤波)。然后再根据滤波后信号的最高频率max设定采样间隔。2.2.2量化与误差量化是对信号采样点取值进行数字化转换的过程。量化结果以一定位数的数字近似表示信号在采样点的取值。由于模/数转换器的位数有限，模/数转换器的输出只能表达一系列具有一定间隔的电平。当模拟信号在采样点上的取值落在两个相邻电平之间时，就要舍入到相近的一个电平上，我们把这一过程称之为量化。若设模/数转换器的位数为N，采用二进制编码，转换器转换的电压范围为V，则相邻电平之间的增量为21NV(2.2.4)2.2.3窗函数或泄漏理论上信号的长度是无限的，但任何观测信号都是在有限时间段内进行观测的。因此，信号采样过程须使用窗函数，将无限长信号截断成为有限长度的信号。从理论上看，截断过程就是在时域将无限长信号乘以有限时间宽度的窗函数窗函数或泄漏最简单的窗函数是矩形窗T|t||t|=T.T|t|=tw0501)(TTT=W)sin(2)(w(t)tT01(a)矩形窗函数-TW()0(b)矩形窗函数幅频曲线Tπ-2TT无限带宽泄漏窗函数或泄漏泄漏与截断长度、所使用的窗函数等有关。不进行信号截断就没有泄漏误差另外，使用不同的窗函数泄漏大小也不同。泄漏取决于窗函数频谱的旁瓣。如果窗函数的旁瓣小，相应的泄漏也小其它窗函数三角窗Tt||T|t||t|T=tw011)(22/)2/sin()(TT=TWw(t)tT01(a)三角窗函数-TW()0(b)三角窗函数幅频曲线T2-T2T其它窗函数汉宁窗Tt|T|t|Tt=tw|0)cos(2121)(2)(11)sin()(TT=Ww(t)tT0(c)汉宁窗函数-TW()0(d)汉宁窗函数幅频曲线T1-T1T几种窗函数频谱的比较几种窗函数频谱的比较名称特点应用矩形窗Rectangle矩形窗使用最多，习惯上不加窗就是使信号通过了矩形窗。这种窗的优点是主瓣比较集中，缺点是旁瓣较高，并有负旁瓣，导致变换中带进了高频干扰和泄漏，甚至出现负谱现象。频率识别精度最高，幅值识别精