一、每题6分，共10题。
1、 试叙述信号分析的不确定原理，并以高斯信号为例解释相关概念。
2、 相对于傅里叶变换，短时傅里叶变换有何特点？窗口应满足什么条件？
3、 相对于信号的谱图，wvd 有何缺点？
4、 什么是小波变换的恒Q 性质？试由此简要说明小波变换的时频分析特点。
5、 试给出能保持信号能量边缘特性的和不能保持信号能量边缘特性的时频变换的例
子。
6、 什么是连续信号的Gabor 展开？实际利用Gabor 展开分析信号时，是采用临界采样
还是过采样？说明理由。
7、 简要叙述Cohen 类时频分布对核函数的要求。
8、 什么是小波包？简要说明如何选取“最佳小波包”。
9、 什么是信号的抽取？什么是信号的插值？提取前以及插值后分别要进行滤波，请给
出各自需要滤波的原因。
10、 小波提升方案与第一代小波构造方法的主要区别是什么？简述小波提升方案
的优点。 二、1
()0t T x t t T ?≤?=?>??,求WVD
三、设(,)x A θτ为信号()x t 的模糊函数，试证明
(,)(0,0)x x A A θτ≤
四、叙述Mallat 分解算法，并给出其二通道滤波器组实现示意图。
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 算法求功率谱估计。
一、 单 项选择题 1. 序列x(n)=Re(e jn π/12 )+I m (e jn π/18 ),周期为( )。 A. 18π B. 72 C. 18π D. 36 2. 设C 为Z 变换X(z)收敛域内的一条包围原点的闭曲线，F(z)=X(z)z n-1 ,用留数法求X(z)的反变换时( )。 A. 只能用F(z)在C 内的全部极点 B. 只能用F(z)在C 外的全部极点 C. 必须用收敛域内的全部极点 D. 用F(z)在C 内的全部极点或C 外的全部极点 3. 有限长序列h(n)(0≤n ≤N-1)关于τ= 2 1 -N 偶对称的条件是( )。 A. h(n)=h(N-n) B. h(n)=h(N-n-1) C. h(n)=h(-n) D. h(n)=h(N+n-1) 4. 对于x(n)= n )21(u(n)的Z 变换，( )。 A. 零点为z=21，极点为z=0 B. 零点为z=0，极点为z=21 C. 零点为z=21，极点为z=1 D. 零点为z=2 1 ，极点为z=2 5、)()(101n R n x =，)()(72n R n x =，用DFT 计算二者的线性卷积，为使计算量尽可能的少，应使DFT 的长度N 满足 。 A.16>N B.16=N C.16
信号时频分析技术及matlab仿真 电路与系统王冠军 201128013926153 摘要:本文介绍了时频分析的一些基础理论,对短时傅里叶变换Wigner-Ville分布做了简单介绍，运用MATLAB语言实现了旨在构造一种时间和频率的密度函数,以揭示信号中所包含的频率分量及其演化特性的wigner-ville分布。并对时频分析方法的优缺点进行了分析。 关键词:时频分析短时傅里叶变换wigner-ville分布 1 引言 基于Fourier变换的传统信号处理技术从信号频域表示及能量的频域分布的角度揭示了信号在频域的特征。但Fourier变换是一种整体变换，只能为人们提供信号在时域或频域的全局特性而无法了解信号频谱随时间变化的情况。因此，需要使用一种时间和频率的联合函数来表示信号，这种表示简称为信号，也就是信号的时频分析。 2 时频分析方法 信号时频分析主要研究非平稳信号或时变信号的频谱含量是怎样随时间变化的。时频分析是当今信号处理领域的一个主要研究热点，目前常用时频分析方法主要有短时傅里叶变换、Gabor展开、小波变换、Wigner-Ville分布。本文主要介绍了短时傅里叶变换和Wigner-Ville分布两种分析方法。 2．1 短时傅立叶变换STFT 从历史上看,信号的时频分析用的最多的是短时傅立叶变换,这种变换的基本思想是用一个窗函数乘时间信号,该窗函数的时宽足够窄,使取出的信号可以被看成是平稳的,然后进行的傅立叶变换可以反映该时宽中的频谱变化规律,如果让窗随时间轴移动,可以得到信号频谱随时间变化的规律。对于时变信号,了解不同时刻附近的频域特征是至关重要的。因此,人们采用时间—频率描述时变信号,将一维的时域信号映射到一个二维的时域平面,全面反映观测信号的时频联合特征。短时傅立叶变换反映了这一思想,对于时变信号,采用某一滑动窗函数截取信号,并认为这些信号是准平稳的,然后,再分别对其进行傅立叶变换,构成时变信号的时变谱。短时傅立叶变换是一种常用的时—频域分析方法,其基本思想
学院________________班级_____________学号________姓名______ 现代信号处理技术试题 一、选择题（下面各题中只有一个答案是正确的，请将正确答案的序号写在每 小题的（）上；每小题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 -±=
1.3 时频分布及其性质 1.3.1 单分量信号与多分量信号 从物理学的角度看，信号可以分为单分量信号和多分量信号两类，而时－频分布的一个主要优点就是能够确定一个信号是单分量的还是多分量的。所谓单分量信号就是在任一时间只有一个频率或一个频率窄带的信号。一般地，单分量信号看上去只有一个山峰（如图 1.2.2），图中所示的是信号)()()(t j e t A t s ?=的时－频表示，在每一个时间，山峰的峰值有明显的不同。如果它是充分局部化的，那么峰值就是瞬时频率；山峰的宽度就是瞬时带宽。一般地，如果)(t z 是信号)(cos )()(t t a t s φ=的解析信号，)(f Z 是)(t z 对应的频谱， 图1.2.2 单分量信号时－频表示及其特征 则其瞬时频率定义如下： )]([arg 21)(t z dt d t f i π= （1.2.1） 与瞬时频率对偶的物理量叫做群延迟，定义如下： )]([arg 21)(f Z dt d f g πτ= （1.2.2） 而多分量信号是由两个（或多个）山峰构成, 每一个山峰都有它自己不同的瞬时 频率和瞬时带宽。（如图1.2.3所示）。 图1.2.3 多分量信号时－频表示及特征
1.3.2 时－频分布定义 Fourier 变换的另一种形式 ?∞ ∞ --=dt e t s f S ft j π2)()( ?∞ ∞ -=df e f S t s tf j π2)()( Cohen 指出，尽管信号)(t z 的时－频分布有许多形式，但不同的时－频分布只是体现 在积分变换核的函数形式上，而对于时－频分布各种性质的要求则反映在对核函数的约束条件上，因此它可以用一个统一形式来表示，通常把它叫做Cohen 类时－频分布，连续时间信号)(t z （)(t z 为连续时间信号)(t s 的解析信号）的Cohen 类时－频分布定义为 ττφτττπdudvd e v u z u z f t P vu f vt j ) (2*),()2 1()21(),(-+-∞ ∞ -∞ ∞ -∞ ∞ --+=?? ? （1.3.1） 式中),(v τφ称为核函数。原则上，核函数可以是时间和频率两者的函数，但常用的核函数与时间和频率无关，只是时延τ和频偏v 的函数，即核函数具有时、频移不变性。这个定义提供了全面理解任何一种时－频分析方法的通用工具，而且能够在信号分析中将信号的一种时－频表示及其性质同另一种时－频表示及其性质联系在一起。进一步可将（1.3.1）简记为 ττφττπdvd e v v A f t P f vt j z )(2),(),(),(+-∞ ∞ -∞ ∞ -? ? = （1.3.2） 式中),(v A z τ是双线性变换（双时间信号）)2 ()2(),(*τ τ τ-+ =t z t z t k z 关于时间t 作 Fourier 反变换得到的一种二维时－频分布函数，称为模糊函数，即 dt e t z t z v A tv j z πτ ττ2*)2 ()2(),(-+=?∞ ∞- （1.3.3） 因为Cohen 类时－频分布是以核函数加权的模糊函数的二维Fourier 变换，所以Cohen 类 时－频分布又称为广义双线性时－频分布。 两个连续信号)(t x ，)(t y 的互时－频分布定义为： ???∞ ∞-∞ ∞--+-∞ ∞ --+= ττφτττπdudvd e v u y u x f t P vu f vt j xy ) (2*),()2 1()21(),( ? ? ∞ ∞-∞ ∞ -+-=dv d e v v A f tv j xy ττφττπ)(2),(),( （1.3.4） 式中 du e u y u x v A vu j xy πτ ττ2*)2 ()2(),(?∞ ∞--+= （1.3.5） 是)(t x 和)(t y 的互模函数。
现代信号处理大作业 姓名：潘晓丹 学号：0140349045 班级：A1403492
作业1 LD 算法实现AR 过程估计 1.1 AR 模型 p 阶AR 模型的差分方程为： )()()(1 n w i n x a n x p i i =-+ ∑=，其中)(n w 是均值为0的白噪声。 AR 过程的线性预测方法为：先求得观测数据的自相关函数，然后利用Yule-Walker 方程递推求得模型参数，再根据公式求得功率谱的估计。 Yule-Walker 方程可写成矩阵形式： ??????? ? ????????= ??????? ? ?? ???? ????????????? ??? ??--+-+--000)()2()1(1) 0() 2()1()()2()0()1()2()1()1()0() 1()() 2()1()0(2 σp a a a r p r p r p r p r r r r p r r r r p r r r r p p p xx xx xx xx xx xx xx xx xx xx xx xx xx xx xx xx 1.2 LD 算法介绍 Levinson-Durbin 算法可求解上述问题，其一般步骤为： 1) 计算观测值各自相关系数p j j r xx ,,1, 0),( =；)0(0xx r =ρ；i=1； 2) 利用以下递推公式运算： ) 1(1,...,2,1),()()()() ()()(2 1111 1 1 1 i i i i i i i i i i i j xx i xx i k i j j i a k j a j a k i a j i r j a i r k -=-=--==-?+ -=-----=-∑ρρρ 3) i=i+1，若i>p ，则算法结束；否则，返回(2)。 1.3 matlab 编程实现 以AR 模型：xn=12xn-1-12xn-2+w(n)为例，Matlab 程序代码如下： clear; clc;
研究生“现代信号处理”课程大型作业 （以下四个题目任选三题做） 1. 请用多层感知器（MLP ）神经网络误差反向传播（BP ）算法实现异或问题（输入为[00;01;10;11]X T =，要求可以判别输出为0或1），并画出学习曲线。其中，非线性函数采用S 型Logistic 函数。 2. 试用奇阶互补法设计两带滤波器组(高、低通互补)，进而实现四带滤波器组；并画出其频响。滤波器设计参数为：F p ＝1.7KHz , F r ＝2.3KHz , F s ＝8KHz , A rmin ≥70dB 。 3. 根据《现代数字信号处理》（姚天任等，华中理工大学出版社，2001）第四章附录提供的数据(pp.352-353)，试用如下方法估计其功率谱，并画出不同参数情况下的功率谱曲线： 1） Levinson 算法 2） Burg 算法 3） ARMA 模型法 4） MUSIC 算法 4. 图1为均衡带限信号所引起失真的横向或格型自适应均衡器(其中横向FIR 系统长M =11), 系统输入是取值为±1的随机序列)(n x ，其均值为零；参考信号)7()(-=n x n d ；信道具有脉冲响应： 1 2(2)[1cos( )]1,2,3()20 n n h n W π-?+=?=???其它 式中W 用来控制信道的幅度失真(W = 2～4, 如取W = 2.9,3.1,3.3,3.5等)，且信道受到均 值为零、方差001.02 =v σ(相当于信噪比为30dB)的高斯白噪声)(n v 的干扰。试比较基 于下列几种算法的自适应均衡器在不同信道失真、不同噪声干扰下的收敛情况(对应于每一种情况,在同一坐标下画出其学习曲线): 1） 横向/格-梯型结构LMS 算法 2） 横向/格-梯型结构RLS 算法 并分析其结果。
1、已知X a (t) 2COS (2 f o t)式中f o =1OOH Z,以采样频率f s =400Hz 对X a (t)进行采样，得 到采样信号X a (t)和时域离散信号X(n)，试完成下面各题： (1) 写出X a (t)的傅里叶变换表示式 X a (j )； (2) 写出X a (t)和x(n)的表达式； (3 )分别求出X a (t)的傅里叶变换和x(n)的傅里叶变换。 解:( 1) j t j t X a (j ) X a (t)e j dt 2cos( o t)e j dt 3、在时域对一有限长的模拟信号以 4KHZ 采样，然后对采到的N 个抽样做N 点DFT ，所得离 散谱线的间距相当于模拟频率 100HZ 。某人想使频率能被看得清楚些，每 50HZ 能有一根谱 线，于是他用8KHZ 采样，对采到的2N 个样点做2N 点DFT 。问：他的目的能达到吗？ 答：不能，因为他忽略了数字频率和模拟频率的区别。 提高采样频率f s ，N 固然大了，数字频率(单位圆)上的样点数确实增加了，但从 模拟频率谱看，样点一点也没有变得密集，这是因为数字频率 2总是对应模拟频率 f s 。 2 f s f s 采样频率由f s 到2 f s 增加一倍，N 也增加一倍，但模拟频率的采样间隔 s s 100Hz 2N N 2 )，不能提 2N 高模拟频率的分辨 (e j 0 t e j 0 t 上式中指数函数和傅里叶变换不存在, X a (j ) 2 [ ( ) (2) x a (t ) X a (t) (t )e j t dt 引入奇异函数 )] 函数，它的傅里叶变换可以表示成: nT) 2cos( 0nT) (t nT) n 2cos( 0nT), 2、用微处理器对实数序列作谱分析, 以下各参数： (1) x(n) 最小记录时间 (2) (3) (4) 解:( 1) T pmin T max N min 要求谱分辨率F 最大取样时间 最少采样点数 在频带宽度不变的情况下将频率分辨率提高一倍的 已知 (2) F 50Hz 1 F 1 T pmin 1 0.02s 50 1 (3) N min 1 s min T P T 0.02s 0.5 10 3s (4) 辩率提高1倍(F 变成原来的12) T p 0.04s N min ~ T 0.5 10 s 频带宽度不变就意味着采样间隔 5OHZ ，信号最高频率1KHz,是确定 N 值。 3 0.5ms 103 40 T 不变，应该使记录时间扩大一倍为 0.04s 实频率分 80 一 2 一点也没有变。所以，增大采样频率，只能提高数字频率的分辨率 (一 N
1.总结学过的滤波器设计方法，用matlab 仿真例子分析不同设计方法的滤波器的性能及适应场合。 答： 1.1模拟低通滤波器的设计方法 1.1.1 Butterworth 滤波器设计步骤： ⑴.确定阶次N ① 已知Ωc 、Ωs 和As 求Butterworth DF 阶数N ② 已知Ωc 、Ωs 和Ω=Ωp()的衰减 Ap 求Butterworth DF 阶数N ③ 已知Ωp 、Ωs 和 Ω=Ωp 的衰减Ap 和As 求Butterworth DF 阶数N 3dB p Ω≠-/10 /1022(/)10 1,(/)101p s A A N N p c s c ΩΩ=-ΩΩ=-则：
⑵.用阶次N 确定 根据公式： 在左半平面的极点即为的极点，因而 1.1.2 切比雪夫低通滤波器设计步骤： ⑴.确定技术指标 归一化： ⑵.根据技术指标求出滤波器阶数N 及： ⑶.求出归一化系统函数 其中极点由下式求出： ()a H s 2,2N ()()a a H s H s -()a H s ,2,,N p Ωp αs Ωs α/1p p p λ=ΩΩ=/s s p λ=ΩΩε0.12 10 1δε=-p δα=
或者由和S 直接查表得 2.数字低通滤波器的设计步骤： （1） 确定数字低通滤波器的技术指标：通带截止频率p ω、通带最大衰减系数 p α、 阻带截止频率ω、阻带最小衰减系数s α。 （2）将数字低通滤波器的技术指标转换成模拟低通滤波器的技术指标。 巴特沃斯： 切比雪夫： N ()a H p /s s p λ=ΩΩ0.1210 1δ ε=-p δα=
中南大学考试试卷 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 滤波器的说法中正确的是（ ）
1、已知0()2cos(2)a x t f t π=式中0f =100HZ,以采样频率s f =400Hz 对()a x t 进行采样，得 到采样信号?()a x t 和时域离散信号()x n ，试完成下面各题： （1）写出()a x t 的傅里叶变换表示式()a X j Ω； （2）写出()a x t 和()x n 的表达式； （3）分别求出()a x t 的傅里叶变换和()x n 的傅里叶变换。 解：（1）000()()2cos()()j t j t a a j t j t j t X j x t e dt t e dt e e e dt ∞∞-Ω-Ω-∞-∞∞ Ω-Ω-Ω-∞Ω==Ω=+??? 上式中指数函数和傅里叶变换不存在，引入奇异函数δ函数，它的傅里叶变换可以表示成：00()2[()()]a X j πδδΩ=Ω-Ω+Ω+Ω （2） 00?()()()2cos()()()2cos(),a a n n x t x t t nT nT t nT x n nT n δδ∞∞=-∞=-∞=-=Ω-=Ω-∞<<∞ ∑∑ 2、用微处理器对实数序列作谱分析，要求谱分辨率50F Hz ≤，信号最高频率1KHz,是确定以下各参数： （1）最小记录时间min p T （2）最大取样时间max T （3）最少采样点数min N （4）在频带宽度不变的情况下将频率分辨率提高一倍的N 值。 解：（1）已知50F Hz ≤ min 110.0250 p T s F = == （2） max 3 min max 1110.52210s T ms f f ====? （3） min 30.02400.510p T s N T s -===? （4）频带宽度不变就意味着采样间隔T 不变，应该使记录时间扩大一倍为0.04s 实频率分辩率提高1倍（F 变成原来的12） min 30.04800.510p T s N T s -===? 3、在时域对一有限长的模拟信号以4KHZ 采样，然后对采到的N 个抽样做N 点DFT ，所得离散谱线的间距相当于模拟频率100HZ 。某人想使频率能被看得清楚些，每50HZ 能有一根谱线，于是他用8KHZ 采样，对采到的2N 个样点做2N 点DFT 。问：他的目的能达到吗？ 答：不能，因为他忽略了数字频率和模拟频率的区别。 提高采样频率s f ，N 固然大了，数字频率（单位圆）上的样点数确实增加了，但从模拟频率谱看，样点一点也没有变得密集，这是因为数字频率π2总是对应模拟频率s f 。 采样频率由s f 到2s f 增加一倍，N 也增加一倍，但模拟频率的采样间隔Hz N f N f s s 10022== 一点也没有变。所以，增大采样频率，只能提高数字频率的分辨率)222(N N ππ→ ，不能提高模拟频率的分辨率。
南邮研究生“现代信号处理”期末课程大作业 （四个题目任选三题做） 1. 请用多层感知器（MLP ）神经网络误差反向传播（BP ）算法实现异或问题（输入为[00;01;10;11]X T =，要求可以判别输出为0或1），并画出学习曲线。其中，非线性函数采用S 型Logistic 函数。 2. 试用奇阶互补法设计两带滤波器组(高、低通互补)，进而实现四带滤波器组；并画出其频响。滤波器设计参数为：F p ＝1.7KHz , F r ＝2.3KHz , F s ＝8KHz , A rmin ≥70dB 。 3. 根据《现代数字信号处理》（姚天任等，华中理工大学出版社，2001）第四章附录提供的数据(pp.352-353)，试用如下方法估计其功率谱，并画出不同参数情况下的功率谱曲线： 1） Levinson 算法 2） Burg 算法 3） ARMA 模型法 4） MUSIC 算法 4. 图1为均衡带限信号所引起失真的横向或格型自适应均衡器(其中横向FIR 系统长M =11), 系统输入是取值为±1的随机序列)(n x ，其均值为零；参考信号)7()(-=n x n d ；信道具有脉冲响应： 1 2(2)[1cos( )]1,2,3()20 n n h n W π-?+=?=???其它 式中W 用来控制信道的幅度失真(W = 2～4, 如取W = 2.9,3.1,3.3,3.5等)，且信道受到均 值为零、方差001.02 =v σ(相当于信噪比为30dB)的高斯白噪声)(n v 的干扰。试比较基 于下列几种算法的自适应均衡器在不同信道失真、不同噪声干扰下的收敛情况(对应于每一种情况,在同一坐标下画出其学习曲线): 1） 横向/格-梯型结构LMS 算法 2） 横向/格-梯型结构RLS 算法 并分析其结果。
数字信号处理在软件无线电A/D中的应用 摘要:讨论了数字信号处理对软件无线电发展的影响及其在软件无线电中的应用。并对在几种软件无线电结构中的应用作了讨论。而A／D变换器是软件无线电的关键器件，本文主要介绍数字信号处理在A/D技术中得应用。在介绍了软件无线电的概述和原理．转后介绍A／D变换器的研究现状和存在的问题，最后讨论了解决的办法并得出结论。 关键词：数字信号处理，软件无线电，AD，中频，射频。 一.引言 1992年 5 月，在美国电信系统会议上首次明确提出了软件无线电的概念。其中心思想就是以一个通用、标准、模块化的硬件平台为依托，通过软件编程来实现无线电台的各种功能，从基于硬件、面向用途的电台设计方法中解放出来。其主要特点有：尽可能多地采用dsp（数字信号处理）技术；开放程度高；适应性强；空中接口可下载。其目的是为了实现不同通信频段，不同的词制方式和数据编码方式的特殊军事电台之间的互相通信，以及延长电台的使用寿命。 近几年来，随着个人通信的迅速发展，在移动通信系统中，多种通信体系并存，比较有代表性的是美国的窄带CDMA，日本的宽带CDMA和欧洲的基于GSM 的TDMA，由于受到各自利益的驱使，他们不可能统一标准，因此为了密切跟踪发展的趋势，延长设备的使用寿命，蜂窝基站的灵活性和兼容性变得十分关键。而采用软件无线电可以顺利的解决标准兼容和灵活性的问题。
软件无线电是指将硬件作为无线通信的基本平台，把尽可能多的无线及个人通信功能 用软件来实现，使得整个系统具有多频带通信，多标准兼容，可重新通过软件再配置等特点，具有很大的灵括性和兼容性，这是继模拟到数字，固定到移动之后无线通信领域的又一次重大突破，被称为第三代移动通信。 本文将首先概述数字信号处理技术在软件无线电中的几个关键技术；接着给出软件无线电的基本概念，并说明这些技术在软件无线电中的作用；然后文中主要讲述AD技术，介绍目前的发展状况和存在问题，比较现在的一些AD技术方法，发现不足，并给出一定的技术解决方案，最后是结论。
2011年的题（大概） P29采样、频率混叠，画图说明 P33列举时域参数（有量纲和无量纲），说明其意义与作用 P37～自相关互相关及作用（举例说明） P51～蝶形算法 P61频谱细化过程，如何复调制 P67Hilbert 变换过程，瞬时频率 循环平稳信号，调频调幅信号边频带的分析 小波双尺度方程 P128下方的图 第六章三种连续小波的原理性质及应用 P157算法图示 P196图7.1.1和图7.1.2 P219EMD 基本流程 P230端点效应的处理 2012年1月9日现代信号处理试题（无敌回忆版） 一、必选题 1.请说明基函数在信号分解与特征提取中的作用。 2.什么是信号的相关分析？试举一例说明其工程应用。 3.什么是倒频谱？倒频谱的量纲单位是什么？如何利用倒频谱实现时域信号解卷积？ 4.解释尺度函数和小波函数的功能，并给出小波分解三层和小波包分解三层的频带划分示意图。 5.试举例说明将任意2种信号处理方法相结合的特征提取技术及其故障诊断工程应用案例。 二、选答题（7选5） 1.请列出你认为重要的小波基函数两种性质，说明理由。 2.解释机械信号在离散化过程中产生的频率混叠现象及其原因？在实践中如何避免发生频率混叠现象？ 3.试说明旋转机械故障诊断中二维全息谱的原理，工频全息谱椭圆较扁说明转子系统存在什么状态现象？ 4.以五点序列为例，给出预测器系数为N=2，更新器系数为2=-N 时的第二代小波分解图。 5.给出经验模式分解（EMD ）的基本过程，并分析出现端点效应的原因与两种减弱或消除端点效应的措施。 6.给出循环平稳信号的定义，并列出机械设备循环平稳信号的特点。 7.根据你的学习体会，谈谈实现故障定量诊断的重要性，并举例说明某一种故障定量诊断方法。
AR 模型的功率谱估计BURG 算法的分析与仿真 钱平 （信号与信息处理 S101904010） 一．引言 现代谱估计法主要以随机过程的参数模型为基础，也可以称其为参数模型方法或简称模型方法。现代谱估计技术的研究和应用主要起始于20世纪60年代，在分辨率的可靠性和滤波性能方面有较大进步。目前，现代谱估计研究侧重于一维谱分析，其他如多维谱估计、多通道谱估计、高阶谱估计等的研究正在兴起，特别是双谱和三谱估计的研究受到重视，人们希望这些新方法能在提取信息、估计相位和描述非线性等方面获得更多的应用。 现代谱估计从方法上大致可分为参数模型谱估计和非参数模型谱估计两种。基于参数建摸的功率谱估计是现代功率谱估计的重要内容，其目的就是为了改善功率谱估计的频率分辨率，它主要包括AR 模型、MA 模型、ARMA 模型，其中基于AR 模型的功率谱估计是现代功率谱估计中最常用的一种方法，这是因为AR 模型参数的精确估计可以通过解一组线性方程求得，而对于MA 和ARMA 模型功率谱估计来说，其参数的精确估计需要解一组高阶的非线性方程。在利用AR 模型进行功率谱估计时，必须计算出AR 模型的参数和激励白噪声序列的方差。这些参数的提取算法主要包括自相关法、Burg 算法、协方差法、 改进的协方差法，以及最大似然估计法。本章主要针对采用AR 模型的两种方法：Levinson-Durbin 递推算法、Burg 递推算法。 实际中，数字信号的功率谱只能用所得的有限次记录的有限长数据来予以估计，这就产生了功率谱估计这一研究领域。功率谱的估计大致可分为经典功率谱估计和现代功率谱估计，针对经典谱估计的分辨率低和方差性能不好等问题提出了现代谱估计，AR 模型谱估计就是现代谱估计常用的方法之一。 信号的频谱分析是研究信号特性的重要手段之一，通常是求其功率谱来进行频谱分析。功率谱反映了随机信号各频率成份功率能量的分布情况，可以揭示信号中隐含的周期性及靠得很近的谱峰等有用信息，在许多领域都发挥了重要作用。然而，实际应用中的平稳随机信号通常是有限长的，只能根据有限长信号估计原信号的真实功率谱，这就是功率谱估计。 二．AR 模型的构建 假定u(n)、x(n)都是实平稳的随机信号，u(n)为白噪声，方差为 ，现在，我们希望建立AR 模型 的参数和x(n)的自相关函数的关系，也即AR 模型的正则方程(normal equation)。 由 )}()]()({[)}()({)(1 n x m n u k m n x E m n x n x E m p k k x a r ++-+-=+=∑= )()()(1 m k m m r r a r xu x p k k x +--=∑= (1) 由于u(n)是方差为 的白噪声，有 ?? ?=≠=-0 00)}()({2 m m m n x n u E σ （2） 由Z 变换的定义， ，当 时，有h(0)=1。综合（1）及（2）两式， ???????=-≥--=∑∑==0)(1)()(1 2 1 m k m k m m p k x k p k x k x r a r a r σ （3） 在上面的推导中，应用了自相关函数的偶对称性。上式可写成矩阵式：
研究生“现代信号处理”课程大型作业 (以下四个题目任选三题做 1. 请用多层感知器(MLP 神经网络误差反向传播(BP 算法实现异或问题(输入为 [00;01;10;11]X T =,要求可以判别输出为0或1 ,并画出学习曲线。其中,非线性函数采用S 型Logistic 函数。 2. 试用奇阶互补法设计两带滤波器组(高、低通互补,进而实现四带滤波器组;并画出其频响。滤波器设计参数为:F p =1.7KHz , F r =2.3KHz , F s =8KHz , A rmin ≥70dB 。 3. 根据《现代数字信号处理》(姚天任等,华中理工大学出版社,2001第四章附录提供的数据(pp.352-353,试用如下方法估计其功率谱,并画出不同参数情况下的功率谱曲线: 1 Levinson 算法 2 Burg 算法 3 ARMA 模型法 4 MUSIC 算法 4. 图1为均衡带限信号所引起失真的横向或格型自适应均衡器(其中横向FIR 系统长M =11, 系统输入是取值为±1的随机序列(n x ,其均值为零;参考信号7((-=n x n d ;信道具有脉冲响应: 12(2[1cos(]1,2,3(20 n n h n W π-?+=?=???其它 式中W 用来控制信道的幅度失真(W = 2~4, 如取W = 2.9,3.1,3.3,3.5等,且信道受到均
值为零、方差001.02=v σ(相当于信噪比为30dB的高斯白噪声(n v 的干扰。试比较基 于下列几种算法的自适应均衡器在不同信道失真、不同噪声干扰下的收敛情况(对应于每一种情况,在同一坐标下画出其学习曲线: 1 横向/格-梯型结构LMS 算法 2 横向/格-梯型结构RLS 算法 并分析其结果。 图1 横向或格-梯型自适应均衡器 参考文献 [1] 姚天任, 孙洪. 现代数字信号处理[M]. 武汉: 华中理工大学出版社, 2001 [2] 杨绿溪. 现代数字信号处理[M]. 北京: 科学出版社, 2007 [3] S. K. Mitra. 孙洪等译. 数字信号处理——基于计算机的方法(第三版[M]. 北京: 电子工
一、 基本概念填空 1、 统计检测理论是利用 信号 与 噪声 的统计特性等信息来建立最佳判决的数学理论。 2、 主要解决在受噪声干扰的观测中信号有无的判决问题 3、 信号估计主要解决的是在受噪声干扰的观测中，信号参量 和 波形 的确定问题。 4、 在二元假设检验中，如果发送端发送为H 1，而检测为H 0，则成为 漏警 ，发送端发送H 0，而检测为H 1，则称为 虚警 。 5、 若滤波器的冲激响应时无限长，称为 IIR 滤波器，反之，称为 FIR 滤波器 6、 若滤波器的输出到达 最大信噪比 成为 匹配 滤波器；若使输出滤波器的 均方估计误差 为最小，称为 维纳 滤波器。 7、 在参量估计中，所包含的转换空间有 参量空间 和 观测空间 8、 在小波分析中，小波函数应满足 ∫φφ(tt )ddtt =0+∞?∞ 和 ∫|φφ(tt )|ddtt =1+∞ ?∞ 两个数学条件。 9、 在小波的基本概念中，主要存在 F (w )=∫ff (tt )ee ?ii ii ii ddtt +∞?∞和f(t)=12ππ∫FF (ww )ee ii ii ii ddww +∞?∞ 两个基本方程。(这个不确定答案，个人感觉是) 10、 在谱估计中，有 经典谱估计 和 现代谱估计 组成了完整的谱估计。 11、 如果系统为一个稳定系统，则在Z 变换中，零极点的分布
应在单位圆内，如果系统为因果系统，在拉普拉斯变换中， 零极点的分布应在左边平面。 二、问题 1、在信号检测中，在什么条件下，使用贝叶斯准则，什么条 件下使用极大极小准则？什么条件下使用Neyman-Pearson准 则？ 答：先验概率和代价函数均已知的情况下，使用贝叶斯准则，先验概率未知，但可选代价函数时，使用极大极小准则，先验 概率和代价函数均未知的情况下，使用Neyman-Pearson准则。 2、在参量估计中，无偏估计和渐进无偏估计的定义是什么？ 答：无偏估计：若估计量的均值等于被估计量的均值（随机变 量），即E?θθ??=EE(θθ)或等于被估计量的真值（非随机参 量）E?θθ??=θθ，则称θθ?为θ的无偏估计。 渐进无偏估计：若lim NN→∞EE?θθ??=EE(θ ),称θθ?为θ的渐进无偏估计。 3、卡尔曼滤波器的主要特征是什么？ 答：随机过程的状态空间模型，用矩阵表示，可同时估计多参 量，根据观测数据，提出递推算法，便于实时处理。 4、在现代信号处理中，对信号的处理通常是给出一个算法， 对一个算法性能的评价，应从那些方面进行评价。 答：算法的复杂度，算法的稳定性和现有算法的比较，算法的 运算速度、可靠性、算法的收敛速度。
P29采样、频率混叠，画图说明 将连续信号转换成离散的数字序列过程就是信号的采样。 它包含了离散和量化两个主要步骤。 若采样间隔Δt 太大，使得平移距离2π/Δt 过小。移至各采样脉冲函数对应频域序列点上的频谱X(ω)就会有一部分相互重叠， 由此造成离散信号的频谱与原信号频谱不一致，这种现象称为混叠。 P33列举时域参数（有量纲和无量纲），说明其意义与作用。 有量纲参数指标包括方根幅值、平均幅值、均方幅值和峰值四种。 无量纲参数指标包括了波形指标、峰值指标、脉冲指标和裕度指标。 偏斜度指标S 表示信号概率密度函数的中心偏离正态分布的程度，反映信号幅值分布相对其均值的不对称性。 峭度指标K 表示信号概率密度函数峰顶的陡峭程度，反映信号波形中的冲击分量的大小。 P37～自相关互相关及作用（举例说明） 相关，就是指变量之间的线性联系或相互依赖关系。 信号x (t )的自相关函数： 信号中的周期性分量在相应的自相关函数中不会衰减，且保持了原来的周期。因此，自相关函数可从被噪声干扰的信号中找出周期成分。 在用噪声诊断机器运行状态时，正常机器噪声是由大量、无序、大小近似相等的随机成分叠加的结果，因此正常机器噪声具有较宽而均匀的频谱。当机器状态异常时，随机噪声中将出现有规则、周期性的信号，其幅度要比正常噪声的幅度大得多。 依靠自相关函数就可在噪声中发现隐藏的周期分量，确定机器的缺陷所在。 （如：自相关分析识别车床变速箱运行状态，确定存在缺陷轴的位置；确定信号周期。） 互相关函数： 互相关函数的周期与信号x(t)和y(t)的周期相同，同时保留了两个信号的相位差信息φ。可在噪音背景下提取有用信息；速度测量；板墙对声音的反射和衰减测量等。 （如：利用互相关分析测定船舶的航速；探测地下水管的破损地点。P42） P51～蝶形算法 FFT 的基本思想是把长度为2的正整数次幂的数据序列｛x k ｝分隔成若干较短的序列作DFT 计算，用以代替原始序列的DFT 计算。然后再把他们合并起来。得到整个序列｛x k ｝DFT 。(图示N=8时FFT) t t x t x T R T T x d )()(1lim )(0 ? ±=∞ →ττt t y t x T R T T xy d )()(1lim )(0 ? +=∞ →ττ x 0 x 1x 2x 3x 4x 5x 6 x 7x 0x 4x 6x 3x 5 x 0x 4x 2x 6x 1x 5x 3x 7x 0x 4x 2x 6x 1x 5x 3x 7x'0 x'4 x' 2x'6 x'1 x'5x'3 x'7 -1 -1 -1 -1 -1 -1 -1-1 -1 -1 -1 -1 X 0 X 1 X 2X 3 X 4X 5X 6X 7 x 7x 1x 2N W N W N W 0N W 0N W N W N W 1 N W 1 N W 1 N W 0N 0N W 2N W 3