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

    Bifurcation analysis of visual angle model with anticipated time and stabilizing driving behavior

    2022-08-01 06:01:16XueyiGuan管學(xué)義RongjunCheng程榮軍andHongxiaGe葛紅霞
    Chinese Physics B 2022年7期
    關(guān)鍵詞:榮軍紅霞

    Xueyi Guan(管學(xué)義), Rongjun Cheng(程榮軍), and Hongxia Ge(葛紅霞)

    Faculty of Maritime and Transportation,Ningbo University,Ningbo 315211,China

    Keywords: visual angle,bifurcation analysis,anticipated time,stabilizing driving behavior,TDGL and mKdV equations

    1. Introduction

    In recent years, the demand for transportation has maintained a continuous growth, which has brought huge challenges to the smooth traffic and traffic safety. Facing the traffic problems,numerous of scholars proposed a series of traffic flow models to explain the complex traffic phenomena. These models[1–9]mainly model the traffic flow from three modules:macroscopic,microscopic and mesoscopic. For example,Sunet al.[1]proposed a car-following model, one of the microscopic models of vehicular traffic, which is concentrating on individual vehicles,and studying the characteristics of the following vehicles in the entire fleet by studying the influence of the changes in the movement state of the leading vehicle on the following vehicles, meanwhile the theoretical analysis and simulation of the improved model not only can provide deep exploring into the properties of the model but also help in better understanding the complex phenomena observed in real traffic fields. In addition,we can also subdivide these models into the car-following models (e.g., Sunet al.,[1]Geet al.[2]and Jianget al.[3]), the cellular automation models (Tianet al.[4]), the gas kinetic models (Nagataniet al.[5]), the hydrodynamic lattice models (Redhuet al.[6]and Changet al.[7])and the continuum models(Chenget al.[8]and Zhaiet al.[9])according to different modeling theories. Bifurcation, as an important research content of non-linear science,mainly studies the instability of the system solution caused by parameter changes and the changing behavior of the number of solutions.Studying the bifurcation in traffic flow will have important theoretical implications in analyzing traffic phenomena, traffic modeling,traffic control,etc.

    Facing numerous of traffic congestions, a number of mathematical models have been proposed to demonstrate the car-following behavior under various conditions with the development of urban transportation. Generally speaking,these models are mostly based on the stimulus-response framework,which was put forward at the General Motors research laboratories (Chandleret al.[10]and Gaziset al.[11]) at the first time. According to the framework,each driver responds when giving a stimulus. And the relationship between response and stimulus is as follows:

    In the framework, the stimulus could be mainly divided into four aspects: driver’s behavior factors, vehicle factors, road factors and other factors. However, when it comes to carfollowing model equation,the stimulus could be related to the location of a large number of cars or their speed, or maybe some other coefficient. The reaction which demonstrates the propulsion of follower has been transformed into vehicle’s speeds or accelerations. In different forms of car-following models, the reason for the different forms of vehicle motion equations lies in the different hypotheses of stimulus properties. The stimulus is probably composed of the speed, the speed difference between two consecutive vehicles,the headway, etc. Following are some classic car-following models that consider different stimulus.

    Gazis–Herman–Rothery (GHR)[11]car-following model is perhaps the most recognized model in the research of traffic flow, which is the foundation of the car-following model proposed by Hermanet al.[12]and Yanget al.,[13]displayed as

    whereλis the sensitivity of the relative velocity, Δvn(t-τ)represents the relative velocity of thenth car at timet-τ,τis reaction time of driver,in addition,andλhas some functional forms as follows:

    where Δxn(t)represents distance between two consecutive vehicles, Δxcritical(t) is the threshold value obtained from actual data analysis,vn(t)represents the velocity ofnth vehicle,whileA,A1andA2are constant. Moreover, Lee[16]believed that drivers make decisions relying on a period of stimulation rather than a transient stimulus,and for the reason he proposed an improved GM model with considering driver’s memory,which is the foundation of the extended car-following model considering driver’s memory and jerk proposed by Jinet al.[17]

    For whether Chandleret al.[10]or Gaziset al.[11]mentioned earlier, they all assumed that both cars are driving at the same speed and the headway is acceptable, but this is not consistent with the actual situation. To overcome this shortcoming,Newell[18]proposed a nonlinear car-following model by using the headway as a stimulus, and Helly[19]assumed that every driver will adjust the following distance to a desired distance according to his driving characteristics, and drivers strive to minimize the relative speed along with the distance between the actual headway and the required headway,which is defined as

    whereDn(t) is driver’s desired following distance. In view of the earlier work, Helly[19]calculated the final calibration results whereτ=0.5 s,A1=0.5 s-1,A2=0.125 s-2, andDn(t)=20+vn(t-0.5). And Helly’s improvement also laid the foundation for the innovation of Koshiet al.[20]and Xinget al.[21]For example, Koshiet al.[20]combined Helly’s model with GHR’s model and did a nonlinear extension, but it still has some shortcomings in horizontal curvature effect.

    Then in 1995, Bando[22]proposed the optimal velocity(OV) model, which describes some nonlinear phenomenon and car-following behavior in traffic, and the specific expression shows as following:

    whereαdonates the coefficient of sensitivity andVrepresents the optimal velocity that drivers prefer based on the condition of the headway.

    Optimal velocity model can reflect some actual traffic flow characteristics such as the disappearance of stability of continuous traffic flow, the queuing of the traffic and the entire process from the generation to the elimination of traffic congestion very well, but when facing high acceleration values, it has certain defects. Helbing and Tilch[23]compiled a large amount of empirical data to overcome this shortcoming and proposed a generalized force (GF) model. The same as GHR model, GF model considered the effect of relative speed, too, but the model ignored the situation where the velocity of the follower is equal to or less than the velocity of the leader. Then, Jianget al.[24]proposed the full velocity difference(FVD)model,which considers the positive relative velocity. Therefore, the FVD model can describe the actual traffic flow better compared to the OV and GF models during the car-following process. Its formula reads

    in whichV(θn(t))represents the relative velocity andλis its sensitivity coefficient. In 2019, Yanet al.[25]improved the FVD by taking the optimal velocity difference and the electronic throttle angle difference of the preceding vehicle into account so as to explain the current complex traffic situation.

    In the model, driver’s response is mainly coming from the stimulus of both speed difference and headway. However,since Michaels[26]introduced psychology-related theories into the car-following model for the first time,a series of follow-up studies (e.g., Winsum[27]) showed that psychological characteristics are indeed related to the stability of traffic flow. Its main performance is that visual angle can be used to replace headway together with angular velocity, which can be used to replace speed difference in some engineering car-following models. Therefore,Andersen and Sauer[28]proposed the driving by visual angle (DVA) model by taking visual angle into account based on Helly’s[19]model.

    Combined the FVD with DVA model,Jinet al.[29]developed the visual angle model(VAM)which takes both positive velocity difference and visual angle into account defined as

    whereV(θn(t)) represents the optimal velocity according to the visual angle of the driver, andθn(t) represents visual angle observed by the driver of carnat timet,as shown in Fig.1.

    As we all know, car-following behavior is a humanrelated process, and inserting visual angle into the FVD can access to 3D parameters information directly, which is closer to reality. However, in order to reduce gasoline consumption and stop-and-go while driving,many drivers always try to keep the vehicle as smooth as possible, which means that drivers will control the headway far enough to accelerate, but in order to avoid the repeated deceleration, the driver may drive smoothly at the original speed rather than accelerating. Therefore, in addition to the influence of visual angle and relative speed,driver’s stabilizing driving behavior can also affect the driving behavior, and for this reason, based on the VAM, in this paper we propose an improved car-following model by taking the effect of visual angle,anticipated time and driver’s stabilizing driving behavior into account.

    Bifurcation is defined as the process in which the nonlinear dynamic parameters of the system continuously change and reach a certain critical point,which causes the global nature of the system to change. In 1978, Huilgol[30]firstly introduced the bifurcation theory to the nonlinear dynamics of vehicle systems,which aroused the enthusiasm of traffic scholars around the world to study the stability of vehicle systems.With the continuous development of bifurcation analysis theory,Igarashiet al.[31]and Oroszet al.[32]selected different optimal velocity functions to study the bifurcation phenomenon in the OV model,and obtained the similar conclusions,which was later extended to the OVM with a single fixed delay by Yamamotoet al.[33]performed Hopf bifurcation analysis for the dissipative system with asymmetric interaction including the OVM. It is generally found that a bifurcation structure is related to the multiple traffic jams via the transitions of emerging clusters(i.e.,traveling waves). In order to further explore the relationship between traffic state and bifurcation,in 2004,Gasseret al.[34]conducted a series of car-following models with bifurcation analyses. They found that the instability of the neutral stable solution increased in optimal velocity models is normally owing to the happen of the Hopf bifurcation.Jinet al.[35]conducted a bifurcation analysis for the FVD model to study the effect of bifurcation behavior on traffic congestion and stability of the system. More recently, Zhanget al.[36]studied the traditional velocity difference model considering time delays. They found that bifurcation happened corresponding to the emergence of traffic congestion in continuous traffic flow. Nowadays,a great progress on traffic flow has been made. However, it is still essential to explore and grasp the nonlinear characteristic of traffic information so as to eliminate the traffic jams.

    Through the review of the development of traffic flow theory,we can find that few studies have comprehensively considered the driver’s visual angle,anticipated time and driver’s stable driving behavior to explore the impact of these characteristics on car-following behavior. While the article focuses on three human factors, comprehensively considers the driver’s visual angle, anticipated time and stable driving behavior to explore the impact of driver behavior factors on car-following behavior. Then the influence of various factors on the stability of traffic flow is obtained through the linear, nonlinear analysis and numerical simulation. Finally, the rationality of the improved model is verified by parameter calibration.

    The specific contribution of this paper are as follows:

    (1)Providing theoretical support for the establishment of intelligent transportation systems. In the process of developing self-driving vehicles,robots can be used to simulate human driving,and the optimal human psycho-physiological characteristics are set during the robot programming,so as to achieve the maximum safety and optimal efficiency of the traffic flow.

    (2)Providing references for traffic management systems.In terms of traffic safety, the optimal driver’s psychologicalphysical characteristics obtained from the study can be used as the reference standard to strictly review the driver’s driving qualification.

    The following is the framework of the article. First of all,we make an overview of the development of the traffic flow in Section 1.Then we take visual angle,anticipated time and stabilizing driving behavior into account and put forward a new visual angle model in Section 2. In Section 3,the neutral stability curve is carried out through linear stability analysis. In Section 4 we conduct bifurcation analysis for the improved model and study the influence of the three driver’s behavior on Hopf bifurcation. Through the numerical calculation, the time-dependent Ginzburg–Landau (TDGL) equation and the modified Korteweg–de Vries (mKdV) equation are given respectively in Sections 5 and 6. Then in Section 7, parameter calibration and numerical simulations are carried out and numerical results are discussed. Finally, we obtain the conclusions in Section 8.

    2. The extended car-following model

    All of the previous car-following models are developed based on the ideal state that the driver’s immediate response to traffic conditions is absolutely rational,meanwhile,drivers have a great ability to perceive the traffic information,such as headway, velocity and acceleration. However, there are few psychological studies can prove this. Although some basic factors affecting the driving behavior,such as time delays,inaccuracies,and anticipated time have been taken into account in the car-following model, the headway and the relative velocity always play the most direct role in driver’s perception.The improved model in this paper is proposed based on the theory,considering from a psychological perspective,that numerous of studies have showed that psychological effects can affect driving behavior, and the results of these studies can be applied to car-following models. The most essential research results are that visual information is strongly related to the driver’s driving.

    Considering both the physical and psychological characteristics of follower, we insert a visual angle of follower’s driver into the car-following model rather than the headway and the difference of speed,whose relationship is displayed in Fig.1,

    whereθn(t)is defined as the visual angle which is measured by the driver of thenth car at timet,Δxn(t)=xn+1(t)-xn(t)is the headway between the preceding (n+1)th car and the followingnth car,andln-1andwn-1represent the length and the width of the(n-1)th car,respectively.

    Fig.1. Car-following model considering the visual angle.

    In view of the above knowledge, we developed an improved car-following model considering the visual angle and anticipated time with stabilizing driving behavior.The dynamical equation is proposed as follows:

    whereαis known as the sensitivity coefficient of drivers to the difference between the current speed and desired speed.vn(t)denotes the speed of carnat timet,εandγare defined as the sensitivity parameter of time-delayed velocity difference and steady-state control effects respectively, andτis the anticipated time. Furthermore,Vrepresents the optimal velocity whose mathematical expression is defined as follows:

    whereve=maxV(x)>0 is called the desired speed.

    Equation (7) proved that the visual angle stimulus can really affect the response of driver to follower. Since carfollowing behavior is closely related to human decisionmaking and response,it needs to take human’s limits on judging headway and velocity stimulus into account. Therefore,considering driver’s visual angle factor can fit to the real traffic flow better.

    3. Linear stability analysis

    In this section, we conduct the prerequisites for the stability of the improved model by linear analysis. For a uniform traffic flow,all vehicles have the same characteristics,including everyone in the unified system is moving at the optimal velocityV(θ0)and the same headwayh,and the visual angleθn(t)remains the same. Obviously,we can obtain

    Integrating the two formulas, we can obtain the neutral stability conditions

    With small disturbances,the neutral stability condition of the uniform traffic flow is as follows:

    Whenαis satisfied with the condition, the traffic is stable with small disturbances.Otherwise,the traffic will become unstable.

    Fig.2. Neutral stability curve with different γ,w,ε and τ.

    Figure 2 depicts the neutral stability curve which describes the functional relationship between the headway and sensitivity coefficient space. As we can see,every neutral stability curve has a critical point (hc,αc) with differentγ,w,εandτ, respectively. The stability graph is separated into two parts by the critical curve, among which, one represents stability region,and the other represents instability region. Each point on the line represents a critical point for the stability to instability region. Whenα >αc, regardless of the headway,traffic flow will always remain stable, otherwise congestion may occur. In other words, the traffic flow will keep going above the neutral stability curve and congestion is impossible to happen. While under the curve,the traffic flow density increased along with the increased possibility of congestion.As shown in Figs. 2(a), 2(c) and 2(d), when fixing other parameters at a proper value,the vertices of the curve gradually decrease with the increases ofγ,εandτrespectively. However,they play different roles for traffic flow stability. Instead,F(xiàn)ig.2(b)describes the instability regions increase with the increasing ofwgradually when fixingε=0.2 s-1,γ=0.2 m-1andτ=0.5 s. From this, we conclude thatγ,εandτhave a positive effect in the traffic stability,which is opposite tow.

    4. Bifurcation analysis and simulation

    In this section,we use the theoretical analysis to explore the bifurcation behavior of the new model first of all so as to study the correspondence between the critical point of Hopf bifurcation and instability, and then conduct the numerical simulation to display the impact of various parameters on the bifurcation behavior intuitively. Equation(7)can be rewritten as

    Through analyzing Eq.(17),we obtain the linearized equation as follows:

    Use the Hopf theorem to find the Hopf bifurcation pointL*.It can be found that when a smooth family of eigenvaluesλ=ψ±iωappears, the Hopf bifurcation condition is satisfied,and the bifurcation occurs.

    Substitutingλ=ψ±iωinto Eq.(19),one obtains

    wherek ∈{1,2,...,N}.

    Separating the real and imaginary parts of Eq. (22), the result shows

    When Hopf bifurcation happens, there are only a pair of imaginary rootsλ1,2=±iω. In Fig. 3, all of other (N-2)eigenvalues are located in the left complex plane when a pair of pure imaginary roots appears at the first time(fork=1 andk=N-1).

    Fig.3. Eigenvalues of the system for N=5,α =1 s-1,ε =0.05 s-1,and τ =0.4 s.

    which is just contrary to the result shown in Eq.(15).

    Next, differentiation of the characteristic Eq.(27)brings about the essential condition for the Hopf bifurcation,i.e.,

    Obviously, from Eq. (27), the Hopf conditionψ(β1)<0 is satisfied absolutely. Therefore,the Hopf bifurcation occurs in our system. We can obtain the critical lengthL*from the following equation:

    The spatial variation of Hopf’s relation(28)is displayed in the three-dimensional diagram of Fig. 4. The yellow surface is graph of

    while the blue surface is the spatial variation curve of

    The intersection of these two surfaces corresponds to the critical point when Hopf bifurcation occurs.

    In order to display and compare the effect ofτ,εandγon Hopf bifurcation more intuitively, Fig. 5 shows the two-dimensional(L,N)-plane corresponding to Fig.4,respectively,which reveals the instability regions of the system,and in which the first pair of eigenvalues have positive real parts drawn in yellow. The Hopf curve consists of the set of all intersections(L,N)of the first pair of eigenvalues and the second pair of eigenvalues, that is, the boundary between the blue plane and the yellow plane in Fig. 4. It is uncovered that the sensitivity anticipated timeτ,εandγplay a different role on the occurrence of Hopf bifurcation. Figures 5(a) and 5(b) indicate the critical point (LC,NC) of Hopf bifurcation from(4.1,3.8)to(5.1,8.9)because of the increasing of anticipated timeτfrom 1 to 2 and fixing the sensitivity coefficientε=0.05 s-1andγ=0.8 m-1.In Figs.5(c)and(d),the critical point(LC,NC)of Hopf bifurcation changes from(4.1,3.8)to(5.4,10.8)with the increasing ofεwhile keeping anticipated timeτandγunchanged. It displays the instability area decreases because of the increasing ofε. In Figs.5(e)and 5(f),the critical point (LC,NC) of Hopf bifurcation changes from(5.8,11.6) to (3.9,1.8) with the decrease of the sensitivityγand fixing the anticipated timeτ=0.4 s andε=0.05 s-1.It shows the instability area increases with the decrease ofγfrom 0.8 m-1to 0.65 m-1. Therefore, parameterτ,εandγall can inhibit the occurrence of bifurcation,it means they all contribute to the system’s stability despite the different effect in the critical point region of Hopf bifurcation.

    Fig.5. Side view for α=1 s-1,ve=1 m/s with different values γ,ε and τ: (a)τ=1 s,ε=0.05 s-1,γ=0.8 m-1;(b)τ=2 s,ε=0.05 s-1,γ =0.8 m-1; (c) τ =1 s, ε =0.05 s-1, γ =0.8 m-1; (d) τ =1 s, ε =0.1 s-1, γ =0.8 m-1; (e) τ =0.4 s, ε =0.05 s-1, γ =0.8 m-1; (f)τ =0.4 s,ε =0.05 s-1,γ =0.65 m-1.

    5. The TDGL eqnarray

    In this section, we carry out the nonlinear analysis of Eq. (7) so as to obtain the TDGL equation that fits the improved model traffic behavior best near the critical point(hc,αc). The slow scales for space variablen, time variabletand constantbare introduced,and the slow variablesXandTare defined as follows:

    While Δxn(t)is defined as follows:

    wherehcstands for safety distance. By inserting Eqs.(29)and(30)to Eq.(31),along with the Taylor expanding to the fifthorder ofμ,the nonlinear partial differential equation is shown as follows:

    and substitute Eq. (35) into Eq. (34), the TDGL equation is transformed as

    where?Φ(S)/?Srepresents the function derivative. In addition to the simplest solutionS=0,another two stable solutions have also been shown,one of which is a uniform solution

    and the other is the kink solution

    wherex0is constant. Through analyzing Eq.(39),we can obtain a coexistence curve when following formula relations are met:

    Inserting Eq.(35)into Eq.(40),we can obtain the coexistence curve with the original parameters

    The spinodal line is based on the following conditions:

    By inserting Eq.(35)into Eq.(42),the relationship of the spinodal line is as follows:

    Based on the condition?φ/?S=0 and Eq.(43),we can calculate the critical point containing the original parameters

    6. The mKdV equation

    All ofgiare shown in Table 1.

    Table 1. The coefficients gi.

    In the table

    In order to export the regularized equation,we give the following definition:

    So,the standard mKdV equation withO(μ)correction term is obtained as

    Then,assuming thatR′(X,T′)=R′0(X,T′)+μR′1(X,T′)and considering theO(μ) correction term, in order to obtain the propagation velocitycfor the kink solution,it can be based on the solvability condition as

    where

    and the general velocitycobtained by calculation is as follows:

    whereV′′′<0, andτrepresents the reciprocal of the driver’s sensitivity coefficient.

    7. Simulation

    7.1. Parameter calibration

    In this part, 558 groups datasets of the 59th vehicle in the NGSIM database including velocity,acceleration and other factors are used to calibrate the optimal parameters. Considering that the modified dataset includes vehicle data in a noncar-following state, and existing the lane-changing behavior,therefore we make the following restrictions on car-following behavior.

    (1)Car-following behavior should occur in the same lane,in other words, the data bar of Lane-ID should be the same value.

    (2) Except for the first car, the selected leading car and following car,namely,preceding and following are not 0.

    (3) In order to ensure that the vehicle stays in the carfollowing state, the headway should not be too large. Otherwise,the behavior change of the previous car will hardly affect the behavior of the target car.

    Filtering by the above conditions,we can obtain the carfollowing data about the 59th car. Following is a partial data display.

    Table 2. Car-following partial data of the 59th car.

    According to the model in this article,we can determine the parameters to be calibrated in this article and the input variables of the car-following model, then using the least square method to perform nonlinear optimization of the calibration parameters. The specific objective function expression we adopt is proposed by Ossenet al.,[37]whose core idea is as follows:

    whereNis the data group number,jdonates the serial number of data,andrepresent actual acceleration and simulated acceleration respectively. Meanwhile, as an evaluation index describing the fitness of the parameter calibration, the smaller the PI, the closer the measured value is to the simulated value, that is, the better the applicability of the fitting result.

    Moreover, both models select parameter combinations optimized by genetic algorithm and within 5% of the error range as their analysis objects.Then normalize each parameter and systematically cluster the optimized processing results of each parameter,the comparison of each parameter distribution of the clustering results of the two methods in Table 3.

    As we can see in Table 3 that considering the anticipated time and stabilizing driving behavior can decrease the influence ofαwith the same data, which means anticipated time and stabilizing driving behavior can eliminate or weaken the adverse effects of optimal velocity in unified traffic flow effectively. In addition, the smallerhand the increase ofvmaxindicate that the minimum safety headway requirement is smaller in the stable traffic flow,and the maximum velocity that each vehicle can reach in the traffic flow has increased,which means the extended VAM is easier to steady than the VAM.

    Table 3. Calibrated parameters.

    Figure 6 displays the fitness between the simulations and observations. It is not difficult to see that the red and pink curve maintain roughly the same trend with the blue line,where the pink and red lines denote the fitted value of extended VAM and VAM respectively, and the blue line represents the actual observation, indicating that the initial calibration coefficient set of the VAM and the extended VAM both can effectively predict the movement trend of the actual traffic flow during the acceleration,uniform velocity and deceleration phases.However,comparing the red and pink lines,obviously,the fitted value set of the extended VAM is closer to the observation in most data despite larger deviations in a few inflection points. And synthesizing the error value of all the data, it is found that the PI value is reduced from 0.6624 to 0.6237 after considering the gain and time-delay feedback control, which means that the controlled visual angle model can better reflect the complex actual flow.

    7.2. Numerical simulation

    In this section, numerical simulation of the improved model has been conducted so as to prove the obtained conclusions by using discrete equations. The discrete process is as follows:

    In Fig. 7, the improved car-following model carries out a series of simulations while controlling other parameters unchanged. The patterns in Figs. 7(a)–7(d) correspond to the variablew(w= 1.6,1.8,2.0,2.2) respectively. Especially,whenw=1.6 m,the traffic flow has the highest stability that means the car width of the front car has a certain influence on the follower’s visual angle. As we can see that giving a small perturbation can induce initial stable traffic flow into an oscillating density wave. It can be clearly seen from the curvature change in Fig. 7 that the amplitude of the kink–antikink solution gradually increases with the increasing of parameterw.In other words,taking the vehicle width of the leader into account will interfere with the follower’s perspective,which will inevitably affect the stability of the system and increase the possibility of congestion at the same time. The wider the vehicle’s width,the more it can interfere with the driver’s perspective, and the greater the possibility of the system’s transition to instability region and traffic jams.

    Figure 8 displays that the headway profiles and density waves for the preceding vehicle widthwatt=6000 s corresponding to panels in Fig. 7 respectively. Obviously, as the width of the vehicle increases, the amplitude of the headway gradually increases,which means considering the width of the leader will restrain the stability of the traffic flow and prolong the time for the unstable traffic flow to stabilize.

    Fig.7. Space–time evolution of the headway from with the t=5000 s to 6000 s.

    Fig.8. Headway profiles of the density wave at t=6000 s correspond to Fig.6 respectively.

    Fig.9. Space–time evolution of the headway from with the t=4400 s to 4500 s.

    In Fig.9,we explore the effect of the anticipated time on the stability of traffic flow in the improved car following model through a series of simulations, where Figs. 9(a)–9(d) correspond toε(ε=0.1,0.3,0.5,0.7),respectively. It can be obviously seen from the change curvature in Fig.8 that considering the anticipated time contributes to stable movement of the system. Specially,whenε=0.7 s-1,the density wave change curvature is zero, at this time, the traffic flow steadily travels with a constant headway,and the entire traffic flow system has the highest stability. Therefore,we conclude that considering the driving characteristics of subject driver can effectively improve the stability of the traffic flow. The proposed new model considers the anticipated time and the subject driver uses the received original traffic information combines with his own prediction of the front traffic state to shorten the time responds to the front traffic situation, which provides a comprehensive improvement to the entire car-following behavior.

    Figure 10 displays the headway profiles and density waves for anticipated time att=4500 s corresponding to panels in Fig. 8 respectively. From Figs. 10(a)–10(d), considering the speed difference,the traffic congestion conditions can be effectively alleviated,and as the anticipated time increases,the amplitude of the density wave gradually decreases,which means the traffic flow tends to move with a stable headway.The improved model has improved to a greater extent in eliminating interference. From Fig.10,we obtain that parameterεcan promote the stability of traffic flow.

    Fig.10. Headway profiles of the density wave at t=4500 s correspond to Fig.8 respectively.

    Figure 11 mainly depicts the hysteresis loop relationships between the velocity and density, which are shown with deferentγ(γ=0.05, 0.1, 0.15, 0.2)respectively. As we can see the size of the loop gradually decrease with the value increasesγin Fig. 11. This also means that considering driver’s stabilizing driving behavior can effectively ease traffic congestion and promote the stability of the system at the same time.

    Fig. 11. The hysteresis loops relations of the velocities v1(t) and the density ρ1(t)=1/Δx1(t)respectively.

    8. Conclusion

    In this paper, we take the anticipated time and stabilizing driving behavior into the visual angle model, combining anticipated time with stabilizing driving behavior to explore driver’s sensitivity of change of angle information. The relationship of neutral stability is obtained by linear stability analysis. Meanwhile, bifurcation analysis for the improved carfollowing model considering the visual angle,anticipated time and stabilizing driving behavior is investigated. Hopf bifurcation analysis mainly consists of two modules, one is theoretical analysis and the other is numerical simulation. Through the Hopf bifurcation analysis we can obtain the critical points.Results show that bifurcation analysis results are highly coincide with the stability condition of linear stability analysis.Furthermore,from the three-dimensional bifurcation diagram,it is found that the critical points of Hopf bifurcation are located on the boundary between the stability region and the instability region. Hopf bifurcation leads to conversion of stable traffic flow to irregular traffic flow. In addition, the TDGL equation is obtained by nonlinear analysis,and then the coexistence curve,metastable curve and critical point are described according to the obtained thermodynamic potential.We obtain the mKdV equation of the model in the same way. Finally,the parameter calibration and numerical simulation results display that the anticipated time and stabilizing driving behavior both play a positive role on traffic flow, which is opposite to the visual angle. Based on the conclusions, the article will make contributions on providing theoretical support for the establishment of intelligent transportation systems and references for traffic management systems.

    Acknowledgements

    Project supported by the Natural Science Foundation of Zhejiang Province, China (Grant Nos. LY22G010001,LY20G010004), the Program of Humanities and Social Science of Education Ministry of China (Grant No. 20YJA630008), the National Key Research and Development Program of China-Traffic Modeling, Surveillance and Control with Connected & Automated Vehicles (Grant No. 2017YFE9134700), and the K.C. Wong Magna Fund in Ningbo University,China.

    猜你喜歡
    榮軍紅霞
    請(qǐng)你幫個(gè)忙
    Traffic flow prediction based on BILSTM model and data denoising scheme
    An extended smart driver model considering electronic throttle angle changes with memory
    Stabilization strategy of a car-following model with multiple time delays of the drivers?
    《烏鴉喝水》中的“想”
    A Study of Combination of English Language Teaching and Context
    大東方(2018年1期)2018-05-30 01:27:23
    讓動(dòng)作“活”起來(lái)
    “大愛(ài)有聲 義工行動(dòng) 關(guān)愛(ài)榮譽(yù)軍人”慰問(wèn)活動(dòng) 走進(jìn)廣東省第一榮軍醫(yī)院
    大社會(huì)(2016年6期)2016-05-04 03:42:20
    “光的直線傳播”“光的反射”練習(xí)
    接力生命花季無(wú)聲
    ——記王臺(tái)中學(xué)初一158班郭榮軍同學(xué)舍己救人的英雄事跡
    91精品伊人久久大香线蕉| 国产国拍精品亚洲av在线观看| 一级,二级,三级黄色视频| av有码第一页| 老司机影院毛片| a级片在线免费高清观看视频| 国产乱来视频区| 美女国产视频在线观看| 大话2 男鬼变身卡| 精品99又大又爽又粗少妇毛片| 在线观看www视频免费| 熟女人妻精品中文字幕| 热re99久久国产66热| 大片免费播放器 马上看| 日日爽夜夜爽网站| 啦啦啦视频在线资源免费观看| 久久久久久久精品精品| 黄片播放在线免费| 欧美人与性动交α欧美精品济南到 | 免费av不卡在线播放| 日韩精品免费视频一区二区三区 | 五月玫瑰六月丁香| 亚洲精品视频女| 熟女av电影| 成年女人在线观看亚洲视频| 99热这里只有精品一区| 曰老女人黄片| 国产av国产精品国产| 国产精品国产三级国产av玫瑰| 一区二区三区精品91| 岛国毛片在线播放| 久久这里有精品视频免费| 国产精品国产三级国产专区5o| 午夜福利影视在线免费观看| 亚洲成人手机| 国产精品欧美亚洲77777| 日韩强制内射视频| 国产亚洲欧美精品永久| 亚洲怡红院男人天堂| 日本黄大片高清| 成人亚洲精品一区在线观看| 久久毛片免费看一区二区三区| 国产一区有黄有色的免费视频| 一本大道久久a久久精品| 视频中文字幕在线观看| 2018国产大陆天天弄谢| 美女主播在线视频| 在线看a的网站| 国产日韩欧美在线精品| 亚洲欧美中文字幕日韩二区| 亚洲av电影在线观看一区二区三区| 精品一区二区免费观看| 免费观看在线日韩| 不卡视频在线观看欧美| 大又大粗又爽又黄少妇毛片口| 妹子高潮喷水视频| 国产精品嫩草影院av在线观看| 黄色毛片三级朝国网站| 综合色丁香网| 青青草视频在线视频观看| 色视频在线一区二区三区| 成人漫画全彩无遮挡| 最新中文字幕久久久久| 欧美日韩精品成人综合77777| 亚洲精品日韩在线中文字幕| 在线看a的网站| 久久久久久久久久久久大奶| 肉色欧美久久久久久久蜜桃| 黄色一级大片看看| 国产午夜精品一二区理论片| 精品一区二区三卡| 亚洲人与动物交配视频| 一二三四中文在线观看免费高清| 亚洲精品国产色婷婷电影| 精品久久国产蜜桃| 亚洲婷婷狠狠爱综合网| 婷婷色av中文字幕| 在线观看人妻少妇| 飞空精品影院首页| 日本欧美国产在线视频| 成人18禁高潮啪啪吃奶动态图 | 三级国产精品欧美在线观看| 97在线人人人人妻| 一区二区三区免费毛片| 国产精品久久久久久精品古装| 亚洲精品aⅴ在线观看| 久久久久国产网址| 看十八女毛片水多多多| 国产伦理片在线播放av一区| 午夜福利视频在线观看免费| 久久99热6这里只有精品| 特大巨黑吊av在线直播| 午夜精品国产一区二区电影| 亚洲国产欧美日韩在线播放| 超碰97精品在线观看| 久久精品夜色国产| a 毛片基地| 国产成人aa在线观看| 在线 av 中文字幕| 国产成人a∨麻豆精品| 亚洲国产精品一区三区| 视频区图区小说| 免费人妻精品一区二区三区视频| 老司机影院毛片| 熟女av电影| 夜夜爽夜夜爽视频| 亚洲人与动物交配视频| 色视频在线一区二区三区| 午夜免费鲁丝| 9色porny在线观看| 人人妻人人添人人爽欧美一区卜| av在线播放精品| 精品久久久精品久久久| 国产成人a∨麻豆精品| 一级毛片我不卡| 欧美精品亚洲一区二区| 亚洲人成网站在线播| 亚洲精品美女久久av网站| 男的添女的下面高潮视频| 日韩大片免费观看网站| 狂野欧美激情性xxxx在线观看| 国语对白做爰xxxⅹ性视频网站| 91精品国产九色| 精品人妻在线不人妻| 成年人免费黄色播放视频| 在线天堂最新版资源| 91aial.com中文字幕在线观看| 亚洲精品自拍成人| 国产一区二区三区av在线| 免费高清在线观看视频在线观看| 亚洲精品自拍成人| 久久久午夜欧美精品| freevideosex欧美| 色94色欧美一区二区| 尾随美女入室| 又黄又爽又刺激的免费视频.| 亚洲图色成人| 国产色爽女视频免费观看| 啦啦啦啦在线视频资源| 蜜桃国产av成人99| 国产精品一二三区在线看| 美女cb高潮喷水在线观看| 久久韩国三级中文字幕| 国产精品三级大全| 日韩伦理黄色片| 3wmmmm亚洲av在线观看| 亚洲精品自拍成人| 日本-黄色视频高清免费观看| 亚洲中文av在线| 国产精品人妻久久久久久| 精品久久久噜噜| 欧美精品高潮呻吟av久久| 日本欧美视频一区| 极品少妇高潮喷水抽搐| 青春草国产在线视频| 亚洲欧美一区二区三区国产| 人妻人人澡人人爽人人| 天美传媒精品一区二区| 18禁裸乳无遮挡动漫免费视频| 国产精品成人在线| 日韩人妻高清精品专区| 不卡视频在线观看欧美| 女人久久www免费人成看片| 国产免费又黄又爽又色| 亚洲国产精品国产精品| 国模一区二区三区四区视频| 午夜激情久久久久久久| 国产一区二区在线观看av| 91午夜精品亚洲一区二区三区| 亚洲无线观看免费| 久久久欧美国产精品| 国产毛片在线视频| 熟女av电影| 不卡视频在线观看欧美| 国产成人午夜福利电影在线观看| tube8黄色片| 久久国产精品男人的天堂亚洲 | 精品午夜福利在线看| 91国产中文字幕| 欧美人与善性xxx| 精品亚洲成a人片在线观看| 亚洲av综合色区一区| 国产亚洲一区二区精品| 欧美精品亚洲一区二区| av不卡在线播放| 国产精品99久久99久久久不卡 | 啦啦啦在线观看免费高清www| 一级毛片电影观看| 精品久久久精品久久久| 热99国产精品久久久久久7| 亚洲国产精品成人久久小说| 亚洲美女搞黄在线观看| 精品亚洲乱码少妇综合久久| 久久久久久久国产电影| h视频一区二区三区| 国产精品99久久久久久久久| 满18在线观看网站| 日本wwww免费看| 中国三级夫妇交换| 日韩制服骚丝袜av| 男女国产视频网站| 国产精品熟女久久久久浪| 久久久a久久爽久久v久久| 欧美97在线视频| 精品熟女少妇av免费看| 国产黄色免费在线视频| 亚洲av.av天堂| 国产成人av激情在线播放 | 久久久精品区二区三区| 黄色欧美视频在线观看| 国产av精品麻豆| 日韩中字成人| 大片电影免费在线观看免费| 最新中文字幕久久久久| 婷婷色综合大香蕉| 亚洲国产av影院在线观看| 亚洲伊人久久精品综合| 欧美日韩视频精品一区| 搡女人真爽免费视频火全软件| 一本色道久久久久久精品综合| 久久久久久久久久久久大奶| 男女免费视频国产| 91国产中文字幕| 国产成人免费观看mmmm| 亚洲内射少妇av| 校园人妻丝袜中文字幕| 好男人视频免费观看在线| 国产 精品1| 麻豆乱淫一区二区| 精品亚洲成国产av| 亚洲无线观看免费| 天堂中文最新版在线下载| 亚洲精品第二区| 国产亚洲一区二区精品| 99久国产av精品国产电影| 久久久久精品性色| 超碰97精品在线观看| 亚洲国产欧美日韩在线播放| 亚洲精品中文字幕在线视频| 99久久人妻综合| 午夜福利视频在线观看免费| 久久女婷五月综合色啪小说| 母亲3免费完整高清在线观看 | 久久99一区二区三区| 汤姆久久久久久久影院中文字幕| 人人妻人人澡人人看| 国产欧美另类精品又又久久亚洲欧美| 一级二级三级毛片免费看| 春色校园在线视频观看| 人妻少妇偷人精品九色| 国产日韩欧美视频二区| 特大巨黑吊av在线直播| 亚洲欧美色中文字幕在线| 一级毛片aaaaaa免费看小| 国产毛片在线视频| 国产精品一区二区三区四区免费观看| 色94色欧美一区二区| av网站免费在线观看视频| 男的添女的下面高潮视频| 91aial.com中文字幕在线观看| 色吧在线观看| 少妇的逼好多水| 国产精品一区二区在线观看99| 寂寞人妻少妇视频99o| 亚洲av福利一区| 精品一区二区三区视频在线| 一区二区日韩欧美中文字幕 | 在线播放无遮挡| 亚洲欧洲国产日韩| 午夜福利影视在线免费观看| 久久精品久久久久久久性| 国产成人一区二区在线| 99精国产麻豆久久婷婷| 女人精品久久久久毛片| 婷婷色综合www| 日韩精品免费视频一区二区三区 | 丝袜喷水一区| 少妇高潮的动态图| 男人操女人黄网站| 卡戴珊不雅视频在线播放| 欧美日韩亚洲高清精品| 人人妻人人澡人人爽人人夜夜| 99热国产这里只有精品6| 亚洲婷婷狠狠爱综合网| 欧美成人午夜免费资源| 高清午夜精品一区二区三区| 人人妻人人澡人人看| 久久国产亚洲av麻豆专区| 国产老妇伦熟女老妇高清| 黄色欧美视频在线观看| 亚洲精品aⅴ在线观看| 国产欧美另类精品又又久久亚洲欧美| 卡戴珊不雅视频在线播放| 午夜激情福利司机影院| 日韩电影二区| 欧美日韩一区二区视频在线观看视频在线| 久久久久精品久久久久真实原创| 中文字幕av电影在线播放| 99热全是精品| av视频免费观看在线观看| 午夜福利网站1000一区二区三区| 欧美日韩视频精品一区| 国产日韩一区二区三区精品不卡 | 国产精品国产三级国产专区5o| videos熟女内射| 午夜福利,免费看| 91精品三级在线观看| av网站免费在线观看视频| 中国三级夫妇交换| 欧美性感艳星| 日本欧美国产在线视频| 18禁观看日本| 午夜视频国产福利| 你懂的网址亚洲精品在线观看| 精品人妻偷拍中文字幕| 亚洲欧美成人综合另类久久久| 亚洲内射少妇av| 日本欧美国产在线视频| 天堂俺去俺来也www色官网| videos熟女内射| 国产欧美亚洲国产| 日本av免费视频播放| xxx大片免费视频| 精品一区二区三区视频在线| av在线app专区| 中文字幕人妻丝袜制服| 搡女人真爽免费视频火全软件| 成人漫画全彩无遮挡| 在现免费观看毛片| 亚洲图色成人| 国产成人午夜福利电影在线观看| 天堂中文最新版在线下载| 成人手机av| 少妇的逼好多水| 黄色毛片三级朝国网站| 97精品久久久久久久久久精品| 精品亚洲成国产av| 色5月婷婷丁香| 纵有疾风起免费观看全集完整版| 看非洲黑人一级黄片| 伦理电影免费视频| a级毛片黄视频| 热re99久久精品国产66热6| 在线观看人妻少妇| 国产白丝娇喘喷水9色精品| 国产熟女午夜一区二区三区 | 蜜桃国产av成人99| 国产不卡av网站在线观看| 日韩av在线免费看完整版不卡| 日韩强制内射视频| 亚洲成人一二三区av| 一个人看视频在线观看www免费| 这个男人来自地球电影免费观看 | 最近最新中文字幕免费大全7| 最近手机中文字幕大全| tube8黄色片| 中文字幕人妻丝袜制服| 十分钟在线观看高清视频www| 久久久久久久精品精品| 建设人人有责人人尽责人人享有的| 一个人看视频在线观看www免费| 一区二区av电影网| 99九九在线精品视频| 秋霞伦理黄片| 如日韩欧美国产精品一区二区三区 | 观看av在线不卡| 一本久久精品| 在线观看免费视频网站a站| 蜜臀久久99精品久久宅男| 在线观看免费高清a一片| 十八禁高潮呻吟视频| 边亲边吃奶的免费视频| 国产视频首页在线观看| 免费日韩欧美在线观看| 美女中出高潮动态图| 亚洲精品日韩在线中文字幕| 亚洲色图综合在线观看| 日日爽夜夜爽网站| 久久精品久久久久久久性| 日韩熟女老妇一区二区性免费视频| 亚洲熟女精品中文字幕| 亚洲高清免费不卡视频| 国产成人午夜福利电影在线观看| 成人黄色视频免费在线看| 丰满少妇做爰视频| 赤兔流量卡办理| 桃花免费在线播放| 建设人人有责人人尽责人人享有的| 国产白丝娇喘喷水9色精品| 亚洲综合色惰| 亚洲精品日本国产第一区| 国产精品不卡视频一区二区| 菩萨蛮人人尽说江南好唐韦庄| 水蜜桃什么品种好| 亚洲人成77777在线视频| 青青草视频在线视频观看| 亚洲av在线观看美女高潮| 欧美精品一区二区大全| 亚洲精品久久午夜乱码| 免费黄频网站在线观看国产| 国产黄色免费在线视频| 日本爱情动作片www.在线观看| 国产亚洲精品第一综合不卡 | 欧美97在线视频| 这个男人来自地球电影免费观看 | 一二三四中文在线观看免费高清| 国产欧美日韩一区二区三区在线 | 久久青草综合色| 人妻系列 视频| videosex国产| 国产欧美另类精品又又久久亚洲欧美| 99热6这里只有精品| 蜜桃在线观看..| 久久99一区二区三区| 18禁在线播放成人免费| 在线亚洲精品国产二区图片欧美 | 99久久中文字幕三级久久日本| 少妇猛男粗大的猛烈进出视频| 久久久久久久国产电影| 你懂的网址亚洲精品在线观看| 精品一区在线观看国产| 久久久久网色| 日韩视频在线欧美| 美女脱内裤让男人舔精品视频| 国产精品国产三级国产av玫瑰| 亚洲四区av| 蜜桃在线观看..| 久久精品国产自在天天线| 在线观看美女被高潮喷水网站| 亚洲av综合色区一区| 国产高清不卡午夜福利| 日韩伦理黄色片| 国产免费视频播放在线视频| 欧美成人午夜免费资源| 国产熟女欧美一区二区| 成人国产av品久久久| 国产精品蜜桃在线观看| 人人妻人人爽人人添夜夜欢视频| 国产精品国产三级国产av玫瑰| 性色av一级| 国产免费一级a男人的天堂| 久久97久久精品| 考比视频在线观看| 极品人妻少妇av视频| 国产乱来视频区| 欧美激情国产日韩精品一区| 十分钟在线观看高清视频www| 国产精品一区www在线观看| www.av在线官网国产| 一本—道久久a久久精品蜜桃钙片| 亚洲国产色片| 99re6热这里在线精品视频| 一本色道久久久久久精品综合| 一本大道久久a久久精品| 久久 成人 亚洲| 午夜影院在线不卡| 国产极品天堂在线| 在线观看免费高清a一片| 天天影视国产精品| 久久精品久久久久久噜噜老黄| videossex国产| 日韩一区二区视频免费看| 男的添女的下面高潮视频| 久久影院123| 一个人免费看片子| 免费高清在线观看视频在线观看| 毛片一级片免费看久久久久| 中文字幕亚洲精品专区| 人成视频在线观看免费观看| 久久久久久久精品精品| 欧美97在线视频| 国产熟女午夜一区二区三区 | 久久亚洲国产成人精品v| 成人国产av品久久久| 国产有黄有色有爽视频| 成人18禁高潮啪啪吃奶动态图 | 精品酒店卫生间| 中文字幕最新亚洲高清| 欧美成人午夜免费资源| 久久久久久久久久久久大奶| 亚洲精品久久午夜乱码| 国产精品麻豆人妻色哟哟久久| 欧美日韩成人在线一区二区| 99久久精品国产国产毛片| 在线精品无人区一区二区三| 国产成人免费观看mmmm| 色吧在线观看| 国产69精品久久久久777片| 岛国毛片在线播放| 九九久久精品国产亚洲av麻豆| 亚洲成人手机| 亚洲精品一二三| 日韩大片免费观看网站| 国产69精品久久久久777片| 久久久久久久久久人人人人人人| 菩萨蛮人人尽说江南好唐韦庄| 日韩欧美精品免费久久| 欧美97在线视频| 国产熟女午夜一区二区三区 | 老女人水多毛片| 18在线观看网站| 观看美女的网站| 精品人妻一区二区三区麻豆| 久久热精品热| 中文字幕av电影在线播放| 少妇精品久久久久久久| 欧美日韩在线观看h| 一本大道久久a久久精品| 不卡视频在线观看欧美| 精品久久久久久电影网| 亚洲精品中文字幕在线视频| 亚洲精品aⅴ在线观看| 男女高潮啪啪啪动态图| 九九爱精品视频在线观看| 亚洲精品亚洲一区二区| 成年人午夜在线观看视频| 狠狠婷婷综合久久久久久88av| 免费大片18禁| 99久久精品一区二区三区| 日韩大片免费观看网站| 久久久久久久国产电影| 亚洲精品久久久久久婷婷小说| 国产精品一区www在线观看| 九九在线视频观看精品| 一区二区三区四区激情视频| 欧美变态另类bdsm刘玥| 国产男女内射视频| 国产日韩欧美在线精品| 十分钟在线观看高清视频www| 久久久久久久国产电影| 免费人妻精品一区二区三区视频| 日韩亚洲欧美综合| 色94色欧美一区二区| 青春草亚洲视频在线观看| 免费看av在线观看网站| 我的老师免费观看完整版| 热99国产精品久久久久久7| 中文字幕人妻丝袜制服| 国产在视频线精品| 欧美三级亚洲精品| 99久国产av精品国产电影| 国产成人精品无人区| av女优亚洲男人天堂| 国产免费福利视频在线观看| 春色校园在线视频观看| 99国产综合亚洲精品| 一二三四中文在线观看免费高清| 黄色一级大片看看| 黄色怎么调成土黄色| 尾随美女入室| 2018国产大陆天天弄谢| 亚洲av成人精品一区久久| 日本wwww免费看| 18禁动态无遮挡网站| 成人黄色视频免费在线看| 亚洲精品第二区| 日韩av不卡免费在线播放| 在线免费观看不下载黄p国产| 国产永久视频网站| 黄色视频在线播放观看不卡| 永久免费av网站大全| 日本欧美国产在线视频| 成年美女黄网站色视频大全免费 | 男女无遮挡免费网站观看| 大片电影免费在线观看免费| 午夜影院在线不卡| 中文精品一卡2卡3卡4更新| 国产色婷婷99| 在线观看三级黄色| 久久这里有精品视频免费| av一本久久久久| 亚洲精品久久成人aⅴ小说 | 亚洲国产精品成人久久小说| 在线看a的网站| 大香蕉97超碰在线| 欧美xxⅹ黑人| 国产老妇伦熟女老妇高清| 国产免费一级a男人的天堂| 水蜜桃什么品种好| 欧美日韩成人在线一区二区| 丝袜在线中文字幕| 99热国产这里只有精品6| 国产深夜福利视频在线观看| 久久久久久伊人网av| 久久久a久久爽久久v久久| 国产成人免费观看mmmm| 乱人伦中国视频| 久久精品国产自在天天线| 久久精品久久久久久久性| 久久99蜜桃精品久久| 国产乱来视频区| 18在线观看网站| 午夜激情av网站| 美女大奶头黄色视频| 久久久久网色| 寂寞人妻少妇视频99o| 人妻人人澡人人爽人人| 色网站视频免费| 成人无遮挡网站| .国产精品久久| 18禁在线无遮挡免费观看视频| 亚洲精品美女久久av网站| 久久影院123| 日韩成人伦理影院| 免费大片18禁| 欧美激情国产日韩精品一区| 日本爱情动作片www.在线观看| 伊人久久精品亚洲午夜| 校园人妻丝袜中文字幕| 看免费成人av毛片| 亚洲色图综合在线观看| 五月开心婷婷网| 精品卡一卡二卡四卡免费| 人人妻人人澡人人爽人人夜夜| 王馨瑶露胸无遮挡在线观看|