现代信号处理大作业
现代信号处理作业

信号时频分析技术及matlab仿真电路与系统王冠军 201128013926153摘要:本文介绍了时频分析的一些基础理论,对短时傅里叶变换Wigner-Ville分布做了简单介绍,运用MATLAB语言实现了旨在构造一种时间和频率的密度函数,以揭示信号中所包含的频率分量及其演化特性的wigner-ville分布。
并对时频分析方法的优缺点进行了分析。
关键词:时频分析短时傅里叶变换wigner-ville分布1 引言基于Fourier变换的传统信号处理技术从信号频域表示及能量的频域分布的角度揭示了信号在频域的特征。
但Fourier变换是一种整体变换,只能为人们提供信号在时域或频域的全局特性而无法了解信号频谱随时间变化的情况。
因此,需要使用一种时间和频率的联合函数来表示信号,这种表示简称为信号,也就是信号的时频分析。
2 时频分析方法信号时频分析主要研究非平稳信号或时变信号的频谱含量是怎样随时间变化的。
时频分析是当今信号处理领域的一个主要研究热点,目前常用时频分析方法主要有短时傅里叶变换、Gabor展开、小波变换、Wigner-Ville分布。
本文主要介绍了短时傅里叶变换和Wigner-Ville分布两种分析方法。
2.1 短时傅立叶变换STFT从历史上看,信号的时频分析用的最多的是短时傅立叶变换,这种变换的基本思想是用一个窗函数乘时间信号,该窗函数的时宽足够窄,使取出的信号可以被看成是平稳的,然后进行的傅立叶变换可以反映该时宽中的频谱变化规律,如果让窗随时间轴移动,可以得到信号频谱随时间变化的规律。
对于时变信号,了解不同时刻附近的频域特征是至关重要的。
因此,人们采用时间—频率描述时变信号,将一维的时域信号映射到一个二维的时域平面,全面反映观测信号的时频联合特征。
短时傅立叶变换反映了这一思想,对于时变信号,采用某一滑动窗函数截取信号,并认为这些信号是准平稳的,然后,再分别对其进行傅立叶变换,构成时变信号的时变谱。
现代信号处理大作业王成志1

《现代信号处理》大作业姓名:王成志学号:1140349078一. L D 迭代算法的matlab 实现1.1 Levinson-Durbin 算法介绍功率谱估计大致可以分为经典谱估计和现代功率谱估计,经典谱估计方法存在着以下三点缺陷:(1)数据加窗或自相关加窗,都隐含着假定在窗外未观测到的数据或自相关系数为零,该假设不切实际。
(2)要性能好往往需要较长的数据,但实际数据长度有限(3)窗函数容易造成谱的模糊。
采用AR 模型的现代谱估计方法可以克服这些不足。
其中LD 递推算法可以在计算机上方便实现。
LD 递推算法具体计算步骤如下:(1) Yule-Walker 方程的矩阵形式(1)所示:⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎣⎡=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡--+-----001)0()2()1()()1()1()0()1()()2()1()0(2,1,σk k k xx xx xx xx x xx xx xx xx xx xx xx a a r k r k r k r k r r r r k r r r r 系数矩阵xx Hxx R R =,为Hermitian 矩阵,对角线上元素相同,即为Topliez 矩阵。
(2) P-1阶Yule-Walker 方程为:21111(0)(1)(1),1(1)(0)(2)0,1(1)(2)(0)0x x x p p x x x p x x x R R R p a R R R p a p R p R p R σ-----⎡⎤⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥-⎢⎥⎢⎥⎢⎥=⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥---⎢⎥⎣⎦⎣⎦⎣⎦ 其中,2211{()}p p E e l σ--=为误差功率。
写成联立方程:2111,0,0()0,1,,1p pp k xk m a R m k m p σ---=⎧=-=⎨=-⎩∑ 取共轭得:21**11,0,0()0,1,,1p pp kxk m aR m k m p σ---=⎧=-=⎨=-⎩∑变量替换,并利用*()()x x R l R l =得:21*11,10,1()0,0,,2p pp p kx k m p aR m k m p σ-----=⎧=--=⎨=-⎩∑ 表示成矩阵:*1*1210(0)(1)(1),10(1)(0)(2),2(1)(2)(0)1x x x p x x x p p x x x R R R p a p R R R p a p R p R p R σ-----⎡⎤⎡⎤⎡⎤-⎢⎥⎢⎥⎢⎥--⎢⎥⎢⎥⎢⎥=⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥--⎢⎥⎣⎦⎣⎦⎣⎦ 求解得:*.1,1,,0,,p k p k p p p k a a K a k p ---=+=22*1p p p p K σσ-=+∆ 2210p p p K σ-=∆+,p p p K a =222*22111[][1]p p p p p p p K K K σσσσ---=+-=-(3) 当k=1时,即一阶递推为:⎥⎦⎤⎢⎣⎡=⎥⎦⎤⎢⎣⎡⎥⎦⎤⎢⎣⎡-01)0()1()1()0(211,1σa R R R R x x x x求解可得:)1()0()0()1( ,11,1211,10,1x x x x R a R R R a a +=-==σ(4) 对于2≥p 时,递推为:10,≡p a , *,1,1,k p p p k p k p aK a a ---+=, ]1[2212p p p K -=-σσ 21,-∆-==p pp p p a K σ∑-=--+=∆11,1)()(p k x kp x p k p R ap R矩阵R x 已知,可得到各阶AR 模型系数为:)0())1(1( ,)0()1()1(2111xx xx xx r a r r a -=-=ρ11111)()()()(--=--∑-+-=∆-=k k l xx k xx k kk l k r l a k r k a ρρ1,,2,1)()()()(*11-=-+=--k i i k a k a i a i a k k k k12))(1(--=k k k k a ρρ1.2实验结果(1) 输入p=3,rr = [70,60,50,40] 时,求得AR 模型估计参数为:a =1.0000 -0.8571 0 0 1.0000 -0.5275 -0.3846 0 1.0000 -0.7572 -0.6996 0.5972 各阶求得的方差为:sigma = 18.5714 15.8242 10.18013阶时,a 3 (1)= -0.7572 a 3 (2)= -0.6996 a 3 (3)= -0.5972(2) 输入p=5,rr = [30,45,26,33,47,43]时,AR 模型估计参数为:a =1.0000 -1.5000 0 0 0 0 1.0000 0.2800 -1.1867 0 0 0 1.0000 0.8227 -1.3147 -0.4573 0 0 1.0000 1.9708 1.9858 -2.5226 -2.5105 0 1.0000 1.0869 1.0977 -1.8235 -1.8166 0.3521 各阶求得的方差为: sigma =37.5000 15.3067 12.1054 64.1881 56.23165阶时, a 5 (1)= 1.0869 a 5(2)= 1.0977 a 5(3)= -1.8235 a 5(4)= -1.8166 a 5(5)= 0.3521二. 一维平稳信号由两个高斯信号叠加而成12241122()()[exp(())exp(())]22z t t t j t t t j t αααωωπ=--++--+,其中12,t t >12ωω>,分别求出()z t 的WV 分布及其模糊函数,画出二者的波形图,指出并分析其信号项和交叉项。
现代信号处理大作业题目+答案

研究生“现代信号处理”课程大型作业(以下四个题目任选三题做)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]. 北京: 电子工业出版社, 2006[4] S.Haykin, 郑宝玉等译. 自适应滤波器原理(第四版)[M].北京: 电子工业出版社, 2003[5] J. G. Proakis, C. M. Rader, F. Y. Ling, etc. Algorithms for Statistical Signal Processing [M].Beijing: Tsinghua University Press, 2003一、请用多层感知器(MLP)神经网络误差反向传播(BP)算法实现异或问题(输入为[00;01;10;11],要求可以判别输出为0或1),并画出学习曲线。
西南交大现代信号处理作业

现代信号处理作业1.(5″)证明下面定理:任何一个无偏估计子方差的下界叫作Cramer-Rao 下界 定理:令1(,,)N x x x =为一样本向量,(|)f x θ是x 的条件密度,若ˆθ是θ的一个无偏估计子,且(|)/f x θθ∂∂存在,则221ˆˆvar()()[ln (|)]E E f x θθθθθ=-≥∂∂式中ˆln (|)()()f x K θθθθθ∂=-∂。
其中()K θ是θ的某个不包含x 的正函数。
2.(10″)Wiener 滤波是信号处理中最常用和基础的波形估计工具之一,对其在自己研究领域的应用情况进行一个简单综述。
3.(5″)二阶滑动平均过程由2()()1(1)2(2),{()~(0,)}x n w n b w n b w n w n N σ=+-+-定义,式中2(0,)N σ表示正态分布,其均值为零、方差为2σ。
求x(n)的功率谱。
4.(20″)信号的函数表达式为:()sin(2100)sin(2300)()sin(2200)()()x t t t A t t dn t n t πππ=++++,其中,A(t)为一随时间变化的随机过程,dn(t)为经过390-410Hz 带通滤波器后的高斯白噪声,n(t)为高斯白噪声,采样频率为1kHz ,采样时间为2.048s 。
(1) 利用现代信号处理知识进行信号的谱估计; (2) 利用现代信号处理知识进行信号的频率提取; (3) 分别利用Wiener 滤波和Kalman 滤波进行去噪; (4) 利用Wigner-Ville 分布分析信号的时频特征。
5.(10″)附件中表sheet1 为某地2008年4月28日凌晨12点至2008年5月4日凌晨12点的电力系统负荷数据,采样时间间隔为1小时,利用ARMA 方法预测该地5月5日的电力系统负荷,并给出预测误差(5月5日的实际负荷数据如表sheet2)。
1、定理:令1(,,)N x x x =为一样本向量,(|)f x θ是x 的条件密度,若参数估计ˆθ是真实参数θ的一个无偏估计子,且(|)/f x θθ∂∂、22(|)/f x θθ∂∂存在,则ˆθ的均方误差所能达到的下界(称为Cramer-Rao 下界)等于Fisher 信息的导数,即:221ˆˆvar()()[ln (|)]E E f x θθθθθ=-≥∂∂ (1-1)不等式中等号成立的充分必要条件是:ˆln (|)()()f x K θθθθθ∂=-∂ (1-2) 其中()K θ是θ的某个正函数,与样本1(,,)N x x x =无关。
现代信号及处理作业

现代信号及其应用雷达信号处理技术发展综述姓名:******班级****** 学号******日期:2015/4/10摘要本文简单介绍了雷达发展、雷达工作原理、雷达分类、航迹、航迹处理、滤波器、滤波器工作原理以及滤波器在航迹处理中的应用。
ABSTRACTThis paper simply introduces the development of radar, the working principle of radar, the classification of radar,the track, the tracking processing, application of filter, the filter and filter in the working principle of the track processing.目录1前言 (4)2雷达 (4)2.1雷达意义简述 (4)2.2雷达分类 (4)2.3.1雷达原理 (4)2.3.2影响雷达监测目标的影响因素 (4)2.3.3探测范围 (4)2.3.4探测距离 (5)2.3.5发射距离 (5)2.3.6发射频率 (5)3航迹 (5)3.1航迹形成过程 (5)3.2航迹处理 (5)3.3.1滤波器 (5)3.3.2滤波器的工作原理 (6)3.3.3滤波器在航迹处理过程中的应用 (6)4总结 (7)雷达信号处理技术发展综述1 前言自第二次世界大战以来,雷达技术得到了迅猛的发展,雷达家族不断扩大。
雷达的用途早已超出了检测与定位的范畴。
特别是为适应现代战争和现代经济发展的需要研制的不同用途的雷达,更是种类繁多,五花八门,本文下面简单介绍雷达的发展史以及滤波器在航迹处理过程中的应用的部分内容。
2 雷达2.1 雷达意义简述雷达为英文Rader 一词的译音,来源于词组Radio Detection and Ranging ,意为无线电监测与定位,即雷达是利用无线电方法来发现目标并测定目标位置参数的一种电子系统。
现代信号处理试题

1、已知X a (t) 2COS (2f 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 tj tX a (j )X a (t)e j dt 2cos( o t)e j dt3、在时域对一有限长的模拟信号以 4KHZ 采样,然后对采到的N 个抽样做N 点DFT ,所得离散谱线的间距相当于模拟频率100HZ 。
某人想使频率能被看得清楚些,每50HZ 能有一根谱线,于是他用8KHZ 采样,对采到的2N 个样点做2N 点DFT 。
问:他的目的能达到吗? 答:不能,因为他忽略了数字频率和模拟频率的区别。
提高采样频率f s ,N 固然大了,数字频率(单位圆)上的样点数确实增加了,但从模拟频率谱看,样点一点也没有变得密集,这是因为数字频率2总是对应模拟频率 f s 。
2 f sf s采样频率由f s 到2 f s 增加一倍,N 也增加一倍,但模拟频率的采样间隔s s100Hz2N N 2),不能提 2N高模拟频率的分辨(e j 0t e j 0t上式中指数函数和傅里叶变换不存在, X a (j ) 2 [ ()(2)x a (t )X a (t) (t)e j t dt引入奇异函数)]函数,它的傅里叶变换可以表示成:nT)2cos( 0nT) (tnT)n2cos( 0nT),2、用微处理器对实数序列作谱分析,以下各参数:(1)x(n) 最小记录时间 (2) (3) (4) 解:( 1)Tpmin T max N min要求谱分辨率F最大取样时间 最少采样点数 在频带宽度不变的情况下将频率分辨率提高一倍的已知(2)F 50Hz1 F 1 T pmin10.02s 50 1(3)N min1 s minT PT 0.02s0.5 10 3s(4)辩率提高1倍(F 变成原来的12)T p0.04sN min~T 0.5 10 s频带宽度不变就意味着采样间隔 5OHZ ,信号最高频率1KHz,是确定N 值。
现代信号处理作业

现代信号处理作业现代信号处理课程作业1.做⼀个⽹络检索,简述现代信号处理技术的主要特征和技术特点,并阐述信号处理在实际⼯程中的应⽤情况代信号处理技术的主要特征和技术特点:1)精度⾼:在模拟系统的电路中,元器件精度要达到10-3以上已经不容易了,⽽数字系统17位字长可以达到10-5的精度,这是很平常的?例如,基于离散傅⾥叶变换的数字式频谱分析仪,其幅值精度和频率分辨率均远远⾼于模拟频谱分析仪?2) 灵活性强:数字信号处理采⽤了专⽤或通⽤的数字系统,其性能取决于运算程序和乘法器的各系数,这些均存储在数字系统中,只要改变运算程序或系数,即可改变系统的特性参数,⽐改变模拟系统⽅便得多?3) 可以实现模拟系统很难达到的指标或特性:例如:有限长单位脉冲响应数字滤波器可以实现严格的线性相位;在数字信号处理中可以将信号存储起来,⽤延迟的⽅法实现⾮因果系统,从⽽提⾼了系统的性能指标;数据压缩⽅法可以⼤⼤地减少信息传输中的信道容量?4)可以实现多维信号处理:利⽤庞⼤的存储单元,可以存储⼆维的图像信号或多维的阵列信号,实现⼆维或多维的滤波及谱分析等?信号处理在实际⼯程中的应⽤情况:数字信号处理是利⽤计算机或专⽤计算机或专⽤处理设备,以数据形式对信号进⾏采集,变换,滤波,估值,增强,压缩,识别等处理,以得到符合⼈们需要的信号形式?数字信号处理是以众多科学为理论基础的,他所涉及的范围及其⼴泛?DSP 技术应⽤到我们的⽣活的每⼀个⾓落,从军⽤到民⽤,从航空航天到⽣产⽣活,都越来越多地使⽤DSP. DSP技术在航空⽅⾯,主要⽤于雷达和声纳信号处理;在通信⽅⾯,主要⽤于移动电话,IP电话,ADSL和HFC的信号传输;在控制⽅⾯,主要⽤于电机控制,光驱和硬盘驱动器;在测试/测量⽅⾯,主要⽤于虚拟仪器,⾃动测试系统,医疗诊断等;在电⼦娱乐⽅⾯,主要⽤于⾼清晰度电视,机顶盒,家庭影院,DVD 等应⽤;还有数字相机,⽹络相机等等都应⽤了SP技术?同时,SOC芯⽚系统,⽆线应⽤,嵌⼊式DSP都是未来DSP的发展⽅向和趋势?可以说,没有DSP就没有对互联⽹的访问,也不会有多媒体,也没有⽆线通信?因此DSP仍将是整个半导体⼯业的技术驱动⼒?现在,DSP应⽤领域不断拓宽,其涵盖⾯包括宽带Internet接⼊业务,下⼀代⽆线通信系统的发展,数字消费电⼦市场,汽车电⼦市场的发展等诸多多⽅⾯?现代数字信号处理器是执⾏⾼速数字信号系统的IC电路,它恰好适合多媒体信息化社会需求,迅速发展壮⼤?如今,世界电⼦器件市场上,各种各样的DSP器件已相当丰富?⼤⼤⼩⼩封装形式的DSP器件,已⼴泛⽤于各种产品的⽣产领域,⽽且DSP的应⽤领域仍在不断的扩⼤,发展速度异常?2?简述信号的频率分析技术及其应⽤,阐述实现精细频率分析的实现⽅法?考虑到数字信号分析中,虽然提⾼信号的采样频率可以改善信号分析的频率分辨率,但是提⾼信号的采样频率通常需要付出额外的硬件代价,往往受制于可实现性与成本问题⽽难以实现?因此,就需要使⽤频谱细化技术在尽可能低的采样频率下提⾼数字信号分析的频率分辨率的措施?频谱细化的基本思路是对信号频谱中的某⼀频段进⾏局部放⼤,也即在某⼀频率附近局部增加谱线密度,实现选带频段分析?频谱细化技术在⽣产实践和科学研究中获得了⽇益⼴泛的应⽤?例如,齿轮箱的故障诊断要求准确分辨齿轮各阶啮合振动的主频和边频等,其频谱图上的频率间隔很细,但频率分布⼜较宽,为了识别谱图的细微结构,就必须对信号进⾏细化分析;直升机?坦克?巡航导弹的声⾳具有显著的⾮平稳性,为了得到准确的时延量,信号的取样不能太长,⽽FFT计算的频谱存在栅栏效应?因此必须采⽤有效的⽅法对频谱进⾏细化,这样才能保证⾜够的相关计算精度;在⽆线电通信信号和其他的实际⼯程信号的分析中,为了获取更⾼的测量精度和实时检测能⼒,需要对信号频谱进⾏细化分析,以提供有⽤信息?因此对频谱细化技术的研究受到普遍重视,也是当前信号处理技术研究中的⼀个⼗分活跃的课题?常见的经典⽅法有:复调制细化法?Chirp-Z变换?FFT+FT细化法?DFT补零法等很多⽅法?复调制细化法:⼜称为选带频率细化选带频谱分析,是20世纪70年代发展起来的?其传统的分析步骤为:移频(复调制)低通滤波器重抽样--FFT及谱分析频率成分调整,因其物理概念⾮常明确,所以⼀直沿⽤⾄今?FFT+FT细化法:该⽅法的原理本质是将连续傅⾥叶变换经过将积分化成求和?时域离散化和时域截断为有限长三个步骤变换得到时间离散?频率连续的特殊傅⾥叶变换形式?FFT+FT连续细化分析傅⾥叶变换法先⽤FFT做全景谱,再对指定的⼀个频率区间进⾏细化计算:先确定频率分辨率,再确定计算频率序列,最后⽤FT连续谱分析⽅法进⾏实部和虚部计算,合成幅值谱和相位谱? Chirp-Z变换:最早提出于1969年,CZT是⼀种在Z平⾯上沿着螺旋线轨道计算有限时宽的Z变换⽅法?基本原理是在折叠频率范围内任意选择起始频率和频率分辨率在这有限带宽⾥对样本信号进⾏Z变换这与频谱校正⽅法中的FFT + FT 连续细化分析傅⾥叶变换法的基本原理是⼀样的?3、通过⽹络检索,对弱信号检测技术进⾏调研,分析⼀下现代弱信号检测的⽅法微弱信号检测(WeakSignalDetection)是⼀门新兴的技术学科,应⽤范围遍及光?电?磁?声?热?⽣物?⼒学?地质?环保?医学?激光?材料等领域?其仪器已成为现代科学研究中不可缺少的设备?微弱信号检测技术是采⽤电⼦学?信息论?计算机及物理学的⽅法,分析噪声产⽣的原因和规律,研究被测信号的特点与相关性,检测被噪声淹没的微弱有⽤信号?微弱信号检测的⽬的是从强噪声中提取有⽤信号,或⽤⼀些新技术和新⽅法来提⾼检测系统输出信号的信躁⽐?信号处理系统的信躁⽐改善等于输⼊(⽩)躁声带宽与系统的躁声等效带宽之⽐?因此,减少系统的躁声等效宽度便可以提⾼系统的输出信躁⽐?对于信躁⽐⼩于1的被躁声淹没的信号,只要信号处理系统的躁声等效带宽做得很⼩,就可以将信号(或信号携带的信息)从躁声中提取出来,这就是通常的微弱信号检测的指导思想之⼀?现代弱信号检测的⽅法和原理窄带滤波法: 使⽤窄带滤波器,滤掉宽带躁声只让窄带宽信号通过(仅有极少量窄带躁声通过)?窄带滤波法能减少躁声对有⽤信号的影响?滤除掉通频带以外躁声,提⾼信号的信躁⽐?但是,由于⼀般滤波器的中⼼频率不稳定,不能满⾜更⾼的滤除躁声的要求?双路消躁声法:由于信号与躁声性能完全不同,信号⼀般为⼀些变化规律已知的量,⽽躁声是⼀些随机量满⾜统计规律?当随机性的躁声从两路到达加法器时,极性正好相反,经过加法器相加后把躁声消掉?只有少数强躁声才通过阀值电路⽽产⽣本底计数,根据统计规律?本底计数时间较长时为恒定值?故可以先测出它,然后从总计数中把它减得到信号计数?这种⽅法只能检测到微弱的正弦信号是否存在,⽽不能复现信号波形?同步累积法:利⽤信号的重复性,躁声的随机性,对信号进⾏重复累积(⼏次),使SNIR提⾼,但需耗费时间?锁定接收法(频域分析法) :锁定检测法是利⽤互相关原理,使输⼊待测的周期信号与频率相同的参考相关器中实现互相关,从⽽将深埋在躁声中的周期信号携带的信息检测出来?相关检测法: 相关检测技术是应⽤信号周期性和噪声随机性的特点,通过⾃相关或互相关运算,达到去除躁声检测出信号的⼀种技术?由于信号和躁声是相互独⽴的过程,根据相关函数和互相关函数的定义,信号只与信号本⾝相关与躁声不相关??取样积分法:取样积分(或信号平均)法是将待测的重复信号逐点多次取样并进⾏同步积累,从⽽达到从噪声中恢复信号波形的⽅法?取样积分也采⽤同步相关检测的原理和⽅法,实现从噪声中提取信号,但它的参考信号只在窗⼝持续期间与被测信相关,每周相关时间很短,此外它的相移也是在很慢的变化?取样积分由单点取样积分与多点取样积分两种?4.利⽤MATLAB产⽣出⼀个线性调频信号(chirp信号),采样频率=8000Hz,持续时间1s,起始频率=500Hz,终⽌频率=1300Hz,给出其时域波形图,请利⽤短时FFT分析函数对数据进⾏时间-频率分析,观测频率随时间的变化情况分析结果:00.10.20.30.40.50.60.70.80.91-1-0.50.51时间t/s幅度线性调频信号Time F r e q u e n c y 线性调频信号的STFT 频谱图50010001500200025003000350015. 研究⼀下利⽤⾃相关实现含噪声的正弦信号检测⽅法,并利⽤MATLAB 进⾏验证:答:相关函数的应⽤很⼴,例如,噪声中信号的检测?信号中隐含周期性信号的检测,信号相关性的检测等?设信号)(n f 由正弦信号) (n x 加均值为零的⽩噪声)(n s 所组成,即)()()(n s n x n f +=;那么)(n f 的⾃相关为∑∞=++++=0)]()()][()([1)(n m n s m n x n s n x N m R=)()()()(m R m R m R m R ss sx xs xx +++其中)(m R xs 和)(m R sx 分别是正弦信号)(n x 和⽩噪声)(n s 的互相关?⽩噪声是随机的,和信号)(n x 应⽆相关性,所以)(m R xs 和)(m R sx 应趋近于零?⽩噪声)(n s 的⾃相关函数)(m R ss 主要在n=0处有值,当0||>n 时,衰减很快?由于)(n x 是周期函数,那么)(m R xx 将呈周期变化,从⽽揭⽰出隐含在)(m R xx 中的周期性?由于)(n x 总为有限长,所以这些峰值将是逐渐衰减的,且)(m R xx 的最⼤延迟应⼩于数据长度?01002003004005006007008009001000-4-224含噪声时域正弦信号01002003004005006007008009001000-0.500.5⾃相关检测出的正弦信号6. 简述⼩波滤波的原理,并利⽤MATLAB 中的⼩波⼯具进⾏⼀个⼩波滤波练习,给出计算结果,并进⾏分析答 :信号去噪是信号处理领域的⼀个经典问题,传统的去噪⽅法主要是线性滤波和⾮线性滤波,例如中值滤波和Wiener 滤波等?⼩波变换具有下列良好特性:①低熵性②多分辨率特性③去相关性④选基灵活性?⼩波在信号去噪领域已经取得越来越⼴泛的应⽤?阈值去噪的⽅法是⼀种较好的⼩波去噪法?阈值去噪⽅法的思想就是对⼩波分解后的个层系数中模⼤于和⼩于某阈值的系数进⾏处理,然后对处理完的⼩波系数再进⾏反变换,重构出经过去噪的信号?01002003004005006007008009001000-11原始信号01002003004005006007008009001000-22含噪信号01002003004005006007008009001000-202去噪后的信号。
现代信号处理作业

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 阶数N3dB p Ω≠-/10/1022(/)101,(/)101p s A A N N p c s c ΩΩ=-ΩΩ=-则:⑵.用阶次N 确定 根据公式:在左半平面的极点即为的极点,因而1.1.2 切比雪夫低通滤波器设计步骤: ⑴.确定技术指标归一化: ⑵.根据技术指标求出滤波器阶数N 及:⑶.求出归一化系统函数 其中极点由下式求出:()a H s 1,2,2N()()a a H s H s -()a H s 2,,N p Ωp αs Ωs α/1p p p λ=ΩΩ=/s s p λ=ΩΩε0.12101δε=-p δα=或者由和S 直接查表得2.数字低通滤波器的设计步骤:〔1〕 确定数字低通滤波器的技术指标:通带截止频率pω、通带最大衰减系数pα、阻带截止频率ω、阻带最小衰减系数s α。
〔2〕将数字低通滤波器的技术指标转换成模拟低通滤波器的技术指标。
巴特沃斯:切比雪夫:N ()a H p /ss p λ=ΩΩ0.12101δε=-p δα=〔3〕把模拟滤波器变换成数字滤波器,即把模拟滤波器的系数)(S H 映射成数字滤波器的系统函数)(z H 。
实现系统传递函数s 域至z 域映射有脉冲响应不变法和双线性映射两种方法。
〔〕脉冲响应不变法。
按照技术要求设计一个模拟低通滤波器,得到模拟低通滤波器的传输函数()s H a 转换成数字低通滤波器的系统函数H(z)。
设模拟滤波器的传输函数为()s H a ,相应的单位冲激响应是()t h a ,()s H a =LT[()t h a ],LT[.]代表拉氏变换,对()t h a 进展等间隔采样,采样间隔为T ,得到()nT h a ,将h(n)= ()nT h a 作为数字滤波器的单位取样响应,那么数字滤波器的系统函数H(z)便是h(n)的Z 变换。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
广泛采用。LMS 算法的基本思想:调整滤波器自身参数,使滤波器的输出信号 与期望输出信号之间的均方误差最小。这样系统输出为有用信号的最佳估计。 最小均方(LMS)自适应算法[3]就是以已知期望响应和滤波器输出信号之间误 差的均方值最小为准的, 依据输入信号在迭代过程中估计梯度矢量,并更新权系 数以达到最优的自适应迭代算法。LMS 算法是一种梯度最速下降方法,其显著 的特点和优点是它的简单性, 这种算法不需要计算相应的相关函数,也不需要进 行矩阵运算。 对于横向型滤波结构, 其误差为: e (n) = d (n) −y (n) 均方误差 ε 表示为:ε=E[ e2(n)]=E[d(n)−y(n)]2 对于横向结构的滤波器,有 ε= E[d(n)] +W (n)RW (n) −2 W (n)P
式中:x(k)为自适应滤波器的输入,y(k)为自适应滤波器的输出,d(k)期望输 出信号,e(k)为误差,wi 为滤波器的加权系数,μ 为收敛系数,M 为滤波器阶数。
收敛因子 μ 与滤波器阶数 M 和输入信号的功率都有关系。为使系统收敛, 在输入同一信号的情况下,μ 的取值应该和滤波器的阶数成反比,且应根据不同 的滤波器阶数取不同的步长,这样才能保证有最佳的信号处理结果:当 M 一定 时,μ 是唯一影响 LMS 算法收敛速度的参数,并且随输入信号功率的变化而变 化。Μ 值的选取不能过大,μ 值过大时,在自适应的过程中会引入较大的梯度噪 声,过渡过程将会出现振荡,不能收敛。如果 μ 值太小,虽然梯度噪声降低了, 但是收敛速度较慢。所以对 μ 值要折中考虑。 表 3.1 列出了 LMS 算法的流程。
2 自适应滤波器原理
自适应滤波器与普通滤波器有两个重要区别: 一是自适应滤波器的滤波参数 是可变的[2],它能够随着外界信号特性的变化而动态地改变参数,保持最佳滤波 状态。 自适应滤波器除了普通滤波器的硬件设备以外还有软件部分,即自适应算 法。二是自适应算法决定了自适应滤波器如何根据外界信号的变化来调整参数。 自适应算法的好坏直接影响滤波的效果。 所谓的自适应滤波, 就是利用前一时刻以获得的滤波器参数的结果,自动的 调节现时刻的滤波器参数,以适应信号和噪声未知的或随时间变化的统计特性, 从而实现最优滤波。 自适应滤波器实质上就是一种能调节其自身传输特性以达到 最优的维纳滤波器。自适应滤波器不需要关于输入信号的先验知识,计算量小, 特别适用于实时处理。 自适应滤波器实质上就是一种能调节其自身传输特性以达到最优的维纳滤 波器。自适应滤波器的特性变化是由自适应算法通过调整滤波器系数来实现的。 一般而言,自适应滤波器由两部分组成,一是滤波器结构,二是调整滤波器系数 的自适应算法。 自适应滤波器是由参数可调的数字滤波器 (或称为自适应处理器) 和自适应算法两部分组成。 参数可调的数字滤波器可以是无限冲击响应数字滤波 器(IIR)或有限冲击响应数字滤波器(FIR) ,还可以是格型数字滤波器。图 2.1 给出了自适应滤波器的一般结构,图中 d ( n )为期望响应,x ( n)为自适应滤波器 的输入,( n )为自适应滤波器的输出,e ( n )为估计误差。自适应滤波器的滤波器 系数受误差信号 e ( n )控制,根据 e ( n )的值和自适应算法自动调整。
T 2 T T
(3.1) (3.2) (3.3)
其中,R(n)=x(n) x (n)是 N × N 的自相关矩阵,它是输入信号采样值间的自 相关矩阵。P= d (n) x(n) 为互相关矢量,代表理想信号 d (n) 与输入矢量 x(n) 的 相关性。 在均方误差 ε 达到最小时,得到最佳权系数 W*=[w*0, w*1,…,w*N-1]T。 它满足下列方程 / W n | W (n)= W*=0 即 RW*−P=0 佳值满足 W*= R-1P (3.6) 在有些应用中,把输入信号的采样值分成相同的段(每段为一帧),再求出 R 和 P 的估计值, 以得到每帧的最佳权系数。 这种方法称为块对块的自适应算法。 如语音信号的线性预测编码 LPC,就是把语音信号分成帧进行处理的。 当输入信号和噪声的统计特性未知或输入过程的统计特性发生变化时, 自适 应滤波器自动调整自身参数以满足某种最佳准则要求。根据不同的准则,产生不 同的自适应算法,但主要有两种基本算法:最小均方误差(LMS)算法和递推最小 二乘(RLS)算法。由 Widow 和 Hoff 提出的最小均方误差算法,因其具有计算最 小、易于实现等优点而在实践中被广泛采用。 LMS 算法是一个随机的递推算法,它是用一个带噪声的梯度估计来代替最 陡下降法中的真实度[4]。LMS 算法应包括以下三个方程:
重庆邮电大学研究生堂下考试答卷
2015—2016 学年第二学期
考试科目
现代信号处理
姓
名
年
级
专
业ቤተ መጻሕፍቲ ባይዱ
2016 年 6 月 24 日
基于 LMS 的自适应噪声对消器的设计与仿真
摘要:自适应信号处理的理论和技术已经成为人们经常使用的语音去噪技术,而 Matlab 为其提供了更为方便快捷的方法来对语音信号进行消噪处理。 通过论述自 适应滤波器及 LMS 算法的工作原理, 在对自适应滤波器相关理论研究的基础上, 重点研究了 LMS 自适应滤波算法,并对 LMS 自适应算法进行了分析,设计了 一个二阶自适应噪声对消器,用 Matlab 对其进行了仿真和实现。 关键词:自适应滤波器;LMS 算法;噪声对消
滤波器是一种信号处理系统,它能够过滤或抑制输入信号中的干扰信号,提 取有用信号。滤波器之所以能够滤波,是因为它对不同频率的信号有不同增益, 能将某些频率的信号放大, 而另外的一些频率信号得到抑制。若希望设计的滤波 器能够最大限度的滤除干扰信号, 这样就要设计最佳的频率响应特性和最佳的滤 波参数。 最佳滤波参数和输入信号的特性相关,最佳滤波参数必须根据输入信号 的特性确定。比如必须根据干扰信号和有用信号的频谱设计滤波器的频率响应。 我们不可能设计一个对任何信号都是最佳的滤波器, 这样当输入信号的特性未知 或者随着时间缓慢变化时, 一般的 FIR 和 IIR 两种具有固定滤波系数的滤波器无 法实现最优滤波。解决的办法是引入自适应滤波器,以跟踪信号和噪声的变化。 应用自适应滤波器主要有两种情形:一是输入信号的特性是不变的,但是未 知的。对于这种情形,最佳滤波参数是固定的。在这种情况下,要求自适应滤波 器的参数尽快收敛到最佳滤波参数。一般把参数收敛过程称为“学习”过程。二 是输入信号的特性是“缓慢”变化的。这里的“缓慢”变化是相对于信号幅度变 化而言的。在这种情况下,最佳滤波参数也是“缓慢”变化的,这样要求自适应 滤波参数能尽快“反应”过来,跟随信号特性的变化而改变。这个过程一般称为 “跟踪”过程。 文献给自适应滤波器下的定义[1]为:自适应滤波器是这样的处理器,它在输 入信号特性未知或者输入信号特性变化时,能够调整自己的参数,以满足某种最 佳滤波准则的要求。在有用信号的传输过程中,通常会受到干扰或噪声的污染。 利用滤波技术可以在信号很微弱或用常规的方法无法检测的噪声场中提取所需 要的有用的信号,同时抑制干扰与噪声信号,以便更有效地利用原始信号。自适 应滤波器是相对固定滤波而言的。 固定滤波器是假定输入信号中的有用成分和希 望去掉的成分各占有不同的频带,因此其希望被滤掉的频率是固定的,而其滤波 的频率则是自动适应输入信号而变化的。 在没有任何关于信号和噪声先验知识的 条件下, 自适应滤波器利用前一时刻已获得的滤波器参数来自动调节现时刻的滤 波器参数,以适应信号和噪声未知或随机变化的统计特性,从而实现最优滤波。 在实际情况中, 由于信号和噪声的统计特性常常未知或无法获知,为自适应滤波 器提供了广阔的应用空间。本文主要讨论自适应滤波器在噪声对消方面的原理、 算法及仿真。 根据环境的改变, 使用自适应算法来改变滤波器的参数和结构,这样的滤波 器称为自适应滤波器。 自适应滤波器的系数是由自适应算法更新的时变系数,即 其系数自动连续地适应于给定信号,以获得期望的响应。自适应滤波器的最重要 的特征在于它能够在未知环境中有效工作,并能够跟踪输入信号的时变特征。本 文在理解 LMS 算法实质的基础上对 LMS 算法在自适应噪声对消器中的应用进 行了仿真实现,同时对其收敛性进行了简单分析。
x(n)
参数可调的数字滤波器
y(n) _
d(n) +
自适应算法
∑
e(n)
图 2.1 自适应滤波器原理
3 自适应滤波器实现 LMS 算法
要使自适应滤波器自动调整自身参数,得到有效的输出,则它必须满足某种 最佳准则要求。不同的准则,可以产生不同的自适应算法。目前主要有两种基本 的算法:最小均方误差(LMS)算法和递推最小二乘(RLS)算法。由 Widow 和 Hoff 提出的最小均方误差(LMS)算法,具有计算量小、易于实现等优点而在实践中被
1 绪论
在日常生活中,人们经常受到各种噪声的干扰。例如,在有线电话、无线通信 中回波是不可避免的。回波的存在严重影响了通信的质量。在电视电话系统中, 由本地扩声系统产生的声反馈引起的回音造成了再生混响,影响了语音的清晰度, 严重时会产生自激嚣叫,妨碍整个系统的正常工作。 各种封闭空间的噪声如厂房、 汽车内的噪声等对人体也会产生不利的影响。长期在有噪声的环境中工作,将危 害人的听力、思维、生理和心理。在嘈杂的环境下工作,人们很容易疲乏、反应 迟钝、工作效率降低,并具容易心情烦躁,在噪声的刺激下,人们的注意力不容易 集中,工作容易出错,影响工作速度和工作质量,并且很容易产生错误的判断、进 行错误的操作,降低了生产效率。 在生活中,噪声的存在也很大程度上影响了人们 的休息和放松,降低了生活质量。在如今这个人们不断追求工作效率、生活质量 的年代里,如何有效地消除和抑制噪声己成为人们研究的一个热门课题。 噪声消除是现代信号处理的核心问题之一, 通常实现最优滤波的滤波器为维 纳滤波器与卡尔曼滤波器, 它们均要求已知信号和噪声的先验知识,但在许多实 际应用中往往无法预先得知。1965 年美国斯坦福大学建成了第一个自适应噪声 抵消(ANC)系统,之后随着计算机技术与集成电路技术的进步,新的自适应算 法不断涌现出来, 自适应噪声抵消在理论和应用上都得到了很大的发展。随着科 学技术的不断发展,特别是当前数字通信日益广泛的应用,高性能、高速度、大容 量的数字通信对通信技术和系统提出了更高的要求,这就要求自适应滤波器具有 高性能、高稳定性、高收敛速度以及更加宽广的适用范围。自适应噪声抵消技术 是一种能够消除背景噪声影响的现代信号处理技术。应用自适应噪声抵消技术, 可在未知外界干扰源特征、 传递途径不断变化、背景噪声和被测对象信号相似的 情况下,有效地消除外界噪声的干扰,提高信号传输中的信噪比。这一技术可为 动态信号在测试环境不太理想的工作现场作测试分析和故障诊断提供了有效的 方法和依据, 具有一定的理论和应用价值。自适应噪声抵消技术是基于自适应滤 波原理的一种扩展, 因此, 在研究自适应噪声抵消技术前先要掌握一般自适应滤 波器的设计原理。