利用Hilbert变换提取信号瞬时特征的算法实现
研究了在工程中如何通过算法来实现利用 Hilbert 变换提取信号的瞬时特征。深入地分析了如何利用数值微分法提高提取瞬时频率特征的精度。最后,给出了一种可行的算法,并通过实验验证了这种方法可以在工程中有效地提取信号的瞬时频率特征。84微机发展第13卷①H(x)=y;H(x)=y;(i=0,1,…n)(j=0,1(11)②在每个小区间/x1,x1+17i=0,1,…,n-1)上由相关定理知:当划分的小区间的长度趋于零时H(x)是三次多项式。s(x)及其一至三阶导数分别一致收敛到f(x)及其一至可以写出分段三次 Hermite插值函数的分段表达式:三阶导数。所以用三次样条插值函数去近似表达用离散值(x)=(1+2x-x过+)2v;+表示的原函数,具有较高的可靠性。3)两种插值的比较挨尔米特 Hermite插值较三次样I-i,1+2条插值具有较好的稳定性与收敛性,但它只能休让各段曲线在连接点上的连续性,而不能保证整条曲线在这些点上y+1Ditl的光滑性。而有时不仅要求曲线连续,而且要求曲线的曲X/(i=0,18)率也连续即要求分段插值函数具有连续的一阶导数,埃H(x)的导数为尔米特 Hermite插值此时就不能满足上述要求6次样条插值较埃尔米特 Hermite插值具有较好的H(x光滑程度,用三次样条插值函数求数值导数比用埃尔米特+2(x-x2(xHermite插值可靠性大,但计算比较复杂,二者的区别见图h2yV+17, h(i-0,12)三次样条插值。已知函数y=f(x)在区间/a,b上的n+1个节点上的值y=f(x;)(i=0,1,…m),求插值函数s(x),使(i=0,1图4 Hermite插值与三次样条插值的比较图2在每个小区间x,x+1(=0.1.…n-1)上利用埃尔米特 Hermite插值得到的2FSK信号的瞬时s(x是三次多项式,记为s(x频率见图5,利用二次样条插值得到的该信号的瞬时频率③3(x)在la,b/上二阶连续可徵。见图6。数s(x)称为f(x)的三次样条插值函数可以利用节点处的二阶导数值为参数,也可以利用节点处的导数值为参数求三次样条插值涵数的表达式。若利用节点处的一阶导数值为参数,求得的三次样条插值函数的表达式为(x)=M-1x-x-)36 h6 hMihi5 DEMeN5a亩pai66hx∈[x;,x+17,b-x+1-x,S"(x)=M图5由 Hermite插值提取图6由三次样条插值提取(j=0,1的2FSK信号的瞬时频率的2FSK信号的瞬时频率对s(x)进行求导,利用S(x)在节点处一阶导数连从图5、图6可以看出利用三次样条插值得到的瞬时续的性质结合边界条件求解出参数M,把求得的参数代频率可以准确反映出信号具有的的摒时频率特征而利用入公式(10),即得三次样条插值函数的s(x)分段表示式。埃尔米特 Hermite插值得到的瞬时频率与信号具有的瞬s;(x)的导数为时频率特征不符。这是因为利用数值微分法求瞬时频率插值以后喫进行求导。三次样条插值函数具有连续的二阶M2 hiM; 2 hj导数,因而具有较好的光滑程度,符合求导条件,所以可以J+1-h(M2+1-M/)准确求出信号的瞬时频率;而埃尔米特 Hernite插值.不够光滑,虽能保证插值多项式收敛于原函数,但不能保证插x Elx,x;+1 h,=xi+I-x, S(xj )=M;值多顷式的导数收敛于原函数的导数,所以求得的值与信o1994-2010ChinaAcademicJournalElectronicPublishingHouse.Allrightsreservedhttp://www.cnki.net第6期刘慧婷等:利用 Hilbert变提取信号瞵时特征的算法实现号实际的瞬时频率值不符。实验结果和理论分析结果是(1) Hilbert变换只能近似应用于窄带信号,即形如纹的(t)=a(1)cosu+6(1)),其中>>B(B为信号带2.3.3结论宽)的信号。但实际应用中,存在许多非窄带信号, Hilbcrt利用数值微分法求瞬时频率ω(t)的步骤可以归纳变换对这些信号无能为力为:首先通过三次样条插值得到分段多项式p(1),(2)对于任意给定时刻,通过 Hilbert变换运算后的结pp(抄);然后分别对分段多项式p(t),Pp()关于变量t果只能存在一个频率值,即只能处理任何时刻为单一频率进行求导,得到pd(,ppd(t);最后求出每一时刻t所对的信号。这显然不合理,因为在实东中同一信号会含有多应的导数值,即求得t(t,u(t)。再把求得的值代入公种频率成分式(6)就完成了提取瞬时频率ω(1)的过程。求解结果见(3)对信号进行 Hilbert变换时,信号的两端会出现严图7重的端点效应。提取某些信号瞬时特征所得的瞬时频率在局部出现了负数,端点效应是造成负频率的一个原因而端点效应可以通过利用特征波对原有数据序列进行延拓的方法来解决,具体解决办法将在今后讨论。尽管目前出现了EMD担论4,其目的是将不满足Hibt变换的信号进行分解得到若干个IMF( intrinsic mode function),然后进行 Hilbert运算,达到提取信号瞬时特征的目的。该理论开辟了信号处理的新空间。但它还不够成熟还需喫进一步的完善和研究图7利用数值微分法提取信号的瞬时频率特征参考文献从图7可以看出,以三次烊条指值进行的数值微分可[]黄长蓉. Hilbert变换及其应用[J].成都气象学院学报以准确岀提取岀信号的瞬时频率特征。199,14(3):273-276.[2]杨小牛,楼A义,徐建良.软件无线电原理与应用[M].北3结束语京:电子工业出版社,2001在工程中, Hilbert变换使得我们对短信号和复杂信号[3]丁丽妤.数值计算方法[M].北京:北京理工大学出版社,的摒时特征的提取成为可能特别是对瞬时频率特征提1997取,在工程中具有十公重要的意义。文中讨论的利用三次[4] Huang N e. The empirical mode decomposition and the hilbert样条插值进行数值徵分以提取瞬时特征的方法是可行的,spectrum for nonlinear and nor stationary time series anal ysis但还存在着如下问题。[].Proc.R.soc.Lond.A,1998,454:903-995(上接第81页)218994。例22(x)=(1-2siny=223101075一般的(A算法计算了120代,求到的最大值为454176.219。154370083改进的α算法计算了34代,求到的最大值为1048575.875。改进后的αA算法收敛速度(指迭代次数)比一般GA算法几乎快了一个数量级,精度也提高了不少,特别是例2的最大值提高一倍多,速度提高这么快是未曾料到的y=74958参考文献+4X Axl Thla[1]陈国良.遗传算法及其应用[M]·北京:人民邮电出版社,图2函数2的图像1996一般GA算法计算了20代,求到的最大值为[2]袁亚湘,孙文瑜.最优化理沦与方法[M]北京:科学出版社,19991.218983[3]张铃,张钹·遗传算法杋理的硏究[J]·软件学报,改进(A算法计算了5代,求到的最大值为2000,11(7):945952o1994-2010ChinaacAdemicJournalElectronicPublishingHouse.Allrightsreservedhttp://www.cnki.net
- 2020-12-05下载
- 积分:1
仿真软件saber说明书 很详细
介绍saber仿真的很好的书籍 建议按步骤学习SABER电气系统培训手册Saber ElectricalSystems Workshop练习指南北京才略科技有限公司TEL:010)82673952/82673953SABER电气系统培训手册介绍Saber Electrical Systems介绍本课程的目的是使用户熟悉 Saber模拟器。课程为 Saber功能培训:这一部分通过相关的简单电路和系统集中讲解如何使用 Saber模拟器的各个功能。内容用户应从以下内容出发:如何应用 Saber改善电路和系统设计如何通过 Saber进行由上至下和由下至上设计如何完成不同类型的分析如何使用 Saber模型如何使用 Saber模型库如何查阅帮助如出错如何解决背景要求需具备基本的工程知识熟练的计算机操作·非必需不需要具备仿真经验北京才略科技有限公司TEL:010)82673952/82673953SABER电气系统培训手册相关说明相关说明Saber book, Saber在线帮助系统,描述了 Saber sketch和 Cosmos Scope的特性以及 Saber的一般功能,比如菜单的使用和打印。还提供了每个 Saber命令和 Saber guide界面的详细信息Sorh8 ing Saberi讲解如何在 Saber sketch中生成电路设计图并如何应用 Saber分析Getting Started Using Saber with the Frameway integrations帮助您分析两个示范电路,并使您熟悉在 Cadence和 Mentor图表的框架环境下主要的使用过程。Analyzing Designs Using saber描述了如何应用 Saber获取图表、模拟设计和优化参数值Saber Design Esamples示范如何应用 Saber对设计图进行仿真和分析。北京才略科技有限公司TEL:010)82673952/82673953SABER电气系统培训手册惯例Saber Electrical systems惯例培训手册采用以下惯例ButtonButton这样的字体用来在用户界面上突出显示按键的描述和编号。ComputerComputer”这样的字体用来突出屏幕上输入内容(即您在命令行或某区域输入内容)。Dialog Menu Form“ DialogMenu Forn”这样的字体表示对话框标题土级菜单标题及表格标题。DocnammeNonAme”这样的字体用来指示印刷手写体标题(同 Computer字样)FieldNameFieldName”这样的字体用来突出区域名称。FilenameFilename”这样的字体显示路径称或目录名称。Menuchoice“ Menu choice”这样的体用来突出菜单路径,如引导选择Fie>openDesign单击迅速地按动并松开鼠标键按键并保持按鼠标键,并不松开。双击连续两次快速按动松开鼠标键。北京才略科技有限公司TEL:010)82673952/82673953SABER电气系统培训手册Saber Electrical SystemsLab#1-DC分析在第一个练习,您将对DC(工作点)分析进行基础的了解。这是你在大部分设计中将要作的基本分析。打开RLC设计UNIX用户1.找到 Saber_ Training_ Files/Saber_ Electrical_Training/ Feature Labs目录2.输入: sketchWindows用户:1. itf Start >Programs Synopsys >Saben)XX> Saber Sketch所有用户(从 Saber Sketch下拉菜单中1.选择Fie>Open> Design2左键单击并按住鼠标在菜单上打开下拉萊单。滑动光标选择相应子菜单并展开。点击下拉菜单中部可固定蕖单,即使松开鼠标键,菜单也不会弹回。3.浏览RLC自录(在 Saber_Training Files/Saber_ Electrical_TrainingFeature Labs路径下)。4.双击 Open Design中的exrl文件名(如果文件扩展名可见,选择带有 ai sch扩展名的文件)。RLC示意图显示如下:25mwimid Ivoutv pulse initial: 0ouselu北京才略科技有限公司TEL:010)82673952/82673953SABER电气系统培训手册Lab#1-DC分析从此处,您开始在 Saber Sketch中设计。5.点击Show/ Hide Saber Guide按键,该按键使您可以进λ Saber guide仿真,在整个设计过程中都将用到该按键。(按键位于Saber sketch图标栏的右侧)。Saber guide图标栏出现。工作点分析1.点击 Operating Point按键,幸2.点击OK,接受默认。生成列表,确定工作点(在生成列表前如提示是否保存,回答yes)。如有错误,将有信息提示。点击 Simulation Transcript按键 Smdl,打开 Saber Guide transcrip窗口您可以通过 Saber guide Transcript窗口监控 Saber guide命令进程。也可显示完成一次分析的执行时间。3.当分析完成,从 Saber sketch下拉菜单中选择Results> Operating Point Report2.显示分析结果。4.单击OK接受工作点报告缺省值Report Tool弹出,显示分析结果。注意所有的显示值都为0。检查该结果是否正确,查看驱动滤波器源电压的初始值。该示意图显示电压初始值为0,脉冲值为1。表示在0时刻电源电压为0,所以该结果正确为得到DC分析的非零值,您可以改变电源的初始值,例如,您可以将初始值设为1,脉冲值置为0。这样您将得到棉反的波形改变输人电压并重新分析1.改变示例中电源的初始电压a,在 initial:0处的0附近单击左键b.通过箭头键将光标置于0的右侧(如果需要)c,单击删除键删除0,在该位置键入1d.单击 Return{ Enter)或鼠标左键2.改变示例中的脉冲值a.在puse:1处的0附近单击左键b.通过箭头键将光标置于1的右侧(如果需要)c.单击删除键删除1,在该位置键入0d.单击 Return( Enter)或鼠标左键北京才略科技有限公司TEL:010)82673952/82673953SABER电气系统培训手册Saber Electrical Systems您已完成下面两项工作改变了示例中电源的初始电压及脉冲值。当您编辑列表并输入 Saber时,这些值会随之自动动态改变,这就意味着模拟器已经接收了新值,您不需要保存或重新生成网表(您可以在 Saber Guide Transcript窗口中看到这些值的变化)。现在重新启动DC分析,除了相应的变化,您可以看到 Saber执行分析时有很多选项。3.从 Saber Guide下拉菜单栏中,选择Analyses Operating Point DC Operating Point.*4.在 Operating Point分析中,分析结束后,执行以下操作,工作点报告可以自动显示在Saber guide transcriptt窗口内a.点击 Display After Analysis旁边的Yes按钮b.点击OKc.在 Saber Guide Transcript窗口内查看结果结果显示1V的输入电压产生0909V输出电压。你理解这个结果吗?为什么?为了便于结合给定电压与实际设计示意图进行分析,可按如下方法在 Saber sketch i中显示DC分析结果:标注分析结果1.从 Saber Sketch下拉菜单栏中选择 Results> Back Annotation2.在后注释栏点击OK注意示例中仿真电压如何出现关闭Lab#11.按以下步骤关闭报告及报告工具栏a.选择File>Cose关闭报告,关闭前不保存 de. rpt文件b,.选择File> Close Window关闭报告工具2.恢复电压源初始值,为下一练习作准备(即将初始值还原为0,脉冲还原为1)请告知培训人员您已经完成了第一部分的内容。北京才略科技有限公司TEL:010)82673952/82673953SABER电气系统培训手册Saber Electrical SystemsLab#2一时域分析在本节,您将对RLC滤波器进行时域(瞬态)分析,以确定其脉冲输入响应。继续采用您上一节建立的RLC电路。瞬态分析1.点击 Transient Analysis按钮,2.在时域瞬态分析中,输入数值:a. End Time: Jomb. Time Step: 0. Iu选择该值是因为脉冲升降时间是1毫秒,初始时间长应置为110脉冲时间)Run DC Analysis First: Yesd. Plot After Analysis: Yes-Open Onl3.点击OK执行分析完成瞬态分析, OsmoscOpe.信号管理器和图形窗们命自动打开Cosmos scope的使用当完成分析后, OsmoscOpe自动打开(因为分析中您选择了Yes. Open Only to Plot AfterAnalysis)。同样生成了两个窗口,如下所示信号管理器 Signal Manager:信号管理器显示当前激活状态下了图形文件名称(图形文件包含仿真数据在 OsmoscOpe中可见)。gnal ManagerFe Polle SionalsSignal FiterOpen ClothesPlotsClose Plotfdles( 1]ex die. ac ai plDisplay PottiesSetupMatch All图形文件窗口 Plot file window:图形文件窗口显示相应图形文件中的信号名称。北京才略科技有限公司TEL:010)82673952/82673953
- 2020-12-03下载
- 积分:1