
In this work, we studied the interleaving distance between merge trees from a combinatorial point of view. In the first part of the paper, we used a particular type of matching between trees to obtain a novel formulation of the distance. This formulation unveiled a link connecting the interleaving distance and edit distances between merge trees, which was of great interest due to the recursive decomposition properties and constrained formulations with polynomial time algorithms that these distances often enjoyed. In the second part of the paper, we built on this connection by applying tools from edit distances to obtain a constrained formulation of the interleaving distance and a recursive procedure which allowed us to find algorithms for upper and lower bounds of the interleaving distance. We implemented those algorithms and used them to test another upper bound presented by other authors, and tackled some simulations and case studies. Motivated by the literature on edit distances, we believe that our novel formulation could lead to novel heuristics to compute the interleaving distance and that applying our recursive scheme to the constrained interleaving distance would produce polynomial time upper bounds of practical relevance.
Citation: Matteo Pegoraro. A graph-matching formulation of the interleaving distance between merge trees[J]. AIMS Mathematics, 2025, 10(6): 13025-13081. doi: 10.3934/math.2025586
[1] | Zhoujin Cui, Xiaorong Zhang, Tao Lu . Resonance analysis and time-delay feedback controllability for a fractional horizontal nonlinear roller system. AIMS Mathematics, 2024, 9(9): 24832-24853. doi: 10.3934/math.20241209 |
[2] | P. Pirmohabbati, A. H. Refahi Sheikhani, H. Saberi Najafi, A. Abdolahzadeh Ziabari . Numerical solution of full fractional Duffing equations with Cubic-Quintic-Heptic nonlinearities. AIMS Mathematics, 2020, 5(2): 1621-1641. doi: 10.3934/math.2020110 |
[3] | Nasser. A. Saeed, Amal Ashour, Lei Hou, Jan Awrejcewicz, Faisal Z. Duraihem . Time-delayed control of a nonlinear self-excited structure driven by simultaneous primary and 1:1 internal resonance: analytical and numerical investigation. AIMS Mathematics, 2024, 9(10): 27627-27663. doi: 10.3934/math.20241342 |
[4] | Rashid Ali, Rabia Imtiaz, Moin-ud-Din Junjua, Fuad A. Awwad, Emad A. A. Ismail, Ahmed S. Hendy . Multiplicity of kink and hump structures in Van der Waals gas system. AIMS Mathematics, 2025, 10(5): 12493-12518. doi: 10.3934/math.2025564 |
[5] | Sunil Kumar, Jai Bhagwan, Lorentz Jäntschi . Numerical simulation of multiple roots of van der Waals and CSTR problems with a derivative-free technique. AIMS Mathematics, 2023, 8(6): 14288-14299. doi: 10.3934/math.2023731 |
[6] | Gamal M. Ismail, Maha M. El-Moshneb, Mohra Zayed . A modified global error minimization method for solving nonlinear Duffing-harmonic oscillators. AIMS Mathematics, 2023, 8(1): 484-500. doi: 10.3934/math.2023023 |
[7] | Hany Bauomy, A. T. EL-Sayed, A. M. Salem, F. T. El-Bahrawy . The improved giant magnetostrictive actuator oscillations via positive position feedback damper. AIMS Mathematics, 2023, 8(7): 16864-16886. doi: 10.3934/math.2023862 |
[8] | Changdev P. Jadhav, Tanisha B. Dale, Vaijanath L. Chinchane, Asha B. Nale, Sabri T. M. Thabet, Imed Kedim, Miguel Vivas-Cortez . On solutions of fractional differential equations for the mechanical oscillations by using the Laplace transform. AIMS Mathematics, 2024, 9(11): 32629-32645. doi: 10.3934/math.20241562 |
[9] | Guanghua Wei, Xiaorong Zhang, Zhoujin Cui . Primary resonance analysis of a novel fractional-order coronary artery model. AIMS Mathematics, 2025, 10(6): 14996-15011. doi: 10.3934/math.2025672 |
[10] | Guijun Xing, Huatao Chen, Zahra S. Aghayan, Jingfei Jiang, Juan L. G. Guirao . Tracking control for a class of fractional order uncertain systems with time-delay based on composite nonlinear feedback control. AIMS Mathematics, 2024, 9(5): 13058-13076. doi: 10.3934/math.2024637 |
In this work, we studied the interleaving distance between merge trees from a combinatorial point of view. In the first part of the paper, we used a particular type of matching between trees to obtain a novel formulation of the distance. This formulation unveiled a link connecting the interleaving distance and edit distances between merge trees, which was of great interest due to the recursive decomposition properties and constrained formulations with polynomial time algorithms that these distances often enjoyed. In the second part of the paper, we built on this connection by applying tools from edit distances to obtain a constrained formulation of the interleaving distance and a recursive procedure which allowed us to find algorithms for upper and lower bounds of the interleaving distance. We implemented those algorithms and used them to test another upper bound presented by other authors, and tackled some simulations and case studies. Motivated by the literature on edit distances, we believe that our novel formulation could lead to novel heuristics to compute the interleaving distance and that applying our recursive scheme to the constrained interleaving distance would produce polynomial time upper bounds of practical relevance.
Fractional calculus is an ancient mathematical topic from the 17th century. Although Fractional calculus is a mathematical subject with a history of more than 300 years, its application in physics and engineering did not attract wide attention until recent decades. In the past twenty years, significant progress has been made in fractional calculus and fractional order dynamical systems [1,2,3,4]. Previous studies have shown that fractional calculus is a mathematical tool, which is applicable to anomalous social and physical systems with non-local, frequency- and history-dependent properties, as well as intermediate states such as soft materials that are neither ideal solids nor ideal fluid [5,6,7]. Moreover, applications of fractional calculus have been reported in many areas such as signal processing, image processing, automatic control and robotics [8,9,10,11]. Research has even applied fractional derivatives to specific fields of research: psychological and life sciences [12,13]. In a word, virtually no area of classical analysis has been left untouched by fractional calculus. These examples and many other similar samples perfectly clarify the importance of consideration and analysis of dynamical systems with fractional-order models.
For the reason that nonlinear vibration systems exist in a large number of engineering technology and scientific practice, the research of nonlinear vibration systems has attracted the attention of many engineering and technical personnel and researchers. As a combination of two typical nonlinear systems, the nonlinear part of the Duffing-van der Pol system contains not only the nonlinear damping term for maintaining self-excited vibration of the van der Pol system, but also the cubic nonlinear restoring force term of the Duffing system, which has particularly rich dynamic characteristics. The Duffing-van der Pol oscillator is also a typical system to simulate some physical phenomena, and its simple nonlinear structure has led to an in-depth study of its dynamic behavior, which can be used as a research model in many fields, such as engineering, electronics, biology and neurology [14,15,16,17]. Furthermore, the fractional Duffing-van der Pol oscillator is considered to be one of the most commonly used equations in various system designs, including dynamics, biology, seismology, etc. [18,19]. Therefore, it is of practical significance to further study the dynamic behavior of such systems.
Many scientists have done a lot of research on the nonlinear dynamic behavior of the Duffing-van der Pol equation, including integer and fractional order cases [20,21,22,23,24,25]. Specially, in [22], a three-dimensional autonomous Van der Pol-Duffing type oscillator is proposed. Using the Routh-Hurwitz stability criterion and the linear stability analysis of equilibrium points, it was found that there is a Hopf bifurcation in the three-dimensional autonomous Van der Pol-Duffing oscillator. In [23], the damping characteristics of two Duffing-van der Pol oscillators with damping terms described by fractional derivative and time delay are studied by the residue harmonic balance method. A Duffing-van der Pol oscillator having fractional derivatives and time delays is investigated by the residue harmonic method in [24]. The dynamical properties of fractional order Duffing-van der Pol oscillators are studied in [25], and the amplitude-frequency response equation of the primary resonance is obtained by the harmonic balance method.
According to current research results, many studies focus on only Duffing-van der Pol equations with lower power terms under different excitations, while relatively few studies focus on the system with high power nonlinear terms. The Duffing-van der Pol systems with high power nonlinear terms are more worthy of study, because in the process of simplifying the actual model to the Duffing equation or van der Pol equation, in order to reduce the error between the system and the actual process, the nonlinear terms need to be approximately expressed as higher powers by power series to study the richer dynamic characteristics of the system. In addition, for the reason that the nonlinear dynamic behavior of the Duffing-van der Pol oscillator is very rich in bifurcation, chaos, primary resonance, superharmonic or subharmonic resonance, stochastic resonance and other vibration characteristics, the dynamic characteristics of the fractional Duffing-van der Pol oscillator are worthy of further exploration and research. It should be noted here that the primary resonance of a fractional Duffing-van der Pol oscillator can be studied by perturbation methods such as harmonic balance method or averaging method. Although a multiple scale method has been commonly used in the study of integer-order models, they are rarely involved in the study of fractional-order related problems. Therefore, it is also very necessary to use the multiple scale method to study the resonance characteristics of fractional Duffing-van der Pol oscillators, for example, to understand how the system parameters affect the periodic solutions resulted from resonance of the fractional Duffing-van der Pol oscillator.
The paper is organized as follows: In Section 2, the approximately analytical solution is obtained by the multiple scale method. The amplitude frequency equation is obtained and the stability conditions for the steady-state solution are given, while numerical simulations of the time history and phase diagrams are also provided. In Section 3, the bifurcations of primary resonance are investigated and the influence of parameters on fractional-order models is investigated. In Section 4, the comparisons of the approximate analytical solution and numerical solution are fulfilled. Finally, concluding remarks and implications in Section 5 close the paper.
The fractional Duffing-van der Pol oscillator with quintic nonlinear-restoring force is given by a second-order non-autonomous differential equation as follows
¨x−μ(1−x2)Dqtx+αx+λx3+βx5=fcosωt, | (2.1) |
where x is the position coordinate and is a function of time t, Dqtx is the q-order derivative of x with respect to time, μ is a positive fractional damping coefficient, ω0 is the natural frequency, α, λ and β are coefficients of nonlinear restoring forces, the external excitation have amplitude and frequency parameters which are f and ω, respectively. There are many definitions available for the fractional-order derivative, and in this study, Dqtx with 0<q≤1 is the Caputo's fractional derivative of x(t) described by
Dqtx(t)=1Γ(1−q)∫t0˙x(τ)(t−τ)−qdτ, |
in which Γ(z) is Gamma function satisfying Γ(z+1)=zΓ(z). Under certain conditions, primary resonance with limited amplitude occurs in the fractional Duffing-van der Pol oscillator (2.1) when ω=ω0def=√α or ω≈ω0. That is, the primary resonance means the excitation frequency is close to the natural one.
In order to find an approximate resonant solution of Eq (2.1), the multiple scale method will be used, for which a small parameter ε is required. Let Tn=εnt,(n=0,1) be two independent time variables, an approximate solution of Eq (2.1) with small amplitudes in three time scales can be represented by
x=x0(T0,T1)+εx1(T0,T1)+⋯. | (2.2) |
With such a small parameter ε, in order to balance the effect of the nonlinearity, damping and excitation, the following scaling is used
μ→εμ,λ→ελ,β→εβ,f→εf. | (2.3) |
Then, Eq (2.1) becomes
¨x−εμ(1−x2)Dqtx+ω20x+ελx3+εβx5=εfcosωt. | (2.4) |
The derivatives with respect to t can be expressed in terms of the new scaled times Tn using chain rule as follows:
ddt=D0+εD1+⋯, | (2.5a) |
d2dt2=D20+2εD0D1+ε2D21+⋯, | (2.5b) |
Dqt=Dq0+qεDq−10D1+⋯, | (2.5c) |
in which Dq0=∂q∂Tq0, Dn=∂∂Tn, D2n=∂2∂T2n and n=0,1.
For the primary resonance with ω=ω0:=√α or ω≈ω0, a detuning parameter σ describing the nearness of ω to ω0 is introduced by
ω=ω0+εσ, | (2.6) |
then ωt=ω0T0+σT1. Substituting (2.2), (2.5b) and (2.5c) into (2.4) leads to the following equation
(D20+2εD0D1+ε2D21)(x0+εx1)−εμ[1−(x0+εx1)2]⋅(Dq0+qεDq−10D1)(x0+εx1)+ω20(x0+εx1)+ελ(x0+εx1)3+εβ(x0+εx1)5=εfcos(ω0T0+σT1). |
Equating the coefficients of the same power of ε, a set of linear differential equations are obtained
O(ε0):D20x0+ω20x0=0, | (2.7) |
O(ε1):D20x1+ω20x1=−2D0D1x0+μ(1−x20)Dq0x0−λx30−βx50+fcos(ω0T0+σT1), | (2.8) |
from which x0 and x1 can solved one-by-one, respectively. In this way, the resonant solution x is dominated by x0 and collected by εx1.
In order to obtain the primary resonance solutions, it is necessary to eliminate the secular term when solving the Eqs (2.7) and (2.8) one by one, so as to find the amplitude-frequency response equation, from which the influence of various parameters on the resonance solution can be studied.
The general solution of Eq (2.7) is of the form,
x0=A(T1)eiω0T0+¯A(T1)e−iω0T0, | (2.9) |
where A(T1) and ¯A(T1) are unknown functions, ¯A(T1) denotes the complex conjugate of A(T1).
To solve Eq (2.8), for sufficient large t≫1, the qth-order (0≤q≤1) derivative of eiΩt is approximated with
DqteiΩt≈−∞DqteiΩt=(iΩ)qeiΩt, | (2.10) |
according to the "Short Memory Principle" [13,26], where
−∞Dqtx(t)=1Γ(1−q)∫t−∞˙x(τ)(t−τ)−qdτ. |
Substituting Eqs (2.9) and (2.10) into Eq (2.8) and using
cos(ω0T0+σT1)=ei(ω0T0+σT1)+e−i(ω0T0+σT1)2, |
the right-hand of Eq (2.8) becomes
[−2iω0D1A+μ(1−2AˉA)A(iω0)q−3λA2ˉA−10βA3ˉA2+f2eiσT1]eiω0T0+NST+cc, |
where NST stands for the terms that do not produce secular terms, cc denotes the complex conjugate of the preceding terms.
In order that x1 is periodic, the secular terms with eiω0T0 must be zero, namely
−2iω0D1A+μ(1−2AˉA)A(iω0)q−3λA2ˉA−10βA3ˉA2+f2eiσT1=0. | (2.11) |
To analyze the solution of Eq (2.11), it is convenient to express A(T1) in the polar form as
A(T1)=a(T1)2eiθ(T1). | (2.12) |
In addition, the Euler formula gives
iq=(eiπ/2)q=eiqπ/2=cosqπ2+isinqπ2. |
Let φdef=σT1−θ. By separating the real and imaginary parts of Eq (2.11), the differential equations governing amplitude a(T1) and φ(T1) of A(T1) are as follows respectively
D1a=μ(1−a22)a2ωq−10sinqπ2+f2ω0sinφ, | (2.13a) |
aD1φ=σa+μ(1−a22)a2ωq−10cosqπ2−3λ8ω0a3−5β16ω0a5+f2ω0cosφ. | (2.13b) |
We are interested in the steady solutions satisfying D1a=0, D1φ=0 in (2.13a) and (2.13b), namely
μ(2−a2)4aωq−10sinqπ2=−f2ω0sinφ, | (2.14a) |
σa+μ(2−a2)4aωq−10cosqπ2−3λ8ω0a3−5β16ω0a5=−f2ω0cosφ. | (2.14b) |
Because sin2φ+cos2φ=1, by eliminating φ from Eqs (2.14a) and (2.14b), the amplitude-frequency response equation can be obtained as follows
[(4μ(2−a2)ωq0sinqπ2)2+(16σω0+4μ(2−a2)ωq0cosqπ2−6λa2−5βa4)2]a2=(8f)2. | (2.15) |
The real solution a of Eq (2.15) determines the dominant part x0 of the primary resonance response amplitude. Similar results have been discussed in [21] for the case when q=1, but the calculation error in the derivation process of this paper leads to inaccurate amplitude-frequency response equation.
The fractional Duffing-van der Pol oscillator with negative acceleration feedback considered is as follows
¨x−μ(1−x2)Dqtx+ω20x+λx3+βx5=fcosωt−u¨x, | (2.16) |
in which −u is the feedback gain reflecting the magnitude of a feedback force. Similar to the previous discussion, the method of multiple time scales is used to obtain a uniformly valid, asymptotic expansion of the solution for Eq (2.16). Let Tn=εnt,(n=0,1) be two independent time variables, an approximate solution of Eq (2.16) with small amplitudes in three time scales can be represented by
x=x0(T0,T1)+εx1(T0,T1)+⋯. | (2.17) |
With such a small parameter ε, in order to balance the effect of the nonlinearity, damping and excitation, the following scaling is used
μ→εμ,λ→ελ,β→εβ,f→εf,u→εu. | (2.18) |
Then, Eq (2.16) becomes
¨x−εμ(1−x2)Dqtx+ω20x+ελx3+εβx5=εfcosωt−εu¨x. | (2.19) |
Substituting (2.17), (2.5b) and (2.5c) into (2.19) leads to the following equation
(D20+2εD0D1+ε2D21)(x0+εx1)−εμ[1−(x0+εx1)2]⋅(Dq0+qεDq−10D1)(x0+εx1)+ω20(x0+εx1)+ελ(x0+εx1)3+εβ(x0+εx1)5=εfcos(ω0T0+σT1)−εu(D20+2εD0D1+ε2D21)(x0+εx1). |
Similar to the previous discussion, the following differential equations about amplitude and phase can be obtained,
D1a=μ(1−a22)a2ωq−10sinqπ2+f2ω0sinφ, | (2.20a) |
aD1φ=σa+μ(1−a22)a2ωq−10cosqπ2−3λ8ω0a3−5β16ω0a5+f2ω0cosφ+uω0a2. | (2.20b) |
The steady solutions satisfying D1a=0, D1φ=0 in (2.20a) and (2.20b), namely
μ(2−a2)4aωq−10sinqπ2=−f2ω0sinφ, | (2.21a) |
σa+μ(2−a2)4aωq−10cosqπ2−3λ8ω0a3−5β16ω0a5+uω0a2=−f2ω0cosφ. | (2.21b) |
By eliminating φ from Eqs (2.21a) and (2.21b), the amplitude-frequency response equation can be obtained as follows
[(4μ(2−a2)ωq0sinqπ2)2+(16σω0+4μ(2−a2)ωq0cosqπ2−6λa2−5βa4+8uω20)2]a2=(8f)2. | (2.22) |
It should be noted that a fractional dynamical system does not have an exact periodic solution, see [27,28], and it may have one if the fractional derivative is defined over (−∞,t] [28]. When the asymptotic behaviors of the fractional system are addressed, a fractional derivative over [0,t] with large t equals approximately the corresponding one over (−∞,t], according to the "Short Memory Principle" of fractional derivatives [26]. Thus, a periodic solution of a fractional system should be understood as an asymptotic solution with transients before a finite time neglected. All the periodic solutions mentioned in this paper are understood in this sense.
To determine the stability of the steady resonant solutions, suppose that (a,φ)=(a∗,φ∗) is a steady solution of Eqs (2.20a) and (2.20b), let
Δa=a−a∗,Δφ=φ−φ∗ | (2.23) |
are perturbations which are assumed to be small compared to a∗ and φ∗. Substituting Eq (2.23) into Eqs (2.20a) and (2.20b), using Eqs (2.21a) and (2.21b), then the linearized differential equations governing Δa and Δφ are
D1Δa=2−3(a∗)24μωq−10sinqπ2Δa+f2ω0cosφ∗Δφ, | (2.24a) |
D1Δφ=[σa∗+2−3(a∗)24a∗μωq−10cosqπ2−9λa∗8ω0−25β(a∗)38ω0+uω02a∗]Δa−f2ω0a∗sinφ∗Δφ. | (2.24b) |
Let M=σ+2−(a∗)24μωq−10cosqπ2−3λ8ω0(a∗)2−5β16ω0(a∗)4+uω02, N=σ+2−3(a∗)24μωq−10cosqπ2−9λ8ω0(a∗)2−25β8ω0(a∗)4+uω02. According to Eqs (2.21a) and (2.21b), the characteristic equation can be written as
|2−3(a∗)24μωq−10sinqπ2−λ−a∗M1a∗N2−(a∗)24μωq−10sinqπ2−λ|=0. |
By expanding the determinant, we can get
λ2+[(a∗)2−1]μωq−10sinqπ2λ+Λ=0, | (2.25) |
in which Λ=[2−3(a∗)2][2−(a∗)2]16(μωq−10sinqπ2)2+MN.
Thus, the steady solution (a,φ)=(a∗,φ∗) is asymptotically stable if and only if the real parts of all roots of the characteristic Eq (2.25) are negative.
In this section, time history diagrams are used to briefly illustrate the impact of feedback control on the model. Figure 1 shows the time history of the system without controller. The various parameters of the system in Figure 1 are ω0=1, q=0.9, μ=0.01, λ=0.04, β=0.01 and f=0.06. It can be observed that the amplitude fluctuates significantly in the first 150s and the steady-state amplitude of the system is about 0.28. The phase diagram shows a limit cycle, indicating that there is no chaos in the system. Different initial conditions were tested and it was found that the steady-state amplitude of the system is not sensitive to the initial conditions.
Figure 2 shows the time history of the system with feedback controller. The various parameters of the system in Figure 2 are ω0=1, q=0.9, μ=0.01, λ=0.04, β=0.01, f=0.06 and u=0.1. From Figure 2, it can be seen that the steady-state amplitude of the system has decreased to approximately 0.115, this indicates the effectiveness of the controller. The steady-state value appears after 150s without control, while it appears before 50 seconds with control, which also reflects the effect of adding control.
Rewritten Eqs (2.15) and (2.22) in terms of the original system parameters as follows
[(4μ(2−a2)ωq0sinqπ2)2+(16(ω−ω0)ω0+4μ(2−a2)ωq0cosqπ2−6λa2−5βa4)2]a2=(8f)2, | (3.1) |
and
[(4μ(2−a2)ωq0sinqπ2)2+(16(ω−ω0)ω0+4μ(2−a2)ωq0cosqπ2−6λa2−5βa4+8uω20)2]a2=(8f)2. | (3.2) |
The effect of different parameters on the amplitude of the resonant solutions is investigated numerically. Here, the difference between the model with quintic nonlinear term and the model without quintic power is given by the amplitude-frequency response curve, as shown in Figure 3.
It can be seen from the figure that the bending degree, resonance region and jump phenomena of the model have obvious changes after considering the quintic nonlinear terms, indicating that it is very important to consider the role of higher order terms for studying the nonlinear characteristics of the model.
This section discusses the influence of parameters on the hardening/softening characteristics and resonance peak value based on Eq (3.1). As shown in Figures 4–7, typical hardening/softening characteristics can be observed also in the curves of the resonant amplitudes with respect to the changes of λ, β, q, μ and f, respectively.
The influence of the nonlinear terms on the amplitude-frequency response curve can be identified with the changes of λ, β. Typical softening characteristic can be observed in Figure 4 with λ=−0.04 and β=−0.01, and the frequency response curve bends to low frequency. When λ=0.04 and β=0.01, hardening characteristic can be observed and the frequency response curve bends to high frequency.
The effect of the fractional-order q on the amplitude-frequency curves is shown in Figure 5. If the fractional order q is closer to 0, the fractional differential term degenerates to linear stiffness, while if q is closer to 1, it degenerates to linear damping. It is further known that the smaller the order q is, the larger the maximum amplitude is. The integer-order model is the lowest. Compared with the integer order case when q=1, the bending degree, resonance peak and resonance region of the amplitude-frequency curve of the fractional-order model change accordingly with the decrease of the fractional-order q.
Figure 6 shows the effect of the damping coefficient μ on the amplitude-frequency curves. With the increase of μ, the nonlinear jump of the system weakens and the resonance amplitude of the system decreases. In other words, as μ increases, the unstable portions decrease.
The effect of the excitation amplitude f on the amplitude-frequency curves can be seen in Figure 7. With the increase of the excitation amplitude f, the frequency response curve will generate a multivalued region, the nonlinear jump of the system is more obvious and the resonance range and resonance amplitude of the system increase.
Figures 5–7 present parameter effects on resonant amplitudes with respect to excitation frequency. With fixed parameter values, all the figures exhibit typical characteristics of hardening spring. Let us take a different view on the resonant solutions. Now, the excitation frequency is fixed, and then with the change of a parameter, the amplitude of the primary resonant solution changes in different ways. In Figures 8–10, the images reflect characteristics similar to hardening spring, where the peak of amplitude is bent to the right as f varies, as λ varies or as β varies, respectively.
Through Eq (3.2), this section discusses the influence of parameters on the resonant solutions under feedback control. Figure 11 considers the influence of the gain value of the feedback control on the vibration state of the system. The red line in Figure 11 shows the primary resonance amplitude frequency response curve of the system without feedback control. It can be seen from Figure 11 that compared with the amplitude frequency curve of the primary resonance of the system without feedback control, after adding the feedback control, the resonance curve of the system moves to the left, and the resonance domain of the system changes. When the amplitude of the value of the gain of the control is increased, the deviation of the amplitude frequency response curve becomes larger.
Figure 12 shows the influence of the frequency parameters on the vibration state of the system. The closer the external excitation frequency ω is to the primary resonance frequency ω0, the smaller the offset of amplitude frequency curve is, and the smaller the change of resonance domain is. On the contrary, the more the phase difference, the larger the offset of the amplitude frequency curve. The peak amplitude of the primary resonance in Figure 12 does not change, but changes the position where it appears.
According to Eq (3.1), the primary resonance amplitude-frequency response curve of the system can be drawn. For comparison, this paper adopts the power series method introduced in references [3,29], and its calculation formula is
Dqtn[y(tn)]≈h−qn∑j=0Cqjy(tn−j), | (4.1) |
where tn=nh is the sample points, h is the sample step and Cqj is the fractional binomial coefficient with the iterative relationship as
Cq0=1,Cqj=(1−1+qj)Cqj−1. | (4.2) |
According to Eqs (4.1) and (4.2), the numerical scheme for Eq (2.4) can be expressed as
x(tn)=y(tn−1)h−n∑j=1C1jx(tn−j), | (4.3a) |
y(tn)={fcos(ωtn)−αx(tn)−λx3(tn)−βx5(tn)+μ[1−x2(tn)]z(tn−1)}h−n∑j=1C1jy(tn−j), | (4.3b) |
z(tn)=y(tn)h1−q−n∑j=1C1−qjz(tn−j). | (4.3c) |
The numerical amplitude-frequency curve marked with circle in Figure 13, where the stepsize of time is h=0.005, and the total computation time is 100s with the first 25s neglected. It shows that the resonant amplitude calculated from Eq (3.1) is in good agreement with the numerical results, especially when the ω≈ω0.
We study the primary resonance and feedback control of the fractional Duffing-van der Pol oscillator with quintic nonlinear-restoring force. The approximate analytical solution and the amplitude-frequency equation are obtained by the multiple scale method. The stability condition for steady-state solution is obtained based on the Lyapunov theory. The influence of parameters on system characteristics was studied using the amplitude frequency response equation and phase plane technology, and bifurcation analysis was conducted to verify the stability of the system. It can also be concluded that fractional orders and damping coefficients are very important in the fractional Duffing-van der Pol oscillator. For example, larger fractional orders and larger damping coefficients can reduce the effective amplitude of resonance and change the resonance frequency.
The authors declare they have not used Artificial Intelligence (AI) tools in the creation of this article.
The author expresses gratitude to the anonymous reviewers for their helpful comments and suggestions, as well as to the financial support from High Level Talent Research Launch Fund of Jiangsu Second Normal University (No. 928201/058).
The author declares that the research was conducted in the absence of any conflict of interest.
[1] | H. Adams, J. Bush, N. Clause, F. Frick, M. Gómez, M. Harrison, et al., Gromov-hausdorff distances, borsuk-ulam theorems, and vietoris-rips complexes, (2022), arXiv: 2301.00246. |
[2] | H. Adams, T. Emerson, M. Kirby, R. Neville, C. Peterson, P. Shipman, et al., Persistence images: A stable vector representation of persistent homology, J. Mach. Learn. Res., 18 (2017), 1–35. |
[3] |
P. K. Agarwal, K. Fox, A. Nath, A. Sidiropoulos, Y. S. Wang, Computing the gromov-hausdorff distance for metric trees, ACM T. Algorithms, 14 (2018), 20. https://doi.org/10.1145/3185466 doi: 10.1145/3185466
![]() |
[4] |
T. Akutsu, D. Fukagawa, A. Takasu, T. Tamura, Exact algorithms for computing the tree edit distance between unordered trees, Theor. Comput. Sci., 412 (2011), 352–364. https://doi.org/10.1016/j.tcs.2010.10.002 doi: 10.1016/j.tcs.2010.10.002
![]() |
[5] |
O. Al-Jowder, E. K. Kemsley, R. H. Wilson, Detection of adulteration in cooked meat products by mid-infrared spectroscopy, J. Agric. Food Chem., 50 (2002), 1325–1329. https://doi.org/10.1021/jf0108967 doi: 10.1021/jf0108967
![]() |
[6] |
M. Balasubramanian, E. L. Schwartz, The isomap algorithm and topological stability, Science, 295 (2002), 7. https://doi.org/10.1126/science.295.5552.7a doi: 10.1126/science.295.5552.7a
![]() |
[7] |
U. Bauer, C. Landi, F. Mémoli, The reeb graph edit distance is universal, Found. Comput. Math., 21 (2021), 1441–1464. https://doi.org/10.1007/s10208-020-09488-3 doi: 10.1007/s10208-020-09488-3
![]() |
[8] | U. Bauer, E. Munch, Y. S. Wang, Strong Equivalence of the interleaving and functional distortion metrics for reeb graphs, 31st International Symposium on Computational Geometry (SoCG 2015), Schloss Dagstuhl, Leibniz-Zentrum für Informatik, 2015,461–475. https://doi.org/10.4230/LIPIcs.SOCG.2015.461 |
[9] | K. Beketayev, D. Yeliussizov, D. Morozov, G. H. Weber, B. Hamann, Measuring the distance between merge trees, In: Topological methods in data analysis and visualization III, Cham: Springer, 2014,151–165. https://doi.org/10.1007/978-3-319-04099-8_10 |
[10] | N. Berkouk, Algebraic homotopy interleaving distance, In: Geometric science of information, Cham: Springer, 2021,656–664. https://doi.org/10.1007/978-3-030-80209-7_70 |
[11] |
H. B. Bjerkevik, M. B. Botnan, M. Kerber, Computing the interleaving distance is NP-hard, Found. Comput. Math., 20 (2020), 1237–1271. https://doi.org/10.1007/s10208-019-09442-y doi: 10.1007/s10208-019-09442-y
![]() |
[12] | P. Bubenik, Statistical topological data analysis using persistence landscapes, J. Mach. Learn. Res., 16 (2015), 77–102. |
[13] | R. Cardona, J. Curry, T. Lam, M. Lesnick, The universal ℓp-metric on merge trees, 38th International Symposium on Computational Geometry (SoCG 2022), Schloss Dagstuhl, Leibniz-Zentrum für Informatik, 2022, 1–24. https://doi.org/10.4230/LIPIcs.SoCG.2022.24. |
[14] | F. Chazal, D. Cohen-Steiner, M. Glisse, L. J. Guibas, S. Y. Oudot, Proximity of persistence modules and their diagrams, In: Proceedings of the twenty-fifth annual symposium on computational geometry, New York: Association for Computing Machinery, 2009,237–246. https://doi.org/10.1145/1542362.1542407 |
[15] |
F. Chazal, B. T. Fasy, F. Lecci, A. Rinaldo, L. Wasserman, Stochastic convergence of persistence landscapes and silhouettes, J. Comput. Geom., 6 (2015), 140–161. https://doi.org/10.20382/jocg.v6i2a8 doi: 10.20382/jocg.v6i2a8
![]() |
[16] |
J. Curry, H. B. Hang, W. Mio, T. Needham, O. B. Okutan, Decorated merge trees for persistent topology, J. Appl. Comput. Topology, 6 (2022), 371–428. https://doi.org/10.1007/s41468-022-00089-3 doi: 10.1007/s41468-022-00089-3
![]() |
[17] | J. Curry, W. Mio, T. Needham, O. B. Okutan, F. Russold, Convergence of leray cosheaves for decorated mapper graphs, (2023), arXiv: 2303.00130. |
[18] |
V. de Silva, E. Munch, A. Patel, Categorified reeb graphs, Discrete Comput. Geom., 55 (2016), 854–906. https://doi.org/10.1007/s00454-016-9763-9 doi: 10.1007/s00454-016-9763-9
![]() |
[19] | V. de Silva, E. Munch, A. Stefanou, Theory of interleavings on categories with a flow, Theor. Appl. Categ., 33 (2018), 583–607. |
[20] |
B. Di Fabio, C. Landi, The edit distance for reeb graphs of surfaces, Discrete Comput. Geom., 55 (2016), 423–461. https://doi.org/10.1007/s00454-016-9758-6 doi: 10.1007/s00454-016-9758-6
![]() |
[21] |
H. Edelsbrunner, D. Letscher, A. Zomorodian, Topological persistence and simplification, Discrete Comput. Geom., 28 (2002), 511–533. https://doi.org/10.1007/s00454-002-2885-2 doi: 10.1007/s00454-002-2885-2
![]() |
[22] | H. Edelsbrunner, J. Harer, Persistent homology–a survey, In: Surveys on discrete and computational geometry: twenty years later, Providence: American Mathematical Society, 2008,257–282. https://doi.org/10.1090/conm/453 |
[23] | Y. Elkin, V. Kurlin, The mergegram of a dendrogram and its stability, 45th International Symposium on Mathematical Foundations of Computer Science (MFCS 2020), Schloss Dagstuhl, Leibniz-Zentrum für Informatik, 2020, 1–13. https://doi.org/10.4230/LIPIcs.MFCS.2020.32 |
[24] | K. H. Esbensen, D. Guyot, F. Westad, L. P. Houmoller, Multivariate data analysis: in practice: An introduction to multivariate data analysis and experimental design, 5 Eds., AS: Camo Process, 2002. |
[25] |
E. F. Touli, Y. S. Wang, FPT-Algorithms for computing gromov-hausdorff and interleaving distances between trees, J. Comput. Geom., 13 (2022), 89–124. https://doi.org/10.20382/jocg.v13i1a4 doi: 10.20382/jocg.v13i1a4
![]() |
[26] |
M. Febrero, P. Galeano, W. González-Manteiga, Outlier detection in functional data by depth measures, with application to identify abnormal Nox levels, Environmetrics, 19 (2008), 331–345. https://doi.org/10.1002/env.878 doi: 10.1002/env.878
![]() |
[27] | F. Ferraty, P. Vieu, Nonparametric functional data analysis: theory and practice, New York: Springer, 2006. https://doi.org/10.1007/0-387-36620-2 |
[28] |
E. Gasparovic, E. Munch, S. Oudot, K. Turner, B. Wang, Y. S. Wang, Intrinsic interleaving distance for merge trees, La Matematica, 4 (2025), 40–65. https://doi.org/10.1007/s44007-024-00143-9 doi: 10.1007/s44007-024-00143-9
![]() |
[29] | A. Hatcher, Algebraic topology, Cambrige: Cambridge University Press, 2001. |
[30] | K. Hirata, Y. Yamamoto, T. Kuboyama, Improved max snp-hard results for finding an edit distance between unordered trees, In: Combinatorial pattern matching, Berlin: Springer, 2011,402–415. https://doi.org/10.1007/978-3-642-21458-5_34 |
[31] | E. Hong, Y. Kobayashi, A. Yamamoto, Improved methods for computing distances between unordered trees using integer programming, In: Combinatorial optimization and applications, Cham: Springer, 2017, 45–60. https://doi.org/10.1007/978-3-319-71147-8_4 |
[32] |
T. Jiang, L. S. Wang, K. Z. Zhang, Alignment of trees–an alternative to tree edit, Theor. Comput. Sci., 143 (1995), 137–148. https://doi.org/10.1016/0304-3975(95)80029-9 doi: 10.1016/0304-3975(95)80029-9
![]() |
[33] |
L. Kanari, A. Garin, K. Hess, From trees to barcodes and back again: theoretical and statistical perspectives, Algorithms, 13 (2020), 335. https://doi.org/10.3390/a13120335 doi: 10.3390/a13120335
![]() |
[34] |
M. Lesnick, The theory of the interleaving distance on multidimensional persistence modules, Found. Comput. Math., 15 (2015), 613–650. https://doi.org/10.1007/s10208-015-9255-y doi: 10.1007/s10208-015-9255-y
![]() |
[35] | S. M. Lane, Categories for the working mathematician, 2Eds., New York: Springer, 1998. |
[36] |
J. S. Marron, J. O. Ramsay, L. M. Sangalli, A. Srivastava, Functional data analysis of amplitude and phase variation, Statist. Sci., 30 (2015), 468–484. https://doi.org/10.1214/15-STS524 doi: 10.1214/15-STS524
![]() |
[37] | J. W. Milnor, M. Spivak, R. Wells, Morse theory, Princeton: Princeton University Press, 1963. https://doi.org/10.1515/9781400881802 |
[38] | E. Munch, B. Wang, Convergence between categorical representations of Reeb space and mapper, In: 32nd International Symposium on Computational Geometry (SoCG 2016), Schloss Dagstuhl, Leibniz-Zentrum für Informatik, 2016, 1–16. https://doi.org/10.4230/LIPIcs.SoCG.2016.53 |
[39] | M. Pegoraro, A persistence-driven edit distance for graphs with abstract weights, (2024), arXiv: 2304.12088v3. |
[40] | M. Pegoraro, A finitely stable edit distance for functions defined on merge trees, (2024), arXiv: 2108.13108. |
[41] | M. Pegoraro, A finitely stable edit distance for merge trees, (2024), arXiv: 2111.02738. |
[42] | M. Pegoraro, P. Secchi, Functional data representation with merge trees, (2024), arXiv: 2108.13147. |
[43] |
M. Pont, J. Vidal, J. Delon, J. Tierny, Wasserstein distances, geodesics and barycenters of merge trees, IEEE T. Vis. Comput. Gr., 28 (2021), 291–301. https://doi.org/10.1109/TVCG.2021.3114839 doi: 10.1109/TVCG.2021.3114839
![]() |
[44] |
C. R.-Carreño, J. L. Torrecilla, M. C. Berrocal, P. Manchón, A. Suárez, scikit-fda: A Python package for functional data analysis, J. Stat. Softw., 109 (2024), 1–37. https://doi.org/10.18637/jss.v109.i02 doi: 10.18637/jss.v109.i02
![]() |
[45] | J. O. Ramsay, B. W. Silverman, Functional data analysis, In: International encyclopedia of the social and behavioral sciences, 2 Eds., New York: Springer, 2005,514–518. https://doi.org/10.1016/B978-0-08-097086-8.42046-5 |
[46] |
S. M. Selkow, The tree-to-tree editing problem, Inform. Process. Lett., 6 (1977), 184–186. https://doi.org/10.1016/0020-0190(77)90064-3 doi: 10.1016/0020-0190(77)90064-3
![]() |
[47] |
P. Smith, V. Kurlin, Generic families of finite metric spaces with identical or trivial 1-dimensional persistence, J. Appl. Comput. Topology, 8 (2024), 839–855. https://doi.org/10.1007/s41468-024-00177-6 doi: 10.1007/s41468-024-00177-6
![]() |
[48] |
R. Sridharamurthy, T. B. Masood, A. Kamakshidasan, V. Natarajan, Edit distance between merge trees, IEEE T. Vis. Comput. Gr., 26 (2020), 1518–1531. https://doi.org/10.1109/TVCG.2018.2873612 doi: 10.1109/TVCG.2018.2873612
![]() |
[49] |
R. Sridharamurthy, V. Natarajan, Comparative analysis of merge trees using local tree edit distance, IEEE T. Vis. Comput. Gr., 29 (2021), 1518–1530. https://doi.org/10.1109/TVCG.2021.3122176 doi: 10.1109/TVCG.2021.3122176
![]() |
[50] | R. D. Tuddenham, M. M. Snyder, Physical growth of California boys and girls from birth to eighteen years, Publ. Child. Dev. Univ. Calif., 1 (1954), 183–364. |
[51] | L. M. Sangalli, P. Secchi, S. Vantini, V. Vitelli, Functional clustering and alignment methods with applications, Communications in Applied and Industrial Mathematics, 1 (2010), 205–224. |
[52] | F. Wetzels, C. Garth, A deformation-based edit distance for merge trees, 2022 Topological Data Analysis and Visualization (TopoInVis), Oklahoma City, USA, 2022, 29–38. https://doi.org/10.1109/TopoInVis57755.2022.00010 |
[53] |
F. Wetzels, H. Leitte, C. Garth, Branch decomposition-independent edit distances for merge trees, Comput. Graph. Forum, 41 (2022), 367–378. https://doi.org/10.1111/cgf.14547 doi: 10.1111/cgf.14547
![]() |
[54] |
K. Z. Zhang, A constrained edit distance between unordered labeled trees, Algorithmica, 15 (1996), 205–222. https://doi.org/10.1007/BF01975866 doi: 10.1007/BF01975866
![]() |
1. | Guohui Li, Ruiting Xie, Hong Yang, Study on fractional-order coupling of high-order Duffing oscillator and its application, 2024, 186, 09600779, 115255, 10.1016/j.chaos.2024.115255 | |
2. | M.N. Abd El-Salam, Rageh K. Hussein, Reducing the Primary Resonance Vibrations of a Cantilever Beam Using a Proportional Fractional-Order Derivative Controller, 2025, 13, 2227-7390, 1886, 10.3390/math13111886 | |
3. | Marina V. Shitikova, Elena I. Osipova, Combinational and Primary Simultaneous Resonance at Nonlinear Vibrations of a Fractionally Damped Plate, 2025, 20, 2224-3429, 79, 10.37394/232011.2025.20.10 |