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

    Heterogeneous dual memristive circuit: Multistability,symmetry,and FPGA implementation?

    2021-12-22 06:46:10YiZiCheng承亦梓FuHongMin閔富紅ZhiRui芮智andLeiZhang張雷
    Chinese Physics B 2021年12期
    關(guān)鍵詞:張雷

    Yi-Zi Cheng(承亦梓), Fu-Hong Min(閔富紅), Zhi Rui(芮智), and Lei Zhang(張雷)

    School of Electrical and Automation Engineering,Nanjing Normal University,Nanjing 210023,China

    Keywords: memristive circuit,chaos,multistability,FPGA implementation

    1. Introduction

    Memristor,the 4th basic circuit element,[1]which reveals the relationship between flux and charge, has great application potentials in the fields of neural network,[2]nonlinear circuit design,[3]image encryption,[4]etc. In recent years, various novel memristor-based circuits were widely built, and the complex dynamical behaviors of such memristive systems have been also studied. For example, a Duffing oscillator with memristor was experimentally realized in Ref. [5],and a Chua’s circuit with memristor, which shows extreme multistability and other complex dynamics, was proposed in Ref. [6]. In Ref. [7], a discrete memristor model was defined and applied to H′enon map. An Sr0.95Ba0.05TiO3(SBT)memristor is a new kind of physical memristor, and the complex behaviors of memristive circuit based on SBT was studied in Ref. [8]. A variety of autonomous or non-autonomous oscillators based on memristors expanded the research areas for these nonlinear systems.[9–11]Compared with the conventional chaotic circuits,[12]the memristor-based circuits can easily procduce diverse nonlinear motions,such as self-excited chaotic behaviors,[13]hidden attractors,[14]hyperchaos,[15]and anti-monotonic properties.[16]Because the memristive circuit is sensitive to initial values, the multistability, a typical behavior of such circuits, has become a research hotspot in recent years. The multistability in circuits with single or multiple memristors has been reported by a large number of researches.[17–21]For example, the serial-parallel memristor was introduced into the Chua’s circuit, and the influence of memristor polarity and series-parallel structure on mutistability phenomenon was explored and concluded.[2]The memristive Chua’s circuit based on active band-pass filter (BPF)was studied,in which the stable distribution of the equilibrium points was divided to confirm the trace of memristive circuit starting from different initial states.[21]And then,the extreme multistability of infinite coexistence attractors was observed,which provides the guidance for further exploring the relationship between the equilibrium points and the multistability.

    To make in depth analysis of the multistability, a new modeling method called state variable mapping (SVM) was proposed.[22]However, the SVM has its disadvantages: it is suitable only for the normalized state equation and those memristive systems with only simple nonlinear terms. To improve SVM, Chenet al.proposed a method of combining the incremental integral change with linear state variable,[23]which was named hybrid state variable incremental integral(HSVII).As the dimensionality reduction methods mentioned above still cannot directly simplify the model, in the light of flux–charge analysis method(FCAM)[24]which offers an effective way to analyze the infinite many dynamics caused by memristor, a dimensionality reduction analysis was carried out in the memeristive circuits in Ref.[25]. However,because of the complexity of the equation derivation and state variable selection, the application of FCAM in memristive systems has rarely been reported.

    In this paper, a heterogeneous DMC with two different memristors, based on Chua’s circuit, is constructed, and its nonlinear dynamical behaviors are investigated from multiple perspectives. This DMC can be used for encrypting an image.The rest of this paper is organized as follows. In Section 2,the DMC schematics are presented, and their mathematical models in the voltage–current domains and the flux–charge domains are established respectively. In Section 3, the symmetrical bifurcation and the parameter mappings with specific parameters are discussed by the 5th-order model(5OM)in the volt–ampere domain. Subsequently in Section 4, the dimensionality reduction model, 3rd-order model (3OM) in flux–charge domain,is derived by FCAM,and its hardware circuit is realized through field-programmable gate array (FPGA).Then,the multistability is illustrated through attraction basins and phase diagrams in 3OM and 5OM.Finally,some conclusions are drawn in Section 6.

    2. Circuit modeling

    In this section, a memristive circuit with two different memristors is proposed based on Chua’s circuit. The voltage–ampere model of heterogeneous dual memristive circuit is derived by Kirchhoff’s law and the volt–ampere relationship,and the flux–charge model is subsequently obtained through flux–charge analysis method to make it more controllable.

    2.1. Voltage–ampere model

    The heterogeneous DMC is constructed based on memristive Chua’s circuit[26]as shown in Fig. 1, where the circuit consists of two capacitorsC1andC2, an inductorL, a linear negative conductance-G, a third-order nonlinear flux–controlled memristorW, and a charge-controlled memristorMdescribed by absolute value. For the rationality of the mathematical model,the flux–controlled memristorW(?)and the capacitorC1are connected in parallel, while the chargecontrolled memristorM(q)and the inductorLare in series.

    Fig.1. Heterogeneous DMC.

    According to Kirchhoff’s law and the volt–ampere relationship, the state equation of the DMC can be derived into five first-order differential equations:

    where the flux-controlled memristor is described asW(?)=a+3b?2,the charge-controlled memristorM(q)=?c+d|q|,v1andv2represent the voltage of capacitorC1andC2,ithe current through the parallel branch of the charge-controlled memristor and the inductanceL,?the internal flux of the cubic memristor,qthe internal charge of the quadratic memristor. Assumingv1=x,v2=y,i=z,?=w,q=v,1/C1=α,1/C2=β,1/L=η,G=γ,we can obtain 5OM from Eq.(1)as follows:

    whereW(w)=a+3bw2andM(v)=?c+d|v|. For convenience,the system parameters are taken as

    2.2. Flux–charge model

    Through the flux–charge method, the increment of flux and charge from 0 totare defined as?(t;0)=?(t)??(0)andq(t;0) =q(t)?q(0), respectively, where?(0,0) = 0 andq(0,0)=0, the flux increments are?M1(t;0),?M2(t;0),?C1(t;0),?C2(t;0),?L(t;0), and?G(t;0). The charge increments areqM1(t;0),qM2(t;0),qC1(t;0),qC2(t;0),qL(t;0),andqG(t;0),and the reference directions are shown in Fig.2.

    Fig.2. DMC with reference directions.

    where the constitutive relationships of the elements can be obtained as follows:

    The initial values of the above elements are expressed asqC1(0) =C1uC1(0),qC2(0) =C2uC2(0),?L(0) =LiL(0),anduC1(0),uC2(0),iL(0) are equivalent to the initial values of the first three equations in Eq. (1). The model of fluxcontrolled memristor can thus be expressed asqM1(t;0) =f[?M1(t;0)+?M1(0)]?qM1(0), and the charge-controlled memristor?M2(t;0) =h[qM2(t;0)+qM2(0)]??M2(0), in whichqM1(0)=f[?M1(0)],?M2(0)=h[qM2(0)]. According to the actual circuit, the mathematical model of two memristors can be described as

    Here,we seta=?0.3,b=0.6,c=0.03,andd=0.02.

    According to Kirchhoff’s law and the circuit model in Fig.2,it follows that

    and the third-order differential equation

    Substituting Eq.(9)into Eq.(10),we obtain in the flux–charge domain

    Compared with the 5OM, the 3OM converts the initial values into the system parameters,thereby improving the controllability of the circuit.

    3. Symmetry sensitive to parameters

    In this section,the symmetric dynamical behavior which depends on the system parameters is analyzed. Firstly, based on the volt–ampere model (2), the symmetric coexistence bifurcation behaviors under special conditions are studied.Then, aiming at the dimensionality reduction model(12), the motions with corresponding parameters are simulated,and the differences between the above two models are compared.

    3.1. Symmetry in voltage–current domain

    For the parameters in Eq. (3) and the initial conditions(±10?9,0,0,0,0),the bifurcation diagrams and Lyapunov exponent spectrumsversus γandηare calculated and depicted in Fig.3,where the blue represents the coexisting bifurcation diagram with (10?9,0,0,0,0), red denotes the coexisting bifurcation diagram with(?10?9,0,0,0,0),two superposed bifurcation trajectories ofxmaxcan be observed. As the variation trend of Lyapunov exponent spectrums under two opposite initial conditions are almost the same,only the Lyapunov exponent spectrums at positive initial condition are exhibited in Fig.3.

    As can be seen from Fig.3,there are four symmetric motions, that is the following: ‘SP (stable point)’, ‘period(periodic state)’,‘CM(complex motion)’,and‘LP(large period)’.Forγincreasing in Fig. 3(a) andγ ∈(0.45,0.85), the system moves from SP to period, then to CM through perioddoubling bifurcation, and finally jumps to LP. Forηincreasing in Fig.3(b)andη ∈(0.65,1.95),the motion is exactly the opposite, starting from LP, then moving into CM and entering into period by reverse period-doubling bifurcation,finally dropping to SP. It should be noted that in the CM area, there appears mainly chaos with different periodic windows in the middle. Although the phase trajectories of the conventional period-1 and the large period are all limit cycles, while their Lyapunov exponents are different. When the large period appears,the minimum Lyapunov exponent will suddenly drop to an abnormally small value. In order to distinguish it between“period”and“LP”,the minimum Lyapunov exponent is used as a criterion.

    Fig.3. Bifurcation and Lyapunov exponent spectrum versus(a)γ and(b)η.

    To intuitively show the distribution characteristics of the dynamical behaviors for varying parametersγandη, the parameter mappings (γ,α) and (η,α) are depicted respectively in Fig. 4, where four colors are used to mark four motions: ‘blue-SP’, ‘orange-period’, ‘purple-CM’, and ‘yellow-LP’.From Fig.4,we can see that the symmetric characteristic still exist within the attraction basins,and it is not affected by other parameters. The whole evolution of “SP–period–CM–LP”, in positive or reverse order, can be fully presented for 1.4≤α ≤6. Ifαis less than 1.4,no matter how other parameters are changed,the circuit always goes through two states,SP and period,in which the system is extremely stable. Phase diagrams of each motion are given in Fig. 5 withα=2.66 and(10?9,0,0,0,0)as verification. When the specific chaotic signals need to be used,it is necessary to avoid fixing the parameter range within this range. It can been seen thatγandηhave good symmetry, for theGand 1/Lhave opposite trends in circuit evolution while other parameters have little effect on the symmetry.

    Fig.4. Parameter mappings with(10?9,0,0,0,0)on(a)γ–α plane and(b)η–α plane.

    Fig.5. Coexisting attractors with α =2.66 and(10?9,0,0,0,0).

    3.2. Symmetry in flux–charge domain

    After the dimensionality reduction by FCAM, the initial values are converted into the system parameters in flux–charge model(12), which can make the in-depth analysis of the heterogeneous dual memristive circuit. The original system (2)has symmetrical characteristics,which should not be changed by the dimensionality reduction model(12). To validate this,the bifurcation trajectory and Lyapunov exponent spectrums of 3OM are calculated and illustrated in Figs. 6(a) and 6(b),

    where the initial condition is (0,0,0), and the parameters are from Eq.(13).

    From Fig.6,we know that the diagrams,such as symmetric bifurcation and four motions,look similar to those in Fig.3,which means that Eq. (12) is consistent with system (1). Forγ1∈(0.45,0.5),g1=±10?9, andg2=g3=0, as shown in Fig.6(a),the corresponding maximum Lyapunov exponentL1is less than zero,and the system is in the stable state. With increasingγ1∈(0.5,0.85), the system enters into the period 1,then changes to period 2 through the period-doubling bifurcation,subsequently moves into complex motion through saddle junction bifurcation,in which several periodic windows exist,finally,suddenly jumps into a large period withL1returning to zero,andL2showing a steep drop. However,forηincreasing,η1∈(0.65,1.95), the movement trend of Fig. 6(b) is almost opposite to that of Fig.6(a).The system first moves from LP to chaos,and then enters into multi-period from the saddle junction bifurcation. Finally,it goes into period 1 through reverse period-doubling bifurcation and then remains stable.

    Comparing with the above analysis, it can be seen that the five-dimensional model in thev–idomain has been transformed into a three-dimensional model in the?–qdomain via FCAM.Although the system dimension decreases and the mathematical model changes,the dynamical behavior with the corresponding parameters does not change, and the special symmetric coexistence bifurcation still exists.

    Fig. 6. Bifurcation diagrams and Lyapunov exponent spectra for (a) γ1 ∈(0.45, 0.85)and(b)η1 ∈(0.65, 1.95).

    4. Multistability sensitive to initial conditions

    In this section, the multistability of the system with the initial values is analyzed based on the system model(Eqs.(2)and (13)), the similarities and differences among the multistable behaviors are studied through bifurcation diagrams and attraction basins.

    4.1. Bifurcation structures with initial conditions

    The different initial values each as a single variable will cause the system to exhibit bifurcation structures. The voltage–current model of the dual memristive system in Eq.(2)is described by the fifth-order equation,which has five initial conditions:x(0),y(0), andz(0) corresponding to the voltage ofC1,C2, and the current ofL, named “conventional initial value”;w(0)andv(0)corresponding to the internal variables ofM1,M2,as“memristive initial value”. To compare the influences of two kinds of initial values on multistability, the conventional initial values are firstly discussed.

    Fig. 7. Bifurcation diagrams with (a) (x(0), 0, 0, 0, 0), x(0)∈(?1.3,1.3), (b) (0, y(0), 0, 0, 0), y(0)∈(?0.5,0.5); and (c) (0, 0,z(0), 0, 0),z(0)∈(?0.7,0.7),respectively.

    Assume that the system parameters are from expression(3),then the bifurcation diagrams of the conventional initial values will be calculated and depicted in Fig. 7, where we can see that although the initial terms and the ranges of its variation are different,the system motion types look similar,i.e., the motion evolves from the LP into the chaos, then back to the LP. The structure of chaotic band in Fig. 7(a) is almost the same as that in Fig. 7(b), and the chaotic band in Fig.7(c)is like a mirror image of the structure of chaotic band in Fig. 7(a) or Fig. 7(b). The difference in the chaotic zone will cause chaotic attractors of different topological structures.Considering that three conventional initial values represent the voltage or current values of different components,it is believed that the difference in the structure of the chaotic zone is caused by the nature of each dynamic element. In addition, the corresponding Lyapunov exponents are also calculated and illustrated,where the curves are similar to those in Fig.8,since the Lyapunov exponents of the LP and chaos are (0,?,?,?,?)and(+,0,?,?,?),respectively.

    Fig. 8. Lyapunov exponent spectra with: (a) (x(0), 0, 0, 0, 0), x(0)∈(?1.3,1.3); (b) (0, y(0), 0, 0, 0), y(0)∈(?0.5,0.5); and (c) (0, 0,z(0), 0, 0),z(0)∈(?0.7,0.7),respectively.

    4.2. Multistability via voltage–current model

    In this subsection, the multistable characteristics of the system are revealed through the basin of attraction corresponding to the voltammetry model. The initial terms with different properties are selected as bivariate combinations to obtain the attraction basins as shown in Fig. 9, where dark blue, green and red color respectively represent the coexisting motions as“stable point”, “period”, and “complex motion” respectively.As the motions corresponding to the conventional initial values are similar, the motion ofxis taken for example. The periodic states include a variety of coexisting periodic states,such as coexisting left and right small period and large periodic limit cycles.

    As can be seen from Fig.9,the multistable diagrams symmetric with respect to the origin are observed,which is consistent with the system symmetric characteristic and the motion states withinx(0)∈(?1.3,1.3) in Figs. 9(a), 9(b), and 9(d)are the same as in Fig.9(a). The motion in Fig.9(a)shows a“reverse S”shape,and the state distribution in Fig.9(b)shows some distortions with respect to the origin. The chaotic attractors in Figs. 9(a) and 9(b) are mainly distributed near the origin, and the motions switch quickly. On the contrary, the distribution in Fig.9(d)is in the strip-type,i.e.,the system has a stable and wide chaotic band. Distinguishingly, the stable points disappear and the color boundary is fuzzy,which means that the system is more sensitive to the memristive initial values. From Fig.9(c), the symmetry of the motion distribution about the origin disappears in the attraction basin,meanwhile,the color boundary is more blurred and a stable fixed point appears,which indicates that the system that relies on the change of the initial value of dual memristance has more complex dynamical behavior.

    Fig.9. Attraction basins in voltage–current domain with(a)(x(0), 0, 0, w(0), 0);(b)(x(0), 0, 0, 0, v(0));(c)(0, 0, 0, w(0), v(0));and(d)(x(0), 0, z(0), 0, 0).

    In summary, the dual memristive system in Eq. (2) is more sensitive to the initial value of memristive resistance.When the conventional initial value is changed,the nonlinear behavior of the system is relatively simple,and it cannot fully reflect the multistable characteristics of the system. However,the change of the initial conditions of the memristive element causes the system to produce a variety of and even an infinite number of different types of attractors, which presents better multistable characteristics.

    4.3. Multistability via flux–charge model

    To validate the consistency between multistabilities before and after dimensionality reduction, the attraction basin in flux–charge domain is calculated by 3OM (12) and illustrated in Figs.10,12(a),and 13(a),in which eight coexistence attractor types are observed as listed in Table 1. Figure 10 shows that the symmetry with respect to the origin still exists in Fig.10(a),but the motion distribution is distorted about the origin in Fig.10(b). The initial values of circuit elements have different effects on the motion distribution,e.g., the initial value of memristor can make the system produce more coexisting structures. With the conventional initial values varying,the chaotic behaviors are observed within a larger parameter range,but the phenomenon of multistability cannot be fully reflected. Comparing with the charge-controlled memristorM2,the motion distribution of the flux-controlled memristorM1has a wide parameter range and concentrates near the origin,so the attractor structure is more stable.

    Fig.10. Two-dimensional initial value plane attraction basin for(a)uC1(0)–qM2(0)and(b)?M1(0)–qM2(0).

    Table 1. Types of attractors in different attraction domains.

    Fig.11. Chaotic attractors in planes of(a)x–z and(b)y–z.

    The existence of multistability needs to be verified by the trajectory of phase plane.The chaotic attractors of 3OM under parameters in Eq.(13)are given in Fig.11 with initial conditions(0,0,0),g1=10?9,andg2=g3=0,and the coexisting attractors are also shown in Figs. 12 and 13, where different kinds of attractors can be captured by attraction basins. In Fig. 12(a), the coexistence behavior presents an ‘S’ type. In Fig.13(a),the motion distribution looks like bands and strips,which is consistent with that observed in voltage–current domain. The evolution process of the coexisting attractor corresponding to Fig. 12(a) is given in Fig. 12(b). With the increase ofuC1(0), the system trajectory changes from a large period (type 3) to a left-side cycle 1 limit cycle (type 2), and then enters into the left-side chaotic state (type 6) after multiple period bifurcation. With the increase ofuC1(0), the left chaotic attractor changes into the right chaotic attractor(type 7), then it becomes the right cycle 1 (type 1) and eventually evolves back into a large cycle limit cycle(type 3)after multiple inverse cycle bifurcations. Similarly, the coexistent attractor change process corresponding to theuC2(0)–iL(0)twodimensional initial value plane is shown in Fig.13(b). WhenuC2(0) increases from?1 to 1, the system first stabilizes in the large-cycle state(type 3), then changes into a small cycle on the left(cycle 2,type 4),and then enters into the complex motion after bicyclic cycle bifurcation has evolved into a left chaotic attractor(type 6). Subsequently,repeating the foregoing process, the attractor structure changes again into a right chaotic attractor (type 7), and then enters into a small period state (right period 2, type 5) and eventually jumps to a large period(type 3). It should be pointed out that there is a transition state of cycle 1(types 1, 2)in the movement state of the system from large cycle to small cycle.

    Fig.12. Evolution process of corresponding coexisting attractors in attraction basin: (a)?M1(0)versus uC1(0)and(b)coexisting attractors.

    Fig.13. Evolution process of corresponding coexisting attractors in the attraction basin: (a)uC2(0)versus iL(0)and(b)coexisting attractors.

    Therefore, in contrast with 5OM (in Eq. (2)), by using 3OM(in Eq.(12)),more precise and effective discussion of infinite many coexisting attractors can be obtained for the transform from initial conditions to system parameters.The motion distribution trend of the circuit on the two-dimensional initial plane before and after dimension reduction are roughly the same,in which the multistable phenomenon shows the consistency.

    4.4. FPGA implementation of coexisting multistability

    In this subsection, equation (12) is transformed by the FCAM, which is beneficial to the study of multistability that depends on the initial values,and its hardware memristive circuit is realized via FPGA. Analog equivalent circuit is also a way to realize the chaotic system,[27]but the setting of initial conditions is not an easy task in practice.Such a digital experiment platform can set system parameters or initial values more quickly and accurately,and is suitable for implementing memristive chaotic circuits with higher precision requirements.[28]The digital implementation of the 3OM (12) is given as follows:

    The iteration stephis set to be 0.01, andKi j(i,j=1,2,3,4)are expressed as follows:

    Fig. 14. FPGA implementation, showing (a) connection diagram, (b)chaotic attractor in x–z plane.

    The physical connection diagram and the chaotic attractor are shown in Fig. 14. The trajectory is corresponding to that of in Fig.11(a),which verifies the correctness of numerical results. Subsequently, to reduce the multistability in the heterogeneous dual memristive circuit, a series of symmetrical coexistence attractors is presented in Fig.15,which physically realizes the unique multistability phenomenon as those shown in Figs.12 and 13. Therefore,the FPGA-based digital implementation of 3OM has high stability,good accuracy,and certain portability and universality.

    Fig.15. Coexisting attractors,showing(a)left-right-period 2,(b)left-right-period 3,and(c)left-right-chaos.

    5. Conclusions

    In this paper, a dual-memristor-based Chua’s circuit,which shows good symmetry and multistability,is constructed by introducing two memristors with different structures. To analyze its dynamical behaviors, such as the system mechanism and the influences of the initial conditions,two analysis models,i.e.,5OM in the voltage–current domain and 3OM in the flux–charge domain, are compared with each other. The symmetric bifurcation behaviors of the DMC are first investigated in the 5OM via parameter mappings to confirm the existence of symmetry under some special system parameters.By using the 3OM, a similar symmetry with corresponding parameters is also illustrated, which partly confirms the effectiveness and correctness of FCAM in the proposed circuit.Moreover,the multistabilities in two models are observed and compared to show the coexistence of complex multiple attractors in DMC. The reliability of FCAM is validated by both numerical simulations and FPGA implementation. For engineering applications,the analysis model selection shall depend on the practical situation.

    猜你喜歡
    張雷
    漲渡湖濕地冬韻
    A new stage of the Asian laser-induced breakdown spectroscopy community
    In-plane optical anisotropy of two-dimensional VOCl single crystal with weak interlayer interaction
    Measurement and analysis of species distribution in laser-induced ablation plasma of an aluminum–magnesium alloy
    黃科院田世民、呂錫芝、張雷入選水利青年拔尖人才
    人民黃河(2022年4期)2022-04-07 09:03:16
    張雷詠
    登銅雀臺
    Dynamic and inner-dressing control of four-wave mixing in periodically-driven atomic system?
    關(guān)于“見元る”的“自發(fā)”與“可能”
    耳邊不停的嘶喊
    女性天地(2018年3期)2018-03-27 09:55:22
    一级毛片我不卡| 两个人视频免费观看高清| 久久久午夜欧美精品| 亚洲电影在线观看av| 久久久久免费精品人妻一区二区| 最近最新中文字幕大全电影3| 精品久久久久久久久亚洲| 色尼玛亚洲综合影院| 国产精品,欧美在线| 免费观看在线日韩| 夜夜夜夜夜久久久久| 草草在线视频免费看| 99国产极品粉嫩在线观看| 日韩在线高清观看一区二区三区| 色av中文字幕| 热99re8久久精品国产| 午夜福利在线观看吧| 免费观看人在逋| 两性午夜刺激爽爽歪歪视频在线观看| 美女内射精品一级片tv| 大型黄色视频在线免费观看| 亚洲精品一卡2卡三卡4卡5卡| 久久午夜亚洲精品久久| 长腿黑丝高跟| 国产三级在线视频| 久久鲁丝午夜福利片| 亚洲不卡免费看| 一本一本综合久久| 亚洲熟妇中文字幕五十中出| 亚洲成人中文字幕在线播放| 熟妇人妻久久中文字幕3abv| 国产高清三级在线| 蜜桃亚洲精品一区二区三区| 中文字幕免费在线视频6| 麻豆av噜噜一区二区三区| 精品无人区乱码1区二区| 久久精品国产亚洲网站| 国产中年淑女户外野战色| 1024手机看黄色片| 久久九九热精品免费| 中国美女看黄片| 国产成人福利小说| 91午夜精品亚洲一区二区三区| 国产精品乱码一区二三区的特点| 大又大粗又爽又黄少妇毛片口| 尤物成人国产欧美一区二区三区| 午夜精品国产一区二区电影 | 91久久精品国产一区二区成人| 亚洲无线在线观看| 国产av不卡久久| 亚洲av免费在线观看| 免费不卡的大黄色大毛片视频在线观看 | 国产精品日韩av在线免费观看| 国产欧美日韩精品亚洲av| 国产欧美日韩精品亚洲av| 熟妇人妻久久中文字幕3abv| 国产视频内射| 色视频www国产| 午夜福利成人在线免费观看| 乱码一卡2卡4卡精品| 国产色婷婷99| 免费av观看视频| 联通29元200g的流量卡| 色尼玛亚洲综合影院| 亚洲天堂国产精品一区在线| 午夜福利在线观看吧| 色噜噜av男人的天堂激情| 天天躁夜夜躁狠狠久久av| 国产 一区 欧美 日韩| 久久6这里有精品| 噜噜噜噜噜久久久久久91| 国产成人91sexporn| 国产三级中文精品| 亚洲精品一区av在线观看| 蜜臀久久99精品久久宅男| 精品人妻视频免费看| 无遮挡黄片免费观看| 亚洲精品成人久久久久久| 最新中文字幕久久久久| 日本色播在线视频| 亚洲最大成人手机在线| 91麻豆精品激情在线观看国产| 神马国产精品三级电影在线观看| 中文亚洲av片在线观看爽| 国产精品乱码一区二三区的特点| 午夜福利成人在线免费观看| av视频在线观看入口| 久久99热这里只有精品18| 99精品在免费线老司机午夜| videossex国产| 亚洲熟妇中文字幕五十中出| 一a级毛片在线观看| 久久久久久九九精品二区国产| 免费无遮挡裸体视频| 亚洲在线观看片| 国产精品久久久久久久久免| 超碰av人人做人人爽久久| 国产精品国产三级国产av玫瑰| 亚洲成人久久性| 欧美日韩综合久久久久久| 99久久九九国产精品国产免费| 欧美最新免费一区二区三区| 1000部很黄的大片| 国产乱人视频| 亚洲真实伦在线观看| 日韩精品有码人妻一区| 我的女老师完整版在线观看| 国产精品日韩av在线免费观看| 男女之事视频高清在线观看| 亚洲四区av| 午夜激情福利司机影院| 日本撒尿小便嘘嘘汇集6| aaaaa片日本免费| 久久久久九九精品影院| 免费观看人在逋| 一夜夜www| 久久精品综合一区二区三区| 久久久精品欧美日韩精品| 亚洲av第一区精品v没综合| 久久久久国内视频| av视频在线观看入口| 舔av片在线| 丝袜美腿在线中文| 嫩草影院入口| 久久天躁狠狠躁夜夜2o2o| 国产亚洲av嫩草精品影院| 最后的刺客免费高清国语| 国产亚洲av嫩草精品影院| 欧美一区二区精品小视频在线| 人妻制服诱惑在线中文字幕| 久久久久国产精品人妻aⅴ院| 你懂的网址亚洲精品在线观看 | 伦理电影大哥的女人| 啦啦啦观看免费观看视频高清| 亚洲四区av| 校园人妻丝袜中文字幕| 国产高清视频在线观看网站| 亚洲一区二区三区色噜噜| 你懂的网址亚洲精品在线观看 | 波多野结衣高清作品| 97在线视频观看| 色播亚洲综合网| 午夜免费激情av| 精品无人区乱码1区二区| 亚洲人成网站在线播| 久久久久久久亚洲中文字幕| 伊人久久精品亚洲午夜| 在线观看一区二区三区| 国产老妇女一区| 午夜爱爱视频在线播放| 国产v大片淫在线免费观看| 成年免费大片在线观看| 亚洲人与动物交配视频| 1024手机看黄色片| 毛片女人毛片| 色av中文字幕| 看片在线看免费视频| 黄色日韩在线| 欧美+亚洲+日韩+国产| 久久久久久伊人网av| or卡值多少钱| 六月丁香七月| 99热全是精品| 国产精品久久电影中文字幕| 精品一区二区三区av网在线观看| 免费一级毛片在线播放高清视频| h日本视频在线播放| 日韩精品青青久久久久久| 亚洲精品一区av在线观看| 国产一区二区在线观看日韩| 国产精品久久久久久亚洲av鲁大| avwww免费| 国产黄a三级三级三级人| 两个人视频免费观看高清| 插阴视频在线观看视频| 国产人妻一区二区三区在| 秋霞在线观看毛片| 校园春色视频在线观看| 亚洲av成人精品一区久久| 久久亚洲国产成人精品v| 亚洲成人av在线免费| 亚洲激情五月婷婷啪啪| 国产私拍福利视频在线观看| 午夜亚洲福利在线播放| 国产真实乱freesex| 精品99又大又爽又粗少妇毛片| 成人综合一区亚洲| 亚洲国产精品合色在线| 91午夜精品亚洲一区二区三区| 久久久久久伊人网av| av在线天堂中文字幕| 午夜精品在线福利| 国产久久久一区二区三区| 男女那种视频在线观看| av在线老鸭窝| 国产伦精品一区二区三区四那| 黄色欧美视频在线观看| 国产成人一区二区在线| 国语自产精品视频在线第100页| 国产免费一级a男人的天堂| 精品熟女少妇av免费看| 成年版毛片免费区| 国产亚洲精品综合一区在线观看| 国产精品乱码一区二三区的特点| 久久欧美精品欧美久久欧美| 不卡视频在线观看欧美| 亚洲电影在线观看av| 神马国产精品三级电影在线观看| 香蕉av资源在线| 综合色av麻豆| 久久这里只有精品中国| 国产精品日韩av在线免费观看| 久久亚洲精品不卡| 欧美一区二区国产精品久久精品| 亚洲av.av天堂| 国产一区二区三区在线臀色熟女| 色吧在线观看| 久久久久久久午夜电影| 国内揄拍国产精品人妻在线| 亚洲精品国产成人久久av| 精品久久久久久久久av| 精品免费久久久久久久清纯| 久久国内精品自在自线图片| 欧美成人免费av一区二区三区| 99热这里只有是精品50| 国产探花极品一区二区| 精品人妻偷拍中文字幕| 搡女人真爽免费视频火全软件 | 久久人人爽人人爽人人片va| 嫩草影院入口| 麻豆av噜噜一区二区三区| 亚洲人成网站高清观看| 国产成人freesex在线 | 欧美zozozo另类| 乱系列少妇在线播放| 九九爱精品视频在线观看| 久久久久久久久中文| 三级男女做爰猛烈吃奶摸视频| 午夜a级毛片| 国产国拍精品亚洲av在线观看| 99在线视频只有这里精品首页| 国产伦精品一区二区三区四那| 91久久精品电影网| 色吧在线观看| 18禁在线无遮挡免费观看视频 | 国产毛片a区久久久久| 色吧在线观看| 久久6这里有精品| 亚洲人与动物交配视频| 最新在线观看一区二区三区| 欧洲精品卡2卡3卡4卡5卡区| 国产高清三级在线| 美女被艹到高潮喷水动态| 精品少妇黑人巨大在线播放 | 亚洲熟妇熟女久久| 国国产精品蜜臀av免费| 国产精品女同一区二区软件| 性插视频无遮挡在线免费观看| 简卡轻食公司| 国产单亲对白刺激| 老师上课跳d突然被开到最大视频| 看片在线看免费视频| 国内精品一区二区在线观看| 在线观看66精品国产| 国产三级在线视频| 99久久精品热视频| 九九久久精品国产亚洲av麻豆| 国产亚洲91精品色在线| 偷拍熟女少妇极品色| 日韩亚洲欧美综合| 久久久a久久爽久久v久久| 成人亚洲欧美一区二区av| 最新中文字幕久久久久| 18禁在线播放成人免费| 秋霞在线观看毛片| 午夜免费男女啪啪视频观看 | 99久久无色码亚洲精品果冻| 免费看日本二区| 91狼人影院| 97碰自拍视频| 一个人免费在线观看电影| 成人二区视频| 国产精品亚洲美女久久久| 九九在线视频观看精品| 日产精品乱码卡一卡2卡三| 久久婷婷人人爽人人干人人爱| 欧美性猛交╳xxx乱大交人| 国产亚洲av嫩草精品影院| 国产在视频线在精品| 亚洲性久久影院| 亚洲乱码一区二区免费版| 国产成人福利小说| 少妇猛男粗大的猛烈进出视频 | 中文亚洲av片在线观看爽| 黄片wwwwww| 最近中文字幕高清免费大全6| 精品无人区乱码1区二区| 亚洲三级黄色毛片| 欧美日韩综合久久久久久| 久久99热这里只有精品18| av卡一久久| 极品教师在线视频| 亚洲精华国产精华液的使用体验 | 精品一区二区三区人妻视频| av免费在线看不卡| 性色avwww在线观看| 国产成人影院久久av| 人妻少妇偷人精品九色| 国产视频内射| 国产色爽女视频免费观看| 国产色爽女视频免费观看| 国产黄片美女视频| 九九在线视频观看精品| 看非洲黑人一级黄片| 久久这里只有精品中国| 狠狠狠狠99中文字幕| 亚洲高清免费不卡视频| 极品教师在线视频| 日韩大尺度精品在线看网址| 亚洲色图av天堂| 亚洲精品日韩av片在线观看| 69人妻影院| 国产精品一二三区在线看| 天堂av国产一区二区熟女人妻| 国产乱人视频| 我要看日韩黄色一级片| 99热6这里只有精品| 国产亚洲精品av在线| 亚洲真实伦在线观看| 男女之事视频高清在线观看| 女的被弄到高潮叫床怎么办| 老熟妇乱子伦视频在线观看| 丰满的人妻完整版| 69av精品久久久久久| a级毛色黄片| 欧美最黄视频在线播放免费| 日本黄色视频三级网站网址| 欧美bdsm另类| 舔av片在线| 最近2019中文字幕mv第一页| 亚洲av电影不卡..在线观看| 成人美女网站在线观看视频| 又黄又爽又免费观看的视频| 国产乱人视频| 久久久久久伊人网av| 丝袜美腿在线中文| 久久精品国产99精品国产亚洲性色| 精品国内亚洲2022精品成人| 在线观看免费视频日本深夜| 91av网一区二区| 亚洲一级一片aⅴ在线观看| av在线天堂中文字幕| 午夜爱爱视频在线播放| 亚洲在线观看片| 人人妻人人澡人人爽人人夜夜 | 国产成人福利小说| 久久午夜福利片| 青春草视频在线免费观看| 国产精品一区二区三区四区免费观看 | 小说图片视频综合网站| 日韩欧美三级三区| 日本欧美国产在线视频| 成人欧美大片| 日韩精品青青久久久久久| a级毛片a级免费在线| 久久精品夜色国产| 日韩欧美在线乱码| 午夜影院日韩av| 两个人视频免费观看高清| 色尼玛亚洲综合影院| 精品少妇黑人巨大在线播放 | 成人特级av手机在线观看| 亚洲国产色片| 日韩大尺度精品在线看网址| 欧美+亚洲+日韩+国产| 亚洲熟妇熟女久久| 在线免费观看不下载黄p国产| 91在线观看av| 国产精品福利在线免费观看| 最近视频中文字幕2019在线8| 国产单亲对白刺激| 免费无遮挡裸体视频| 高清日韩中文字幕在线| 国产人妻一区二区三区在| 中国美白少妇内射xxxbb| 欧美日韩精品成人综合77777| 国产精品久久久久久av不卡| 最近最新中文字幕大全电影3| 一进一出抽搐gif免费好疼| 国产精品日韩av在线免费观看| 亚洲av电影不卡..在线观看| 插逼视频在线观看| 久久精品影院6| 久久国内精品自在自线图片| 国产精品久久久久久久久免| 波多野结衣巨乳人妻| 亚洲国产精品sss在线观看| 成人午夜高清在线视频| 国产人妻一区二区三区在| 中文字幕人妻熟人妻熟丝袜美| 中出人妻视频一区二区| 精品日产1卡2卡| 久久久久久久久中文| 你懂的网址亚洲精品在线观看 | 狠狠狠狠99中文字幕| 国产欧美日韩一区二区精品| 中文字幕人妻熟人妻熟丝袜美| 午夜日韩欧美国产| 一进一出抽搐gif免费好疼| 国产一区亚洲一区在线观看| 国内精品宾馆在线| 日韩大尺度精品在线看网址| 国产成人精品久久久久久| 日韩欧美免费精品| 在线播放国产精品三级| 1000部很黄的大片| 岛国在线免费视频观看| 男女那种视频在线观看| 此物有八面人人有两片| 69人妻影院| 天堂网av新在线| 男女下面进入的视频免费午夜| 禁无遮挡网站| 18+在线观看网站| 最新中文字幕久久久久| 午夜精品在线福利| 久久久久久久久久成人| 日韩欧美一区二区三区在线观看| 亚洲成人久久爱视频| 少妇熟女aⅴ在线视频| 午夜a级毛片| 亚洲精品久久国产高清桃花| 亚洲av成人精品一区久久| 天天躁夜夜躁狠狠久久av| 深夜精品福利| av女优亚洲男人天堂| 晚上一个人看的免费电影| 日日摸夜夜添夜夜添小说| 两个人的视频大全免费| 久久久精品欧美日韩精品| 晚上一个人看的免费电影| 黄色日韩在线| 免费不卡的大黄色大毛片视频在线观看 | 亚洲人成网站在线播放欧美日韩| 日日撸夜夜添| 亚洲欧美中文字幕日韩二区| 99精品在免费线老司机午夜| 在线播放无遮挡| 亚洲天堂国产精品一区在线| 免费av不卡在线播放| 久久久久久久久大av| 乱系列少妇在线播放| 色5月婷婷丁香| 高清日韩中文字幕在线| 秋霞在线观看毛片| 国产精品1区2区在线观看.| 99久久中文字幕三级久久日本| 午夜激情福利司机影院| 国产精品久久久久久久电影| 精品一区二区三区视频在线观看免费| 国产黄色小视频在线观看| 麻豆一二三区av精品| 久久精品国产清高在天天线| 热99re8久久精品国产| 99久久九九国产精品国产免费| 国产精品日韩av在线免费观看| 精品日产1卡2卡| ponron亚洲| 亚洲高清免费不卡视频| 久久久久久久久中文| 国产精品免费一区二区三区在线| 色视频www国产| 国产一级毛片七仙女欲春2| 国产黄色小视频在线观看| 淫妇啪啪啪对白视频| 日韩欧美精品v在线| 91久久精品国产一区二区成人| 久久精品国产99精品国产亚洲性色| 欧美日韩国产亚洲二区| 国产午夜精品久久久久久一区二区三区 | 欧美国产日韩亚洲一区| 免费观看在线日韩| 国产成年人精品一区二区| 日韩欧美精品免费久久| 又爽又黄a免费视频| 在线观看66精品国产| 1000部很黄的大片| 久久国产乱子免费精品| 床上黄色一级片| 两个人视频免费观看高清| 中出人妻视频一区二区| 又黄又爽又刺激的免费视频.| 91在线观看av| 99久久久亚洲精品蜜臀av| 欧美bdsm另类| 日日摸夜夜添夜夜爱| 午夜激情福利司机影院| av黄色大香蕉| 天天一区二区日本电影三级| 婷婷色综合大香蕉| 精品无人区乱码1区二区| 国产精品国产高清国产av| 久久人人爽人人爽人人片va| 成年av动漫网址| 久久国产乱子免费精品| 韩国av在线不卡| 不卡一级毛片| 午夜老司机福利剧场| 中国美白少妇内射xxxbb| 国产男人的电影天堂91| 高清毛片免费看| 五月伊人婷婷丁香| 97超级碰碰碰精品色视频在线观看| 国产一级毛片七仙女欲春2| 免费搜索国产男女视频| 日韩国内少妇激情av| 久久久精品94久久精品| 精品一区二区三区av网在线观看| 高清日韩中文字幕在线| 午夜精品国产一区二区电影 | 色综合亚洲欧美另类图片| 欧美中文日本在线观看视频| 一进一出抽搐gif免费好疼| av卡一久久| 真人做人爱边吃奶动态| 日韩欧美三级三区| 干丝袜人妻中文字幕| 免费观看在线日韩| 欧美性猛交黑人性爽| 久久久久久久久久久丰满| 亚洲乱码一区二区免费版| 麻豆国产av国片精品| 日本成人三级电影网站| 久久综合国产亚洲精品| 天堂网av新在线| 国产伦精品一区二区三区视频9| 麻豆精品久久久久久蜜桃| 国产午夜福利久久久久久| 欧美在线一区亚洲| 91在线观看av| 亚洲精品粉嫩美女一区| 最近最新中文字幕大全电影3| 国产不卡一卡二| 最新在线观看一区二区三区| 欧美一区二区精品小视频在线| 亚洲精品国产av成人精品 | 精品久久国产蜜桃| 直男gayav资源| 欧美日本亚洲视频在线播放| 1024手机看黄色片| 免费av不卡在线播放| 两个人视频免费观看高清| 18+在线观看网站| 九色成人免费人妻av| 精品人妻一区二区三区麻豆 | 搡老岳熟女国产| 麻豆国产av国片精品| av.在线天堂| 毛片女人毛片| 在线天堂最新版资源| 国模一区二区三区四区视频| 变态另类成人亚洲欧美熟女| 欧美xxxx黑人xx丫x性爽| 国语自产精品视频在线第100页| 岛国在线免费视频观看| 久久久久国产精品人妻aⅴ院| 亚洲aⅴ乱码一区二区在线播放| 国产探花极品一区二区| 免费大片18禁| 国产精品99久久久久久久久| 一夜夜www| 国产精品国产高清国产av| 最近中文字幕高清免费大全6| 免费看光身美女| 亚洲七黄色美女视频| 一个人看视频在线观看www免费| 成人漫画全彩无遮挡| 免费一级毛片在线播放高清视频| 亚洲国产欧洲综合997久久,| 丝袜喷水一区| 亚洲综合色惰| 一卡2卡三卡四卡精品乱码亚洲| 国产精品福利在线免费观看| 国产日本99.免费观看| 大型黄色视频在线免费观看| 午夜爱爱视频在线播放| 99久久久亚洲精品蜜臀av| 亚洲国产精品成人久久小说 | 色5月婷婷丁香| 深爱激情五月婷婷| 欧美日韩国产亚洲二区| 超碰av人人做人人爽久久| 国产精品99久久久久久久久| 日本a在线网址| 久久鲁丝午夜福利片| 亚洲七黄色美女视频| 亚洲av五月六月丁香网| 男女视频在线观看网站免费| 又黄又爽又免费观看的视频| 久久韩国三级中文字幕| 伦理电影大哥的女人| 久久午夜福利片| 欧美日韩一区二区视频在线观看视频在线 | 嫩草影院新地址| 色综合色国产| 亚洲五月天丁香| 你懂的网址亚洲精品在线观看 | 国产亚洲精品久久久久久毛片| 亚洲aⅴ乱码一区二区在线播放| 国产国拍精品亚洲av在线观看| 亚洲在线自拍视频| 国产成人91sexporn| 日本在线视频免费播放| 亚洲av.av天堂| 精品久久久久久久久久久久久|