现代信号处理复习题


    现代信号处理复习题 1、已知式中=100HZ,以采样频率=400Hz对进行采样,得到采样信号和时域离散信号,试完成下面各题: (1)写出的傅里叶变换表示式; (2)写出和的表达式; (3)分别求出的傅里叶变换和的傅里叶变换。 解:(1) 上式中指数函数和傅里叶变换不存在,引入奇异函数函数,它的傅里叶变换可以表示成: (2) 2、用微处理器对实数序列作谱分析,要求谱分辨率,信号最高频率1KHz,是确定以下各参数: (1)最小记录时间 (2)最大取样时间 (3)最少采样点数 (4)在频带宽度不变的情况下将频率分辨率提高一倍的N值。 解:(1)已知 (2) (3) (4)频带宽度不变就意味着采样间隔不变,应该使记录时间扩大一倍为0.04s实频率分辩率提高1倍(变成原来的) 3、在时域对一有限长的模拟信号以4KHZ采样,然后对采到的N个抽样做N点DFT,所得离散谱线的间距相当于模拟频率100HZ。某人想使频率能被看得清楚些,每50HZ能有一根谱线,于是他用8KHZ采样,对采到的2N个样点做2N点DFT。问:他的目的能达到吗? 答:不能,因为他忽略了数字频率和模拟频率的区别。 提高采样频率 , 固然大了,数字频率(单位圆)上的样点数确实增加了,但从模拟频率谱看,样点一点也没有变得密集,这是因为数字频率总是对应模拟频率 。采样频率由到2 增加一倍,也增加一倍,但模拟频率的采样间隔 一点也没有变。所以,增大采样频率,只能提高数字频率的分辨率 ,不能提高模拟频率的分辨率。 4、在A/D变换之前和D/A变换之后都要让信号通过一个低通滤波器,他们分别起什么作用? 解:在 变换之前让信号通过一个低通滤波器,是为了限制信号的最高频率,使其满足当采样频率一定时,采样频率应大于等于信号最高频率2倍的条件。此滤波器亦称为“抗折叠”滤波器。 在 变换之后都要让信号通过一个低通滤波器,是为了滤除高频延拓谱,以便把抽样保持的阶梯形输出波平滑化,故又称为“平滑”滤波器。 5、已知,分析其因果性和稳定性。 解: 的极点为, (1) 收敛域,对应的系统是因果系统,但由于收敛域不包含单位圆,因此是不稳定系统。单位脉冲响应,这是一个因果序列,但不收敛。 (2) 收敛域,对应的系统是非因果且不稳定系统。其单位脉冲响应,这是一个非因果且不收敛的序列。 (3) 收敛域,对应的系统是一个非因果系统,但由于收敛域包含单位圆,因此是稳定系统。其单位脉冲响应,这是一个收敛的双边序列。 6、什么叫做数字滤波器?FIR和IIR的比较和各自的设计方案? 答:所谓数字滤波器,是指输入、输出均为数字信号,通过一定运算关系改变信号所含频率成分的相对比例或者滤除某些频率成分的器件。 FIR:有限脉冲响应滤波器 IIR:无限脉冲响应滤波器 ★ IIR极点可存在与单位圆的任何地方,有较强的幅度选择性,但相位特性差。 FIR相位呈线性,但幅度特性需高阶才可调节的较好。 ★ FIR计算不产生振荡,误差影响小,可以采用FFT算法。 IIR有稳定问题,有限字长可能产生振荡,同阶递归算法速度受到限制。 ★ IIR可用模拟滤波器成果,得到有效的封闭式公式,设计工作量小,要求低。 FIR仅窗函数有公式,但无显式表达通、阻带,需要计算机辅助设计。 ★ IIR设计已规格化,频率特性为分段常数的滤波器。 FIR主要适应特殊应用,且高阶IIR不易达到指标的滤波器。 IIR数字滤波器设计 ★直接设计: 原型变换(由一低通经过频率变形设计低通、高通、带通、带阻等) 频域设计(零、极点配置;幅度平方函数), 时域设计(帕德(Pade)逼近;波形形成) ★ 优化技术设计(依据一定的优化准则进行设计) ★ FIR数字滤波器设计 ★线性相位:   零点的镜像存在。 偶对称: 奇对称: ★窗函数(时域加权平均):矩形,三角,余弦,布莱克曼(Blackman)系列,凯塞(Kaiser)系列,高斯 ★频率取样:在H(z)的单位圆上等分取样(是否带初相) ★优化技术设计:(依据一定的优化准则进行设计) 8:长度为N=10的两个有限长序列 作图表示、和(圆周卷积),循环卷积区间长度L=10。 解:、和分别如题3解图()、()、()所示 9:若序列是因果序列,其傅里叶变换的实部为,求序列的及其傅里叶变换 。 解: 10、什么是宽平稳随机过程?什么是严平稳随机过程?它们之间有什么联系? 答:若一个随机过程的数学期望与时间无关,而其相关函数仅与有关,则称这个随机过程是宽平稳的或广义平稳的。所谓严平稳随机过程是指它的任何n维分布函数或概率密度函数与时间起点无关。严平稳的随机过程一定是宽平稳的,反之则不然。 取模的平方 FFT 观测数据x(n) 15、如图所示: 1/N (1)在描述随机信号的频率特性时为什么不用信号的傅里叶变换而改用功率谱估计? (2)观察上述框图,说出这是哪一种经典功率谱估计的方法,并写出描述估计关系式。 (3)根据维纳-辛钦定理及相关估计方法写出另一种经典功率谱估计描述估计关系式,结合框图或关系式说明上述框图所示方法的优点。 (4)两种经典功率谱估计都有一个致命的缺点,请简要说明并写出常用的改进方法的名称。 解:1.对于随机信号,其傅里叶变换并不存在,因此转向研究其功率谱。 2.图中所示的是周期图法 3. 周期图法简单,不用估计自相关函数,且可以用FFT进行计算。 4.经典谱估计得致命缺点是频率分辨率低,其原因是傅里叶变换域是无限大,而用作估计的观察数据只有有限个,认为剩余的数据为0,造成系统偏差。改进的方法有:1.平均周期法2.窗函数法3.修正的周期图求平均法。 16、如图所示的RC电路,若输入电压的功率谱密度为X(),求输出电压的功率谱密度 Y()。 R C Y() X() 解:RC电路系统的频率响应函数为 H() = = H()= 由线性系统的输出谱密度与输入谱密度之间的关系可得: Y() = H()* X()= 17、已知LTI系统的传输函数为h(t),输入是实平稳随机过程X(t),输出是Y(t),求三者间的关系? 解:平稳随机过程经过LTI系统输出还是平稳随机过程,所以 其中是卷积运算。 18、常用的自适应滤波理论与算法有哪些? 从理论上讲,自适应滤波问题没有惟一的解。为了得到自适应滤波器及其应用系统,可以采用各种不同的递推算法,这些自适应算法都有各自的特点,适用于不同场合。 常用的自适应滤波理论与算法有: (1)、基于维纳滤波理论的方法。 (2)、基于卡尔曼滤波理论的方法。 (3)、基于最小均方误差准则的方法。 (4)、基于最小二乘准则的方法。 19、简述自适应信号处理技术的应用 自适应滤波处理技术可以用来检测平稳的和非平稳的随机信号。常应用于: (1)、自适应滤波与逆滤波。 (2)、系统辨识。 (3)、自适应均衡。 (4)、自适应回波抵消。 (5)、自适应噪声抵消与谱线增强。 (6)、自适应谱估计。 (7)、自适应波束形成。 (8)、自适应神经智能信息处理。 (9)、盲自适应信号处理。 20、设为一随机电报信号,其样本函数如图1所示,取+1,-1概率相等,在时间间隔内波形变号次数服从参数为的泊松分布,即: 求的自相关函数。 1 0 -1 t 解: 在时间间隔内可能变号偶次,,将同时取+1或-1,若变号奇次,,将异号。 当时, 显然当时, 21、设,其中和是相互独立的随机变量,在(0,)均匀分布,试讨论的平稳性和各态历经性。 解: 所以,是广义平稳过程。 因此,具有均值各态历经性。 因此,不具有自相关函数各态历经性。 22、从最速下降法出发: 其中, 是第j+1个抽样时刻的滤波器权矢量, 控制收敛稳定性和速率, 是误差-性能曲面的真实梯度,推导自适应噪声消除的 Widrow-Hopf 的LMS算法。 解答: 梯度矢量▽,初级输入与刺激输入的互相关P以及初级输入的自相关R之间的关系为: = 在LMS算法中,使用的瞬时估计,则有 =-2+2=-2+2 (1) =- 其中 用(1)式替换最速下降法的梯度,我们得到基本的Widrow-Hopf 的LMS算法: 其中 = 23、自适应滤波器的特点及应用范围 答案:由于滤波器的参数可以按照某种准则自动地调整到满足最佳滤波的要求;实现时不需要任何关于信号和噪声的自相关特性,尤其当输入统计特性变化时,自适应滤波器都能调整自身的参数来满足最佳滤波的需要,即具有学习和跟踪的性能。当符合下面几个情况时都可以应用自适应滤波(1)需要滤波器特性变化以自适应改变的情况时(2)当信号和噪声存在频谱重叠时(3)噪声占据的频谱是时变或未知。例如电话回声对消,雷达信号处理,导航系统,通信信道均衡和生物医学信号增强。 24、已知输入信号向量U(n)的相关矩阵R及期望响应信号的互相关向量P分别为 R= , 且已知期望响应的平均功率为 E{}=30 (1) 计算维纳滤波的最优权向量; (2) 推导误差性能面的表达式; (3) 计算最小均方误差。 答案:(1)由可得 == (2)假设在时刻,输入信号为,横向滤波器输出信号 (3)把代入得: 25、怎样判断随机过程是宽平稳随机过程?并证明随机过程是宽平稳过程,其中,Y , Z是相互独立的随机变量,且。 答:(1)如果满足,如下条件: (a)是二阶矩过程; (b)对任意,常数; (c)对任意,。则判定是宽平稳随机过程。 (2) 证明: 因为Y,Z是相互独立的随机过程,且,所以 = =常数 ,只与时间间隔有关,与时间起点无关。 所以,是宽平稳随机过程。 26、若为均方连续的实平稳随机过程,则其自相关函数具有那些常用性质?在计算其功率谱时有什么作用? 答:(1)具有如下常用性质: (a) (b)=,是实偶函数; (c)||; (d)若是周期为T的周期函数,即=,则; (e)若是不含周期分量的非周期过程,当时,与相互独立,则 。 (2)若,根据辛钦—维纳定理 = 自相关函数和功率谱是一对傅里叶变换对。 27、从随机过程的平稳性上考虑,卡尔曼滤波的适用范围? 答案:卡尔曼滤波不仅适用于平稳随机过程,同样也适用于非平稳随机过程。 28、简述基于卡尔曼滤波理论的方法。 答: 为使自适应滤波器能工作在工作在平稳或非平稳的环境,可以借助于卡尔曼滤波器来推导自适应滤波算法。卡尔曼滤波是线性无偏最小方差递推滤波,估计性能是最优的,且递推计算形式又能适应实时处理的需要对于一个线性动态系统的卡尔曼滤波问题,可以用状态方程和测量方程来描述,前者以状态矢量刻画系统的动态,后者描述系统中的测量误差。 假设研究离散线性动态系统的N维参数的状态矢量为,M维观察数据的测量矢量为,通常矢量和都是随机变量,由他们表示系统模型的状态方程和测量分别为 (1) (2) 其中为系统在和时刻的状态转移矩阵,为已知的测量矩阵 系统动态噪声和观察噪声的统计特性为 , (3) , (4) (5) “H”表示共轭转置;当n=k,=1,当n≠k,=0;噪声矢量和统计独立的。根据观察数据的测量矢量…,可求出系统状态的线性无偏最小方差估计。当时,这种最佳估计问题成为卡尔曼滤波;当时,则称为最优预测;两者之间存在密切的关系。 29、设有两个线性时不变系统如图所示,它们的频率响应函数分别为和。若两个系统输入同一个均值为零的平稳过程,它们的输出分别为、。问如何设计和才能使、互不相关。 解答: 其中,上式表明与的互相关函数只是时间函数的函数。由 故当设计两个系统的频率响应函数的振幅频率特性没有重叠时,则=0,从而有=0=,即与互不相关。 30、什么叫白噪声? 答: 白噪声是指功率谱密度在整个频域内均匀分布的噪声。 所有频率具有相同能量的随机噪声称为白噪声。理想的白噪声具有无限带宽,因而其能量是无限大 31、设连续时间随机信号以等时间间隔取样后得到离散时间随机序列。为使是白色随机序列,讨论应满足什么条件? 解:设的自相关函数和功率谱分别为和,的自相关序列和功率谱分别为 和。 根据定义,有 ① ② 按照定义,自相关序列为 因此就是以周期对的等间隔取样。 根据上式和的定义,有 == ③ 另一方面,根据S(e)的定义,有 = ④ 将式③表示成无限个积分之和,其中每个积分都在长为的区间上进行,即= 变量置换:,即得 = 交换求和与积分的次序,考虑到对于所有整数和有,因此= 再将代入式中,得 = ⑤ 比照式⑤与式④,得 = ⑥ 或 = ⑦ 如果是白色随机序列,则它的自相关序列应当是一个幅度为的冲激序列,即 上式代入式②,得 ⑧ 由式⑥、⑧有 , 上式表明,若要是白色随机序列,则要求 =常数。 32、用雷达测量地球和月球之间的距离,测量过程用下列方程描述 其中是均值为零,方差为的白噪声序列,它表示测量误差。为了提高测量精度,现采用以下两种滤波器分别对进行处理,试比较其方差的大小。 滤波器1 滤波器2 式中。 解答:两个滤波器的系统函数分别为 因此两滤波器的直流增益,无直流失真。两滤波器输出噪声平均功率为 因此 33、简述经典功率谱估计与现代功率谱估计的差别。 信号的频谱分析是研究信号特性的重要手段之一,通常是求其功率谱来进行频谱分析。功率谱反映了随机信号各频率成份功率能量的分布情况,可以揭示信号中隐含的周期性及靠得很近的谱峰等有用信息,在许多领域都发挥了重要作用。然而,实际应用中的平稳随机信号通常是有限长的,只能根据有限长信号估计原信号的真实功率谱,这就是功率谱估计。 功率谱估计分为经典谱估计和现代谱估计。经典谱估计是将数据工作区外的未知数据假设为零,相当于数据加窗,主要方法有相关法和周期图法;现代谱估计是通过观测数据估计参数模型再按照求参数模型输出功率的方法估计信号功率谱,主要是针对经典谱估计的分辨率低和方差性能不好等问题提出的。 34、两个联合平稳信号和的互相关函数为: 其中为单位阶跃函数。求互功率谱密度和。 解:直接查傅氏变换表,得 利用互谱密度的性质有 = = 35、观测信号为,其中有信号是为恒量平稳序列,其统计特征已求得为 噪声是零均值白噪声,且与有用信号不相关,即 求维纳滤波器? 解:观测的自相关函数为 观测有与有用信号之间的互相关函数为: ] 则维纳-霍甫方程式为: 由此得维纳滤波器为: 故滤波输出为: 36、平稳信号形成的滤波器为 其中为方差为0.64的零均值白噪声。观测信号为 其中是零均值单位(即方差为1)白噪声。求维纳滤波器(其脉冲响应序列的项数指定为2)。 解:由的形成滤波器得其频率特性为 的功率谱密度为: 上式的收敛域必是包括平面单位圆在内的环域,因为只有这样谱密度才存在。则的自相关函数为上式的傅氏反变换即双边反变换: 容易验证,符合上述收敛域的原函数为 其中右边序列的变换成。的自相关函数为 另外,与的互相关函数为 当维纳滤波器的脉冲响应序列为2项时,维纳-霍夫方程式为 把算得的相关函数代入上式,得 这就是维纳滤波器的脉冲响应函数。 38、设观测量由有用信号和与不相关的零均值白噪声相加而成,即 且已估计出它们的相关函数分别为 (…,-1,0,1…) 求非因果维纳滤波器的频率特性。 解: 又有 故有 最后可得,非因果维纳滤波器的频率特性为 41、令x(t)是一个是不变的标量随机变量,它在加性高斯白噪声v(t)中被观测,即y(t)=x(t)+v(t)为观测数据,若用kalman滤波器自适应估计x(t),试设计kalman滤波器。 (1) 构造离散时间的状态空间方程 (2) 求出状态变量想x(k)的更新公式。 解:x(t)是一个时不变的随机变量,故x(t)关于时间t的一阶导数等于0,即有=0 这就是连续时间的状态方程。观测方程为y(t)=x(t)+v(t) 令x(t)是一个具有均值、方差的随机变量,记做x~(,);加性观测噪声v(t)的均值为0,方差为,记做v(t)~(0,)。现在,用T=1作为采样间隔,对x(t)和v(t)等离散化,则离散时间的状态空间模型: x(n+1)=x(n) y(n)=x(n)+v(n) 式中,加性观测噪声v(n)~N(0,). 上述状态空间模型的kalman滤波算法如下: g(n)= (n+1)=(n)+g(n)[y(n)-(n)] K(n+1,n)=K(n,n-1)[1-g(n)]=g(n) g(n)的一般表达式: 由于K(1,0)=E{|x(1)—E{x(1)}|}=故当n=1时有 g(1)= K(2,1)=g(1) 当n=2时,则有g(2)== K(3,2)=g(2)= n=3时,有 g(3)== K(4,3)=g(3)= 若令 g(n-1)= K(n,n-1)=g(n-1)= 则 g(n)== k(n+1,n)=g(n)= 42、卡尔曼滤波和维纳滤波的关系及存在的问题。 答:卡尔曼滤波有一个过渡过程,而在稳态下与维纳滤波有相同的结果,是因为它们都是以最小均方误差为准则的线性估计器。卡尔曼滤波与维纳滤波中解决最佳滤波的方法也不相同。维纳滤波是用频域及传递函数的方法,而卡尔曼滤波是用时域及状态变量的方法,在理论上是维纳滤波的推广和发展,特别是在处理多变量系统,时变线性系统及非线性系统的最佳滤波等领域,提供了一种比较有效的方法,克服了基于频域处理所遇到的困难。困难包括:维纳滤波要求平稳,而卡尔曼滤波则不要求;他容许初始时间不是负无穷大,这在很多情况下是有实际意义的;卡尔曼滤波的另一个不同点是把状态或信号过程的产生看成是白噪声激励有限维数系统的输出;此外,维纳滤波要求过程的自相关函数和互关函数的简单(先验)知识,而卡尔曼滤波则要求时域中状态变量及信号产生过程的详细知识。卡尔曼滤波在时域上采用线性递推形式对观测值进行处理,能实时地给出系统状态的最优估计,并突破了单维输入和输出的限制。卡尔曼滤波算法的这些优点使它在信号和信息系统中得到了比较广泛的应用。卡尔曼滤波算法在具体应用中也存在一些实际问题,包括:   (1)模型误差和数值发散。 即使能够获得精确的模型,也常会因精确模型太复杂,维数过高而与实时处理必须减少计算量及尽量简化模型的要求相矛盾。近似或化简的模型与精确模型之间存在误差,模型误差必然会给滤波带来影响,严重时还会造成滤波结果不收敛。 2)实时要求。   影响卡尔曼滤波算法的实时性主要是状态维数n和增益矩阵的计算,它们往往有很大的计算量。一般在计算中采取某些措施,例如应用定常系统新算法或在精度损失允许的情况下尽量减少维数等到措施,从而减少计算量以满足实时滤波的要求。 43、卡尔曼滤波的特点 卡尔曼滤波具有以下的特点: 答:(1) 算法是递推的状态空间法采用在时域内设计滤波器的方法,因而适用于多维随机过程的估计;离散型卡尔曼算法适用于计算机处理。 (2) 用递推法计算,不需要知道全部过去的值,用状态方程描述状态变量的动态变化规律,因此信号可以是平稳的,也可以是非平稳的, 即卡尔曼滤波适用于非平稳过程。 (3) 卡尔曼滤波采取的误差准则仍为估计误差的均方值最小。 44、令,……是一个具有概率密度函数的正态分布得到的随机观测样本,试确定均值和方差的最大似然估计。 解:似然函数是均值和方差二者的函数,故有 从而有 分别求关于和的偏导,然后令偏导为零,得到 从可以解出 将其代入可解得 由样本均值和样本方差是无偏的,因此均值最大似然估计为无偏估计,而方差的最大似然估计则是有偏的。 45、Burg递推较levenson递推法有什么优势,并写出Burg递推法求解AR模型参数的递推公式。 答:(1)列文森(levenson)递推法需要先由信号的观测数据估计自相关函数,这是它的缺点,而伯格(Burg)递推法则由信号观测数据直接计算AR模型的参数,Burg递推法利用levenson递推公式,导出前向预测误差和后向预测误差,并按照它们最小的原则求出,从而避免求自相关函数这一难题。 (1) Burg递推法求AR模型参数的递推公式如下: ① ② ③ ④ ⑤ ⑥ ⑦ ⑧ ⑨ ⑩ 46、简述正交性原理信义,维纳滤波原理及其结论。 解:(1)正交性原理: 文字表述:使代价函数最小化的充分必要条件是估计误差与输入正交。 (1) 定义输入向量 则其相关矩阵为 式中使用了自相关函数的性质,类似地,输入与期望响应的互相关向量为 综上式子可以将winer-Hopf方程组写成紧凑的矩阵形式 ① 式中表示横向滤波器的最优抽头权向量: 由矩阵方程①,立即可以得到最优抽头权向量的解为 满足这个关系的离散时间横向滤波器称为维纳滤波器。 47、关于维纳滤波器的两个主要结论: ①维纳滤波器最优抽头权向量的计算需要已经以下统计量:(1)输入向量的自相关矩阵;(2)输入向量与期望响应的互相关向量。 ②维纳滤波器实际上是无约束优化最优滤波问题的解。 48、已经信号的四个观察数据为分别用自相关法和协方差法估计AR(1)模型参数。 解:自相关法: 协方差法: 49、假定是一个满足差分方程式 的AR()过程,且该过程是在一与独立的加性观测白噪声中观测的,即,其中的方差为,求的功率谱。 解:由已知差分方程式可得的谱密度 当与互相独立时, 故的功率谱 所以 50、分别解释“滤波”和“预测”。 解:用当前的和过去的观测值来估计当前的信号y(n)=(n)称为滤波;用过去的观测值来估计当前的或将来的信号 ,N≥0,称为预测。 51、介绍维纳滤波和卡尔曼滤波解决问题的方法。 解:维纳滤波是根据全部过去观测值和当前观测值来估计信号的当前值,因此它的解形式是系统的传递函数H(Z)或单位脉冲响应h(n);卡尔曼滤波是用当前一个估计值和最近一个观测值来估计信号的当前值,它的解形式是状态变量值。 维纳滤波只适用于平稳随机过程,卡尔曼滤波就没有这个限制。设计维纳滤波器要求已知信号与噪声的相关函数,设计卡尔曼滤波要求已知状态方程和量测方程。 h(n) 52、已知下图中x(n)=s(n)+w(n),且与统计独立,其中的自相关序列为=,w(n)是方差为1的单位白噪声,试设计一个N=2的维纳滤波器来估计s(n),并求最小均方误差。 x(n)=s(n)+w(n) 维纳滤波的输入和输出的关系 解:依题意,已知信号的自相关和噪声的自相关为:=,=δ(w),代入,得 解得:h(0)=0.451,h(1)=0.165 将上式结果代入式E=- ,求得最小均方误差: E=-(m)=1-h(0)-0.6h(1)=0.45 1. 某独立观测序列其均值为m,方差为。现有两种估计算法: 算法A:均值估计为,算法B:均值估计为 请对这两种估计算法的无偏性和有效性进行讨论。(12分) 答:算法A:均值估计为,则 ,, 均值估计是无偏估计 1.设是离散时间平稳随机过程,证明其功率谱。 证明:将通过冲激响应为的LTI离散时间系统,设其频率响应为 输出随机过程的功率谱为 输出随机过程的平均功率为 当频率宽度时,上式可表示为 由于频率是任意的,所以有 5、假定输入信号{x(t)}是一个零均值的高斯白噪声,其功率谱为,且线性系统的冲激响应为 求输出y(t)=x(t)*h(t)的功率谱及协方差函数。 解:由题知,系统的传递函数为 有此得 由输出功率谱与输入功率谱、系统函数之间的关系,得 输出的协方差函数为功率谱的傅里叶反变换,故有 6、BT谱估计的理论根据是什么?请写出此方法的具体步骤。 答:(1)相关图法又称BT法,BT谱估计的理论根据是:通过改善对相关函数的估计方法,来对周期图进行平滑处理以改善周期图谱估计的方差性能。 (2)此方法的具体步骤是: ①给出观察序列,估计出自相关函数: ②对自相关函数在(-M,M)内作Fourier变换,得到功率谱: 式中,一般取,为一个窗函数,通常可取矩形窗。 可见,该窗函数的选择会影响到谱估计的分辨率。 7、对于连续时间信号和离散时间信号,试写出相应的维纳-辛欣定理的主要内容。 答:(1)连续时间信号相应的维纳-辛欣定理主要内容: 连续时间信号的功率谱密度与其自相关函数满足如下关系: (2)离散时间信号相应的维纳-辛欣定理主要内容: 离散时间信号的功率谱密度与其自相关函数满足如下关系: 12、AR谱估计的基本原理是什么?与经典谱估计方法相比,其有什么特点? 答:(1)AR谱估计的基本原理是:阶的AR模型表示为: 其自相关函数满足以下YW方程: 取,可得到如下矩阵方程: 在实际计算中,已知长度为N的序列,可以估计其自相关函数,再利用以上矩阵方程,直接求出参数及,于是可求出的功率谱的估计值。 13、已知信号模型为s(n)=s(n-1)+w(n),测量模型为x(n)=s(n)+v(n),这里w(n)和v(n)都是均值为零的白噪声,其方差分别为0.5和1,v(n)与s(n)和w(n)都不相关。现设计一因果IIR维纳滤波器处理x(n),以得到对s(n)的最佳估计。求该滤波器的传输函数和差分方程。 解:根据信号模型和测量模型方程可看出下列参数值:a=1,c=1,Q=0.5,R=1。将它们代入Ricatti方程Q=P-a2RP/(R+c2P) 得 0.5=P-P/(1+P) 解此方程得P=1或P=-0.5,取正解P=1。 再计算维纳增益G和参数f:G=cp/(R+c2P) =1/ (1+1) =0.5 f=Ra/(R+c2P) =1/ (1+1) =0.5 故得因果IIR维纳滤波器的传输函数和差分方程分别如下: Hc(z)=G/(1-fz-1)=0.5/(1-0.5z-1) (n)=0.5(n-1)+0.5x(n) 14、简述AR模型功率谱估计步骤。 步骤1:根据N点的观测数据uN(n)估计自相关函数,得,m=0,1,2,…,p, 即 步骤2:用p+1个自相关函数的估计值,通过直接矩阵求逆或者按阶数递推的方法(如Levinson-Durbin算法),求解Yule-Walker方程式,得到p阶AR模型参数的估计值 和 步骤3:将上述参数代入AR(p)的功率表达式中,得到功率谱估计 ,即 3、已知输入信号向量u(n)的相关矩阵及数学期望响应信号d(n)的互相关向量分别为 且已知期望相应d(n)的平均功率为E{d2(n)}=30。 (1)计算维纳滤波器的权向量。(2)计算误差性能面的表达式和最小均方误差。 解:(1)根据维纳霍夫方程Rω0=p 得 ω0=R-1p (2)误差性能面的表达式为J(ω)=σ2d-pHω-ωHp+ωHRω 最小均方误差值为将ω0代入上面的误差性能面表达式得 Jmin=σ2d-pHω-ωHp+ωHRω =σ2d-pHω0 =30-14=16。 5、在测试某正弦信号的过程中叠加有白噪声v(n),即测试结果为 设计一个长为N=4的有限冲激响应滤波器,对x(n)进行滤波后得到s’(n),它与s(n)的误差的均方值最小,求该滤波器的冲激响应和估计误差的平均功率。 解答:已知 设 取 7、设信号的自相关序列为 观测信号为式中是方差为0.45的零均值白噪声,它与相互统计独立。试设计一个长为N=3的FIR滤波器来处理,使其输出与的差的均方值最小。 解答: 已知所以 10.用观测数据(y(n),y(n-1))自适应估计随机变量x(n).已知Ryy=[1 0.4;0.4 1],为保证收敛,µ的值应限制在什么范围?若Ryy=[1 0.8;0.8 1],则自适应滤波器的收敛速度将会更快还是更慢? 解: 11、已知满足AR(2)模型,即满足如下差分方程: 其中是均值为零、方差为的白噪声。试用自相关函数来表示系数、。 答案:AR模型的正则方程式可以表示为 和 将带入上面两式为: 和 可以解得 30、白噪声是现代信号处理中常用的一种随机信号,请从时域和频域两个角度对其加以阐述。 答:设为实值平稳过程,若它的均值为零, 在时域中,其自相关函数仅在0点有一个冲击函数,在其他点均为0; 在频域中,谱函数在所有频率范围内为非零的常数,则称X(t)为白噪声过程。 33、 一个差分滤波器的输出为:y(n)=x(n)+x(n-1), n=1,2…… 令x(n)的功率谱为1/(1+f2),试求差分滤波器输出y(n)的功率谱密度。 解:由题意知,系统的传递函数H(z)=1+z-1,故 所以系统输出y(t)的功率谱密度为: 39、已知输入信号向量的相关矩阵及与期望响应信号的互相关向量分别为 ,且已知期望响应的平均功率为。 (1) 计算维纳滤波器的权向量。 (2) 计算误差性能面的表达式和最小均方误差。 解: 40、离散时间的二阶AR过程由差分方程 描述,式中是一零均值、方差为的白噪声。证明的功率谱为 解: 由AR过程的功率谱公式知(式1) 式中 = = =(式2) 将(式2)代入(式1)中可得: 证毕。 45、设随机序列,其中是两两互不相关的随机变量且 ,序列被称作白噪声。验证白噪声序列是平稳序列。 解:显然均值函数为常数,当时,因为不相关,所以 当 所以,只是时间差的函数,序列是平稳的 46、若序列x(n)为实因果序列,h(0)=1,其傅氏变换的虚部为H1(ejω)=—sinω,求序列h(n)及其傅氏变换H(ejω)。 解:因为H1(ejω)=—sinω=—( ejω—e-jω)=0(n) e-jω h0(n)= =-δ(n+1)+δ(n-1), h(n)= = 所以 h(n)=δ(n)+δ(n-1); H(ejω)=1+ e-jω 49、简述AR模型功率谱估计的方法 答:(1)根据N点的观测数据估计自相关函数,得,即 (2)用个自相关函数的估计值,通过直接矩阵求逆或者按阶数递推的方法(如Levinson-Durbin算法),求解Yule-Walker方程 得到p阶AR模型的参数估计值和。 (3)将上述参数带入AR(p)的功率谱表达式中,得到功率谱估计式,即 50、简述LMS算法 答:(1)初始化, 权向量: 估计误差: 输入向量: (2)对 权向量的更新: 期望信号的估计: 估计误差: (3)令,转到(2) 63、请写出ARMA的数学模型表达式,并画出该模型的电路框图。 答:(1)ARMA的数学模型表达式: 式中,为常数, (2)该模型的电路框图如下所示: 100、AR谱估计中的虚假谱峰是怎样产生的?怎样避免产生虚假谱峰?谱线分裂的原因是什么?相位怎样影响谱线的位置?如何减少这种影响?怎样避免谱线分裂?噪声对AR谱估计有什么影响?怎样减少这种影响? 答:1.如果自相关函数的取样值的估计没误差,模型参数的估计在理论上应该为:,但是实际上自相关函数的估计是有误差的,这使得对大于P的值有,相应的也会产生n-p个额外的极点在单位圆附近,这样就会形成虚假谱峰。2.要避免产生虚假谱峰,要求模型的技术不宜过高,最高不超过N/2,其中N是观测数据长度。3.如果估计的随机过程是由一个正弦信号叠加噪声构成的,那么对某些算法会观察到AR谱估计中存在两个靠的很近的谱峰,似乎在随机过程中还存在着另一个正弦信号,这种现象叫谱线分裂。原因有高信噪比,正弦分量初相位为π/4的奇数倍,数据序列的长度为正弦分量的1/4周期的奇数倍,估计的AR参数数目与数据的个数相比占有较大的百分比,虚假谱峰常伴随着谱线分裂现象,这与观测数据长度太短有关。4.AR谱估计中谱峰出现的位置与正弦信号的初相位有着很密切的关系,谱峰位置对相位的依赖性随着观测数据长度的增加而减少,一般认为谱峰对相位的依赖是由正弦信号的正负频率成分之间的相互作用引起的。5.解决的方法是用解析信号代替真实信号,对解析信号进行欠抽样,并利用复数据进行AR谱的估计。这样就不要求有复共轭极点,而模型的阶数可减少一半。另一种方法是对Burg算法的反射系数估计公式进行修改,用一实值窗函数加权。6.通过增加观测数据的长度,或者利用向前和向后最小二乘算法和Marple递推算法。7.AR谱估计方法对观测噪声比较敏感,噪声会使谱峰展宽,从而导致分辨率下降,而且会使谱峰偏离正确的位置。8.减少噪声对AR谱估计的影响一般使用四种方法:一、采用ARMA谱估计方法;二、对观测数据进行滤波减小噪声;三、采用高阶AR模型;四、补偿自相关函数或者反射系数估计中噪声的影响。 104、设某积分电路输入输出之间满足以下关系,式中,T为积分时间。并设输入输出都是平稳过程。求证输出功率谱密度为 (提示:,而,是矩形方波。) 证明:因为 所以 所以 所以 所以 111、Wiener滤波器是现代信号滤波处理的经典,其核心在于考察滤波器输入输出信号之间的关系,请用恰当的数学模型对其加以描述。 滤波器的理想输出为s(t+a) 估计误差为e(t)=s(t+a)-y(t) 估计误差的平方为: 而 代入上式,两边取数学期望,得到均方误差: 其中,Rs s(t)的自相关函数 Rx x(t)=s(t)+n(t)的自相关函数 Rs,x s(t)和x(t)之间的互相关函数 若信号s(t)和噪声n(t)不相关,且噪声均值为零,即E[n(t)]=0,则有: 维纳滤波就是希望求出最优h(u),使得最小。 本文档由香当网(https://www.xiangdang.net)用户上传

    下载文档到电脑,查找使用更方便

    文档的实际排版效果,会与网站的显示效果略有不同!!

    需要 2 香币 [ 分享文档获得香币 ]

    下载文档

    相关文档

    西南交大现代信号处理部分答案

    西南交大现代信号处理部分答案题1:(1) 是随变化的随机信号,因此=.所以谐波信号的均值为 =由于谐波信号的均值等于零,故其方差等于二阶矩,既有 所以x(t)的方差为 谐波信号的自相关函数又所...

    8个月前   
    103    0

    测试技术与信号处理复习题(客观题)

    第 1、3 章一、判断题.1. 对多次测量的数据取算数平均值,就可以减少随机误差的影响。 ( √ )2. 传感器的线性范围越宽,表明其工作量程越大。 ( √ )3. 一台仪器的重复性很好,但测...

    2年前   
    864    0

    信号处理上机

    XX大学硕士研究生 2014/2015 学年上机作业考试科目: 近代信号处理 课程编号: 姓名: ...

    2年前   
    210    0

    数字信号处理复习总结

    信号是信息的载体,以某种函数的形式传递信息。这个函数可以是时间域、频率域或其它域,但最基础的域是时域。

    2年前   
    400    0

    《数字信号处理》习题集(附答案)

    第一章 数字信号处理概述简答题:1. 在A/D变换之前和D/A变换之后都要让信号通过一个低通滤波器,它们分别起什么作用?答:在A/D变化之前让信号通过一个低通滤波器,是为了限制信号的最高频率,...

    8个月前   
    95    0

    现代公文写作与处理:通报

    现代公文写作与处理:通报  一、通报的含义   通报是党政机关、企事业单位和社会团体把工作情况、经验教训、典型事例及有教育、指导、规范意义的事件向下级或公众发布的周知性公文。适用于表彰先进,批...

    8年前   
    453    0

    2017年现代公文写作与处理:公文处理

    现代公文写作与处理:公文处理  公文的处理  一、公文处理的概念   公文处理是包括公文拟制、办理、管理、立卷、归档在内的一系列衔接有序的工作。  二、公文处理的作用   公文处理具有十分重要...

    6年前   
    298    0

    现代公文写作与处理:讲话稿

    现代公文写作与处理:讲话稿  一、 讲话稿的概念   讲话稿也称“发言稿”,是讲话人在各种会议、集会、仪式上发表讲话而准备的发言文稿。它能帮助讲话人准确表达意思,避免讲话的随意性。   二、话...

    11年前   
    344    0

    现代公文的写作与处理:公文概述

    现代公文的写作与处理:公文概述  作者:陈华平  王太钧 孙杰  公文的涵义与特点   一、公文的涵义   公文是“公务文书”的简称。   从广义上讲,公文是指党政机关、社会团体、企事业单位为...

    9年前   
    743    0

    现代控制理论复习题库

    现代控制理论复习题库一、选择题1. 下面关于建模和模型说法错误的是( C )。A.无论是何种系统,其模型均可用来提示规律或因果关系。B.建模实际上是通过数据、图表、数学表达式、程序、逻...

    2个月前   
    0    0

    数字信号的处理与滤波设计课程设计

    — — 装 订 线 — —数字信号的处理与滤波 摘要: 数字信号处理是利用计算机或专用处理设备,以数值计算的方法,对信号进行采集、滤波、增强、压缩、估值和识别等加工处理,...

    2年前   
    221    0

    《数字信号处理》期末试题库(有答案)

    。一. 填空题1、一线性时不变系统,输入为 x (n)时,输出为 y(n) ;则输入为 2x(n)时,输出为 2y(n) ;输入为 x(n-3 )时,输出为 y(n-3) 。2、从奈奎斯特采...

    1年前   
    514    0

    《数字信号处理》第三版课后习题答案

    数字信号处理课后答案 1.2 教材第一章习题解答 1. 用单位脉冲序列及其加权和表示题1图所示的序列。 解: 2. 给定信号: (1)画出序列的波形,标上各序列的值; (2...

    4年前   
    891    0

    信号电缆、电源接地处理技术总结报告

    信号电缆、电源接地处理技术总结报告 2013年,乌西车间党总支在段、段党委的领导下,以制约车间安全生产的薄弱环节为切入点,全力开展党员攻关活动,处理了大量电源电缆接地隐患,取得了良好的效果。...

    9年前   
    8639    0

    噪声中正弦信号的现代法频谱分析

    实验报告一、实验名称噪声中正弦信号的现代法频谱分析二、实验目的通过对噪声中正弦信号的现代法频谱分析,来理解和掌握现代谱估计的基本概念,以及学会应用现代谱估计以及改进后的方法。三、基本原理1.参...

    1个月前   
    0    0

    2017年现代公文写作与处理:调查报告

    现代公文写作与处理:调查报告  调查报告概述  一、调查报告的含义  调查报告是根据某种特定的需要,有计划地对典型事物、社会问题或工作情况等,进行认真调查研究后形成的书面报告形式的公文。  二...

    6年前   
    194    0

    现代公文写作与处理:公文要素之公文结构

    现代公文写作与处理:公文要素之公文结构  公文结构的概念和作用   (一)公文结构的概念   公文结构是指公文内容的组合和构造。换句话说,也就是如何按照主题的需要,安排材料、组织成文的方式,比...

    8年前   
    711    0

    现代公文的写作与处理:公文要素之公文主题

    现代公文的写作与处理:公文要素之公文主题  公文主题  公文主题,亦称之为公文主旨,是公文的首要和核心的要素。它是公文的灵魂,在文章中起到统率作用,它决定和制约材料的取舍、结构的安排、语言的组...

    8年前   
    757    0

    2017年现代公文写作与处理:述职报告

    现代公文写作与处理:述职报告  述职报告概述  一、述职报告的含义   述职报告是领导干部向选举或任命机构、上级领导机关、主管部门及本单位的群众汇报自己在一定时期内履行职务责任的书面报报告,是...

    5年前   
    265    0

    现代公文写作与处理:公文要素之公文语言

    现代公文写作与处理:公文要素之公文语言  一、公文语言的特点   (一)准确性   要使公文语言做到准确无误,要注意以下几点:   1 .概念准确,避免歧义。   2 .措词得体,防止误用。 ...

    10年前   
    890    0

    文档贡献者

    文***品

    贡献于2022-12-26

    下载需要 2 香币 [香币充值 ]
    亲,您也可以通过 分享原创文档 来获得香币奖励!
    下载文档

    该用户的其他文档