Vehicular Ad-Hoc Using Matlab
Vehicular Ad-Hoc Using Matlab
s
3- The segment number is the headway rounded to +.
Fig. 4-5 and Fig. 4-6 show the difference between headway-based and distance-based
segmentation with vehicles running at different speeds.
Headway (sec)
Headway (sec)
49
Fig. 4-5. Distance-based segmentation
Fig. 4-6. Headway-based segmentation
Fig. 4-5 shows a 3-lane highway with three vehicles running at different speeds,
(30,60,120 Km/h) with reference to distance (meter). This figure is a real snap-shot image.
Fig. 4-6 shows the same situation after calculating the headway for each vehicle to produce an
imaginary calculated image. This image reveals that headway-based segmentation mimics
dangerous situations better than distance-based, that it puts the 120Km/h-vehicle into the first
segment, which is identical to the intuitive analysis of the situation.
From now on, figures of the highway will be of two types: the first one is figures with
reference to distance (meter) (e.g Fig. 4-5) which is a real snap-shot image. The second one is
that with reference to headway (sec) (e.g. Fig. 4-6) which is an imaginary calculated image
based on the location and speed of each vehicle.
Headway (sec)
CW [0-3] [4-7] [8-11] [9-12]
Segment S1 S2 S3 S4
Distance (meter)
120 Km/h
30 Km/h
CW [0-3] [4-7] [8-11] [9-12]
Segment S1 S2 S3 S4
50
4.4.1 Discussion
- Assuming multiple lanes highway:
A question that may arise if the analysis is for multiple lanes is do vehicles in different
lanes are prone to the same danger as vehicles running in the same lane following the
abnormal one? Studies found that some drivers avoid obstacles by steering rather than by
braking or even perform the both [14]. It is found that the response time for steering is about
0.3 sec faster than that for breaking. This is the cause that we should consider that the
abnormal vehicle may often use steering in conjunction with breaking and that the danger area
is not only the same lane of the abnormal vehicle, but also adjacent lanes.
- Assuming a single lane highway:
A question that may arise if the analysis is for a single lane is will a far fast vehicle
overtakes a near slow one?
Let us study this problem quantitatively; assume that there are two vehicles following the
abnormal one in a single lane highway as indicated in the figure.
Fig. 4-7. Assuming a single lane highway
The first vehicle is running in a speed S1, and the distance to the abnormal vehicle is d1
meters and H1 secs, where
E1 =
J1
S1
The second vehicle is running in a speed S2, and the distance to the first vehicle is d2 meters
and H2 secs, where
H1 & d1 H2 & d2
Ht & dt
Headway (sec)
S1 S2
51
E2 =
J2
S2
and the distance to the abnormal vehicle is dt meters and Ht secs, where
Et =
Jt
S2
The condition stated in the question is
Et < E1
J1 +J2
S2
<
J1
S1
J1
S2
+E2 <
J1
S1
J1 _
1
S1
-
1
S2
] > E2
J1 _
S2 -S1
S1 - S2
] > E2
J1 > E2_
S1 - S2
S2 -S1
]
E1 =
d1
S1
E1 > E2[
S2
S2-S1
For example, assume that S1 = 80km/h (22.2m/s), S2 = 120km/h (33.3m/s) and H2 is the
average minimum headway (1 sec) [41]
E1 > S sccs
Thus, if the two vehicles are running at 80, 120km/h (in developed countries, such a speed
difference is not expected to happen in the same highway lane), this situation may happen at a
minimum of three seconds away from the abnormal vehicle; i.e. a relaxed situation. Then,
either the fast vehicle will slow down (not to hit the slower one) or try to pass it and truly
become more threatened than the slower one and the one that should logically reply with the
CTB (or ACK) message.
For short, the stated situation may happen only if
- The speed difference is very large
- and the two vehicles are still far away from the abnormal vehicle
52
With step-2 modification, the algorithm elects the nearest-in-time vehicle with a plain
uniform headway-based segmentation method.
In this step, we emphasized that any VANET protocol should expect that vehicles are
running at different speeds and it should cope with these variations efficiently.
4.5 Step-3: Non-uniform Segmentation (Headway Model)
What modification that can be applied to the proposed protocol if we highlighted that, "the
distribution of vehicle-headway along highways has never been uniform"?
The analysis of the distribution of vehicle-headways is classified under the Headway
Model. The concept of Headway Model will be discussed in Sec 4.5.1, before continuing
with its impact on the protocol development with Sec 4.5.2.
4.5.1 Headway Model
The Headway Model is a mathematical equation that describes the average naturalistic
headway that drivers tend to leave apart. This model is fundamental in any traffic engineering
application because it provides a laboratorial method of generating vehicles in any traffic flow
simulator.
Traffic engineering researchers introduced many headway models trying to mimic realistic
situations. Some of these models are: the negative exponential distribution, the shifted
exponential distribution, the gamma distribution, the lognormal distribution and the semi-
poisson distribution.
Luttinen (1996) [32] ,in his PhD thesis, gave a good study comparing each model with
empirical data. Fig. 4-8 shows the comparison between the probability density function (pdf)
of these distributions and the empirical density function (edf) with reference to 10-sec
headway. The headway distribution shows that drives cannot hold a zero sec headway and the
average headway expected is about 1.5 sec. The headway model of Fig. 4-8 is measured at a
traffic density of 1439 vehicle/hour.
He stated that the model that gives the best goodness-of-fit results when compared with
the edf is the Semi-Poisson distribution, and that it is the recommended distribution for use in
53
simulators with adequate computational facilities. The Semi-Poisson distribution is the one
used throughout this work.
Fig. 4-8. Sample Headway models
Studying the headway model shows that, broadcasting protocols that depend on the time
headway (only this contribution at this time) have a high robustness at different traffic
volumes. This claim is clear comparing the headway model at different traffic volumes (Fig.
4-9 [32]) and with the simulation in Sec 5.5.
Fig. 4-9. Headway at different traffic volumes
a. 1439 veh/h b. 795 veh/h c. 328 veh/h
d. Lognormal Distribution e. Semi-Poisson
a. Negative Exponential b. Shifted Exponential c. Gamma Distribution
54
Although nothing has been changed in the developing protocol, this section satisfies a new
goal, namely High robustness at different traffic volumes.
4.5.1.1 The Semi-Poisson Distribution
The probability density function (pdf) of the Semi-Poisson distribution is recalled here;
t t
e
t
p e
t
p p t f
) 1 (
) (
) , (
) 1 (
) (
) (
) , , , | (
1
1 0 ; 0 , , ; 0 > p t Eq. 4-1
Where, is the Gamma Function and is the Incomplete Gamma Function. Parameters p
and changes according to the desired traffic volume; however, parameters and do
not correlate with the traffic volume. Parameters chosen in (Fig. 4-10) are considered at a
traffic volume of approximately 330 veh/h [32].
Fig. 4-10. Semi-Poisson Headway Model
Notwithstanding that there is a common verbal instruction used in driving manuals around
the world to maintain a minimum headway of 2 seconds [9], Authors in [41] found that the
0 2 4 6 8 10
0
0.05
0.1
0.15
0.2
0.25
0.3
0.35
Headway (sec)
D
e
n
s
i
t
y
Semi-Poisson Distribution at 330 veh/h
P = 0.14
[ = 2.929
= 5.488
0 = 0.06
55
naturalistic average minimum headway for drivers is about 0.66 second as indicated in Fig.
4-10, and it is insensitive to the driving speeds. This 0.66 sec is a very strict timing taking into
account that the best driver response time is about 0.70 to 0.75 sec [14].
4.5.2 Protocol Improvement
The headway model can dramatically change the segmentation algorithm. It can be a basis
for a non-uniform segmentation where the width of each segment is adapted to give any
required distribution of collision probability where a minimum probability of collision leads
directly to a minimum latency (Fig. 4-11).
Fig. 4-11. Non-uniform headway-based segmentation
Without loss of generality, assume that there are only 2 vehicles in the transmission range
of the abnormal vehicle. The headway between the abnormal vehicle and the first following
one is
1
X sec, and the headway between the two following vehicles is
2
X sec as shown in
Fig. 4-12. It is clear that both and are random variables with a Semi-Poisson
probability distribution function identical to Fig. 4-10. According to the probability theory,
both variables can be replaced with a common variable X , having the same pdf.
Fig. 4-12. Study area of the analytical solution
1
X
2
X
i
l
f
l
Study area
Headway (sec)
1
X
2
X
Headway (sec)
Segment S1 S2 S3 S4
56
Assuming that, the highway is only one lane and the CW is static ( = 1 ) i.e. there is no
contention or backoff, all vehicles within one segment will wait for exactly one TS before
initiating the CTB transmission i.e. for a successful broadcast, segments should include
exactly one vehicle.
To study the collision probability in one of the segments, we assume that the segment
under consideration is in-between any arbitrary headways l
and l
]
sec (as shown in Fig.
4-12). Although the headway is a contentious variable and the probabilities should have been
in the integral form, we list here the discrete forms (i.e. after discritization) to be identical to
the ones used in the optimization program that will be discussed in Sec 4.5.3.
There will be a collision in the CTB message only if there were more than one vehicle in
the segment, that the probability of collision (P
C
) is the probability that X
1
lies inbetween
l
& l
]
and X
2
is smaller than ( l
]
-X
1
)
=
< = =
f
i
l
l x
f C
x l X P x X P P ) ( ) ( Eq. 4-2
For completeness and check of correctness, probabilities of the other conditions that may
happen in the study area are:
The probability of success broadcast (P
B
) (i.e. only one node in the segment) is the probability
that X
1
lies inbetween l
& l
]
and X
2
is larger than ( l
]
-X
1
)
=
> = =
f
i
l
l x
f B
x l X P x X P P ) ( ) ( Eq. 4-3
The probability of idle (P
I
) (i.e. there is no nodes in the segment) is the probability that X
1
is
larger than l
]
.
) (
f I
l X P P > = Eq. 4-4
The probability (P
0
) (one of the prior nodes has captured the broadcast phase) is the
probability that X
1
is smaller than l
.
) (
0 i
l X P P < = Eq. 4-5
As a check of correctness, we note that, these probabilities sum to one as shown in Fig. 4-13.
57
Fig. 4-13. Probabilities associated with an arbitrary segment
In the analysis, we are only concerned with the probability of collision and latency (i.e. P
C
and P
B
). Hence, P
C
and P
B
are rewritten below with a normalization factor so that their
summation equals one i.e. the P
C
and P
B
given that the first vehicle lies within the segment.
=
<
< <
=
=
f
i
l
l x
f
f i
C
x l X P
l X l P
x X P
P ) (
) (
) (
Eq. 4-6
=
>
< <
=
=
f
i
l
l x
f
f i
B
x l X P
l X l P
x X P
P ) (
) (
) (
Eq. 4-7
Eq. 4-6 is the one that will be used in the analytical analysis of the proposed protocol.
Although we assumed only two vehicles for simplification, the stated equations are
applicable to any number of following vehicles because only the fastest two ones are the main
cause of collision. Following the same procedures, the analysis is expandable to as multiple
lanes as practically needed.
0 2 4 6 8 10
0
0.05
0.
0.15
0.
0.25
0.
0.35
Headway (sec)
D
e
n
s
i
t
y
pdf of
1
X
I
P
0
P
B C
P P +
58
Now, we can state the design problem as; within a single communication range (10 sec as
recommended by DSRC [8]), find the best points of segmentation that result in P
C
s that are
linearly increasing with a minimum slop as shown in Fig. 4-14.
Fig. 4-14. Suggested Distribution of Collisions
There are two reasons behind minimizing the slop instead of the absolute minimum:
- Intuitively, vehicles in the first segments are more threatened to danger than those in the last
segments. Hence, they should take a higher weight in the minimization process. Each vehicle
is exposed to a danger that is inversely proportional to the time before collision (a vehicle
collision) which is the same as the headway time. Fortunately, segmentation points are in
headway time and we will get a linearly increasing latency through a linearly increasing
C
P
.
- The other reason is a traffic concept that if there were no vehicles in the first segment, we
can expect that the traffic is moderate or low, and let later segments be of wider width.
0 2 4 6 8 10 12
Headway (sec)
P
r
o
b
a
b
i
l
i
t
y
o
f
C
o
l
l
i
s
i
o
n
s
(
r
e
q
u
i
r
e
d
)
The required probability of collision
Suggested
distribution
of P
C
Find the
minimum
59
4.5.3 Analytical Results
Using Matlab commercial program, a simple program was developed to compute the Semi-
Poisson equation (Sec 4.5.1.1) and to perform the minimization process as indicated above for
different number of segments ranging from (4) to (10) segments. The best points of
segmentation are listed in Table 4-1.
Table 4-1 Best segmentation points for 330 vehicle /h (in headway sec)
10-seg = [1.01 1.74 2.51 3.34 4.26 5.27 6.36 7.52 8.75 10];
9-seg = [1.12 1.93 2.79 3.73 4.79 5.96 7.23 8.58 10];
8-seg = [1.25 2.15 3.12 4.22 5.47 6.85 8.35 10];
7-seg = [1.42 2.45 3.59 4.93 6.47 8.18 10];
6-seg = [1.62 2.81 4.18 5.86 7.81 10];
5-seg = [1.86 3.27 5.03 7.28 10];
4-seg = [2.15 3.88 6.38 10];
We note that the width of each segment is monotonically increasing as indicated earlier;
with an upper bound to 10 sec. For example, for (6) segments, the width of each segment is
{1.19 - 1.37 - 1.68 - 1.95 - 2.19}. The width of the first segment is {1.62 sec} as an indication
of the effect of the minimum headway (as stated in 4.5.1.1).
Using the above results, the probability of collision are computed as in Eq. 4-6 and the
results are as shown in Fig. 4-15.
There is a note in this figure and all subsequent figures; values (here, P
C
) are measured as
an average on the segment, and are plotted on the last point of it. For example, the value of P
C
of the last segment of 4-seg case is an average over the range {6.38 10} and is plotted at
headway of 10 sec.
It is clear that the program did its job correctly, and it indicates that increasing the number
of segments leads directly to decreasing the overall average probability of collisions.
60
Fig. 4-15. Analytical calculation of Pc for best segmentation
4.6 Step-4: Application Adaptive (Modes of Operation)
Although the majority of VANET applications require message broadcasting, each
application has its unique flavor and requires a special treatment. The difference between
them is which of the following vehicles should have the highest priority to respond first,
either with a reply to the source vehicle or a relay to the following vehicles. This modification
enables changing the order of priority according to the application requirement. The RTB
message should contain a field for the Mode which will inform other vehicles of the category
of the broadcasting application so that the following nodes arrange priority accordingly.
Without loss of generality, we propose only four modes covering major applications.
4.6.1 Mode 0 Basic Broadcasting
The zero mode is the original basic mode, where broadcasting is omnidirectional with no
intended vehicle nor acknowledgment. This mode is still useful in VANET environment
especially in case of the status message, where, as recommended by DSRC [8], every
0 2 4 6 8 10 12
0
0.05
0.1
0.15
0.2
0.25
0.3
0.35
Headway (sec)
P
r
o
b
a
b
i
l
i
t
y
o
f
C
o
l
l
i
s
i
o
n
Pc for Non-uniform seg at 330 v/h (analytical)
10-seg
9-seg
8-seg
7-seg
6-seg
5-seg
4-seg
61
vehicle should broadcast its position, speed, direction of travel, and acceleration every
300 ms, and this transmission is intended for all vehicles within 10-sec travel time.
Fig. 4-16. Mode 0 Basic Broadcasting
4.6.2 Mode 1 The Furthest Following Vehicle
The intended vehicle in this mode is the furthest one following the transmitting vehicle
(Fig. 4-17). This mode is suitable to be used with dissemination protocol for applications like
Traffic Information, and Work Zone Warning. With these applications, the broadcast is
required to be delivered to the physically furthest node; that we recommend the regular
uniform distance-based protocols (e.g. The Smart Broadcasting Protocol) to be used in this
mode. Acknowledgment is recommended to be with a basic ACK packet.
Fig. 4-17. Priority arrangement of mode 1
4.6.3 Mode 2 The Nearest-in-time Following Vehicle
The intended vehicle, in this mode, is the nearest-in-time one running behind the
transmitting vehicle (Fig. 4-18). This mode is suitable to be used with reliable protocols for
Distance (meter)
Segment S4 S3 S2 S1
62
all public-safety related applications like Cooperative Collision Warning and Stop Light
Assistant. As indicated above, the proposed non-uniform headway-based protocol is the
superior in this mode. Acknowledgment is recommended to be with the same message setting
the ACK flag. This same message ACK is to compensate collisions at far range nodes due to
hidden node problem (Sec 4.3.1).
Fig. 4-18. Priority arrangement of mode 2
4.6.4 Mode 3 The Furthest Leading Vehicle
The intended vehicle, in this mode, is the furthest one leading the transmitting vehicle.
This mode is suitable for emergency applications like Approaching Emergency Vehicle
either it was an ambulance or a police car. In this case, the headway is identical to distance
because the speed is constant (headway is measure with reference to the speed of the vehicle
that comes later). However, with headway-based protocols, we can implement a non-uniform
segmentation based on headway studies. Acknowledgment is recommended to be with a basic
ACK packet.
Fig. 4-19. Priority arrangement of mode 3
Headway (sec)
Segment S1 S2 S3 S4
Headway (sec)
Segment S1 S2 S3 S4
63
Although these four modes cover major application categories, any other mode can be
added according to the application requirements.
4.7 The Proposed Algorithm
Although we already presented every modification and its impact on performance, this
section is devoted for implementation aspects.
The suggested Wave Short Message (WSM) frame format is the same as that indicated in
Fig. 2-9 with a new 8-bit field for Mode and ACK. This field can be a part of the Provider
Service Identifier, but as indicated in Fig. 4-20, it is recommended to be a part of the WSM
data, not to confuse with the P1609.3 standard [20].
1 1 1 1 1 4 2 1 variable
WSM
Version
Security
Type
Channel
Number
Data
Rate
Tx
Power
Level
Provider
Service
Identifier
WSM
Length
Mode
&
ACK
WSM
Data
Fig. 4-20. The suggested WSM frame format
4.7.1 Algorithm of the Transmitting node
Actions of the transmitting vehicle:
In case of an OBU has a message to broadcast, the MAC layer of the system has to proceed
with the following (Fig. 4-21),
1- It sends an RTB message including its MAC address, current location, current speed,
message propagation direction and the mode of operation.
2- It waits for a valid CTB message within SIFS+N+1 time-slots (assuming N segments).
If it received a valid CTB, then it should send the unencrypted broadcast to with intended
receiver as that indicated in the CTB message. Otherwise (if not), repeat the procedure from
the beginning (as long as the application requires).
64
Fig. 4-21. Actions of the transmitting MAC
4.7.2 Algorithm of Other Vehicles
Actions of the other vehicles/nodes MAC:
Upon receiving of an RTB message, other nodes proceed with the following algorithm
(Fig. 4-22),
1- Set the NAV to be SIFS+N+2 time-slots so that nodes will not start a new session until
the end of the current broadcast.
2- Check the broadcasting mode field.
3- Compare the geographical coordinates of the transmitting vehicle with their own, and
obtain its relative position. If the receiving vehicle is in the opposite driving direction or not
in the message propagation direction, ignore the message and go to end. Otherwise, if the
receiving vehicle is in the message propagation direction, continue to Step 4.
4- Compute the headway in seconds (or distance in meter for mode 2) then determines its
segment number with reference to the operating mode. Widths of each segment are
implemented according to Table 4-1.
send RTB
received CTB
within
SIFS+N+1 Ts?
again?
start
end
yes
no
no
yes
continue
session
65
Fig. 4-22. Actions of other vehicles
5- Assuming that the segment number equals S
i
where ( i <= N ) and ( i ) is the segment
number. Set the back-off counter to be equal to ( i-1 ).
So that, nodes wait for the SIFS then decrement the back-off counter by one in each time slot
while listening to the channel for any valid CTB message, if locked with a valid CTB
message then the node should exit the contention phase and listen for the incoming broadcast.
Set NAV to
SIFS+N+2 TS
start
end
yes
Is an affected
Vehicle?
no
Find segment
number Si'
Send CTB
Received CTB
within
SIFS+ i-1 Ts?
no
yes
Check mode
Continue session
66
The node that reaches zero initiates a CTB message including its MAC address and continues
the session with the transmitting node.
This completes the analytical analysis of the proposed protocol. The next chapter will
further investigate the protocol with some simulation programs.
67
Chapter 5
Simulation Results
The only possible method of evaluating new VANET protocols is by simulation (Sec 2.4).
In order to study the performance of the proposed protocol, it is required to implement two
new models;
- A mobility model that supports generating of vehicles with positions that change
according to random speeds, keeping the headway difference as a random with the
Semi-Poisson distribution.
- A new network model where decisions are taken based on the current location, speed
and direction of travel of vehicles i.e. any simulation programs that implements the
mobility model before stating the simulation of the network model will not satisfy our
needs.
According to the complexity of both tasks, we preferred using the well-known Matlab
commercial program for being popular, intuitive and easy to use. Matlab offers a full control
of all simulation parameters.
5.1 Performance Metrics
The performance metrics used to validate the proposed protocol are:
Latency: the total time required measured from the first attempt to broadcast to the
complete of the broadcasting phase.
Collision Probability: the average probability of collision in the ACK messages in each
segment of the transmission range.
68
5.2 Measurement Methodology
Latency is measured in the simulation program using four distinct delay sources (Fig. 5-1);
1- Contention starting time: that is equal to DIFS + RTB transmission time
2- Success broadcasting time: in case of a following node replied with a valid CTB
message, the time required for the rest of the broadcasting phase, equals to SIFS + CTB
transmission time + SIFS + message transmission time + SIFS + ACK transmission time.
3- Collision time: in case of a collision in the CTB packet, the time wasted will be equal to
DIFS + RTB transmission time + SIFS + CTB transmission time + SIFS. Note that a collision
destroys a complete phase.
4- Waiting time: the time taken by a node to decide either it will reply with a CTB message
or not, each instance of this time equals to a single time slot.
Fig. 5-1. RTB/CTB/data/ACK timeline
Note that: event (1) and (2) must happen once per any broadcasting phase, however, event
(3) and (4) may happen with a variable number according to the protocol design and the
network condition. There is a tradeoff between events (3) and (4) in all broadcasting protocols
that depends on segmentation of the transmission area (e.g. [12] and [30]); increasing the
number of segments leads to smaller number of vehicles in each segment, hence, decreases
collision probability while increases waiting times, and vice versa.
Assuming that both CW
min
and CW
max
equal to one (i.e. there is no contention or random
backoff), collision probability is measured in the simulation program by dividing the number
of broadcast phases that happen to have more than one node in the first non-empty segment
by the total number of broadcast phases.
SIFS
DIFS
SIFS
SIFS
RTB DATA Transmitter
Receiver
CTB ACK
69
5.3 Simulation Parameters
Table 5-1 summarizes parameters that are taken during simulation. These parameters are
quoted from the 802.11p [22] standard. It is assumed that the length of ACK message is the
same as the original broadcast. The ACK message is a mere repeat of the original broadcast
setting the ACK field, which is considered as a compensation to the expected collisions at far
range nodes (Sec. 4.3.1).
Table 5-1 Matlab parameters
Time-Slot 16 s CTB 14 bytes
SIFS 32 s Messages 512 bytes
DIFS 64 s ACK 512 bytes
RTB 20 bytes Data rate 3 Mbps
5.4 Random Number Generator
The problem of a new mobility model is treated with a random number generator to be
used as the headway between vehicles of the same lane. The histogram of one of these
variables is shown in Fig. 5-2 (at 2000 sample) indicating the flexibility of Matlab
environment in the simulation. This headway distribution is considered at a traffic volume of
approximately 330 veh/h to be identical to the one used for the analytical analysis (Fig. 4-10).
5.5 Simulation Results
Using these random variables, a simulation program was conducted for estimating the
probability of collisions and the average latency within each segment of the communication
range (10 sec). The width of each segment is taken according to Table 4-1. The probability of
collision is shown in Fig. 5-3 while the average latency is shown in Fig. 5-4.
For clearness, we have to repeat the following note: In these figures, values are measured
as an average on the segment, and are plotted on the last point of it.
Comparing simulation results with analytical results (Fig. 4-15) proves the correctness of
the analytical analysis.
70
Fig. 5-2. Histogram of one of the variables
Fig. 5-3. Simulated calculation of Pc for best segmentation
0 2 4 6 8 10 12
0
0.05
0.1
0.15
0.2
0.25
0.3
0.35
headway (sec)
P
r
o
b
a
b
i
l
i
t
y
o
f
C
o
l
l
i
s
i
o
n
Pc for non-uniform seg at 330v/h (simulation)
10 seg
9 seg
8 seg
7 seg
6 seg
5 seg
4 seg
0 2 4 6 8 10
0
10
20
30
40
50
60
headway (sec)
N
u
m
b
e
r
o
f
s
a
m
p
l
e
s
Histogram of X
71
Fig. 5-4. Simulated calculation of latency at best segmentation
In Fig. 5-4, the curves are close to each other, and thus, the curves were plotted only for
cases with even number of segments. The average latency associated with each segment
reveals that, the case of 6-seg gives the minimum latency (best performance) before over-
segmentation begins to take place with 8 and 10 segments.
5.6 Robustness at Different Traffic Volumes
The performance of almost all previously published protocols, changes drastically with
changing the traffic volume. However, the proposed protocol possesses unique robustness at
different traffic volumes. In this section, we will study the performance of the same protocol
under different traffic volumes. The headway distribution at traffic volume of 330veh/h and
1300 veh/h (very low vs. very high) are shown in Fig. 5-5 [32].
It could be seen that increasing the traffic volume results in increasing the ratio of short
headways and decreasing that of long headways.
0 2 4 6 8 10 12
2850
2900
2950
3000
3050
headway (sec)
L
a
t
e
n
c
y
(
s
e
c
)
Latency for non-uni seg at 330veh/h (simulation)
10 seg
8 seg
6 seg
4 seg
72
Fig. 5-5. Headway distribution at 330v/h and 1300v/h
Assuming a traffic volume of 1300 veh/h, we implemented the same simulation program
with the same segmentation table (Table 4-1) which was extracted for the best operation
under 330 veh/h. We compare the performance of this table with a new table calculated
especially for the new traffic volume (Table 5-2). By this comparison, we answer the question
should we change the segmentation table according to the current traffic volume or not?
Table 5-2 Best segmentation points for 1300 vehicle /h (in headway sec)
10-seg = [ 1.16 1.93 2.72 3.55 4.45 5.43 6.5 7.63 8.81 10];
9-seg = [ 1.28 2.13 3.01 3.95 4.98 6.13 7.36 8.67 10];
8-seg = [ 1.42 2.36 3.35 4.42 5.64 7 8.46 10];
7-seg = [ 1.59 2.66 3.79 5.07 6.57 8.22 10];
6-seg = [ 1.79 3.01 4.34 5.95 7.84 10];
5-seg = [ 2.04 3.47 5.14 7.33 10];
4-seg = [ 2.34 4.06 6.42 10];
Fig. 5-6 and Fig. 5-7 show the PC and latency at 1300veh/h for 6-seg computed with,
- Segmentation points that result in best operation under a traffic volume of 330v/h
- Newly calculated segmentation points for best operation under a traffic volume of 1300v/h
0 2 4 6 8 10
0
0.05
0.1
0.15
0.2
0.25
0.3
0.35
0.4
0.45
0.5
Headway (sec)
D
e
n
s
i
t
y
Headway model for 330 V/h and 1300 V/h
1300 V/h
330 V/h
73
Fig. 5-6. P
C
for 6-seg at 1300v/h
Fig. 5-7. Latency for 6-seg at 1300v/h
0 2 4 6 8 10 12
0
0.05
0.1
0.15
0.2
0.25
Headway (sec)
P
r
o
b
a
b
i
l
i
t
y
o
f
c
o
l
l
i
s
i
o
n
Pc for 6-seg at 1300 V/h
Best for 1300 V/h
Best for 330 V/h
0 2 4 6 8 10 12
2850
2900
2950
3000
3050
3100
Headway (sec)
L
a
t
e
n
c
y
(
s
e
c
)
Average latency for 6-seg at 1300 V/h
Best for 1300 V/h
Best for 330 V/h
74
The difference in P
C
is in a range of only 0.025 and the difference in latency is in a range
of 10sec. Hence, the answer to the above question is simply no. This protocol possesses
high robustness at different traffic volumes that we do not need to recalculate the table of best
segmentation points for each traffic volume.
5.7 Protocol Comparison
The last section of this chapter is about comparing the performance of the proposed
protocol with previously published ones. We performed the same simulation analysis to
compare the performance of the proposed protocol with The Urban Multihop Broadcast
Protocol (UMB) [30] and The Smart Broadcasting Protocol (SB) [12] protocols. The
difficulty with this comparison is that the proposed protocol is the first one to use the concept
of headway-based segmentation. To accomplish this comparison, we assumed that the speed
is uniform (i.e., neglecting the effect of the proposed headway-based segmentation). The
objective of this comparison becomes only to study the effect of non-uniform segmentation
on the performance. This is to illustrate that the proposed protocol has uniquely succeeded in
achieving a linearly increasing latency with a minimum slope. Despite the neglecting of
headway, the proposed protocol has two new features: it uses slightly different positions for
segmentation and a reversed order of priority.
Fig. 5-8 shows the collision probability per each segment with respect to each protocol,
while Fig. 5-9 shows the average latency per segment.
Fig. 5-8 shows the superior performance of our protocol with respect to P
C
while lines of
both SB and UMB coincide on each other. Fig. 5-9 shows that both versions (uniform and
non-uniform segmentation) of our protocol perform better than SB and UMB. Note that
the non-uniform segmentation mimics danger situation better than uniform segmentation,
where the latency is required to be directly proportional to the headway.
This section showed that, even if we neglected the effect headway-based segmentation (the
core contribution), the proposed protocol still outperforms previously presented ones.
75
Fig. 5-8. Probability of Collision (protocol comparison)
Fig. 5-9. Latency (protocol comparison)
0 2 4 6 8 10 12
0.02
0.04
0.06
0.08
0.1
0.12
0.14
0.16
headway (sec)
P
r
o
b
a
b
i
l
i
t
y
o
f
C
o
l
l
i
s
i
o
n
Comparison of Pc at 330 V/h
6-seg non-uni
6-seg uni
6-seg SB
6-seg UMB
0 2 4 6 8 10 12
2850
2900
2950
3000
headway (sec)
L
a
t
e
n
c
y
(
s
e
c
)
Comparison of average latency at 330 V/h
6-seg non-uni
6-seg uni
6-seg SB
6-seg UMB
76
Chapter 6
Conclusion
In this research, we introduced a novel broadcasting protocol in VANET environments
with these new features:
The first protocol to use the concept of headway-based segmentation and to include
effects of human behaviors in its design with the headway model.
Non-uniform segmentation achieving a unique a minimum slope linearly increasing
latency distribution.
Unique robustness at different speeds and traffic volumes rooted to the headway
robustness at different traffic volume variations. For example the latency difference between
the traffic volume of 330veh/h and 1300veh/h is in a range of 10sec.
Superior minimum latency for public safety applications.
Application adaptability with special multi-mode operations.
Considered offering a solution to applications never discussed in literature, like
Approaching Emergency Vehicle.
77
Appendix A - List of Co-authored Publications
[1] VANET-DSRC Protocol for Reliable Broadcasting of Life Safety Messages, In Proc.
of the 7
th
IEEE Int. Symposium on Signal Processing and Information Technology ISSPIT07,
Cairo, EGYPT, pp. 104-109, Dec.2007.
Available: http://ieeexplore.ieee.org/xpl/freeabs_all.jsp?arnumber=4458046
[2] Integrated intra-vehicle VANET system for increasing the road safety, in Proc. of
the Global Knowledge Forum NOOR-2008, Almadina, KSA, June 2008.
[3] A Novel Headway-Based Vehicle-to-Vehicle Multi-Mode Broadcasting Protocol
accepted for publication in Proc. of the 68
th
IEEE Vehicular Technology Conf. VTC2008-
Fall, Calgary, Alberta, CANADA, 2124 September 2008.
78
Appendix B - Word-Wide VANET Projects
USA:
DSRC http://www.leearmstrong.com/DSRC/DSRCHomeset.htm
http://grouper.ieee.org/groups/scc32/index.html
California PATH http://www-path.eecs.berkeley.edu/
DynaMIT http://mit.edu/its/index.html
ITS Research Program http://www.ivhs.washington.edu/
CITranS http://citrans.pti.psu.edu/
CISR http://www.cisr.gwu.edu/
TIDE http://www.njtide.org/
ODU VANET http://www.cs.odu.edu/~vanet/index.html
Europe:
AIDE http://www.aide-eu.org/
Car to Car http://www.car-2-car.org/
CarTALK 2000 http://www.cartalk2000.net/
CVIS http://www.cvisproject.org/
FleetNet http://www.et2.tu-harburg.de/fleetnet/index.html
Invent http://www.invent-online.de/
Network on Wheels http://www.network-on-wheels.de/
SEVECOM http://www.sevecom.org/
SAFESPOT http://www.safespot-eu.org/pages/page.php
SmartPark http://smartpark.epfl.ch/
WATCH-OVER http://www.watchover-eu.org/
WILLWARN http://www.prevent-ip.org/en/prevent_subprojects/
safe_speed_and_safe_following/willwarn/
Japan:
ITS Consortium http://www.its-jp.org/english/
79
Appendix C - VANET Simulation Programs
Network Simulators
NS-2 http://www.isi.edu/nsnam/ns/
GloMoSim http://pcl.cs.ucla.edu/projects/glomosim/
QualNet http://www.scalable-networks.com/
OPNet http://www.opnet.com/
NCTUns http://nsl.csie.nctu.edu.tw/nctuns.html
MATLAB http://www.mathworks.com/
VANET Mobility Generators
VanetMobiSim http://vanet.eurecom.fr/
CanuMobiSim http://canu.informatik.uni-stuttgart.de/mobisim/
Joint Mobility and network simulators for VANET
TraNS http://trans.epfl.ch/
MOVE http://www.csie.ncku.edu.tw/~klan/move/index.htm
80
Appendix D - MATLAB Scripts
Program set #1:
%Thi s f unct i on gener at es and pl ot s t he Semi - Poi sson Di st r i but i on
%x must be equal t o ei t her 330 or 1300
%- - - - - - - - - - - - - - - - - - - - - - - - - - - - -
x=i nput ( ' Ent er t he r equi r ed vehi cl e di st r i but i on 330 or 1300: ' ) ;
i f ( x==330)
P=0. 14; %par amet er s f or appr oxi mat e 330 veh/ h
B=2. 929;
A=5. 488;
T=0. 06;
el sei f ( x==1300)
P=. 64; %par amet er s f or appr oxi mat e 1300 veh/ h
B=2. 929;
A=5. 488;
T=. 18;
el se
er r or ( ' t he val i d i nput s ar e onl y 330 or 1300' ) ;
end
f =[ ] ;
f or t =0: 0. 01: 9. 99 %1000 sampl e i n 10 sec
r =P*( B*t ) ^( A- 1) / gamma( A) *B*exp( - B*t ) +( 1-
P) *gammai nc( A*t , B) / gamma( B) *( 1+T/ A) ^B*T*exp( - T*t ) ;
f =[ f r ] ;
end
f =f / ( . 01*t r apz( f , 2) ) ; %nor mal i zi ng t he sumt o 1
pl ot ( [ 0: . 01: 9. 99] , f , ' - b' )
xl abel ( ' headway ( sec) ' )
yl abel ( ' Densi t y %' )
81
Program set #2:
%Thi s f unct i on cal cul at es t he best poi nt s of segment at i on
%- - - - - - - - - - - - - - - - - - headway pdf - - - - - - - - - - - - - - - -
x=i nput ( ' Ent er t he r equi r ed vehi cl e di st r i but i on 330 or 1300: ' ) ;
i f ( x==330)
P=0. 14; %par amet er s f or appr oxi mat e 330 veh/ h
B=2. 929;
A=5. 488;
T=0. 06;
f 0=0. 33864;
el sei f ( x==1300)
P=. 64; %par amet er s f or appr oxi mat e 1300 veh/ h
B=2. 929;
A=5. 488;
T=. 18;
f 0=0. 79708;
el se
er r or ( ' t he val i d i nput s ar e onl y 330 or 1300' ) ;
end
f =[ ] ;
f or t =0: 0. 01: 29. 99 %3000 sampl e i n 10 sec
r =P*( B*t ) ^( A- 1) / gamma( A) *B*exp( - B*t ) +( 1-
P) *gammai nc( A*t , B) / gamma( B) *( 1+T/ A) ^B*T*exp( - T*t ) ;
f =[ f r ] ;
end
f =f / f 0; %nor mal i zi ng t he summat i on t o 1
%%- - - - - - - - - - - - - - - - Run f or al l
opt i ons = opt i mset ( ' Di spl ay' , ' of f ' ) ; %opt i mi zat i on par amet er s
H_new10=0;
Ns=10; %number of segment s
Best _H1=f mi nsear ch( @opt mm, 1, opt i ons, f , Ns) ; %mi n. f unc. on t he f i r st poi nt
H_new10=ext r act _H( Best _H1, f , Ns) ; %t he compl et e set of poi nt s
di sp ( ' done f or 10- seg' ) ;
Pc10=ext r act _Pc( H_new10, f ) ; %Pc f or t he set of poi nt s
H_new9=0;
Ns=9; %same f unct i on used f or 9- seg
Best _H1=f mi nsear ch( @opt mm, 1, opt i ons, f , Ns) ;
H_new9=ext r act _H( Best _H1, f , Ns) ;
di sp ( ' done f or 9- seg' ) ;
Pc9=ext r act _Pc( H_new9, f ) ;
H_new8=0;
Ns=8; %same f unct i on used f or 8- seg
Best _H1=f mi nsear ch( @opt mm, 1, opt i ons, f , Ns) ;
H_new8=ext r act _H( Best _H1, f , Ns) ;
di sp ( ' done f or 8- seg' ) ;
Pc8=ext r act _Pc( H_new8, f ) ;
H_new7=0;
82
Ns=7; %same f unct i on used f or 7- seg
Best _H1=f mi nsear ch( @opt mm, 1, opt i ons, f , Ns) ;
H_new7=ext r act _H( Best _H1, f , Ns) ;
di sp ( ' done f or 7- seg' ) ;
Pc7=ext r act _Pc( H_new7, f ) ;
H_new6=0;
Ns=6; %same f unct i on used f or 6- seg
Best _H1=f mi nsear ch( @opt mm, 1, opt i ons, f , Ns) ;
H_new6=ext r act _H( Best _H1, f , Ns) ;
di sp ( ' done f or 6- seg' ) ;
Pc6=ext r act _Pc( H_new6, f ) ;
H_new5=0;
Ns=5; %same f unct i on used f or 5- seg
Best _H1=f mi nsear ch( @opt mm, 2, opt i ons, f , Ns) ;
H_new5=ext r act _H( Best _H1, f , Ns) ;
di sp ( ' done f or 5- seg' ) ;
Pc5=ext r act _Pc( H_new5, f ) ;
H_new4=0;
Ns=4; %same f unct i on used f or 4- seg
Best _H1=f mi nsear ch( @opt mm, 2, opt i ons, f , Ns) ;
H_new4=ext r act _H( Best _H1, f , Ns) ;
di sp ( ' done f or 4- seg' ) ;
Pc4=ext r act _Pc( H_new4, f ) ;
di sp ( ' The best poi nt s of segment at i on ar e' )
H_new10
H_new9
H_new8
H_new7
H_new6
H_new5
H_new4
f i gur e
hol d
pl ot ( H_new10, Pc10, ' - bo' )
pl ot ( H_new9, Pc9, ' - gx' )
pl ot ( H_new8, Pc8, ' - r +' )
pl ot ( H_new7, Pc7, ' - c*' )
pl ot ( H_new6, Pc6, ' - ms' )
pl ot ( H_new5, Pc5, ' - yd' )
pl ot ( H_new4, Pc4, ' - kv' )
l egend( ' 10- seg' , ' 9- seg' , ' 8- seg' , ' 7- seg' , ' 6- seg' , ' 5- seg' , ' 4-
seg' , ' Locat i on' , ' best ' )
xl abel ( ' Headway ( sec) ' )
yl abel ( ' Pr obabi l i t y of Col l i si on' )
83
f unct i on D=opt mm( H1, f , Ns)
%Thi s f unct i on i s used by t he Mat l ab opt i mi zat i on pr ogr am
%I t accept s t he f i r st poi nt of segment at i on ( H1) and t he used
%di st r i but i on ( f ) and t he number of segment s( Ns)
%I t t r i es t o f i nd a l i near l y i ncr easi ng sl op st ar t i ng at H1 and ends at
%t he exat l y 10
H=zer os( 1, Ns) ;
H( 1) =r ound( H1*100) ;
Px=sum( f ( 1: H( 1) ) ) ;
Pc=0; %f i nd sl op of Pc f or f i r st segment
f or x=1: H( 1) - 1
Pc=Pc+( f ( x) / Px) *sum( f ( 1: H( 1) - x) ) ;
end
S=Pc/ H( 1) ;
f or i =1: Ns- 1;
f or Next _H=H( i ) +1: 2000 %l oop on f i ndi ng next H
Px=sum( f ( H( i ) : Next _H) ) ;
Pc=0;
f or x=H( i ) : Next _H- 1
Pc=Pc+( f ( x) / Px) *sum( f ( 1: Next _H- x) ) ;
end
i f ( Pc/ Next _H) >= S %f i nd next H havi ng same sl ope
H( i +1) =Next _H;
br eak
end
end
end
D=abs( H( Ns) - 1000) ;
f unct i on Pc=ext r act _Pc ( H, f )
%ext r ac Pc f or Best _H
L=l engt h( H) ;
H=r ound( H. *100) ;
H=[ 1 H] ;
Pc=zer os( 1, L) ;
f or i =1: L
Px=sum( f ( H( i ) : H( i +1) ) ) ;
f or x=H( i ) : H( i +1)
Pc( i ) =Pc( i ) +( ( f ( x) / Px) *sum( f ( 1: ( H( i +1) - x) ) ) ) ;
end
end
Pc=Pc/ 100;
f unct i on H=ext r act _H( H1, f , Ns)
%ext r act al l poi nt s of segment at i on f r omt he f i r st one
H=zer os( 1, Ns) ;
84
H( 1) =r ound( H1*100) ;
Px=sum( f ( 1: H( 1) ) ) ;
Pc=0; %f i nd sl op of Pc f or f i r st segment
f or x=1: H( 1) - 1
Pc=Pc+( f ( x) / Px) *sum( f ( 1: H( 1) - x) ) ;
end
S=Pc/ H( 1) ;
f or i =1: Ns- 1;
f or Next _H=H( i ) +1: 2000 %l oop on f i ndi ng next H
Px=sum( f ( H( i ) : Next _H) ) ;
Pc=0;
f or x=H( i ) : Next _H- 1
Pc=Pc+( f ( x) / Px) *sum( f ( 1: Next _H- x) ) ;
end
i f ( Pc/ Next _H) >= S %f i nd next H havi ng same sl ope
H( i +1) =Next _H;
br eak
end
end
end
H=H/ 100;
85
Program set #3:
%Thi s f unct i on comput e t he aver age PC and l at ency f or t he best poi nt s of
%segment at i on i n 330 and 1300 v/ h
x=i nput ( ' Ent er t he r equi r ed vehi cl e di st r i but i on 330 or 1300: ' ) ;
i f ( x==330)
seg10=[ 1. 01 1. 74 2. 51 3. 34 4. 26 5. 27 6. 36 7. 52
8. 75 10. 05] ;
seg9 =[ 1. 12 1. 93 2. 79 3. 73 4. 79 5. 96 7. 23 8. 58
10. 02] ;
seg8 =[ 1. 25 2. 15 3. 12 4. 22 5. 47 6. 85 8. 35 9. 96] ;
seg7 =[ 1. 42 2. 45 3. 59 4. 93 6. 47 8. 18 10. 06] ;
seg6 =[ 1. 62 2. 81 4. 18 5. 86 7. 81 10. 04] ;
seg5 =[ 1. 86 3. 27 5. 03 7. 28 10. 02] ;
seg4 =[ 2. 15 3. 88 6. 38 9. 95] ;
el sei f ( x==1300)
seg10=[ 1. 16 1. 93 2. 72 3. 55 4. 45 5. 43 6. 5 7. 63 8. 81
10] ;
seg9 =[ 1. 28 2. 13 3. 01 3. 95 4. 98 6. 13 7. 36 8. 67 10] ;
seg8 =[ 1. 42 2. 36 3. 35 4. 42 5. 64 7 8. 46 10] ;
seg7 =[ 1. 59 2. 66 3. 79 5. 07 6. 57 8. 22 10] ;
seg6 =[ 1. 79 3. 01 4. 34 5. 95 7. 84 10] ;
seg5 =[ 2. 04 3. 47 5. 14 7. 33 10] ;
seg4 =[ 2. 34 4. 06 6. 42 10] ;
el se
er r or ( ' t he val i d i nput s ar e onl y 330 or 1300' ) ;
end
d10=zer os( 50, l engt h( seg10) ) ; %pr e al l ocat i on f or bet t er speed
c10=zer os( 50, l engt h( seg10) ) ;
d9=zer os( 50, l engt h( seg9) ) ; %d: del ay
c9=zer os( 50, l engt h( seg9) ) ; %c: col l i si ons
d8=zer os( 50, l engt h( seg8) ) ;
c8=zer os( 50, l engt h( seg8) ) ;
d7=zer os( 50, l engt h( seg7) ) ;
c7=zer os( 50, l engt h( seg7) ) ;
d6=zer os( 50, l engt h( seg6) ) ;
c6=zer os( 50, l engt h( seg6) ) ;
d5=zer os( 50, l engt h( seg5) ) ;
c5=zer os( 50, l engt h( seg5) ) ;
d4=zer os( 50, l engt h( seg4) ) ;
c4=zer os( 50, l engt h( seg4) ) ;
f or i =1: 50
[ x1] =headway( x) ;
[ d10( i , : ) c10( i , : ) ] =MMB( x1, seg10) ;
[ d9( i , : ) c9( i , : ) ] =MMB( x1, seg9) ;
[ d8( i , : ) c8( i , : ) ] =MMB( x1, seg8) ;
[ d7( i , : ) c7( i , : ) ] =MMB( x1, seg7) ;
[ d6( i , : ) c6( i , : ) ] =MMB( x1, seg6) ;
[ d5( i , : ) c5( i , : ) ] =MMB( x1, seg5) ;
[ d4( i , : ) c4( i , : ) ] =MMB( x1, seg4) ;
end
86
Del ay10=( mean( d10) ) ' ;
Pc10=( mean( c10) ) ' ;
Del ay9=( mean( d9) ) ' ;
Pc9=( mean( c9) ) ' ;
Del ay8=( mean( d8) ) ' ;
Pc8=( mean( c8) ) ' ;
Del ay7=( mean( d7) ) ' ;
Pc7=( mean( c7) ) ' ;
Del ay6=( mean( d6) ) ' ;
Pc6=( mean( c6) ) ' ;
Del ay5=( mean( d5) ) ' ;
Pc5=( mean( c5) ) ' ;
Del ay4=( mean( d4) ) ' ;
Pc4=( mean( c4) ) ' ;
f i gur e
hol d
pl ot ( seg10, Del ay10, ' - bo' )
pl ot ( seg9, Del ay9, ' - gx' )
pl ot ( seg8, Del ay8, ' - r +' )
pl ot ( seg7, Del ay7, ' - c*' )
pl ot ( seg6, Del ay6, ' - ms' )
pl ot ( seg5, Del ay5, ' - yd' )
pl ot ( seg4, Del ay4, ' - kv' )
l egend( ' 10 seg' , ' 9 seg' , ' 8 seg' , ' 7 seg' , ' 6 seg' , ' 5 seg' , ' 4
seg' , ' Locat i on' , ' best ' )
xl abel ( ' headway ( sec) ' )
yl abel ( ' Lat ency ( usec) ' )
f i gur e
hol d
pl ot ( seg10, Pc10, ' - bo' )
pl ot ( seg9, Pc9, ' - gx' )
pl ot ( seg8, Pc8, ' - r +' )
pl ot ( seg7, Pc7, ' - c*' )
pl ot ( seg6, Pc6, ' - ms' )
pl ot ( seg5, Pc5, ' - yd' )
pl ot ( seg4, Pc4, ' - kv' )
l egend( ' 10 seg' , ' 9 seg' , ' 8 seg' , ' 7 seg' , ' 6 seg' , ' 5 seg' , ' 4
seg' , ' Locat i on' , ' best ' )
xl abel ( ' headway ( sec) ' )
yl abel ( ' Pr obabi l i t y of Col l i si on' )
f unct i on [ x1] =headway( x)
%Headway Model and r andomnumber gener at or
%- - - - - - - - - - - - - - - - - - - - - - - - - - - - -
i f ( x==330)
P=0. 14; %par amet er s f or f or appr oxi mat e 330 veh/ h f or t est i ng of MMB
B=2. 929;
A=5. 488;
87
T=0. 06;
el sei f ( x==1300)
P=. 64; %par amet er s f or f or appr oxi mat e 1300 veh/ h f or t est i ng of MMB
B=2. 929;
A=5. 488;
T=. 18;
end
f =zer os( 1, 101) ;
i =1;
f or t =0: 0. 1: 10
f ( i ) =P*( B*t ) ^( A- 1) / gamma( A) *B*exp( - B*t ) +( 1-
P) *gammai nc( A*t , B) / gamma( B) *( 1+T/ A) ^B*T*exp( - T*t ) ;
i =i +1;
end
f =f / t r apz( f , 2) ;
%- - - - - - - - - - - - - - - - - - - - - - - - - - - - -
%Randomnumber gener at or
%- - - - - - - - - - - - - - - - - - - - - - - - - - - - -
f =f *10;
x1=zer os( 2000, 2) ;
y=zer os( 101, 5) ;
f or j =1: 2
i =1;
whi l e i <2000
r =10*r and;
c=cei l ( 10*r ) ;
i f y( c, j ) < 200*f ( c)
y( c, j ) =y( c, j ) +1;
x1( i , j ) =r ;
i =i +1;
end
end
x1( : , j ) =shuf f l e( x1( : , j ) ) ;
end
x1( : , 2) =x1( : , 1) +x1( : , 2) ;
f unct i on Y = shuf f l e( X)
%t hi s f unct i on i s used t o shuf f l e t he r andomnumber
Y=zer os( 1, l engt h( X) ) ;
f or I =1: l engt h( X)
A=cei l ( l engt h( X) *r and) ;
Y( I ) =X( A) ;
X( A) =[ ] ;
end
f unct i on [ d, c] =MMB( x1, seg)
%t he pr oposed pr ot ocol
%I t accept s t he r andomnumber and t he segment at i on
%- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
88
T0=64+50. 86; %DI FS+RTB RTB=20byt e/ 3Mbps
TB=32+35. 604+32+1302. 1+32+1302. 1; %SI FS+CTB+SI FS+Dat a+SI FS+ACK
CTB=14byt e/ 3Mbps message=512byt e / 3Mbps
TI =16; %si ngl e t i me- sl ot
TC=64+50. 86+32+35. 604; %DI FS+RTB+SI FS+CTB
L_seg=l engt h( seg) ;
seg=[ 0 seg] ;
Count er =zer os( 1, L_seg) ; %same l engt h as f
Del ay=zer os( 1, L_seg) ;
C_col l =zer os( 1, L_seg) ;
f or xc1=1: 2000 %l engt h of x1
t =T0;
vehi cl e=x1( xc1, : ) ;
f or i =1: L_seg
n=sum( ( vehi cl e>seg( i ) ) & ( vehi cl e<seg( i +1) ) ) ;
i f ( n==0)
t =t +TI ;
el sei f ( n==1)
t =t +TB;
Del ay( i ) =Del ay( i ) +t ;
Count er ( i ) =Count er ( i ) +1;
br eak
el sei f ( n>1)
t =t +TC;
Del ay( i ) =Del ay( i ) +t ;
C_col l ( i ) =C_col l ( i ) +1;
br eak
end
end
end
d=Del ay. / Count er ;
c=C_col l . / ( Count er +C_col l ) ;
89
Pr ogr amset #4:
%Thi s f unct i on i s t p compar e t he per f or mance of t he pr oposed pr ot ocol wi t h
%ot her publ i shed ones
seg6=[ 1. 62 2. 81 4. 18 5. 86 7. 81 10. 04] ;
segu=[ 10/ 6 20/ 6 30/ 6 40/ 6 50/ 6 10] ;
d6=zer os( 200, l engt h( seg6) ) ;
c6=zer os( 200, l engt h( seg6) ) ;
du=zer os( 200, l engt h( segu) ) ;
cu=zer os( 200, l engt h( segu) ) ;
ds=zer os( 200, l engt h( segu) ) ;
cs=zer os( 200, l engt h( segu) ) ;
dU=zer os( 200, l engt h( segu) ) ;
cU=zer os( 200, l engt h( segu) ) ;
f or i =1: 200 %r edo t he pr ogr amf or 200 t i mes
[ x1] =headwaySB;
[ d6( i , : ) c6( i , : ) ] =MMB( x1, seg6) ;
[ du( i , : ) cu( i , : ) ] =MMB( x1, segu) ;
[ ds( i , : ) cs( i , : ) ] =SB( x1, segu) ;
[ dU( i , : ) cU( i , : ) ] =UMB( x1, segu) ;
end
Del ay6=( mean( d6) ) ' ;
Pc6=( mean( c6) ) ' ;
Del ayu=( mean( du) ) ' ;
Pcu=( mean( cu) ) ' ;
Del ays=( mean( ds) ) ' ;
Pcs=( mean( cs) ) ' ;
Del ayU=( mean( dU) ) ' ;
PcU=( mean( cU) ) ' ;
f i gur e
hol d
pl ot ( seg6, Del ay6, ' - ms' )
pl ot ( segu, Del ayu, ' - kd' )
pl ot ( segu, Del ays, ' - c*' )
pl ot ( segu, Del ayU, ' - gx' )
l egend( ' 6- seg non- uni ' , ' 6- seg uni ' , ' 6- seg SB' , ' 6- seg
UMB' , ' Locat i on' , ' best ' )
xl abel ( ' headway ( sec) ' )
yl abel ( ' Lat ency ( usec) ' )
f i gur e
hol d
pl ot ( seg6, Pc6, ' - ms' )
pl ot ( segu, Pcu, ' - kd' )
pl ot ( segu, Pcs, ' - c*' )
pl ot ( segu, PcU, ' - gx' )
l egend( ' 6- seg non- uni ' , ' 6- seg uni ' , ' 6- seg SB' , ' 6- seg
UMB' , ' Locat i on' , ' best ' )
xl abel ( ' headwaySB ( sec) ' )
yl abel ( ' Pr obabi l i t y of Col l i si on' )
90
f unct i on [ x1] =headwaySB( )
%Headway Model
%- - - - - - - - - - - - - - - - - - - - - - - - - - - - -
P=0. 14; %par amet er s f or f or appr oxi mat e 330 veh/ h f or t est i ng of MMB
B=2. 929;
A=5. 488;
T=0. 06;
f =zer os( 1, 101) ;
i =1;
f or t =0: 0. 1: 10
f ( i ) =P*( B*t ) ^( A- 1) / gamma( A) *B*exp( - B*t ) +( 1-
P) *gammai nc( A*t , B) / gamma( B) *( 1+T/ A) ^B*T*exp( - T*t ) ;
i =i +1;
end
f =f / t r apz( f , 2) ;
%t =0: 0. 1: 38;
%pl ot ( t , f )
%- - - - - - - - - - - - - - - - - - - - - - - - - - - - -
%Randomnumber gener at or
%- - - - - - - - - - - - - - - - - - - - - - - - - - - - -
f =f *10;
x1=zer os( 2000, 10) ;
y=zer os( 101, 10) ;
f or j =1: 10
i =1;
whi l e i <2000
r =10*r and;
c=cei l ( 10*r ) ;
i f y( c, j ) < 200*f ( c)
y( c, j ) =y( c, j ) +1;
x1( i , j ) =r ;
i =i +1;
end
end
x1( : , j ) =shuf f l e( x1( : , j ) ) ;
end
x1( : , 2) =x1( : , 1) +x1( : , 2) ; %f or meach 10 val ues ser i al l y
x1( : , 3) =x1( : , 2) +x1( : , 3) ;
x1( : , 4) =x1( : , 3) +x1( : , 4) ;
x1( : , 5) =x1( : , 4) +x1( : , 5) ;
x1( : , 6) =x1( : , 5) +x1( : , 6) ;
x1( : , 7) =x1( : , 6) +x1( : , 7) ;
x1( : , 8) =x1( : , 7) +x1( : , 8) ;
x1( : , 9) =x1( : , 8) +x1( : , 9) ;
x1( : , 10) =x1( : , 9) +x1( : , 10) ; %f our l ane
f unct i on Y = shuf f l e( X)
Y=zer os( 1, l engt h( X) ) ;
f or I =1: l engt h( X)
A=cei l ( l engt h( X) *r and) ;
91
Y( I ) =X( A) ;
X( A) =[ ] ;
end
f unct i on [ d, c] =MMB( x1, seg)
%The pr oposed pr ot ocol
%- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
T0=64+50. 86; %DI FS+RTB RTB=20byt e/ 3Mbps
TB=32+35. 604+32+1302. 1+32+1302. 1; %SI FS+CTB+SI FS+message+SI FS+ACK
CTB=14byt e/ 3Mbps message=512byt e / 3Mbps
TI =16; %si ngl e t i me- sl ot
TC=64+50. 86+32+35. 604; %DI FS+RTB+SI FS+CTB
L_seg=l engt h( seg) ;
seg=[ 0 seg] ;
Count er =zer os( 1, L_seg) ; %same l engt h as f
Del ay=zer os( 1, L_seg) ;
C_col l =zer os( 1, L_seg) ;
f or xc1=1: 2000 %l engt h of x1
t =T0;
vehi cl e=x1( xc1, : ) ;
f or i =1: L_seg
n=sum( ( vehi cl e>seg( i ) ) & ( vehi cl e<seg( i +1) ) ) ;
i f ( n==0)
t =t +TI ;
el sei f ( n==1)
t =t +TB;
Del ay( i ) =Del ay( i ) +t ;
Count er ( i ) =Count er ( i ) +1;
br eak
el sei f ( n>1)
t =t +TC;
Del ay( i ) =Del ay( i ) +t ;
C_col l ( i ) =C_col l ( i ) +1;
br eak
end
end
end
d=Del ay. / Count er ;
c=C_col l . / ( Count er +C_col l ) ;
f unct i on [ d, c] =SB( x1, seg)
%The Smar t pr ot ocol
%- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
T0=64+50. 86; %DI FS+RTB RTB=20byt e/ 3Mbps
TB=32+35. 604+32+1302. 1+32+1302. 1; %SI FS+CTB+SI FS+message+SI FS+ACK
CTB=14byt e/ 3Mbps message=512byt e / 3Mbps
92
TI =16; %si ngl e t i me- sl ot
TC=64+50. 86+32+35. 604; %DI FS+RTB+SI FS+CTB
L_seg=l engt h( seg) ;
seg=[ 0 seg] ;
Count er =zer os( 1, L_seg) ; %same l engt h as f
Del ay=zer os( 1, L_seg) ;
C_col l =zer os( 1, L_seg) ;
f or xc1=1: 2000 %l engt h of x1
t =T0;
vehi cl e=x1( xc1, : ) ;
f or i =L_seg+1: - 1: 2
n=sum( ( vehi cl e<seg( i ) ) & ( vehi cl e>seg( i - 1) ) ) ;
i f ( n==0)
t =t +TI ;
el sei f ( n==1)
t =t +TB;
Del ay( i - 1) =Del ay( i - 1) +t ;
Count er ( i - 1) =Count er ( i - 1) +1;
br eak
el sei f ( n>1)
t =t +TC;
Del ay( i - 1) =Del ay( i - 1) +t ;
C_col l ( i - 1) =C_col l ( i - 1) +1;
br eak
end
end
end
d=Del ay. / Count er ;
c=C_col l . / ( Count er +C_col l ) ;
f unct i on [ d, c] =UMB( x1, seg)
%The UMB pr ot ocol
%- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
T0=64+50. 86; %DI FS+RTB RTB=20byt e/ 3Mbps
TB=32+35. 604+32+1302. 1+32+1302. 1; %SI FS+CTB+SI FS+message+SI FS+ACK
CTB=14byt e/ 3Mbps message=512byt e / 3Mbps
TI =16; %si ngl e t i me- sl ot
TC=64+50. 86+32+35. 604; %DI FS+RTB+SI FS+CTB
L_seg=l engt h( seg) ;
seg=[ 0 seg] ;
Count er =zer os( 1, L_seg) ; %same l engt h as f
Del ay=zer os( 1, L_seg) ;
C_col l =zer os( 1, L_seg) ;
f or xc1=1: 2000 %l engt h of x1
93
t =T0;
vehi cl e=x1( xc1, : ) ;
f or i =L_seg+1: - 1: 2
n=sum( ( vehi cl e<seg( i ) ) & ( vehi cl e>seg( i - 1) ) ) ;
i f ( n==1)
t =t +( i - 1) *TI ;
t =t +TB;
Del ay( i - 1) =Del ay( i - 1) +t ;
Count er ( i - 1) =Count er ( i - 1) +1;
br eak
el sei f ( n>1)
t =t +( i - 1) *TI ;
t =t +TC;
Del ay( i - 1) =Del ay( i - 1) +t ;
C_col l ( i - 1) =C_col l ( i - 1) +1;
br eak
end
end
end
d=Del ay. / Count er ;
c=C_col l . / ( Count er +C_col l ) ;
94
Appendix E - References
[1] ARIB official homepage http://www.arib.or.jp/english/
[2] N. Balon, and J. Guo, Increasing Broadcast Reliability in Vehicular Ad Hoc Networks,
In Proc. of the 3
rd
ACM Int. Workshop on Vehicular Ad Hoc Networks VANET'06, NY,
USA, pp. 104-105, 2006.
[3] Bibliography on Secure Vehicular Communications
http://bbcr.uwaterloo.ca/~rxlu/sevecombib.htm
[4] B. Parno, and A. Perrig Challenges in securing vehicular networks, In Proc. of the Int.
Workshop on Hot Topics in Networks (HotNets-IV), 2005.
[5] Assessment of the applicability of cooperative vehicle-highway automation systems to
bus transit and intermodal freight: case study, California Partners for Advanced Transit
and Highways (PATH), 2004.
http://repositories.cdlib.org/cgi/viewcontent.cgi?article=1629&context=its/path
[6] Car-to-Car communication consortium http://www.car-2-car.org/
[7] CEN official homepage http://www.cen.eu/cenorm/homepage.htm
[8] Dedicated Short Range Communications (DSRC) Home,
http://www.leearmstrong.com/DSRC/DSRCHomeset.htm
http://grouper.ieee.org/groups/scc32/index.html
[9] "The official highway code," Department of Transport (DOT), Driving Standards Agency
(DSA), The Stationery Office, 2007.
[10] Statewide Plan, Department of Transportation (DOT), Intelligent Transportation
Systems (ITS), September 2005.
[11] Projects Book, Department of Transportation (DOT), Intelligent Transportation
Systems (ITS), January 2000.
[12] E. Fasolo, A. Zanella, and M. Zorzi, An Effective Broadcast Scheme for Alert Message
Propagation in Vehicular Ad hoc Networks, In Proc. of the IEEE Int. Conf. on
Communications ICC'06, vol. 9, pp. 3960-3965, 2006.
[13] GM advanced technology, Vehicle-to-Vehicle technology,
http://media.gm.com:8221/us/gm/en/technology/advanced_technology/safety_telematics/
V2V_communications.htm
[14] M. Green, How long does it take to stop? Methodological analysis of driver perception-
brake times, In Transportation Human Factors, vol.2, pp.195-216, 2000.
[15] H. Alshaer and E. Horlait, An optimized adaptive broadcast scheme for inter-vehicle
communication, In Proc. of the 61
st
Int. IEEE Vehicular Technology Conf. VTC05,
vol.5, pp. 28402844, 2005.
[16] L. Huang, A. Arora, and TH. Lai, Reliable MAC layer multicast in IEEE 802.11
wireless networks, In Proc. of the Int. Conf. on Parallel Processing ICPP'02, USA, 2002.
95
[17] IEEE 802.15.4a, Wireless medium access control (MAC) and physical layer (PHY)
specifications for low-rate wireless personal area networks (WPANs), amendment 1: add
alternate PHYs, August 2007.
[18] IEEE P1609.1, Trial-use standard for wireless access in vehicular environments
(WAVE) - resource manager, Draft D17, July 2006.
[19] IEEE P1609.2, Trial-use standard for wireless access in vehicular environments
(WAVE) - security services for applications and management messages, Draft D7, April
2006.
[20] IEEE P1609.3, Trial-use standard for wireless access in vehicular environments
(WAVE) - networking services, Draft D22, January 2007.
[21] IEEE P1609.4, Trial-use standard for wireless access in vehicular environments
(WAVE) - multi-channel operation, Draft D08, July 2006.
[22] IEEE P802.11p, Amendment 3: wireless access in vehicular environments (WAVE),
Draft D0.26, January 2006.
[23] IEEE P802.11-REVma
TM
/D7.0, Wireless LAN Medium Access Control (MAC) and
Physical Layer (PHY) Specifications, Revision of 802.11-1999, June 2006.
[24] IEEE P802.11s, Amendment: ESS Mesh Networking, Draft D1.03, April 2007.
[25] Internet Protocol Version 6 (IPv6) Specification, IETF Request for Comments, RFC
2460.
[26] Intelligent transportation systems (ITS) Home, http://www.its.dot.gov
[27] ITS Consortium http://www.its-jp.org/english
[28] J. Guo and N. Balon, "Vehicular ad hoc networks and dedicated short-range
communication," Book Chapter, 2006.
[29] J. Harri, F. Filali, and C. Bonnet, Mobility models for vehicular ad hoc networks: a
survey and taxonomy, The Institut Eurecom Department of Mobile Communications,
FRANCE, 2007.
[30] G. Korkmaz, E. Ekici, F. zgner, and U. zgner, "Urban multi-hop broadcast protocol
for inter-vehicle communication systems," In Proc. of the 1
st
ACM Int. Workshop on
Vehicular Ad Hoc Networks VANET'04, NY, USA, pp. 76-85, 2004.
[31] Largest Tiny Network Yet, http://webs.cs.berkeley.edu/800demo/
[32] T. Luttinen, Statistical analysis of vehicle time headways, Teknillinen korkeakoulu,
pp.155-172, 1996
[33] M. Heissenbttel, T. Braun, M. Wlchli, and T. Bernoulli, Optimized stateless
broadcasting in wireless multi-hop networks, in IEEE Infocom 2006, Barcelona, April
2006.
[34] M. Torrent-Moreno, D. Jiang, and H. Hartenstein, Broadcast reception rates and effects
of priority access in 802.11-based vehicular ad-hoc networks, in Proc. of the of the 1
st
ACM Int. Workshop on Vehicular Ad Hoc Networks VANET04, pp.10-18, 2004.
96
[35] M. Heron, D. Hoyert, J. Xu, C. Scott, and B. Tejada-Vera Deaths: Preliminary Data for
2006, Division of vital statistics, National Vital Statistics Reports, Vol.56, N.16, 2008.
[36] Fcc allocates spectrum in 5.9 ghz range for intelligent transportation systems uses,
News, Federal Communications Commission, 1999.
http://www.fcc.gov/Bureaus/Engineering_Technology/News_Releases/1999/nret9006.ht
ml
[37] S.Y. Ni, Y.C. Tseng, Y.S. Chen, and J.P. Sheu, "The broadcast storm problem in a
mobile ad hoc network," In Proc. of the 5th ACM/IEEE Int. Con. on Mobile Computing
and Networking MobiCom'99, NY, USA, pp. 151-162, 1999.
[38] Official IEEE 802.11 working group project timelines.
http://grouper.ieee.org/groups/802/11/Reports/802.11_Timelines.htm
[39] K. Rmer, and F. Mattern, The design space of wireless sensor networks, In Proc. of
the IEEE Wireless Communications and Networking Conf., vol.11, No.6, pp. 54-61,
2004.
[40] Status of Project IEEE 802.11 Task Group p.
http://grouper.ieee.org/groups/802/11/Reports/tgp_update.htm
[41] M. Taieb-Maimon, and D. Shinar, Minimum and Comfortable Driving Headways:
Reality versus Perception, In Human Factors, The Journal of the Human Factors and
Ergonomics Society, vol.43. pp.159-172, 2001.
[42] K. Tang, M. Gerla, MAC reliable broadcast in ad hoc networks, In Proc. of
Communications for Network-Centric Operations, Creating the Information Force, IEEE
Military Communications Conf. vol.2, pp.1008-1013, 2001.
[43] London Congestion Charging Technology Trials, Transport for London, Feb 2005.
http://www.tfl.gov.uk/assets/downloads/technology-trials.pdf
[44] USA National Transportation Statistics 2007, Bureau of Transportation Statistics,
USA, 2007.
[45] J. Xie, A. Das, S. Nandi, and A.K. Gupta, Improving the reliability of IEEE 802.11
broadcast scheme for multicasting in mobile ad hoc networks, In Proc. of the IEEE
Wireless Communications and Networking Conf., vol.1, pp.126-131. 2005.
[46] Q. Xu, T. Mak, J. Ko, and R. Sengupta, "Vehicle-to-vehicle safety messaging in DSRC,"
In Proc. of the 1
st
ACM Int. Workshop on Vehicular Ad Hoc Networks VANET'04, NY,
USA, pp.19-28, 2004.
[47] X. Yang, L. Liu, N.H. Vaidya, and F. Zhao, "A vehicle-to-vehicle communication
protocol for cooperative collision warning," In Proc. of the 1
st
Int. Conf. on Networking
and Services MOBIQUITOUS04, pp.114-123, 2004.
[48] A. Zanella, G. Pierobon, and S. Merlin, On the limiting performance of broadcast
algorithms over unidimensional ad-hoc radio networks, in Proc. of 7
th
Int. Symposium
on Wireless Personal Multimedia Communications, WPMC04, Abano Terme, Italy,
2004.
. ) - (
) VANETs (
.
) IEEE802.11 (
(IEEE802.11p)
.
.
: ) ACK (
.
) ACK ( .
.
) - (
/
2008
: -
. /
) (
. /
) (
. /
) (
/
) (
: -
. /
) (
/
) (
/
) (
) - (
/
2008