本征时间尺度分解(Intrinsic Time-Scale Decomposition)

1.在诊断领域应用。

本征时间尺度分解能够很好地处理振动信号,分解出振动信号最具代表性的特征信号,在如滚动轴承故障诊断上,采取本征时间尺度分解作为信号处理手段来进行故障诊断能够大大提高诊断性能,

2.数学原理。

本征时间尺度分解最早由 Frei Osorio 在本征时间尺度分解——非平稳信号的时频能
量分析和实时滤波 一文提出,旨在将复杂信号分解为成多个相互独立的固有旋转分量和一
个单调函数或常数之和。具体计算下下:
给定一个待分解的复杂信号 F ( t ), 定义一个操作过程 P 以提取 F ( t ) 中的基线成分( P 相当于是提取低通信号的算子),当去除 F ( t ) 中的基线信号后,剩余的信号就成为了固有旋转分 量。一次分解过程可以描述为式(1)
{\color{Magenta} F(t)=PF(t)+(1-P)F(t)=P(t)+Q(t)}(1)
式中:
P ( t ) ——基线信号或者低通信号;
Q ( t ) ——固有旋转分量。
具体实现步骤如下:
Step 1:确定分析的实值信号 F ( t ),为了便于定义,τ0=0 是实值信号 F ( t ) 的第一个样本, 设在对应时刻τ k k =1,2,..., n )下的极值 F ( τ k )
Step 2:假设 P ( t ) H ( t ) 定义在 [0, τ k ] 之间, F ( t ) 的时间 t [0, τ k+2 ] 取值,基线提取算子 P 定义为两个极点之间的区间( τ k , τ k+1 ] 上的分段线性函数,如式(2)所示。
{\color{Magenta} PF(t)=P(t)=P(\tau_{k} )+\frac{P(\tau_{k+1})-P(\tau_{k} )}{F(\tau_{k+1})-F(\tau_{k} )}[F(t)-F(\tau_{k})],t\in (\tau_{k},\tau_{k+1}]}(2)
式中:
{\color{Magenta} P(\tau_{k+1})=\alpha\left \{ F(\tau_{k})+\frac{\tau_{k+1}-\tau_{k}}{\tau_{k+2}-\tau_{k}} \left [ F(\tau_{k+2})- F(\tau_{k}) \right ]\right \}+(1-\alpha ) F(\tau_{k+1}),k=1,2,...,n-2}
线性增益 \alpha取值为 0< \alpha <1在计算过程中常取 0.5
Step 3:将基线信号(低通信号) P ( t ) 重新当作新原信号,重复上述步骤,直至基线信
号(低通信号)为一个单调函数或常函数。因此复杂信号 F ( t ) 可分解为式(3)
{\color{Magenta} P(t)=PF(t)+QF(t)=(P+Q)F(t)+QF(t)}
{\color{Magenta} \Rightarrow \left [ Q(1+P)+P^{2} \right ]F(t)}
{\color{Magenta} \Rightarrow (Q\sum_{k=0}^{r-1}P^{k}+P^{r})F(t)}
{\color{Magenta} \Rightarrow Q_{t}^{1}+Q_{t}^{2}+...+Q_{t}^{r}+P_{t}^{r}}(3)
式中:
Q ——固有旋转分量提取算子;
Q_{t}^{i}——第 i 次分离出来的固有旋转分量;
P_{t}^{r}——为单调函数或常函数。
matlab编程见作者资源部分。
  • 1
    点赞
  • 5
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 1
    评论
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

L_2018

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值