资源预览内容
第1页 / 共44页
第2页 / 共44页
第3页 / 共44页
第4页 / 共44页
第5页 / 共44页
第6页 / 共44页
第7页 / 共44页
第8页 / 共44页
第9页 / 共44页
第10页 / 共44页
亲,该文档总共44页,到这儿已超出免费预览范围,如果喜欢就下载吧!
资源描述
第六章 数字信号处理,概述数字信号处理的基本步骤信号数字化出现的问题离散傅里叶变换快速傅里叶变换(FFT)其他谱分析技术FFT算法的应用思考题(作业),概述,什么是信号处理、信号分析信号处理和信号分析没有明确的界限。通常把研究信号的构成和特征值称为信号分析,把信号经过必要的变换以获得所需信息的过程称为信号处理信号处理可用模拟信号处理系统和数字信号处理系统来实现模拟信号处理系统由一系列能实现模拟运算的电路,诸如模拟滤波器、乘法器、微分放大器等环节组成。模拟信号处理也作为任何数字信号处理的前奏,例如滤波、限幅、隔直、解调等预处理。数字处理之后也常需作模拟显示、记录。数字信号处理是用数字方法处理信号,它既可在通用计算机上借助程序来实现,也可以用专用信号处理机来完成。数字信号处理具有稳定、灵活、快速、高效、应用范围广、设备体积小、重量轻等优点,在各行业中得到了广泛的应用。,概述,数字信号处理应用领域生物医学工程、声学、声纳、雷达、宇航、地震、语音通讯等领域在医学领域,心电、脑电分析、超声、CT(电子计算机体层扫描仪)、核磁共振技术等,都是近年来医学上的重大科技成果。可以说,它们的基本原理都是传感与计算机技术的有机结合,是依据数字信号处理的数据或图形,获取人体内器官的工作状态或体内异物是否存在、形状、大小、位置等信息。在军事上,利用声纳、雷达、卫星系统,实现军事目标的探测,这同样是利用数字信号的存贮或处理,以获取有关军事设施、兵力部署,或者飞机、舰艇的航行位置、速度、方向等信息。,概述,机械、冶金、建筑、交通、电力等部门,数字信号处理的应用亦得到了迅速的发展。例如,大型旋转机械的监测与故障诊断系统,是将振动、声音、位移、温度等物理量,通过传感器转换为电信号,输入计算机系统,对信号处理加工变换,得到一系列特征参数,以监测系统的工作状态即便是普通的家用电器,如冰箱、电风扇、空调器、抽油烟机等,运用数字信号处理技术,找出振动与噪声源,提出对策,这对于提高产品质量亦是一种科学手段。数字信号处理技术的发展除了在通用计算机上发展各种数字信号处理软件外,还发展了有专用硬件的数字信号处理机。在运算速度、分辨能力、功能等方面,都显示出优越性。采用数字信号分析技术,对于1024采样点进行A/D转换,仅需4-15s,FFT运算仅需250ms,较快的只需数ms。因此,数字信号处理具有很强的实时能力,这使旋转机械的动态谱分析,以及高速运动物体的实时分析与监视均成为可能。此外,数字信号处理的分辨能力在高频段(50kHz),可达25Hz;在超低频段可达0.0025Hz。,概述,数字信号处理的主要内容频谱分析包含有相关与统计分析,其数学运算的核心是离散傅里叶变换与快速傅里叶变换。数字滤波包含了无限冲激响应滤波(IIR)与有限冲激响应滤波(FIR)等。,6-1 数字信号处理的基本步骤,一、信号预处理把信号变成适于数字处理的形式,以减轻数字处理的困难。电压幅值调理,以便适宜于采样。希望电压幅值峰峰值足够大,以便充分利用A/D转换器的精确度。如12位的A/D转换器,其参考电压为5V。由于2124096,故其末位数字的当量电压为2.5mv。若信号电平较低,转换后二进制数的高位都为0,仅在低位有值,其转换后的信噪比将很差。若信号电平绝对值超过5V,则转换中又将发生溢出,这是不允许的。必要的滤波,以提高信噪比,并滤去信号中的高频噪声。隔离信号中的直流分量(如果所测信号中不应有直流分量)。 如原信号经过调制,则应先行解调。,6-1 数字信号处理的基本步骤,二、A/D转换模拟信号经采样、量化并转化为二进制数的过程。,6-1 数字信号处理的基本步骤,三、运算处理数字信号处理器或计算机对离散的时间序列进行运算处理。计算机只能处理有限长度的数据,所以首先要把长时间的序列截断,对截取的数字序列有时还要人为地进行加权(乘以窗函数)以成为新的有限长的序列。对数据中的奇异点(由于强干扰或信号丢失所引起的数据突变)应予以剔除。对温漂、时漂等系统性干扰所引起的趋势项(周期大于记录长度的频率成分)也应予以分离。如有必要,还可以设计专门的程序来进行数字滤波。然后把数据按给定的程序进行运算,完成各种分析。,6-2 信号数字化出现的问题,一、概述信号数字化过程包括:采样:模拟信号离散信号量化:离散信号数字信号(时间序列)截断(加窗):无限长时间序列有限长时间序列DFT(离散傅里叶变换):连续频谱离散频谱每一步骤都可能引起信号和其蕴含信息的失真。,6-2 信号数字化出现的问题,1、采样:用一个等时距的周期脉冲序列s(t)(采样函数)乘以原模拟信号x(t)。时距Ts称为采样周期,fs=1/Ts称为采样频率。,信号x(t)及其频谱|X( f )|,采样信号s(t)及其频谱|S(f)|,采样后的信号x(t)s(t)及其频谱|X(f)*S( f )|,Ts太大即fs 太小时,平移后的图形会交迭!,6-2 信号数字化出现的问题,2、截取(加窗处理)计算机只能对有限长的序列进行运算,对长时间序列进行截短,相当于对采样后的信号进行加窗处理(加矩形窗),设窗宽为T,则t T 时,视x(t)=0。窗内数据点数(序列长度)N=T/Ts 。,进入计算机的信号:x(t)s(t)w(t),其频谱为X( f )*S( f )*W( f ),是频域连续函数。因W( f )的旁瓣引起了新频谱的皱波!,时窗函数及其幅频谱,有限长离散信号及其幅频谱,6-2 信号数字化出现的问题,3、频谱的离散化x(t)s(t)w(t)的频谱是一连续的频率函数,但计算机只能处理数字信号,用离散傅里叶变换(DFT)对离散的x(t)s(t)w(t)进行傅里叶变换后,输出的是离散的频率序列,即对X( f )*S( f )*W( f )进行了频域的采样处理。相当于在频域中乘上了频域采样函数D( f )。,计算机输出的频率序列X( f )p对应的时域函数x(t)p既不是原来的时域函数x(t),也不是x(t)s(t),而是一个周期函数。与原信号有很大差别!,频域采样函数及其时域函数,DFT后的频谱及其时域函数(t)p,6-2 信号数字化出现的问题,小结从以上过程看到,原来希望获得模拟信号x(t)的频域函数X( f ),但由于输入计算机的数据却是序列长为N的离散采样后信号x(t)s(t)w(t),计算机输出的是X( f )p。X( f )p已非X( f ),而是用X( f )p来近似X( f )。处理过程中的每一个步骤:采样、截断、DFT计算都会引起失真或误差,必须充分注意。好在工程上不仅关心有无误差,而更重要的是了解误差的具体数值,以及是否能以经济、有效的手段提取足够精确的信息。只要概念清楚,处理得当,就可以利用计算机有效地处理测试信号,完成在模拟信号处理技术中难以完成的工作。,6-2 信号数字化出现的问题,二、时域采样、混叠和采样定理1、采样把连续时间信号变成离散时间序列的过程,即按一定时间间隔读取连续信号瞬时值的过程。设采样长度为T,采样间隔为Ts,则采样得到的离散时间序列可表示为x(n)=x(nTs)=x(n/fs),n=0,1,2,N-1,N=T/Ts。采样频率fs=1/Ts。若Ts很小,fs很高,当T一定时,N很大,计算工作量很大;若N一定,而fs很高时,T很小,只能处理很短的时间历程。,6-2 信号数字化出现的问题,2、混叠现象图a)中,设1、2、3、4为四个连续的采样点,此时,信号A、B、C正好在这四个采样点的瞬时值是相等的,所以,用这些采样点的值无法分清信号A、B、C的区别。图b)是用过大的Ts对两个不同频率的正弦波采样的结果,得到一组相同的采样值,无法辨别两者的差别。采样后,将高频信号误认为某种相应的低频信号的现象称为混叠现象。造成这种现象的原因是采样频率太低而采样间隔太大。,6-2 信号数字化出现的问题,从频域来看,信号经时域采样之后成为离散信号,新信号的频域函数是一周期函数,周期为1/Ts=fs。若Ts太大, fs太低,平移距离1/Ts过小,那么移至各采样脉冲所在处的频谱X( f )就有一部分相互交叠,新合成的X( f )*S( f )图形与原X( f )不一致,这种现象称为混叠。发生混叠后,改变了原来频谱的部分幅值,这样,就不可能从离散的采样信号x(t)s(t)准确地恢复原来的时域信号x(t)。,6-2 信号数字化出现的问题,3、采样定理因为原信号的频谱X( f )是f的偶函数,并以f=0为对称轴,X( f )*S( f )又是以fs为周期的函数,故有混叠时,必出现在f=fs/2的两侧,称fs/2为折叠频率。为避免频率混叠,应使fs=1/Ts2fh(如图),fh为x(t)的最高频率分量。对不满足要求的信号,采样之前,应通过模拟低通滤波器滤去高频成分,为满足上述条件创造条件。这种处理称为抗混叠预处理。把一个没有混叠的频谱通过中心频率为0,带宽为(fs/2)的理想低通滤波器,就可完整地取出原信号的频谱,从而有可能从此频谱中准确地恢复原模拟信号x(t)。为避免混叠以使采样处理后仍有可能准确地恢复原信号,采样频率fs必须大于信号最高频率fh的两倍,即fs2fh,这就是采样定理。由于理想低通滤波器不存在,故采样频率常取为(34)fc, fc为滤波器的截止频率。,6-2 信号数字化出现的问题,三、量化和量化误差采样是对时间坐标的离散化,使连续的模拟信号变成了离散信号,但计算机只能处理数字信号,量化则是对幅值坐标的离散化,使离散信号变成数字信号。量化是从一组有限个离散电平中取一个来代表信号采样点的实际幅值电平。当A/D转换器的位数为b,动态范围为D(如5V)时,相邻两量化电平之间的差为: x=D/2(b-1) (字长第一位用作符号位)量化误差为(x/2,+x/2),当b一定时,量化误差是一定的。A/D转换的位数越多,转化精度越高,但转换速率下降,成本也显著增加。,6-2 信号数字化出现的问题,四、截断与泄漏计算机只能处理有限长度的时域信号,因此,对长时间历程信号必须截断,亦即必须对实际信号乘以时域上有限的矩形窗函数,得:,由于W( f )是一个无限宽的sinc函数,所以,既使x(t)是带限信号,在截断后也必然成为无限带宽的信号,这种信号的能量在频率轴上分布扩展的现象称为泄漏。同时,由于截断后信号带宽变为无限宽,因此,无论采样频率多高,信号总是不可避免地出现混叠而导致一些误差。为减小截断的影响,常采用其他的时窗函数来对所截取的时域信号进行加权处理,6-2 信号数字化出现的问题,五、时域周期延拓和栅栏效应x(t)经过时域采样、截断、频域采样之后,所得到的信号x(t)s(t)w(t)*d(t)是一个周期信号。原信号x(t)不一定是周期信号,但经频率离散化之后,将其改造成了周期信号,这一现象称为时域周期延拓。对一函数实行采样,其效果有如透过栅栏的缝隙看外景一样,只有落在缝隙前的少数景象被看到,其余景象都被栅栏档住,视为0,这种现象被称为栅栏效应。时域采样和频域采样均有栅栏效应。对时域采样,若满足采样定理要求,栅栏效应不会有什么影响。而对频域采样,“档住”的频率成分可能是重要的或具有特征的成分,使整个处理失去意义,6-2 信号数字化出现的问题,六、频率分辨力、整周期截取频率采样间隔f 是频率分辨力的一项指标, f 越小,被档住的频率成分就越少,频率分辨力就越高。 由 f =fs/N=1/T (fs采样频率;T时窗宽度) 可知,频率分辨力与计算工作量是矛盾的。fs由采样定理决定, fs决定后,为减小f 而必须增大N,从而急剧增加计算工作量。解决这一矛盾的途径有两条:一是在DFT的基础上,采用频率细化技术(ZOOM),提高感兴趣的局部频段的分辨力;二是采用DFT之外的其他方法得到频谱。整周期截取是指对周期信号用数字化处理时,截取信号的长度T必须为信号周期的整数倍。只有这样,才能获得精确的频谱,以及使周期延拓后的信号和原信号完全重合,否则,波形和频谱都会发生畸变,
网站客服QQ:2055934822
金锄头文库版权所有
经营许可证:蜀ICP备13022795号 | 川公网安备 51140202000112号