• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    A class of two-dimensional rational maps with self-excited and hidden attractors

    2022-03-12 07:47:50LiPingZhang張麗萍YangLiu劉洋ZhouChaoWei魏周超HaiBoJiang姜海波andQinShengBi畢勤勝
    Chinese Physics B 2022年3期
    關(guān)鍵詞:張麗萍劉洋海波

    Li-Ping Zhang(張麗萍) Yang Liu(劉洋) Zhou-Chao Wei(魏周超)Hai-Bo Jiang(姜海波) and Qin-Sheng Bi(畢勤勝)

    1Faculty of Civil Engineering and Mechanics,Jiangsu University,Zhenjiang 212013,China

    2School of Mathematics and Statistics,Yancheng Teachers University,Yancheng 224002,China

    3College of Engineering,Mathematics and Physical Sciences,University of Exeter,Exeter EX4 4QF,UK

    4School of Mathematics and Physics,China University of Geosciences,Wuhan 430074,China

    Keywords: two-dimensional rational map, hidden attractors, multi-stability, a line of fixed points, chaotic attractor

    1. Introduction

    Recently,hidden attractor has become an attractive topic for researchers in non-linear sciences, and has been found in many practical systems, such as the Chua system,[1-3]the drilling system,[4]and the automatic control system with piecewise-linear non-linearity.[5]They have been explored intensively in continuous dynamical systems with different structures of equilibria (see Ref. [6]). As it is known, if a dynamical system with a given set of parameters has more than one attractor according to its initial conditions,this phenomenon is called multi-stability. Multi-stability has been studied extensively in the literature since it exists in many areas, such as physics, chemistry, biology and economics (see,e.g.,Refs.[7-11]).

    Many researchers have studied the self-excited and hidden attractors in discrete-time maps because of their broad applications (see e.g., Refs. [12-24]). Among these works, a few studies[15-19]have identified the attractors with particular types of fixed points. For example, Jianget al.[15,16]explored the hidden chaotic attractors with no fixed point and a single stable fixed point in a class of two-dimensional and three-dimensional maps, respectively. Jianget al.[17]studied the hidden chaotic attractors in a class of two-dimensional chaotic maps with closed curve fixed points. Luo[19]analyzed the complexity and singularity of a class of discrete systems with infinite-fixed-points by using the local analysis methods.Also,hidden attractors have been found in the non-linear maps in a fractional form.[20-25]In Ref. [20], Ouannaset al.studied a fractional map with no fixed point that exhibits hidden attractors. Hadjabiet al.proposed two new classes of twodimensional fractional maps with closed curve fixed points and studied their dynamics in Ref. [21]. Great efforts[26-30]have also been made for searching and controlling hidden attractors. For example,Dudkowskiet al.[26,27]presented a new method to locate hidden and coexisting attractors in non-linear maps based on the concept of perpetual points. In Ref. [28]Danca and Fe?can used impulsive control to suppress chaos and produced hidden attractors in a one-dimensional discrete supply and demand dynamical system. Then, Danca and Lampart[29]proposed an algorithm to identify hidden attractors in maps and numerically studied the dynamics of a heterogeneous Cournot oligopoly model exhibiting self-excited and hidden attractors. In Ref. [30], Zhanget al.put forward a method to distinguish hidden and self-excited attractors by constructing random bifurcation diagrams. They employed the linear augmentation method to control hidden attractors and multi-stability in a class of two-dimensional maps.

    On the other hand, the extreme multi-stability of which the map exhibits infinite many coexisting attractors is of great interest of many researchers.[31-33]In Ref. [31], Zhanget al.constructed a new class of two-dimensional maps by introducing a sine term and presented infinitely many coexisting attractors in the map. Baoet al.proposed a new two-dimensional hyper-chaotic map with a sine term. They investigated the parameter-dependent and initial-boosting bifurcations for the map with line and no fixed points in Ref. [32]. In Ref. [33],Konget al.proposed a novel two-dimensional hyper-chaotic map with two sine terms and investigated the map’s conditional symmetry and attractor growth. Among different maps,the maps with rational fraction (see e.g., Refs. [34-41]) are complex and it is challenging for investigation. In Ref. [35],Luet al.investigated the complex dynamics of a new rational chaotic map. Changet al.constructed a new two-dimensional rational chaotic map and studied its tracking and synchronization in Ref. [36]. Elhadj and Sprott proposed a new twodimensional rational map that presents a quasi-periodic route to chaos in Refs. [37,38]. In Ref. [39], Somarakis and Baras investigated the complex dynamics of the two-dimensional rational map proposed in Ref.[37]analytically and numerically by studying the strange attractors,bifurcation diagrams,periodic windows,and invariant characteristics. Chenet al.studied the boundedness of the attractors and the corresponding estimation of absorbing set of the Zeraoulia-Sprott map[37]analytically in Ref.[40]. In Ref.[41], Ouannaset al.considered the dynamics,control and synchronization of the rational maps in a fractional form based on the Rulkov,[34]Chang,[36]and Zeraoulia-Sprott[37]rational maps. However, to the best of authors’ knowledge, the current research work on the hidden attractors in rational maps is very limited, which is the motivation of the present work.

    This paper aims to explore several simple chaotic rational maps exhibiting self-excited and hidden attractors by performing an exhaustive computer search.[13,14]The main focus of this paper is as follows: (i)A new class of two-dimensional rational maps with self-excited and hidden attractors are developed. (ii) The chaotic attractors are analyzed numerically by using the basins of attraction,the Lyapunov exponent spectrum(Les)and the Kaplan-Yorke dimension. (iii)By varying parameters,the map can display four types of fixed points,i.e.,no fixed point, one single fixed point, two fixed points and a line of fixed points. Furthermore,it can exhibit different types of coexisting attractors.

    The rest of this paper is structured as follows. In Section 2,the mathematical model of this class of two-dimensional rational maps is formulated, and the existence and stability of their fixed points are studied. In Section 3, the dynamical behaviors of the rational map are investigated by using various numerical analysis tools. Finally,concluding remarks are drawn in Section 4.

    2. System model

    Inspired by Refs. [31-33], we construct a new class of two-dimensional rational maps in the following form:

    stand for the determinant of the Jacobian matrix and the trace of the Jacobian matrix,respectively. Eigenvalues ofJ,λ1,andλ2are termed as multipliers of the fixed point. Denote the numbers of multipliers of the fixed point(x*,y*)lying inside,on and outside the unit circle{λ ∈C:|λ|=1}byn-,n0,andn+, respectively. The fixed point is stable if the roots of the characteristic equationλ1andλ2satisfy that|λ1,2|<1,where|·|refers to the modulus of a complex number.

    Definition 1 (Definition 2.10 in Ref. [42]) A fixed point(x*,y*) is called hyperbolic ifn0=0, that is, if there is no eigenvalue of the Jacobian matrix evaluated at this fixed point on the unit circle. Otherwise, the fixed point is called nonhyperbolic,that is,there is at least one eigenvalue of the Jacobian matrix evaluated at the fixed point on the unit circle.

    To distinguish the hidden and self-excited attractors of the rational map(1),the following definition is introduced.

    Definition 2[43]An attractor is called a hidden attractor if its basin of attraction does not intersect with the small neighborhoods of equilibria(fixed points)of the system(map);otherwise,it is called a self-excited attractor.

    Five cases can be divided from Eq.(2)as follows.

    Ifa/=0,b=0,c/=0, there are two different cases, i.e.,d=0 andd/=0.

    Case I:a/=0,b=0,c/=0,d=0(a line of fixed points).

    Whena/=0,b=0,c/=0, andd=0, there is a line of fixed points(0,e),e ∈R. The Jacobian matrix of the rational map(1)at the fixed points(0,e)can be rewritten as

    The multipliers of the fixed points, eigenvalues ofJ,λ1=1 andλ2=1+a/(1+e2). Thus,according to definition 1,these fixed points are all non-hyperbolic.Then these fixed points are critical stable.

    Case II:a/=0,b=0,c/=0,d/=0(no fixed point A).

    Whena/=0,b=0,c/=0, andd/=0, Eq. (2) has none solution. So there is no fixed point.

    Ifa/=0,b/=0,andc/=0,thenx*=-d/c,-ad/[c(1+(y*)2)]+b=0,i.e.,(y*)2=(ad/bc)-1. So ifa/=0,c/=0,andb/=0,three cases can be obtained by considering the size relation between(ad/bc)and 1.

    Case III:a/=0,b/=0,c/=0,(ad/bc)<1(no fixed point B).

    Whena/=0,b/=0,c/=0, and(ad/bc)<1, Eq.(2)has none solution. So there is no fixed point.

    Case IV:a/= 0,b/= 0,c/= 0, (ad/bc) = 1 (one fixed point).

    Whena/=0,b/=0,c/=0, and(ad/bc)=1, Eq.(2)has two equal solutions,i.e.,x*=-d/c,y*=0. So there is only one fixed point. The Jacobian matrix of the rational map(1)at the fixed point(-d/c,0)can be expressed as

    The multipliers of the fixed points can also be obtained by solving the characteristic equation of the Jacobian matrix(8),which will be analyzed in the following section.

    3. Dynamical behaviors of the rational map

    There are several ways to construct a bifurcation diagram,for example, by using random or fixed initial value, and forward or backward. In the random bifurcation diagram, many initial values are selected randomly in a range for each bifurcation parameter value. So the random bifurcation diagram can be used to exhibit all possible attractors. In the forward(backward)bifurcation diagram,one initial value is chosen first and the method of “follow the attractor” is used, that is, the last steady state is used for the initial state of next value of the increasing(decreasing)bifurcation parameter. In the bifurcation diagram calculated by using the fixed initial value,only a fixed initial value is chosen for all the values of the bifurcation parameter.

    In this section, the dynamical behaviors of the rational map (1) will be investigated in five cases in accordance with different types of fixed points. In each case,the random bifurcation diagram[30]is drawn firstly to show the possible attractors of the rational map(1). If no multi-stability is observed,the forward(backward)bifurcation diagram will be presented to show the complex behaviors of the rational map(1). Moreover, the Lyapunov exponent spectrum (Les) and the Lyapunov (Kaplan-Yorke) dimension (Dky) of the attractors of the rational map(1)are computed by using the Wolf methods given in Refs.[14,44].

    3.1. Case I:a line of fixed points

    By random searching,chaotic attractors were found with several coefficients satisfying the condition thata/=0,b=0,c/=0, andd=0, i.e., there exist a line of fixed points. For example,whena=-3,b=0,c=0.1,andd=0,the rational map(1)can present chaotic attractors.

    In order to show the complex dynamics of the rational map (1) with a line of fixed points, random bifurcation diagram and Lyapunov exponent spectrum diagram of the map were plotted in Fig. 1 by usingaas a branching parameter and fixing (b,c,d) as (0,0.1,0), where 500 initial values were randomly chosen within [-5,5] for each value of the parametera. The steady states after transience marked by black dots were shown in Fig.1(a),and the largest Lyapunov exponents (Le1) were indicated by black dots in Fig. 1(b).Since there were many fixed points, many Lyapunov exponents of the map for each parameter value were presented.To see the largest Lyapunov exponents clearly, the smallest Lyapunov exponents were omitted. Some samples of phase portraits of the rational map (1) with different values of the parameterawere presented in Fig. 2. As can be seen from Fig.1,when-4≤a ≤-1,the rational map(1)shows a line of fixed points (Fig. 2(a)). Asadecreases from-1 to-2, a period-2 solution(Fig.2(b))emerges. Whena=-2.667,the period-2 solution bifurcates to a period-8 solution(Fig.2(c)).Whena=-2.892, there is a reverse period-doubling bifurcation, and the period-8 solution bifurcates to a period-16 solution (Fig. 2(d)), and then becomes multiple-piece chaos(Fig. 2(e)) and two-piece chaos (Fig. 2(f)) after a reverse period-doubling cascade. Thereafter the map experiences a small window of periodic solutions (Fig. 2(g)). There is coexistence of stable fixed points,periodic solutions and chaotic solutions(Fig.2(h)). Finally, the two-piece chaotic attractors are merged into one-piece chaotic attractors(Fig.2(i)),which terminate to emerge ata=-3.38.

    Fig. 1. Random bifurcation diagrams of (a) x, and (b) the largest Lyapunov exponent (Le1) of the rational map (1) calculated for a ∈[-4,-1]and (b,c,d)=(0,0.1,0). In Fig. 1(b), the horizontal line denotes the zero value of the largest Lyapunov exponent.

    Fig. 2. Phase portraits of the solutions of the rational map (1) calculated at(b,c,d)=(0,0.1,0) and (a) a=-1.95 (stable fixed points), (b) a=-2.1(stable fixed points and period-2 solution),(c)a=-2.7(stable fixed points and period-8 solution), (d) a=-2.9 (stable fixed points and period-16 solution),(e)a=-2.98(stable fixed points and eight-piece chaotic solution),(f) a = -3.05 (a line of fixed points and two-piece chaotic solution), (g)a = -3.11 (stable fixed points and period-12 solution), (h) a = -3.135(stable fixed points, period-24 solution and four-piece chaotic solution), (i)a=-3.37(stable fixed points and one-piece chaotic solution),respectively.

    To show the stability of a line of fixed points,we plotted the basin of attraction of the rational map(1)whena=-3.05,b=0,c=0.1,andd=0 as shown in Fig.3.The black and red dots represent the chaotic attractor and the stable fixed points,respectively. The basins of the chaotic attractor, the stable fixed points and unbounded solutions were painted cyan,yellow and white, respectively. From Fig. 3, although the line of fixed points are all non-hyperbolic, many fixed points are stable while there exist unstable fixed points. Since the knowledge about the line of fixed points is not very helpful for the localization of attractors,the attractors with a line of fixed points can also be considered as hidden attractors.

    Fig.3. Basins of attraction of the rational map(1)when a=-3.05, b=0,c=0.1,d=0 in the region{(x,y)|x ∈[-30,30],y ∈[-3,3]}.The unbounded basin of attraction which is the set of initial points going into the region({(x,y)||x|+|y|>100}) is shown in white. The chaotic attractor, the stable fixed points and unstable fixed points are denoted by black,red and blue dots,respectively. The basin of chaotic attractors and stable fixed points are shown in cyan and yellow,respectively.

    3.2. Case II:no fixed point A

    By random searching,chaotic attractors were found with several coefficients satisfying the condition thata/=0,b=0,c/=0,andd/=0,i.e.,there exists no fixed point. So the attractors are hidden. For example,whena=-3,b=0,c=1,andd=0.1,the rational map(1)has chaotic attractors.

    To show the complex dynamics of the rational map (1)with no fixed point, forward bifurcation diagram and Lyapunov exponent spectrum diagrams of the map were plotted in Fig.4 by usingaas a varying parameter and fixing(b,c,d)as (0,1,0.1), where the initial value was assigned as (1,0),and the final state at the end of each iteration of the parameter was used as the initial state for the next iteration of computation. The steady states after transience denoted by black dots were presented in Fig.4(a),and the largest Lyapunov exponent (Le1), the smallest Lyapunov exponent (Le2) and the Lyapunov(Kaplan-Yorke)dimension(Dky)were indicated by red,blue and black lines in Fig.4(b),respectively.From Fig.4,the Lyapunov exponent spectrum diagram is in good agreement with the bifurcation diagram. Figure 5 presents some samples of phase portraits of the solutions for the rational map(1). It can be seen from Fig.4 that,whena=-2.361,the rational map(1)shows a hidden period-2 solution(Fig.5(a)).Asadecreases to-2.521, the rational map(1)experiences a reverse period-doubling bifurcation, and the hidden period-2 solution bifurcates to a hidden period-4 solution (Fig. 5(b)).Whena=-2.741, another reverse period-doubling bifurcation is encountered,and this hidden period-4 solution converts into a hidden period-8 solution (Fig. 5(c)). Ata=-2.798,this hidden period-8 solution transforms into a hidden period-16 solution(Fig.5(d)),and then becomes multiple-piece chaos after a reverse period-doubling cascade. After that, a small window of hidden period-6 solutions (Fig. 5(e)) and hidden period-12 solutions are recorded, and the rational map (1)goes into chaotic states again ata=-2.965 through a reverse period-doubling cascade. Finally,whena=-3.144,the hidden two-piece chaotic attractors are jointed together into hidden one-piece chaotic attractors(Fig.5(f)),which disappear ata=-3.148.

    Fig. 4. Forward bifurcation diagram of (a) x, and (b) Lyapunov exponent spectrum(Les)and Lyapunov(Kaplan-Yorke)dimension(Dky)of the rational map(1)calculated for a ∈[-3.2,-2.3]and(b,c,d)=(0,1,0.1)using the initial value(1,0). The largest Lyapunov exponent(Le1), the smallest Lyapunov exponent (Le2) and the Lyapunov (Kaplan-Yorke) dimension (Dky)are indicated by red,blue and black lines,respectively.The dashed horizontal line denotes that the zero value of the Lyapunov exponent and the Lyapunov(Kaplan-Yorke)dimension.

    Fig. 5. Phase portraits of the solutions of the rational map (1) calculated at (b,c,d)=(0,1,0.1) and (a) a=-2.361 (hidden period-2 solution), (b)a=-2.600(hidden period-4 solution),(c)a=-2.780(hidden period-8 solution),(d)a=-2.817(hidden period-16 solution),(e)a=-2.946(hidden period-6 solution), (f) a=-3.144 (one-piece hidden chaotic solution), respectively.

    3.3. Case III:no fixed point B

    Whena=-3,b=1,c=-0.1, andd=-0.1, the condition thatb/=0 and(ad/bc)<1 is fulfilled, i.e., there is no fixed point. So the attractors are hidden. Ifb=1,c=-0.1,d=-0.1, anda <1, then (ad/bc)<1, so there is no fixed point and the attractors are hidden.

    Fig. 6. Backward bifurcation diagram of (a) x, and (b) Lyapunov exponent spectrum(Les)and Lyapunov(Kaplan-Yorke)dimension(Dky)of the rational map(1)calculated for a ∈[-3.14,-2.44]and(b,c,d)=(1,-0.1,-0.1)using the initial value (7,0). The largest Lyapunov exponent (Le1), the smallest Lyapunov exponent (Le2) and Lyapunov (Kaplan-Yorke) dimension(Dky)are denoted by red,blue and black lines,respectively. The dashed horizontal line denotes the zero value of the Lyapunov exponent and the Lyapunov(Kaplan-Yorke)dimension.

    Fig. 7. Phase portraits of the solutions of the rational map (1) calculated at(b,c,d)=(1,-0.1,-0.1)and(a)a=-2.459(hidden period-2 solution),(b)a=-2.508(hidden period-4 solution),(c)a=-2.718(hidden period-8 solution),(d)a=-2.776(hidden period-16 solution),(e)a=-2.926(hidden period-6 solution), (f) a=-3.133 (one-piece hidden chaotic solution),respectively.

    In order to show the complex dynamics of the rational map (1) with no fixed point, backward bifurcation diagram and Lyapunov exponent spectrum diagram of the map were depicted in Fig. 6 by usingaas a bifurcation parameter and fixing(b,c,d)as(1,-0.1,-0.1). The initial value was set as(7,0)and the final state at the end of each iteration of the parameter was used as the initial state for the next iteration of computation. The steady states after transience were shown in Fig.6(a),and the largest Lyapunov exponent(Le1),the smallest Lyapunov exponent (Le2) and the Lyapunov (Kaplan-Yorke)dimension(Dky)were denoted by red,blue and black lines in Fig. 6(b), respectively. From Fig. 6, the Lyapunov exponent diagram agrees well with the bifurcation diagram.Figure 6 presents some samples of phase portraits of the solutions for the rational map(1).It is apparent from Fig.6 that the bifurcation procedure of case III is very similar to that of Case II. Whena=-2.459, a hidden period-2 solution (Fig. 7(a))exists. Asadecreases to-2.504, there is a reverse perioddoubling bifurcation, leading the hidden period-2 solution to a hidden period-4 solution (Fig. 7(b)). Whena=-2.706,the occurrence of another reverse period-doubling bifurcation leads the hidden period-4 solution to a hidden period-8 solution (Fig. 7(c)). Ata=-2.760, this hidden period-8 solution becomes a hidden period-16 solution (Fig. 7(d)),and then changes into hidden multiple-piece chaos after a reverse period-doubling cascade. Hereafter, the map exhibits a small window of hidden period-6 solutions (Fig. 7(e)) and hidden period-12 solutions,and bifurcates into chaos again ata=-2.937 via a reverse period-doubling cascade.Finally,the hidden two-piece chaotic attractors are combined into hidden one-piece chaotic attractors(Fig.7(f)),which cease to exist ata=-3.158.

    3.4. Case IV:a single fixed point

    By random searching, several coefficients satisfying the condition thatb/=0,(ad/bc)=1 and yielding chaotic attractors of the rational map(1)were found. In this case,the rational map(1)has a single non-hyperbolic fixed point.Ifa=-3,b=1,c=0.6, andd=-0.2, the non-hyperbolic fixed point is (1/3,0) and the map shows a two-piece chaotic attractor with the initial value chosen as (0,0). Figures 8(a) and 8(b)show the statesxandyof the chaotic attractor of the rational map (1) as a function of the stepk, respectively. Figure 8(c)presents the phase portrait of the chaotic attractor of the rational map(1). Based on our numerical computation,Lyapunov exponent spectrum (Les) of the chaotic attractor are 0.2182,-0.1304,and its Lyapunov(Kaplan-Yorke)dimension(Dky)is 2,which proves the chaotic property of the rational map(1).

    The basin of attraction of the rational map (1) was depicted in Fig. 9. The cyan and white regions correspond to the basin of the chaotic attractor and the unbounded solution,respectively. The black dots and the blue dot represent the chaotic attractor and the non-hyperbolic fixed point, respectively. From Fig.9,the non-hyperbolic fixed point lies inside the chaotic attractor’s basin,so the initial values chosen from any small punctured neighborhood of the fixed point will tend to the chaotic attractor. However, the basin of the chaotic attractor exhibits to be striped.

    Fig. 8. Time histories of (a) x, (b) y, and (c) phase portrait of the chaotic attractor of the rational map (1) calculated by using the parameter (a,b,c,d)=(-3,1,0.6,-0.2)and the initial value(0,0).

    Fig. 9. Basins of attraction of the rational map (1) when a=-3, b=1,c = 0.6, and d = -0.2 in the region {(x,y)|x ∈[-10,10],y ∈[-10,10]}.The chaotic attractor and the non-hyperbolic fixed point are denoted by black dots and a blue dot, respectively. The basin of the chaotic attractor and the unbounded solution are shown in cyan and white,respectively.

    3.5. Case V:two fixed points

    Ifb/=0 and (ad/bc)>1, the rational map (1) has two fixed points. By numerical exploration,the stabilities of these fixed points can be classified into two cases, i.e., one stable and one unstable fixed points,and two unstable fixed points.In the following,we will discuss these cases in detail.

    3.5.1. Case VA: one stable and one unstable fixed points

    Ifa=-3,b=1,c=0.1,andd=-0.1,thenb/=0 and(ad/bc)>1,i.e.,there are two fixed points. Ifb=1,c=0.1,d=-0.1,anda <-1,then(ad/bc)>1,there are two fixed points.

    In order to show the complex dynamics of the rational map (1) with two fixed points, random bifurcation diagram and Lyapunov exponent spectrum diagram of the map were presented in Fig. 10 by varyingain the range [-4,-1] and fixing (b,c,d) as (1,0.1,-0.1), where 500 initial states were randomly selected in the interval[-5,5]for each value of the parametera. The steady states after transience marked by black dots were shown in Fig. 10(a). The largest Lyapunov exponent (Le1) and the smallest Lyapunov exponent (Le2)were indicated by red and blue dots in Fig.10(b),respectively.According to Fig. 10, there is a good agreement between the Lyapunov exponent spectrum diagram and the bifurcation diagram. When-4<a <-1,there are one unstable fixed point(UFP) and one stable fixed point (SFP), represented by blue dashed lines and red lines, lying in down and up branches,respectively. Asaincreases to-1, the stable and unstable branches intersect and the stable fixed point becomes unstable via a saddle-node bifurcation. Figure 11 presents some samples of phase portraits of the solutions for the rational map (1). For-2.361<a <-1, there only exists one stable fixed point(Fig.11(a)). Whena=-2.361,the map(1)shows a coexisting period-2 solution (Fig. 11(b)). Asadecreases to-2.521, a reverse period-doubling bifurcation occurs, and the coexisting period-2 solution loses stability,leaving a coexisting period-4 solution (Fig. 11(c)). Whena=-2.741, another reverse period-doubling bifurcation arises, and this coexisting period-4 solution bifurcates to a coexisting period-8 solution(Fig.11(d)),and then becomes a coexisting multiplepiece chaotic attractor (Fig. 11(e)) through a reverse perioddoubling cascade. Some small range in which there are more than three coexisting attractors are observed, e.g., the coexistence of fixed point, period-6 solution and chaotic attractor(Fig.11(f)). After that,a cascade of period-doubling bifurcations is observed and the state of the map(1)goes to the coexisting period-8 solution(Fig.11(g)),which evolves into chaos again as a result of a reverse period-doubling cascade. Finally,whena=-3.14,the two-piece chaotic attractors(Fig.11(h))are formed into one-piece chaotic attractors(Fig.11(i)).When-4<a <-3.38,there only exists one stable fixed point.

    Fig. 10. Random bifurcation diagram of (a) x, and (b) Lyapunov exponent spectrum(Les)of the rational map(1)calculated for a ∈[-4,-1]and(b,c,d)=(1,0.1,-0.1). The stable fixed points(SFP),unstable fixed points(UFP) and attractors expect for the stable fixed points are denoted by red lines,dashed lines and black dots,respectively. The largest Lyapunov exponent (Le1) and the smallest Lyapunov exponent (Le2) are indicated by red and blue dots,respectively. In Fig.10(b),the horizontal line denotes the zero value of the Lyapunov exponent.

    Fig.12. Basins of attraction of the rational map(1)when(a)a=-2.98,(b)a=-3.06 and b=1,c=0.1,d=-0.1 in the region{(x,y)|x ∈[-25,25],y ∈[-2,2]}. The chaotic attractor, the period-6 attractor, the stable fixed point and the unstable fixed point are represented by black,magenta,red and blue dots,respectively. The basins of the chaotic attractors,the period-6 attractor,the stable fixed point and the unbounded solution are shown in cyan,green,yellow and white,respectively.

    Fig.11. Phase portraits of the solutions of the rational map(1)calculated at(b,c,d)=(1,0.1,-0.1)and(a)a=-1(a fixed point),(b)a=-2.24(fixed point and period-2 solution), (c) a=-2.56 (fixed point and period-4 solution),(d)a=-2.79(fixed point and period-8 solution),(e)a=-2.94(fixed point and chaotic solution),(f)a=-2.98(one fixed point,period-6 solution and chaotic solution),(g)a=-3.04(fixed point and period-8 solution),(h)a=-3.06(fixed point and two-piece chaotic solution),(i)a=-3.14(fixed point and one-piece chaotic solution),respectively.

    To show the hidden and self-excited attractors of the rational map(1),we computed the basin of attraction of the map whena=-2.98 anda=-3.06,b=1,c=0.1,d=-0.1,as demonstrated in Fig.12,respectively. The chaotic attractor,the period-6 attractor, the stable fixed point and the unstable fixed point were represented by black, magenta, red and blue dots, respectively. The basins of the chaotic attractors,the period-6 attractor, the stable fixed points and the unbounded solutions were colored in cyan, green, yellow and white, respectively. From Fig.12(a), the period-6 attractor is self-excited since one unstable fixed point stays in the basin of the period-6 attractor. However,the chaotic attractor is hidden since the basin of the chaotic attractor is not connected with any small neighborhoods of the fixed points. It can be concluded from Fig.12(b)that,the chaotic attractor is self-excited since one unstable fixed point lies in the basin of the chaotic attractor, i.e., the basin of the chaotic attractor is connected with the small neighborhoods of the unstable fixed point.

    3.5.2. Case VB: two unstable fixed points

    Ifa=-1,b=-1,c=1,andd=2,then we haveb/=0 and(ad/bc)>1,i.e.,there exist two fixed points. Ifa=-1,c=1,d=2 and-2<b <0, then (ad/bc)>1. Thus the rational map(1)admits two fixed points. In the following,we chosea=-1,c=1,d=2, and selectedbas a branching parameter. Figure 13 presents the fixed points of the rational map (1) calculated forb ∈[-2.1,-0.01] and (a,c,d)=(-1,1,2). The stable and unstable fixed points were indicated by red and blue lines,respectively.For-2.1<b <-2,there is an unstable fixed point. Then the rational map(1)experiences a saddle-node bifurcation atb=-2, which yields one stable fixed point and one unstable fixed point. When the parameterbincreases to-1.867, the stable fixed point loses stability and becomes unstable through a saddle-node bifurcation.When the parameterbincreases further, this unstable fixed point gains stability and becomes stable atb=-0.133 via a saddle-node bifurcation. Thus there are two unstable fixed points forb ∈[-1.866,-0.132].

    Fig. 13. The fixed points of the rational map (1) calculated for b ∈[-2.1,-0.001]and(a,c,d)=(-1,1,2).The stable and unstable fixed points are indicated by red and blue lines,respectively.

    Fig. 14. Random bifurcation diagram of (a) x, and (b) Lyapunov exponent spectrum(Les)of the rational map(1)calculated for b ∈[-1.1,-0.001]and(a,c,d)=(-1,1,2) using the initial value randomly chosen in the interval[-5,5]. The largest Lyapunov exponent (Le1) and the smallest Lyapunov exponent(Le2)are indicated by red and blue dots,respectively. The dashed horizontal line denotes the zero value of the Lyapunov exponent.

    Fig.15. Phase portraits of the solutions of the rational map(1)calculated at(a,c,d)=(-1,1,2)and(a)b=-0.58(period-11 solution), (b)b=-0.55(period-12 solution), (c) b = -0.46 (period-13 solution), (d) b = -0.40(period-14 solution), (e) b = -0.37 (period-15 solution), (f) b = -0.34(period-16 solution), (g) b = -0.32 (period-17 solution), (h) b = -0.31(period-18 solution), (i) b = -0.27 (period-19 solution), (j) b = -0.26(period-20 solution), (k) b = -0.24 (period-21 solution), (l) b = -0.23(period-22 solution),respectively.

    For showing the complex dynamics of the rational map (1) with two fixed points, random bifurcation diagram and Lyapunov exponent spectrum diagram of the map were illustrated in Fig. 14 by takingbas a control parameter and fixing (a,c,d) as (-1,1,2), where 500 initial states were randomly selected in the interval [-5,5] for each value of the parameterb. The steady states after transience were depicted in Fig. 14(a). The largest Lyapunov exponent (Le1)and the smallest Lyapunov exponent (Le2) were indicated by red and blue dots in Fig.14(b),respectively. From Fig.14,the Lyapunov exponent diagram coincides with the bifurcation diagram. Figure 15 presents some samples of phase portraits of the solutions for the rational map (1). It is apparent from Fig. 14 that, there exist many chaotic regimes and periodic windows. For-0.58≤b <0, the period of the periodic solutions increases by one from eleven, as displayed in Figs. 15(a)-15(l), which can be considered as periodadding bifurcation. And the length of the chaotic regimes becomes shorter as the parameteraincreases. For-1.2<b <-0.58, several period-doubling bifurcations and inverse period-doubling bifurcations are observed.Finally,the chaotic attractors disappear atb=-1.075. Although there exist two unstable fixed points, the rational map (1) may also exhibit hidden attractors.

    Fig. 16. Basins of attraction of the rational map (1) when (a) b=-0.745,(b)b=-0.58 and a=-1,c=1,d=2,respectively. The chaotic attractors,the period-11 attractor,the stable fixed point and the unstable fixed point are denoted by black, magenta, red and blue dots, respectively. The basins of the chaotic attractors, the period-11 attractor, the stable fixed point and the unbounded solution are shown in cyan,green,yellow and white,respectively.

    To distinguish the hidden and self-excited attractors of the rational map (1), we drew the basin of attraction of the map(1)forb=-0.745 andb=-0.58,a=-1,c=1,d=2,as demonstrated in Fig.16. The chaotic attractor, the period-6 attractor, the stable fixed point and the unstable fixed point were represented by black, magenta, red and blue dots, respectively. The basins of the chaotic attractors, the period-11 attractor, the stable fixed points and the unbounded solutions were colored in cyan, green, yellow and white, respectively.From Fig.16(a), the basins of the period-11 attractor and the chaotic attractor are fractal,so the basins of the chaotic attractor and period-11 attractor are connected with the small neighborhoods of the fixed points.Hence the period-11 attractor and the chaotic attractor are all self-excited. It can be found from Fig.16(b)that,the period-11 attractor is self-excited since one unstable fixed point lies in the basin of the period-11 attractor,i.e.,the basins of the period-11 attractor is connected with the small neighborhoods of the unstable fixed point. We also verified the class of attractors according to the definition of hidden and self-excited attractors by exploiting the numerical tools as follows. Firstly, we obtained two fixed points. Then we randomly chose initial points in a very small neighborhood(<0.001) of these fixed points. Finally, some points could tend to these self-excited attractors.

    4. Conclusions

    A new class of two-dimensional rational maps with different types of fixed points was introduced and studied in this paper. The existence and stabilities of fixed points of the rational map were discussed. Several numerical analysis tools were employed to demonstrate the rich and complex dynamics of the rational map. Both self-excited and hidden attractors were shown and explored in the rational map. In addition,multi-stability, especially the hidden multi-stability, was further investigated. The proposed rational map can be applied to secure network communications, such as data and image encryption.[32]Future works will focus on the investigation of high-dimensional rational maps with self-excited and hidden attractors.

    Acknowledgments

    Project supported by the National Natural Science Foundation of China (Grant Nos. 11672257, 11772306,11972173, and 12172340) and the 5th 333 High-level Personnel Training Project of Jiangsu Province of China (Grant No.BRA2018324).

    猜你喜歡
    張麗萍劉洋海波
    搏浪
    汽車(chē)ABS控制仿真分析
    又見(jiàn)劉洋
    海峽姐妹(2020年6期)2020-07-25 01:26:10
    Three dimensional nonlinear shock waves in inhomogeneous plasmas with different size dust grains and external magnetized field
    山清水秀
    劉洋作品
    藝術(shù)家(2019年9期)2019-12-17 08:28:19
    說(shuō)海波
    劉洋 藏石欣賞
    寶藏(2018年8期)2018-08-31 07:28:00
    劉洋作品
    張麗萍 勿忘初心 立己達(dá)人
    最近最新免费中文字幕在线| 一进一出抽搐gif免费好疼| 韩国av一区二区三区四区| 直男gayav资源| 成人无遮挡网站| 狠狠狠狠99中文字幕| 有码 亚洲区| 一区福利在线观看| 两个人的视频大全免费| a级一级毛片免费在线观看| 蜜桃久久精品国产亚洲av| 久久久国产成人精品二区| 午夜爱爱视频在线播放| 一边摸一边抽搐一进一小说| 一个人免费在线观看电影| 国产高清有码在线观看视频| 狠狠狠狠99中文字幕| 亚洲欧美日韩无卡精品| 午夜福利在线观看免费完整高清在 | 日韩欧美精品免费久久| 免费看美女性在线毛片视频| 亚洲精品色激情综合| 国产免费av片在线观看野外av| 无遮挡黄片免费观看| 两个人的视频大全免费| 久久香蕉精品热| 亚洲va日本ⅴa欧美va伊人久久| 成人一区二区视频在线观看| 干丝袜人妻中文字幕| 无人区码免费观看不卡| 麻豆成人午夜福利视频| 天天躁日日操中文字幕| 校园人妻丝袜中文字幕| 99九九线精品视频在线观看视频| 国产精品人妻久久久影院| 男女下面进入的视频免费午夜| 欧美日韩精品成人综合77777| 国产真实伦视频高清在线观看 | 免费人成视频x8x8入口观看| 成人午夜高清在线视频| 91久久精品国产一区二区成人| 男女视频在线观看网站免费| 久久6这里有精品| 天美传媒精品一区二区| 露出奶头的视频| 国产精品自产拍在线观看55亚洲| 成人特级av手机在线观看| 国产一区二区激情短视频| 久久精品影院6| 国产精品,欧美在线| 国产蜜桃级精品一区二区三区| 亚洲真实伦在线观看| 国产乱人视频| 身体一侧抽搐| 国产一区二区三区视频了| av在线蜜桃| 老熟妇仑乱视频hdxx| 亚洲第一区二区三区不卡| 国产视频一区二区在线看| 黄色日韩在线| 女同久久另类99精品国产91| 一个人看的www免费观看视频| 在线天堂最新版资源| av在线蜜桃| 午夜福利在线在线| 麻豆久久精品国产亚洲av| 悠悠久久av| 免费在线观看影片大全网站| 99视频精品全部免费 在线| 99久久精品一区二区三区| 精品国内亚洲2022精品成人| 色综合婷婷激情| 天堂av国产一区二区熟女人妻| 日韩大尺度精品在线看网址| 麻豆一二三区av精品| 亚洲国产日韩欧美精品在线观看| 老司机深夜福利视频在线观看| 国产精品乱码一区二三区的特点| 国产精品乱码一区二三区的特点| 亚洲va日本ⅴa欧美va伊人久久| 亚洲va日本ⅴa欧美va伊人久久| 免费无遮挡裸体视频| 免费av不卡在线播放| 亚洲av五月六月丁香网| 麻豆成人av在线观看| 白带黄色成豆腐渣| 在线观看免费视频日本深夜| 小说图片视频综合网站| 色综合色国产| 免费av不卡在线播放| 亚洲电影在线观看av| 1000部很黄的大片| 国产午夜福利久久久久久| 国产高潮美女av| 人人妻,人人澡人人爽秒播| 在线观看av片永久免费下载| 日本免费一区二区三区高清不卡| 亚洲真实伦在线观看| av天堂在线播放| 免费看a级黄色片| avwww免费| 国产91精品成人一区二区三区| 国产男人的电影天堂91| 日韩强制内射视频| 成年女人永久免费观看视频| 亚洲av一区综合| 嫩草影院新地址| 丝袜美腿在线中文| 亚洲精品影视一区二区三区av| 成人无遮挡网站| 简卡轻食公司| 久久精品国产亚洲网站| 女的被弄到高潮叫床怎么办 | av福利片在线观看| av福利片在线观看| 国产精品综合久久久久久久免费| 22中文网久久字幕| 12—13女人毛片做爰片一| 老熟妇仑乱视频hdxx| 午夜免费成人在线视频| 在线观看一区二区三区| 国产精品国产三级国产av玫瑰| 国产av不卡久久| 乱人视频在线观看| 91久久精品国产一区二区成人| 俺也久久电影网| 国产精品美女特级片免费视频播放器| 成年女人毛片免费观看观看9| 国产av一区在线观看免费| 精品一区二区三区av网在线观看| 窝窝影院91人妻| 欧美精品啪啪一区二区三区| 欧美最新免费一区二区三区| 亚洲中文字幕日韩| 人妻夜夜爽99麻豆av| 亚洲一区二区三区色噜噜| 午夜福利在线在线| 亚洲av.av天堂| 99九九线精品视频在线观看视频| 亚洲精品久久国产高清桃花| 可以在线观看的亚洲视频| 蜜桃久久精品国产亚洲av| 国产欧美日韩精品一区二区| 免费人成在线观看视频色| 日韩欧美在线乱码| 波野结衣二区三区在线| 国产高潮美女av| 一级a爱片免费观看的视频| www.色视频.com| 亚洲熟妇中文字幕五十中出| 舔av片在线| 性插视频无遮挡在线免费观看| 日本欧美国产在线视频| 成年免费大片在线观看| 日本在线视频免费播放| 午夜精品一区二区三区免费看| 亚洲性久久影院| 欧美最新免费一区二区三区| 国产精品一区二区性色av| 精品人妻偷拍中文字幕| 波多野结衣高清无吗| 久久久久久伊人网av| 日韩强制内射视频| 97热精品久久久久久| 欧美xxxx性猛交bbbb| 人妻久久中文字幕网| 国内精品久久久久精免费| 在线播放无遮挡| 亚洲电影在线观看av| 99国产极品粉嫩在线观看| 亚洲综合色惰| 最近视频中文字幕2019在线8| 精品久久久久久久末码| 最近在线观看免费完整版| 成人国产一区最新在线观看| 少妇丰满av| 亚洲熟妇熟女久久| 我要搜黄色片| 俺也久久电影网| 一级av片app| 日日夜夜操网爽| 欧美激情久久久久久爽电影| ponron亚洲| 欧美精品国产亚洲| 搡老妇女老女人老熟妇| 亚洲四区av| 亚洲成人久久爱视频| 国产又黄又爽又无遮挡在线| 熟妇人妻久久中文字幕3abv| 亚洲色图av天堂| 亚洲精品亚洲一区二区| 88av欧美| а√天堂www在线а√下载| 亚洲av免费高清在线观看| 性欧美人与动物交配| 免费人成在线观看视频色| 熟女电影av网| 亚洲成人久久性| 国产主播在线观看一区二区| 麻豆国产97在线/欧美| 91在线精品国自产拍蜜月| 麻豆成人午夜福利视频| 亚洲avbb在线观看| 18+在线观看网站| 国产亚洲精品久久久久久毛片| 亚洲自拍偷在线| 麻豆国产97在线/欧美| 亚洲国产精品成人综合色| 简卡轻食公司| 中文字幕人妻熟人妻熟丝袜美| 成人高潮视频无遮挡免费网站| 日韩一本色道免费dvd| 香蕉av资源在线| 美女xxoo啪啪120秒动态图| 毛片女人毛片| 婷婷丁香在线五月| 亚洲精华国产精华液的使用体验 | 成年人黄色毛片网站| 麻豆av噜噜一区二区三区| 亚洲精华国产精华液的使用体验 | 99精品久久久久人妻精品| 如何舔出高潮| 亚洲三级黄色毛片| 熟女电影av网| 久久99热这里只有精品18| 成人午夜高清在线视频| 欧美3d第一页| 国产探花在线观看一区二区| 三级男女做爰猛烈吃奶摸视频| 亚洲经典国产精华液单| 老师上课跳d突然被开到最大视频| 乱系列少妇在线播放| 国产私拍福利视频在线观看| 国产视频内射| 亚洲自拍偷在线| 亚洲自偷自拍三级| 淫秽高清视频在线观看| 国产欧美日韩一区二区精品| 黄色丝袜av网址大全| 变态另类成人亚洲欧美熟女| 日韩精品有码人妻一区| 狠狠狠狠99中文字幕| 欧美中文日本在线观看视频| 亚洲无线观看免费| av在线观看视频网站免费| 国产精品爽爽va在线观看网站| 最近在线观看免费完整版| 午夜视频国产福利| 亚洲美女黄片视频| 欧美zozozo另类| 久久精品国产鲁丝片午夜精品 | 一个人观看的视频www高清免费观看| 国内精品久久久久久久电影| 国产探花在线观看一区二区| 不卡一级毛片| 黄色女人牲交| 嫩草影院入口| 亚洲不卡免费看| 国产精品98久久久久久宅男小说| 少妇人妻一区二区三区视频| 嫩草影院精品99| 在线看三级毛片| 少妇的逼水好多| 麻豆久久精品国产亚洲av| 亚洲av.av天堂| 18禁裸乳无遮挡免费网站照片| 亚洲欧美精品综合久久99| 国产伦一二天堂av在线观看| 嫩草影视91久久| 国产一区二区在线av高清观看| 在现免费观看毛片| 十八禁国产超污无遮挡网站| 国产高清不卡午夜福利| 日韩精品有码人妻一区| 搡老妇女老女人老熟妇| 美女 人体艺术 gogo| 国产激情偷乱视频一区二区| 国产视频一区二区在线看| 日本色播在线视频| 精品久久国产蜜桃| 国内精品宾馆在线| 观看免费一级毛片| 国产一区二区三区视频了| 日本精品一区二区三区蜜桃| 国产淫片久久久久久久久| 国产成人a区在线观看| 熟女电影av网| 深夜精品福利| 亚洲最大成人手机在线| 一进一出好大好爽视频| 精品乱码久久久久久99久播| 免费人成视频x8x8入口观看| 午夜福利成人在线免费观看| 国内精品一区二区在线观看| 3wmmmm亚洲av在线观看| 亚洲狠狠婷婷综合久久图片| 午夜精品久久久久久毛片777| 人妻丰满熟妇av一区二区三区| 永久网站在线| 波多野结衣高清作品| 久久午夜亚洲精品久久| 九色国产91popny在线| 国产精品久久久久久久电影| 搡女人真爽免费视频火全软件 | 亚洲中文字幕一区二区三区有码在线看| 日韩欧美三级三区| 真人一进一出gif抽搐免费| 男女那种视频在线观看| 我的老师免费观看完整版| 十八禁网站免费在线| 在现免费观看毛片| 国内精品久久久久精免费| 国产一区二区三区视频了| 久久久久久久久久黄片| 俺也久久电影网| 啦啦啦观看免费观看视频高清| 两个人的视频大全免费| 长腿黑丝高跟| 久9热在线精品视频| 日韩大尺度精品在线看网址| 波多野结衣高清无吗| 乱人视频在线观看| 欧美成人免费av一区二区三区| 看免费成人av毛片| 午夜福利欧美成人| 午夜福利在线观看吧| 99热这里只有是精品50| 成人av一区二区三区在线看| 亚洲av免费在线观看| 美女cb高潮喷水在线观看| 精品一区二区三区av网在线观看| 精品一区二区三区视频在线观看免费| 亚洲av美国av| 麻豆成人午夜福利视频| 99久久精品一区二区三区| 免费看a级黄色片| 日韩欧美免费精品| 老司机午夜福利在线观看视频| 网址你懂的国产日韩在线| 国产免费av片在线观看野外av| 亚洲精品久久国产高清桃花| 国产综合懂色| 我的老师免费观看完整版| 九色成人免费人妻av| 夜夜夜夜夜久久久久| 日韩欧美在线乱码| 日韩欧美国产一区二区入口| 九九热线精品视视频播放| 一级黄片播放器| 中国美白少妇内射xxxbb| 人人妻,人人澡人人爽秒播| 制服丝袜大香蕉在线| 久久精品国产亚洲网站| 国产探花极品一区二区| 五月伊人婷婷丁香| 啪啪无遮挡十八禁网站| 色尼玛亚洲综合影院| 夜夜看夜夜爽夜夜摸| 少妇的逼好多水| 中文字幕精品亚洲无线码一区| 看黄色毛片网站| 午夜福利在线在线| 成人国产麻豆网| 成人无遮挡网站| 99久国产av精品| 国产午夜精品论理片| 国产精品福利在线免费观看| 国产黄色小视频在线观看| 亚洲avbb在线观看| 国产视频一区二区在线看| 亚洲内射少妇av| 少妇人妻精品综合一区二区 | 久久欧美精品欧美久久欧美| 午夜亚洲福利在线播放| 亚洲图色成人| 国产免费一级a男人的天堂| 国产精品自产拍在线观看55亚洲| 女的被弄到高潮叫床怎么办 | 夜夜爽天天搞| 99热这里只有是精品50| 两个人的视频大全免费| 变态另类丝袜制服| 干丝袜人妻中文字幕| 黄色配什么色好看| 免费不卡的大黄色大毛片视频在线观看 | 欧美丝袜亚洲另类 | 嫩草影院精品99| 一卡2卡三卡四卡精品乱码亚洲| 在线观看免费视频日本深夜| 香蕉av资源在线| 无人区码免费观看不卡| 又紧又爽又黄一区二区| a级毛片a级免费在线| 尤物成人国产欧美一区二区三区| 日本免费一区二区三区高清不卡| 制服丝袜大香蕉在线| 别揉我奶头 嗯啊视频| 99久久久亚洲精品蜜臀av| 天美传媒精品一区二区| 精品人妻一区二区三区麻豆 | 国产在视频线在精品| 色哟哟哟哟哟哟| 深爱激情五月婷婷| 国产av一区在线观看免费| 午夜福利高清视频| 在线观看美女被高潮喷水网站| 亚洲国产精品久久男人天堂| 国产免费一级a男人的天堂| 我要搜黄色片| 淫妇啪啪啪对白视频| av在线观看视频网站免费| 在线观看午夜福利视频| 小蜜桃在线观看免费完整版高清| 最新中文字幕久久久久| 别揉我奶头~嗯~啊~动态视频| 精品日产1卡2卡| 欧美日韩瑟瑟在线播放| 日韩欧美三级三区| 一级av片app| 91久久精品电影网| 狂野欧美激情性xxxx在线观看| 久久草成人影院| 国产精品综合久久久久久久免费| 日韩av在线大香蕉| 99热精品在线国产| 少妇丰满av| 久9热在线精品视频| 国产白丝娇喘喷水9色精品| 国产成人福利小说| 亚洲图色成人| 神马国产精品三级电影在线观看| 国产精品久久视频播放| 可以在线观看毛片的网站| 一区福利在线观看| 久久久国产成人免费| 亚洲国产精品合色在线| 国产一区二区在线av高清观看| 成熟少妇高潮喷水视频| 亚洲 国产 在线| 一夜夜www| 啦啦啦观看免费观看视频高清| 美女免费视频网站| 国产成人aa在线观看| 小说图片视频综合网站| 国产精品女同一区二区软件 | 久久这里只有精品中国| 麻豆成人午夜福利视频| 网址你懂的国产日韩在线| 国产一区二区在线av高清观看| 蜜桃久久精品国产亚洲av| 日本与韩国留学比较| 最近在线观看免费完整版| 欧美性感艳星| 国产不卡一卡二| 亚洲国产色片| 亚洲经典国产精华液单| 欧美性猛交╳xxx乱大交人| av在线老鸭窝| 精华霜和精华液先用哪个| 亚洲成人久久爱视频| 麻豆成人av在线观看| 一夜夜www| 久久久久久久久中文| 搡女人真爽免费视频火全软件 | 亚洲国产欧美人成| 精品午夜福利视频在线观看一区| 婷婷精品国产亚洲av在线| 一边摸一边抽搐一进一小说| 亚洲最大成人中文| 色综合站精品国产| 老师上课跳d突然被开到最大视频| 久久人人精品亚洲av| 嫁个100分男人电影在线观看| 欧美日韩瑟瑟在线播放| 91久久精品国产一区二区三区| 乱码一卡2卡4卡精品| 91午夜精品亚洲一区二区三区 | 国产精品一区二区三区四区免费观看 | 日韩,欧美,国产一区二区三区 | 天美传媒精品一区二区| 成人综合一区亚洲| 亚洲av.av天堂| 中文字幕免费在线视频6| 日本与韩国留学比较| 能在线免费观看的黄片| 熟女电影av网| 69av精品久久久久久| 欧美日本视频| 国产成人福利小说| 国产色婷婷99| 美女黄网站色视频| 女人被狂操c到高潮| 国产高潮美女av| 国产精华一区二区三区| 日本成人三级电影网站| 搡老妇女老女人老熟妇| 亚洲精华国产精华液的使用体验 | av专区在线播放| 亚洲av.av天堂| 老司机福利观看| 自拍偷自拍亚洲精品老妇| 国内精品宾馆在线| 国产熟女欧美一区二区| 毛片一级片免费看久久久久 | 91午夜精品亚洲一区二区三区 | 国产一级毛片七仙女欲春2| 两个人视频免费观看高清| 亚洲精华国产精华液的使用体验 | 亚洲成av人片在线播放无| 欧美日韩中文字幕国产精品一区二区三区| 观看免费一级毛片| 一进一出好大好爽视频| 国产单亲对白刺激| 亚洲欧美日韩无卡精品| 在线看三级毛片| 日日啪夜夜撸| 少妇人妻精品综合一区二区 | 一个人观看的视频www高清免费观看| 欧美激情国产日韩精品一区| or卡值多少钱| 麻豆一二三区av精品| 欧美精品国产亚洲| 日本一本二区三区精品| 久久久成人免费电影| 亚洲美女搞黄在线观看 | 国产精品av视频在线免费观看| 91麻豆av在线| 日韩欧美精品v在线| x7x7x7水蜜桃| 22中文网久久字幕| 成人性生交大片免费视频hd| 九九在线视频观看精品| 成人无遮挡网站| 在线观看66精品国产| 欧美成人a在线观看| 搞女人的毛片| 国产伦人伦偷精品视频| 亚洲午夜理论影院| 狂野欧美激情性xxxx在线观看| 人妻丰满熟妇av一区二区三区| 五月玫瑰六月丁香| 国产高清三级在线| 久久久久久久久大av| 国内久久婷婷六月综合欲色啪| 在线观看一区二区三区| 干丝袜人妻中文字幕| 国产午夜精品久久久久久一区二区三区 | 免费在线观看影片大全网站| 97热精品久久久久久| 熟女电影av网| 欧美中文日本在线观看视频| 女同久久另类99精品国产91| 国产精品久久视频播放| 亚洲精品一卡2卡三卡4卡5卡| 亚洲国产精品成人综合色| 欧美不卡视频在线免费观看| 亚洲aⅴ乱码一区二区在线播放| 天天一区二区日本电影三级| 日本免费a在线| 尾随美女入室| 超碰av人人做人人爽久久| 欧美激情国产日韩精品一区| 在线观看美女被高潮喷水网站| 国产成人一区二区在线| 中文字幕熟女人妻在线| 午夜日韩欧美国产| 99九九线精品视频在线观看视频| 91狼人影院| 日本免费a在线| 欧美激情国产日韩精品一区| 91在线精品国自产拍蜜月| 欧美不卡视频在线免费观看| 欧美黑人巨大hd| 老熟妇乱子伦视频在线观看| 中亚洲国语对白在线视频| 精品不卡国产一区二区三区| 18禁裸乳无遮挡免费网站照片| 成年女人看的毛片在线观看| 久9热在线精品视频| 久久久久久九九精品二区国产| 国产精品美女特级片免费视频播放器| 欧美一区二区精品小视频在线| 国产中年淑女户外野战色| 亚洲乱码一区二区免费版| 男人狂女人下面高潮的视频| 欧美日韩精品成人综合77777| 一卡2卡三卡四卡精品乱码亚洲| 亚洲美女视频黄频| 老女人水多毛片| 亚洲精品色激情综合| 欧美潮喷喷水| 国产精品爽爽va在线观看网站| 日日摸夜夜添夜夜添av毛片 | 在线国产一区二区在线| 日本免费a在线| 亚洲av五月六月丁香网| a在线观看视频网站| 不卡一级毛片| 久99久视频精品免费| 国产男靠女视频免费网站| 精品午夜福利视频在线观看一区| 婷婷六月久久综合丁香| 精品久久久噜噜| 国产成人福利小说| 中文资源天堂在线| 69av精品久久久久久| 国产精品亚洲美女久久久| 69av精品久久久久久| 熟女电影av网| 国产精品一区www在线观看 | 1024手机看黄色片| 久久久午夜欧美精品| 极品教师在线视频| 久久精品夜夜夜夜夜久久蜜豆| 国产色爽女视频免费观看|