Next Article in Journal
Numerical Integration Schemes Based on Composition of Adjoint Multistep Methods
Next Article in Special Issue
On the Numerical Treatment of the Temporal Discontinuity Arising from a Time-Varying Point Mass Attachment on a Waveguide
Previous Article in Journal
A Mask-Based Adversarial Defense Scheme
Previous Article in Special Issue
Modular Stability Analysis of a Nonlinear Stochastic Fractional Volterra IDE
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

On a Hypothetical Model with Second Kind Chebyshev’s Polynomial–Correction: Type of Limit Cycles, Simulations, and Possible Applications

by
Nikolay Kyurkchiev
1,2,*,† and
Anton Iliev
1,2,†
1
Faculty of Mathematics and Informatics, University of Plovdiv Paisii Hilendarski, 24, Tzar Asen Str., 4000 Plovdiv, Bulgaria
2
Institute of Mathematics and Informatics, Bulgarian Academy of Sciences, Acad. G. Bonchev Str., Bl. 8, 1113 Sofia, Bulgaria
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Algorithms 2022, 15(12), 462; https://doi.org/10.3390/a15120462
Submission received: 2 November 2022 / Revised: 30 November 2022 / Accepted: 30 November 2022 / Published: 6 December 2022
(This article belongs to the Special Issue Computational Methods and Optimization for Numerical Analysis)

Abstract

:
In this article, we explore a new extended Lienard-type planar system with “corrections” of the second kind Chebyshev’s polynomial U n . The number and type of limit cycles are also studied. The discussion on the y ( t ) —component of the solution of the Lienard system is connected to searching for the solution of the synthesis of filters and electrical circuits. Numerical experiments, depicting our outcomes using CAS MATHEMATICA, are presented.

1. Introduction

Hilbert [1] proposed 23 mathematical problems, of which the second part of the 16th one is to find the maximal number of limit cycles and their relative locations for polynomial vector fields. A large part of the scientific articles is concerned with the Lienard system x = y , y = g ( x ) + ϵ f ( x ) , where ϵ is a small parameter and f ( x ) and g ( x ) are polynomials of degree m and n.
The Melnikov function [2] for the Lienard system [3]
d x d t = y ϵ a 1 x + a 2 x 2 + + a 2 n + 1 x 2 n + 1 d y d t = x
is defined as
M ( α , μ ) = 2 π α 2 a 1 2 + 3 8 a 3 α 2 + + 2 n + 2 n + 1 a 2 n + 1 2 2 n + 2 α 2 n .
The Melnikov polynomial is defined as
P ( r 2 , n ) = 1 2 π r 2 M ( r , μ ) .
Perko and collaborators [4,5] ensure the needful information about radii and the number of limit cycles. More precisely, the system (1) for small enough  ϵ 0 has at most nlimit cycles asymptotic to circles of radii r j , j = 1 , 2 , , n as ϵ 0 if and only if then-th degree polynomial in r 2 ,
P ( r 2 , n ) = a 1 2 + 3 8 a 3 r 2 + + 2 n + 2 n + 1 a 2 n + 1 2 2 n + 2 r 2 n
has n positive roots r 2 = r j 2 , j = 1 , 2 , , n .
Denote by U n the second kind Chebyshev’s polynomial. The polynomials participate in the solution to some of the practical tasks [6]. In this article, we explore a new extended Lienard-type planar system with the polynomial U n . The number and type of limit cycles is also investigated. The discussion (in Section 2.3) is on y ( t ) —the component of the solution of Lienard system that is connected to the solution of some practical issues such as the synthesis of “adaptive” filters and electrical circuits. Numerical examples, confirming the results we obtained through CAS MATHEMATICA, are received.

2. Main Results–Simulations

2.1. Extended Lienard-Type Planar System

The next model will be investigated:
d x d t = y ϵ U n ( x ) d y d t = x
where ϵ > 0 and U n ( x ) for n = 3 , 5 , 7 , 9 , is the Chebyshev’s polynomial of the second kind. For example, for n = 5 , 7 , 9 , 13 we have (see Figure 1)
U 5 ( x ) = 32 x 5 32 x 3 + 6 x U 7 ( x ) = 128 x 7 192 x 5 + 80 x 3 8 x U 9 ( x ) = 512 x 9 1024 x 7 + 672 x 5 160 x 3 + 10 x
Under stated circumstances it can be indicated that the Lienard’s-type system (1) has a limit cycle. The proof of existence of the limit cycle is based on the check of the circumstances in Lienard’s theorem [3], and we will skip it. The process for user-initiated coefficient ϵ = 0.006 and n = 5 , with the model (5) for x 0 = 0.6 , y 0 = 0.4 is shown in Figure 2. For ϵ = 0.001 and n = 5 and x 0 = 0.7 , y 0 = 0.3 see Figure 3. The simulation for ϵ = 0.006 and n = 7 , with the model (5) for x 0 = 0.6 , y 0 = 0.4 , is shown in Figure 4. For ϵ = 0.003 and n = 7 ; x 0 = 0.5 , y 0 = 0.5 , see Figure 5.

2.2. The New Model in the Light of Melnikov’s Considerations

The case n = 5 . Let us explore the model
d x d t = y ϵ ( 32 x 5 32 x 3 + μ x ) d y d t = x
where μ > 0 , ϵ > 0 . The following is correct:
Proposition 1. 
The Lienard-type system (6) for n = 5 , and for all small enough ϵ 0 for μ = 7.2 has a limit cycle with multiplicity–2.
Proof. 
For the Melnikov polynomial in r 2 (see Figure 6), we have:
P ( r 2 , 2 ) = μ 2 12 r 2 + 10 r 4 .
The roots are
6 36 5 μ 10 ; 3 5 + 1 10 36 5 μ .
The radii r i ; i = 1 , 2 for some values of μ are tabulated in Table 1.
Evidently, for example μ = 7.2 we obtain a limit cycle with multiplicity–2 (see Table 1). This ends the proof of the proposition. The solution of the system (6) for ϵ = 0.001 , n = 5 , μ = 7.2 with x 0 = 0.7 , y 0 = 0.3 is visualized on Figure 7.
The case n = 7 . Let us explore the model
d x d t = y ϵ ( 128 x 7 192 x 5 + 80 x 3 μ x ) d y d t = x
where μ > 0 , ϵ > 0 . The following is correct:
Proposition 2. 
The Lienard-type system (8) for n = 7 , and for all small enough ϵ 0 for μ = 7.00805257 has a simple limit cycle and limit cycle with multiplicity–2.
Proof. 
The Melnikov polynomial in r 2 (see Figure 8) gives:
P ( r 2 , 3 ) = μ 2 + 30 r 2 60 r 4 + 35 r 6 .
The radii r i ; i = 1 , 2 , 3 for some values of μ are tabulated in Table 2.
For an proper choice of the parameter μ (for example μ = 7.00805257 ), we have a limit cycle and limit cycle with multiplicity 2 (see Table 2). The solution of the system (8) for ϵ = 0.001 , n = 7 , μ = 7.00805257 with x 0 = 0.6 , y 0 = 0.4 is visualized in Figure 9.
The case n = 9 . Consider the model
d x d t = y ϵ ( 512 x 9 1024 x 7 + 672 x 5 160 x 3 + μ 1 x ) d y d t = x
where μ 1 > 0 , ϵ > 0 . The following is correct:
Proposition 3. 
The Lienard-type system (10) for n = 9 , and for all small enough ϵ 0 for μ 1 = 9.0482691 , has 2 simple limit cycles and limit cycle with multiplicity–2.
The proof of existence of limit cycles is based on the ideas given in Propositions 1 and 2, and we will omit it here.
We note that for the Melnikov polynomial in r 2 (see Figure 10) we obtain:
P ( r 2 , 4 ) = μ 2 60 r 2 + 210 r 4 280 r 6 + 126 r 8 .
The solution of the system (10) for ϵ = 0.001 , n = 9 , μ 1 = 9.0482691 with x 0 = 0.5 , y 0 = 0.5 is visualized in Figure 11.

2.3. Related Problems and Possible Applications

Let us explore the following Lienard system:
d x d t = y d y d t = g ( x ) + ϵ f ( x ) y
where 0 ϵ 1 .
Numerical examples
Example 1. 
The solution of the system (12) for ϵ = 0.001 , g ( x ) = U 7 ( x ) , f ( x ) = x x 3 + x 5 1 7 x 7 (see oscillator model considered in [7]) is visualized in Figure 12.
Example 2. 
The solution of the system (12) for ϵ = 0.001 , g ( x ) = U 9 ( x ) , f ( x ) = x x 3 + x 5 1 7 x 7 is visualized on Figure 13.
We will explicitly note that the y ( t ) -components of the differential systems discussed above can be used successfully in modeling and approximating functions and point sets in the field of signal theory, filters, and other characteristics.
Applications
We will consider 2 typical examples.
(A) Example 1 shows that the component y ( t ) of the solution of system (12) (with x 0 = 1.1 , y 0 = 0.1 ) can be used for the modeling and synthesis of electric circuits of the type shown in Figure 14 at an appropriately selected interval.
(B) A number of studies are known on the design of in some sense optimal filters, such as the exponential optimal filters of Tadmore and Tanner [8,9] (for the intrinsic properties of these filters concerning Hausdorff metrics, see [10]).
It is easy to take into account that the alteration of the variable t with t = a cos θ in the y-component of the solution of the system (12) leads to a diagram characteristic of a filter.
By varying the parameter a, a variety of interesting models are obtained (see Figure 15).

3. Concluding Remarks

The reader can pin down the adequate approximation results for some n. For example, for fixed U 13 ( x ) the Lienard system through Melnikov’s approach is defined as
d x d t = y ϵ μ x 448 x 3 + 4032 x 5 15360 x 7 + 28160 x 9 24576 x 11 + 8192 x 13 d y d t = x
The Melnikov polynomial is of the form (see Figure 16):
P ( r 2 , 6 ) = μ 2 168 r 2 + 1260 r 4 4200 r 6 + 6930 r 8 5544 r 10 + 1716 r 12
and the following statement is valid:
Proposition 4. 
The Lienard-type system for μ = 13.25125164 has four simple limit cycles: 0.254682 , 0.521069 , 0.61736 , and 0.992178 and limit cycle 0.874307 with multiplicity–2.
The topics in this paper are connected to solving some practical issues such as the synthesis of antennas and electrical circuits. For example, a simple circuit for testing the memristor model and corresponding current–voltage characteristics is given in [11], see Figure 17. For other results, see [11,12,13,14,15,16].
Intriguing radiation diagrams can be received with a proper selection of the functions f ( x ) and g ( x ) arising in the Lienard system. Interesting outcomes are obtained in [7,17,18,19,20].
Remark 1. 
The study of dynamical systems includes bifurcation theory with branch catastrophe theory [21]. Arnold [22] discussed the catastrophes of the ADE classification, because of their relation with the Lie groups. Consider the following model in the light of Zeeman’s approach [23]:
d x d t = c ( F ( x ) y ) d y d t = 1 c x
with c > 0 and
F ( x ) = p x 448 x 3 + 4032 x 5 15360 x 7 + 28160 x 9 24576 x 11 + 8192 x 13 .
For the model the catastrophe surface ( x , y , p ) = F ( x ) y is depicted in Figure 18.

3.1. Numerical Issues Connected to the Investigation of the Roots of the Melnikov Polynomial

Let D ( λ ) be the polynomial with real and multiple roots λ 1 , , λ m with the multiplicities s i , i = 1 , , m , where i = 1 m s i = n , i.e.,
D ( λ ) = j = 1 m ( λ λ j ) s j .
The next iteration process can be used [24]:
λ i k + 1 = λ i k s i H ( λ i k ) l i m s l λ i k λ l k Δ l R , k , i = 1 , , m ; k = 0 , 1 , 2 , .
where
H ( λ i k ) = D ( λ i k ) D ( λ i k ) Δ t R , k = s i H ( λ i k ) q t m s q λ t k λ q k Δ q R 1 , k 1 , t = 1 , 2 , , m Δ t 0 , k = 0 , t = 1 , 2 , , m ; k = 0 , 1 , 2 , .
The next theorem is fulfilled
Theorem 1 
([24]). Let 0 < q < 1 , d = min i j | λ i λ j | and c > 0 be numbers such that
d > c ( 2 + 3 q n / s ) , c 2 ( n s ) ( d c ) ( d 2 c s q n / s ) 1 ( n / s + 1 ) ( n / s 1 ) c 2 q 2 ( d c ) ( d 2 c s q n / s ) 1 < 1 , s = max 1 j m s j ; s = min 1 j m s j .
If the initial approximations { λ i 0 } , i = 1 , 2 , , m of the zeros { λ i } ; i = 1 , 2 , , m satisfy the inequalities | λ i 0 λ i | c q , i = 1 , 2 , , m , then the inequalities
| λ i k λ i | c q ( 2 R + 3 ) k , i = 1 , 2 , , m ; k = 0 , 1 , 2 ,
hold true.
For other results see [25,26,27,28,29].

3.2. The Level Curves

For more details of the existing important results on the generalized polynomial Lienard differential systems and the limit-cycle bifurcations of some generalized polynomial Lienard systems, see [30,31,32,33,34,35,36,37,38,39,40,41,42,43,44,45,46,47].
Consider the class of Lienard polynomial systems of the type
d x d t = y d y d t = U 7 ( x ) + ϵ ( a x + b x 2 + c x 4 + d x 6 ) y
where 0 ϵ < 1 ; (a) U 7 ( x ) = 128 x 7 192 x 5 + 80 x 3 8 x is the Chebyshev polynomial of the second kind and a , b , c , d are bounded parameters. Without going into details, we will note some interesting level curves. For ϵ = 0 , system (14) is a Hamiltonian system with Hamiltonian
H ( x , y ) = y 2 2 16 x 8 + 32 x 6 20 x 4 + 4 x 2 .
The level curves L h i = { H ( x , y ) = h i } are depicted in Figure 19. The case b) U 9 ( x ) = 512 x 9 1024 x 7 + 672 x 5 160 x 3 + 10 x . The Hamiltonian of system (14) ( ϵ = 0 ) is
H ( x , y ) = y 2 2 5 x 2 + 40 x 4 112 x 6 + 128 x 8 256 x 10 5 .
The level curves L h i = { H ( x , y ) = h i } are depicted in Figure 20.
The outcomes received in this paper are based on the following algorithms:
(i) For the arbitrary order n fixed by the user—the production of the extended oscillator;
(ii) For the existence of the boundary cycle—the automatic checking of the rules in Lienard’s theorem;
(iii) For determining their number and type—receiving sure evaluations for the roots of the polynomial and the high-order Melnikov polynomials P ( r 2 , n ) ;
(iv) For embedding polynomial-type rectification factors—generating dynamic models.
Software for study and depiction have also been produced, written in CAS Mathematica. These modules are improving similar ones realized in computer algebraic systems designed for scientific calculations. The offered modules are only part of the more general project for investigating nonlinear models.

Author Contributions

Conceptualization, N.K.; methodology, N.K.; software, A.I.; validation, N.K. and A.I.; formal analysis, N.K.; investigation, N.K.; resources, A.I.; data curation, N.K.; writing—original draft preparation, N.K.; writing—review and editing, A.I.; visualization, A.I.; supervision, N.K.; project administration, N.K.; and funding acquisition, A.I. All authors have read and agreed to the published version of the manuscript.

Funding

This work has been accomplished with the financial support of grant no BG05M2OP001-1.001-0003, financed by the Science and Education for Smart Growth Operational Program (2014–2020) and co-financed by the European Union through the European Structural and Investment Funds.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors would like to express their gratitude to the anonymous reviewers for their helpful comments and to the editors for suggestions that improved the quality and presentation of the paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Hilbert, D. Mathematical problems (M.Newton Transl.). Bull. Am. Math. Soc. 1901, 8, 437–479. [Google Scholar] [CrossRef] [Green Version]
  2. Melnikov, V.K. On the stability of a center for time–periodic perturbation. Trudy Moskovskogo Matematicheskogo Obshchestva 1963, 12, 3–52. [Google Scholar]
  3. Lienard, A. Etude des oscillations entretenues. Revue Generale de e’Electricite 1828, 23, 901–912, 946–954. [Google Scholar]
  4. Blows, T.; Perko, L. Bifurcation of Limit Cycles from Centers and Separatrix Cycles of Planar Analytic Systems. SIAM Rev. 1994, 36, 341–376. [Google Scholar] [CrossRef]
  5. Perko, L. Differential Equations and Dynamical Systems; Springer: New York, NY, USA, 1991. [Google Scholar]
  6. Kyurkchiev, N.; Andreev, A. Approximation and Antenna and Filters Synthesis. Some Moduli in Programming Environment Mathematica; LAP LAMBERT Academic Publishing: Saarbrucken, Germany, 2014; 150p, ISBN 978-3-659-53322-8. [Google Scholar]
  7. Kyurkchiev, V.; Kyurkchiev, N. On an extended relaxation oscillator model: Number of limit cycles, simulations. I. Commun. Appl. Anal. 2022, 26, 19–42. [Google Scholar]
  8. Tadmor, E.; Tanner, J. Adaptive filters for piecewise smooth spectral data. IMA J. Numer. Alg. 2005, 25, 535–547. [Google Scholar] [CrossRef] [Green Version]
  9. Tanner, J. Optimal filter and mollifier for piecewise smooth spectral data. Math. Comput. 2006, 75, 767–790. [Google Scholar] [CrossRef] [Green Version]
  10. Kyurkchiev, N. Some Intrinsic Properties of Tadmor-Tanner Functions: Related Problems and Possible Applications. Mathematics 2020, 8, 1963. [Google Scholar] [CrossRef]
  11. Mladenov, V. A unified and open LTSPICE memristor model library. Electronics 2021, 10, 1594. [Google Scholar] [CrossRef]
  12. Mladenov, V.; Kirilov, S. A Neural Synapse Based on Ta2O5 Memristor. In Proceedings of the 2021 17th International Workshop on Cellular Nanoscale Networks and their Applications (CNNA), Catania, Italy, 29–30 September 2021. [Google Scholar] [CrossRef]
  13. Mladenov, V.; Kirilov, S. A Simplified Tantalum Oxide Memristor Model, Parameters Estimation and Application in Memory Crossbars. Technologies 2022, 10, 6. [Google Scholar] [CrossRef]
  14. Mladenov, V.M.; Zaykov, I.D.; Kirilov, S.M. Application of a Nonlinear Drift Memristor Model in Analogue Reconfigurable Devices. In Proceedings of the 2022 26th International Conference Electronics, Palanga, Lithuania, 13–15 June 2022; pp. 1–6. [Google Scholar]
  15. Mladenov, V. Advanced Memristor Modeling–Memristor Circuits and Networks; MDPI: Basel, Switzerland, 2019. [Google Scholar]
  16. Mladenov, V. Analysis and Simulations of Hybrid Memory Scheme Based on Memristors. Electronics 2018, 7, 289. [Google Scholar] [CrossRef] [Green Version]
  17. Kyurkchiev, V.; Iliev, A.; Rahnev, A.; Kyurkchiev, N. A technique for simulating the dynamics of some extended relaxation oscillator models. II. Commun. Appl. Anal. 2022, 26, 43–59. [Google Scholar]
  18. Kyurkchiev, V.; Iliev, A.; Rahnev, A.; Kyurkchiev, N. Another extended polynomial Lienard systems: Simulations and applications. III. Int. Electron. J. Pure Appl. Math. 2022, 16, 55–65. [Google Scholar]
  19. Kyurkchiev, V.; Iliev, A.; Rahnev, A.; Kyurkchiev, N. Investigations on some polynomial Lienard-type systems: Number of limit cycles, simulations. Int. J. Differ. Equ. Appl. 2022, 21, 117–126. [Google Scholar]
  20. Kyurkchiev, V.; Kyurkchiev, N.; Iliev, A.; Rahnev, A. On Some Extended Oscillator Models: A Technique for Simulating and Studying Their Dynamics; Plovdiv University Press: Plovdiv, Bulgaria, 2022; ISBN 978-619-7663-13-6. [Google Scholar]
  21. Beltrami, E. Mathematics for Dynamic Modeling; Academic Press: Boston, MA, USA, 1987. [Google Scholar]
  22. Arnold, V. Geometrical Methods in the Theory of Ordinary Differential Equation; Springer: Berlin/Heidelberg, Germany, 1988. [Google Scholar]
  23. Zeeman, E. Catastrophe Theory. Selected Papers 1972–1977; Addison-Wesley: Reading, MA, USA, 1977. [Google Scholar]
  24. Kyurkchiev, N.; Andreev, A.; Popov, V. Iterative methods for the computation of all multiple roots of an algebraic polynomial. Annu. Univ. Sofia Fac. Math. Mech. 1984, 78, 178–185. [Google Scholar]
  25. Proinov, P.; Vasileva, M. On the convergence of high-order Ehrlich-type iterative methods for approximating all zeros of a polynomial simultaneously. J. Inequalities Appl. 2015, 2015, 336. [Google Scholar] [CrossRef]
  26. Proinov, P.D.; Vasileva, M.T. Local and Semilocal Convergence of Nourein’s Iterative Method for Finding All Zeros of a Polynomial Simultaneously. Symmetry 2020, 12, 1801. [Google Scholar] [CrossRef]
  27. Kyncheva, V.; Yotov, V.; Ivanov, S. Convergence of Newton, Halley and Chebyshev iterative methods as methods for simultaneous determination of multiple zeros. Appl. Numer. Math. 2017, 112, 146–154. [Google Scholar] [CrossRef]
  28. Petkovic, I.; Herceg, D. Computer methodologies for comparison of computational efficiency of simultaneous methods for finding polynomial zeros. J. Comput. Appl. Math. 2020, 368, 112513. [Google Scholar] [CrossRef]
  29. Kanno, S.; Kjurkchiev, N.; Yamamoto, T. On some methods for the simultaneous determination of polynomial zeros. Jpn. J. Appl. Math. 1995, 13, 267–288. [Google Scholar] [CrossRef]
  30. Llibre, J.; Valls, C. Global centers of the generalized polynomial Lienard differential systems. J. Differ. Equ. 2022, 330, 66–80. [Google Scholar] [CrossRef]
  31. Chen, H.; Lie, Z.; Zhang, R. A sufficient and necessary condition of generalized polynomial Lienard systems with global centers. arXiv 2022, arXiv:2208.06184. [Google Scholar]
  32. He, H.; Llibre, J.; Xiao, D. Hamiltonian polynomial differential systems with global centers in the plane. Sci. China Math. 2021, 48, 2018. [Google Scholar]
  33. Smale, S. Mathematical problems for the next century. Math. Intell. 1998, 20, 7–15. [Google Scholar] [CrossRef]
  34. Zhao, Y.; Liang, Z.; Lu, G. On the global center of polynomial differential systems of degree 2k + 1. In Differential Equations and Control Theory; CRC Press: Boca Raton, FL, USA, 1996; 10p. [Google Scholar]
  35. Andronov, A.A.; Leontovich, E.A.; Gordon, I.I.; Maier, A.G.; Gutzwiller, M.C. Qualitative Theory of Second Order Dynamic Systems; Wiley: New York, NY, USA, 1973. [Google Scholar]
  36. Hale, J.K. Ordinary Differential Equations; Wiley: New York, NY, USA, 1980. [Google Scholar]
  37. Garcia, I.A. Cyclicity of Nilpotent Centers with Minimum Andreev Number. 2019. Available online: https://repositori.udl.cat/handle/10459.1/67895 (accessed on 1 October 2022).
  38. Sun, X.; Xi, H. Bifurcation of limit cycles in small perturbation of a class of Lienard systems. Int. J. Bifurc. Chaos 2014, 24, 1450004. [Google Scholar] [CrossRef]
  39. Asheghi, R.; Bakhshalizadeh, A. On the distribution of limit cycles in a Lienard system with a nilpotent center and a nilpotent saddle. Int. J. Bifurc. Chaos 2016, 26, 1650025. [Google Scholar] [CrossRef]
  40. Zaghian, A.; Kazemi, R.; Zangenech, H. Bifurcation of limit cycles in a class of Lienard system with a cusp and nilpotent saddle. UPB Sci. Bull. Ser. A 2016, 78, 95–106. [Google Scholar]
  41. Gaiko, V.; Vuik, C.; Reijm, H. Bifurcation Analysis of Multi–Parameter Lienard Polynomial System. IFAC-PapersOnLine 2018, 51, 144–149. [Google Scholar] [CrossRef]
  42. Cai, J.; Wei, M.; Zhu, H. Nine limit cycles in a 5-degree polynomials Lienard system. Complexity 2020, 2020, 8584616. [Google Scholar] [CrossRef]
  43. Xu, W.; Li, C. Limit cycles of some polynomial Lienard system. J. Math. Anal. Appl. 2012, 389, 367–378. [Google Scholar] [CrossRef]
  44. Xu, W. Limit cycle bifurcations of some polynomial Lienard system with symmetry. Nonlinear Anal. Differ. Equ. 2020, 8, 77–87. [Google Scholar] [CrossRef]
  45. Hou, J.; Han, M. Melnikov functions for planar near–Hamiltonian systems and Hopf bifurcations. J. Shanghai Norm. Univ. (Nat. Sci.) 2006, 35, 1–10. [Google Scholar]
  46. Han, M.; Yang, J.; Tarta, A.; Gao, Y. Limit cycles near homoclinic and heteroclinic loops. J. Dyn. Differ. Equ. 2008, 20, 923–947. [Google Scholar] [CrossRef]
  47. An, Y.; Han, M. On the number of limit cycles near a homoclinic loop with a nilpotent singular point. J. Differ. Equ. 2015, 258, 3194–3247. [Google Scholar] [CrossRef]
Figure 1. U n ( x ) for n = 5 , n = 7 and n = 9 .
Figure 1. U n ( x ) for n = 5 , n = 7 and n = 9 .
Algorithms 15 00462 g001
Figure 2. The system solutions ( ϵ = 0.006 ; n = 5 ).
Figure 2. The system solutions ( ϵ = 0.006 ; n = 5 ).
Algorithms 15 00462 g002
Figure 3. The system solutions ( ϵ = 0.001 ; n = 5 ).
Figure 3. The system solutions ( ϵ = 0.001 ; n = 5 ).
Algorithms 15 00462 g003
Figure 4. The system solutions ( ϵ = 0.006 ; n = 7 ).
Figure 4. The system solutions ( ϵ = 0.006 ; n = 7 ).
Algorithms 15 00462 g004
Figure 5. The system solutions ( ϵ = 0.003 ; n = 7 ).
Figure 5. The system solutions ( ϵ = 0.003 ; n = 7 ).
Algorithms 15 00462 g005
Figure 6. P ( r 2 , 2 ) for n = 5 and μ = 7.2 .
Figure 6. P ( r 2 , 2 ) for n = 5 and μ = 7.2 .
Algorithms 15 00462 g006
Figure 7. (a) The solutions of the system (6) ( ϵ = 0.001 ; n = 5 ; μ = 7.2 ); (b) the portrait ( x , y ) ; (c) 3D plot.
Figure 7. (a) The solutions of the system (6) ( ϵ = 0.001 ; n = 5 ; μ = 7.2 ); (b) the portrait ( x , y ) ; (c) 3D plot.
Algorithms 15 00462 g007
Figure 8. P ( r 2 , 3 ) for n = 7 and μ = 7.00805257 . The roots are: 0.409106 —simple and 0.879462 —with multiplicity 2.
Figure 8. P ( r 2 , 3 ) for n = 7 and μ = 7.00805257 . The roots are: 0.409106 —simple and 0.879462 —with multiplicity 2.
Algorithms 15 00462 g008
Figure 9. (a) The solutions of the system (8) ( ϵ = 0.001 ; n = 7 ; μ = 7.00805257 ); (b) the portrait ( x , y ) ; (c) 3D plot.
Figure 9. (a) The solutions of the system (8) ( ϵ = 0.001 ; n = 7 ; μ = 7.00805257 ); (b) the portrait ( x , y ) ; (c) 3D plot.
Algorithms 15 00462 g009
Figure 10. P ( r 2 , 4 ) for n = 9 and μ 1 = 9.0482691 . The simple roots are 0.338954 , 0.987437 and root 0.75243 —with multiplicity 2.
Figure 10. P ( r 2 , 4 ) for n = 9 and μ 1 = 9.0482691 . The simple roots are 0.338954 , 0.987437 and root 0.75243 —with multiplicity 2.
Algorithms 15 00462 g010
Figure 11. The solutions of the system (10) ( ϵ = 0.001 ; n = 9 ; μ 1 = 9.0482691 ).
Figure 11. The solutions of the system (10) ( ϵ = 0.001 ; n = 9 ; μ 1 = 9.0482691 ).
Algorithms 15 00462 g011
Figure 12. The solutions of the system (12) (Example 1).
Figure 12. The solutions of the system (12) (Example 1).
Algorithms 15 00462 g012
Figure 13. The solutions of the system (12) (Example 2).
Figure 13. The solutions of the system (12) (Example 2).
Algorithms 15 00462 g013
Figure 14. Application for modeling and analysis of electrical circuits and signal processes.
Figure 14. Application for modeling and analysis of electrical circuits and signal processes.
Algorithms 15 00462 g014
Figure 15. The model y ( θ ) for x 0 = 0.7 , y 0 = 0.3 in ( π 2 , π 2 ) for (a) a = 0.4 ; (b) a = 0.49 ; and (c) a = 0.5 .
Figure 15. The model y ( θ ) for x 0 = 0.7 , y 0 = 0.3 in ( π 2 , π 2 ) for (a) a = 0.4 ; (b) a = 0.49 ; and (c) a = 0.5 .
Algorithms 15 00462 g015
Figure 16. The Melnikov polynomial P ( r 2 , 6 ) .
Figure 16. The Melnikov polynomial P ( r 2 , 6 ) .
Algorithms 15 00462 g016
Figure 17. A simple circuit for testing the memristor model and corresponding current–voltage characteristics [11].
Figure 17. A simple circuit for testing the memristor model and corresponding current–voltage characteristics [11].
Algorithms 15 00462 g017
Figure 18. The catastrophe surface ( x , y , p ) = F ( x ) y ; p = 0.1 , 10 , 20 .
Figure 18. The catastrophe surface ( x , y , p ) = F ( x ) y ; p = 0.1 , 10 , 20 .
Algorithms 15 00462 g018
Figure 19. The level curves (the case a).
Figure 19. The level curves (the case a).
Algorithms 15 00462 g019
Figure 20. Level curves (the case b).
Figure 20. Level curves (the case b).
Algorithms 15 00462 g020
Table 1. The radii r i ; i = 1 , 2 for some values of μ .
Table 1. The radii r i ; i = 1 , 2 for some values of μ .
μ r 1 r 2
50.5180130.965227
60.5958620.919211
70.7071070.83666
7.20.7745970.774597
Table 2. The values of μ that proceed in 3 positive zeros of the Melnikov polynomial (Theorem 2); for μ = 7.00487909525 : a simple limit cycle and limit cycle with multiplicity–2.
Table 2. The values of μ that proceed in 3 positive zeros of the Melnikov polynomial (Theorem 2); for μ = 7.00487909525 : a simple limit cycle and limit cycle with multiplicity–2.
μ r 1 r 2 r 3
7.50.4333910.8082630.934435
7.30.4232550.8268620.922735
7.10.4134690.8515450.904544
7.0088052570.4091060.8784620.879469
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kyurkchiev, N.; Iliev, A. On a Hypothetical Model with Second Kind Chebyshev’s Polynomial–Correction: Type of Limit Cycles, Simulations, and Possible Applications. Algorithms 2022, 15, 462. https://doi.org/10.3390/a15120462

AMA Style

Kyurkchiev N, Iliev A. On a Hypothetical Model with Second Kind Chebyshev’s Polynomial–Correction: Type of Limit Cycles, Simulations, and Possible Applications. Algorithms. 2022; 15(12):462. https://doi.org/10.3390/a15120462

Chicago/Turabian Style

Kyurkchiev, Nikolay, and Anton Iliev. 2022. "On a Hypothetical Model with Second Kind Chebyshev’s Polynomial–Correction: Type of Limit Cycles, Simulations, and Possible Applications" Algorithms 15, no. 12: 462. https://doi.org/10.3390/a15120462

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop