现代信号处理大作业

合集下载

现代信号处理作业

现代信号处理作业

信号时频分析技术及matlab仿真电路与系统王冠军 201128013926153摘要:本文介绍了时频分析的一些基础理论,对短时傅里叶变换Wigner-Ville分布做了简单介绍,运用MATLAB语言实现了旨在构造一种时间和频率的密度函数,以揭示信号中所包含的频率分量及其演化特性的wigner-ville分布。

并对时频分析方法的优缺点进行了分析。

关键词:时频分析短时傅里叶变换wigner-ville分布1 引言基于Fourier变换的传统信号处理技术从信号频域表示及能量的频域分布的角度揭示了信号在频域的特征。

但Fourier变换是一种整体变换,只能为人们提供信号在时域或频域的全局特性而无法了解信号频谱随时间变化的情况。

因此,需要使用一种时间和频率的联合函数来表示信号,这种表示简称为信号,也就是信号的时频分析。

2 时频分析方法信号时频分析主要研究非平稳信号或时变信号的频谱含量是怎样随时间变化的。

时频分析是当今信号处理领域的一个主要研究热点,目前常用时频分析方法主要有短时傅里叶变换、Gabor展开、小波变换、Wigner-Ville分布。

本文主要介绍了短时傅里叶变换和Wigner-Ville分布两种分析方法。

2.1 短时傅立叶变换STFT从历史上看,信号的时频分析用的最多的是短时傅立叶变换,这种变换的基本思想是用一个窗函数乘时间信号,该窗函数的时宽足够窄,使取出的信号可以被看成是平稳的,然后进行的傅立叶变换可以反映该时宽中的频谱变化规律,如果让窗随时间轴移动,可以得到信号频谱随时间变化的规律。

对于时变信号,了解不同时刻附近的频域特征是至关重要的。

因此,人们采用时间—频率描述时变信号,将一维的时域信号映射到一个二维的时域平面,全面反映观测信号的时频联合特征。

短时傅立叶变换反映了这一思想,对于时变信号,采用某一滑动窗函数截取信号,并认为这些信号是准平稳的,然后,再分别对其进行傅立叶变换,构成时变信号的时变谱。

现代信号处理大作业王成志1

现代信号处理大作业王成志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 分布及其模糊函数,画出二者的波形图,指出并分析其信号项和交叉项。

信号处理原理与应用答案

信号处理原理与应用答案

信号处理原理与应用答案【篇一:《现代信号处理理论与应用》作业题】txt>(2)纸质版(手写),认真,不得抄袭。

《现代信号处理理论与应用》作业题(1)推导《现代信号处理》(张贤达版)p37页公式(2-5-8)。

(2)推导矢量参数估计的cramer-rao不等式,并讨论等号成立条件。

(3)令观测样本为xi?s??i,(i?1,?,n)其中{?i}是一高斯白噪声,其均值为零,方差为1。

证明:s的极大似然估计是无偏的和一致的。

(4)若信号满足:s1(t)?a1co?st?a2co2s?t???apcosp?ts2(t)?b1sin?t?b2sin2?t???bpsinp?t观测信号为x(t)?s1(t)?s2(t)?n(t),n(t)是均值为0,均方差为1的高斯白噪声。

计算a1,a2,?,ap,b1,b2,?,bp的最小二乘估计。

(5) 输入信号x(t)为高斯-马尔可夫信号s(t)和噪声n(t)的叠加,信号和噪声假定不相关,其功率谱分别为ss(?)?3和sn(?)?1。

计算?分别取0,+1和-1时,物理不可实现21??维纳滤波器的冲激响应和最小均方误差。

(6) 设信号满足的状态方程及观测方程分别为:?11?xk?1???xk?uk 01 ??yk?1?[1,0]xk?1?nk?1{uk,k?0}和{nk?1,k?0}是均值为零的高斯白噪且与初始状态x0独立,并有?50??10?,var{n}?r?2?(?1)k?1,初始状态的方程矩阵。

vx0??var{uk}?q??k?1k?1???010??01?计算卡尔曼滤波增益m(k)。

(7)接收信号为x(t)?s(t)?n(t),其中n(t)是功率谱密度为n0的高斯白噪。

信号为 2t???2s(t)??2e,t?0。

??0,t?0求匹配滤波器传输函数及其脉冲响应,该匹配滤波器是物理可实现吗?有无可能将它变为物理可实现?若可能,求出滤波器的传递函数及其输出信噪比最大时刻,若不可能,说明理由。

现代信号处理作业

现代信号处理作业

现代信号处理作业现代信号处理课程作业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. 请用多层感知器(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),并画出学习曲线。

信号处理大作业

信号处理大作业
(4)分析特点:自相关的频谱分析与FFT分析类似,经过自变换的信号可以直接与原始信号进行对比分析,并能直观的看出其相似程度,且自相关图像在0点的值最大;之后变小,
若信号中有周期成分,则自相关函数也有周期性。例题中正弦信号的自相关函数为余弦函数。
FFT变换能有效降低时域信号的信噪比阀值,即具有良好的信号识别能力。能够看出信号幅值、相位的变化情况。而自相关变换却不可以。
s=s0+v;
subplot(2,1,1),plot(t,s),title('输入信号')
y=filtfilt(bz,az,s);
y1=filter(bz,az,s);
subplot(2,1,2),plot(t,y,t,y1,':'),title('输出信号'),xlabel('时间/s')
legend('filtfilt','filter')
bz =
1.0e-003 *
0.0604 0.3020 0.6039 0.6039 0.3020 0.0604
az =
1.0000 -3.9823 6.4281 -5.2458 2.1610 -0.3591
频率特性图像为2-1
图2-1数字滤波器的频率特性图形
Chebyschev滤波器源程序:
Fs=1000; %采样频率
string=['Boxcar,N=',num2str(N)];
h=hd.*(B)'; %得到FIR数字滤波器
[H,w]=freqz(h,[1],1024); %求其频率响应
db=20*log10(abs(H)+eps); %得到幅值

现代信号处理大作业

现代信号处理大作业

现代信号处理大型作业一.试用奇阶互补法设计两带滤波器组(高、低通互补),进而实现四带滤波器组;并画出其频响。

滤波器设计参数为:F p =1.7KHz , F r =2.3KHz , F s =8KHz , A rmin ≥70dB 。

(一)、分析与通常的滤波器相比,互补滤波器具有优良的结构特性和结构特性,具有较低的噪声能量和系数敏感性,其定义如下:一组滤波器H 12(),(),.......()Z H Z H Z n 如果满足下式:He Kjw k n(),==∑110<w<2π 则称这组滤波器为幅度互补滤波器;如果满足下式:He kjw k n()=∑=121, 0<w<2π则称这组滤波器为功率互补滤波器,同时互补滤波器还应该满足:Hz A z kk n()()=∑=1其中A(z)为全通函数,适当的选择全通函数,可以使两带函数具有所需要的低通和高通特性。

(二)、设计步骤(1) 对Fp 、Fr 进行预畸);();(''FsFrtg FsFptg r p ∏=Ω∏=Ω(2) 计算'''*r p c ΩΩ=Ω,判断'c Ω是否等于1,即该互补滤波器是否为互补镜像滤波器(3) 计算相关系数⎪⎩⎪⎨⎧-==+++=+-=-=ΩΩ=--=偶数)N 为(;21奇数)N 为 (;;lg /)16/1lg(;150152;1121;1;;])110)(110[(1213090500''02'''211-min1.0min1.0i i u q k N q q q q q k k q k k k k rp Ar Ap;)2cos()1(21))12(sin()1(21)1(21'2∑∑∞=∞=+-++-=Ωm mm m m m m i u Nm q u Nm q q ππ;42⎥⎦⎤⎢⎣⎡=N N;221N N N -⎥⎦⎤⎢⎣⎡=;)/1)(1(2'2'k k v i i i Ω-Ω-=12'1212,1;12N i v i i i =Ω+=--α 22'22,1;12N i v iii =Ω+=β (4) 互补镜像滤波器的数字实现;22i ii A αα+-=;22iii B ββ+-=1221,1;1)(N i ZA Z A Z H i i i =++=∏--22212,1;1)(N i ZB Z B Z Z H i i i =++=∏--- )];()([21)(21Z H Z H Z H L +=(三)、程序与结果 1. 二带滤波器组 (1) 源程序: clear; clf;Fp=1700;Fr=2300;Fs=8000; Wp=tan(pi*Fp/Fs); Wr=tan(pi*Fr/Fs); Wc=sqrt(Wp*Wr); k=Wp/Wr;k1=sqrt(sqrt(1-k^2)); q0=0.5*(1-k1)/(1+k1);q=q0+2*q0^5+15*q0^9+150*q0^13; N=11;N2=fix(N/4); M=fix(N/2); N1=M-N2; for jj=1:M a=0;for m=0:5a=a+(-1)^m*q^(m*(m+1))*sin((2*m+1)*pi*jj/N);%N is odd, u=j end ab=0;for m=1:5b=b+(-1)^m*q^(m^2)*cos(2*m*pi*jj/N); end bW(jj)=2*q^0.25*a/(1+2*b);V(jj)=sqrt((1-k*W(jj)^2)*(1-W(jj)^2/k)); endfor i=1:N1alpha(i)=2*V(2*i-1)/(1+W(2*i-1)^2); endfor i=1:N2beta(i)=2*V(2*i)/(1+W(2*i)^2); endfor i=1:N1a(i)=(1-alpha(i)*Wc+Wc^2)/(1+alpha(i)*Wc+Wc^2); endfor i=1:N2b(i)=(1-beta(i)*Wc+Wc^2)/(1+beta(i)*Wc+Wc^2); endw=0:0.0001:0.5;LP=zeros(size(w));HP=zeros(size(w));for n=1:length(w)z=exp(j*w(n)*2*pi);H1=1;for i=1:N1H1=H1*(a(i)+z^(-2))/(1+a(i)*z^(-2)) ;endH2=1/z;for i=1:N2H2=H2*(b(i)+z^(-2))/(1+b(i)*z^(-2));endLP(n)=abs((H1+H2)/2);HP(n)=abs((H1-H2)/2);endplot(w,LP,'b',w,HP,'r');hold on;xlabel('digital frequency');ylabel('amptitude');(2)运行结果:见图1图1 二带数字滤波器组2.四带滤波器组(1)源程序:clf;Fp=1700;Fr=2300;Fs=8000;Wp=tan(pi*Fp/Fs);Wr=tan(pi*Fr/Fs);Wc=sqrt(Wp*Wr);k=Wp/Wr;k1=sqrt(sqrt(1-k^2));q0=0.5*(1-k1)/(1+k1);q=q0+2*q0^5+15*q0^9+150*q0^13;N=11;N2=fix(N/4);M=fix(N/2);N1=M-N2;for jj=1:Ma=0;for m=0:5a=a+(-1)^m*q^(m*(m+1))*sin((2*m+1)*pi*jj/N); % N is odd, u=jendb=0;for m=1:5b=b+(-1)^m*q^(m^2)*cos(2*m*pi*jj/N);endW(jj)=2*q^0.25*a/(1+2*b);V(jj)=sqrt((1-k*W(jj)^2)*(1-W(jj)^2/k));Endfor i=1:N1alpha(i)=2*V(2*i-1)/(1+W(2*i-1)^2);endfor i=1:N2beta(i)=2*V(2*i)/(1+W(2*i)^2);endfor i=1:N1a(i)=(1-alpha(i)*Wc+Wc^2)/(1+alpha(i)*Wc+Wc^2);endfor i=1:N2b(i)=(1-beta(i)*Wc+Wc^2)/(1+beta(i)*Wc+Wc^2);endw=0:0.0001:0.5;LLP=zeros(size(w));LHP=zeros(size(w));HLP=zeros(size(w));HHP=zeros(size(w));for n=1:length(w)z=exp(j*w(n)*2*pi);H1=1;for i=1:N1H1=H1*(a(i)+z^(-2))/(1+a(i)*z^(-2)) ;endH21=1;for i=1:N1H21=H21*(a(i)+z^(-4))/(1+a(i)*z^(-4)) ;H2=1/z;for i=1:N2H2=H2*(b(i)+z^(-2))/(1+b(i)*z^(-2));endH22=1/(z^2);for i=1:N2H22=H22*(b(i)+z^(-4))/(1+b(i)*z^(-4));endLP=((H1+H2)/2);HP=((H1-H2)/2);LLP(n)=abs((H21+H22)/2*LP);LHP(n)=abs((H21-H22)/2*LP);HHP(n)=abs((H21+H22)/2*HP);HLP(n)=abs((H21-H22)/2*HP);endplot(w,LLP,'b',w,LHP,'r',w,HLP,'k',w,HHP,'m')hold onxlabel('digital frequency');ylabel('amptitude');(2)运行结果:见图2图2 四带数字滤波器组二、根据《现代数字信号处理》第四章提供的数据,试用如下方法估计其功率谱,并画出不同参数情况下的功率谱曲线:1)Levison算法2)Burg算法3) ARMA 模型法 4) MUSIC 算法 1 Levinson 算法Levinson 算法用于求解Yule-Walker 方程,是一种按阶次进行递推的算法,即首先以AR (0)和AR (1)模型参数作为初始条件,计算AR (2)模型参数;然后根据这些参数计算AR (3)参数,等等,一直到计算出AR (p )模型参数为止,需要的运算量数量级为2p ,其中p 为AR 模型的阶数。

现代信号处理试题及答案总结汇编

现代信号处理试题及答案总结汇编

P29采样、频率混叠,画图说明将连续信号转换成离散的数字序列过程就是信号的采样。

它包含了离散和量化两个主要步骤。

若采样间隔Δt 太大,使得平移距离2π/Δt 过小。

移至各采样脉冲函数对应频域序列点上的频谱X(ω)就会有一部分相互重叠,由此造成离散信号的频谱与原信号频谱不一致,这种现象称为混叠。

P33列举时域参数(有量纲和无量纲),说明其意义与作用。

有量纲参数指标包括方根幅值、平均幅值、均方幅值和峰值四种。

无量纲参数指标包括了波形指标、峰值指标、脉冲指标和裕度指标。

偏斜度指标S 表示信号概率密度函数的中心偏离正态分布的程度,反映信号幅值分布相对其均值的不对称性。

峭度指标K 表示信号概率密度函数峰顶的陡峭程度,反映信号波形中的冲击分量的大小。

P37~自相关互相关及作用(举例说明)相关,就是指变量之间的线性联系或相互依赖关系。

信号x (t )的自相关函数:信号中的周期性分量在相应的自相关函数中不会衰减,且保持了原来的周期。

因此,自相关函数可从被噪声干扰的信号中找出周期成分。

在用噪声诊断机器运行状态时,正常机器噪声是由大量、无序、大小近似相等的随机成分叠加的结果,因此正常机器噪声具有较宽而均匀的频谱。

当机器状态异常时,随机噪声中将出现有规则、周期性的信号,其幅度要比正常噪声的幅度大得多。

依靠自相关函数就可在噪声中发现隐藏的周期分量,确定机器的缺陷所在。

(如:自相关分析识别车床变速箱运行状态,确定存在缺陷轴的位置;确定信号周期。

)互相关函数:互相关函数的周期与信号x(t)和y(t)的周期相同,同时保留了两个信号的相位差信息φ。

可在噪音背景下提取有用信息;速度测量;板墙对声音的反射和衰减测量等。

(如:利用互相关分析测定船舶的航速;探测地下水管的破损地点。

P42)P51~蝶形算法FFT 的基本思想是把长度为2的正整数次幂的数据序列{x k }分隔成若干较短的序列作DFT 计算,用以代替原始序列的DFT 计算。

现代信号分析课程大作业.

现代信号分析课程大作业.

课程名称:现代信号分析考试形式:□专题研究报告□论文大作业□综合考试学生姓名:学号:评阅人:时间:年月目录一、FFT算法的特点。

......................................... - 3 -二、平稳随机过程............................................. - 8 -三、功率谱估计.............................................. - 10 -四、自适应滤波器(LMS和RLS)............................... - 16 -五、维纳滤波器.............................................. - 23 -六、FIR维纳滤波器.......................................... - 25 -七、AR/MA滤波器............................................ - 28 -八、卡尔曼滤波器............................................ - 31 -九、小波分析方法及应用...................................... - 37 -十、其他现代信号分析方法.................................... - 43 -现代信号分析- 3 -一、FFT 算法的特点。

题目:请举例说明FFT 算法的特点计算离散傅里叶变换的一种快速算法,简称FFT 。

快速傅里叶变换是1965年由J.W.库利和T.W.图基提出的。

采用这种算法能使计算机计算离散傅里叶变换所需要的乘法次数大为减少,特别是被变换的抽样点数N 越多,FFT 算法计算量的节省就越显著。

根据对序列分解与选取方法的不同而产生了FFT 的多种算法,基本算法是基2DIT 和基2DIF 。

现代信号处理作业

现代信号处理作业

现代信号处理作业2014年6月Cramer-Rao 不等式的证明定义:设),,,(21n x x x θθ=为参数θ的估计子,参数θ的估计子θ 的偏差定义为该估计子误差的期望值,即θθθθθ-=-=}{}{)(E E b def则估计子θ 称为无偏估计子。

若变差)(θ b 等于零或者θθ=}{ E ,即估计子的期望值等于真实参数,则成估计子θ 的渐进无偏估计子。

无偏性反映了参数估计量的取值在真值θ周围摆动程度。

一个参数往往具有不止一个无偏估计子,因此引入了估计量的有效性这一概念来判断偏估计的优劣,以方差大小来衡量无偏估计子的有效性。

假定参数θ存在无偏估计子1θ 和2θ ,若1θ 具有比2θ 更大的方差,即)var()var(21θθ >,则2θ 对于1θ 相对有效。

对于任意参数θ,我们自然希望得到其更为有效的无偏估计子,即要求方差值尽量小。

究竟无偏估计量的方差能够小到什么程度,是否有下界?如果有,是否有估计子能达到这一下届?这就是本次要讨论的Cramer-Rao 下界。

定义:参数θ的估计子θ 的均方误差)(2θ M 定义为该估计子与真实参数的误差平方的期望值,即)()var(}){()(222θθθθθ b E M +=-=定理:(Cramer-Rao 不等式)令),,,(21n x x x x =为样本向量。

若参数估计θ是真实参数θ的无偏估计,并且θθ∂∂)|(x f 和22)|(θθ∂∂x f 存在,则θ 的均方误差所能达到的下界成为Cramer-Rao 下界等于Fisher 信息的倒数,即)(1}){()(22θθθθJ E M ≥-= 不等式中等号成立的充分必要条件是))(()|(ln θθθθθ-=∂∂ K x f 。

其中)(θK 是θ的某个正函数,并与样本n x x x ,,,21 无关。

定义:品质函数的的方差成为Fisher 信息,用)(θJ 表示,定义为)}|(ln {})]|(ln {[)(222θθθθθx f E x f E J ∂∂-=∂∂=证明Cramer-Rao 不等式:由假设条件可得θθ=}{ E ,则可得0)|()(}{=-=-⎰∞∞-dx x f E θθθθθ(1) 对式(1)两边求偏导,得0)|()(}{=-∂∂=-∂∂⎰∞∞-dx x f E θθθθθθθ 即可得 0)|()()|(=∂∂-+-⎰⎰∞∞-∞∞-dx x f dx x f θθθθθ(2) 另,由复合函数求导法可得)|()|(ln )|(θθθθθx f x f x f ⎥⎦⎤⎢⎣⎡∂∂=∂∂ (3) 又有 ⎰∞∞-=1)|(dx x f θ (4)将式(3)和式(4)代入到式(2),得1))(|()|(ln =-⎥⎦⎤⎢⎣⎡∂∂⎰∞∞-dx x f x f θθθθθ 可改写为 []1)|()()|()|(ln =-⎥⎦⎤⎢⎣⎡∂∂⎰∞∞-dx x f x f x f θθθθθθ (5) 由Cauchy-Schwartz 不等式,可得1)()|()|()|(ln 2≥-⎥⎦⎤⎢⎣⎡∂∂⎰⎰∞∞-∞∞-dx x f dx x f x f θθθθθθ 则 ⎰⎰∞∞-∞∞-⎥⎦⎤⎢⎣⎡∂∂≥-dx x f x f dx x f )|()|(ln 1)()|(2θθθθθθ (6) 由Cauchy-Schwartz 不等式等号成立的条件可知:当且仅当)|())(()|()|(ln θθθθθθθx f K x f x f -=∂∂ ,不等式可取等号。

北京邮电大学-现代信号处理-作业

北京邮电大学-现代信号处理-作业

2012111203班AaronHwang 学号:20121406191.2 设 x(n) =5cos(0.25二n), n =0,1^|,15,为有限长序列。

(1) 计算16点DFT并画出幅度谱序列解:程序代码如下n=0:15;x=5*cos(0.25*pi* n);figure(1);stem( n, x);xlabel(' n');ylabel('x( n)');X=fft(x);X=abs(X);figure(2);stem( n,X); xlabel('k');ylabel('X(K)');所得图像如下(2) 在给序列后面补16个零后,计算32点DFT并画出DFT畐度谱序列解:程序代码如下n=0:31;n1=0:15;x1=5*cos(0.25*pi* n1);x=[x1 zeros(1,16)];figure(1);stem( n,x);xlabel(' n');ylabel('x (n)');X=fft(x);X=abs(X);figure (2) ;stem( n,X);stem( n,Xk);xlabel('k');ylabel('X(k)');所得图像如下:(4)依据DTFT 与 DFT 之间的关系,解释补零操作对 DFT 的影响。

解:通过补零操作和DFT 点数的增加可以更加准确的得出原信号的 DFS 图像。

1.5 ( 1)设H 1(z) = 1 - r o e -1扯z ,为单零点系统的传递函数,求其幅度响应HMe^)的解析式。

令r 。

=0.95「0 =0.25二,画出幅度响应的草图。

解:由已知H i (z) = 1 - r °e,可知响应为|H i (e j)|=、1 r 。

22r °cos( s •)(3)把DFT 的点数扩大为64,然后重复(2) 解:程序代码如下 n=0:63; n1=0:15; x1=5*cos(0.25*pi* n1) x=[x1 zeros(1,48)]; figure(1);stem( n,x);xlabel(' n');ylabel('x (n)'); X=fft(x); X=abs(X);figure (2) ;stem( n,X);stem( n, X);xlabel('k');ylabel('X(k)');所得图像如下:程序代码如下r0=0.95;w0=0.25*pi;w=li nspace(0,2*pi,100);H=1-rO*exp(-j*wO)*exp(-j*w);plot(w,abs(H));title('H (jw)');所得图像如下:(2)设H2(Z)= --------- [竹为单极点系统的传递函数,求其幅度响应H2(e尬)的解析1 -r0e式。

信号分析处理大作业报告+程序

信号分析处理大作业报告+程序

1.设计方案如下①利用MATLAB中的wavread命令来读入语音信号,将它赋值给某一向量。

再将该向量看作一个普通的信号,对其进行FFT变换实现频谱分析,再依据实际情况对它进行滤波。

对于波形图与频谱图(包括滤波前后的对比图)都可以用MATLAB画出。

②由于音频信号是连续且长度未知,故可以采用N阶低通滤波器。

滤掉低频部分的噪音,剩下的就是原信号了。

③将去噪后的信号写成wav格式的文件可以使用wavwrite函数。

2. 步骤①录制一段歌曲,采用Matlab工具对此音频信号用FFT作谱分析。

②录制一段音频信号并命名为信xinhao1.wav存放在文件夹中。

③使用wavread函数读出此信号。

④用函数FFT进行傅里叶变换,得到波形图,幅值图,频谱图。

⑤加入一个随机高斯噪声,将原始信号与噪声叠加产生加噪之后的声音文件,得到xinhao2.wav文件。

⑥通过N阶低通滤波器对噪声语音滤波,在Matlab中,FIR 滤波器利用函数filter对信号进行滤波,得到xinhao3.wav文件。

首先通过MATLAB工具编程获取音频文件的原始信号波形,原信号幅值和原始信号频谱图如下:然后通过加一个高斯噪声对其分析可得加噪声后信号波形,加噪声后幅值和加噪声后信号频谱图如下:最后再通过N阶低通滤波器对噪声信号滤波,在Matlab中,FIR 滤波器利用函数filter对信号进行滤波,从而得到滤波后信号波形,滤波后幅值和滤波后信号频谱图:程序[x]=wavread('C:\Users\h\Desktop\xinhao1.wav');X=fft(x,2048);figure(1)fs=abs(X);plot(fs);xlabel('HZ');ylabel('|Y(d)|');subplot(2,2,1);plot(x);xlabel('HZ');ylabel('|Y(d)|');title('原始信号波形');subplot(2,2,2);plot(X);xlabel('HZ');ylabel('|Y(d)|');title('原始信号幅值');subplot(2,2,3);plot(fs);xlabel('HZ');ylabel('|Y(d)|');title('原始信号频谱');figure(2)N=length(x); %计算原始语音信号的长度y1=0.05*randn(N,1); %加上一个高斯随机噪声x1=x+y1;x2=fft(x1,2048);mt=abs(x1);plot(mt);xlabel('HZ');ylabel('|Y(d)|');subplot(2,2,1);plot(x1);xlabel('HZ');ylabel('|Y(d)|');title('加噪声后的波形');subplot(2,2,2);plot(x2);xlabel('HZ');ylabel('|Y(d)|');title('加噪声后的幅值');subplot(2,2,3);plot(mt);xlabel('HZ');ylabel('|Y(d)|');title('加噪声后的频谱');ht=43000;bits=16;wavwrite(x1,ht,bits,'C:\Users\wentao\h\xinhao2.wav');%将加噪声的信号保存figure(3)N=15;wc=0.3;[b,a]=butter(N,wc);x3=fft(x);fp=abs(x3);y2=filter(b,a,x);Y1=fft(y2);subplot(2,2,1);plot(y2);xlabel('HZ');ylabel('|Y(d)|');title('滤波后信号的波形');subplot(2,2,2);plot(Y1);xlabel('HZ');ylabel('|Y(d)|');title('滤波后信号的幅值');subplot(2,2,3);plot(fp);xlabel('HZ');ylabel('|Y(d)|');title('滤波后信号的频谱');wavwrite(y2,ht,bits,'C:\Users\h\Desktop\xinhao3.wav');%将滤波之后保存。

现代信号处理大作业

现代信号处理大作业

1. Levinson-Durbin 算法1.1 Levinson-Durbin 算法P 阶AR 模型的差分方程为:1x()()()pi i n a x n i w n =+-=∑,其中()w n 是均值为0的白噪声。

AR 过程的线性预测方法为:先求得观测数据的自相关函数,然后利用Yule-Walker 方程递推求得模型参数,再根据公式求得功率谱的估计。

LD 递推算法是一种采用AR 模型的现代谱估计方法。

LD 递推算法具体计算步骤如下:(1) Yule-Walker 方程的矩阵形式如下所示:2,1,1(0)(1)(2)()(1)(0)(1)(1)0()(1)(2)(0)0xx xx xx xx k xx xx xx x k k xx xx xx xx r r r r k a r r r r k a r k r k r k r σ---⎡⎤⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥--+⎢⎥⎢⎥⎢⎥=⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥--⎣⎦⎣⎦⎣⎦系数矩阵Hxx xx R R =为Hermitian 矩阵,其中k=1,2……p, (i=1…k)表示k 阶时的预测系数,通过归纳递推求解出p 阶的预测系数。

(2) 当p=1时,即一阶递推为:⎥⎦⎤⎢⎣⎡=⎥⎦⎤⎢⎣⎡⎥⎦⎤⎢⎣⎡-01)0()1()1()0(211,1σa R R R R x x x x求解可得:1,01,1211,1(1)1, (0)(0)(1)x x x x R a a R R a R σ==-=+(3) 当2≥p 时,递推公式如下:10,≡p a , *,1,1,kp p p k p k p a K 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矩阵x R 已知,可得到各阶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 MATLAB 算法实现function [a_p var_p] = Levinson_Durbin(x,p) N = length(x); for ii=1:NRxx(ii) = x(1:N-ii+1)*(x(ii:N))'/N; end a(1)=1;a(2)=-Rxx(2)/Rxx(1);for k=1:p-1 % Levinson-Durbin algorithm var(k+1) = Rxx(0+1)+a(1+1:k+1)*Rxx(1+1:k+1)';reflect_coefficient(k+1+1) = -a(0+1:k+1)*(fliplr(Rxx(2:k+1+1)))'/var(k+1); var(k+1+1) = (1-(reflect_coefficient(k+1+1))^2)*var(k+1); a_temp(1) = 1; for kk=1:ka_temp(kk+1) = a(kk+1)+reflect_coefficient(k+1+1)*a(k+1-kk+1); enda_temp(k+1+1) = reflect_coefficient(k+1+1); a = a_temp; enda_p = a; % prediction coeffecients var_p = var(p+1); % prediction error power1.3 仿真结果(1) 当=2p 时,仿真结果如下:图1 p=2时仿真结果图预测系数:222[(0),(1),(2)]a a a =[1,-0.5236,0.5094]; 误差功率:0.9790(2) 当=10p 时,仿真结果如下:图2 p=10时仿真结果图预测系数:101010[(0),(1),(2),]a a a = [1,-0.4956,0.5058,-0.0057,0.0373,-0.0161,0.0276,0.0285,-0.0221,0.0212,-0.0179] 误差功率:0.9963(3) 当=50p 时,仿真结果如下:图3 p=50时仿真结果图预测系数:50505050[(0),(1),(2),,(50)]a a a a =[1,-0.4958,0.5031,…,0.0078,0.0191]误差功率:0.96761.4 结果分析从p=2的仿真结果可以看出估计得到的功率谱与原始功率谱基本吻合,且曲线平滑没有毛刺;p=10时,采用LD 算法进行估计后,得到的功率谱产生振荡;p=50时,得到的功率谱产生较大的振荡。

现代数字信号处理习题

现代数字信号处理习题

()u n 是离散时间平稳随机过程,证明其功率谱()w 0S ≥。

证明:将()u n 通过冲激响应为()h n 的LTI 离散时间系统,设其频率响应()w H 为()001,w-w w 0,w-w wH w⎧<∆⎪=⎨>∆⎪⎩ 输出随机过程()y n 的功率谱为()()()2y S w H w S w =输出随机过程()y n 的平均功率为()()()002011r 022w wy y w w S w dw S w dw πππ+∆-∆==⎰⎰当频率宽度w 0∆−−→时,上式可表示为()()()01r 00y S w w π=∆≥由于频率0w 是任意的,所以有()w 0S ≥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 zn z n ∈-=-=--ξξα步骤3 一步预测误差自相关矩阵NN 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 新息过程自相关矩阵MM H C n Q n C n n P n C n A *2)()()1,()()(∈+-= 步骤5 卡尔曼增益MN H C n A n C n n P n K *1)()()1,()(∈-=- 或)()()()(12n Q n C n P n K H-= 步骤6 状态估计1*1)()()|(ˆ)|(ˆN n n C n n K n x n x∈+=-αξξ步骤7 状态估计自相关矩阵 NN 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 HH +---= 步骤8 重复步骤1-7,进行递推滤波计算 4、经典谱估计方法:直接法:又称为周期图法,它把随机序列x(n)的N 个观测数据视为一能量有限的序列,直接计算x(n)的离散傅里叶变换,得到X(k), 然后再取其幅值的平方,并除以N ,作为序列x(n)的真实功率普估计自相关法 :1949年,Tukey 根据Wiener —Khintchine 定理提出了对有限长数据进行谱估计的自相关法,即利用有限长数据估计自相关函数,再对该自相关函数球傅立叶变换,从而得到谱的估计。

西南交大现代信号处理作业

西南交大现代信号处理作业

现代信号处理作业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 =无关。

现代信号处理及其应用作业

现代信号处理及其应用作业

成绩:现代信号处理及其应用题目:现代信号处理在智能交通系统中的应用学号:121044408姓名:何俊霖2015年6月10日现代信号处理在智能交通系统中的应用摘要:如今科技高速发展,信息处理技术也日新月异,信息处理技术的发展将会影响人们的生活。

本文主要论述了现代信号中的图像信号处理技术在智能交通管理系统中应用的热点和关键问题。

针对智能交通管理系统的多个应用实例,如城市交通智能监控、电子警察系统及流量监测和交通诱导系统等,本文给出了应用视频技术的解决方案,并与传统的交通管理方式作了对比。

关键字:图像信号处理;交通管理系统;视频监控引言在世界经济飞速发展的今天,交通在国民经济中占据着越来越来重要的地位,交通的现代化程度成为衡量一国一城现代化水平的重要标志。

然而,交通的现代化却带来了许多问题:道路的拥挤、事故的增加交通环境的恶化等等。

这就需要有强有力的交通管理手段,而旧有的依靠人为管理的管理系统由于效率低、成本高,显然己不适应交通现状。

第一代交通自动管理系统采用无线遥感、红外、雷达等方式来获取环境信息,并传达给指挥人员,由指挥人员对交通实施管理,采用这种方式比纯人工方式有了巨大进步,减轻了交警的负担, 缓解了交警警力人员不足的问题。

但是采用这种方式也有其缺点,指挥人员接收到的只是信号、数据,对现场状况无法有更深入的了解,决策的正确性更多地依赖于经验,因此还不能说是真正的自动化、智能化。

得益于人工智能技术和硬件技术的飞速发展, 图像和视频技术已经广泛应用于新一代,智能交通系统( IST )中。

人工智能技术使系统的判定更可靠、算法更优化。

硬件技术的发展使我们可以获得更加高质量的图像,并为我们系统的实时性提供了保证。

据统计,正常人从外界获取的信息中75%以上是通过眼睛获得的,而司机则达到了90%以上。

第一代交通自动管理系统是利用微波、遥感等方式获取信息的, 大约相当于利用人的听觉、触觉来代替人的视觉,效果形同盲人摸象。

现代信号处理作业

现代信号处理作业
fpe(m+1)=zaosheng1(m+1)*(N+m+2)/(N-m-2);
if (fpe(m)-fpe(m+1))<1e-6 %用FPE准则定阶
break
end
m=m+1;
end
end
disp(['根据模型求出的自相关序列为:',num2str(R)])
disp(['题目给出的AR模型参数为:',num2str(a1)])
R1 = [R(3),R(2),R(1);R(4),R(3),R(2);R(5),R(4),R(3)];
r = [-R(4);-R(5);-R(6)];
A=linsolve(R1,r);%估计模型参数。
A
% 2)估计ARMA模型中的MA参数。
AA=[1,A(1),A(2),A(3)];
v =filter(AA,1,x);
结果与分析:
根据模型求出的自相关序列为:4.9377 4.3287 4.1964 3.8654 3.6481 3.4027 3.1919 2.986 2.797
题目给出的AR模型参数为:-0.58333 -0.375 0.041667
估计出的AR模型参数为:-0.56844 -0.36884 0.039549
(3)由于已经估计出ARMA模型的所有参数,根据ARMA模型与功率谱等价的关系可以算出功率谱,并用plot函数画出图形。
在这里,噪声方差为一,各阶参数都已经估计出,进而画出图形。
结果与分析:
A =
-0.2630
-0.1870
0.0537
B =
0.4696
0.0317
由结果可知,求解结果与已知系数相近,但存在一定误差,当观测的数据个数不断增大时,结果误差会不断减小,一直到不再用明显变化。

清华大学研究生 现代信号处理大作业__1up

清华大学研究生 现代信号处理大作业__1up

现代信号处理 大作业-1大作业-1描述:求下列假设的最佳判决方式,并求解每种错误的概率。

(提示,可以假设K 非常大,并利用中心极限定理将大量同分布随机变量之和近似为高斯。

)200211222:[](0,):[](0,)=1:[](0,)H r k N H r k N k K H r k N σσσ::L : ,,解:根据题意2022210011(|)exp [](2)2KK k p r H r k πσσ=⎡⎤=-⎢⎥⎣⎦∑ 2122211111(|)exp [](2)2K K k p r H r k πσσ=⎡⎤=-⎢⎥⎣⎦∑2222212211(|)exp [](2)2K K k p r H r k πσσ=⎡⎤=-⎢⎥⎣⎦∑ 假设判决为0H ,则应满足以下条件:00110011(|)()(|)()(|)()(|)()p r H p H p r H p H p r H p H p r H p H ≥⎧⎨≥⎩假设先验概率为0121()=()=()=3p H p H p H ,并假设0120σσσ≥≥>,代入上式得:212210012222100211exp []2211exp []22K K k KK k r k r k σσσσσσσσ==⎧⎡⎤⎛⎫⎛⎫⎪≥-⎢⎥ ⎪ ⎪⎪⎝⎭⎝⎭⎣⎦⎨⎡⎤⎛⎫⎛⎫⎪≥-⎢⎥⎪ ⎪⎪⎝⎭⎝⎭⎣⎦⎩∑∑ 进一步化简可得到:2101220122012202(ln ln )[]1122(ln ln )[]1122K k Kk K r k K r k σσσσσσσσ==-⎧≥⎪-⎪⎪⎨-⎪≥⎪-⎪⎩∑∑即 21020122220102(ln ln )(ln ln )[]max ,11112222Kk K K r k σσσσσσσσ=⎛⎫ ⎪-- ⎪≥ ⎪-- ⎪⎝⎭∑。

根据上述同样的思路,可以得到判决为1H ,需要满足:21021122221201(ln ln )(ln ln )[]11112222K k K K r k σσσσσσσσ=--≤≤--∑判决为2H ,需要满足:22021122220212(ln ln )(ln ln )[]min ,11112222Kk K K r k σσσσσσσσ=⎛⎫⎪-- ⎪≥ ⎪-- ⎪⎝⎭∑ 下面为了易于表述,记211[]()Kk r k f r K==∑。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

4
4.5 x 10
5
³ ³ ³ ³ ³ ³ ³ ³ ³ ³ ³ ³ ³ ³ ³
Yn=fft(Y,L); FY=abs(Yn(1:1:L/2)); figure(3) subplot(211) plot(f,FY); grid on; xlabel('频率(f)');ylabel('幅度(FY)'); title('加噪前音频语音波形的频域图'); Y1n=fft(Y1,L); FY1=abs(Y1n(1:1:L/2)); subplot(212); plot(f,FY1); grid on; xlabel('频率(f)');ylabel('幅度(FY1)'); title('加噪后音频语音波形的频域图');
fp=100;fz=200;fc=150;fs=2000; wc=2*pi*fc/fs;%将模拟滤波器技术指标转换为数字滤波器技术指标 M=ceil(3.44*fs/100);%求窗函数项数M alpha=(M-1)/2; n=[0:1:(M-1)]; m=n-alpha+eps; hd=sin(wc*m)./(pi*m);%理想冲击响应 w_ham=(hamming(M))';%M项哈明窗 h=hd.*w_ham;%加哈明窗,计算脉冲响应 [H,w]=freqz(h,1);%求频谱 mag=abs(H);%求实部 db=20*log(mag); subplot(221);stem(hd);title('理想冲击响应'); ylabel('hd[n]'); subplot(222);stem(w_ham);title('哈明窗'); ylabel('w[n]'); subplot(223);stem(h);title('实际冲击响应'); ylabel('h[n]'); subplot(224); plot(w/(2*pi),db);title('衰减幅度');ylabel('H[w]');
MATLAB程序
运行结果
0.15 0.1 理想冲击响应 1 0.8 哈明窗
hd[n]
0.05 0 -0.05 0 20 40 60 80
w[n]
0.6 0.4 0.2 0 0 20 40 衰减幅度 60 80
0.15 0.1 0.05 0 -0.05 0
实际冲击响应
100 0
H[w]
20 40 60 80
sf=filter(h,1,s); %使用filter函数对信号进行滤波 subplot(211); plot(t,sf); %滤波后信号的图像 xlabel('时间/s');ylabel('幅度');title('信号滤波后的时域图 '); subplot(212); Fsf=fft(sf,512); %滤波后信号频域图 AFst=abs(Fsf); %信号频域图的幅值 f=(0:255)*fs/512; % 频率采样 plot(f,AFst(1:256)) %滤波后的信号频域图 xlabel('频率/赫兹'); ylabel('幅度'); title('信号滤波后的频 域图');
MATLAB程序
运行结果
1 0.5 幅 度 ( Y) 0 -0.5 -1 加噪前音频语音波形的时域图
0
0.5
1
1.5
2
时 间 ( t)
2.5
3
3.5
4
4.5 x 10
5
3 2 幅 度 ( Y) 1 0 -1 -2 -3 0 0.5 1 1.5
加噪后音频语音波形的时域图
2
时 间 ( t)
2.5
3
3.5
MATLAB程序
运行结果
4 2
幅度
信号滤波前时域图像
0 -2 -4 0 0.05 0.1 0.15 0.2 0.25
300 200
幅度
时 间 /秒 信号滤波前频域图
100 0
0
100
200
300
400 500 600 频 率 /赫 兹
700
800
900
1000
³ ³ ³ ³
³ ³ ³ ³ ³ ³
现代信号处理新方法
时间:2016.5.3(第9周)
实验一
(1)蒲丰投针实验计算PI值
³ ³ ³ ³ ³
在平面上画有一组间距为a的平行线,将一根长 度为l(l<a)的针任意掷在这个平面上,求此 针与平行线中任一条相交的概率 投掷n次 相交次数为m 相交概率p=m/n=2l/(a* π) π=2l/(a*p)
MATLAB程序

运行结果
1 0.9 0.8 0.7 0.6 0.5 0.4 0.3 0.2 0.1 0 0 0.2 0.4 0.6 0.8 1
n
1000 5000 10000 50000 100000
PI
3.1560 3.1368 3.1600
3.1392 3.1475
n=10000运行结果图
实验结论
50 0 -50 H[w] -100 -150 -200
衰减幅度
h[n]
0.2 0.1 0 -0.1 -0.2 0 5 10 15 20 25 30 35 40 45 50
-250
0
0.05
0.1
0.15
0.2
0.25
0.3
0.35
0.4
0.45
0.5
³ ³ ³ ³ ³ ³ ³ ³ ³ ³ ³ ³ ³ ³ ³ ³
³
³ ³
³ ³ ³ ³ ³ ³ ³ ³ ³ ³ ³ ³ ³ ³ ³ ³ ³ ³ ³ ³ ³
fp=1500;fz=2000;fc=1750;fs=7000; wc=2*pi*fc/fs;%将模拟滤波器技术指标转换为数字滤波器技术指标 M=ceil(3.44*fs/500);%求窗函数项数M alpha=(M-1)/2; n=[0:1:(M-1)]; m=n-alpha+eps; hd=sin(wc*m)./(pi*m);%理想冲击响应 w_ham=(hamming(M))';%M项哈明窗 h=hd.*w_ham;%加哈明窗,计算脉冲响应 [H,w]=freqz(h,1);%求频谱 mag=abs(H);%求实部 db=20*log(mag); figure(1) subplot(221);stem(hd);title('理想冲击响应'); ylabel('hd[n]'); subplot(222);stem(w_ham);title('哈明窗'); ylabel('w[n]'); subplot(223);stem(h);title('实际冲击响应'); ylabel('h[n]'); subplot(224); plot(w/(2*pi),db);title('衰减幅度');ylabel('H[w]');
h[n]
-100 -200 -300 0 0.2 0.4 0.6 0.8
叠加的正弦信号
频率分别为f1=50HZ;f2=300HZ;f3=700HZ 采样频率fp=2000HZ
³
³
³ ³ ³ ³ ³ ³ ³ ³ ³ ³ ³ ³ ³ ³ ³ ³ ³
f1=50;f2=300;f3=700;%待滤波正弦信号频率 fs=2000;%采样频率 subplot(211); t=0:1/fs:0.25;%定义时间范围和步长 s=sin(2*pi*f1*t)+sin(2*pi*f2*t)+sin(2*pi*f3*t); plot(t,s);%滤波前的信号图像 xlabel('时间/秒'); ylabel('幅度') title('信号滤波前时域图像'); subplot(212); Fs=fft(s,512);%将信号变换到频域 AFs=abs(Fs);%信号频域图的幅值 f=(0:255)*fs/512;%频率采样 plot(f,AFs(1:256));%滤波前的信号频域图 xlabel('频率/赫兹'); ylabel('幅度') title('信号滤波前频域图')
MATLAB程序
运行结果
0.5 0.4 0.3 hd[n] w[n] 0.2 0.1 0 -0.1 -0.2 0 5 10 15 20 25 30 35 40 45 50 0.2 0 0.6 0.4 理想冲击响应 1 0.8 哈明窗
0
5
10
15
20
25
30
35
40
45
50
0.5 0.4 0.3
实际冲击响应
MATLAB程序
运行结果
2 1
幅度
信号滤波后的时域图
0 -1 -2 0 0.05 0.1 0.15 0.2 0.25
300 200
幅度
时 间 /s 信号滤波后的频域图
100 0
0
100
200
300
400 500 600 频 率 /赫 兹
700
800
900
1000
语音加噪滤波
数字低通滤波器 通带边缘频率fp=1500HZ; 阻带边缘频率fz=2000HZ; 采样频率fs=7000HZ 加入噪声为高斯白噪声
MATLAB程序
运行结果
1500 幅 度 ( FY) 1000 500 0 加噪前音频语音波形的频域图
0
500
相关文档
最新文档