High-Order Synchrosqueezed Chirplet Transforms for Multicomponent Signal Analysis

Yi-Ju Yen,De-Yan Lu,Sing-Yuan Yeh,Jian-Jiun Ding,Chun-Yen Shen

2024-05-11

Abstract:This study focuses on the analysis of signals containing multiple components with crossover instantaneous frequencies (IF). This problem was initially solved with the chirplet transform (CT). Also, it can be sharpened by adding the synchrosqueezing step, which is called the synchrosqueezed chirplet transform (SCT). However, we found that the SCT goes wrong with the high chirp modulation signal due to the wrong estimation of the IF. In this paper, we present the improvement of the post-transformation of the CT. The main goal of this paper is to amend the estimation introduced in the SCT and carry out the high-order synchrosqueezed chirplet transform. The proposed method reduces the wrong estimation when facing a stronger variety of chirp-modulated multi-component signals. The theoretical analysis of the new reassignment ingredient is provided. Numerical experiments on some synthetic signals are presented to verify the effectiveness of the proposed high-order SCT.

Numerical Analysis,Signal Processing

What problem does this paper attempt to address?

### What problem does this paper attempt to solve?

This paper aims to solve the problem of signal analysis with multiple components and cross - instantaneous frequencies (IF). Specifically, the existing synchrosqueezing modulation transform (SCT) performs poorly when dealing with strong frequency - modulated signals because of its inaccurate estimation of the instantaneous frequency (IF). This paper proposes a high - order synchrosqueezing modulation transform (HSCT) to improve the performance of existing methods in the face of complex multi - component signals.

#### Background and problem description

1. **Challenges in multi - component signal analysis**:
- In many practical applications (such as radar systems, seismic wave detection, and medical signals, etc.), the observed signals are usually oscillatory and composed of multiple components.
- The non - stationary characteristics of these signals make it difficult for traditional time - frequency analysis tools to effectively decompose and analyze these signals.

2. **Limitations of existing methods**:
- Methods such as the short - time Fourier transform (STFT) and empirical mode decomposition (EMD) need to meet specific separation conditions, that is, the instantaneous frequencies of different components must be spaced far enough apart.
- The synchrosqueezing transform (SST) and its high - order versions, although they improve the resolution of the time - frequency representation, are still limited by the separation conditions.
- The synchrosqueezing modulation transform (SCT) performs poorly when dealing with strong frequency - modulated signals because of its inaccurate estimation of the instantaneous frequency, especially in the presence of high - order modulation.

3. **The core of the problem**:
- When the signal contains high - order frequency - modulated modulation, existing methods cannot accurately estimate the instantaneous frequency and its derivatives, resulting in a decrease in the accuracy of signal decomposition and analysis.

#### Proposed solution

To overcome the above problems, this paper proposes a new high - order synchrosqueezing modulation transform (HSCT). The main improvements include:

1. **Introduction of high - order derivative estimation**:
- The new method not only estimates the first - and second - order derivatives of the phase function but also introduces the estimation of the third - order derivative, thus more accurately capturing the high - order modulation characteristics of the signal.
- The specific formulas are:
\[
\theta_0=\frac{1}{\pi i}\partial_t\left(\frac{\partial_t\left(\frac{\partial_t T(g)f}{T(g)f}\right)}{2 + \partial_t\left(\frac{q_1 f}{q_0 f}\right)}\right)
\]
\[
\mu_0+\theta_0 t=\frac{1}{2\pi i}\left(\frac{2q_0 f+\partial_t q_1 f}{\partial_t\left(\frac{\partial_t T(g)f}{T(g)f}\right)}-\frac{\partial_{tt}\left(\frac{\partial_t T(g)f}{T(g)f}\right)q_1 f}{2\left(q_0 f\right)^2}+\frac{\partial_t q_1 f q_0 f - q_1 f\partial_t q_0 f}{q_0 f}\right)
\]
\[
\omega_0+\mu_0 t+\frac{\theta_0}{2}t^2=\frac{1}{2\pi i}\left(\frac{\partial_t T(g)f}{T(g)f}-(\mu_0+\theta_0 t)\frac{T(tg)f}{T(g)f}-\frac{\theta_0}{2}\frac{T(t^2 g)f}{T(g)f}\right)
\]

2. **Improvement of signal decomposition accuracy**:
- The new method reduces...