频谱分析的思想

频谱分析的思想
频谱分析的思想

谱可能是信号处理中最重要的概念之一了。信号处理与数学分析之间最大的区别之一在于信号处理经常要在时域与频域之间的转换,这种转换称为傅里叶变换,也常叫做频谱分析。在信号处理中,我们经常会遇到频谱、功率谱、能量谱、幅度谱、相位谱等各式各样的关于谱的概念。

那么,谱的概念最初是从哪里得来的呢?和牛顿有关系吗?我们的故事,正要从1666这个科学史上的奇迹年开始。

牛顿疾步走到窗前,并拉上窗帘,实验室变得完全黑暗。然后快速回到他的实验桌边,急切地开始他的实验。虽然他对实验的结果已经确信无疑,但他还是要用事实来证明这个结果。为这个实验,他已经准备了很长时间。

这是公元1669年,牛顿刚刚成为剑桥大学卢卡斯席位的特聘教授。也许是这个席位与生俱来的贵族传统,也许是牛顿的声名太显赫,这个教授席位在日后的科学史中,不仅是剑桥的荣耀,不仅是英国科学家的荣耀,更记载着人类对自然探索的光荣。成为卢卡斯教授后,牛顿开始的第一项研究即是光学,为此甚至不惜推后对微积分理论的进一步完善。实际上,不管是光学,还是微积分,基础都来源于他1665至1666年在老家躲避瘟疫的开创性工作。而对于光学,最开始是因为他发现了一种几乎是完美的研磨透镜的新方法。

利用这些透镜,他发现当白光通过透镜的时候会产生很多颜色的光,也就是色差现象。如何来解释这一现象呢?牛顿自己的理论是白光由各种颜色的光组成,透镜本身并不能产生颜色,他们仅仅是将光的组成部分进行了分解。但是在严格的实验验证之前,他还不能贸然发表他的这些观点。上述实验或许将证明他的理论。

牛顿在仔细检查了一下实验设备。主要设备是两个三棱镜:一个用于将白光分解为各种不同颜色的光,另一个希望能将各种颜色的光恢复成白光。之前,他已经非常仔细地研磨这些三棱镜,保证即便实验不成功也不会是棱镜的原因。一切准备就绪后,他小心地举起光源,并将其对准第一个三棱镜。然后他看到了他所预想的现象:在两个棱镜之间是七色的彩虹,在第二个棱镜之后,又是白光。然后,又经过多次的反复实验,每次都得到了预想的结果。于是他确信,对这个实验,即便是皇家学会那些最挑剔的批评者也将无话可说了。

牛顿认识到,白光包含了各种颜色的光。他将这些通常用肉眼看不见的各种种颜色的光,看做是幽灵一般。于是,在他的拉丁文手稿中用了specter这个词。后来,经过逐步的演化,人们用spectrum来表示谱:意思是彩虹中的各种颜色。这就是谱(spectrum)的由来。

牛顿的下一步工作本来应该是去深入认识这些不同颜色的光源于其频率的不同。然而很遗憾的是,即便是牛顿这样那个时代最杰出的科学家,也终于没能迈出这一步,主要是因为他坚信光的粒子性而不是波动性。在此之后,光的波粒大战还将延续几百年,并深刻影响着物理学的进程,这已经是题外话了。

将频率与谱联系并统一起来,才构筑了信号处理中我们通常所说的频谱的概念。

功率信号在时间域上是无限的、不收敛,所以无法直接做傅立叶变换。如果对时间T内的信号做傅立叶变换,T在趋于无穷,其实也就是得到了功率信号的频谱,其模的平方也就是功率谱了。

如果这个信号不是确定信号,而是随机信号,那功率普的计算为其自相关函数的傅立叶变换。

不过在实际实现中,通过一段随机信号的采样来计算出其自相关函数,然后做傅立叶变换得到的功率谱,其实和把它看成一段确知信号,做傅立叶变换再取模平方得到的功率谱是一样的。这个我实验过

对功率信号,其傅立叶变换不存在,这句话应该不确切。对于有些能量无限信号,如正弦,是不满足傅氏变换的条件的,但是引入delta函数后,还是可以有傅立叶变换。

对于随机信号,从时间研究是没多大意义(除了大数据),所以有了维纳-辛钦定理。

从物理含义上分析更容易理解一些:

1)时间信号与功率和能量;2)频谱与功率谱和能量谱,它们其实是对同一个信号的不同角度的看法而已。

从时间的角度来看,一个信号有实时的时间信号x(t),对应地有其功率和能量。而对于功率信号来说,其能量是无限大的,功率是有限大的;对于能量信号来说,能量有限,功率为0。从频率的角度来看,x(f),其功率谱就相当于时间角度的功率,能量与时间角度的能量的含义是相同的。

对一个时域信号进行傅里叶变换,就可以得到信号的频谱,信号的频谱由两部分构成:幅度谱和相位谱。这个关系倒还是简单。那么,什么是功率谱呢?什么又是能量谱呢?功率谱或能量谱与信号的频谱有什么关系呢?

要区分功率谱和能量谱,首先要清楚两种不同类型的信号:功率信号和能量信号。我们从一个具体的物理系统来引出能量信号和功率信号的概念。已知阻值为R的电阻上的电压和电流分别为v(t)和i(t),则此电信号的瞬时功率为:p(t) = v2(t)/R = i2(t)R。我们在做定性分析时,为了方便起见,通常设定电阻R为1,从而就得到归一化的功率值。做定量计算时候可以通过去归一化,即将实际电阻值代入既可以得到实际功率。将上面对电信号功率的思考方法做一个抽象,对信号x(t)定义其瞬时功率为|f (t)|2,在时间间隔(-T/2T/2)内的能量为:

E=int(|f (t)|2 ,-T/2,T/2) (1)

该间隔内的平均功率为:p = E/T (2)

当且仅当f(t)在所有时间上对应所得到的能量不为0且有限时,该信号为能量信号,即(1)

式中的T趋于无穷大的时候E为有限。典型的能量信号如方波信号、三角波信号等。但是有些信号不满足能量信号的条件,如周期信号和能量无限的随机信号,此时就需要用功率来描述这类信号。当且仅当x(t)在所有时间上对应所得到的功率不为0且有限时,该信号为功率信号,即(2)式中的T趋于无穷大的时候p为有限。一般来说,周期信号和随机信号是功率信号,而非周期的确定信号是能量信号。将信号区分为能量信号和功率信号可以简化对各种信号和噪声的数学分析。还有一类信号其功率和能量都是无限的,如f(t) = t,这类信号很少会用到。

了解了信号可能是能量信号也可能是功率信号,或者非能量信号非功率信号(这类信号实际太少,基本可不考虑),就可以很好地理解功率谱和能量谱。对于能量信号,常用能量谱来描述。所谓的能量谱,也称为能量谱密度,是指用密度的概念表示信号能量在各频率点的分布情况。也即是说,对能量谱在频域上积分就可以得到信号的能量。能量谱是信号幅度谱的模的平方,其量纲是焦/赫。对于功率信号,常用功率谱来描述。所谓的功率谱,也称为功率谱密度,是指用密度的概念表示信号功率在各频率点的分布情况。也就是说,对功率谱在频域上积分就可以得到信号的功率。从理论上来说,功率谱是信号自相关函数的傅里叶变换。因为功率信号不满足傅里叶变换的条件,其频谱通常不存在,维纳-辛钦定理证明了自相关函数和傅里叶变换之间对应关系。在工程实际中,即便是功率信号,由于持续的时间有限,可以直接对信号进行傅里叶变换,然后对得到的幅度谱的模求平方,再除以持续时间来估计信号的功率谱。

对确定性的信号,特别是非周期的确定性信号,常用能量谱来描述。而对于随机信号,由于持续期时间无限长,不满足绝对可积与能量可积的条件,因此不存在傅立叶变换,所以通常用功率谱来描述。周期性的信号,也同样是不满足傅里叶变换的条件,常用功率谱来描述,这些在前面已经有所说明。只有如单频正弦信号等很少的特殊的信号,在引入delta函数之后,才可以求解信号的傅里叶变换。

对于用功率谱描述的随机信号而言,白噪声是一个特例。根据定义,白噪声是指功率谱密度在整个频域内均匀分布的噪声。严格地说,白噪声只是一种理想化模型,因为实际噪声的功率谱密度不可能具有无限宽的带宽,否则它的功率将是无限大,是物理上不可实现的。然而,白噪声在数学处理上比较方便,因此它是系统分析的有力工具。一般,只要一个噪声过程所具有的频谱宽度远远大于它所作用系统的带宽,并且在该带宽中其频谱密度基本上可以作为常数来考虑,就可以把它作为白噪声来处理。例如,热噪声和散弹噪声在很宽的频率范围内具有均匀的功率谱密度,通常可以认为它们是白噪声。

频谱分析仪的使用方法

频谱分析仪的使用方法(第一页) 13MHz信号。一般情况下,可以用示波器判断13MHz电路信号的存在与否,以及信号的幅度是否正常,然而,却无法利用示波器确定13MHz电路信号的频率是否正常,用频率计可以确定13MHz电路信号的有无,以及信号的频率是否准确,但却无法用频率计判断信号的幅度是否正常。然而,使用频谱分析仪可迎刃而解,因为频谱分析仪既可检查信号的有无,又可判断信号的频率是否准确,还可以判断信号的幅度是否正常。同时它还可以判断信号,特别是VCO信号是否纯净。可见频谱分析仪在手机维修过程中是十分重要的。 另外,数字手机的接收机、发射机电路在待机状态下是间隙工作的,所以在待机状态下,频率计很难测到射频电路中的信号,对于这一点,应用频谱分析仪不难做到。 一、使用前须知 在使用频谱分析仪之前,有必要了解一下分贝(dB)和分贝毫瓦(dBm)的基本概念,下面作一简要介绍。 1.分贝(dB) 分贝是增益的一种电量单位,常用来表示放大器的放大能力、衰减量等,表示的是一个相对量,分贝对功率、电压、电流的定义如下: 分贝数:101g(dB) 分贝数=201g(dB) 分贝数=201g(dB) 例如:A功率比B功率大一倍,那么,101gA/B=10182’3dB,也就是说,A功率比B功率大3dB, 2.分贝毫瓦(dBm) 分贝毫瓦(dBm)是一个表示功率绝对值的单位,计算公式为: 分贝毫瓦=101g(dBm) 例如,如果发射功率为lmw,则按dBm进行折算后应为:101glmw/1mw=0dBm。如果发射功率为40mw,则10g40w/1mw--46dBm。 二、频谱分析仪介绍 生产频谱分析仪的厂家不多。我们通常所知的频谱分析仪有惠普(现在惠普的测试设备分离出来,为安捷伦)、马可尼、惠美以及国产的安泰信。相比之下,惠普的频谱分析仪性能最好,但其价格也相当可观,早期惠美的5010频谱分析仪比较便宜,国产的安泰5010频谱分析仪的功能与惠美的5010差不多,其价格却便宜得多。 下面以国产安泰5010频谱分析仪为例进行介绍。 1.性能特点 AT5010最低能测到2.24uv,即是-100dBm。一般示波器在lmv,频率计要在20mv以上,跟频谱仪比相差10000倍。如用频率计测频率时,有的频率点测量很难,有的频率点测最不准,频率数字显示不

ASK--FSK--PSK频谱特性分析

分析ASK 、FSK 、PSK 调制信号的频谱特性 ASK(Amplitude-shift Keying):幅移键控 ASK 指的是振幅键控方式。在二进制数字调制中每个符号只能表示0和1(+1或-1)。但在许多实际的数字传输系统中却往往采用多进制的数字调制方式。与二进制数字调制系统相比,多进制数字调制系统具有如下两个特点: 第一:在相同的信道码源调制中,每个符号可以携带log2M 比特信息,因此,当信道频带受限时可以使信息传输率增加,提高了频带利用率。但由此付出的代价是增加信号功率和实现上的复杂性。 第二,在相同的信息速率下,由于多进制方式的信道传输速率可以比二进制的低,因而多进制信号码源的持续时间要比二进制的宽。加宽码元宽度,就会增加信号码元的能量,也能减小由于信道特性引起的码间干扰的影响等。 ASK 这种调制方式是根据信号的不同,调节正弦波的幅度。幅度键控可以通过乘法器和开关电路来实现。载波在数字信号1或0的控制下通或断,在信号为1的状态载波接通,此时传输信道上有载波出现;在信号为0的状态下,载波被关断,此时传输信道上无载波传送。那么在接收端我们就可以根据载波的有无还原出数字信号的1和0。对于二进制幅度键控信号的频带宽度为二进制基带信号宽度的两倍。 设S(t)频谱为S(ω),S2ASK(t)频谱为: 21 ()[()()] 2ASK c c S w s w w s w w =++- 2ASK 信号的频谱是将数字基带频谱中心搬移到载频处,带宽为基带带宽的两倍;又由 ()() n s n s t a g t nT =-∑ 可知,基带信号是由若干基本脉冲组成的, 因而基带信号的带宽完全由基本脉冲带宽决定。2ASK 信号的带宽取决于基带基本脉冲的带宽,是基本脉冲带宽的两倍。设矩形脉冲: 1,||/2()()() 20,s s t T T f t g t f t ≤?=?=-??其它 对其傅里叶变换得()f t 频谱为: sin(/2) ()/2S wT F w W =

9种谱校正方法及matlab代码

9种谱校正方法及matlab 程序代码 采样间隔归一化成1T ?=,采样长度为N .这样FFT 离散谱线为0,1)i X i N =-(,相应的频率分辨率2/(1/)N f N ωπ?=?=. 设FFT 离散谱线局部极高谱线为m (为了数学上简洁,假定从0开始,注意在MATLAB 环境下数组实际操作的是从1开始),记频偏量δωδω=?. 我们需要使用谱线m 和与之相邻一条次高谱线,记这连续两条谱线中左边一条序号为M (当次高谱线在m 左侧时1M m =-,反之M m =). 下面列出若干算法的δ计算公式 1. 加矩形窗的精确谱校正[1] i i i X U jV =+ 111()sin()()cos()M M M M opt M M V V M U U M K U U ωω+++-?+-?=- 1211cos()sin()cos()sin()opt M M opt M M K M Z V U M K M Z V U M ωωωωωω++-???=+?????-?+???=+???+??? 2121 cos()cos()()Z M Z M M m Z Z ωωωδ?+?-?=+-- 2. 加矩形窗情形,采用解析单频模型的幅值比校正[1, 2] 11||()|||| M M M X M m X X δ++=+-+ 3. 加汉宁窗情形,采用解析单频模型的幅值比校正[1, 2] 112||||()|||| M M M M X X M m X X δ++-=+-+ 4. 加矩形窗情形,采用解析单频模型的复比值校正[3] 1 1Re ()M M M X M m X X δ++??=+- ?-?? 5. 加汉宁窗情形,采用解析单频模型的复比值校正[3] 112()M M M M X X M m X X δ+++=+-- 6. 加矩形窗情形,采用解析单频模型的复合复比值校正[3]

AdobeAudition系列教程二频谱分析仪

Adobe Audition系列教程(二):频谱分析仪 频谱分析仪是研究信号频谱特征的仪器,在电子技术一日千里的今天,是研究、开发、调试维修中的有力武器。现代频谱分析仪都趋向于智能化,虚拟仪器技术广泛应用,有些就是以专用的计算机系统为核心设计的。其结果是结构大大简化、性能飞速提高。当然专业的频谱分析仪就比示波器更加昂贵了,业余爱好者更难用上。不过不必灰心,我们可以充分利用Adobe Audition的频谱分析功能,让你拥有精确频谱分析仪的美梦成真! 1. 频谱显示模式 Adobe Audition本身有一种“频谱显示”模式。先打开一段波形,或用《妙用Adobe Audition:数字存储示波器》一文介绍的方法录制一段波形,即可进行频谱分析。这里我们新建一段20秒的对数扫频信号(本文大多选用直接建立的波形,以便了解信号原始波形的标准频谱特征),然后选择“View=>Spectral View”(视图=>频谱),如图1,或点击快捷工具栏的“Toggle between Spectral and Waveform views”(切换频谱视图/波形视图)按扭,即可将波形以频谱显示的方式显示出来,如图2。扫频的频谱显示见图3。 图1

图2 图3 可以看到,横轴为时间,纵轴为频率指示。每个时刻对应的波形频谱都被显示出来了,可以看到扫描速度是指数增加的,即将频率轴取对数时扫描速度是线性的。如图中光标处18秒处频谱指示约11KHz。实际上频谱指示的颜色是代表频谱能量的高低的,颜色从深蓝到红再到黄,指示谱线电平由低到高的变化。这实际上跟地图的地形鸟瞰显示是比较相似的,看图4频谱复杂变化的声音频谱就更容易理解这点了。 图4

用FFT对信号作频谱分析 实验报告

实验报告 实验三:用FFT 对信号作频谱分析 一、 实验目的与要求 学习用FFT 对连续信号和时域离散信号进行谱分析的方法,了解可能出现的分析误差及其原因,以便正确应用FFT 。 二、 实验原理 用FFT 对信号作频分析是学习数字信号处理的重要内容,经常需要进行分析的信号是模拟信号的时域离散信号。对信号进行谱分析的重要问题是频谱分辨率D 和分析误差。频谱分辨率直接和FFT 的变换区间N 有关,因为FFT 能够实现的频率分辨率是2π/N ,因此要求2π/N 小于等于D 。可以根据此式选择FFT 的变换区间N 。误差主要来自于用FFT 作频谱分析时,得到的是离散谱,而信号(周期信号除外)是连续谱,只有当N 较大时,离散谱的包络才能逼近连续谱,因此N 要适当选择大一些。 三、 实验步骤及内容(含结果分析) (1)对以下序列进行FFT 分析: x 1(n)=R 4(n) x 2(n)= x 3(n)= 选择FFT 的变换区间N 为8和16两种情况进行频谱分析,分别打印出幅频特性曲线,并进行讨论、分析与比较。 【实验结果如下】: n+1 0≤n ≤3 8-n 4≤n ≤7 0 其它n 4-n 0≤n ≤3 n-3 4≤n ≤7 0 其它 n

实验结果图形与理论分析相符。(2)对以下周期序列进行谱分析: x4(n)=cos[(π/4)*n]

x5(n)= cos[(π/4)*n]+ cos[(π/8)*n] 选择FFT的变换区间N为8和16两种情况进行频谱分析,分别打印出幅频特性曲线,并进行讨论、分析与比较。 【实验结果如下】: (3)对模拟周期信号进行频谱分析: x6(n)= cos(8πt)+ cos(16πt)+ cos(20πt) 选择采样频率Fs=64Hz,FFT的变换区间N为16、32、64三种情况进行频谱分析,分别打印出幅频特性曲线,并进行讨论、分析与比较。 【实验结果如下】:

倒频谱

倒频谱 倒频谱分析是一种二次分析技术,是对功率谱的对数值进行傅立叶逆变换的结果。其计算公式为: 该分析方法受传感器的测点位置及传输途径的影响小,能将原来频谱图上成族的边频带谱线简化为单根谱线,以便提取、分析原频谱图上肉眼难以识别的周期性信号。但是进行多段平均的功率谱取对数后,功率谱中与调制边频带无关的噪声和其他信号也都得到较大的权系数而放大,降低了信噪比。 本书分7章介绍了振动信号现代分析的理论方法与技术应用。第1章在介绍振动信号的基本概念、方法后,给出了振动信号的时频分析方法。第2章介绍了近十几年来发展和应用的小波分析方法及其工程实现技术。第3章介绍了谐波小波分析的理论方法与应用技术。第4章介绍了振动信号的Hilbert—Huan9变换分析方法,及其与其他方法相结合的发展和应用。第5章给出了振动信号分形分析方法及其与其他方法相结合的研究。第6章介绍了微弱振动信号的混沌识别方法。第7章介绍了振动信号的盲源分离技术与方法。每一章都给出了若干振动信号分析的应用实例。本书可作为具有一定振动分析基础的大学毕业生、研究生和工程技术人员的学习参考书,也可以作为航空航天类、大机械类、力学类、动力工程类专业研究生“机械振动学”课程的后续专业课教学参考书。目录 第1章振动信号的时频分析方法 1.1幅值域分析法 1.2振动信号的时差域分析方法1.2.1信号预处理 1.2.2相关分析 1.3傅里叶变换 1.3.1连续傅里叶变换1.3.2离散傅里叶变换 1.4振动信号的频域分析方法 1.4.1频谱与频谱分析1.4.2经典谱估计方法 1.4.3倒频谱分析 1.4.4加窗与细化分析 1.5多相干分析技术 1.5.1相干函数 1.5.2频段上非独立输入信号的优先级排序1.5.3多相干分析 1.6时频分析 1.6.1时频分析的基本概念 1.6.2信号的时频表示和相平板 1.6.3时频分析的窗函数 1.7短时傅里叶变换 1.7.1短时傅里叶变换的概念 1.7.2离散短时傅里叶变换 1.7.3短时傅里叶变换在振动信号分析中的应用参考文献第2章振动信号的小波分析方法 2.1小波分析的基本概念 2.1.1小波与小波函数 2.1.2从傅里叶变换到小波变换 2.2二进离散小波变换 2.2.1二进小波变换 2.2.2二进小波的构造 2.2.3数字信号的二进小波变换 2.3多分辨分析与正交小波 2.3.1多分辨分析 2.3.2正交小波基 2.3.3Mallat塔形算法 2.3.4多分辨分析的工程实现技术 2.4二进小波变换在振动信号分析中的应用 2.4.1离散数字信号的二进小波变换过程2.4.2信号的频带分离 2.4.3奇异信号检测 2.4.4带噪转子信号的小波消噪2.4.5行驶车辆实际振动信号的小波分析……第3章振动信号的谐波小波分析第4章振动信号分析的Hilbert-Huang变换方法第5章振动信号分析的分形方法第6章微弱振动信号的混沌识别第7章振动信号的盲源分离参考文献

频谱分析报告仪地使用方法

频谱分析仪的使用方法 13MHz信号。一般情况下,可以用示波器判断13MHz电路信号的存在与否,以及信号的幅度是否正常,然而,却无法利用示波器确定13MHz电路信号的频率是否正常,用频率计可以确定13MHz电路信号的有无,以及信号的频率是否准确,但却无法用频率计判断信号的幅度是否正常。然而,使用频谱分析仪可迎刃而解,因为频谱分析仪既可检查信号的有无,又可判断信号的频率是否准确,还可以判断信号的幅度是否正常。同时它还可以判断信号,特别是VCO信号是否纯净。可见频谱分析仪在手机维修过程中是十分重要的。 另外,数字手机的接收机、发射机电路在待机状态下是间隙工作的,所以在待机状态下,频率计很难测到射频电路中的信号,对于这一点,应用频谱分析仪不难做到。 一、使用前须知 在使用频谱分析仪之前,有必要了解一下分贝(dB)和分贝毫瓦(dBm)的基本概念,下面作一简要介绍。 1.分贝(dB) 分贝是增益的一种电量单位,常用来表示放大器的放大能力、衰减量等,表示的是一个相对量,分贝对功率、电压、电流的定义如下: 分贝数:101g(dB) 分贝数=201g(dB) 分贝数=201g(dB) 例如:A功率比B功率大一倍,那么,101gA/B=10182’3dB,也就是说,A功率比B功率大3dB, 2.分贝毫瓦(dBm) 分贝毫瓦(dBm)是一个表示功率绝对值的单位,计算公式为: 分贝毫瓦=101g(dBm) 例如,如果发射功率为lmw,则按dBm进行折算后应为:101glmw/1mw=0dBm。如果发射功率为40mw,则10g40w/1mw--46dBm。 二、频谱分析仪介绍 生产频谱分析仪的厂家不多。我们通常所知的频谱分析仪有惠普(现在惠普的测试设备分离出来,为安捷伦)、马可尼、惠美以及国产的安泰信。相比之下,惠普的频谱分析仪性能最好,但其价格也相当可观,早期惠美的5010频谱分析仪比较便宜,国产的安泰5010频谱分析仪的功能与惠美的5010差不多,其价格却便宜得多。 下面以国产安泰5010频谱分析仪为例进行介绍。 1.性能特点 AT5010最低能测到2.24uv,即是-100dBm。一般示波器在lmv,频率计要在20mv以上,跟频谱仪比相差10000倍。如用频率计测频率时,有的频率点测量很难,有的频率点测最不准,频率数字显示不稳定,甚至测不出来。这主要足频率计灵敏度问题,即信号低于20mv频率计就无能为力了,如用示波器测量时,信号5%失真示波器看不出来,在频谱仪上万分之一的失真都能看出来。

信号的频谱分析

实验三信号的频谱分析 方波信号的分解与合成实验 一、任务与目的 1. 了解方波的傅立叶级数展开和频谱特性。 2. 掌握方波信号在时域上进行分解与合成的方法。 3. 掌握方波谐波分量的幅值和相位对信号合成的影响。 二、原理(条件) PC机一台,TD-SAS系列教学实验系统一套。 1. 信号的傅立叶级数展开与频谱分析 信号的时域特性和频域特性是对信号的两种不同的描述方式。对于一个时域的周期信号f(t),只要满足狄利克莱条件,就可以将其展开成傅立叶级数: 如果将式中同频率项合并,可以写成如下形式: 从式中可以看出,信号f(t)是由直流分量和许多余弦(或正弦)分量组成。其中第一项A0/2是常数项,它是周期信号中所包含的直流分量;式中第二项A1cos(Ωt+φ1)称为基波,它的角频率与原周期信号相同,A1是基波振幅,φ1是基波初相角;式中第三项A2cos(Ωt+φ2)称为二次谐波,它的频率是基波的二倍,A2是基波振幅,φ2是基波初相角。依此类推,还有三次、四次等高次谐波分量。 2. 方波信号的频谱 将方波信号展开成傅立叶级数为: n=1,3,5… 此公式说明,方波信号中只含有一、三、五等奇次谐波分量,并且其各奇次谐波分量的幅值逐渐减小,初相角为零。图3-1-1为一个周期方波信号的组成情况,由图可见,当它包含的分量越多时,波形越接近于原来的方波信号,还可以看出频率较低的谐波分量振幅较大,它们组成方波的主体,而频率较高的谐波分量振幅较小,它们主要影响波形的细节。

(a)基波(b)基波+三次谐波 (c)基波+三次谐波+五次谐波 (d)基波+三次谐波+五次谐波+七次谐波 (e)基波+三次谐波+五次谐波+七次谐波+九次谐波 图3-1-1方波的合成 3. 方波信号的分解 方波信号的分解的基本工作原理是采用多个带通滤波器,把它们的中心频率分别调到被测信号的各个频率分量上,当被测信号同时加到多路滤波器上,中心频率与信号所包含的某次谐波分量频率一致的滤波器便有输出。在被测信号发生的实际时间内可以同时测得信号所包含的各频率分量。本实验便是采用此方法,实验中共有5路滤波器,分别对应方波的一、三、五、七、九次分量。 4. 信号的合成 本实验将分解出的1路基波分量和4路谐波分量通过一个加法器,合成为原输入的方波信号,信号合成电路图如图3-1-2所示。 图3-1-2 三、内容与步骤 本实验在方波信号的分解与合成单元完成。 1. 使信号发生器输出频率为100Hz、幅值为4V的方波信号,接入IN端。 2. 用示波器同时测量IN和OUT1端,调节该通路所对应的幅值调节电位器,使该通路输出方波的基波分量,基波分量的幅值为方波信号幅值的4/π倍,频率于方波相同并且没有相位差.(注意:出厂时波形调节电位器已调到最佳位置,其波形基本不失真,基本没有相位差。若实验中发现存在波形失真或有相位差的现象,请适当调节波形调节电位器,使波形恢复正常。) 3. 用同样的方法分别在OUT3、OUT5、OUT7、OUT9端得到方波的三、五、七、九此谐波分量(注意其他谐波分量各参数应当满足式3-1-1所示)。 4. 完成信号的分解后,先后将OUT1与IN1、OUT3与IN2、OUT5与IN3、OUT7与IN4、OUT9与IN5连接起来,即进行谐波叠加(信号合成),分别测量(1)基波与三次谐波;(2)基波、三次谐波与五次谐波;(3)基波、三次谐波、五次谐波与七次谐波;(4)基波、三次谐波、五次谐波、七次谐波与九次谐波合成后的波形。并分别保

基于频域的校正方法及实验设计

2016届毕业(设计)论文 题目基于频域的校正方法及实验设计 专业班级过程自动化 学号 1204160134 学生姓名于春明 第一指导教师陈杰 指导教师职称 学院名称电气信息院 完成日期: 2016年 5月 20日

基于频域的校正方法及实验设计 Correction Method and Experimental Design based on Frequency Domain 学生姓名于春明 第一指导教师陈杰

摘要 在经典控制理论中,系统校正设计,就是在给定的性能指标下,对于给定的对象模型,确定一个能够完成系统满足的静态与动态性能指标要求的控制器,即确定校正器的结构与参数。串联校正控制器的频域设计方法中,使用的校正器有超前校正器、滞后校正器、滞后-超前校正器等。超前校正设计方法的特点是校正后系统的截止频率比校正前的大,系统的快速性能得到提高,这种校正设计方法对于要求稳定性好、超调量小以及动态过程响应快的系统被经常采用。滞后校正设计方法的特点是校正后系统的截止频率比校正前的小,系统的快速性能变差,但系统的稳定性能却得到提高,因此,在系统快速性要求不是很高,而稳定性与稳态精度要求很高的场合,滞后校正设计方法比较适合。滞后-超前校正设计是指既有滞后校正作用又有超前校正作用的校正器设计。它既具有了滞后校正高稳定性能、高精确度的好处,又具有超前校正响应快、超调小的优点,这种设计方法在要求较高的场合经常被采用。 关键词:

ABSTRACT In classical control theory, system design correction, that is, given the performance for a given object model to determine a controller to complete the system meets the static and dynamic performance requirements, namely to determine the structure and corrector parameter. Frequency Domain controller series corrected using correction has lead corrector corrector lag lag - lead Corrector like. Features advanced design correction method is to correct the system cut-off frequency than the fast performance is improved before the correction is large, the system, this correction method for design requires good stability, small overshoot and dynamic process of rapid response systems are often use. Lag compensation design approach is characterized by system after correcting the cutoff frequency than the fast performance before correction is small, the system is deteriorated but stable performance of the system has improved, therefore, the system speed requirements are not high, and the stability and steady-state high precision of the occasion, lag compensation design method is more suitable. Lag - lead Corrector design means there are both lagging corrective action ahead of corrective action calibration designs. It has a lag correction high stability, high accuracy advantages, but also has the leading correction fast response, small overshoot of the advantages of this design approach in demanding situations are often used. Keywords:

Adobe-Audition-系列教程(二):频谱分析仪

Adobe Audition系列教程(二):频谱分析仪 频谱分析仪是研究信号频谱特征的仪器,在电子技术一日千里的今天,是研究、开发、调试维修中的有力武器。现代频谱分析仪都趋向于智能化,虚拟仪器技术广泛应用,有些就是以专用的计算机系统为核心设计的。其结果是结构大大简化、性能飞速提高。当然专业的频谱分析仪就比示波器更加昂贵了,业余爱好者更难用上。不过不必灰心,我们可以充分利用AdobeAudition的频谱分析功能,让你拥有精确频谱分析仪的美梦成真!? 1. 频谱显示模式? Adobe Audition本身有一种“频谱显示”模式。先打开一段波形,或用《妙用Adobe Audition:数字存储示波器》一文介绍的方法录制一段波形,即可进行频谱分析。这里我们新建一段20秒的对数扫频信号(本文大多选用直接建立的波形,以便了解信号原始波形的标准频谱特征),然后选择“View=>Spectral View”(视图=>频谱),如图1,或点击快捷工具栏的“Toggle between Spectral and Waveform views”(切换频谱视图/波形视图)按扭,即可将波形以频谱显示的方式显示出来,如图2。扫频的频谱显示见图3。 图1

图2 图3 可以看到,横轴为时间,纵轴为频率指示。每个时刻对应的波形频谱都被显示出来了,可以看到扫描速度是指数增加的,即将频率轴取对数时扫描速度是线性的。如图中光标处18秒处频谱指示约11KHz。实际上频谱指示的颜色是代表频谱能量的高低的,颜色从深蓝到红再到黄,指示谱线电平由低到高的变化。这实际上跟地图的地形鸟瞰显示是比较相似的,看图4频谱复杂变化的声音频谱就更容易理解这点了。 图4

倒谱分析

倒谱分析 (1).倒频谱的数学描述 倒频谱函数CF(q)(power cepstrum)其数学表达式为: (2.6) CF(q)又叫功率倒频谱,或叫对数功率谱的功率谱。工程上常用的是式(2.6)的开方形式,即: (2.7) C0(q)称为幅值倒频谱,有时简称倒频谱。 倒频谱变量q的物理意义 为了使其定义更加明确,还可以定义: (2.8) 即倒频谱定义为信号的双边功率谱对数加权,再取其傅里叶逆变换,联系一下信号的自相关函数: 看出,这种定义方法与自相关函数很相近,变量q与τ在量纲上完全相同。 为了反映出相位信息,分离后能恢复原信号,又提出一种复倒频谱的运算方法。若信号x(t)的傅里叶变换为X(f): (2.9) x(t)的倒频谱记为: (2.10) 显而易见,它保留了相位的信息。

倒频谱与相关函数不同的只差对数加权,目的是使再变换以后的信号能量集中,扩大动态分析的频谱范围和提高再变换的精度。还可以解卷积(褶积)成分,易于对原信号的分离和识别。 (2).倒频谱的应用 分离信息通道对信号的影响 图2.26对数功率谱关系图。在机械状态监测和故障诊断中,所测得的信号,往往是由故障源经系统路径的传输而得到的响应,也就是说它不是原故障点的信号,如欲得到该源信号,必须删除传递通道的影响。如在噪声测量时,所测得之信号,不仅有源信号而且又有不同方向反射回来的回声信号的混入,要提取源信号,也必须删除回声的干扰信号。若系统的输入为x(t),输出为y(t),脉冲响应函数是h(t),两者的时域关系为: y(t)=x(t)*h(t) 频域为: Y(f)=X(f)*H(f)或Sy(f)=Sx(f)*|H(f)|2 对上式两边取对数,则有: (2.11) 式(2.72)关系如图(2.26)所示,源信号为具有明显周期特征的信号,经过系统特性logGk(f)的影响修正,合成而得输出信号logGy(f)。

安捷伦glenB 频谱分析仪使用说明简介

Agilent E4402B ESA-E Series Spectrum Analyzer 使用方法简介 宁波之猫 2009-6-17

目录

1简介 Agilent ESA-E系列是能适应未来需要的Agilent中性能频谱分析仪解决方案。该系列在测量速度、动态范围、精度和功率分辨能力上,都为类似价位的产品建立了性能标准。它灵活的平台设计使研发、制造和现场服务工程师能自定义产品,以满足特定测试要求,和在需要时用新的特性升级产品。该产品

采用单键测量解决方案,并具有易于浏览的用户界面和高速测量的性能,使工程师能把较少的时间用于测试,而把更多的时间用在元件和产品的设计、制作和查错上。 2.面板 操作区 1.观察角度键,用于调节显示,以适于使用者的观察角度。 2.Esc键,可以取消输入,终止打印。 3.无标识键,实现左边屏幕上紧挨的右边栏菜单的功能。 4.Frequency Channel(频率通道)、Span X Scale(扫宽X刻度)和Amplitude Y scale(幅度Y 刻度)三个键,可以激活主要的调节功能(频率、X轴、Y轴)并在右边栏显示相应的菜单。 5.Control(控制)功能区。 6.Measure(测量)功能区。 7.System(系统)功能区。 8.Marker(标记)功能区。 9.软驱和耳机插孔。 10.步进键和旋钮,用于改变所选中有效功能的数值。 11.音量调节。 12.外接键盘插口。 13.探头电源,为高阻抗交流探头或其它附件提供电源。 14.Return键,用于返回先前选择过的一级菜单。 15.Amptd Ref Out,可提供-20dBm的50MHz幅度参考信号。 16.Tab(制表)键,用于在界限编辑器和修正编辑器中四处移动,也用于在有File菜单键所访问对话 框的域中移动。 17.信号输入口(50Ω)。在使用中,接50ΩBNC电缆,探头上必须串联一隔直电容(30PF左右,陶瓷 封装)。探头实物:

Adobe-Audition-系列教程(二):频谱分析仪

AdobeAudition系列教程(二):频谱分析仪 频谱分析仪是研究信号频谱特征的仪器,在电子技术一日千里的今天,是研究、开发、调试维修中的有力武器。现代频谱分析仪都趋向于智能化,虚拟仪器技术广泛应用,有些就是以专用的计算机系统为核心设计的。其结果是结构大大简化、性能飞速提高。当然专业的频谱分析仪就比示波器更加昂贵了,业余爱好者更难用上。不过不必灰心,我们可以充分利用AdobeAudition的频谱分析功能,让你拥有精确频谱分析仪的美梦成真! 1. 频谱显示模式 AdobeAudition本身有一种“频谱显示”模式。先打开一段波形,或用《妙用Adobe Audition:数字存储示波器》一文介绍的方法录制一段波形,即可进行频谱分析。这里我们新建一段20秒的对数扫频信号(本文大多选用直接建立的波形,以便了解信号原始波形的标准频谱特征),然后选择“View=>Spe ctral View”(视图=>频谱),如图1,或点击快捷工具栏的“Toggle between Spectral and Waveform views”(切换频谱视图/波形视图)按扭,即可将波形以频谱显示的方式显示出来,如图2。扫频的频谱显示见图3。 图1

图2 图3 可以看到,横轴为时间,纵轴为频率指示。每个时刻对应的波形频谱都被显示出来了,可以看到扫描速度是指数增加的,即将频率轴取对数时扫描速度是线性的。如图中光标处18秒处频谱指示约11KHz。实际上频谱指示的颜色是代表频谱能量的高低的,颜色从深蓝到红再到黄,指示谱线电平由低到高的变化。这实际上跟地图的地形鸟瞰显示是比较相似的,看图4频谱复杂变化的声音频谱就更容易理解这点了。

ASK--FSK--PSK频谱特性分析

分析ASK 、FSK 、PSK 调制信号的频谱特性 ASK(Amplitude-shift Keying):幅移键控 ASK 指的是振幅键控方式。在二进制数字调制中每个符号只能表示0和1(+1或-1)。但在许多实际的数字传输系统中却往往采用多进制的数字调制方式。与二进制数字调制系统相比,多进制数字调制系统具有如下两个特点: 第一:在相同的信道码源调制中,每个符号可以携带log2M 比特信息,因此,当信道频带受限时可以使信息传输率增加,提高了频带利用率。但由此付出的代价是增加信号功率和实现上的复杂性。 第二,在相同的信息速率下,由于多进制方式的信道传输速率可以比二进制的低,因而多进制信号码源的持续时间要比二进制的宽。加宽码元宽度,就会增加信号码元的能量,也能减小由于信道特性引起的码间干扰的影响等。 ASK 这种调制方式是根据信号的不同,调节正弦波的幅度。幅度键控可以通过乘法器和开关电路来实现。载波在数字信号1或0的控制下通或断,在信号为1的状态载波接通,此时传输信道上有载波出现;在信号为0的状态下,载波被关断,此时传输信道上无载波传送。那么在接收端我们就可以根据载波的有无还原出数字信号的1和0。对于二进制幅度键控信号的频带宽度为二进制基带信号宽度的两倍。 设S(t)频谱为S(ω),S2ASK(t)频谱为: 21 ()[()()] 2ASK c c S w s w w s w w =++- 2ASK 信号的频谱是将数字基带频谱中心搬移到载频处,带宽为基带带宽的两倍;又由 ()() n s n s t a g t nT =-∑ 可知,基带信号是由若干基本脉冲组成的, 因而基带信号的带宽完全由基本脉冲带宽决定。2ASK 信号的带宽取决于基带基本脉冲的带宽,是基本脉冲带宽的两倍。设矩形脉冲: 1,||/2()()() 20,s s t T T f t g t f t ≤?=?=-??其它 对其傅里叶变换得()f t 频谱为:

频谱仪的简单操作使用方法

R3131A频谱仪简单操作使用方法 一.R3131A频谱仪简介。 R3131A频谱仪是日本ADV ANTEST公司的产品,用于测量高频信号,可测量的频率范围为9K—3GHz。对于GSM手机的维修,通过频谱仪可测量射频电路中的以下电路信号, (维修人员可以通过对所测出信号的幅度、频率偏移、干扰程度等参数的分析,以判断出故障点,进行快速有效的维修): 1.手机参考基准时钟(13M,26M等); 2.射频本振(RFVCO)的输出频率信号(视手机型号而异); 3.发射本振(TXVCO)的输出频率信号(GSM:890M—915M;DCS:1710—1785M); 4.由天线至中频芯片间接收和发射通路的高频信号; 5.接收中频和发射中频信号(视手机型号而异)。 面板上各按键(如图-1所示)的功能如下: A区:此区按键是其他区功能按键对应的详细功能选择按键,例如按下B区的FREQ 键后,会在屏幕的右边弹出一列功能菜单,要选择其中的“START”功能就可通过按下其对 (图-1) B区:此区按键是主要设置参数的功能按键区,包括:FREQ—中心频率; SPAN—扫描频率宽度;LEVEL—参考电平。此区中按键只需直接按下对应键输入数值及单位即可。 C区:此区是数字数值及标点符号选择输入区,其中“1”键的另一个功能是“CAL(校

准)”,此功能要先按下“SHIFT(蓝色键)”后再按下“1”键进行相应选择才起作用; “-”是退格删除键,可删除错误输入。 D 区:参数单位选择区,包括幅度、电平、频率、时间的单位,其中“Hz ”键还有“ENTER(确认)”的作用。 E 区:系统功能按键控制区,较常使用的有“SHIFT ”第二功能选择键,“SHIFT+CONFIG(PRESET )”选择系统复位功能,“RECALL ”调用存储的设置信息键,“SHIFT+RECALL(SA VE )”选择将设置信息保存功能。 F 区:信号波形峰值检测功能选择区。 G 区:其他参数功能选择控制区,常用的有“BW ”信号带宽选择及“SWEEP ”扫描时间选择,“SWEEP ”是指显示屏幕从左边到右边扫描一次的时间。 显示屏幕上的信息(如图-2所示)。 二.一般操作步骤。[“ ”表示的是菜单面板上直接功能按键,“ ” 表 示单个菜单键的详细功能按键(在显示屏幕的右边)]: 1) 按Power On 键开机。 2) 每次开始使用时,开机30分钟后进行自动校准,先按 Shift+7(cal ) ,再选择 cal all 键,校准过程中出现“Calibrating ”字样,校准结束后如通过则回复校准前状态。校准过程约进行3分钟。 3) 校准完成后首先按 FREQ 键,设置中心频率数值,例如需测中心频率为902.4M 的信

实验:典型信号频谱分析报告

实验3.2 典型信号频谱分析 一、 实验目的 1. 在理论学习的基础上,通过本实验熟悉典型信号的波形和频谱特征,并 能够从信号频谱中读取所需的信息。 2. 了解信号频谱分析的基本方法及仪器设备。 二、 实验原理 1. 典型信号及其频谱分析的作用 正弦波、方波、三角波和白噪声信号是实际工程测试中常见的典型信号,这些信号时域、频域之间的关系很明确,并且都具有一定的特性,通过对这些典型信号的频谱进行分析,对掌握信号的特性,熟悉信号的分析方法大有益处,并且这些典型信号也可以作为实际工程信号分析时的参照资料。本次实验利用DRVI 快速可重组虚拟仪器平台可以很方便的对上述典型信号作频谱分析。 2. 频谱分析的方法及设备 信号的频谱可分为幅值谱、相位谱、功率谱、对数谱等等。对信号作频谱分析的设备主要是频谱分析仪,它把信号按数学关系作为频率的函数显示出来,其工作方式有模拟式和数字式二种。模拟式频谱分析仪以模拟滤波器为基础,从信号中选出各个频率成分的量值;数字式频谱分析仪以数字滤波器或快速傅立叶变换为基础,实现信号的时—频关系转换分析。 傅立叶变换是信号频谱分析中常用的一个工具,它把一些复杂的信号分解为无穷多个相互之间具有一定关系的正弦信号之和,并通过对各个正弦信号的研究来了解复杂信号的频率成分和幅值。 信号频谱分析是采用傅立叶变换将时域信号x(t)变换为频域信号X(f),从而帮助人们从另一个角度来了解信号的特征。时域信号x(t)的傅氏变换为: 式中X(f)为信号的频域表示,x(t)为信号的时域表示,f 为频率。 3. 周期信号的频谱分析 周期信号是经过一定时间可以重复出现的信号,满足条件: dt e t x f X ft j ?+∞ ∞--=π2)()(

多种频谱校正方法及matlab代码

多种频谱校正方法 采样间隔归一化成1T ?=,采样长度为N .这样FFT 离散谱线为0,1)i X i N =-(,相应的 频率分辨率2/(1/)N f N ωπ?=?=.设FFT 离散谱线局部极高谱线为m (为了数学上简洁,假定从0开始,注意在MATLAB 环境下数组实际操作的是从1开始),记频偏量δωδω=?.我们需要使用谱线m 和与之相邻一条次高谱线,记这连续两条谱线中左边一条序号为M (当次高谱线在m 左侧时1M m =-,反之M m =). 下面列出若干算法的δ计算公式 1.加矩形窗的精确谱校正[1] i i i X U jV =+111()sin()()cos() M M M M opt M M V V M U U M K U U ωω+++-?+-?=-1211cos()sin()cos()sin()opt M M opt M M K M Z V U M K M Z V U M ωωωωωω++-???=+???? -?+???=+???+??? 2121 cos()cos()()Z M Z M M m Z Z ωωωδ?+?-?=+--2.加矩形窗情形,采用解析单频模型的幅值比校正[1,2]11||()||||M M M X M m X X δ++= +-+3.加汉宁窗情形,采用解析单频模型的幅值比校正[1,2] 112||||()|||| M M M M X X M m X X δ++-=+-+4.加矩形窗情形,采用解析单频模型的复比值校正[3] 11Re ()M M M X M m X X δ++??=+- ?-?? 5.加汉宁窗情形,采用解析单频模型的复比值校正[3] 112()M M M M X X M m X X δ+++=+--6.加矩形窗情形,采用解析单频模型的复合复比值校正[3]

频谱分析(完整版)

Matlab 信号处理工具箱 帮助文档 谱估计专题 翻译:无名网友 & Lyra 频谱分析 Spectral estimation (谱估计)的目标是基于一个有限的数据集合描述一个信号的功率(在频率上的)分布。功率谱估计在很多场合下都是有用的,包括对宽带噪声湮没下的信号的检测。 从数学上看,一个平稳随机过程n x 的power spectrum (功率谱)和correlation sequence (相关序列)通过discrete-time Fourier transform (离散时间傅立叶变换)构成联系。从normalized frequency (归一化角频率)角度看,有下式 ()()j m xx xx m S R m e ωω∞ -=-∞ = ∑ 注:()()2 xx S X ωω=,其中()/2 /2 1lim N j n n N n N X x e N ωω→∞=-=∑πωπ-<≤。其matlab 近似为X=fft(x,N)/sqrt(N),在下文中()L X f 就是指matlab fft 函数 的计算结果了 使用关系2/s f f ωπ=可以写成物理频率f 的函数,其中s f 是采样频率 ()()2/s jfm f xx xx m S f R m e π∞ -=-∞ = ∑ 相关序列可以从功率谱用IDFT 变换求得: ()()()/2 2//2 2s s s f jfm f j m xx xx xx s f S e S f e R m d df f πωπ π ωωπ--= =? ? 序列n x 在整个Nyquist 间隔上的平均功率可以表示为 ()()() /2 /2 02s s f xx xx xx s f S S f R d df f π π ωωπ--= =? ? 上式中的 ()()2xx xx S P ωωπ= 以及()() xx xx s S f P f f =

频谱分析

2.1频谱分析原理 时域分析只能反映信号的幅值随时间的变化情况,除单频率分量的简单波形外,很难明确提示信号的频率组成和各频率分量大小,而频谱分析能很好的解决此问题。由于从频域能获得的主要是频率信息,所以本节主要介绍频率(周期)的估计与频谱图的生成。 2.2.1DFT与FFT 对于给定的时域信号y,可以通过Fourier变换得到频域信息Y。Y可按下式计算 式中,N为样本容量,Δt = 1/Fs为采样间隔。 采样信号的频谱是一个连续的频谱,不可能计算出所有的点的值,故采用离散Fourier变换(DFT),即 式中,Δf = Fs/N。但上式的计算效率很低,因为有大量的指数(等价于三角函数)运算,故实际中多采用快速Fourier变换(FFT)。其原理即是将重复的三角函数算计的中间结果保存起来,以减少重复三角函数计算带来的时间浪费。由于三角函数计算的重复量相当大,故FFT能极大地提高运算效率。 2.2.2 频率、周期的估计 对于Y(kΔf),如果当kΔf = 时,Y(kΔf)取最大值,则为频率的估计值,由于采样间隔的误差,也存在误差,其误差最大为Δf / 2。 周期T=1/f。 从原理上可以看出,如果在标准信号中混有噪声,用上述方法仍能够精确地估计出原标准信号的频率和周期,这个将在下一章做出验证 2.2.3 频谱图 为了直观地表示信号的频率特性,工程上常常将Fourier变换的结果用图形的方式表示,即频谱图。 以频率f为横坐标,|Y(f)|为纵坐标,可以得到幅值谱;

以频率f为横坐标,arg Y(f)为纵坐标,可以得到相位谱; 以频率f为横坐标,Re Y(f)为纵坐标,可以得到实频谱; 以频率f为横坐标,Im Y(f)为纵坐标,可以得到虚频谱。 根据采样定理,只有频率不超过Fs/2的信号才能被正确采集,即Fourier 变换的结果中频率大于Fs/2的部分是不正确的部分,故不在频谱图中显示。即横坐标f ∈[0, Fs/2] 2.5.运行实例与误差分析 为了分析软件的性能并比较时域分析与频域分析各自的优势,本章给出了两种分析方法的频率估计的比较,分析软件的在时域和频域的计算精度问题。2.5.1标准正弦信号的频率估计 用信号发生器生成标准正弦信号,然后分别进行时域分析与频域分析,得到的结果如图 4所示。从图中可以看出,时域分析的结果为f = 400.3702Hz,频域分析的结果为f = 417.959Hz,而标准信号的频率为400Hz,从而对于标准信号时域分析的精度远高于频域分析的精度。 2.5.2 带噪声的正弦信号的频率估计 先成生幅值100的标准正弦信号,再将幅值50的白噪声信号与其混迭,对最终得到的信号进行时域分析与频域分析,结果如图 5所示,可以看出,时域分析的结果为f = 158.9498Hz,频域分析的结果为f = 200.391Hz,而标准信号的频率为200Hz,从而对于带噪声的正弦信号频域分析的精度远高于时域分析的精度。 2.5.3 结果分析与结论

相关文档
最新文档