局部特征尺度分解方法及其分解能力研究_杨宇
一种新的时频分析方法——局部特征尺度分解
2 O 2 年 1
湖
南
大
学
学
报 (自 然 科 学 版 )
Vo _ 9, . l 3 No 6
6 月
J u n l f Hu a ie st ( t r lS i n e ) o r a n n Un v r i Na u a c e c s o y
Th na y i e u t v e e a l ss r s ls ha e d mons r t d t e v ld t h wo d c t a e h a i iy oft e t e ompo ii n m e ho . M or ov r sto t ds e e ,LCD s i s upe i o EM D n c rort i ompu a i n le fc e y a d r s rc i n nd e f c s n a ii t to a fi inc n e t i to ofe fe t .I dd ton,LCD sa s p id i l o a ple t a t d a no i o o lr be rng a h na y i e ulsf o t c ua a tvi a i n sgn lha e f — o f ul i g s sf r r le a i nd t e a l s sr s t r m he a t lf ul br to i a v ur t e r e h fe tv ne s o h r p ov d t e e f c i e s fLCD. Ke r s: a t i gn s s l a c r c e itc s a e d c y wo d f ul d a o i ; oc l ha a t rs i— c l e ompo ii n; i t i i c l c sto n rnsc s a e omp e ;r lr on nt ole be rn a i gs
自适应特征尺度分解方法及其应用
方法都 具 有 一 定 的 局 限 性 。 局 部 特 征 尺 度 分 解 ( , l o c a l c h a r a c t e r i s t i c- s c a l e d e c o m o s i t i o n L C D) p [ 方法是最近提 出 的 一 种 自 适 应 时 频 分 析 方 法 4]。 该方法在定义瞬时频率具有物理意义的单分量信 号— — — 内禀尺度分 量 ( , i n t r i n s i c s c a l e c o m o n e n t p 基 础 上, 通过提取数据本身的特征尺度参 I S C) 数, 可自适应地将复 杂 信 号 分 解 为 若 干 个 I S C分 量之和 。 相 对 于 EMD 和 LMD 方 法 , L C D 方法 抑制 端 点 效 应 和 模 态 混 叠 等 方 面 具 在分解速度 、
[ 1 0] 方法和 A MD) C L C D 方 法。 对 具 有 碰 摩 故 障 的转子振动位移 信 号 的 分 析 结 果 验 证 了 A C D的
·3 1 8 4·
中国机械工程第 2 6 卷第 2 3期2 0 1 5年1 2 月上半月
自适应特征尺度分解方法及其应用
吴占涛 程军圣 杨 宇
湖南大学汽车车身先进设计制造国家重点实验室 , 长沙 , 4 1 0 0 8 2
摘要 : 针对局部特征尺度分解 ( 存在的模态混叠问题和其在均值曲线定义方面存在的不足 , 在 L C D) 充分借鉴经验模态分解 ( 和L 对L C D 方法研究的基础上 , EMD) C D 等此类基于筛分的信号分解方法的 定义了一种新的瞬时频率具有物理意义的单分量信号 — — — 内禀致密尺度分量 ( , 并提出了一 思路 , I C C) 种新的自适应信号分解方法 — — —自 适 应 特 征 尺 度 分 解( 方 法。同 时, 给 出 了I A C D) C C 分 量 评 价 准 则, 通过对 A 选取 C D 每阶筛分中由不同均值曲 线 和 致 密 系 数 取 值 得 到 的 一 组 不 同 的 分 解 分 量 进 行 对 比 , 从而保障最终分解效果优于 L 最优分量作为该阶筛分的 I C C 分量 , C D 方 法 分 解 效 果。 对 仿 真 信 号 的 分析结果证实了 A 方法的分解效果优于 、 、 总体平均经验模态分解 ( 和自主致密局 C D EMD L C D E EMD) 部特征尺度分解 ( 方法的分解效果 ; 对实验数据的分析结果验证了 A 从而为旋转 A C L C D) C D 的有效性 , 机械故障诊断提供了一种新的方法 。 关键词 : 自适应特征尺度分解 ; 局部特征尺度分解 ; 经验模态分解 ; 内禀致密尺度分量 ; 故障诊断 中图分类号 : : / TH 1 6 5; TH 9 1 1. 7 DO I 1 0. 3 9 6 9 . i s s n. 1 0 0 4 - 1 3 2 X. 2 0 1 5. 2 3. 0 1 2 j A d a t i v e C h a r a c t e r i s t i c s c a l e D e c o m o s i t i o n M e t h o d a n d I t s A l i c a t i o n s - p p p p Wu Z h a n t a o C h e n J u n s h e n Y a n Y u g g g , k e L a b o r a t o r o f A d v a n c e d D e s i n a n d M a n u f a c t u r e f o r V e h i c l e B o d S t a t e y y g y , , H u n a n U n i v e r s i t C h a n s h a 4 1 0 0 8 2 y g : , A b s t r a c t a e r B a s e d o n t h e r e s e a r c h o f L C D t h i s b o r r o w e d i d e a s f r o m t h e s i f t i n b a s e d s i n a l p p g g m e t h o d s s u c h a s EMD a n d L C D, a n d t h e n a n e w m o n o - c o m o n e n t s i n a l w i t h d e c o m o s i t i o n r e s e n t e d p g p p , , h s i c a l l r o b l e m m e a n i n f u l i n s t a n t a n e o u s f r e u e n c i e s i . e . I C C, f o r r e s t r a i n i n t h e m o d e m i x i n p y y p g q g g , o f L C D a n d i m r o v i n t h e d e f i c i e n c i n t e r m s o f t h e m e a n c u r v e d e f i n i t i o n. F i n a l l a n e w s i n a l d e p g y y g - , , r o o s e d . M e a n w h i l et i v e n. I n c o m o s i t i o n m e t h o dA C D w a s h e e v a l u a t i o n c r i t e r i o n o f I C C w a s a l s o p p g p , t h e s i f t i n f o r s e a r a t i n c e r t a i n o r d e r c o m o n e n t s a s e t o f I C C s w a s o b t a i n e d b u s i n d i f r o c e d u r e g p g p y g p - f e r e n t m e a n c u r v e s a n d c o m a c t c o e f f i c i e n t s . T h e o t i m a l I C C f o r t h i s o r d e r s i f t i n w o u l d b e s e l e c t e d p p g u a r a n t e e d t h e c a n d i d a t e I C C s u s i n t h e e v a l u a t i o n c r i t e r i o n o f I C C,w h i c h A C D o u t e r f o r m s f r o m g g p L C D. T h e s i m u l a t i o n r e s u l t s i n d i c a t e t h a t t h e d e c o m o s i t i o n e f f e c t o f A C D i s b e t t e r t h a n t h a t o f EMD, p L C D, e n s e m b l e e m i r i c a l m o d e d e c o m o s i t i o n( E EMD) a n d a u t o n o m o u s c o m a c t l o c a l c h a r a c t e r i s t i c - p p p ( ) , s c a l e d e c o m o s i t i o nA C L C D e x e r i m e n t a l r e s u l t s o f t h e e x e r i m e n t a l s i n a l s s h o w i t s v a l i d i t . p p p g y r o v i d e d . T h u s a n e w w a f o r f a u l t d i a n o s i s o f r o t a t i n m a c h i n e r i s p y g g y : ; K e w o r d s a d a t i v e c h a r a c t e r i s t i c - s c a l e d e c o m o s i t i o n( A C D) l o c a l c h a r a c t e r i s t i c - s c a l e d e c o m o y p p p - ( ; ( ; ( ; s i t i o n L C D) e m i r i c a l m o d e d e c o m o s i t i o n EMD) i n t r i n s i c c o m a c t - s c a l e c o m o n e n t I C C) f a u l t d i p p p p - a n o s i s g
基于LCD和局部Hilbert边际能量谱的滚动轴承故障特征提取方法
基于LCD和局部Hilbert边际能量谱的滚动轴承故障特征提取方法杨潘;李东文;曾鸣【期刊名称】《机械强度》【年(卷),期】2015(37)1【摘要】在定义局部Hilbert边际能量谱的基础上,提出了一种基于局部特征尺度分解(Local characteristic-scale decomposition,LCD)和局部Hilbert边际能量谱的滚动轴承故障特征提取方法。
采用LCD方法对滚动轴承原始振动信号进行分解得到若干内禀尺度分量(Intrinsic scale component,ISC),然后对各个ISC分量进行Hilbert解调得到信号的时频分布。
根据信号时频分布中能量分布确定频率段的下限和上限频率,从而得到相应的局部Hilbert边际能量谱,计算该频率段内信号的能量并将其作为故障特征参数。
实验分析结果表明,所提出的方法能有效地提取滚动轴承故障特征信息。
【总页数】4页(P9-12)【关键词】局部特征尺度分解;局部Hilbert边际;能量谱;能量;滚动轴承;故障诊断【作者】杨潘;李东文;曾鸣【作者单位】湖南大学机械与运载工程学院;湖南大学汽车车身先进设计制造国家重点实验室【正文语种】中文【中图分类】TH165.3【相关文献】1.局部Hilbert边际能量谱在滚动轴承故障诊断中的应用 [J], 杨恭勇;周小龙;李家飞;梁秀霞2.基于LCD-Hilbert谱奇异值和QRVPMCD的滚动轴承故障诊断方法 [J], 杨宇;何知义;潘海洋;程军圣3.基于Hilbert边际谱的滚动轴承故障诊断方法 [J], 杨宇;于德介;程军圣4.基于广义形态学滤波和 Hilbert边际谱的滚动轴承故障诊断 [J], 崔宝珍;马泽玮;李会龙;王珊5.基于Hilbert边际谱和IPSO-SVDD的滚动轴承故障诊断 [J], 高峰;申江江;曲建岭;袁涛;何勰;余路因版权原因,仅展示原文概要,查看原文内容请购买。
总体局部特征尺度分解及ELM的滚动轴承故障诊断
量( I n t r i n s i c S e M e C o mp o n e n t , I S C ) ; 其次, 根据分解后 I S C分量计算时域指标、 能量、 相对熵, 利用特征评估法提取敏感特征 ; 最后, 将敏感特征向量输入极限学习机 ( E x t r e m e L e a r n i n g Ma c h i n e , E L M) 进行训练与测试 , 从而识别滚动轴承的故障类型。 对 实验信号的分析表明, 该方法能够有效的诊 断出滚动轴承不 同的工作状 态, 且效果较局部特征尺度分解方法好 。
关键词 : 滚动轴承 ; 故障诊断 ; 极 限学 习机 ; 总体局部特征尺度分解 ; 特征评估法 中图分类号 : T Hl 6 文献标识码 : A 文章编号 : 1 0 0 1 — 3 9 9 7 ( 2 0 1 7 ) 0 2 — 0 2 2 6 — 0 5
F a u l t Di a g n o s i s o f B a l l Be a r i n g B a s e d o n EL CD a n d EL M
机 械 设 计 与 制 造
22 6 Ma c h i n e r y De s i g n & Ma n u f a c t ur e
第 2期 2 0 1 7年 2月
总体局部特征尺度分解及 E L M 的滚动轴承故 障诊 断
董素鸽 , 胡代 弟, 葛明涛
( 郑州大学 西亚斯 国际学院 电子工程系 , 河南 郑州 4 5 1 1 5 0 )
A b s t r a c t : Ai mi n g讲t h e n o s t ti a o n a r y c h a r a c t e r i s t i c o fag e a r f a u l t v i b r t a i o n s i g n a 1 . i t p r o p o s e s a r e c o g n i t i o n m e t h o d b me d o n
基于EMD的奇异值分解技术在滚动轴承故障诊断中的应用_杨宇
式中 , 当 n 为偶数时 , J = n /2, 当 n 为奇数时 , J = (n + 1) /2。 这样 , 通过 EMD 分解 , 滚动轴承振动信号 x ( t) 的特征就可以由初始特征向量矩阵 A 和 B 来刻画 , 对 初始特征向量矩阵 A 和 B 进行奇异值分解 , 就可以提 取滚动轴承振动信号的特征 。 本文提出的滚动轴承故障诊断方法流程图如图 3 所示 。
n
(3) 式可得到 : x (t) = ∑ ci (t) + rn (t) i =1
(4 )
式中 , rn (t) 为残余函数 , 代表信号的平均趋势 。 而各 M F 分量 c1 (t), c2 (t)… cn (t) I 分别包含了信号从高到 低不同频率段的成分 , 每一频率段所包含的频率成分 都是不 同的 , 且随 信号 本 身的 变 化而 变 化 。 通 常 , EMD 方法分解出来的前几个 I M F 分量集中了原信号 中主要的信息 , 从这个角度上讲 , EMD 方法就是一种 新的主成分分析方法 。 图 1 为滚动轴承 外圈有缺陷时的时域加速度振 动信号 。 图 2 为其 E M D 分解图 , 共有 22 个 I MF 分 量 , 由于篇幅原因 , 此处只给出了前 10 个 I M F 分量 。 从图中可以 看出 , EMD 把信号 分解成了若干 个 I MF 分量之和 , 不同的 I M F 分量包含了不同的时间尺度 , 可以使信号的特征在不同的 分辨率下显示出来 。 从 图 2 中还可以看出 , E M D 分解结果的前几个 I M F 分量 包含了原始信号的主要信息 。
DO I 牶 牨 牥 牣 牨 牫 牬 牰 牭 牤 j牣 cnki 牣 jvs牣 牪 牥 牥 牭 牣 牥 牪 牣 牥 牥 牬
基于LMD和AR模型的转子系统故障诊断方法
Fa tDig o i pr a h f r Ro o s e s ul a n ss Ap o c o t rSy t m
Ba e n LM D nd A R o e sdo a M dl
YANGYu ,YANG i in L — a g。CHENG u — h n x J ns e g
f tt y L b rtr fAd a cdDein a dMa ua tr o hceB d 。Hu a i,C ag h Hu a 4 0 8 , ia S aeKe a o aoyo vn e s n n fcu efrVe il o y g n nUnv h n sa, n n 10 2 Chn )
s o h t b h EM D n h wn t a ot a d LM D e h a pp i d o t e r t r s t m a t d a no i e f c i ey m t od c n be a le t h o o ys e f ul i g ss fe tv l .
文 章编 号 : 6 4 2 7 ( 0 0 0 — 0 4 0 1 7 — 9 4 2 1 ) 90 2 5
基 于 LM D 和 AR 模 型 的 转 子 系 统 故 障 诊 断 :- h法
杨 宇 , 丽 湘 , 军 圣 杨 程
( 湖南 大 学 汽 车 车身 先 进 设 计制 造 国家 重 点 实 验 室 , 南 长沙 湖 408) 1 0 2
a nd AR mod l sp op s d e wa r o e .Fis l r ty,by u i D t o sng IM me h d,t e v b a i i a fr o y t mswa — h i r ton sgn lo ot rs s e sde
一种新的非平稳信号分析方法——局部特征尺度分解法
L MD方 法 将一个 单分 量 的调 幅一 频信 号看 成 是其 调 本 身 的 包 络 信 号 和 一 个 纯 调 频 信 号 的乘 积 , P 即 F
( r d c u cin 简称 P ) 量 , P ou t n t , f o F分 因此 能 自适 应 地 将一 个复 杂信 号分 解 为若 干 个瞬 时频 率具 有物 理意 义 的P F分 量之 和 。 MD方法 虽然 避免 了E L MD 方法 中的过 包 络 、 欠包 络 以及 由希 尔 伯 特变 换 而 产 生 的 负 频 率等 问题 , L 但 MD 本 身也 存在 迭代 计 算 量大 、 频 率混 淆 、 点 效应 等 问题 。2 0 端 0 6年 , r ri Mak G F e 提 出 了另 外一 种 自适应 时频 分析 方法 —— 本 征 时间 尺 度 分 解 方 法 (nr s i —c l d c mp st n, Iti i t n c mesae eo o io i
波 基 , 乏 自适 应 性L 。 缺 2 ] 19 9 8年 , 国 国 家 宇 航 局 的 Nod n E Hu n 美 re a g
等 创造性 地 提 出了处理 非平 稳信 号 的新方 法—— 希 尔 伯特 黄 变换 方法 ( let a gta so m, 称 Hi r Hu n rn fr 简 b — HHT) 州] [ 。该 方 法 能 自适 应 地 将 一个 复 杂 信 号分 3
引 言
经典 的傅 里 叶变换 方法 只能 处理 线性 和平 稳信 号 , 自然界 中的大部 分信 号是 非线 性 和非平 稳 的 。 而
S t mi h提 出 了另一 种新 的 时频 分析 方法 —— 局 部均 值分 解 ( o a me nd cmp s in 简 称 L L cl a eo o i o , t MD)1 [ ,
一种基于FastText_的恶意代码家族分类方法
word2vec 首先 生 成 词 典ꎬ 再 将 此 词 典 中 的 所 有
API 函数映射为词向量ꎮ 所生成的词向量仅是样
本语义到向量的简单映射ꎬ局限于 API 函数的顺
序特征且依赖词典ꎬ缺乏对词典外的 API 函数表
征能力ꎮ 恶意代码随着混淆技术的发展ꎬ执行程
98% ꎮ 王栋等
[9]
基于 VGG 模型构建了一维卷积
分类网络模型 ID ̄CNN ̄IMIRꎬ相较于其他深度学
习模型具有更好的性能ꎬ准确率可达 98 94% ꎮ 静
态分析的缺点是容易受到混淆的影响ꎬ可直接分
析原始程序代码的语法和语义ꎬ进而提取到更深
层次的特征ꎮ
在 针 对 API 序 列 的 静 态 检 测 中ꎬ 常 用
指具有相似功能、相同来源和不同演进程度的恶
意代码集合
[2]
ꎬ其分类研究可以视为对不同恶意
在恶意代码家族分类中ꎬ受文献[10] 启发ꎬ本
文对收集的 API 数据进行向量化ꎮ 由于卷积神经
代码的相似性、关联性判断ꎬ有助于研究人员发现
网络( CNN) 和长短时记忆( LSTM) 网络在特征提
同类型的变异代码并快速了解恶意代码感染策
proposed. The FastText framework is utilized to achieve multidimensional vector transformation and
precise expression of code sample API sequencesꎬand one ̄dimensional convolution and long short
基于三次Hermite插值的局部特征尺度分解方法
E-mail: pansea1989@ 通讯作者:潘孟春(1963-),男,教授。
析有窗口 Fourier 变换、Wigner-Wille 分布和小波变 换 ,但上述方法都存在一定的缺陷 ,如窗口 Fourier 变换具有一定的时频窗口,Wigner-Wille 分布有交叉 项的干扰,而小波变换也需要小波基的选择,缺乏自 适应 。 [1] 因此,传统的信号分析方法在应用方面具 有很大的局限性。
Abstract : Since linear transformation is used to obtain baseline signal in intrinsic time scale decomposition (ITD) method, burr and instantaneous frequency distortion will appear in the decomposition results. Therefore, a rational Cubic Hermite interpolation—Local characteristic- scale decomposition (CHLCD) method is presented. In this method, any complex signal can be adaptively decomposed into a sum of several independent rational intrinsic scale components (ISCs), whose instantaneous frequencies have obvious physical meanings. Firstly, the principle of the CHLCD method was analyzed. Then, the detailed steps of CHLCD of signal were given. Finally, a simulation signal was adopted to verify the CHLCD method. Experimental results show that the CHLCD method can effectively decompose signals.
基于局部特征尺度分解和核最近邻凸包分类算法的滚动轴承故障诊断方法
级 分类器 将其 转化 为 多 个 两类 问题 来 处 理 , 而 采用 不 同 的分 类器 结构 ( 如 一对 一 、 一对多、 有 向无 环 图 和 二叉树 结 构 等 ) 也会对 S VM 的分 类 性 能 有 所 影
工
类 问题 且 需 优 化 的参 数 只有 核 参 数 。实 验 分 析 结 果 表 明 , 所 提 出 的方 法 能 有 效 地 提 取 滚 动 轴 承 故 障 特 征 信 息 , 而
程E 叫
学
且 在 小 样 本 的情 况 下 仍 能 准 确 地 对 滚 动 轴 承 的工 作 状 态 和 故 障 类 型 进 行 分 类 。 同 时 , 与支持向量机 ( S u p p o r t v e c ~
t o r ma c h i n e , S VM ) 算 法 的 对 比分 析 结 果 表 明 , KN NC H 算 法 的 分 类 性 能 的 稳 定 性要 高 于 S V M 算法 。
关键词 : 滚动轴承 ; 故 障诊 断 ;局 部 特 征 尺 度 分解 ; 核 最 近 邻 凸包 分 类 算 法 ; 奇 异 值 分 解
响l _ 3 ] 。另 外 , S VM 的分类 性 能还受核 参数 以及 惩 罚 系 数 的影 响 , 需要优 化 的参数 较 多 。针 对这些 问题 ,
种新 的 自适应 时频 分析 方法 _ 6 ] 。该 方法采 用 I T D
中基 线信 号 的构 造 方 法 , 即 以 原始 信 号 任 意 两个 相 邻 的极值 点为跨 度对 信号 进行分 段线 性变 换来 构造
摘要 : 提 出 了一 种 基 于 局 部 特 征 尺 度 分 解 ( L o c a l c h a r a c t e r i s t i c — s c a l e d e c o mp o s i t i o n , L C D) 和核最近 邻凸包 ( Ke r n e l
一种新的非平稳信号分析方法_局部特征尺度分解法_程军圣
一种新的非平稳信号分析方法——局部特征尺度分解法程军圣,郑近德,杨 宇(湖南大学汽车车身先进设计制造国家重点实验室,湖南长沙410082)摘要:在定义瞬时频率具有物理意义的内禀尺度分量(Intrinsic scale component,简称ISC)的基础上,提出了一种新的非平稳信号分析方法——局部特征尺度分解方法(Local char acteristic-scale deco mpo sition,简称L CD),该方法能够自适应地将一个复杂信号分解为若干个瞬时频率具有物理意义的内禀尺度分量之和。
首先对LCD方法的原理进行了分析,然后给出了采用L CD对信号进行分解的详细步骤,最后采用仿真信号对LCD和经验模态分解(Empirica l mode decomposition,简称EM D)方法进行了对比分析,结果表明了LCD方法的有效性及在端点效应、迭代次数和分解时间等方面都优于EM D分解方法。
关键词:局部特征尺度分解;经验模态分解;内禀尺度分量;非平稳信号中图分类号:T N911.7;T H165.3 文献标识码:A 文章编号:1004-4523(2012)02-0215-06引 言经典的傅里叶变换方法只能处理线性和平稳信号,而自然界中的大部分信号是非线性和非平稳的。
由于时频分析方法能同时提供非平稳信号在时域和频域的局部化信息而得到了广泛的应用。
典型的时频分析方法有窗口傅里叶变换、Wigner-Wille分布、小波变换等[1]。
但这些方法都有各自的缺点,如窗口傅里叶变换具有固定的时频窗口,Wigner-Wille分布存在交叉项干扰,而小波变换则需要事先选择小波基,缺乏自适应性[2]。
1998年,美国国家宇航局的No rden E Huang 等创造性地提出了处理非平稳信号的新方法——希尔伯特黄变换方法(H ilbert-Hua ng transform,简称H HT)[3~6]。
该方法能自适应地将一个复杂信号分解为若干个内禀模态函数(Intrinsic mode functio n,简称IM F)之和,并进一步对各个IM F分量进行Hilbert变换,得到瞬时幅值和瞬时频率。
一种基于改进辛几何模态分解的复合故障诊断方法
一种基于改进辛几何模态分解的复合故障诊断方法作者:杨宇程健彭晓燕潘海洋程军圣来源:《湖南大学学报·自然科学版》2020年第02期摘要:针对辛几何模态分解方法分析结果的不确定性,提出一种改进的辛几何模态分解方法. 首先对原时间序列进行相空间变换,获得轨迹矩阵;然后通过辛几何相似变换求得特征值和对应的特征向量,并通过对角平均得到一系列的初始辛几何分量;最后采用层次聚类方法对初始辛几何分量进行自适应重组,进而得到最终的聚类辛几何分量. 实验结果表明:改进的辛几何模态分解方法可以有效地对旋转机械复合故障信号进行特征提取,提高故障诊断的准确性.关键词:层次聚类;改进的辛几何模态分解;故障诊断;旋转机械;信号处理中图分类号:TH113 文献标志码:AAbstract:Aiming at the uncertainty of Symplectic Geometry Mode Decomposition (SGMD)method,an Improved Symplectic Geometry Mode Decomposition (ISGMD) method was proposed. Firstly,the phase space transformation of the original time series was carried out to obtain the trajectory matrix. Then the eigenvalues and the corresponding eigenvectors were obtained by symplectic geometric similarity transformation,and a series of initial symplectic geometric components were obtained by diagonal average. Finally,the hierarchical cluster method was employed to adaptively restructure the initial symplectic geometry components to obtain the final Cluster Symplectic Geometry Component (CSGC). The experimental results indicate that the ISGMD method can effectively extract features from composite fault signals of rotating machinery and improve the accuracy of fault diagnosis.Key words:hierarchical cluster;Improved Symplectic Geometry Mode Decomposition (ISGMD);fault diagnosis;rotating machinery;signal processing齿轮、滚动轴承等零部件作为机械设备的核心元素,其状态好坏直接关系到整个机械设备的运转性能. 当前,针对这种现象学者们已经开展研究,并取得了一定的成果. 但是大多学者只针对某单一故障进行研究,却忽视了复合故障存在的可能性. 当机械设备出现复合故障时,多种故障信号相互耦合以及在强噪声的干扰下使得故障类型分析难度加大[1]. 因此,需要采用合适的信号分析方法提取出有效的故障信息进行分析.经验模态分解(Empirical Mode Decomposition,EMD)[2-3]方法作为一种典型的信号分析方法,可以将待分析信号分解为若干个有意义的内禀模态函数(Intrinsic Mode Function,IMF)之和. IMF分量的频率受采样频率的影响,同时也随着原始信号的变化而变化,通过迭代和极值点包络可以有效地从待诊断信号中提取出含有故障信息的分量信號. 因此EMD方法是一种自适应的信号分析方法,适合用于处理非平稳信号. 局部特征尺度分解(Local Characteristic-scale Decomposition,LCD)[4-5]方法是近几年提出的一种新方法,该方法借鉴EMD的思想通过设置基线信号将信号分解成多个内禀尺度分量(Intrinsic Scale Component,ISC)之和. 与EMD方法相比,该方法运算速度大大提高,信号处理效率得到明显提升,适用于在线监测. 虽然EMD和LCD在处理非平稳信号以及运算效率方面具有一定的优势,但它们同样也存在着许多难以克服的缺点,如端点效应和模态混叠等. 对此,相关学者已经做出大量研究工作,并提出相应的改进方法,抑制了端点效应、模态混叠等现象,但是又带来了新的不足. 例如提出的集合模态分解(Ensemble Empirical Mode Decomposition,EEMD)方法一定程度上抑制了模态混叠现象,但增加了计算量,同时无法完全中和添加的白噪声,使EEMD不具有完备性.近年来,随着辛几何谱分析(Symplectic Geometry Spectrum Analysis,SGSA)[6-8]方法的提出,许多学者对其做了大量研究. 该方法是一种基于辛几何的信号分析方法,其核心是通过辛几何相似变换在不改变原时间序列本质特征的基础上获得相应的辛几何分量. SGSA方法具有较强的抗噪能力,在处理非平稳故障信号中有着较好的效果,但它同时也存在着自身的不足,主要有以下两点:第一,无法自适应地选择嵌入维数构造轨迹矩阵;第二,当进行辛几何重构时,该方法利用周期、频率等相似性来重构出最终的辛几何分量,并需要人为设置终止条件. 这些不足会导致最终的分析结果具有不确定性,严重影响辛几何谱分析方法的推广. 针对第一点不足,Pan等[9]将辛几何模态分解(Symplectic Geometry Mode Decomposition,SGMD)方法应用于故障诊断,采用功率谱密度(Power Spectral Density,PSD)方法自适应地确定嵌入维数,从而消除了人为设置嵌入维数给结果带来的不确定性. 但其辛几何分量的构造依然采用频率和相关系数相似度,同时也需要人为地设置迭代终止条件. 该缺陷将导致SGMD方法所得到的最终分析结果依然存在着一定的不确定性. 针对SGMD方法的这一缺陷,本文提出了一种改进的辛几何模态分解(Improved Symplectic Geometry Mode Decomposition,ISGMD)方法,该方法首先对待分析的时间序列信号进行相空间变换,获得轨迹矩阵;其次将轨迹矩阵进行辛几何相似变换求出特征值,所求特征值对应的特征向量经过重构得到具有原始信号完整结构信息的辛几何分量矩阵;然后对辛几何分量矩阵做对角平均化得到一系列的初始辛几何分量;最后采用层次聚类方法对初始辛几何分量进行自适应重组,进而得到最终的聚类辛几何分量(Cluster Symplectic Geometry Component,CSGC).本文利用層次聚类(Hierarchical Cluster)方法[10-11]不需要设置聚类的阈值就可以得到较为准确的聚类结果这一优点,将层次聚类方法引入到SGMD方法中,其目的是为了实现分量自适应地重组,提高分析结果的准确性. 基于ISGMD是一种非平稳信号处理方法,同时可以将复合故障信号分解成若干个含有不同频率信息的辛几何分量的特点,本文的ISGMD方法可用于旋转机械复合故障诊断.从图19~图22可以看出,图19、图20、图21中提取出了轴承故障特征频率,而齿轮故障特征频率被完全覆盖;图22中同时提取出了齿轮和轴承故障特征频率. 对比可知,ISGMD方法在复合故障诊断方面优于EMD、LCD和SGMD方法,ISGMD方法是一种较好的复合故障诊断方法.4 结论本文提出了一种改进的辛几何模态分解方法,采用了层次聚类方法对SGMD方法进行改进,使其能够自适应地进行分量重组,并将其应用于复合故障诊断中. 通过对复合故障仿真信号和实验信号分析并将ISGMD方法和EMD、LCD、SGMD方法对比,结果表明ISGMD方法能更好地识别复合故障.参考文献[1] 杨宇,罗鹏,程军圣. 广义变分模态分解及其在齿轮箱复合故障诊断中的应用[J]. 中国机械工程,2017,28(9):1069—1073.YANG Y,LUO P,CHENG J S. GVMD and its applications in composite fault diagnosis for gearboxes.[J]. China Mechanical Engineering,2017,28 (9):1069—1073. (In Chinese)[2] YU D J,YANG Y,CHENG J S. Application of time-frequency entropy method based on Hilbert-Huang transform to gear fault diagnosis[J]. Measurement,2007,40 (9):823—830.[3] YANG Z H,YU Z H,CHAO X,et al. Application of Hilbert-Huang transform to acoustic emission signal for burn feature extraction in surface grinding process[J]. Measurement,2014,47:14—21.[4] ZHENG J D,CHENG J S,YANG Y. A rolling bearing fault diagnosis approach based on LCD and fuzzy entropy[J]. Mechanism and Machine Theory,2013,70:441—447.[5] LUO S R,CHENG J S.VPMCD based novelty detection method on and its application to fault identification for local characteristic-scale decomposition[J]. Cluster Computing,2017,20 (4):2955—2965.[6] XIE H B,DOKOS S,SIVAKUMAR B,et al. Symplectic geometry spectrum regression for prediction of noisy time series[J]. Physical Review E,2016,93(1/5):052217.[7] LEI M,WANG Z Z,FENG Z J. A method of embedding dimension estimation based on symplectic geometry[J]. Physics Letters A,2002,303(2):179—189.[8] XIE H B,DOKOS S. A hybrid symplectic geometry spectrum analysis and central tendency measure method for detection of determinism in noisy time series with application to mechanomyography[J]. Chaos,2013,23:023131.[9] PAN H Y,YANG Y,LI X,et al. Symplectic geometry mode decomposition and its application to rotating machinery compound fault diagnosis[J]. Mechanical Systems and Signal Processing,2019,114:192—194.[10] LIN Y C. Minimax problems under hierarchical structures[J]. Journal of Inequalities and Applications,2015(1):1—10.[11] CALÓ D G,MONTANARI A,VIROLI C. A hierarchical modeling approach for clustering probability density functions[J]. Computational Statistics & Data Analysis,2014,71(3):79—91.[12] BONIZZ P,KAREL M H,MESTE O,et al. Singular spectrum decomposition:a new time series decomposition [J]. Advances in Adaptive Data Analysis,2014,6(4):107—109.。
2015年湖南省优秀硕士学位论文名单
口腔鳞状细胞癌癌相关成纤维细胞中FAK 表达及其生物学行为的初 100302 口腔临床医学 步研究 长沙市H1N1流感大流行防控措施的动力学模型研究 乳腺癌靶向长春瑞滨纳米粒构建及体外评价 100400 公共卫生与预防医学 100702 药剂学
COX-2,ABO和NR3C2基因多态性对依那普利咳嗽不良反应及疗效的影 100706 药理学 响 6
作者姓名 邓金 张砚 王俏 韩知秋 彭星 张海 马志丽 郭静楠 刘卫 周江 李娟萍 陈兰 张锡霖 章运生 谢敏敏 赵应香 谭力鹏 李维逸 彭艳 张庆飞 肖文甲
论文名称
学科代码及名称
导师姓名 周平 胡军华 许民利 王坚强 李斌 杨招军 黄桂芳 赵敬哲 王青 聂舟 李继山 蔡青云 刘承斌 李丹 赵小英 何晓晓 张屹 经济法 030301 社会学 030301 社会学 070200 物理学 070302 分析化学 070305 高分子化学与物理 071000 生物学 071010 生物化学与分子生物学 080200 机械工程 080200 机械工程 080501 材料物理与化学 080502 材料学 080503 材料加工工程 080601 冶金物理化学 080603 有色金属冶金 080621 材料冶金 080703 动力机械及工程 080800 电气工程 080900 电子科学与技术 081000 信息与通信工程 081100 控制科学与工程
论文名称 基于电子地图的GPS/INS组合导航技术与算法研究 蛋白质网络的聚类分析方法及评估 节理岩体力学特征及隧道施工力学行为研究 高速列车荷载作用下交叉隧道动力响应特性及影响分区研究 联合地基GPS及空基COSMIC的对流层水汽三维层析 节理岩体采动损伤演化及破坏规律研究 中深孔爆破块度控制及测定方法研究 羟肟酸类捕收剂与细粒锡石的作用机理及浮选研究 视网膜血管增强方法研究 张量分解在大数据分析中的应用研究
基于分段多项式的局部特征尺度分解方法及应用
正交性 、精确性 等具有一 定优 越性 ;由转子碰摩故障诊断表明该方法的有效 性。
关键词 :局部特征尺度分解 ;分段 多项式 ;故障诊断 ;转子 ;内禀尺度分量
中 图分 类 号 :TH165+.3;TH911.7 文 献 标 志 码 :A
DOI:10.13465/j.cnki.jVS.2016.04.037
2.M odern Engineering Training Center, Hunan University,Chmagsha 410082,China;
3.School of Mechanical Engineering,Anhui University of Technology,Ma’anshan 243002,China)
is discontinuous,thereby the decomposition precision is lowered. Piecewise polynomial was used in PPLCD to replace line connection in LCD ,and the interpolation points on the mean curve were computed by the polynomial which was generated
第 35卷第 4期
振 动 与 冲 击
JOURNAL OF VIBRATION AND SH0CK
Vo J.35 No.4 2016
基于分段 多项式的局部特 征尺度分解方法及应 用
吴 占涛 ,程 军圣 ,曾 鸣 ,郑近德
(1.湖南 大学 汽车车身先进设计 制造 国家重点 实验 室 ,长沙 410082;2.湖南大学 现代工程训练 中心 ,长沙 3.安徽工业 大学 机械工程学院 ,安徽 马鞍山 243002)
一种局部特征尺度分解局部均值点求取方法
No.9Sep.2020第9期2020年9月组合机床与自动化加工技术Modular Machine Tool & Automatic Manufacturing Techninue文章编号:1001 -2265(2020)09 -0047 -05DOI : 10.13462//. .nki. mmtamt. 2020. 09. 011一种局部特征尺度分解局部均值点求取方法崔伟成,杨利斌,桑德一(海军航空大学岸防兵学院,山东烟台264001)摘要:局部特征尺度分解是一种基于信号时域局部特征的自适应分解方法,其核心在于准确估计局部均值点。
为保证局部均值点的估计精度,进而提高分解质量,给出了一种基于分段三次样条插值 的局部特征尺度分解均值点求取方法。
估计某一个极值点时刻对应的局部均值点时,首先,找出与 该极值点相邻的4个不同类型的极值点,基于三次样条方法插值出极值点之间的局部包络线;其 次,在包络线上求取极值点时刻对应的镜像极值,点;最后,把原始极值点与对应的镜像极值点取算 术平均值作为局部均值点。
通过对仿真信号和齿轮断齿振动信号的自适应分解发现了该方法能有效提高分解质量。
关键词:局部特征尺度分解;局部均值点;分段插值;三次样条插值中图分类号:TH16 ;TG506文献标识码:AA Computing Algorittm of Local Mean Point of Local Characteristic-scale Decomposition CUI Wei-cheng ,YANG Li-bin ,SANG De-yl(Colleae of Coastal Defense , Naval Aeronautical Uniwrsity , Yantai Shandong 264001, China )Abstract : The local characteristic -scale decomposition (LCD ) is an effective adaptive decompositionmehhod , whosefoundahion eiesin accuraheeshimahion ofhheeocaemean poinhs.In orderho improvehhedel composition quality of standard LCD , a computing algosthm of local mean points was given based on sub section polynomial mtemolaWon. When estimated the local mean point corresponding to an exheme point ,Firstly , the local mean envelope was estimated based on subsection cubic spline inWmolCion according to four adjacent dtferent types of extreme points. Secondly , the meTor extreme point was estimated on thelocal mean envelope. Finally , the mean of the micor point and the exheme point was calculated as the lo- caemean poinh.Theanaey&iofhhe&imueahed dahaand achuaegearing fauehdaha&how&hhahhhemehhod can improvehhequaeihy ofdecompo&ihion e f echiveey.Key words : local charactemstW-scale decomposition & local mean point & subsection inWmoStion ; cubic splineinherpoeahion0引言局部特征尺度分解(Local characte/syc-ocale de composition , LCD )[1-4]是湖南大学程军圣、杨宇等学者在本征时间尺度分解(IntOnsio ame-ocale decomposi tion , ITD ) [5]的基础上创造性地提出的一种有效的自 适应分解方法。
局部特征尺度分解方法在机械振动中的应用
丁i + 2 一丁
( t ) : = ∑/ / S C ( £ ) ( £ J… … … … . J……………… .( ( 5 )
2 L C D方 法 的 改进
上述方法采用线性变换所得到的基线信号不够
(
j ) … … …… ……
( 1 )
光滑,而采用三次样条插值虽会使基线信号具有很 好的光滑性 ,但由于被插值 函数序列的导数只在边
/
其 中的 t ∈( r _ ’ 丁 i + 。 ) 。
将各段 依次连接起来得到基线信号 。 ( 3 )把 日 从原信号 中分离 出来 ,得到剩余信
l。
( 3 )依次连接 日 k 得到基线信号 日。 。
将求 出的 日 代替原 L C D 中用 线性 变 换得 到 的
+ 1 . 7 s i n ( 5 . 5 7 r t ) s i n ( 6 0 7 r t ) … … …… ( 6 )
( 4 )如果 h ( t ) 满足 / S C分量条件 , 则/ S C 就 等于 h , ( t ) ;否则 ,把 h 。 ( t ) 作 为新 的原始信号 , 重复 上 面两个 步骤 ,直 到 h ( t) 满足 / S C分 量 的条 件为止 ,此时 h l k 【 ) 即为第一个 / S C 分量 。 / S C分 量应 满足 以下 条件 : ①在整个数据范 围内,相邻极值点 的符号互 异。这一条件可保证极值点之间的单调性 ,使得数 据 序列 具有单 一 的模态 。
信号再重复以上步骤 ,直到剩余量为一单调函数为
止。
最 后 ,原 时 间序列 被分解 为 T t 个/ S C分 量与 一
个残余量 肛 ( £ ) 之和 ,即
基于局部特征尺度分解的汽轮机故障诊断研究
法袁 在汽轮机实际运行中进行故障诊断遥 同时与 RBF 神经网络袁 BP 人工神经网络和支持向量机的 诊断结果进行了比较袁 结果表明袁 ELM 作为分类器 能准确地对汽轮机的故障进行诊断遥 1 局部特征尺度分解特征提取 1.1 经验模态分解 渊EMD冤
经验模式分解法 渊Empiri-cal Mode Decomposi原 tion袁 EMD冤 可根据振动信号自身的时间特征尺度 进行自适应分解袁 将轴承振动信号的状态信息分解 到不 同的内 蕴 模 式 分 量 渊Intrinsic Mode Function袁 IMF冤 中袁 从而为轴承状态深层次的挖掘奠定了基 础遥 EMD 最主要问题是模式混淆问题袁 它可使各 IMF 分量不能反映信号的真实内涵[4]遥 1.2 小波包变换
领域袁 但传统的神经网络袁 如 BP 神经网络由于采
用梯度下降网络学习算法而存在收敛速度慢袁 易陷
入局部最优等问题遥 又如径向基神经网络 渊RBF冤
是一种前向神经网络袁 其结构与多层前向网络的结
构相似袁 它能将一切问题的特征都变为数字袁 最终
结果会丢失信息[7]遥 再由 Vapnik 提出的支持向量机
小波包分解不仅对低频部分进行分解袁 也可对 高频部分实施分解袁 而且小波包分解能根据信号特 性和分析要求选择相应频带与信号频谱相匹配袁 是 一种比小波分解更为精细的分解方法[5]遥 1.3 局部特征尺度分解算法
LCD 是在本征时间尺度分解法的基础上提出 的袁 本征时间尺度分解法的平均曲线是基于信号本 身的线性变换袁 它明确了信号分量的波形失真遥
基于局部特征尺度分解法的声反馈抑制算法
基于局部特征尺度分解法的声反馈抑制算法叶磊;黄莹;孟木子【摘要】出于抑制扬声系统中声反馈的目的,提出基于新时频分析法——局部特征尺度分解(Local Characteristicscale Decomposition,LCD)的声反馈抑制新算法.通过一次分解,分离出局部特征最明显的内禀尺度分量(Intrinsic Scale Component,ISC)——啸叫单频分量,并以此算法为理论基础,给出了基于DSP的系统设计.对导入的语音信号叠加啸叫单音进行仿真分析,分析结果表明,此方法能有效分离语音与啸叫单音,并且可完好保留语音信息.【期刊名称】《电声技术》【年(卷),期】2015(039)008【总页数】5页(P54-57,62)【关键词】语音信号;局部特征尺度分解法;内禀尺度分量;DSP【作者】叶磊;黄莹;孟木子【作者单位】南京邮电大学电子科学与工程学院,江苏南京210000;南京邮电大学电子科学与工程学院,江苏南京210000;南京邮电大学电子科学与工程学院,江苏南京210000【正文语种】中文【中图分类】TP3911 引言啸叫现象在扩声系统中十分普遍,扩声设备发出的声音被拾音设备接收,构成一个闭环系统,声反馈信号在系统中循环放大,最终系统自激,产生尖锐的啸叫音。
啸叫不仅掩盖了有用信号,还可能对人耳听力造成损伤甚至烧坏扩声设备。
现行的声反馈抑制方法包括:(1)频率均衡法[4]——宽带陷波(缺点:需要专业调音师);(2)反馈抑制器法[5]——窄带陷波(缺点:误抑制人声);(3)调相法、移频法[6](缺点:影响音质);(4)随机相位法[7](缺点:效率不高);(5)调整距离法(缺点:弱适应性);(6)现有的自适应陷波器[8],缺点一是其针对工频50 Hz,二是它的自适应性为小范围[9]——在确定频率左右微调以求得最佳陷波效果。
语音信号同自然界大多数信号一样具有随机性和不稳定性,信号的时域分析、频域分析[1]在处理时变非平稳信号时仍有欠缺。
局部特征尺度分解结合局部均值解调的齿轮故障诊断
局部特征尺度分解结合局部均值解调的齿轮故障诊断崔伟成;许爱强;李伟;孟凡磊【期刊名称】《中国机械工程》【年(卷),期】2016(027)024【摘要】为了准确地进行故障诊断,根据齿轮故障振动信号的多分量调幅调频特征,提出了一种新的解调方法———局部均值解调法,将之与局部特征尺度分解相结合进行齿轮故障诊断。
该诊断方法首先对齿轮振动信号运用局部特征尺度分解,得到若干个内禀尺度分量,然后应用局部均值解调法求取每个分量的调频分量,最后根据瞬时频率的频谱进行故障诊断。
采用仿真信号将局部均值解调法与 Hilbert 解调法、经验调幅调频分解法进行了对比,结果表明,局部均值解调法的精确性更好。
通过齿轮故障振动数据的分析,验证了局部特征尺度分解结合局部均值解调的故障诊断方法的有效性。
【总页数】6页(P3332-3337)【作者】崔伟成;许爱强;李伟;孟凡磊【作者单位】海军航空工程学院,烟台,264001;海军航空工程学院,烟台,264001;海军航空工程学院,烟台,264001;海军航空工程学院,烟台,264001【正文语种】中文【中图分类】TN911.23;TP206.3【相关文献】1.基于局部特征尺度分解与瞬时频率谱的齿轮故障诊断 [J], 崔伟成;许爱强;李伟;孟凡磊2.一种局部特征尺度分解局部均值点求取方法 [J], 崔伟成;杨利斌;桑德一3.局部特征尺度分解方法及其在齿轮故障诊断中的应用 [J], 程军圣;杨怡;杨宇4.基于局部特征尺度分解的经验包络解调方法及其在机械故障诊断中的应用 [J], 程军圣;郑近德;杨宇5.基于局部特征尺度分解和数学形态滤波的齿轮故障诊断方法研究 [J], 齐嘉兴;崔伟成;赵修平因版权原因,仅展示原文概要,查看原文内容请购买。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
j= 1
通过上述步骤 , LCD方法可将一个多分量信号 分解成若干 ISC 分量之和。 与 IT D 方法相比 , L CD 方法在每次迭代过程中都对剩余信号进 行分量判 定 ,并通过多次迭代自适应地获得 ISC 分量。
图 1 I T D方法中基线信号的构造 Fig. 1 Co nstructio n o f ba seline sig na l fo r IT D metho d
1 2
量组成 ,任何两个 ISC 分量之间相互独立 , 这样一个 多分量信号 x ( t )就可以被分解为有限个 ISC 分量之 和 ,其中任何一个 ISC 分量满足以下条件 : ( 1)在整个数据段内 , 任意两个相邻极值点符号 互异 ; ( 2)在整个数据段内 , 其极值点为 X k ( k = 1, 2, … , M ) , 各个极值点相对应的时刻为 f k ( k = 1, 2,… , M ) ,由任意两个极大 (小 )值点 (f k , X k ) , (f k+ 2 , X k+ 2 ) 连 接 形 成的 线 段 在 其 中 间 极 小 ( 大 ) 值 点 (f k+ 1 , k+ 1的函数值 Ak+ 1与该极小 (大 )值 X k+ 1 )相对应时刻 f X k+ 1的比值关系近似不变 。 根据所定义的 ISC分量 , 对信号 x ( t )进行 L CD 方法分解 , 可将其分解为若干个 ISC 分量之和 , 其分 解过程如下: ( 1)确定原始信号 x ( t )的所有极值点 (f k , Xk ) , 利用式 ( 1)~ ( 3)构造基线信号 m 11 ( t ) = Lt ( 4) ( 2)将 m 11 ( t )从原始信号中分离出来 , 得到剩余 信号 h 11 ( t ) = x ( t ) - m 11 ( t ) ( 5) 理想地 ,如果 h 11 ( t )满足 ISC 分量判据 , 则 h11 ( t ) 为信号 x ( t )的第一个分量 ISC1 ( t ) ; 如果 h 11 ( t )不满足 ISC 分量判据 , 将 h 11 ( t )作为原始信号重复步骤 ( 1) , 则循环 i 次得到 h1i ( t ) = h1i - 1 ( t ) - m 1i ( t ) , 使得 h 1i ( t ) 满足 ISC 分量判据 , h 1i 即为信号 x ( t )的第一个分量 ISC1 ( t )。 ( 3)将 ISC1 ( t )从原始信号中分离出来 , 得到 r1 ( t ) = x ( t ) - ISC1 ( t ) ( 6) 将 r1 ( t )作为原始信号重复步骤 ( 1)和 ( 2) , 得到 x ( t )的满足 ISC分量判据的第二个分量 ISC2 ( t ) , 重 复循环 n 次 , 得到信号 x ( t )的 n 个满足 ISC分量判据 的分量 , 直到 rn ( t )为一单调函数或者小于预设阀值 为止 。 这样便可以将 x ( t )分解为 n 个 ISC 分量和一个 剩余函数 rn ( t )之和 , 即
n
- f k ( X k+ - f k
2
- Xk )
( 1)
这样可在点 (f k+ 1 , Ak+ 1 ) 与 (f k+ 1 , X k+ 1 )之间用 线性插值得到基线信号控制点 (f k+ 1 , Lk#43; 1 + ( 1 - a ) X k+
1
第 25 卷第 5 期 2012年 10 月
振 动 工 程 学 报 Journal of Vibratio n Engi neeri ng
V o l. 25 No. 5 Oct. 2012
局部特征尺度分解方法及其分解能力研究
杨 宇 , 曾 鸣 , 程军圣
(湖南大学汽车车身先进设计制 造国家重点实验室 , 湖南 长沙 410082) 摘要 : 在研究内禀时间尺度分 解 ( Int rinsic Time-Scale Decompositio n, I T D)方法 的基础上提出了一 种新的自适应 时频分 析方法 — — 局部特征 尺度分解 ( Lo ca l Char acteristicscale Decompositio n, L CD) 方法 ,该 方法可以 自适应地 将一个 复杂的多分量信号 分解为若干个瞬时 频率具有物 理意义的内 禀尺度分 量 ( Int rinsic Scale Component , ISC ) 之和 。对 LCD方法的基本原理进行了 阐述 ,通过建立分解能力研究模型分别研究了分量信 号的频率比 、幅值比和初 相位差以及单个分量迭代次数对 L CD方法分解能力的影响 ,并给出 了多分量信号可由 L CD方法分解的大致条件 。 关键词 : 自适应时频分析 ; 局部特征尺度分解 ; 内禀时间尺度分解 ; 分解能力 中图分类号 : T H165. 3; T N 911. 6 文献标 识码 : A 文章编号 : 1004-4523( 2012) 05-060208
2 仿真信号分析
2. 1 LCD与 ITD 的对比分析 考察下式所示仿真信号 , 采样率为 2 048 Hz,
1. 2 LCD 方法分解过程 L CD 方法假设任何复杂信号由不同的 ISC 分
604
振 动 工 程 学 报
第 25卷
t∈ [ 0, 1 ]。 x 1( t)= [1 + 0. 5si n( 5 πt ) ]cos [ 300 πt+ 2cos( 10 πt ) ] x 2 ( t ) = 2e- 2t si n( 60 πt ) x ( t )= x 1 ( t )+ x 2 ( t ) ( 8) 仿真信号 x ( t )由调幅调频信号 x 1 ( t )和调幅信 号 x 2 ( t )合成 , 其时域波形如图 2所示。
图 4 L CD方法分解结果 Fig. 4 Deco mposi tion r esult of L CD me tho d 图 2 仿真信号的时域波形 Fig. 2 Time do main w av efo r m of the simulated sig nal
2. 2 LCD与 EMD 的对比分析 本文仅从计算效率方面对 L CD 和 EMD 方法进 行对比分析 。 L CD 和 EM D方法一样 ,算法本身都存 在着端点效应 , 而端点效应会增大基线信号的偏差 , 继而增多迭代次数 , 降低计算效率 。 为了尽量减弱端 点效应对计算效率的影响 , 将仿真信号 x ( t )的仿真 时间延长至区间 [- 1, 2] , 这样就相当于对原始信号 进行了端点波形延拓 , 然后分别采用 LCD 和 EM D 方法对式 ( 8)所示仿真信号进行分解以对比两者的 计算效率 。 对 信号进行 LCD方法分解 , ISC 分量判据借鉴 EM D 方法中的标准差判据 , 通过多种类 型仿真信 号 (如正弦 、调幅、 调频以及调幅调频等信号之间的 相互叠加 )试验得到一个较为合理的标准差范围 : SD < 0. 01 。 LCD方法分解结果如图 5所示。 经过 2 次迭代后满足分量判据条件 , 得到第一个分量 ISC1 , 其运行时间为 0. 031 s; 经过 1 次迭代后得到第二个 分量 ISC2 , 其运行时间为 0. 016 s。 L CD 分解所需总 运行时间为 0. 047 s。
1
]
( 3) Lt k 表示对原始信号的第 k ( k = 1, 2, … , M - 1) 个区间进行线性变换后得到的基线信号段 , 将各个 极值点区间内的基线信号段首尾相连即可得到基线 信号 Lt , 如图 1中的点线所示 。
x (t) =
∑
ISCj ( t ) + rn ( t )
( 7)
之和 。 由于对于脑电波的癫痫信号采用单步迭代分
引 言
自适应时频分析方法的特点主要表现在不需要 对被分析信号的形态特征或者信息做出预测和限制 的前提下 , 可以在对信号进行分解的过程中根据信 号本身的特性自动产生基线信号 , 从而使得分解结 果具有一定的物理意义 。 近年来最具代表性的自适 应时频分析方法是经验模态分解 ( Empi rical M ode [ 1, 2] , 该方法在定义 Decom posi tio n, 简称 EM D)方法 瞬时频率具有物理意义的内禀模态函 数 ( Int ri nsic Mode Functio n,简称 IM F)分量的基础上 , 通过上下 极值点包络线的平均来构造基线信号 , 从而将复杂 的 多分量信号自适应地分解为若干个 IM F 分量之 和。 由此可见 , 自适应时频分析方法的关键在于基线 信号的构造。 文献 [ 3]介绍了一种新的基线信号构造方法 , 即 以原始信号任意两个相邻的极值点为跨度对信号进 行分段线性变换来构造基线信号 , 且由此提出了内 禀时间尺度分解方法 , 并将其应用于脑电波的癫痫 分析。 IT D方法每次只经过单步迭代就得到固有旋 转分量 ( Proper Rot atio n Com po nent, P RC) , 即把基 线信号从原始信号中分离后得到的剩余信号作为 P RC分量 , 并将基线信号作为下一次迭代的原始信 号 , 如此循环下去就可将信号分解为若干 P RC分量
第 5期
杨 宇 ,等 : 局部特征尺 度分解方法及其分解能力研究
603
1 LCD方法
1. 1 ITD 方法基线信号的构造 IT D 方法采用分段的形式对信号任意两个相 邻的极值点之间的数据段进行线性变换而获得基线 信号 [3 ]。 设原始信号 x ( t )的极值为 X k ( k = 1, 2, … , M ) , 其相对应的时刻为 f k ( k = 1, 2, … , M ) , 如图 1中 “· ” 所示 。 由任 意两 个 极大 (小 ) 值点 (f k , X k ) , (f k+ 2 , X k+ 2 ) 连接 形 成 的 线 段 在 其 中 间 极 小 (大 ) 值 点 (f k + 1 , X k+ 1 )相对应时刻 f k+ 1的值为 Ak + 1 = X k + f k+ f k+