现代信号处理汇编资料

1、 现代信号处理 的几个(j )边沿问题 湖 南 工 学 院 曹 才 开共六十八页目录(ml)绪论 1 信号处理学科的地位 2 经典信号处理技术的困境 3 现代信号处理的基本内容现代信号处理的几个边沿问题介绍 1 时谱（倒谱）和功率频谱分析(fnx) 2 基于多尺度估计理论的信号检测 3 基于相关函数处理随机信号 4 实时检测的速度问题 三 结束语共六十八页一 绪 论共六十八页 1948创立的系统论、信息论和控制论三大科学理论，对于信号处理学科的发展起到非常重要的奠基和推动作用。系统论是美国生物学家贝格朗菲创立的，他为确立适用于系统的一般原则做出了重要贡献。信息论是美国数学家仙农建立的，它是现代

2、通信理论的基础，在计算技术、自动控制等方面得到广泛应用。控制论是美国数学家维纳提出的，它促进了通信、计算机和人工智能等方面得到广泛应用。随着大规模集成工艺和计算机技术的飞速发展，近几十年来，信号与系统学科得到惊人的发展。信号处理是信息论的一个分支学科，它的基本概念与分析方法还在不断的发展，其应用范围也在不断的扩大，它在通信、航空与航天、电工及电子电路、机械、声学、地震学、探矿(tn kung)、生物工程、能源、化学等许多领域里起着重要的作用。该学科水平的高低反映一个国家的整体科技水平。1 信号处理学科(xuk)的地位共六十八页 2 经典信号处理技术(jsh)的困境 二十世界60年代以来，由于微

3、电子集成电路技术(jsh)的不断发展，为复杂信号处理的实现提供了可能，极大促进了信号处理向新的领域发展。 随着科学技术的飞速发展，经典信号处理技术越来越力不从心。 (一) 局限性 (1) 假设信号及其背景噪声是高斯的和平稳的; (2) 其对象系统只限于时不变(或缓慢) 、线性、因果、最小相位的系统;共六十八页 (3) 信号分析方法只限于二阶矩特性和傅氏频谱。 (二) 傅里叶变换的困境 在信号分析和故障诊断技术等领域中，以前最为普遍是利用快速(kui s)傅里叶变换 (FFT) 的频域分析法，这种方法虽然能够分辨平稳信号在频域中的位置与大小，但对非线性、非平稳随机信号的检测问题、时域频域变化规律

4、等方面的分析显得力不从心。 这是因为傅里叶变换 把任意信号分解为 共六十八页 由于n=0,1,2,3,., 傅里叶变换算法的基函数是sin (t) 或cos (t) 的频率特性(点通) ，仅能检测平稳信号的整次谐波，不能检测信号的非整次谐波，所以，傅里叶变换算法不能检测非平稳信号的特性信息。 现代信号处理技术在非线性、非平稳随机信号的检测问题、时域频域变化规律等方面(fngmin)优于傅里叶变换的地方，是它在时域和频域同时具有良好的局部化性质，而且由于对高频成分采用逐渐精细的时域和空域取样步长，从而可以聚焦到对象的任意细节。由于这一特点，它能将不同频率组成的混合信号分解成不同频率的块信号，可有

5、效地进行信噪分离、信号特征提取、故障诊断等。 共六十八页 3 现代(xindi)信号处理的基本内容 (1) 统计信号处理 (2) 多维多信道(xn do)信号处理 (3) 非高斯信号处理 (4) 非平稳信号处理共六十八页 (1) 统计(tngj)信号处理 * 参数估计理论(lln) * 波形估计 * 现代频谱分析 * 自适应滤波 * 鲁棒参数估计 * 倒谱分析 * 统计性能分析 * 信号检测共六十八页 (2) 多维多信道(xn do)信号处理 * 二维信号处理的特点和难点 稳定性、谱因子分解、模型参数(cnsh)的非唯一可识别性、二维最大熵法尚无闭式解、多元时间序列分析等。 *多信道信号处理

6、多元时间序列分析(最佳线性预测、多元AR过程的建模、多元ARMA过程的建模等) 共六十八页 (3) 非高斯(o s)信号处理 非参数化双谱估计、非最小相位系统辨识、非因果系统辨识、有色噪声中的谐波恢复、非高斯(o s)噪声中非高斯(o s)信号检测等。 共六十八页 (4) 非平稳(pngwn)信号处理 基于时频分析的信号检测、基于多尺度估计理论的信号检测(小波变换、短时分形变换、分布式系统状态融合估计等)、智能信息处理技术(模糊(m hu)计算技术、人工神经网络)等。 一般来说,智能信息处理可以划分为两大类,一类为基于传统计算机的智能信息处理,另一类为基于神经计算的智能信息处理。 共六十八页

7、二 现代信号处理 的几个边沿问题(wnt)介绍共六十八页 1 时谱（倒谱）和功率(gngl)频谱分析 时谱分析(Cepstrum analysis)是一种(y zhn)非线性信号处理技术，它在语言、图像、和噪声处理领域中都有广泛的应用。时谱可分为两类：复时谱和功率时谱。 共六十八页 (1) 复时谱（Complex cepstrum）分析(fnx) 复时谱（Complex cepstrum）的定义(dngy)为： 共六十八页 由上式可见，复时谱实际上是序列x(n)的Fourier变换的自然对数，再取Fourier逆变换，得到的复时谱仍然是一个序列。也就是说，复时谱是x(n)从时间(shjin)域

8、至频率域、频率域至频率域、频率域至时间(shjin)域的三次变换。共六十八页例1: 设原信号是一个45Hz的正弦波，在传播过程中遇到障碍产生回声，回声振幅衰减为原信号的0.5，并与原信号有0.2s的延迟(ynch)。在某测点测到的信号是原信号和回声信号的叠加。试使用复时谱分析该测点的信号。由于复时谱从复频谱计算得到，不损失相位信息，因此复时谱是可逆的。实时谱过程是不可逆的。 MATLAB 仿真见图1 。 共六十八页 图1 正弦波与回声信号(xnho)叠加的波形和时谱形状共六十八页 (2) 功率(gngl)频谱(不是功率时谱)例2: 求信号 的功率谱。其中f1=50Hz,f2=120Hz,w(t

9、)为白噪声，采样频率(pnl)Fs=1000Hz。（1）信号长度 N=256；（2）信号长度 N=1024。共六十八页 图2 含有噪声信号(xnho)的功率谱(下图)和无噪声信号(xnho)的功率谱(上图)共六十八页 2 基于多尺度估计理论的信号(xnho)检测多尺度系统理论(Multiscale System Theory,MST)也称为(chn wi)多尺度估计理论,或称为(chn wi)多尺度变换。1990.12由法国A.S.Willsky数学家首先提出的。尺度变换:若原信号f(t),则f(at)称为对f(t)的尺度变换,根据f(t)不同,a取值方式不同,出现了许多种尺度变换。多尺度变换

10、:f(a 、b 、c 、)t,不同尺度(或分辨级)是根据信号几何图形不同而设计的,达到提高信号分辨等级之目的。共六十八页 * 两种基本多尺度变换: 小波分析(fnx),短时分形维数算法 (1) 小波分析 由于小波分析具有能够根据分析对象自动调整有关参数的“自适应性”和能够根据观测对象自动“调焦”的特性而广泛应用于各个领域。 Fourier变换不能提取信号中的奇异性和突变点的信息，它只是将这些(zhxi)信息铺开到整个频率轴上。但是小波变换Tg( )是将分解成具有局部特性的小波 : 共六十八页其中小波 是将具有局部特性的小波函数g(t)通过平移(pn y)和尺度变换（放大倍数为1/a）而构成的。

11、参数a具有时间的量纲，也称为小波尺度;f(t)为被处理的信号。 小波函数g(t)称为小波母函数,有多种,以便适应各种非平稳信号的检测。当对信号进行小波变换时，其局部化特性与所选取小波函数有关，因此，要根据信号的特征选择适当的小波母函数才能获得满意的检测效果。 共六十八页 常见(chn jin)的小波函数有: Harr小波函数 墨西哥小波函数 Daubechies小波函数 Morlet小波函数 Meyer小波函数 共六十八页Harr小波函数定义(dngy)为 :共六十八页 为什么叫“小波”? “小”,是指g(t)具有衰减性; “波”,是指g(t)具有波动性,即其振幅(zhnf)在正负相间进行震荡

12、。共六十八页应用举例: 基于小波包变换实时检测电机振动速度信号 小波包变换(WPD)不仅能检测非平稳信号的整次谐波，还能检测信号的非整次谐波，又因为小波变换本身对信号的奇异(qy)点十分敏感，这个特点可以用来跟踪电机振动速度信号。在虚拟仪器(VI：Virtual Instrument) LabVIEW 6.i平台上，基于小波包变换算法设计了VI程序，实现了电机振动速度信号实时检测系统。经过信号处理，该系统还具有信噪分离、测量电机振动功率谱、电机振动信号的时域频域变化规律、电机振动速度信号三维图、伴有噪声的原始振动波形和噪声波等测量功能。 共六十八页电机(dinj)振动自动检测系统的基本原理见图

13、1。 共六十八页 信噪 分离(fnl) 共六十八页 图3 基于小波分析(fnx)测得的电机振动速度信号 共六十八页 图4 基于小波分析测得的电机(dinj)振动功率谱 共六十八页 HP3562A动态信号(xnho)分析仪测得的结果 共六十八页图6为时域频域谱图；图7为电机某点振动(zhndng)速度信号的三维视图。 共六十八页(2) 短时分形(fn xn)维数算法 随着混沌分形理论研究的逐渐(zhjin)深入，应用范围也日益广泛，现在已渗透到图像处理、语音处理、模式识别、人工智能以及非线性电路等信息科学的许多分支。分形理论的研究对象是非线性系统中产生的不光滑和不可微的复杂信号波形，对应的定量参

14、数是分形维数。分形维数描述了系统及其测量时间序列的复杂程度。 共六十八页 分形维数作为图像模式的形态特征已用于图像分析与模式分类、图像生成、内插与计算机仿真、信号滤波、图像压缩编码、分形神经网络乃至于非线性混沌的研究。在工程实践中，人们对于复杂系统的探测往往通过(tnggu)对某时间变量的观测（即时间序列）来实现。但是测量噪声对系统重构及估计分形维数有着不良影响，因而必须对测量信号进行滤波。 共六十八页 分形(fn xn): 把信号波形看成一幅图,对图分割成网格(wn ),这就是分形。在n维空间上的网格是指的一个分割，Rn=U Xj，Xj 之间只相差一个平移或旋转，Xj称为一个格子。例如常用的

15、网格有正方形网格、矩形网格和三角形网格等。网格的尺寸定义为Xj的直径（或宽度），并记为网格，不同的网格格子之间是相似的。 事实上，网格分形是一种简单实用的分形方式，特别适用于计算机数字化处理。 共六十八页短时: 小时间(shjin)区间。 应用举例: 开关电源传导干扰信号的短时分形维数模糊控制滤波 基于短时分形维数的模糊控制滤波方法(fngf)，对开关电源传导干扰信号中的噪声进行滤波。该方法(fngf)提出了网络分形维数和短时分形维数的新算法，并讨论了模糊控制滤波方法(fngf)中的模糊控制参数的选取算法。基于虚拟仪器(VI) LabVIEW 6.i平台上对开关电源传导干扰信号进行实时检测。经

16、过信号处理，该系统还具有信噪分离、测量传导干扰功率谱等功能。结果表明，该方法滤波效果良好。 共六十八页 测量(cling)基本原理 线性阻抗稳定(wndng)网络（LISN） 图1 测量基本原理 共六十八页图2 附加(fji)噪声的 图3 短时分形维数模糊控制电流传导干扰信号 滤波后的电流传导干扰信号 共六十八页 图4 =0.05时的电流(dinli) 图5 =0.20时的电流 传导干扰信号 传导干扰信号 共六十八页模糊控制参数 是短时分形(fn xn)维数的函数 图4为 取固定值0.05测得波形，可见，取值太小，有用信号受损；图5为 取固定值0.20测得波形，可见，取值太大，噪声残留太多。由

17、此可见，固定值的滤波效果总是不太理想，而采用分形模糊(m hu)滤波方法能获得较好的滤波效果。因此，采用分形模糊(m hu)滤波的效果(图3 )优于取固定值时的滤波效果。 共六十八页 3 基于相关函数处理随机(su j)信号 随机信号和确定信号是两类性质完全不同的信号，对它们(t men)的描述、分析和处理方法也不相同。随机信号是一种不能用确定数学关系来描述的信号，无法预测未来时刻精确值的信号，也无法用实验的方法重复再现。 共六十八页 (1) 随机(su j)信号x(t)自相关函数 对于(duy)随机信号x(t)，自相关函数为：若去掉x(t)的均值部分，则相应的自相关函数称为自协方差，即： 共

18、六十八页对于离散随机信号(xnho)序列，x(n)的自相关函数和自协方差分别为： 式中，m为延迟。 共六十八页 (2) 应用举例(j l):地震信号估计 按照前面的计算方法，得到了中国(zhn u)数字地震台网（CDSN）改造后7个台站VHZ通道的功率谱密度估计（图6）。这里需要注意的是，由于功率谱密度不包含相位信息，所以不涉及地震仪器的相位响应。图6为中国数字地震台网记录的2001年昆仑山口西8.1级地震震后VH频段波形数据。 共六十八页图6 中国数字地震台网记录(jl)的2001年昆仑山口西8.1级地震震后VH频段功率谱密度波形共六十八页图7 中国数字地震(dzhn)台网记录的2001年昆

19、仑山口西8.1级地震后垂直分向波形 共六十八页 (2) 应用举例(j l): 设计数字滤波器对白噪声序列滤波 设计一个归一化频率为0.2的FIR数字滤波器，对一个白噪声信号序列进行滤波，对滤波后的信号绘制(huzh)置信区间为0.95的功率谱估计曲线。 共六十八页 图8 数字滤波器对白噪声(zoshng)序列滤波后的功率谱 共六十八页 4 实时检测的速度(sd)问题 在有些数字信号处理系统中，虽然运算简单，例如加法、乘法运算，但这些系统需要很快的实时处理速度(如实时处理图像(t xin)系统，要求采样率高达10MHz/s) ，这时，若用普通的数字处理系统、甚至是数字信号处理器(DSP)专用芯片

20、也很难达到要求，这时可以考虑采用特殊的硬件实现方法 ,例如，用简单的电路配合硬件查表法就可以解决这个问题。 共六十八页 (1) 查表法的基本(jbn)思想 事先把两个数相加(或相减)的结果按照输入数据的组合顺序制成一张表，并把它固化在ROM里。假设输入数据的字长为8位，在执行两个数相加(或相减)的运算时，把输入的两个8位数据组成一个16位地址，用这个地址作为(zuwi)ROM的输入地址，对存储在ROM中的数据表进行寻址，从ROM中就能很快读出相应的运算结果。表中的结果是事先计算好的，因此，这样可以大大加快运算速度。 共六十八页 (2) 基本硬件(yn jin)单元我们知道，在有些数字信号处理系

21、统中，基本硬件单元是加法器、乘法器和延迟器。 下面讨论(toln)用查表法实现这些基本硬件单元。 共六十八页1) 用查表法实现(shxin)加法器 EPROM27512 图1 用查表法实现(shxin)的实数加法器 共六十八页图1是用查表法实现的实数加法器。设A、B是两个8位输入数据，由于(yuy)8位数据所能表示的数的个数为 ，所以，A、B是两个8位数据的组合数为 个。因此可将计算好的加法表格数据存在存储器中，例如EPROM 27512型存储器( )。查表时,把A、B两个8位输入数据组合起来查找表格数据的地址。 共六十八页 两个(lin )复数相加 图2 用查表法实现(shxin)的复数加法

22、器共六十八页 可制作(A+C)及(B+D)两个表，分别用两个EPROM 27512型存储器进行(jnxng)实部和虚部相加运算，如图2所示。 共六十八页 2)用查表法实现(shxin)乘法器 图3 基于(jy)查表法实现乘法器 共六十八页 因为在二进制乘法运算(yn sun)中，都是通过反复相加和移位操作而实现的。用查表法实现乘法器原理与加法器类同，区别在于乘法运算(yn sun)表格中的数据为两数相乘的积。图3为两个实数的乘法器。 共六十八页 3) 单位(dnwi)延迟器 图4 用D触发器实现(shxin)单位延迟器 共六十八页 用D触发器实现单位延迟器如图4所示。D触发器型号为74LS273。当在触发器D端输入信号为 时，D触发器型Q端输出的信号为 ，即延迟了一个单位时间间隔(一个取样周期)。系统的取样速率、延迟器的移位和输出数据均由时钟(shzhng)CLK控制。 共六十八页 (3) 举 例 一阶IIR滤波器(IIR:无限脉冲响应数字滤波器, 即数字滤波器的单位冲激响应包含无限个采样值),对应的差分方程为 式中的乘法运算用查表选取EPROM27512实现；加法运算用查表选取 EPROM27512 实现，但考虑到DAC0832要求原码输入，为了避免前级码制转换的麻烦，所以(suy)也选用查表法。这里把实常数乘法器和加法器合并起来，用一个EPROM实现。