《dsp技术3数字信号处理.ppt》由会员分享，可在线阅读，更多相关《dsp技术3数字信号处理.ppt（57页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、数字信号处理第一部分 离散时间信号与系统1.1 离散时间信号o将连续时间信号x(t)按照时间间隔T抽样，就形成了只有在离散时间点t=nT上才有非零值的信号，称为离散时间信号，其数学模型是将x(t)乘以单位冲激串T(t)：离散时间信号o如果将周期T归一化为1，xnT就成为离散时间序列xn（discrete time series）o离散时间信号尽管在时间上是离散的，但在幅度上仍然是连续变化的，因此仍然是模拟信号，只有经过量化器，也就是将各离散时间点上的信号的幅度归并到有限的若干个幅度电平上，并用数字来表示，才能称其为数字信号（digital signal）。1.2 典型的离散时间信号o单位数字冲
2、激信号n典型的离散时间信号o单位阶跃序列典型的离散时间信号o门序列典型的离散时间信号o正弦序列正弦序列o正弦序列本身就有一个频率f，再加上一个采样频率fs=1/T，就产生了一个f与fs的关系问题，令：o称为圆周频率圆周频率，它是一个相对频率量，也是一个归一化的频率，记：典型的离散时间信号o实际上，任何序列都可以表示成n的移位移位加权和加权和，这是一个非常重要的概念，因为我们在讨论信号通过系统时，只需要讨论单位数字冲激通过系统时的响应，即单位冲激响单位冲激响应应，然后将系统对各移位的冲激的响应叠加，就得到系统的输出。o整个一部信号与系统的理论，就是建立在对信号用不同形式分解，然后再综合的基础上的
3、。典型的离散时间信号o时域分析，就是将信号分解为n的移位加权和o频域分析，就是将信号分解为不同频率的正弦（sin）信号之和o复频域分析，就是将信号分解为不同的复指数之和1.3 离散线性位移不变系统o一个离散时间系统，可以抽象为一个变换或映射，把输入序列xn变换为输出序列yn:o我们来看两个例子：例1：yn=xn+ayn-1o这是一个一阶差分方程，而且是递归的，其框图如下：离散线性位移不变系统o一个离散的时间系统总是由延时器、加法器与数乘器组成的：例2:yn=b0 xn+b1xn-1+b2xn-2o这是一个三点的加权平均器离散线性位移不变系统o如果系统的输入xn是单位冲激n，则输出yn就是系统的
4、单位冲激响应hn。在例2中，若xn=n，则：hn=b0n+b1n+b2no这是一个有限冲激响应系统（FIR：Finite Impulse Response），而例1由于有递归，因而是一个无限冲激响应系统（IIR：Infinite Impulse Response）。系统的线性性o设 x1n y1n x2n y2no如果满足ax1n+bx2n ay1n+by2no则称系统是线性的系统的位移不变性o设 xn yno如果满足 xn-n0 yn-n0o则称系统是位移不变的线性位移不变系统o同时具有线性性与位移不变系统的离散时间系统，称为线性位移不变系统。o只有线性位移不变系统，才可以将表示为单位数字冲
5、激的移位加权和的输入信号中的每一个数字冲激通过系统，然后将其响应相加，得到总的输出，这是讨论问题的基础。1.3 离散信号通过线性位移不变系统o将输入信号xn表示为n的移位加权和：o将其中的每一个冲激xkn-k通过系统，得到相应的响应xkhn-k，然后将所有的输出叠加，得到总的输出：o称为xn与hn的卷和卷和，其性质类似于连续系统中的卷卷积积。卷和示例卷和示例 xn=hn=n+n-1+n-2o求oyn=n+2n-1+3n-2+2n-3+n-4数字信号处理第二部分 Z变换及离散时间系统分析2.1 Z变换的定义o给定一个离散信号xn，n=(-,)，可以直接给出xn的Z变换：oxn的Z变换也可以由拉普
6、拉斯变换推导出来2.2 Z变换与拉普拉斯变换的关系o信号f(t)的拉普拉斯变换的定义o式中，s=+j，因此可以将上式写为：o由此可以看出拉普拉斯变换与傅里叶变换的关系Z变换与拉普拉斯变换的关系o从z=es=e+j=eej可以看出拉氏变换所在的S平面和Z变换所在的Z平面之间的映射，即直角坐标系和极坐标系之间的映射。oS平面的j轴（=0）映射为Z平面的单位圆（r=1）oS平面的左半平面（0）映射为Z平面的单位圆内（r0）映射为Z平面的单位圆外（r1）Z变换与拉普拉斯变换的关系o时域离散，则频域是周期的o在Z平面上，=t=2f/fs，fs是采样频率。o当频率从0fs时，从02 o或者说，当频率从-f
7、s/2变到fs/2时，从-o每当变化2/T时，Z平面上完成一个圆周，从而体现其周期性o当=0时，拉氏变换就演变为傅里叶变换，即在S平面上，傅里叶变换始终是在j轴上进行的o映射到Z平面上，若变量z仅在单位圆上取值，Z变换也就演变成为傅里叶变换2.3 常用信号的Z变换o单位数字冲激 no只有当k=0时，n=1，这时z-k=1o单位阶跃序列un常用信号的Z变换o指数序列e-anTo正弦序列sin(nT)o利用欧拉公式sinx=(1/2j)(ejx-e-jx)得：利用指数Z变换的结果2.4 Z变换的性质o将xn及其Z变换X(z)之间的关系化为xn X(z)o设a，b为常数，则oax1n+bx2n aX
8、1(z)+bX2(z)o这就是我们熟悉的线性性Z变换的性质o移位性oxk+NzNX(z)oxk-N z-NX(z)o由于时移是数字系统的基本运算之一（其余两种是数乘和求和），在时域作单位时延（即延时一个样本），相应的Z变换多一个z-1因子，因此，往往将z-1用来表示单位时延。Z变换的性质o卷和定理ox1n*x2n X1(z)X2(z)o这是一个极为重要的定理，类似于模拟系统的卷积定理，将时域内比较繁琐的卷和运算，转化为Z域内的算术乘法。2.5 离散信号通过线性位移不变系统o对于线性位移不变系统，满足以下关系：o其中，yn=xn*hn Y(z)=X(z)H(z)数字信号处理第三部分 离散时间信号
9、的傅里叶变换3.1 连续时间信号的傅里叶变换3.2 离散时间序列的傅里叶变换(DTFT)o定义o设线性位移不变系统的单位冲激响应为hn，则该系统的频率响应为：o这就是离散时间序列的傅里叶变换，是以2为周期的周期函数，也可以将其看成周期信号H(j)在频域内展开的傅里叶级数，其傅里叶系数即是hno由Z变换的定义，容易得到H(j)=H(z)|z=ejDTFT的性质o即离散时间序列的傅里叶变换就是仅在单位圆上取值的Z变换。o线性性线性性o时移性：时移性：设xn X(j)，则o时域卷积定理时域卷积定理 若yn=xn*hn 则Y(j)=X(j)H(j)o频域卷积定理频域卷积定理 o若yn=xn hn 则Y
10、(j)=X(j)*H(j)o注意：时域是离散的，频域是周期的，故上述积分只在一个周期内（-，）进行DTFT的性质o能量定理（Parseval定理）o信号在时域的总能量等于其在频域的总能量o这个结论是显然的，因为无论是时域的xn，还是频域的X(j)，它们所表征的是同一个信号，其能量当然应该是相同的。o由于离散时间信号的频域表示是周期的，则其频域的总能量是在一个周期内做积分，因此|X(j)|2是信号的能量谱，而|X(j)|2是信号在d这一极小的频带内的能量3.3 离散傅里叶变换（DFT）oDTFT是离散时间序列的傅里叶变换，其时域是离散的，谱是周期的，但还是连续的，仍然不能用计算机来处理，因为计算
11、机要求时域和频域都是离散的，实际上是要求时域和频域都即是离散的，也是周期的。o时域的连续的非周期信号频域连续非周期谱（FT）o时域的周期信号频域离散谱（FST）o时域的离散信号频域周期谱（DTFT）o时域的离散周期信号频域离散周期谱（DFT）o在DFT中，变换的两边都是离散的，从而才是真正能用计算机来做数字信号处理的变换对o两边都是周期的，从而处理可以只在一个周期内进行，这有两个重要的意义：一是所做的处理是有限的（这对计算机来说是必须的）；二是只在一个周期内就可以保留全部的信息（这对准确处理来说是必须的）。信号的离散周期化o工程中的实际信号往往是连续非周期的，为了进行DFT，就必须将其离散周期
12、化o离散化就是采样的过程离散化就是采样的过程o周期化分为两种情况周期化分为两种情况:n若xn是有限长序列，长度为N，则将该N点序列看成是周期信号的一个周期，进行周期化延拓;n若xn是无限长序列，首先将其截为长度为N点的有限长序列，然后再做周期化延拓。截尾的过程，必然带来失真，因此要研究各种截尾的方法，这也就是各种加窗算法的由来.信号的离散周期化DFT的定义o这是一个以N为周期的信号，k表示它是离散的，每经过N点，在单位圆上转一周其中特别注意旋转因子DFT的性质o线性性oDFTax1n+bx2n=aX1k+bX2ko移位特性oDFTxn+m=W-kmXkoDFTxn-m=WkmXko由于xn是周
13、期序列中的一个周期，所以对xn的移位应该是整个序列的移位，即前面的移出去，后面的移进来，移位后仍然是N点的周期序列，这就称为循环移位循环移位循环移位o非周期函数移位 周期函数循环移位DFT的性质o能量定理（Parseval定理）o时域一个周期内的信号的能量和频域内一个周期的能量是相等的，这个结论的正确性是显然的。DFT的性质o时域的循环卷和o式中的n mod N表示以N为模对n求余o做循环卷和时，将一个信号翻摺，是将整个序列翻摺，而不仅仅是将一个周期翻摺o在翻摺后移位时，在一个周期内，有移出就有移入，求和始终是在一个周期内进行、o循环卷和的结果yn也是周期的，周期仍然是N数字信号处理第四部分
14、模拟信号的数字化4.1 Nyquist采样定理o截止频率为fm的带限信号，由相距不大的于1/2fm的均匀间隔上的采样序列唯一的确定。oT=1/2fm称为Nyquist间隔（既不失真的最大采样间隔）ofs=q/T=2fm称为Nyquist频率（既不失真的最低采样频率）o由于理想低通滤波器是不可实现的，实际的低通滤波器都有上升时间和下降时间，因此工程上往往采用fs3fm不满足Nyquist条件会产生频谱混叠4.2 采样o首先让信号通过一个低通滤波器（称为防混叠滤波器），使其成为带限信号o往往将采样频率选择为防混叠滤波器带宽的3倍，以保证无混叠产生o采样脉冲尽可能的窄，以便接近理想的冲激序列，否则，
15、采样信号会产生Sa(x)的包络调制o使用采样/保持器（S/H），使采样得到的电平保持一段时间不变化，以便A/D变换器有充分的时间将采样信号量化4.3 量化o采样后的信号，时间上是离散的，但幅度仍然是连续的，还不是真正的数字信号，必须对其幅度作量化，即将离散时间点上的幅度归入有限数目的等级o这个过程称为A/D变换，目前都是用专门的A/D变换器来完成，不同的A/D变换器，工作方式不同，如逐次逼近，闪烁(flash)，积分，-调制等逐次逼近A/D变换器o数字输出从最高值开始，通过D/A与输入比较，决定该位是0还是1，以2n逐次降低，逐级比较，最后得到从最高位到最低位的数字输出on位转换至少需要n个时
16、钟周期，n越大，即分辨率越高，转换速度就越慢闪烁式A/D转换器o输入信号同时与2n-1个不同的参考电平进行比较o在同一个时钟周期内完成转换，转换速度最快，其精度取决于参考电平的精度与分级的精度4.3 量化噪声oA/D的量化分级越细，精度越高o8-bit 256级o12-bit 4096级o两个量化等级之间的步长q越小，精度越高，因为舍入误差越小，但无论如何还是会引入误差，这种因为量化而引入的误差，称为量化噪声（量化误差）量化噪声o设量化噪声在步长q上式均匀分布的，则因为量化而引入的均方误差为：o其均方根值为：量化噪声o设A/D的满度正向电压为viq=2vi/2n vi=2n-1qS/N=20
17、log10vi/erms =(6.02n+4.77)dBo即每增加1bit，S/N提高约6dBo12bit的A/D S/N70dBo16bit的A/D S/N94dB量化噪声o量化噪声的降低（即S/N的提高），是要付出代价的n转换速度n成本n存储空间与数据吞吐速度o当信号本身的噪声已经大于量化噪声时，再降低量化噪声已经没有任何意义了量化噪声o上述信噪比的推导，是以信号满度为前提的o输入信号大于满度信号vi时，饱和，也就产生大的误差o输入信号太小，量化噪声加大o因此，往往在采样前加入一级低噪声放大器，程控其增益，使输入信号即接近满度，又不饱和，从而达到最好的信噪比4.4 数据采集的其它问题o程控防混叠滤波o程控前置放大器（衰减器）o与传感器的耦合方式o与外部电路的隔离o共模抑制o多通道时，是分时采集，还是同时采集