第一章 绪论
1、 试举例说明信号与信息这两个概念的区别与联系。
信息反映了一个物理系统的状态或特性，是自然界、人类社会和人类思维活动中普遍存在的物质和事物的属性。
信号是传载信息的物理量是信息的表现形式，如文字、语言、图像等。
如人们常用qq 聊天，即是用文字形式的信号将所要表达的信息传递给别人。
2、 什么是信号的正交分解？如何理解正交分解在机械故障诊断中的重要价值？
P9正交函数的定义
信号的正交分解如傅里叶变换、小波分解等，即将信号分解成多个独立的相互正交的信号的叠加。从而将信号独立的分解到不同空间中去，通常指滤波器频域内正交以便于故障分析和故障特征的提取。
傅里叶变换将信号分解成各个正交的傅里叶级数，将信号从时域转换到频域从而得到信号中的各个信号的频率。正交小波变换能够将任意信号（平稳或非平稳）分解到各自独立的频带中；正交性保证了这些独立频带中状态信息无冗余、无疏漏，排除了干扰，浓缩了了动态分析与监测诊断的信息。
3、 为什么要从内积变换的角度来认识常见的几种信号处理方法？如何选择合适的信号处理方法？ 在信号处理各种运算中内积变换发挥了重要作用。内积变换可视为信号与基函数关系紧密程度或相似性的一种度量。对于平稳信号，是利用傅里叶变换将信号从时域变为频域函数实现的方式是信号函数x （t ）与基函数i t e ω 通过内积运算。匹配出信号x （t ）中圆频率为w 的正弦波.而非平稳信号一般会用快速傅里叶变换、离散小波变换、连续小波变换等这些小波变换的内积变换内积运算旨在探求信号x （t ）中包含与小波基函数最相关或最相似的分量。
“特征波形基函数信号分解”旨在灵活运用小波基函数 去更好地处理信号、提取故障特征。用特定的基函数分解信号是为了获得具有不同物理意义的分类信息。
不同类型的机械故障会在动态信号中反应出不同的特征波形，如旋转机械失衡振动的波形与正弦波形有关，内燃机爆燃振动波形是具有钟形包络的高频波；齿轮轴承等机械零部件出现剥落。裂纹等王府机械活塞连杆、气阀磨损缺陷在运行过程中产生的冲击振动呈现出接近单边震荡衰减波形，等等充分利用基函数的各种性质，根据研究对象的特点和需求，选用针对性强的小波基函数，才能合理地解决工程实际问题，融合表征各种不同类型机械状态特征波形的混合基函数，是现代信号处理进行机械动态分析和检测诊断的一个新的研究方向。
4、 对于基函数的各种性质的物理意义如何理解？
1、 正交性——是小波基函数一个非常优良的性质，他保证信号处理时将信息独立化的提取出来。
2、 正则性——在数学上表现为小波函数的光滑性或可微性。
3、 消失矩——小波基函数的消失矩必须具有足够高的阶数，一个小波消失矩为N ，则它的滤波器长
度不能少于2R 。在信号奇异性检测中要求有足够高的消失矩，但不能过高否则会将奇异的信号平滑掉。表示基函数必行光滑性的程度，R 越大越光滑。
4、 紧支性——函数在区间[a ，b]以外恒为零，支撑区间越小，小波局部化能力越强，越有利于信号
)( ,t b a ψ
点的检测。例如谐波小波，它在频域进制具有完全的盒形频谱。
5、 对称性——具有偶对称或奇对称的尺度函数和小波函数在小波变换信号处理时可得到线性相位
和零相移的分析结果。可进行实时性处理。例如谐波小波，其实属部分为偶对称其负数部分为奇对称，可对信号进行实时性处理。
6、 相似性——这是利用小波分型技术分析信号非平稳性和复杂性的理论基础。
7、 冗余度——表示信号通过某种变换后，由逆变换重建原来信号过程中，基函数所包含重建信息
的过剩量。对信号的重构和图像的恢复有意义。
上述性质应用：
傅里叶基函数：时域无紧支性，但频域有优良的对称性、正交性和紧支性，可得到准确的相位幅值频率。
小波函数：1、haar 小波时间局部化能力强，频率弱。
2、shannon 相反。
3、daubechise 紧支性正交小波，应用广泛但没有解析式，只有离散形式，计算量大。不具备严格的对称性。
4、调频高斯小波：可进行连续小波变换。非正交冗余小波。
5、谐波小波：有解析表达式、频域紧支的正交小波，频域有很好的盒性适合旋转机械的
监测诊断
6、laplace 单边衰减震荡，对齿轮轴承等因缺陷在运行中产生的冲击响应以及旋转机械碰
磨、蒸汽激振等故障特征的提取以及模态分析很有效。
7、hermitian 敏感的识别信号的奇异性。
8、第二代小波：获得与信号更好拼配的期望小波基函数。
第二章 信号的时域分析
1、解释理想滤波器的特点。
信号滤波处理是消除或减弱干扰噪声，保留有用信号的过程。理想的滤波器具有以下四个特点
（1） 理想低通滤波器能使信号中低于频率wc 的各频率分量以同样的放大倍数通过，使高于wc
的频率成分减小为0.高通相反。
（2） 理想的低通滤波器具有矩形幅频特性和线性相位特性。
0() 1 () c p c H j j ωωωδφωτωωω?=
（3） 理想低通滤波器物理上是不可实现的。
（4） Wc 越小时，信号失真大。
2、描述实际滤波器的参数有哪些？其物理含义是什么？
1、截止频率wc ——当滤波器幅值等于
0.707时对应的
频率，也称半功率点。
2、通带边缘频率wp 、阻带边缘频率ws
——划分通
带、过渡带和租代的两个指标。
3、允许的波动量
4、衰减——波动的大小。 3、图示说明采样定理的基本原理。实际测试时如何确定采样频率和数据长度？p29
离散信号的频谱相当于将原信号频谱依次平移ws=2*pi/△t 至各采样脉冲函数对应的频域序列点上然后叠加而成，当△t 太大时ws 过小，移至各采样脉冲函数对应的频域序列点上的频谱会有一部分重叠，导致信号与原信号不一致，称为混叠。因此采样频率ws 必须不小于原信号中的最大频率的2倍。从而t πω
?≤或max 12t f ?≤。实际采样中，一般去ws 为最高频率的2.5~4倍。由于测量信号中高频信号往往是由干扰引起的噪声信号或我们不关注的频谱，因此采样前先需对信号进行低通滤波再定采样频率和时间间隔。
数字信号的分辨率包括时间分辨率和频率分辨率。数字信号的时间分辨率即采样间隔?t ，它反映了数字信号在时域中取值点之间的细密程度。数字信号的频率分辨率为?ω =2π / T ，其中T =N ?t 为数字信号的时间跨度，N 为数字信号的长度。频率分辨率表示了数字信号的频谱在频域中取值点之间的细密程度。因此，当采样频率或采样间隔确定后，增大采样点数就可增加信号的时间长度和频率分辨率。
4、 窗函数为什么会导致频谱泄露？试讨论检测两个频率相近幅度不同的信号 ，选择哪种窗函数比
较合适？p30图
理论上任何信号的长度都是无限的，但任何观测信号的长度都是在有限的时间段内进行的。因此，信号采样过程必须使用窗函数，将无限长信号截断成为有限长度的信号。从理论上看，截断过程是在时域将无限长信号乘以有限时间宽度的窗函数。由卷积定理知在频域内则为信号的频谱与窗函数频谱的卷积。由于窗函数的幅频曲线是一个无线带宽的函数，所以即使原信号为有限带宽信号，截断后信号的频谱也必然是无线带宽的。这就说明信号的能量截断后被扩展了。由此可见信号截断必然会带来一定的误差。泄露取决于窗函数的旁瓣，旁瓣越小，相应的泄露越少。
窗函数选择：1、仅要求获得主瓣的频率——矩形窗。
1、 要求幅值精度高，泄露量小——汉宁窗。
这里由于两个信号频率相近但幅值不同，因此在检测过程中要求幅度进度高，应选择汉宁窗。
5、有量纲指标与无量纲指标各有什么优缺点？试举例说明。
有量纲参数指标：平均幅值、方根幅值、均方幅值、峰值四种。
优缺点：不但与机器的状态有关，且与机器的运动参数如转速、载荷等有关。 无量纲参数指标：波形指标、峰值指标、脉冲指标、裕度指标、偏斜度指标、峭度指标。
优缺点：具有对信号幅值和变化率均不敏感的特点，它们与机械的运动条件无关，只依赖于概率密度函数的形状，是一种较好的机械状态监测诊断参数。
例如，偏斜度指标表示信号概率密度函数的中心偏离正态分布的程度，反应信号的幅值分布和相对其均值的不对称性。峭度指标，表示信号概率密度函数封顶的陡峭程度，反应信号波形中冲击分量的大小。
6.结合你自己的研究方向，谈谈如何应用自相关函数与互相关函数。自己设计一个并解释。 相关是指变量之间的线性联系或相互依赖关系。通过反应信号之间的内积或投影大小来刻画。 自相关函数反应了信号自身取值随自变量时间前后变化的相似性。——可从被噪声干扰的信号中找出周期成分。正常机械噪声是由大量无序、大小近似相等的随机成分叠加的结果，因此正常机器噪声具有较宽而均匀的频谱。当机械状态异常时，随机噪声中将出现有规则、周期性的信号，其幅值要大得多。特别对于早期故障，周期信号不明显是尤为重要。Eg 车床变速箱运行状态识别。 互相关函数完整的描述了两信号之间的相关情况或取值依赖关系。P43例子。
第三章 信号的频域分析
1. 谈谈你对信号频谱的物理本质是如何理解的？结合傅里叶变换的性质，试举例说明其重要作用。 在整个时间轴上的非周期信号()x t 是由频率为ω的谐波()d jwt
X e ωω沿频率从-∞到+∞，通过积分叠加得到的。由于对不同的ω，d ω是一样的，所以只需()X ω就能真实地反映不同高频率谐波的振幅和位移变化。
频谱是信号在频域上的重要特征，反映了信号的频率成分以及分布情况。
7个性质（线性叠加（可分离）、时移性质、频移性质、时间伸缩性、时间微分、时间积分性、卷积定理）。
2. 解释机械信号在离散化过程中产生频率混叠现象及其原因？在工程实践中如何避免频率混叠现象？
机械信号离散化过程中，若采样间隔t ?太大，使得平移距离2t π?过小。移至各采样脉冲函数对应频域序列点上的频谱()X ω就会有一部分重叠，由此造成离散信号的频谱与原信号的频谱不一致，这就导致频率混叠现象。
如果信号中的最高频率（截止频率）为max ω，则在选择采样间隔t ?时应保证max t πω?≤，或max 12t f ?≤，其中max f 是信号中的最高频率（Hz ）。在工程实际中选取采样频率时往往留有余地，一般选取采样频率s ω为处理信号最高频率的2.5~4倍。
3. 在进行信号频谱分析时，为何要加窗函数？如果要求频谱分析结果的幅值精度高，泄漏量小，应该选择什么窗函数？为什么？
理论上信号的长度是无限的，但是任何观测信号都是在有限时间段内进行观测的。因此，信号采样过程中须使用窗函数，将无限长信号截断为有线长度的信号。
如果要求幅值精度高，泄漏量小，应选择汉宁窗函数。因为，汉宁窗函数的旁瓣小，因而相应的泄漏量也较小，采样过程中导致的能量泄露小，能获得较高幅值精度。
要求精确获得主瓣的频率则应选择矩形窗函数。
4. 什么是倒频谱？倒频谱的量纲单位是什么？你如何利用倒频谱原理将时域中两个卷积信号转换为倒频域中相应的两个线性相加的倒频谱？
倒频谱就是对功率谱()x S f 的对数值进行傅里叶逆变换。()()()1log p x
C q F S f -= 其中，q 是倒频率，具有与时间自相关函数()x R τ中的自变量τ相同的时间量纲。q 值大者为高倒频率，表示频谱上的低速波动，q 值小者为低倒频谱，表示频谱上的快速波动。
()()()()()()()()
()()()
()()()()()()
()()()
0111d log log log log log log y x h y x h y x h y x h y t x t h t x h t S f S f S f S f S f S f F S f F S f F S f
C q C q C q τττ∞
---=*=-?=??=+?=+?=+? 工程应用：
（1） 机械中齿轮、滚动轴承等出现故障时信号的频谱上会出现难以识别的多簇调制边频带。采
用到频谱分析可以分解和识别故障频率、原因和部位。
（2） 语音和回声分析及解卷积。有利于问题本质的研究
5. 请说明旋转机械故障诊断中二维全息谱的原理。工频全息谱椭圆较扁说明转子系统存在什么状态现象？
全息谱：将机组的振动信号在完成频域转换后，进一步讲频谱上的谱线加以集成而形成的谱图或轴心轨迹。
将转子测量截面上水平和垂直两方向的振动信号作傅里叶变换，从中提取各主要频率分量的频率、幅值和相位。然后按照各主要频率分量分别进行合成，并将合成结果按频率顺序排列在一张谱图上，就得到了二维全息谱。
工频椭圆扁：支承刚度不对称或受力不均；
椭圆较大、较园则说明转子存在不平衡、轴瓦间隙大或转子永久弯曲等。
工频的二倍频圆比较大、较扁，且工频的四倍频椭圆扁，说明转子存在对中不良、受力不均、基础变形等。
工频的二倍频椭圆较大、较园，说转子存在裂纹或其它故障。
第四章循环平稳信号分析
1.给出循环平稳信号的定义，并解释机械设备循环平稳信号的特点。
某阶统计量随时间变化的信号称为非平稳信号，非平稳信号中的一个重要的子类，其统计量随时间按照周期或多周期规律变化，称这类信号为循环平稳信号。机械设备循环平稳信号的特点：（1）正常无故障的机械信号一般是平稳随机信号，统计量基本不随时间变化。（2）故障信号产生周期成分或调制现象，其统计量呈现周期性变化，此时信号成为循环平稳信号。（3）统计量中的某些周期信息反映机械故障的发生。
2.为什么齿轮、轴承等机械设备在故障发生时，其振动信号往往具有循环平稳性？
由于齿轮、轴承等机械设备零件的结构具有对称性，其运动方式为旋转或往复，因此产生的信号中含有大量的周期成分，其以及统计量出现周期性。此外，当轴承、齿轮等复杂机械零部件出现故障时，振动信号是典型的调制信号，及故障特征频率调制机械系统某阶共振频率，其二阶统计量出现周期性。可见，由于机械结构和工作特点，其产生的非平稳信号具有循环平稳性，因此将机械信号视为循环平稳信号是基于客观实际，是对特定的非平稳信号的一种合理化的假设。
3.对于时间序列为整数，试给出其循环自相关函数的算法步骤。（答案见书本86-88页）
4.如何通过循环谱识别调幅信号的调制频率和载波频率？
根据循环自相关函数三维图和循环自相关函数切片图（τ），可见循环自相关函数将载波信息和调制信息划分到了循环频率高低两个不同的频段，循环频率的高频段的信息既含载波信息又含调制信息，循环频率的低频段只含调制信息，根据这两个频段的信息，可以准确的判别载波信息和调制信息。对于调频信号：根据循环自相关函数三维图和循环自相关函数切片图（τ），可见循环自相关函数将载波信息和调制信息划分到了循环频率高低两个不同的频段，循环频率的高频段的信息既含载波信息又含调制信息，循环频率的低频段只含调制信息，根据这两个频段的信息，可以准确的判别载波信息和调制信息。
对于调福信号，载波信息在频域内的值与其自身相等，而在循环频域内的频率信息是其载波频率的2倍。而调制频率在频率域和循环频率域内的值没有变化。利用循环频率与频率之间的相关性，用切片图可以将有用的信息提取出来并进而分析频率信息特征。
第五章非平稳信号处理方法
1. 请结合时频平面划分的不同，对比说明短时傅里叶变换与小波变换时频分辨率的区别？
（1）受到Heisenberg不确定性原理的约束（公式），短时傅里叶变换的时域分辨率与频域分辨率不可能同时任意小，要提高时间分辨率，只能降低频率分辨率，反之要提高频率分辨率，只能降低时间分辨率。时间分辨率和频率分辨率一旦确定，则STFT在整个时频平面上的时频分辨率保持不变。（2）小波变换则相当于通过小波尺度因子和时移因子变化去观测信号。当a减小时，小波函数的时宽减小，频宽增大；反之亦然。因此，小波变换的局部化是变化的，在高频处时间分辨率高，频率分辨率低，而在低频处，时间分辨率低，频率分辨率高，即具有“变焦”的性质，也就是自适应窗的性质。用大尺度可观察信号的总体，用小尺度可观察信号的细节。
1、请结合时频平面划分的不同，对比说明短时傅里叶变换与小波变换时频分辨率的区别？
小波变换和短时傅立叶变换都是线性时频分析方法。
短时傅立叶变换：用一个在时间上可滑移的时窗来进行傅里叶变换，从而实现时间域和频率域上都具有较好的局部性的分析方法。（STFT是一种时窗大小及形状都固定不变的时频局部化分析）其时间和频率分辨率一旦确定，则在整个时频平面上的时频分辨率保持不变。因为信号的频率与其时间周期成反比，因此反映高频成分需要窄时窗，而反映信号低频成分需要寛时窗，短时傅立叶不能同时满足这些要求。短时傅立叶变换只具有单一分辨率的分析，窗函数一旦确定，其窗口的大小和形状便固定，即窗口只能在同一频率中平移而不能伸缩，因而只能实现一定程度的时频局部化，对信号无法做出仔细的分辨。
小波变换的时频窗口在相平面上不是均匀分布的。其窗口是可调时频窗，在高频时使用窄窗口，在低频时则用宽窗口，即以不同的尺度观察信号，以不同的分辨率分析信号。这充分体现了自适应分辨分析的思想，小波变换能够把任何信号映射到由一个母小波伸缩（变换频率）、平移（刻划时间）而形成的一组基函数上去，实现信号在不同频带、不同时刻的合理分离。在时域和频域都具有很好地局部化特性。小波变换较好地解决了时间和频率分辨率的矛盾，被称为“数学显微镜”
2、解释尺度函数和小波函数的功能，并给出小波分解三层和小波包分解三层的频带划分示意图。
小波函数由基本小波经过伸缩和平移形成的函数（a尺度因子，b时移因子）。一般为双窗函数，具有带通性质。小波函数的时移因子决定了它窗口的时频位置，而他的尺度因子不仅影响窗口在频域的位置也影响了窗口的形状。当尺度减小时，小波函数在时域变窄，而在频域不仅窗口增加，而且中心位置也向高频移动。
尺度函数：通常每个基本小波都与一个尺度函数相对应，甚至有许多小波函数是通过尺度函数而生成。与小波函数相反，尺度函数是一个低通滤波器。由于尺度函数的正交性，小波函数的正交性以及尺度函数和小波函数之间的正交性，通过小波函数（带通）和尺度函数（低通）使信号中的低频信号和高频信号分解到相互独立的频带中，信息无冗余、无疏漏的独立化地提取出来。
频带划分见教材128页
3、解释什么是小波基函数的双尺度关系？为什么小波变换能够对信号进行时间-尺度（时-频）分
析？
双尺度关系：教材122页。V
j 中的尺度函数()
j
t
?和
j
W中的小波函数()
j
t
φ均可由
1
V
j+
中的尺
度函数
1()
j t
?
+
给出。尺度函数可以由比它小一级的尺度函数和函数h离散卷积得到。
4、简述Mallat塔形算法的基本原理和特点。
参见教材127页
它是一种无冗余的迭代算法，计算快，占用存储空间少，同时保存了信号原有的全部信息，因而可以通过逆变换恢复全部信号。
基本原理：
算法过程：在分解计算中，进行隔二抽取，将输入序列每隔一个输出一次（例如只取偶数），组成长度缩短一般的新序列。这样每次分解等获得到的逼近信号和细节信号的数据长度是上一次逼近信号数据长度的一半。当J次分解后，逼近信号和细节信号的数据长度缩减为原始信号长度的2j 。在重构计算的每一步中，先在数据之间插零后在参与同低通、带通滤波器系数的运算，结果重构数据长度加倍。
5、简述小波包频带能量监测的基本原理。教材132页
小波包信号分解是正交分解，各个分解频带的信号互相独立，他们无冗余，不疏漏，小波包信号分解遵循能量守恒原理。小波包信号分解是将包括正弦信号在内的任意信号划归到相应的频带中，用每个频带里的能量来反映机械设备的状态。小波包频带能量监测更具有合理性，通过相应频带里能量比例的变化，可对机械设备进行有效的动态分析与监测诊断。
第六章连续小波变换及其工程应用
1.连续小波变换和离散小波变换的各有哪些优点和缺点？
对于离散小波变换：
三种连续小波分别是谐波小波、Laplace小波和Hermitian小波，它们都具有明确的解析表达式。
谐波小波是复小波，在频域紧支，且具有完全“盒形”的频谱。谐波小波分解算法是通过信号的FFT以及IFFT实现的，算法速度快，精度高，因而具有很好的工程应用价值。随着小波层（即j）的变大，谐波小波的频谱宽度倍增而幅值降低。处于不同层的谐波小波总是正交的，同一层内的小波也是正交的。分析频宽从高频到低频是以1/2关系逐渐减小的，对信号的低频部分划分比较细，而高频部分划分比较粗，这说明谐波小波分解是一种小波分解。小波小波具有“锁定”信号相位的能力。
缺点：谐波小波分解系数，低频频带内的数据点数少，高频频带内的数据点数多。
对于离散小波变换：
小波分析中广泛使用的Daubechies类小波与样条小波都是实小波，它们没有明确的解析表达式，对信号的小波分解是通过构造相应的正交滤波器系数运用Mallat快速算法实现
的。
小波分解可以完全表征函数f(t)的全部信息(提供了一个数学上完备的描述)；
小波变换通过选取合适的滤波器，可以极大的减小或去除所提取得不同特征之间的相关性；
小波变换具有“变焦”特性，在低频段可用高频率分辨率和低时间分辨率(宽分析窗口)，在高频段，可用低频率分辨率和高时间分辨率(窄分析窗口)。
小波变换实现上有快速算法(Mallat 小波分解算法)
小波包技术将信号无冗余、无疏漏、正交地分解到独立的频带内，每个频带里信号的能量对于机械动态分析与监测诊断都是十分有用的信息。
2.以谐波小波变换为例，说明如何实现连续小波变换的快速算法？
Newland 快速算法是指谐波小波的分解算法，通过信号的快速傅里叶变换FFT 和快速傅里叶逆变换IFFT 实现。设有离散信号x (r)，r = 0,…, N-1,其中N =n 2 ，其谐波小波分解为s a , s= 0,…, N – 1。令s FFT(x(r))=F ，s a 由s F 经分段、对每一段作IFFT 得到，下两式为其表达式： j j j 1s s N s s a IFFT F s 2,21,,21a a +-=()
=[][]+???[]-=
其中2j 1,0,,log (N/4)=-???.下图表示一数据长度为16的实序列的谐波小波分解示意图。
系数s a 为一复序列，包含信号时频两方面的信息。s a
某一段对应谐波小波分解的一层，
也就是信号特定频段的成分。s a 与j 的关系如下图所示： (N/2)1023456712a ,a a ,a a ,a ,a ,a ,,a ,,,j 1j 0j 1j 2j log (N/4)
-?????????=-==== 图： s a 的分段与j 的关系
3.为什么小波分析与分形理论可以相结合构成小波分形技术？除了谐波小波轴心轨迹的盒维数应用外，请提出其它的振动信号小波分形应用方法。
分形是一门以不规则事物为研究对象、探索复杂性的科学，因此，很自然地被用来描述设备振动信号的不规则形和复杂性。事实上，分形理论和小波分析在自相似性的本质上和认识事物由粗到细的过程中是一致的。小波分形技术原理是通过比较小波分解后不同频带内信号分形维数的大小及其变化来反映信号的不规则形和复杂性，刻画信号的非平稳性。
小波变换、小波包分析是一种基于事物认识过程的多分辨分析方法，如同人们从远到近逐步深化得观测事物那样，首先看到的是总体轮廓，然后注意到结构线条，最后才聚焦于纹理细节。在振动信号处理中，小波变换、小波包分析可以由粗到细逐步给出振动信号在不同尺度下的波形。这种从低分辨到高分辨的过渡原则与分形过程中的从总体到局部、从宏观到微观深化是一致的。
分形：H (t)r (rt)
r,H 0ββ-=> 小波：1
2a,b t b (t)a ()a,b R,a 0a
ψψ--=∈≠ 分形的自相似仿射算子r 与小波变换的伸缩因子a 是作用相同，因而可以说小波和分形都具有自相似性，上述小波变换、小波包分析与分形在认识事物的过程的一致性和都具备的自相似性保证了两者结合的可行性和必然性。
4.自学Laplace 小波和Hermitian 小波，简要说明其特点与工程应用价值。
Laplace 小波：
该小波在复数空间内为螺旋衰减曲线，其实部和虚部与单自由度结构系统的自由衰减响应函数非常相似。紧支性是显而易见，不具备正交性，其频域盒形不好，故滤波特性较差。
Laplace 小波搜寻信号中的单边衰减波形发生的时刻，振荡频率和阻尼比，实现被测对
象的模态参数识别。Laplace 小波相关滤波法能够在强大噪声或其它干扰中准确捕捉到脉冲响应信号，识别出响应波形的参数。可以预测，Laplace 小波相关滤波法在模态识别和设备故障诊断中具有良好的应用前景
Hermitian 小波：
从Hermitian 小波实部和虚部时域波形可以看出，实部实际上就是Mexico-Hat 小波，它是偶函数，在支撑区域内振荡2次。虚部为奇函数，在支撑区域内振荡1.5次，可见只需要少量离散点即可表达Hermitian 小波，因此，它具有很强的时域局部化能力，这恰好是奇异性检测所需要的。Hermitian 小波的傅里叶变换也是实数，对信号进行卷积变换时不会影响信号的相位。
第七章 基于第二代小波变换的信号处理
1. 与经典小波相比，第二代小波的优势哪些？
它不依赖于傅里叶变换，在时域中完成对双正交小波的构造，具有结构化设计和自适应构造的优点； 构造方法灵活，可以通过提升(lifting scheme)改善小波函数的特性，从而构造出具有期望特性的小波； 不再是某一给定小波函数的伸缩和平移，它适合于不等间隔采样问题的小波构造；
算法简单，运算速度快，占用内存少，执行效率高，可以分析任意长度的信号。
2. 图示说明第二代小波分解的基本原理。
具体细节详见教材196-198
3. 简要说明第二代小波尺度函数和小波函数特性。
图202页，从图中可以看出第二代尺度函数和小波函数都具有紧支性和对称性。小波函数的形状与冲击信号的波形非常相似，可有效提取振动波形中的特征分量。当更新器系数与预测器系数不同时，尺度函数和小波函数的支撑区间和光滑性都发生变化，但形状相似。在工程中可根据这一特点灵活选择与喜好特征匹配的预测器和更新器。
4. 简要说明基于第二代小波包分解的滚动轴承损伤定量识别方法。
7.3.2 滚动轴承损伤定量识别方法
s d
滚动轴承和齿轮发生损伤时，在其缺陷部位产生的冲击脉冲激励下，会出现振荡衰减的脉冲响应信号。理想的方法是利用包络解调技术获取故障特征频率对应的幅值，再利用冲击脉冲法（Shock Pulse Method ，SPM ）对解调结果进行定量处理。可定量识别轴承等部件的损伤程度。SPM 给出dB （分贝）的故障等级经验计算公式：
0.6
200020log
SV B N D ?=? 根据n dB 的如下值判断轴承的运行状态：
[0,21)——正常状态，轴承工作状态良好；
[21,35]——轻微故障，轴承有早期损伤；
(35,60]——严重故障，轴承已有明显损伤。 （1） 将原始振动信号按第二代小波算法进行分解和重构。
（2） 将重构得到的各频带信号进行希尔伯特包络解调然后对信号进行STFT ，得到各频带重构信号
的包络谱。
（3） 计算各个频带重构信号的包络谱分贝值
（4） 计算滚动轴承保持架、滚动体、外圈、内圈的故障特征频率。提取各个故障特征频率在第二
代小波包各个频带中对应包络谱幅值的分贝值。
（5） 选取各频带中同一故障频率对应分贝值得最大值，对滚动轴承的损伤进行定量识别。
5. 为什么要进行冗余第二代小波变换？简述其工程应用。
在机电设备运行过程中，故障特征信息往往被淹没在噪声中，因此降噪是故障诊断领域的重要课题。小波阈值降噪是一种常用的降噪方法，该方法采用下抽样运算，会造成重构后信号失真，且计算量大。为了克服上述缺陷，对信号进行冗余第二代小波变换；该方法，先根据被分析信号的特点，选择相应的更新器和预测器；在每次信号分解前对预测器和更新器进行差值补零语段，获得每层小波分解新的冗余预测器和冗余更新器，该方法去掉了剖分部分直接利用新的冗余预测器和更新器对信号进行分解和重构，能够有效的去除信号中的噪声，较好地保留信号的时域特征；（分解信号具有时移不变性，信息冗余性，重构失真小）
去除振动信号中的噪声，提取被噪声淹没的信号故障特征——降噪
第八章 基于EMD 的时频分析方法及其应用
1. 为什么在EMD 分解时会出现端点效应？试给出三种消除或减弱端点效应的措施。
在应用EMD 方法对非平稳信号进行分解时，在数据的两端会产生发散现象，并且这种发散的结果会逐渐向内“污染”整个数据序列而使所得分解结果严重失真，这就是所谓EMD 分解过程中产生的端点效应问题。
包括直接对原始数据进行简单延托的方法、采用神经网络对数据延托法、在端点出按照端点数据变化的“平衡位置”附加两条平行线段的方法、边界波形匹配预测法、极值点延托法、基于AR 模型的时间序列线性预测方法、神经网络、支持向量机等
2. 瞬时频率的定义和物理意义是什么？如何理由瞬时频率进行转子摩擦故障的诊断？
瞬时频率：相位的微分
物理意义：
转子摩擦故障：此时IMF 的瞬时频率在工频上下波动，也就是说工频的分量出现了频率调制现象。对瞬时频率作傅里叶变换即可得到工频和二倍工频。
这可以解释为转子发生周期性碰磨故障，导致转子转动的线速度发生周期性的变化：转子每转动一周，摩擦一次，线速度都将减小一次，摩擦结束以后又回复到正常速度，因此工频振动分量就发生了上述频率调制现象。
摩擦力可以分解为法向力和切向力，正是这个切向力使转子在1号瓦的波动速度在转动一周的过程中变化一次，从而造成了前文所述1号瓦信号的EMD 分解工频分量的调制现象；
3. 与小波分解相比较，试说明EMD 方法的完备性和正交性。
完备性：就是把分解后的各个分量相加就能获得原信号的性质。根据EMD 的基本原理可知，原信号x(t)为若干基本模式分量与剩余分量的和，可见EMD 方法是完备的。
正交性：到目前为止，经验模式分解的正交性在理论上还难以严格地进行证明，只能在分解后在数值上进行检验。
文献[2] 和[11]分别用某一齿轮箱的振动信号和某一风波信号模式分解的正交性进行了检验，结果证明EMD 方法基本上是正交的，或者称是近似正交的。
整体正交性指标：
两个基本模式分量之间的正交性：
HHT 在傅里叶变换表达式中，在某一频率ω处能量的存在，代表一个正弦或余弦在整个时间长度上都存在。而HHT 中，在某一频率ω处能量的存在，仅代表在数据的整个事件长度上，很可能有这样一个频率的振动在局部出现过。
提高局部均值计算精度和速度的方法有：自适应时变滤波发，极值域均值模式分解法，改进的极值域均值模式分解法。
历年题目答案：
2、 对信号处理这本书你学的方法，这些方法的共同点一级各自的特点。
（1） 时域分析法：滤波、放大、统计特征计算、相关性分析等处理。
特点：可有效的提高信噪比，求取信号波形在不同时刻的相似性和相关性，获得反应机械设备运行状态的特征菜蔬，为机械系统动态分析和故障诊断提供有效的信息。
（2） 频域分析法：
1）传统的分析方法：
112011IO (()()/())
T n n i k t i k c t c t x t ++====∑∑∑,220()()IO ()()T
i k i k t i k c t c t c t c t ==+∑
a相干分析：确定输出信号中分量多大程度的来自输入信号x（他）
b离散傅里叶变换：可真实的反映不同频率谐波的振幅和位移的变化情况，但工作量大，计算时间长。
C复调细化：可将一部分需要重点研究的信号细化放大，提高分辨率。
D倒频谱：将两个卷积信号变成两个倒频谱的线性叠加，从而得到清晰的分析结果。
传统分析方法的缺点：分辨率低、频谱能量泄露、需要较长的原始数据。。并且将幅值和相位信息分离，将相位信息完全的忽略。
2）现代的分析方法：AR、MA/ARMA
优点：谱精度高、适合于短数据、适合分析非平稳信号。
缺点：p阶次难确定。大——峰值发生分叉产生虚假峰值。小——对噪声敏感，要求高。
3）全息谱：将信号各频率分量的频率相位和幅值均反映到谱图上更加全面的反映了转子的振动情况。
实验6 数字滤波器的网络结构 一、实验目的： 1、加深对数字滤波器分类与结构的了解。 2、明确数字滤波器的基本结构及其相互间的转换方法。 3、掌握用MA TLAB 语言进行数字滤波器结构间相互转换的子函数及程序编写方法。 二、实验原理： 1、数字滤波器的分类 离散LSI 系统对信号的响应过程实际上就是对信号进行滤波的过程。因此，离散LSI 系统又称为数字滤波器。 数字滤波器从滤波功能上可以分为低通、高通、带通、带阻以及全通滤波器；根据单位脉冲响应的特性，又可以分为有限长单位脉冲响应滤波器（FIR ）和无限长单位脉冲响应滤波器（IIR ）。 一个离散LSI 系统可以用系统函数来表示： M -m -1-2-m m m=0 012m N -1-2-k -k 12k k k=1 b z b +b z +b z ++b z Y(z)b(z)H(z)=== =X(z)a(z) 1+a z +a z ++a z 1+a z ∑∑ 也可以用差分方程来表示： N M k m k=1 m=0 y(n)+a y(n-k)=b x(n-m)∑∑ 以上两个公式中，当a k 至少有一个不为0时，则在有限Z 平面上存在极点，表达的是以一个IIR 数字滤波器；当a k 全都为0时，系统不存在极点，表达的是一个FIR 数字滤波器。FIR 数字滤波器可以看成是IIR 数字滤波器的a k 全都为0时的一个特例。 IIR 数字滤波器的基本结构分为直接Ⅰ型、直接Ⅱ型、直接Ⅲ型、级联型和并联型。 FIR 数字滤波器的基本结构分为横截型（又称直接型或卷积型）、级联型、线性相位型及频率采样型等。本实验对线性相位型及频率采样型不做讨论，见实验10、12。 另外，滤波器的一种新型结构——格型结构也逐步投入应用，有全零点FIR 系统格型结构、全极点IIR 系统格型结构以及全零极点IIR 系统格型结构。 2、IIR 数字滤波器的基本结构与实现 （1）直接型与级联型、并联型的转换 例6-1 已知一个系统的传递函数为 -1-2-3 -1-2-3 8-4z +11z -2z H(z)=1-1.25z +0.75z -0.125z 将其从直接型（其信号流图如图6-1所示）转换为级联型和并联型。
1、对模拟信号（一维信号，是时间的函数）进行采样后，就是 离散 信号，再进行幅度量化后就是 数字信号。 2、若线性时不变系统是有因果性，则该系统的单位取样响应序列h(n)应满足的充分必要条件是 当n<0时，h(n)＝0 。 3、序列)(n x 的N 点DFT 是)(n x 的Z 变换在 单位圆 的N 点等间隔采样。 4、)()(5241n R x n R x ==，只有当循环卷积长度L ≥8 时，二者的循环卷积等于线性 卷积。 5、已知系统的单位抽样响应为h(n)，则系统稳定的充要条件是 ()n h n ∞ =-∞ <∞∑ 6、用来计算N ＝16点DFT ，直接计算需要（N 2）16*16＝256_次复乘法，采用基2FFT 算法，需要__(N/2 )×log 2N ＝8×4＝32 次复乘法。 7、无限长单位冲激响应（IIR ）滤波器的基本结构有直接Ⅰ型，直接Ⅱ型，_级联型_和 并联型_四种。 8、IIR 系统的系统函数为)(z H ，分别用直接型，级联型，并联型结构实现，其中 并联型的运算速度最高。 9、数字信号处理的三种基本运算是：延时、乘法、加法 10、两个有限长序列 和 长度分别是 和 ，在做线性卷积后结果长度是 __N 1＋N 2－1_。 11、N=2M 点基2FFT ，共有 M 列蝶形，每列有N/2 个蝶形。 12、线性相位FIR 滤波器的零点分布特点是 互为倒数的共轭对 13、数字信号处理的三种基本运算是： 延时、乘法、加法 14、在利用窗函数法设计FIR 滤波器时，窗函数的窗谱性能指标中最重要的是___过渡带宽___与__阻带最小衰减__。 16、_脉冲响应不变法_设计IIR 滤波器不会产生畸变。 17、用窗口法设计FIR 滤波器时影响滤波器幅频特性质量的主要原因是主瓣使数字滤波器存在过渡带，旁瓣使数字滤波器存在波动，减少阻带衰减。 18、单位脉冲响应分别为 和 的两线性系统相串联，其等效系统函数时域及频域表 达式分别是h(n)=h1(n)*h2(n), =H1(ej ω)×H2(ej ω)。 19、稳定系统的系统函数H(z)的收敛域包括 单位圆 。 20、对于M 点的有限长序列x(n)，频域采样不失真的条件是 频域采样点数N 要大于时域采样点数M 。
数字信号处理作业
DFT 习题 1. 如果)(~n x 是一个周期为N 的周期序列，那么它也是周期为N 2的周期序列。把)(~ n x 看作周期为N 的周期序列，令)(~1k X 表示)(~n x 的离散傅里叶级数之系数，再把)(~ n x 看作周期为N 2的周期序列，再令)(~2k X 表示)(~n x 的离散傅里叶级数之系数。当然，)(~1k X 是周期性的，周期为N ，而)(~2k X 也是周期性的，周期为N 2。试利用)(~1k X 确定)(~2k X 。（76-4）
2. 研究两个周期序列)(~n x 和)(~n y 。)(~n x 具有周期N ,而)(~ n y 具有周期M 。序列)(~n w 定义为)()()(~ ~~n y n x n w +=。 a. 证明)(~n w 是周期性的，周期为MN 。 b. 由于)(~n x 的周期为N ，其离散傅里叶级数之系数)(~k X 的周期也是N 。类似地， 由于)(~n y 的周期为M ，其离散傅里叶级数之系数)(~k Y 的周期也是M 。)(~n w 的离散傅里叶级数之系数)(~k W 的周期为MN 。试利用)(~k X 和)(~k Y 求)(~k W 。（76-5）
3. 计算下列各有限长度序列DFT （假设长度为N ）: a. )()(n n x δ= b ．N n n n n x <<-=000) ()(δ c ．10)(-≤≤=N n a n x n （78-7） 4. 欲作频谱分析的模拟数据以10千赫速率被取样，且计算了1024个取样的离散傅里叶变换。试求频谱取样之间的频率间隔，并证明你的回答。（79 -10）
数字信号处理模拟试题一 一、单项选择题(本大题共10小题，每小题2分，共20分) 在每小题列出的四个备选项中只有一个是符合题目要求的，请将其代码填写在题后的括号内。错选、多选或未选均无分。 1.在对连续信号均匀采样时，要从离散采样值不失真恢复原信号，则采样角频率Ωs与信号最高截止频率Ωc应满足关系（A ） A.Ωs>2Ωc B.Ωs>Ωc C.Ωs<Ωc D.Ωs<2Ωc 2.下列系统（其中y(n)为输出序列，x(n)为输入序列）中哪个属于线性系统？（D） A.y(n)=y(n-1)x(n) B.y(n)=x(n)/x(n+1) C.y(n)=x(n)+1 D.y(n)=x(n)-x(n-1) 3.已知某序列Z变换的收敛域为5>|z|>3，则该序列为（D ） A.有限长序列 B.右边序列 C.左边序列 D.双边序列 4.实偶序列傅里叶变换是（A ） A.实偶序列 B.实奇序列 C.虚偶序列 D.虚奇序列 5.已知x(n)=δ(n)，其N点的DFT［x(n)］=X(k)，则X(N-1)=（B） A.N-1 B.1 C.0 D.-N+1 6.设两有限长序列的长度分别是M与N，欲通过计算两者的圆周卷积来得到两者的线性卷积，则圆周卷积的点数至少应取（B ） A.M+N B.M+N-1 C.M+N+1 D.2(M+N) 7.下面说法中正确的是（C） A.连续非周期信号的频谱为周期连续函数 B.连续周期信号的频谱为周期连续函数 C.离散非周期信号的频谱为周期连续函数 D.离散周期信号的频谱为周期连续函数 8.下列各种滤波器的结构中哪种不是IIR滤波器的基本结构？（C ） A.直接型 B.级联型 C.频率抽样型 D.并联型 9.下列关于FIR滤波器的说法中正确的是（C） A.FIR滤波器容易设计成线性相位特性
现代数字信号处理复习题 一、填空题 1、平稳随机信号是指：概率分布不随时间推移而变化的随机信号，也就是说，平稳随机信号的统计特性与起始 时间无关，只与时间间隔有关。 判断随机信号是否广义平稳的三个条件是： （1）x(t)的均值为与时间无关的常数：C t m x =)( （C 为常数） ； （2）x(t)的自相关函数与起始时间无关，即：)(),(),(ττx i i x j i x R t t R t t R =+=； （3）信号的瞬时功率有限，即：∞<=)0(x x R D 。 高斯白噪声信号是指：噪声的概率密度函数满足正态分布统计特性，同时其功率谱密度函数是常数的一类噪 声信号。 信号的遍历性是指：从随机过程中得到的任一样本函数,好象经历了随机过程的所有可能状态,因此,用一个 样本函数的时间平均就可以代替它的集合平均 。 广义遍历信号x(n)的时间均值的定义为： ，其时间自相关函数的定义为： 。 2、连续随机信号f(t)在区间上的能量E 定义为： 其功率P 定义为： 离散随机信号f(n)在区间 上的能量E 定义为： 其功率P 定义为： 注意:(1)如果信号的能量0 数字信号处理期末试卷(含答案) 一、单项选择题(在每小题的四个备选答案中，选出一个正确答案，并将正确答案的序号填在括号内。 1.若一模拟信号为带限，且对其抽样满足奈奎斯特采样定理，则只要将抽样信号通过( )即可完全不失真恢复原信号。 A.理想低通滤波器 B.理想高通滤波器 C.理想带通滤波器 D.理想带阻滤波器 2．下列系统（其中y(n)为输出序列，x(n)为输入序列）中哪个属于线性系统？( ) A.y(n)=x 3(n) B.y(n)=x(n)x(n+2) C.y(n)=x(n)+2 D.y(n)=x(n 2) 3.．设两有限长序列的长度分别是M 与N ，欲用圆周卷积计算两者的线性卷积，则圆周卷积的长度至少应取( )。 A ．M+N B.M+N-1 C.M+N+1 D.2(M+N) 4.若序列的长度为M ，要能够由频域抽样信号X(k)恢复原序列，而不发生时域混 叠现象，则频域抽样点数N 需满足的条件是( )。 A.N ≥M B.N ≤M C.N ≤2M D.N ≥2M 5.直接计算N 点DFT 所需的复数乘法次数与( )成正比。 A.N B.N 2 C.N 3 D.Nlog 2N 6.下列各种滤波器的结构中哪种不是FIR 滤波器的基本结构( )。 A.直接型 B.级联型 C.并联型 D.频率抽样型 7.第二种类型线性FIR 滤波器的幅度响应H(w)特点( )： A 关于0=w 、π、π2偶对称 B 关于0=w 、π、π2奇对称 C 关于0=w 、π2偶对称 关于=w π奇对称 D 关于0=w 、π2奇对称 关于=w π偶对称 8.适合带阻滤波器设计的是： （ ） A )n N (h )n (h ---=1 N 为偶数 B )n N (h )n (h ---=1 N 为奇数 C )n N (h )n (h --=1 N 为偶数 数字信号处理上机作业 学院：电子工程学院 班级：021215 组员： 实验一：信号、系统及系统响应 1、实验目的 (1) 熟悉连续信号经理想采样前后的频谱变化关系，加深对时域采样定理的理解。 (2) 熟悉时域离散系统的时域特性。 (3) 利用卷积方法观察分析系统的时域特性。 (4) 掌握序列傅里叶变换的计算机实现方法，利用序列的傅里叶变换对连续信号、离散信号及系统响应进行频域分析。 2、实验原理与方法 (1) 时域采样。 (2) LTI系统的输入输出关系。 3、实验内容及步骤 (1) 认真复习采样理论、离散信号与系统、线性卷积、序列的傅里叶变换及性质等有关内容，阅读本实验原理与方法。 (2) 编制实验用主程序及相应子程序。 ①信号产生子程序，用于产生实验中要用到的下列信号序列： a. xa(t)=A*e^-at *sin(Ω0t)u(t) b. 单位脉冲序列：xb(n)=δ(n) c. 矩形序列： xc(n)=RN(n), N=10 ②系统单位脉冲响应序列产生子程序。本实验要用到两种FIR系统。 a. ha(n)=R10(n); b. hb(n)=δ(n)+2.5δ(n-1)+2.5δ(n-2)+δ(n-3) ③有限长序列线性卷积子程序 用于完成两个给定长度的序列的卷积。可以直接调用MATLAB语言中的卷积函数conv。 conv 用于两个有限长度序列的卷积，它假定两个序列都从n=0 开始。调用格式如下： y=conv (x, h) 4、实验结果分析 ①分析采样序列的特性。 a. 取采样频率fs=1 kHz,，即T=1 ms。 b. 改变采样频率，fs=300 Hz，观察|X(e^jω)|的变化，并做记录(打印曲线)；进一步降低采样频率，fs=200 Hz，观察频谱混叠是否明显存在，说明原因，并记录(打印)这时的|X(e^j ω)|曲线。 程序代码如下： close all;clear all;clc; A=50; a=50*sqrt(2)*pi; m=50*sqrt(2)*pi; fs1=1000; fs2=300; fs3=200; T1=1/fs1; T2=1/fs2; T3=1/fs3; N=100; 一、选择题 2、对于x(n)=n 21??? ??u(n)的Z 变换，( )。 A. 零点为z=21，极点为z=0 B. 零点为z=2 1 ，极点为z=2 C. 零点为z=21，极点为z=1 D. 零点为z=0，极点为z=21 3、()?? ? ??=n A n x π513sin 是一个以( )为周期的序列。 A. 16 B. 10 C. 14 D. 以上都不对，是一个非周期序列 6、序列()1+n δ的波形图为( )。 C B A 7、s 平面的虚轴对应z 平面的( )。 A. 单位圆内 B. 单位圆外 C. 正实轴 D. 单位圆上 8、关于快速傅里叶变换，下述叙述中错误的是( )。 A.相对离散傅里叶变换来说，它不是一种全新的算法 B.nk N W 具有对称、周期和可约性 C.每个蝶形运算的两个输出值仍放回到两个输入所在的存储器中，能够节 省存储单元 D.就运算量来说，FFT 相对DFT 并没有任何减少 9、下列关于FIR 滤波器的说法中正确的是( )。 A. FIR 滤波器不能设计成线性相位 B. 线性相位FIR 滤波器的约束条件是针对()h n C. FIR 滤波器的单位冲激响应是无限长的 D.不管加哪一种窗，对于FIR 滤波器的性能都是一样的 10、幅度量化、时间离散的的信号是( )。 A. 连续时间信号 B. 离散时间信号 C. 数字信号 D. 模拟信号 11、幅值连续、时间为离散变量的信号是（ ）。 A. 连续时间信号 B. 离散时间信号 C. 数字信号 D. 模拟信号 12、右面的波形图代表序列（ ）。 A. ()34-n R B. ()25+n R C. ()25-n R D. ()24-n R 13、序列()??? ??-=ππ6183cos n A n x 的周期为（ ）。 A. 16 B. 10 C. 14 D. 以上都不对，是一个非周期序列 14、从奈奎斯特采样定理得出，要使信号采样后能够不失真还原，采样频率f 与信号最高频率 f h 关系为：（ ）。 A. f ≤2f h B. f ≥2f h C. f ≥f h D. f ≤f h 16、无限长单位冲激响应（IIR ）滤波器的结构是（ ）型的。 A. 非递归 B. 无反馈 C. 递归 D. 不确定 17、已知序列Z 变换的收敛域为｜z ｜<1，则该序列为（ ）。 A.有限长序列 B. 左边序列 C. 右边序列 D.双边序列 18、下面说法中正确的是（ ）。 A. 连续非周期信号的频谱为周期连续函数 B. 连续周期信号的频谱为周期连续函数 C. 离散周期信号的频谱为周期连续函数 D. 离散非周期信号的频谱为周期连续函数 19、利用矩形窗函数法设计FIR 滤波器时，在理想频率特性的不连续点附近形 成的过滤带的宽度近似等于（ ）。 2013《现代信号处理》试题 1. （10分）设观察样本{x i }(i =1,…,n )的分布密度为 22 2exp{}0(,) 0 0x x i xe e x f x x λλλλ+??->?=?≤?? 其中未知参数0λ>．试求λ的极大似然估计。 2. （30分）现代信号处理与传统的数字信号处理相比，一个最大的区别在于处理的信号是统计性的随机信号而不再是确定性信号，请回答下述问题： （1）当研究宽平稳信号时，需要有各态历经性的理论基础来支撑，请对该性质加以 论述。 （2）白噪声是现代信号处理中常用的一种随机信号，请从时域和频域两个角度对其 加以阐述。 （3）为了便于分析和设计，白化滤波器被提了出来，请从其作用和应用两个方面对 其加以阐述。 3. （30分）与传统的数字信号处理相比，现代信号处理另一个最大的区别在于更多的关注信号之间的关系，如相关函数、功率谱密度函数、信噪比等，请回答下述问题： （1）Wiener 滤波器是现代信号滤波处理的经典，其核心在于考察滤波器输入输出信 号之间的关系，请用恰当的数学模型对其加以描述。 （2）功率谱密度是对时域自相关函数进行傅立叶变换得到的结果。请阐述在工程中 对功率谱密度进行测量有何应用？ （3）高阶谱在传统功率谱的基础上发展起来的，请对其概念、特点与具体应用进行 简要介绍。 4. （15分）梯度搜索法的基本原理是什么？Widrow 提出的LMS 算法与基本的梯度法有何不同？试写出Widrow 提出的LMS 算法的基本步骤。 5. （15分）用计算机仿真计算功率谱，用下式生成一个随机序列 ()2cos(2.02)0.5sin(52)()x t t t e t ππ=?+?+ e (t )为白噪声，均值为零，方差为0.1~1（可任选）或为信号的5%~30%（可任选）。 （1）用周期图法求功率谱估计。 （2）用参数模型法求功率谱估计。 （3）采用Burg 算法求功率谱估计。 数字信号处理作业 哈尔滨工业大学 2006.10 DFT 习题 1. 如果)(~n x 是一个周期为N 的周期序列，那么它也是周期为N 2的周期序列。把)(~ n x 看作周期为N 的周期序列，令)(~ 1k X 表示)(~n x 的离散傅里叶级数之系数，再把)(~ n x 看作周期为N 2的周期序列，再令)(~ 2k X 表示)(~n x 的离散傅里叶级数之系数。当然，)(~ 1k X 是周期性的，周期为N ，而)(~ 2k X 也是周期性的，周期为N 2。试利用)(~ 1k X 确定)(~ 2k X 。（76-4） 2. 研究两个周期序列)(~ n x 和)(~ n y 。)(~ n x 具有周期N ,而)(~ n y 具有周期M 。序列 )(~n w 定义为)()()(~ ~~n y n x n w +=。 a. 证明)(~ n w 是周期性的，周期为MN 。 b. 由于)(~n x 的周期为N ，其离散傅里叶级数之系数)(~ k X 的周期也是N 。类似地， 由于)(~n y 的周期为M ，其离散傅里叶级数之系数)(~k Y 的周期也是M 。)(~ n w 的离散傅里叶级数之系数)(~ k W 的周期为MN 。试利用)(~ k X 和)(~ k Y 求)(~ k W 。（76-5） 3. 计算下列各有限长度序列DFT （假设长度为N ）: a. )()(n n x δ= b ．N n n n n x <<-=000)()(δ c ．10)(-≤≤=N n a n x n （78-7） 4. 欲作频谱分析的模拟数据以10千赫速率被取样，且计算了1024个取样的离散傅里叶变换。试求频谱取样之间的频率间隔，并证明你的回答。（79 -10） 一、 填空题（每题2分，共10题） 1、 1、 对模拟信号（一维信号，是时间的函数）进行采样后，就是 信号，再 进行幅度量化后就是 信号。 2、 2、 )()]([ω j e X n x FT =，用)(n x 求出)](Re[ω j e X 对应的序列 为 。 ３、序列)(n x 的N 点DFT 是)(n x 的Z 变换在 的N 点等间隔采样。 ４、)()(5241n R x n R x ==，只有当循环卷积长度L 时，二者的循环卷积等于线性卷积。 ５、用来计算N ＝16点DFT ，直接计算需要_________ 次复乘法，采用基2FFT 算法，需要________ 次复乘法，运算效率为__ _ 。 ６、FFT 利用 来减少运算量。 ７、数字信号处理的三种基本运算是： 。 ８、FIR 滤波器的单位取样响应)(n h 是圆周偶对称的，N=6， 3)3()2(2 )4()1(5.1)5()0(======h h h h h h ，其幅度特性有什么特性？ ，相位有何特性？ 。 ９、数字滤波网络系统函数为 ∑=--= N K k k z a z H 111)(，该网络中共有 条反馈支路。 １０、用脉冲响应不变法将)(s H a 转换为)(Z H ，若)(s H a 只有单极点k s ，则系统)(Z H 稳定的条件是 （取s T 1.0=）。 二、 选择题（每题3分，共6题） 1、 1、 )6 3()(π-=n j e n x ，该序列是 。 A.非周期序列 B.周期 6π = N C.周期π6=N D. 周期π2=N 2、 2、 序列 )1()(---=n u a n x n ，则)(Z X 的收敛域为 。 A. a Z < B. a Z ≤ C. a Z > D. a Z ≥ 3、 3、 对)70()(≤≤n n x 和)190()(≤≤n n y 分别作20点DFT ，得)(k X 和)(k Y ， 19,1,0),()()(Λ=?=k k Y k X k F ，19,1,0)],([)(Λ==n k F IDFT n f ， n 在 范围内时，)(n f 是)(n x 和)(n y 的线性卷积。 A.70≤≤n B.197≤≤n C.1912≤≤n D.190≤≤n 4、 4、 )()(101n R n x =，) ()(72n R n x =，用DFT 计算二者的线性卷积，为使计算量尽可 能的少，应使DFT 的长度N 满足 。 A.16>N B.16=N C.16 学院________________班级_____________学号________姓名______ 现代信号处理技术试题 一、选择题（下面各题中只有一个答案是正确的，请将正确答案的序号写在每 小题的（）上；每小题2分，共20分） 1. 下列四个离散信号，只有（ ）是周期序列。 A.)100sin(n B. n j e 3 C.)30sin()cos(n n +π D.5432π π j j e e + 2.x(n)非零范围为21N n N ≤≤，h(n)的非零范围为43N n N ≤≤，y(n)=x(n)*h(n) 的非零范围为（ ）。 A.4231N N n N N +≤≤+ B. 42311N N n N N +≤≤-+ C. 14231-+≤≤+N N n N N D. 114231-+≤≤-+N N n N N 3.求周期序列[]?? ? ??=k k x 5cos 2~π的DFS 系数为（ ）。 A.[]???==others m m x 09,12~ B. []???==others m m x 09,110~ C. []???==others m m x 0510~ D. []? ??==others m m x 05,15~ 4.序列[]{}210121，，：，，==k k x 的幅度谱和相位谱为（ ） 。 A.()()02cos 42=ΩΩ=Ωφ，j e X B. ()()Ω-=ΩΩ=Ωφ，2 cos 42j e X C. () ()2 -2cos 42πφ+Ω=ΩΩ=Ω，j e X D. ()()Ω-=Ω=Ωφ，4j e X 5.当序列x[k]为实序列，且具有周期偶对称性，则序列的DFT 满足（ ）。 A.X[m]周期共轭对称 B. X[m]虚部为零，实部周期奇对称 C.X[m]实部为零，虚部周期奇对称 D. X[m]虚部为零，实部周期偶对称 6.与512点的DFT 相比，512点的FFT 只需（ ）。 A.1/2的计算量 B.1/100的计算量 C.2倍的计算量 D.1/10的计算量 7.通带和阻带内均有波纹的IIR 滤波器是（ ）。 A.Butterworth B.Chebyshev I C.Chebyshev II D.椭圆 8.M 阶FIR 滤波器具有线性相位的条件是（ ）。 A. ()()n h n h -= B. ()()n M h n h -±= 数字信号处理作业-答案 数字信号处理作业 DFT 习题 1. 如果)(~ n x 是一个周期为N 的周期序列，那么它也是周期为N 2的周期序列。把)(~ n x 看作周期为N 的周期序列，令)(~ 1 k X 表示)(~ n x 的离散傅里叶级数之系数，再把)(~ n x 看作周期为N 2的周期序列，再令)(~2 k X 表示)(~ n x 的离散傅里叶级数之系数。当然，)(~ 1 k X 是周期性的，周期为N ，而)(~ 2 k X 也是周期性的，周期为N 2。试利用)(~ 1k X 确定)(~ 2 k X 。（76-4） 2. 研究两个周期序列)(~ n x 和)(~ n y 。)(~ n x 具有周期N ,而)(~ n y 具有周期M 。序列)(~ n w 定义为)()()(~~ ~ n y n x n w +=。 a. 证明)(~ n w 是周期性的，周期为MN 。 b. 由于)(~ n x 的周期为N ，其离散傅里叶级数之系数)(~k X 的周期也是N 。类似地，由于)(~ n y 的周期为M ，其离散傅里叶级数之系数)(~ k Y 的周期也是M 。)(~n w 的离散傅里叶级数之系数)(~ k W 的周期为MN 。试利用)(~k X 和)(~k Y 求)(~ k W 。（76-5） 3. 计算下列各有限长度序列DFT （假设长度为N ）: a. )()(n n x δ= b ．N n n n n x <<-=0 0)()(δ c ．10)(-≤≤=N n a n x n （78-7） 4. 欲作频谱分析的模拟数据以10千赫速率被取样，且计算了1024个取样的离散傅里叶变换。试求频谱取样之间的频率间隔，并证明你的回答。（79 -10） 一、填空题（每空1分, 共10分） 1．序列()sin(3/5)x n n π=的周期为 。 2．线性时不变系统的性质有 律、 律、 律。 3．对4()()x n R n =的Z 变换为 ，其收敛域为 。 4．抽样序列的Z 变换与离散傅里叶变换DFT 的关系为 。 5．序列x(n)=(1，-2，0，3；n=0，1，2，3), 圆周左移2位得到的序列为 。 6．设LTI 系统输入为x(n) ，系统单位序列响应为h(n)，则系统零状态输出y(n)= 。 7．因果序列x(n)，在Z →∞时，X(Z)= 。 答案： 1.10 2.交换律，结合律、分配律 3. 4 11,01z z z --->- 4. k N j e Z π2= 5.{0，3，1，-2; n=0,1,2,3} 6.()()()y n x n h n =* 7. x(0) 二、单项选择题（每题2分, 共20分） 1．δ(n)的Z 变换是 （ a ） A.1 B.δ(ω) C.2πδ(ω) D.2π 2．序列x 1（n ）的长度为4，序列x 2（n ）的长度为3，则它们线性卷积的长度是 （ c ） A. 3 B. 4 C. 6 D. 7 3．LTI 系统，输入x （n ）时，输出y （n ）；输入为3x （n-2），输出为 （ b ） A. y （n-2） B.3y （n-2） C.3y （n ） D.y （n ） 4．下面描述中最适合离散傅立叶变换DFT 的是 （ d ） A.时域为离散序列，频域为连续信号 B.时域为离散周期序列，频域也为离散周期序列 C.时域为离散无限长序列，频域为连续周期信号 D.时域为离散有限长序列，频域也为离散有限长序列 5．若一模拟信号为带限，且对其抽样满足奈奎斯特条件，理想条件下将抽样信号通过 即可完 全不失真恢复原信号 （ a ） A.理想低通滤波器 B.理想高通滤波器 C.理想带通滤波器 D.理想带阻滤波器 6．下列哪一个系统是因果系统 （ b ） A.y(n)=x (n+2) B. y(n)= cos(n+1)x (n) C. y(n)=x (2n) D.y(n)=x (- n) 7．一个线性时不变离散系统稳定的充要条件是其系统函数的收敛域包括 （ c ） A. 实轴 B.原点 C.单位圆 D.虚轴 1.设()u n 是离散时间平稳随机过程，证明其功率谱()w 0S ≥。 证明：将()u n 通过冲激响应为()h n 的LTI 离散时间系统，设其频率响应()w H 为 ()001,w -w w 0, w -w w H w ???? 输出随机过程()y n 的功率谱为()()()2y S w H w S w = 输出随机过程()y n 的平均功率为()()()00201 1r 022w w y y w w S w dw S w dw π π π+?-?= =?? 当频率宽度w 0???→时，上式可表示为()()()01 r 00y S w w π =?≥ 由于频率0w 是任意的，所以有()w 0 S ≥ 3、已知：状态方程 )()1,()1()1,()(1n n n n x n n F n x ν-Γ+--=观测方程 )()()()(2n n x n C n z ν+= )()]()([111n Q n n E H =νν )()]()([222n Q n n E H =νν 滤波初值 )]0([)|0(0x E x =ξ } )]]0([)0()]][0([)0({[)0(H x E x x E x E P --= 请简述在此已知条件下卡尔曼滤波算法的递推步骤。 解：步骤1 状态一步预测，即 1 *11)|1(?)1,()|(N n n C n x n n F n x ∈--=--∧ ξξ 步骤2 由观测信号z(n)计算新息过程，即 1*11)|(?)()()|(?)()(M n n C n x n C n z n z n z n ∈-=-=--ξξα 步骤3 一步预测误差自相关矩阵 N N H H C n n n Q n n n n F n P n n F n n P *1)1,()1()1,() 1,()1()1,()1,(∈-Γ--Γ+---=- 步骤4 新息过程自相关矩阵M M H C n Q n C n n P n C n A *2)()()1,()()(∈+-= 步骤5 卡尔曼增益M N H C n A n C n n P n K *1)()()1,()(∈-=- 或 )()()()(1 2n Q n C n P n K H -= 步骤6 状态估计 1*1)()()|(?)|(?N n n C n n K n x n x ∈+=-αξξ 步骤7 状态估计自相关矩阵 N N C n n P n C n K I n P *)1,()]()([)(∈--= 或 )()()()]()()[1,()]()([)(2n K n Q n K n C n K I n n P n C n K I n P H H +---= 步骤8 重复步骤1-7，进行递推滤波计算 4、经典谱估计方法： 中南大学考试试卷 2009-- 2010学年 下 学期期末考试试题 时间100分钟 现代信号处理 课程 56 学时 3.5 学分 考试形式： 闭 卷 专业年级： 2009级 总分100分，占总评成绩 70 % 注：此页不作答题纸，请将答案写在答题纸上 一、填空题 （本题28分，每空2分） 1. 一线性时不变系统，输入为 x （n ）时，输出为y （n ） ；则输入为3x （n-2）时，输出为 。 2. 对连续信号采样时，当采样频率fs 确定情况下，一般在采样前进行预滤波，滤除 的频率成分，以免发生频率混叠现象。 3. 有一模拟系统函数5()2 a H s s =+ ，已知采样周期为T ，采用脉冲响应不变法将其转换为数字系统函数H(z)是 。 4. 设采样频率Hz f s 1000=，则当ω为π/2时，信号的模拟角频率Ω为 。 5. 有限长序列x （n ）的X （k ）与)e (X jw 之间的关系： 6. 单位脉冲响应不变法设计IIR DF 时不适合于 滤波器的设计 7.已知FIR 滤波器4321521----++++=z az z z )z (H 具有线性相位，则a ＝ ,单位脉冲响应h （2）＝ 。 8. 已知一6点实序列x （n ）在4个点上DFT 的值为 ：X （0）=1，X(1)=1+j ，X(3)=3,X(4)=2-j;试写出其它两点的DFT 值X （2）= ，X(5)= 。 9.已知线性相位FIR 数字滤波器的零点为/20.5j z e π=，则可判断该系统函数还具有的零点为: 。 10. 已知序列{}()1,3,2,4;0,1,2,3x n n ==，则序列55(())()x n R n -= 。 11. 已知序列x(n)={4，2，3，1，6,5}，X （K ）为其8点DFT ，则X （4）= ，若38()()k Yk W Xk = ,则y(n)=IDFT[Y(k)]= . 二、选择题(10分，每题2分) 1. 已知序列Z 变换的收敛域为｜z ｜<1，则该序列为( )。 A.有限长序列 B.右边序列 C.左边序列 D.双边序列 2. 计算N=2L （L 为整数）点的按时间抽取基-2FFT 需要( )级蝶形运算。 A ．L B.L/2 C.N D.N/2 3. 下列关于FIR 滤波器的说法中正确的是（ ） 《数字信号处理Ⅰ》作业 姓名： 学号： 学院： 2012 年春季学期 第一章 时域离散信号和时域离散系统 月 日 一 、判断： 1、数字信号处理和模拟信号处理在方法上是一样的。（ ） 2、如果信号的取值和自变量都离散，则称其为模拟信号。（ ） 3、如果信号的取值和自变量都离散，则称其为数字信号。（ ） 4、时域离散信号就是数字信号。（ ） 5、正弦序列都是周期的。（ ） 6、序列)n (h )n (x 和的长度分别为N 和M 时，则)n (h )n (x *的长度为N+M 。（ ） 7、如果离散系统的单位取样响应绝对可和，则该系统稳定。（ ） 8、若满足采样定理，则理想采样信号的频谱是原模拟信号频谱以s Ω（采样频率）为周期进行周期延拓的结果。（ ） 9、序列)n (h )n (x 和的元素个数分别为21n n 和，则)n (h )n (x *有（1n n 21-+）个元素。（ ） 二、选择 1、R N (n)和u(n)的关系为（ ）： A. R N (n)=u(n)-u(n-N) B. R N (n)=u(n)+u(n-N) C. R N (n)=u(n)-u(n-N-1) D. R N (n)=u(n)-u(n-N+1) 2、若f(n)和h(n)的长度为别为N 、M ，则f(n)*h(n)的长度为 （ ）： A.N+M B.N+M-1 C.N-M D.N-M+1 3、若模拟信号的频率范围为[0，1kHz],对其采样，则奈奎斯特速率为（ ）： A.4kHz B. 3kHz C.2kHz D.1kHz 4、LTIS 的零状态响应等于激励信号和单位序列响应的（ ）： A.相乘 B. 相加 C.相减 D.卷积 5、线性系统需满足的条件是（ ）： A.因果性 B.稳定性 C.齐次性和叠加性 D.时不变性 6、系统y(n)=f(n)+2f(n-1)(初始状态为0)是（ ）： A. 线性时不变系统 B. 非线性时不变系统 C. 线性时变系统 D. 非线性时变系统 数字信号处理期末试卷(含答案) 填空题（每题2分，共10题） 1、 1、 对模拟信号（一维信号，是时间的函数）进行采样后，就是 信号，再 进行幅度量化后就是 信号。 2、 2、 )()]([ωj e X n x FT =，用)(n x 求出)](Re[ωj e X 对应的序列 为 。 ３、序列)(n x 的N 点DFT 是)(n x 的Z 变换在 的N 点等间隔采样。 ４、)()(5241n R x n R x ==，只有当循环卷积长度L 时，二者的循环卷积等于线性卷积。 ５、用来计算N ＝16点DFT ，直接计算需要_________ 次复乘法，采用基2FFT 算法，需要________ 次复乘法，运算效率为__ _ 。 ６、FFT 利用 来减少运算量。 ７、数字信号处理的三种基本运算是： 。 ８、FIR 滤波器的单位取样响应)(n h 是圆周偶对称的，N=6， 3)3()2(2 )4()1(5 .1)5()0(======h h h h h h ，其幅 度特性有什么特性？ ，相位有何特性？ 。 ９、数字滤波网络系统函数为 ∑=--= N K k k z a z H 111)(，该网络中共有 条反馈支路。 １０、用脉冲响应不变法将)(s H a 转换为)(Z H ，若)(s H a 只有单极点k s ，则系统)(Z H 稳定的条件是 （取s T 1.0=）。 一、 选择题（每题3分，共6题） 1、 1、 )6 3()(π-=n j e n x ，该序列是 。 A.非周期序列 B.周期 6π = N C.周期π6=N D. 周期π2=N 2、 2、 序列)1()(---=n u a n x n ，则)(Z X 的收敛域为 。 A.a Z < B.a Z ≤ C.a Z > D.a Z ≥ 3、 3、 对)70() (≤≤n n x 和)190()(≤≤n n y 分别作20点DFT ，得)(k X 和)(k Y ， 19,1,0),()()( =?=k k Y k X k F ，19,1,0)],([)( ==n k F IDFT n f ， n 在 范围内时，)(n f 是)(n x 和)(n y 的线性卷积。 A.70≤≤n B.197≤≤n C.1912≤≤n D.190≤≤n 4、 4、 )()(101n R n x =，)()(72n R n x =，用DFT 计算二者的线性卷积，为使计算量尽可 能的少，应使DFT 的长度N 满足 。 A.16>N B.16=N C.16 “现代数字信号处理”复习思考题 变换 1.给出DFT的定义和主要性质。 2.DTFT与DFT之间有什么关系？ 3.写出FT、DTFT、DFT的数学表达式。 离散时间系统分析 1.说明IIR滤波器的直接型、级联型和并联型结构的主要特点。 2.全通数字滤波器、最小相位滤波器有何特点？ 3.线性相位FIR滤波器的h(n)应满足什么条件？其幅度特性如何？ 4.简述FIR离散时间系统的Lattice结构的特点。 5.简述IIR离散时间系统的Lattice结构的特点。 采样 1．抽取过程为什么要先进行滤波，此滤波器应逼近什么样的指标？ 维纳滤波 1．画出Wiener滤波器结构，写出平稳信号下的滤波方程，导出Wiener-Hopf方程。 2．写出最优滤波器的均方误差表示式。 3．试说明最优滤波器满足正交性原理，即输出误差与输入信号正交。 4．试说明Wiener-Hopf方程和Yule-Walker方程的主要区别。 5．试说明随机信号的自相关阵与白噪声的自相关阵的主要区别。 6．维纳滤波理论对信号和系统作了哪些假设和限制？ 自适应信号处理 1．如何确定LMS算法的μ值，μ值与算法收敛的关系如何？ 2．什么是失调量？它与哪些因素有关？ 3．RLS算法如何实现？它与LMS算法有何区别？ 4．什么是遗忘因子，它在RLS算法中有何作用，取值范围是多少？ 5．怎样理解参考信号d(n)在自适应信号处理处理中的作用？既然他是滤波器的期望响应，一般在滤波前是不知道的，那么在实际应用中d(n)是怎样获得的，试举两个应用例子来加以说明。 功率谱估计 1.为什么偏差为零的估计不一定是正确的估计？ 2.什么叫一致估计？它要满足哪些条件？ 3.什么叫维拉－辛钦(Wiener-Khinteche)定理？ 4.功率谱的两种定义。 5.功率谱有哪些重要性质？ 6.平稳随机信号通过线性系统时输入和输出之间的关系。 7.AR模型的正则方程(Yule-Walker方程)的导出。 8.用有限长数据估计自相关函数的估计质量如何？ 9.周期图法谱估计的缺点是什么？为什么会产生这些缺点？ 10.改进的周期图法谱估计有哪些方法？它们的根据是什么？ 11.既然隐含加窗有不利作用，为什么改进周期图法谱估计是还要引用各种窗？ 12.经典谱估计和现代谱估计的主要差别在哪里？ 13.为什么AR模型谱估计应用比较普遍？ 14.对于高斯随机过程最大熵谱估计可归结为什么样的模型？ 15.为什么Levison-Durbin快速算法的反射系数的模小于１？ 16.什么是前向预测？什么是后向预测？ 17.AR模型谱估计自相关法的主要缺点是什么？ 18.Burg算法与Levison-Durbin算法的区别有哪些？数字信号处理期末试卷(含答案)全..
数字信号处理上机作业
数字信号处理期末复习题
2013现代信号处理试题
数字信号处理作业+答案讲解
数字信号处理期末试卷(含答案)
现代信号处理技术试题
数字信号处理作业-答案
数字信号处理期末试题及答案(1)
现代数字信号处理习题
《现代信号处理》2011试卷
数字信号处理作业-2012
数字信号处理期末试卷(含答案)
2012《现代数字信号处理》课程复习...
相关主题
文本预览