An improved car-following model considering variable safety headway distance. Yu-han Jia a,*, Jian-ping Wu a

Similar documents
ANOTHER PROOF FOR FERMAT S LAST THEOREM 1. INTRODUCTION

Nonparametric Goodness-of-Fit Tests for Discrete, Grouped or Censored Data 1

Fluids Lecture 2 Notes

Construction of Control Chart for Random Queue Length for (M / M / c): ( / FCFS) Queueing Model Using Skewness

I. Existence of photon

THE MEASUREMENT OF THE SPEED OF THE LIGHT

Observer Design with Reduced Measurement Information

ME260W Mid-Term Exam Instructor: Xinyu Huang Date: Mar

Fuzzy Dynamic Characteristic of Concrete. Material under Impact Loads

Summation Method for Some Special Series Exactly

Lecture 8. Dirac and Weierstrass

ε > 0 N N n N a n < ε. Now notice that a n = a n.

Local Estimates for the Koornwinder Jacobi-Type Polynomials

Interaction Mechanical Analysis between the Lunar Rover Wheel-Leg Foot and Lunar Soil

Certain inclusion properties of subclass of starlike and convex functions of positive order involving Hohlov operator

Optimal Management of the Spare Parts Stock at Their Regular Distribution

The beta density, Bayes, Laplace, and Pólya

Doppler frequency offset estimation and diversity reception scheme of high-speed railway with multiple antennas on separated carriage

One way Analysis of Variance (ANOVA)

Physics 3 (PHYF144) Chap 8: The Nature of Light and the Laws of Geometric Optics - 1

Lesson 4. Thermomechanical Measurements for Energy Systems (MENR) Measurements for Mechanical Systems and Production (MMER)

THE ANALYSIS OF GEOTHERMAL FIELD CHARACTERISTICS IN SICHUAN BASIN OF CHINA

Chapter 4: Angle Modulation

An improved attribute reduction algorithm based on mutual information with rough sets theory

Chapter 8 Hypothesis Testing

COMP26120: Introducing Complexity Analysis (2018/19) Lucas Cordeiro

Journal of Engineering Science and Technology Review 9 (6) (2016s) Research Article. Bai Shi-ye* and Zhu Jie-jiang

Basic Probability/Statistical Theory I

Nonstandard Lorentz-Einstein transformations

EconoQuantum ISSN: Universidad de Guadalajara México

WRITTEN ASSIGNMENT 1 ANSWER KEY

Traffic flow of adaptive cruise control vehicles on highway with an on-ramp

What is a Hypothesis? Hypothesis is a statement about a population parameter developed for the purpose of testing.

After the completion of this section the student. V.4.2. Power Series Solution. V.4.3. The Method of Frobenius. V.4.4. Taylor Series Solution

Analog Filter Synthesis

Bernoulli Numbers. n(n+1) = n(n+1)(2n+1) = n(n 1) 2

B. Maddah ENMG 622 ENMG /20/09

Computer Science 188 Artificial Intelligence. Introduction to Probability. Probability Ryan Waliany

Advanced Sensitivity Analysis of the Semi-Assignment Problem

Research on Fuzzy Clustering Image Segmentation Algorithm based on GA and Gray Histogram Baoyi Wang1, a, Long Kang1, b, Shaomin Zhang1, c

Effect of Magnetic Field on Marangoni Convection in Relatively Hotter or Cooler Liquid Layer

Principal Component Analysis. Nuno Vasconcelos ECE Department, UCSD

(Dependent or paired samples) Step (1): State the null and alternate hypotheses: Case1: One-tailed test (Right)

Calculus 2 TAYLOR SERIES CONVERGENCE AND TAYLOR REMAINDER

Optimization design in Wind Turbine Blade Based on Specific Wind Characteristics

Class #25 Wednesday, April 19, 2018

Digital Signal Processing. Homework 2 Solution. Due Monday 4 October Following the method on page 38, the difference equation

Société de Calcul Mathématique SA Mathematical Modelling Company, Corp.

Michelson's Repetition of the Fizeau Experiment:

The Use of L-Moments in the Peak Over Threshold Approach for Estimating Extreme Quantiles of Wind Velocity

Chapter 13, Part A Analysis of Variance and Experimental Design

= 47.5 ;! R. = 34.0 ; n air =

A first-order phase-transition, a super-cooled fluid, and a glass in a twodimensional. and A. Baram (2)

16th International Symposium on Ballistics San Francisco, CA, September 1996

LINEAR STABILITY ANALYSIS OF A PLANE-POISEUILLE HYDROMAGNETIC FLOW USING ADOMIAN DECOMPOSITION METHOD

Quasi Normal Modes description of transmission properties for Photonic Band Gap structures.

Sx [ ] = x must yield a

CONSIDERATION OF SOME PERFORMANCE OF CONTAINERS FLOW AT YARD

This article appeared in a journal published by Elsevier. The attached copy is furnished to the author for internal non-commercial research and

An Introduction to the Theory of Imprecise Soft Sets

Mass Transfer Chapter 3. Diffusion in Concentrated Solutions

Polynomial Functions and Their Graphs

Effects of Air Humidity on the Performance of a Polymer Insulator under Lightning Induced Voltage Conditions

FIR Filters. Lecture #7 Chapter 5. BME 310 Biomedical Computing - J.Schesser

INF5620 Key exercise 1

REVERSIBLE NON-FLOW PROCESS CONSTANT VOLUME PROCESS (ISOCHORIC PROCESS) In a constant volume process, he working substance is contained in a rigid

Stabilization of a Magnetic Levitation Control System via State-PI Feedback

Spectral Partitioning in the Planted Partition Model

Rational Bounds for the Logarithm Function with Applications

Confidence Interval for Standard Deviation of Normal Distribution with Known Coefficients of Variation

Solving the ZF Receiver Equation for MIMO Systems Under Variable Channel Conditions Using the Block Fourier Algorithm

Quasi Normal Modes description. of transmission properties. for Photonic Band Gap structures.

An Optimized Classification Model for Time-Interval Sequences

Abstract. Fermat's Last Theorem Proved on a Single Page. "The simplest solution is usually the best solution"---albert Einstein

Chapter 11 Output Analysis for a Single Model. Banks, Carson, Nelson & Nicol Discrete-Event System Simulation

ESTIMATION OF MACHINING ERRORS ON GLEASON BEVEL

NEW IDENTIFICATION AND CONTROL METHODS OF SINE-FUNCTION JULIA SETS

Goodness-Of-Fit For The Generalized Exponential Distribution. Abstract

Optimal design of N-Policy batch arrival queueing system with server's single vacation, setup time, second optional service and break down

Review of Grischuk and Sachin Gravitational Wave Generator via Tokamak Physics

SYNTHESIS OF SIGNAL USING THE EXPONENTIAL FOURIER SERIES

4. Optical Resonators

Fixed Point Approximation of Weakly Commuting Mappings in Banach Space

COM BIN A TOR I A L TOURNAMENTS THAT ADMIT EXACTLY ONE HAMILTONIAN CIRCUIT

Finite element analysis for delamination of laminated vibration damping steel sheet

Study on Coal Consumption Curve Fitting of the Thermal Power Based on Genetic Algorithm

Statistical Fundamentals and Control Charts

Closed-form two-step weighted-leastsquares-based. localisation using invariance property of maximum likelihood estimator in multiplesample

ENGI 4421 Confidence Intervals (Two Samples) Page 12-01

THE APPEARANCE OF FIBONACCI AND LUCAS NUMBERS IN THE SIMULATION OF ELECTRICAL POWER LINES SUPPLIED BY TWO SIDES

Double Stage Shrinkage Estimator of Two Parameters. Generalized Exponential Distribution

Comparison Study of Series Approximation. and Convergence between Chebyshev. and Legendre Series

Microscopic traffic flow modeling

MULTIPLE TRAFFIC JAMS IN FULL VELOCITY DIFFERENCE MODEL WITH REACTION-TIME DELAY

STK4011 and STK9011 Autumn 2016

Nonlinear regression

A Question. Output Analysis. Example. What Are We Doing Wrong? Result from throwing a die. Let X be the random variable

PAijpam.eu ON TENSOR PRODUCT DECOMPOSITION

Dr R Tiwari, Associate Professor, Dept. of Mechanical Engg., IIT Guwahati,

Complex Algorithms for Lattice Adaptive IIR Notch Filter

Transcription:

A improved ar-followig model osiderig variable safety headway distae Yu-ha Jia a,*, Jia-pig Wu a ( a Departmet of Civil Egieerig, Tsighua Uiversity, Beijig 00084, Chia) Abstrat: Cosiderig high speed followig o expressway or highway, a improved ar-followig model is develed i this paper by itroduig variable safety headway distae. Stability aalysis of the ew model is arried out usig the otrol theory method. Fially, umerial simulatios are implemeted ad the results show good osistey with theoretial study. Re śume : A l égard du phéomèe de ar-followig, est-à-dire la suite de voiture, àgrade vitesse sur l autoroute, et artile présete u modèle ar-followig basésur la distae de séuritémiimum de variable. Esuite, l aalyse de la stabilité du ouveau modèle e utilisat la théorie de la ommade est effetué, suivie par u grad ombre d expérimetatios simulées mises e oeuvre. Efi, les résultats motret que l étude théorique du ouveau modèle est orrete. Key words: Car-followig model; Variable safety headway distae; Stability aalysis. * Correspodig author: Yu-ha Jia E-mail address: jiajiayuyuhaha@63.om

. Itrodutio Car-Followig theory is a method to desribe how vehiles follow oe aother o a roadway. I the model ars are regarded as disrete ad iteratig partiles without overtakig to aalyze traffi flow harateristis. So far, a variety of models have bee develed, iludig safety distae model, stimulus-respose model ad fuzzy logi based model -8. ew ar-followig model osiderig variable safety headway distae (VSHD). I setio the ew model is develed ad stability aalysis is arried out i setio 3. The umerial simulatio experimets are performed to verify the theoretial study i setio 4. The summary is give i setio 5.. Improved model The typial OVM is preseted as Amog those models, the timal veloity model d x t () [ V ( x( t)) v( t)] dt () (OVM) is well kow for auray ad ratioality 3,4,9. Afterwards the OVM has bee exteded by itroduig geeralized fore, full veloity differee ad lateral effets 0-4. However, the models metioed aot be used to preisely reflet the ar-followig pheomeo i expressway or highway with higher speed ad larger headway (mostly loger tha 0m), i whih ase the hyperboli taget part ( tah( ) ) i where x ad v are the positio ad veloity of the th vehile; x is the headway distae betwee the th ad its leadig vehile; is the sesitivity parameter of the driver; V () is the OV futio desribed as where roadway; vmax V ( x ( t)) [tah( x ( t) h ) tah( h )] () v max is the maximum veloity o a partiular h meas the safety headway distae. OVM remais as a ostat. Therefore the safety However, as otied i the study o expressway headway distae parameter i OVM eeds to be modified. where vehiles have higher speed ad larger headway, the x h part is larger tha o regular road Based o previous work, this paper ivestigates a whih makes the tah( ) remai as. Hee the fol-

lowig vehile gets little or o ifluee despite x () t is hagig, whih otradits the real traffi situatio. dv [ Vew ( x ( t), v ( t)) v ( t)] dt v dy v ( t) v( t) dt, (4) To avoid the metioed problem, the fixed parameter h i OV futio is replaed by VSHD as where y ( t) x ( t) x ( t ). We assume the first vehile is ot iflueed by others ad rus ostatly at hf bv ts h, where t s is the time step uit. The speed v 0, the the steady state is give by physial meaig of h f is that the aeptable safety headway distae of a driver is dyami hagig T ( ), ( ), ( 0) T ew v t y t v V v. (5) based o v level. The liearized system of (9) a be alulated Referrig to previous study 4,5, the veloity aroud steady state (0) as differee betwee the th ad its leadig vehile v () t is itrodued, makig the ew futio as d x [ V ( ( ), ( )) ( )] ew x t v t v t dt v v V x t v t x t h tah( hf )] hf bv ts h max ew ( ( ), ( )) [tah( ( ) f ) 3. Stability aalysis. (3) Aordig to stability aalysis method 5-7, the stable oditio of the modified ar-followig model a be show as d v y ()Λ t v tλ v t dt, (6) ( v t vt ) d y v t vt dt where v t v t v, () partial derivatives ad Vew ( x ( t), v ( t)) Λ v y t y t V v, Vew ( x ( t), v ( t)) Λ x ( t) y tv ew ( v0 ) After Laplae trasformatio, we have s Λ. y tv ew ( v0 ) V s V s Y s ps () s Λ, (7) where V s L( v t ), Y s L( y t ), L () is, the Laplae trasformatio ad the harateristi polyomial is p s s s s Λ s Λ. The 3

the trasfer futio a be obtaied as The suffiiet oditio a be obtaied as s Λ Gs 0 ps s Λ. (8) s Λ s Λ p s s s s s Λ Λ Based o stability theory, the traffi jam will Λ Λ Λ Λ 0 0, whih a be rewritte as, (9) ever happe i this system if G s. ps is stable ad Λ Λ Λ Λ 0. (0) The eutral stability surfae i the spae of As d x () t dt is i positive orrelatio with (Λ,Λ, ) with differet b are show i Fig.. [ V ( x ( t), v ( t)) v ( t)], we have 0. Beause of ew the OV futio harateristi, we have Λ >0, so the oditio for p s to be stable is Λ. The we osider expressed as G( j) G( j) G( j) G s whih a be Λ. Λ ( Λ ) Whe b=0, the model equals the full veloity differee model (FVDM) 4, ad the figure shows a eutral stability lie i Fig.(a). Aordig to Fig.(b) - (d), with higher b, the peak of stability surfae moves to larger x whe v ireases. As a result, to avoid traffi jam at higher speed, driver eeds to be more arefully to keep above the stability surfae. Fig. Neutral stability surfae with differet b: (a) b=0 (FVDM); (b) b=0.; () b=0.3; (d) b=0.5. 4

4. Numerial simulatios Let us osider a 00-vehile system ruig o a expressway without overtakig. I the simulatios, the iitial values are set as v (0) =5m/s, v max h =7m, = =0.5. =0m/s, aeleratio or deeleratio to keep the traffi system i stable situatio. 5. Summary I this paper, we prosed a ew ar-followig Firstly, we assume i the simulatio proess the leadig vehile has a radom flutuatio i aeleratio or deeleratio. Fig.(a) shows the spae-time plot whe b =0 (FVDM) ad we a observe the osillatig headway distaes i the lower vehiles. Fig.(b) represets the ew model with b =0.3, idiatig traffi jam ever happes. model osiderig the variable safety headway distae based o the typial OVM. The stability oditio of ew model is aalyzed by applyig the otrol theory. Fially, umerial simulatios are give, ad the results are osistet with the theoretial study. I olusio, the ew model suppresses the traffi jam i typial OVM with high speed ad larger headway Seodly, we osider the situatio with disturbae where the leadig vehile sts suddely distae. from t=0s to t=5s. The spae-time plot ad veloity behavior of three vehiles are desribed i Fig.3-5. We a see heavy traffi jams whe b =0, ad as b ireases the traffi jams are suppressed. Simulatios illustrate that with a larger headway distae but higher speed, vehiles i ew model are Akowledgemet The researh i this paper was oduted as part of the projet The ourree ad evolutio of traffi gridlok i mega-ity uder storm rai oditios, whih is fuded by Beijig Natural Siee Foudatio (Projet No. 9300). more sesitive to the variatio of leadig vehile s 5

(a) (b) Fig. Spae-time plot of: (a) Typial model (b=0); (b) New model (b=0.3). (a) (b) Fig.3 (a) Spae-time plot; (b) Veloity of the first, 5 th ad 50 th vehiles (b=0). (a) (b) Fig.4 (a) Spae-time plot; (b) Veloity of the first, 5 th ad 50 th vehiles (b=0.05). 6

(a) (b) Fig.5 (a) Spae-time plot; (b) Veloity of the first, 5 th ad 50 th vehiles (b=0.3). 7

[9] T.Nagatai, Phys. Rev. E. 58, 47 (998). Referees: [0] H.X. Ge, X.P. Meg, J. Ma, ad S.M. Lo, Phys. Lett. A. [] T.Q. Tag, J.G. Li, H.J. Huag, ad X.B. Yag, Meas. 48, 377, 9 (0). 63 (04). [] T. Nagatai, Phys. Rev. E. 6, 3564 (000). [] M. Brakstoe ad M. MDoald, Trasp Res F., 8 [] L. C. Davis, Phys. A. 39, 457 (003). (999). [3] D. Helbig ad B. Tilh, Phys. Rev. E. 8, 33 (998). [3] M. Bado, K. Hasebe, A. Nakayama, A. Shibata, ad Y. [4] R. Jiag, Q.S. Wu, ad Z.J. Zhu, Phys. Rev. E. 64, 070 Sugiyama, Phys. Rev. E. 5, 035 (995). (00). [4] M. Bado, K. Hasebe, A. Nakayama, A. Shibata, ad Y. [5] K. Koishi, H. Kokame, ad K. Hirata, Eur. Phys. J. B. 5, Sugiyama, Phys. Rev. E. 58, 549 (998). 75 (000). [5] H.X. Ge, S.Q. Dai, L.Y. Dog, ad Y. Xue, Phys. Rev. E. [6] K. Koishi, H. Kokame, ad K. Hirata, Phys. Rev. E. 60, 70,06634 (004). 4000 (999). [6] M. Papageorgiou, Trasp. Res. A. 3, 33 (998). [7] D. Yag, P. Ji, P. Yu, ad B. Ra, Phys. A. 395, 37 [7] M. Brakstoe, B. Sulta, ad M. MDoald, Trasp. Res. (04). F. 5, 3 (00). [8] J.P. Wu, M. Brakstoe, ad M. MDoald, Trasp Res C.,463 (003). 8

List of Figures Captios Fig. Neutral stability surfae with differet b: (a) b=0 (FVDM); (b) b=0.; () b=0.3; (d) b=0.5. Fig. Spae-time plot of (a) Typial model (b=0); (b) New model (b=0.3) Fig.3 (a) Spae-time plot (b=0); (b) Veloity of the first, 5 th ad 50 th vehiles (b=0) Fig.4 (a) Spae-time plot (b=0.05); (b) Veloity of the first, 5 th ad 50 th vehiles (b=0.05) Fig.5 (a) Spae-time plot (b=0.3); (b) Veloity of the first, 5 th ad 50 th vehiles (b=0.3) 9