Abstract
In this paper, we study the probability of occurrence of phase portraits in the set of random planar homogeneous polynomial vector fields, of degree n. In particular, for \(n=1,2,3,\) we give the complete solution of the problem; that is, we either give the exact value of each probability of occurrence or we estimate it by using the Monte Carlo method. Remarkably is that all but two of these phase portraits are characterized by the index at the origin and by the number of invariant straight lines through this point.
Similar content being viewed by others
Data availability
All data analyzed in this study are included in this article.
References
Argèmi, J.: Sur les points singuliers multiples de systèmes dynamiques dans \(\mathbb{R}^2\). Ann. Mat. Pure Appl. 79, 35–69 (1968)
Arnold, V.I.: Eksperimental’naya matematika (Experimental Mathematics). Fazis, Moscow (2005)
Arnold, V.I., Varchenko, A.N., Gusein-Zade, S.M.: Singularities of Differentiable Maps: Volume I: The Classification of Critical Points Caustics and Wave Fronts. Springer, New York (2012)
Artés, J.C., Llibre, J.: Statistical measure of quadratic vector fields. Resenhas 6, 85–97 (2003)
Artés, J.C., Llibre, J., Schlomiuk, D., Vulpe, N.: From topological to geometric equivalence in the classification of singularities at infinity for quadratic vector fields. Rocky Mt. J. Math. 45, 29–113 (2015)
Artés, J.C., Llibre, J., Schlomiuk, D., Vulpe, N.: Global topological configurations of singularities for the whole family of quadratic differential systems. Qual. Theory Dyn. Syst. 19, 32 (2020)
Bratley, P., Fox, B.L., Schrage, L.E.: A Guide to Simulation, 2d edn. Springer, New York (1987)
Chen-Charpentier, B.M., Stanescu, D.: Epidemic models with random coefficients. Math. Comput. Model. 52, 1004–1010 (2010)
Cima, A., Gasull, A., Mañosa, V.: Stability index of linear random dynamical systems. Electron. J. Qual. Theory Differ. Equ. arXiv:1904.05725 [math.DS] (to appear)
Cima, A., Gasull, A., Torregrosa, J.: On the relation between index and multiplicity. J. Lond. Math. Soc. 57, 757–768 (1998)
Cima, A., Llibre, J.: Algebraic and topological classification of the homogeneous cubic vector fields in the plane. J. Math. Anal. Appl. 747, 420–448 (1990)
Conrad, K.: Probability distributions and maximum entropy. https://kconrad.math.uconn.edu/blurbs/analysis/entropypost.pdf. Accessed 29 Nov 2020
Date, T.: Classification and analysis of two-dimensional real homogeneous quadratic differential equation systems. J. Differ. Equ. 32, 311–334 (1979)
Dumortier, F., Artés, J.C., Llibre, J.: Qualitative Theory of Planar Differential Systems. Springer, Berlin (2006)
Edelman, A., Kostlan, E.: How many zeros of a random polynomial are real? Bull. Am. Math. Soc. 32, 1–37 (1995)
Eisenbud, D., Levine, I.: An algebraic formula for the degree of a \(C^\infty \) map germ. Ann. Math. 106, 19–44 (1977)
Grigoriu, M., Soong, T.: Random Vibration in Mechanical and Structural Systems. Prentice-Hall, Englewood Cliffs (1993)
Leonov, G.A., Burova, I.G., Aleksandrov, K.D.: Visualization of four limit cycles of two-dimensional quadratic systems in the parameter space. Differ. Equ. 49, 1675–1703 (2013)
Llibre, J., Pérez del Río, J.S., Rodríguez, J.A.: Structural stability of planar homogeneous polynomial vector fields: applications to critical points and to infinity. J. Differ. Equ. 125, 490–520 (1996)
Lopes, H., Pagnoncelli, B., Palmeira, C.: Coeficientes aleatórios de equações diferenciais ordinárias lineares. Mat. Univ. 14, 44–50 (2008). (English version [25])
Marsaglia, G.: Choosing a point from the surface of a sphere. Ann. Math. Stat. 43, 645–646 (1972)
Milnor, J.: Topology from a Differentiable Viewpoint. The University Press of Virginia, Charlottesville (1965)
Morgan, B.J.T.: Applied Stochastic Modelling. Arnold Publishers, Oxford (2000)
Muller, M.E.: A note on a method for generationg points uniformly on \(N\)-dimensional spheres. Commun. ACM 2, 19–20 (1959)
Pagnoncelli, B., Lopes, H., Palmeira, C., Palmeira, B.: Random linear systems and simulation. Preprint MAT. 18/06, Pontifícia Universidade Católica de Rio de Janeiro. http://bernardokp.uai.cl/preprint.pdf. Accessed 29 Nov 2020
Stanescu, D., Chen-Charpentier, B.M.: Random coefficient differential equation models for bacterial growth. Math. Comput. Model. 50, 885–895 (2009)
Stoer, J., Bulirsch, R.: Introduction to Numerical Analysis. Springer, New York (2002)
Strogatz, S.H.: Nonlinear Dynamics and Chaos. Westview Press, Cambridge (1994)
Acknowledgements
We want to acknowledge Prof. Maria Jolis for her helpful indications.
We also acknowledge the anonymous reviewer for the careful reading and for addressing us very interesting comments. Some of his/her indications are summarized in Remark 13. We also thank him/her for sharing with us the results of his/her Monte Carlo simulations using the approximation of that remark, which are consistent with our results.
Author information
Authors and Affiliations
Corresponding author
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
The authors are supported by Ministry of Science and Innovation–State Research Agency of the Spanish Government through Grants PID2019-104658GB-I00 (MICINN/AEI, Anna Cima and Armengol Gasull) and DPI2016-77407-P (MICINN/AEI/FEDER, UE, Víctor Mañosa). The first and second authors are also supported by the Grant 2017-SGR-1617 from AGAUR, Generalitat de Catalunya. The third author acknowledges the group’s research recognition 2017-SGR-388 from AGAUR, Generalitat de Catalunya.
Appendix
Appendix
In this appendix we give conditions on the coefficients of polynomial equations of degree 3 or 4 to know its number real zeros zeros and for degree 3 also their signs.
Lemma 15
Set \(p(x)=x^3+ax^2+bx+c,\) \(c_2=ab-9c,\) \(d_2=a^2-3b\) and \(d_3=-27c^2+18abc-4a^3c+a^2b^2-4b^3.\) Assume that \(a\cdot b\cdot c \cdot c_2\cdot d_2\cdot d_3\ne 0.\) Then the following holds:
-
(i)
p has a unique real root of multiplicity one, if and only if \(d_3<0.\) This root is negative (resp. positive) if \(c>0\) (resp. \(c<0\)).
-
(ii)
p has three negative real roots if and only if \(d_3>0,\,c_2>0,\,b>0,c>0.\)
-
(iii)
p has two negative roots and one positive one if and only if either, \(d_3>0,\,c<0,c_2>0\) or \(d_3>0,\,c<0,\,c_2<0,b<0.\)
-
(iv)
p has one negative root and two positive roots if and only if either, \(d_3>0,\,c>0,c_2<0\) or \(d_3>0,\,c>0,\,c_2>0,b<0.\)
-
(v)
p has three positive roots if and only if \(d_3>0,\,c_2<0,\,b>0,c<0.\)
Proof
The proof is based in Sturm’s method which asserts that if \(\big (p_0=p,p_1\ldots ,p_n\big )\) is a Sturm’s sequence of p in [a, b] with \(p(a)\cdot p(b)\ne 0,\) then the number of real zeros of p in (a, b) is \(V(a)-V(b)\) where V(x) is the number of changes of sign in the ordered sequence \(\big (p_0(x),p_1(x),\ldots ,p_n(x)\big ),\) where the zeroes are disregarded. For any polynomial without multiple roots such a sequence always exists, see [27]. For p, without multiple roots and \(d_2\ne 0,\) one Sturm sequence is \(p_0=p,\) \(p_1=p'\) and
The quantity \(d_3\) is the classical discriminant of p and it is known that \(d_3<0\) if and only if p has a real root and two complex ones whereas \(d_3>0\) if and only if p has three real roots (see [11] for instance). We are going to consider the following two tables depending on the sign of \(d_3\) where \(b,c,d_2,c_2\) stands for their respective signs.
\(-\infty \) | 0 | \(\infty \) | \(-\infty \) | 0 | \(\infty \) | ||
---|---|---|---|---|---|---|---|
\(d_3<0\) | \(d_3>0\) | ||||||
\(p_0\) | − | c | \(+\) | \(p_0\) | − | c | \(+\) |
\(p_1\) | \(+\) | b | \(+\) | \(p_1\) | \(+\) | b | \(+\) |
\(p_2\) | \(-d_2\) | \(c_2\) | \(d_2\) | \(p_2\) | \(-d_2\) | \(c_2\) | \(d_2\) |
\(p_3\) | − | − | − | \(p_3\) | \(+\) | \(+\) | \(+\) |
We separate the proof in two cases, depending on the sign of \(d_3.\)
Case 1 Assume that \(d_3<0.\) We observe that \(V(-\infty )=2\) and \(V(+\infty )=1.\) Then
-
The polynomial p has a negative real root if and only if \(V(0)=1.\) It is easy to see that it happens when \(c>0\) and one of the three following conditions hold:
$$\begin{aligned} c_2>0,b>0;\quad c_2<0,b>0;\quad c_2<0,b<0. \end{aligned}$$We observe that \(d_3<0,\,c_2>0,\,b<0,\,c>0\) is not compatible because then \(V(0)-V(+\infty )=3\) but \(V(-\infty )-V(+\infty )=1.\) Summarizing, p has a negative real root and two more complex ones if and only if \(d_3<0,\,c>0.\)
-
The polynomial p has a positive real root if and only if \(V(0)=2.\) It is easy to see that it happens when \(c<0\) and one of the three following conditions hold:
$$\begin{aligned} c_2>0,b>0;\quad c_2>0,b<0;\quad c_2<0,b>0. \end{aligned}$$As before conditions \(d_3<0,c_2<0,b<0,c<0\) are incompatible and then p has a positive real root and two more complex ones if and only if \(d_3<0,\,c<0\) as announced.
Case 2 Assume that \(d_3>0\) and consider the above right-hand side table of signs. We see that in that case \(d_2\) must be positive. Otherwise \(V(-\infty )=1,\,V(\infty )=2\) which is not possible (we can also argue that since \(d_3>0\) implies that p has three real roots, its derivative has two real roots and hence its discriminant which is equal to \(4\,d_2\) has to be positive).
It is straightforward to see that items (ii), (iii), (iv), (v) are equivalent to \(V(0)=0,V(0)=1,V(0)=2,V(0)=3\) respectively and that these number of changes of sign are satisfied exactly when the conditions on \(b,c,c_2\) are the ones stated in the lemma. \(\square \)
Lemma 16
Let \(p(x)=x^4+ax^3+bx^2+cx+d.\) Assume that \(c\cdot d\cdot c_2\cdot c_3\cdot d_2\cdot d_3\cdot d_4\ne 0.\) Then the following holds:
-
(a)
p has two simple real and two complex roots if and only if \(d_4<0.\)
-
(b)
p has four different real roots if and only if \(d_4>0,d_2>0\) and \(d_3>0.\)
-
(c)
p has no real root if and only if \(d_4>0\) and either, \(d_2<0\) or \(d_3<0.\)
Proof
A Sturm sequence of p, is \(p_0=p,\) \(p_1=p',\)
where
We note that \(d_4\) is the discriminant of p and since \(d_4\ne 0\) all its roots are simple.
It is known, see [11] for instance, that \(d_4<0\) if and only if p has two simple real and two complex roots. And that \(d_4>0\) if and only if p has either, four different real roots or no real root. It is very easy to distinguish between these last two possibilities using the Sturm method. Consider the corresponding table when \(d_4>0,\) where again each value stands for its sign.
\(-\infty \) | 0 | \(\infty \) | |
---|---|---|---|
\(p_0\) | \(+\) | d | \(+\) |
\(p_1\) | − | c | \(+\) |
\(p_2\) | \(d_2\) | \(c_2\) | \(d_2\) |
\(p_3\) | \(-d_3\) | \(c_3\) | \(d_3\) |
\(p_4\) | \(+\) | \(+\) | \(+\) |
If p has four real roots then \(V(-\infty )-V(\infty )\) must be four and this only can happen if \(V(-\infty )=4\) and \(V(\infty )=0\) what immediately says that \(d_2\) and \(d_3\) must be positive. If \(d_4>0\) and \(d_2<0\) or \(d_3<0\) then \(V(-\infty )-V(+\infty )=2-2=0\) and the result follows. \(\square \)
Rights and permissions
About this article
Cite this article
Cima, A., Gasull, A. & Mañosa, V. Phase Portraits of Random Planar Homogeneous Vector Fields. Qual. Theory Dyn. Syst. 20, 3 (2021). https://doi.org/10.1007/s12346-020-00437-7
Received:
Accepted:
Published:
DOI: https://doi.org/10.1007/s12346-020-00437-7
Keywords
- Ordinary differential equations with random coefficients
- Planar homogeneous vector fields
- Index
- Phase portraits