Next Article in Journal
A Dual-Perforation Electromagnetic Bandgap Structure for Parallel-Plate Noise Suppression in Thin and Low-Cost Printed Circuit Boards
Previous Article in Journal
A Survey on RF and Microwave Doherty Power Amplifier for Mobile Handset Applications
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

3D Spatial Characteristics of C-V2X Communication Interference

1
School of Electronic Information Engineering, Yangtze Normal University, Chongqing 408100, China
2
College of Communication Engineering, Chongqing University, Chongqing 400044, China
3
College of Big Data and Intelligent Engineering, Yangtze Normal University, Chongqing 408100, China
*
Author to whom correspondence should be addressed.
Electronics 2019, 8(6), 718; https://doi.org/10.3390/electronics8060718
Submission received: 11 May 2019 / Revised: 19 June 2019 / Accepted: 21 June 2019 / Published: 25 June 2019
(This article belongs to the Section Microwave and Wireless Communications)

Abstract

:
In C-V2X (cellular vehicle-to-everything) communication networks, dense spatial reuse of the available radio spectrum is required to achieve efficient spectral usage. Spectrum reuse causes severe network interference where signals from many undesired transmitters are aggregated at a receiver. This paper investigates the 3D spatial characteristics of C-V2X communication interference in the angular domain. A 3D GIDM (Gaussian interference distribution model) is proposed, and the corresponding interference APD (angular power density) is given. Then, the closed-form expressions of some key spatial statistics of interference are derived based on the interference APD. Finally, the closed-form expressions of the probability density function and spatial correlation function of SIR (signal–to–interference ratio) are derived based on the 3D multipath APD model and spatial statistics of the Rice channel. Simulation analysis shows that 3D spatial angular directions have significant effect on these spatial statistics of interference and the spatial correlation function of SIR. The results provide useful insight on the analysis and design of the interference-limited networks.

1. Introduction

C-V2X (cellular vehicle-to-everything) communication is essential for redefining transportation by providing real-time, highly reliable, and actionable information flows to enable safety, mobility, and environmental applications [1]. C-V2X communication networks comprise a huge number of spatially scattered nodes, and the high-density spectrum spatial multiplexing is the main method to improve spectrum utilization. However, C-V2X communication networks become interference-limited due to the frequency reuse. An appropriate interference model is of great significance to the design and performance enhancement of C-V2X communication systems, and it is also a research hotspot in the field of wireless communication. However, the complex and changeable signal propagation environment brings severe challenges to the study of the 3D model and spatial characteristics of C-V2X communication interference.
A lot of work has been done in various ways in the field of wireless communication interference modeling and analysis. The interference modeling and characteristic research results of cellular network, ad hoc network, and D2D (Device–to–Device) communication have important reference significance for the corresponding research of C-V2X communication. A circular interference model that aggregates given interferer topologies to power profiles along circles was proposed in [2] for heterogeneous cellular networks. A statistical model for quantifying the aggregate interference in heterogeneous cellular networks was proposed, and its impact on system performance was evaluated in [3]. Interference from static cellular cooperation using the nearest neighbor model was analyzed in [4]. Interference was modeled using hierarchical spatial clustering of terrain and user density maps in [5]. Interference prediction was investigated in [6] in a mobile ad hoc network with a finite number of nodes by proposing and using a general-order linear model for node mobility. HetNet interference was modeled and analyzed using Poisson cluster processes in [7]. Efficient simulation laws for the tail of the interference in a simple wireless ad hoc network model were provided in [8] considering node locations distributed according to a Poisson point process and various classes of light-tailed fading distributions. Average received interference power of D2D communication in the cellular network was analyzed in [9]. The interference graph construction for D2D communications in cellular networks was studied, and an interference graph construction protocol was proposed in [10]. Mathematical characteristics of the uplink interference region in D2D communications underlaying cellular networks were analyzed in [11]. A cross-layer model was established for D2D communications underlaying cellular network to characterize the realistic interference scenario in [12]. The impact of co-channel interference on the performance of V2V (vehicle–to–vehicle) communication was investigated based on 5.9 GHz real-world measurements conducted at intersections under LoS (line-of-sight) and NLoS (non-line-of-sight) conditions in [13]. The effect of hidden terminal interference on safety-critical traffic in vehicular ad hoc networks was investigated in [14]. A theoretical interference model of V2V communications at an intersection that uses a transmission power control method was proposed in [15]. The influence of interference in an inter-vehicular communication system was studied, and the interfering links were modeled using the multiple scattering radio channel in [16]. Interference signal approximations for system performance predictions and the impact on the IEEE 802.11p system in terms of bit and packet error probability and packet delays were analyzed in [17]. Methodology, experiment configuration, and final results from a set of measurements executed in order to validate the interference model assumptions adopted when simulating wireless networks with NS-3 were described in [18]. A novel network description method named the interference azimuth spectrum was proposed to determine the fundamental second-order small-scale fading statistics in the 2D angular domain in [19] and [20].
Since the heights of the vehicle nodes are different, and the signal is also scattered by the scatterers in the 3D space, the C-V2X communication interference is distributed in the 3D space. However, the model and characteristics of C-V2X communication interference have rarely been investigated so far in the 3D angular domain according to the above analysis. Therefore, this paper mainly investigates spatial characteristics of C-V2X communication interference in the 3D angular domain and obtains the effect of 3D spatial angular directions (i.e., azimuth angle and elevation angle) on spatial statistics of interference and SIR (signal–to–interference ratio). This paper strives to alleviate the current lack of analytical studies by proposing a 3D GIDM (Gaussian interference distribution model), giving the corresponding interference APD (angular power density) and deriving the closed-form expressions of some key spatial statistics of interference and SIR in the 3D angular domain.
The remainder of the paper is organized as follows. Section 2 proposes a 3D GIDM and gives the corresponding interference APD. Section 3 derives some key spatial statistics of interference. Section 4 derives the PDF (probability density function) and SCF (spatial correlation function) of the SIR. Discussions on obtained simulation results are provided in Section 5. Finally, conclusions are mentioned in Section 6.

2. 3D GIDM and Interference APD

2.1. 3D GIDM

Figure 1 shows the 3D GIDM to describe the interference power distribution at the receiver R X . The 3D GIDM is based on the following assumptions and settings:
(1)
the closer the interference signal to the receiving end, the stronger the influence, so the incoming interference power is assumed to obey the Gaussian distribution centered on the receiver in the 3D spatial distance [21];
(2)
since large-scale C-V2X communication networks have many nodes, it is assumed that the number of interference waves tends to infinity, and the multipath channel modeling theory can be used to study the spatial statistical characteristics of interference [19,20].
In Figure 1, I is the interferer if the interferer (base station or other terminal node) and R X have the LoS path, and I is the last scatterer before the arrival of the interference signal if the interferer and R X have no LoS path. α [ 0 , 2 π ] and β [ 0 , π / 2 ] are the azimuth angle and elevation angle of the interference wave, respectively. r is the distance from I to R X . σ is the standard deviation of the Gaussian distribution.
“Small probability events” generally refer to events with a probability of less than 5%, while Gaussian distribution variables with a mean of zero fall below [ 3 σ , 3 σ ] with a probability of less than 0.3%. In practical problems, it is often assumed that the corresponding event does not occur. Basically, the interval [ 3 σ , 3 σ ] can be regarded as the actual possible range of the Gaussian random variable. This is called the 3 σ principle of the Gaussian distribution. Therefore, according to the 3 σ principle of the Gaussian distribution, the maximum effective value of r in Figure 1 is 3 σ .

2.2. Interference APD

Since the incoming interference power is assumed to obey the Gaussian distribution in 3D space, its expression is
f G ( x , y , z ) = A 0 2 π σ exp ( x 2 + y 2 + z 2 2 σ 2 )
where f G ( x , y , z ) is the spatial Gaussian distribution in Cartesian coordinates, and z 0 . A 0 is a constant used to ensure that f G ( x , y , z ) is integrated into 1 in the 3D space. In actual communication, the distributions of interferers (or interference scatterers) in the horizontal and vertical directions are different. Generally, the smaller the elevation angle, the more interferers (or interference scatterers) exist, and the greater the impact the interference causes. Therefore, the standard deviation of the Gaussian distribution needs to be controlled separately in the horizontal and vertical directions. This paper assumes that the standard deviation of the Gaussian distribution in the horizontal direction is σ α , the standard deviation in the vertical direction is σ β , and σ α > σ β . Then the maximum effective value of r in the vertical direction is 3 σ β , and the maximum effective value of r in the horizontal direction is 3 σ α . The 3D standard deviation σ can be expressed as
σ = σ α 2 σ β 2 σ β 2 cos 2 β + σ α 2 sin 2 β
Equation (1) in Cartesian coordinates can be converted to an equation in spherical coordinates by the Jacobian formula [22]:
P G ( r , α , β ) = f G ( x , y , z ) | J ( x , y , z ) | | x = r cos β cos α y = r cos β sin α z = r sin β = f G ( x , y , z ) r 2 cos β = A 0 r 2 2 π σ exp ( r 2 2 σ 2 ) cos β
Then the interference APD can be calculated as
P G ( α , β ) = 0 p G ( r , α , β ) d r = A 0 cos β 2 π σ 0 r 2 exp ( r 2 2 σ 2 ) d r
because
x 2 exp ( a x 2 ) d x = x 2 a exp ( a x 2 ) + a π 4 a 2 erf ( a x )
where erf ( · ) is the error function, erf ( x ) = 2 π 0 x exp ( t 2 ) d t . Let x = r and a = 1 / 2 σ 2 , using the conclusion of Equation (5), and the closed expression of interference APD can be written as
P G ( α , β ) = A 0 cos β 2 π σ [ σ 2 r exp ( r 2 2 σ 2 ) + 2 π σ 3 2 erf ( 2 2 σ r ) ] | | 0 = A 0 σ 2 cos β 2 = A 0 σ α 2 σ β 2 cos β 2 σ β 2 cos 2 β + 2 σ α 2 sin 2 β
Therefore, the interference APD in the vertical direction is
P G ( β ) = 0 2 π A 0 σ α 2 σ β 2 cos β 2 σ β 2 cos 2 β + 2 σ α 2 sin 2 β d α = A 0 π σ α 2 σ β 2 cos β σ β 2 cos 2 β + σ α 2 sin 2 β
The interference APD in the horizontal direction is
P G ( α ) = 0 π 2 A 0 σ α 2 σ β 2 cos β 2 σ β 2 cos 2 β + 2 σ α 2 sin 2 β d β = 1 2 π 0 π 2 P G ( β ) d β = 1 2 π
because the integral of the interference APD is 1, therefore
0 π 2 P G ( β ) d β = 0 π 2 A 0 π σ α 2 σ β 2 cos β σ β 2 cos 2 β + σ α 2 sin 2 β d β = A 0 π σ α 2 σ β [ ( σ α A 3 ) A 2 arctan ( σ β / A 1 ) ( σ α + A 3 ) A 1 arctan ( σ β / A 2 ) ] A 1 A 2 A 3 = 1
where A 1 = σ β 2 + 2 σ α ( σ α σ α 2 σ β 2 ) , A 2 = σ β 2 + 2 σ α ( σ α + σ α 2 σ β 2 ) , A 3 = σ α 2 σ β 2 . So, the expression of A 0 is A 1 A 2 A 3 / π σ α 2 σ β ( σ α A 3 ) A 2 arctan ( σ β / A 1 ) ( σ α + A 3 ) A 1 arctan ( σ β / A 2 ) .

3. Spatial Statistics of Interference

3.1. 3D Shape Factors

According to the definition in [23], the l-th degree and m-th order unnormalized complex spherical harmonic coefficient S l m for the interference APD P G ( α , β ) can be derived as
{ S 0 0 = A 0 π 2 σ α 2 σ β 2 ( σ β + σ α ) ,   S 1 0 = A 0 π σ α 2 σ β 2 ( σ β 2 σ α 2 ) 3 / 2 { σ β 2 σ α 2 + σ α arctan [ 2 σ α 2 σ β 2 2 σ α σ β 2 σ α 2 ] σ α arctan [ σ α σ β 2 σ α 2 ] } S 1 1 = 0 ,   S 2 0 = A 0 π 2 σ α 2 σ β ( σ β 2 σ α ) 12 ( σ β + σ α ) 2 ,   S 2 1 = 0 ,   S 2 2 = 0
Therefore, the 3D shape factors defined in [24] can be derived as
Λ = 1 ( S 1 0 ) 2 + | S 1 1 | 2 ( S 0 0 ) 2 = 1 ( S 1 0 ) 2 ( S 0 0 ) 2
ξ = 3 2 S 2 0 S 0 0 ( S 1 0 ) 2 + 1 2 | S 1 1 | 2 ( S 0 0 ) 2 ( S 1 0 ) 2 | S 1 1 | 2 = 3 S 2 0 S 0 0 / 2 ( S 1 0 ) 2 ( S 0 0 ) 2 ( S 1 0 ) 2
χ = 2 | S 2 1 S 0 0 S 1 0 S 1 1 | ( S 0 0 ) 2 ( S 1 0 ) 2 | S 1 1 | 2 = 0
γ = | S 2 2 S 0 0 ( S 1 1 ) 2 | ( S 0 0 ) 2 ( S 1 0 ) 2 | S 1 1 | 2 = 0
α β 45 max = arg { S 2 1 S 0 0 S 1 0 S 1 1 } = arg { 0 }
α β 0 max = 1 2 arg { S 2 2 S 0 0 ( S 1 1 ) 2 } = 1 2 arg { 0 }
where Λ , ξ , χ , γ , α β 45 max and α β 0 max are the 3D angular spread, elevational constriction, 45° inclined constriction, azimuthal constriction, azimuth of maximum fading at 45° elevation, and azimuth of maximum fading at zero elevation, respectively. Λ [ 0 , 1 ] is a measure of how interference power concentrates about a single direction; ξ [ 0.5 , 1 ] shows how much the interference APD is concentrated on a single elevational cone or along two paths at the same azimuth and opposite elevations; χ = 0 shows the interference APD is horizontal or vertical mirror symmetrical; γ = 0 shows two-path interference APD has no significant difference in the azimuth direction; Equations (15) and (16) indicate that α β 45 max and α β 0 max do not exist.

3.2. Interference Fading Rate Variance

Since the mean derivative of a stationary process is zero, the variance of the fading rate of a received complex voltage is the simplest statistic that measures the fading rate, and it can be derived as [24].
σ I 2 ( α , β ) = 4 π 2 Λ 2 P I 3 λ 2 { 1 + 3 2 [ ξ ( 2 sin 2 β 2 3 ) + χ sin 2 β cos ( α α β 45 max ) + γ cos 2 β cos 2 ( α α β 0 max ) ] } = 4 π 2 Λ 2 P I [ 1 + ξ ( 3 sin 2 β 1 ) ] 3 λ 2
where P I is the total received interference power in a local volume, and λ is the carrier wavelength. According to Equation (17), the maximum and minimum values of σ I 2 ( α , β ) are respectively,
σ I , max 2 = max { 4 π 2 Λ 2 P I [ 1 + ξ ( 3 sin 2 β 1 ) ] 3 λ 2 } = 4 π 2 Λ 2 P I ( 1 + 2 ξ ) 3 λ 2 , β = π 2
σ I , min 2 = min { 4 π 2 Λ 2 P I [ 1 + ξ ( 3 sin 2 β 1 ) ] 3 λ 2 } = 4 π 2 Λ 2 P I ( 1 ξ ) 3 λ 2 , β = 0
For the mobile receiver of C-V2X communication networks, it is often convenient to measure the variance of the fading rate in terms of the change over time. According to [25], the mean-squared time rate-of-change σ V 2 ( α , β ) of a received complex voltage is equal to σ I 2 ( α , β ) times the square of the moving speed v R of the receiver. Since this paper mainly studies the variation of interference and SIR, regardless of the specific amplitude, the moving speed v R of the receiver is set to 1, and σ I 2 ( α , β ) is used instead of σ V 2 ( α , β ) in the derivation of expressions below.

3.3. Interference Level Crossing Rate and Average Fade Duration

Each interference fading signal is assumed to obey Rayleigh distribution or Rice distribution. According to [26,27], the sum of multiple Rayleigh or Rice distributions can approximate a Nakagami-m distribution. Therefore, the total interference signal at the receiver R X obeys the Nakagami-m distribution. Ref. [20] has already derived the expressions of 2D spatial LCR (level crossing rate) and AFD (average fade duration) for a Nakagami-m fading signal, and the fading rate variance is only closely related to the space dimension. Therefore, the 3D LCR is
N I ( ρ N , α , β ) = σ I ( α , β ) π P I m m 1 / 2 Γ ( m ) ρ N 2 m 1 exp ( m ρ N 2 ) = 2 Λ ρ N 2 m 1 3 π [ 1 + ξ ( 3 sin 2 β 1 ) ] 3 λ exp ( m ρ N 2 )
where ρ N = R l / P T is the normalized threshold level, R l is the envelope threshold, m is the Nakagami-m shape factor, and Γ ( · ) is the Gamma function. By definition, the 3D AFD can be derived as
l ¯ I ( ρ N , α , β ) = 1 N I ( ρ N , α , β ) 0 R P m ( r ) d r = 1 N I ( ρ N , α , β ) Γ ( m , m ρ N 2 ) Γ ( m ) = 3 π Γ ( m , m ρ N 2 ) λ 2 π Λ ρ N 2 m 1 Γ ( m ) 1 + ξ ( 3 sin 2 β 1 ) exp ( m ρ N 2 )
where P m ( r ) is the PDF of the Nakagami-m distribution, and Γ ( · , · ) is the incomplete Gamma function.

3.4. Interference Spatial Correlation and Coherence Distance

The interference envelope SCF can be defined as
R I ( r c , α , β ) = Ε [ r ( p 0 ) r ( p 0 + r c p ) ] [ Ε ( r ) ] 2 Ε ( r 2 ) [ Ε ( r ) ] 2
where E ( ) is the statistical expectation operator, p 0 is any initial position, p is the unit direction vector, and r c is the spatial distance. Equation (22) can be approximated by a Mclaurin expansion [28]:
R I ( r c , α , β ) = 1 + n = 1 ( 1 ) n r c 2 n ( 2 n ) ! Ε [ ( d n r d r c n ) 2 ] Ε ( r 2 ) [ Ε ( r ) ] 2 = 1 Ε [ ( d r d r c ) 2 ] 2 { Ε ( r 2 ) [ Ε ( r ) ] 2 } r c 2 + = 1 σ I 2 ( α , β ) 2 V m ( r ) r c 2 +
where V m ( r ) is the variance of Nakagami-m distribution.
At the same time, SCF can also be approximated as an exponential function and its Mclaurin expansion is as follows:
R I ( r c , α , β ) exp [ a ( r c λ ) 2 ] 1 a ( r c λ ) 2 +
Let Equation (23) = Equation (24), and interference SCF can be obtained
R I ( r c , α , β ) exp [ σ I 2 ( α , β ) 2 V m ( r ) r c 2 ] = exp { 2 π 2 Λ 2 m Γ 2 ( m ) [ 1 + ξ ( 3 sin 2 β 1 ) ] 3 λ 2 [ m Γ 2 ( m ) Γ 2 ( m + 1 2 ) ] r c 2 }
If the SCF equals exp ( 1 ) for the spatial distance r c , the distance is defined as the coherence distance D I ,
D I = λ 6 m 2 m π Λ Γ ( m ) m Γ 2 ( m ) Γ 2 ( m + 1 2 ) 1 + ξ ( 3 sin 2 β 1 )

4. Spatial Statistics of SIR

4.1. Probability Density Function of SIR

The useful signal r S at the receiver R X is assumed to obey the Rice distribution, and its average power is expressed as P S . The total interference signal r I at R X obeys the Nakagami-m distribution, and its average power is expressed as P I . Therefore, the SIR can be written as η = r S 2 / r I 2 .
According to [26,29], the square of the Nakagami-m distribution variable r I 2 (i.e., the instantaneous power of r I ) obeys the Gamma distribution G ( a 1 , b 1 ) ,
f ( x ; a 1 , b 1 ) = b 1 ( b 1 x ) b 1 1 exp ( b 1 x ) Γ ( a 1 )
where a 1 = m , b 1 = m / P I . The mean of the Gamma distribution is a 1 / b 1 and the variance is a 1 / b 1 2 .
The Rice distribution can be approximated as the Nakagami-m distribution when K = m s 2 m s / ( m s m s 2 m s ) [26], and K is the Rice K-factor defined as the ratio of the power contributions by the LoS path to the remaining multipath, m s is the Nakagami-m shape factor for Rice distribution. Therefore, the square of the Rice distribution variable also approximately obeys the Gamma distribution G ( a 2 , b 2 ) , and a 2 = m s = ( K + 1 ) 2 / ( 2 K + 1 ) , b 2 = m s / P S = ( K + 1 ) 2 / [ P S ( 2 K + 1 ) ] . The mean and the variance of G ( a 2 , b 2 ) are a 2 / b 2 and a 2 / b 2 2 , respectively.
The ratio of the two Gamma distribution variables G ( a 1 , b 1 ) / G ( a 2 , b 2 ) obeys the compound Gamma distribution β ( a 1 , a 2 , 1 , b 2 / b 1 ) which is a special case of the general beta basic distribution [30], and its expression is
f ( x ; a 3 , b 3 , 1 , c 3 ) = β ( a 3 , b 3 , 1 , c 3 ) = c 3 b 3 x a 3 1 Β ( a 3 , b 3 ) ( c 3 + x ) a 3 + b 3
where Β ( x , y ) = 0 1 t x 1 ( 1 t ) y 1 d t is beta function. Therefore, the PDF of SIR η can be derived as
P ( η ) = ( b 2 / b 1 ) a 2 η a 1 1 Β ( a 1 , a 2 ) ( b 2 / b 1 + η ) a 1 + a 2 = b 1 a 1 b 2 a 2 η a 1 1 Β ( a 1 , a 2 ) ( b 2 + b 1 η ) a 1 + a 2 , { a 1 = m b 1 = m / P I a 2 = ( K + 1 ) 2 / ( 2 K + 1 ) b 2 = ( K + 1 ) 2 / [ P S ( 2 K + 1 ) ]
The k-th order moment, mean and variance of SIR η are respectively [30]
Ε [ η k ] = b 2 k Β ( a 2 k , a 1 + k ) b 1 k Β ( a 1 , a 2 ) , a 2 > k
Ε [ η ] = a 1 b 2 b 1 ( a 2 1 ) , a 2 > 1
V a r [ η ] = a 1 b 2 2 ( a 1 + a 2 1 ) b 1 2 ( a 2 1 ) 2 ( a 2 2 ) , a 2 > 2

4.2. Spatial Correlation Function of SIR

Similar to the interference, the SCF of SIR can be defined as
R ( η c , α , β ) = Ε [ η ( p 0 ) η ( p 0 + η c p ) ] [ Ε ( η ) ] 2 Ε ( η 2 ) [ Ε ( η ) ] 2
where η c is the spatial distance.
According to Equations (30)–(32),
Ε ( η 2 ) [ Ε ( η ) ] 2 = V a r [ η ] = a 1 b 2 2 ( a 1 + a 2 1 ) b 1 2 ( a 2 1 ) 2 ( a 2 2 ) , a 2 > 2
[ Ε ( η ) ] 2 = a 1 2 b 2 2 b 1 2 ( a 2 1 ) 2 , a 2 > 1
For the convenience of derivation, let Φ η ( η c ) = η ( p 0 ) η ( p 0 + η c p ) , Φ S ( η c ) = r S ( p 0 ) r S ( p 0 + η c p ) and Φ I ( η c ) = r I ( p 0 ) r I ( p 0 + η c p ) , then
Ε [ Φ η ( η c ) ] = Ε [ η ( p 0 ) η ( p 0 + η c p ) ] = Ε [ r S 2 ( p 0 ) r I 2 ( p 0 ) r S 2 ( p 0 + η c p ) r I 2 ( p 0 + η c p ) ] = Ε [ Φ S 2 ( η c ) Φ I 2 ( η c ) ]
According to [28], Equation (36) can be approximated as
Ε [ Φ η ( η c ) ] = Ε [ Φ S 2 ( η c ) Φ I 2 ( η c ) ] Ε 2 [ Φ S ( η c ) ] Ε 2 [ Φ I ( η c ) ] + 1 2 [ 2 Φ η ( η c ) Φ S 2 ( η c ) σ Φ S ( η c ) 2 + 2 Φ η ( η c ) Φ I 2 ( η c ) σ Φ I ( η c ) 2 ] = Ε 2 [ Φ S ( η c ) ] Ε 2 [ Φ I ( η c ) ] + 1 2 { 2 Ε 2 [ Φ I ( η c ) ] σ Φ S ( η c ) 2 + 6 Ε 2 [ Φ S ( η c ) ] Ε 4 [ Φ I ( η c ) ] σ Φ I ( η c ) 2 } Ε 2 [ Φ S ( η c ) ] Ε 2 [ Φ I ( η c ) ] + 1 Ε 2 [ Φ I ( η c ) ] σ Φ S ( 0 ) 2 + 3 Ε 2 [ Φ S ( η c ) ] Ε 4 [ Φ I ( η c ) ] σ Φ I ( 0 ) 2
where σ Φ S ( η c ) 2 and σ Φ I ( η c ) 2 are variances of Φ S ( η c ) and Φ I ( η c ) , respectively. Φ S ( 0 ) and Φ I ( 0 ) obey the Gamma distribution, so
σ Φ S ( 0 ) 2 = Ε [ Φ S 2 ( 0 ) ] Ε 2 [ Φ S ( 0 ) ] = a 2 b 2 2
σ Φ I ( 0 ) 2 = Ε [ Φ I 2 ( 0 ) ] Ε 2 [ Φ I ( 0 ) ] = a 1 b 1 2
According to Equation (22),
Ε [ Φ I ( η c ) ] = Ε [ r I ( p 0 ) r I ( p 0 + η c p ) ] = R I ( η c , α , β ) { Ε ( r I 2 ) [ Ε ( r I ) ] 2 } + [ Ε ( r I ) ] 2 = R I ( η c , α , β ) V m + E m 2
where V m and E m are the variance and mean of Nakagami-m distribution, respectively. Similarly, let R S ( η c , α , β ) denote SCF of the useful signal, whose expression is as shown in Equation (16) in our previous paper [31], then
Ε [ Φ S ( η c ) ] = Ε [ r S ( p 0 ) r S ( p 0 + η c p ) ] = R S ( η c , α , β ) { Ε ( r S 2 ) [ Ε ( r S ) ] 2 } + [ Ε ( r S ) ] 2 = R S ( η c , α , β ) V R i c + E R i c 2
where V R i c and E R i c are the variance and mean of Rice distribution, respectively.
Substituting Equations (38)–(41) into Equation (37),
Ε [ Φ η ( η c ) ] [ R S ( η c , α , β ) V R i c + E R i c 2 ] 2 [ R I ( η c , α , β ) V m + E m 2 ] 2 + a 2 b 2 2 [ R I ( η c , α , β ) V m + E m 2 ] 2 + 3 a 1 [ R S ( η c , α , β ) V R i c + E R i c 2 ] 2 b 1 2 [ R I ( η c , α , β ) V m + E m 2 ] 4 = b 2 2 [ R S ( η c , α , β ) V R i c + E R i c 2 ] 2 + a 2 b 2 2 [ R I ( η c , α , β ) V m + E m 2 ] 2 + 3 a 1 [ R S ( η c , α , β ) V R i c + E R i c 2 ] 2 b 1 2 [ R I ( η c , α , β ) V m + E m 2 ] 4
Substituting Equations (34), (35) and (42) into Equation (33), the SCF of SIR η can be obtained
R ( η c , α , β ) = b 1 2 ( a 2 1 ) 2 ( a 2 2 ) a 1 b 2 2 ( a 1 + a 2 1 ) { b 2 2 [ R S ( η c , α , β ) V R i c + E R i c 2 ] 2 + a 2 b 2 2 [ R I ( η c , α , β ) V m + E m 2 ] 2 + 3 a 1 [ R S ( η c , α , β ) V R i c + E R i c 2 ] 2 b 1 2 [ R I ( η c , α , β ) V m + E m 2 ] 4 a 1 2 b 2 2 b 1 2 ( a 2 1 ) 2 }

5. Results and Discussions

This section illustrates and analyses the interference APD, some key spatial statistics of interference, PDF and SCF of SIR described in Section 2, Section 3 and Section 4 through the MATLAB simulation. Unless indicated otherwise, the values of the simulation parameters used to obtain the curves are set to λ = 0.051 m (5.9 GHz), m = 2 , and σ β = 5   m for the general urban C-V2X communication scenarios. For ease of analysis, P I = P S = 1 , v R = 1   m / s .

5.1. Interference APD

Figure 2 demonstrates interference APD for different standard deviations σ α of the Gaussian distribution in the horizontal direction and an interference APD in the scenario where the incoming interference power is assumed to obey the uniform distribution in the 3D space. In the actual urban C-V2X communication, the effective interferer’s height is limited, but it may be far from the receiver R X in the horizontal direction. Therefore, the horizontal interference has a greater influence on the statistical characteristics of the interference, and Figure 2 only considers the impact of different σ α . In the above equations, the unit of the elevation angle is radians. For ease of description and understanding, the unit of the elevation angle in all simulations in Section 5 is degrees, so the interference APD in Figure 2 has some values greater than 1. If the incoming interference power is assumed to obey the uniform distribution in the 3D space, the interference APD can be derived as cos β / 2 π . As shown in Figure 2, the larger σ α is, the more interference power is concentrated in the small elevation angle range. When σ α 400   m , the interference APD curves tend to coincide. According to the 3 σ principle of the Gaussian distribution, this indicates that the interference effects outside the horizontal distance 3 σ α 1200   m are small due to the large path loss and shadow fading. For different σ α , the interference APD curves decrease sharply with the increase of the elevation angle. When the elevation angle is 15°, the interference APD values are 0.072, 0.037, 0.018, and 0.014, respectively; that is, the curves are almost reduced to zero. The interference APD curve under the assumption of uniform distribution decreases extremely slowly with the increase of the elevation angle, which cannot well characterize the law that the influence of the interference in the vertical direction is small in the actual C-V2X communication. The above spatial statistics of interference can provide reference for the design of the antenna arrays and the beamforming technology, so that the antenna and the beam avoid the main interference to improve the C-V2X communication performance.

5.2. Interference Fading Rate Variance

Figure 3 demonstrates the root-mean square (RMS) normalized fading rate variance of interference for different σ α and elevation angles. For facilitating comparison with the 2D GIDM, Figure 3 uses the RMS normalized fading rate variance, which is the root-mean square of the “3D fading rate variance/2D fading rate variance”. The 2D fading rate variance can be derived as σ I 2 D 2 = 2 π 2 P I / λ 2 . As shown in Figure 3, as σ α increases, the RMS normalized fading rate variance curves almost coincide in the range of (0°, 60°) of elevation. The curves differ in the range of (60°, 90°) of elevation, and moreover, the larger σ α is, the smaller the RMS normalized fading rate variance is at the same elevation angle, but when σ α 400   m , the curves tend to be consistent. For example, the average relative deviation of RMS normalized fading rate variances at σ α = 400   m and σ α = 500   m has been reduced to 0.8%, which also shows that the interference effects outside the horizontal distance 3 σ α 1200   m are small. The RMS normalized fading rate variance in Figure 3 is less than or equal to one, which indicates that the received interference power fades slower for 3D GIDM than 2D GIDM. For wireless communication, the smaller interference, and the higher fading speed of interference are good for improving system performance. As shown in Figure 2 and Figure 3, the interference APD is almost zero and the RMS normalized fading rate variances are large near the elevation angle of 15°. Therefore, the elevation angle of 15° is the proper direction of antenna and beam pointing.

5.3. Interference LCR

Figure 4 demonstrates the interference LCR for different normalized thresholds and elevation angles. As shown in Figure 4, the normalized threshold has a decisive effect on the interference LCR, and the elevation angle also has a significant effect on the interference LCR. As the normalized threshold increases, the interference LCR also increases rapidly. As the elevation angle increases, the interference LCR gradually decreases. When the elevation angle approaches 90°, the interference LCR decreases to almost zero, indicating that the interference fading speed becomes slower and the fluctuation also becomes smaller when approaching the vertical direction, and Figure 3 also reflects this characteristic of interference. According to Equation (21), the interference AFD is inversely proportional to the interference LCR at the same normalized interference threshold. Therefore, the variation law of the interference AFD for different normalized thresholds and elevation angles can also be obtained from Figure 4.

5.4. Interference SCF

Figure 5 demonstrates the interference SCF for different spatial distances and elevation angles and whether the spatial distance or the elevation angle has a decisive effect on the interference SCF. In the range of (0°, 60°) of elevation, the effect of the spatial distance is dominant, and when the spatial distance is r c 0.5 λ , the interference SCF approaches zero. However, in the range of (60°, 90°) of elevation, the effect of the elevation angle is dominant, and when the elevation angle approaches 90°, even if the spatial distance is r c = λ , the interference SCF is still large. The SCF is an important indicator to measure the spatial selectivity of interference, which can guide the optimal design of the separation distance of antenna arrays. According to Figure 5, the separation distance of antenna arrays should be changed with the change of the elevation angle, that is, a small separation distance is required at low elevation angles to reduce the influence of interference correlation, and a large separation distance is required at high elevation angles. Similarly, the interference SCF is almost close to zero when r c 0.25 λ near the elevation angle of 15°.

5.5. PDF of SIR

Figure 6 demonstrates the PDF of SIR for different Nakagami-m shape factor m and the Rice K-factor K. According to the Shannon formula, the channel capacity mainly depends on the SIR, which is an important indicator for measuring channel quality. Large-capacity transmission requires a high SIR. Since this paper mainly studies the variation of interference and SIR, regardless of the specific amplitude, the average power of the interference signal and the useful signal is set to P I = P S = 1 . Therefore, the peaks of the PDF curves in Figure 6 are around 0 dB, and the peak position moves slightly with the change of m and K. The larger the K is, the smaller the overall fluctuation and attenuation of the multipath signal. Similarly, the larger the m is, the smaller the overall fluctuation and attenuation of the total interference at the receiver. As shown in Figure 6, the PDF of SIR is directly affected by m and K. The larger the m or K is, the more concentrated the PDF curve is, and the closer the peaks of the PDF curves are to the 0 dB. Moreover, the effects of m and K have a superposition effect.

5.6. SCF of SIR

Figure 7 demonstrates the SCF of SIR for different azimuth angles, elevation angles, and spatial distances. Since Equation (43) requires a 2 = ( K + 1 ) 2 / ( 2 K + 1 ) > 2 , this section sets K = 4 . In addition, the azimuth and elevation angles of the LoS path of the useful signal are set to α 0 = π / 4 and β 0 = π / 6 , respectively. As shown in Figure 7, the separation distance or the elevation angle has a decisive effect on the SCF of SIR, while the azimuth angle has less of an effect. When the separation distance η c < 0.1 λ , the SCF of SIR is small, and it is basically not affected by the 3D direction. As the separation distance increases, the SCF of SIR also increases, which indicates that the SCF of the useful signal decays more slowly than the SCF of the interference signal. The SCF of SIR is large in the range of (0°, 60°) of elevation angle, however, it is small in the range of (60°, 90°) of elevation angle, which indicates that the SCFs of the useful signal and the interference signal also have different fading speeds in different ranges of elevation angle. Figure 6 and Figure 7 describe the variation of the spatial statistics of SIR in the 3D space, which provides theoretical support for the design and performance evaluation of the large-capacity C-V2X transmission technology.

6. Conclusions

In this paper, a 3D GIDM has been proposed to analyze the spatial characteristics of C-V2X communication interference in the angular domain. Then, the closed-form expressions of interference APD, some key spatial statistics of interference, PDF, and SCF of SIR have been derived based on the 3D multipath shape factors theory and spatial statistics of the Rice channel. The interference APD and the key spatial statistics of interference can well characterize the different strengths of co-channel interference of C-V2X communication in different 3D angular directions. Simulation analysis shows that 3D spatial angular directions have a significant effect on these spatial statistics of interference and the SCF of SIR. The results provide useful insight on the design and performance evaluation of the large-capacity transmission technology for the interference- limited networks. In future research, the following work will be very important: (1) the frequency reuse, channel access, and modulation can also deeply affect the performance of C-V2X communication interference and these factors need to be considered in further researches; (2) the difference between 3D model and 2D model of interference in a quantitative manner needs to be investigated by actual signal measurements.

Author Contributions

Formal analysis, S.H. and Y.L.; funding acquisition, D.D., X.J. and Y.T.; investigation, D.D.; methodology, D.D. and X.Z.; software, X.W., Y.T., X.Z., S.H. and Y.L.; Validation, X.W.; writing—original draft, D.D.; writing—review and editing, X.J.

Funding

This research was funded by Fundamental Research Funds for the Central Universities (No. 2018CDXYTX0009, No. 2019CDXYTX0023), the National Natural Science Foundation of China (No. 61571069, No. 61701054, No. 61771080, and No. 61703062), and the Scientific and Technological Research Program of Chongqing Municipal Education Commission (Grant No. KJ1601204).

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Du, D.; Zeng, X.; Jian, X.; Yang, F.; Sun, M. 3-D V2V MIMO Channel Modeling in Different Roadway Scenarios with Moving Scatterers. Prog. Electromagn. Res. 2018, 64, 43–54. [Google Scholar]
  2. Taranetz, M.; Rupp, M. A circular interference model for heterogeneous cellular networks. IEEE Trans. Wirel. Commun. 2016, 15, 1432–1444. [Google Scholar] [CrossRef]
  3. Zhang, T.; Lu, A.; Chen, Y.; Chai, K.K. Aggregate interference statistical modeling and user outage analysis of heterogeneous cellular networks. In Proceedings of the IEEE International Conference on Communications (ICC), Sydney, Australia, 10–14 June 2014; pp. 1260–1265. [Google Scholar]
  4. Giovanidis, A.; Corrales, L.D.A.; Decreusefondf, L. Analyzing interference from static cellular cooperation using the nearest neighbour model. In Proceedings of the 13th International Symposium on Modeling and Optimization in Mobile, Ad Hoc, and Wireless Networks (WiOpt), Mumbai, India, 25–29 May 2015; pp. 576–583. [Google Scholar]
  5. Koufos, K.; Jantti, R. Interference Modelling Using Hierarchical Spatial Clustering of Terrain and User Density Map. In Proceedings of the IEEE 79th Vehicular Technology Conference (VTC Spring), Seoul, Korea, 18–21 May 2014; pp. 1–5. [Google Scholar]
  6. Cong, Y.; Zhou, X.; Kennedy, R.A. Interference prediction in mobile ad hoc networks with a general mobility model. IEEE Trans. Wirel. Commun. 2015, 14, 4277–4290. [Google Scholar] [CrossRef]
  7. Chun, Y.J.; Hasna, M.O.; Ghrayeb, A. Modeling and analysis of HetNet interference using Poisson Cluster Processes. In Proceedings of the IEEE 25th Annual International Symposium on Personal, Indoor, and Mobile Radio Communication (PIMRC), Washington, DC, USA, 2–5 September 2014; pp. 681–686. [Google Scholar]
  8. Torrisi, G.L.; Leonardi, E. Simulating the tail of the interference in a Poisson network model. IEEE Trans. Inform. Theory 2013, 59, 1773–1787. [Google Scholar] [CrossRef]
  9. Huan, C.; Liu, C.; Zheng, W.; Fu, Y. Average received interference power analysis of D2D communication in the cellular network. In Proceedings of the IEEE International Conference on Ubiquitous Wireless Broadband (ICUWB), Nanjing, China, 16–19 October 2016; pp. 1–4. [Google Scholar]
  10. Zhang, Y.; Zheng, J.; Lu, P.S.; Chen, S. Interference Graph Construction for Cellular D2D Communications. IEEE Trans. Veh. Technol. 2016, 66, 3293–3305. [Google Scholar] [CrossRef]
  11. Li, X.; Zhang, W.; Zhang, H.; Li, W. Mathematical characteristics analysis of uplink interference region in D2D communications underlaying cellular networks. In Proceedings of the IEEE Seventh International Conference on Ubiquitous and Future Networks (ICUFN), Sapporo, Japan, 7–10 July 2015; pp. 557–561. [Google Scholar]
  12. Lu, H.; Wang, Y.; Chen, Y.; Liu, K.J. Interference model and analysis on device-to-device cellular coexist networks. In Proceedings of the IEEE Global Conference on Signal and Information Processing (GlobalSIP), Orlando, FL, USA, 14–16 December 2015; pp. 1086–1090. [Google Scholar]
  13. Tchouankem, H.; Lorenzen, T. Measurement-based evaluation of interference in Vehicular Ad-Hoc Networks at urban intersections. In Proceedings of the IEEE International Conference on Communication Workshop (ICCW), London, UK, 8–12 June 2015; pp. 2381–2386. [Google Scholar]
  14. Bastani, S.; Landfeldt, B. The Effect of Hidden Terminal Interference on Safety-Critical Traffic in Vehicular Ad Hoc Networks. In Proceedings of the 6th ACM Symposium on Development and Analysis of Intelligent Vehicular Networks and Applications, Valletta, Malta, 13–17 November 2016; pp. 75–82. [Google Scholar]
  15. Kimura, T.; Saito, H. Theoretical Interference Analysis of Inter-vehicular Communication at Intersection with Power Control. In Proceedings of the 19th ACM International Conference on Modeling, Analysis and Simulation of Wireless and Mobile Systems, Valletta, Malta, 13–17 November 2016; pp. 3–10. [Google Scholar]
  16. Bithas, P.S.; Efthymoglou, G.P.; Kanatas, A.G. Intervehicular communication systems under co-channel interference and outdated channel estimates. In Proceedings of the IEEE International Conference on Communications (ICC), Kuala Lumpur, Malaysia, 23–27 May 2016; pp. 1–6. [Google Scholar]
  17. Tengstrand, S.O.; Fors, K.; Stenumgaard, P.; Wiklundh, K. Jamming and interference vulnerability of IEEE 802.11 p. In Proceedings of the IEEE International Symposium on Electromagnetic Compatibility (EMC Europe), Gothenburg, Sweden, 1–4 September 2014; pp. 533–538. [Google Scholar]
  18. Fuxjaeger, P.; Ruehrup, S. Validation of the NS-3 Interference Model for IEEE802. 11 Networks. In Proceedings of the IEEE 8th IFIP Wireless and Mobile Networking Conference (WMNC), Munich, Germany, 5–7 October 2015; pp. 216–222. [Google Scholar]
  19. Chen, Y.; Mucchi, L.; Wang, R. Angular spectrum and second order statistics of interference in wireless networks. In Proceedings of the IEEE International Conference on Communication, Networks and Satellite (COMNETSAT), Yogyakarta, Indonesia, 3–4 December 2013; pp. 41–45. [Google Scholar]
  20. Chen, Y.; Mucchi, L.; Wang, R.; Huang, K. Modeling Network Interference in the Angular Domain: Interference Azimuth Spectrum. IEEE Trans. Commun. 2014, 62, 2107–2120. [Google Scholar] [CrossRef] [Green Version]
  21. Jiang, H.; Zhang, Z.; Dang, J.; Wu, L. A Novel 3D Massive MIMO Channel Model for Vehicle-to-Vehicle Communication Environments. IEEE Trans. Commun. 2018, 66, 79–90. [Google Scholar] [CrossRef]
  22. Nawaz, S.J.; Qureshi, B.H.; Khan, N.M. A Generalized 3-D Scattering Model for a Macrocell Environment with a Directional Antenna at the BS. IEEE Trans. Veh. Technol. 2010, 59, 3193–3204. [Google Scholar] [CrossRef]
  23. Bell, W.W. Special Functions for Scientists and Engineers; Dover Publications: New York, NY, USA, 2004. [Google Scholar]
  24. Valchev, D.G.; Brady, D. Three-dimensional multipath shape factors for spatial modeling of wireless channels. IEEE Trans. Wirel. Commun. 2009, 8, 5542–5551. [Google Scholar] [CrossRef]
  25. Durgin, G.D.; Rappaport, T.S. Theory of multipath shape factors for small-scale fading wireless channels. IEEE Trans. Antennas Propag. 2000, 48, 682–693. [Google Scholar] [CrossRef]
  26. Nakagami, M. The m-distribution-A general formula of intensity distribution of rapid fading. In Statistical Method of Radio Propagation; Pergamon: Oxford, UK, 1960; pp. 3–36. [Google Scholar]
  27. Babich, F.; Lombardi, G. General Nakagami approximation for sum of Ricean interferers. Electron. Lett. 1998, 34, 23–24. [Google Scholar] [CrossRef]
  28. Papoulis, A.; Pillai, S.U. Probability, Random Variables, and Stochastic Processes; Tata McGraw-Hill Education: New York, NY, USA, 2002. [Google Scholar]
  29. Karagiannidis, G.; Georgopoulos, C.; Kotsopoulos, S. Outage probability analysis for a Nakagami signal in L Nakagami interferers. Eur. Trans. Telecommun. 2001, 12, 145–150. [Google Scholar] [CrossRef]
  30. Dubey, S.D. Compound gamma, beta and F distributions. Metrika 1970, 16, 27–31. [Google Scholar] [CrossRef]
  31. Du, D.; Zeng, X.; Jian, X.; Yu, F.; Miao, L. Analysis of Three-Dimensional Spatial Selectivity for Rician Channel. Radioengineering 2018, 27, 249–255. [Google Scholar] [CrossRef]
Figure 1. The 3D Gaussian interference distribution model (GIDM) of interference.
Figure 1. The 3D Gaussian interference distribution model (GIDM) of interference.
Electronics 08 00718 g001
Figure 2. The interference angular power density (APD).
Figure 2. The interference angular power density (APD).
Electronics 08 00718 g002
Figure 3. The root-mean square normalized fading rate variance of interference.
Figure 3. The root-mean square normalized fading rate variance of interference.
Electronics 08 00718 g003
Figure 4. The interference level crossing rate (LCR).
Figure 4. The interference level crossing rate (LCR).
Electronics 08 00718 g004
Figure 5. The interference spatial correlation function (SCF).
Figure 5. The interference spatial correlation function (SCF).
Electronics 08 00718 g005
Figure 6. The probability density function (PDF) of signal–to–interference ratio (SIR).
Figure 6. The probability density function (PDF) of signal–to–interference ratio (SIR).
Electronics 08 00718 g006
Figure 7. The SCF of SIR.
Figure 7. The SCF of SIR.
Electronics 08 00718 g007

Share and Cite

MDPI and ACS Style

Du, D.; Jian, X.; Wu, X.; Tan, Y.; Zeng, X.; Huang, S.; Li, Y. 3D Spatial Characteristics of C-V2X Communication Interference. Electronics 2019, 8, 718. https://doi.org/10.3390/electronics8060718

AMA Style

Du D, Jian X, Wu X, Tan Y, Zeng X, Huang S, Li Y. 3D Spatial Characteristics of C-V2X Communication Interference. Electronics. 2019; 8(6):718. https://doi.org/10.3390/electronics8060718

Chicago/Turabian Style

Du, Derong, Xin Jian, Xuegang Wu, Yong Tan, Xiaoping Zeng, Shijian Huang, and Yadong Li. 2019. "3D Spatial Characteristics of C-V2X Communication Interference" Electronics 8, no. 6: 718. https://doi.org/10.3390/electronics8060718

APA Style

Du, D., Jian, X., Wu, X., Tan, Y., Zeng, X., Huang, S., & Li, Y. (2019). 3D Spatial Characteristics of C-V2X Communication Interference. Electronics, 8(6), 718. https://doi.org/10.3390/electronics8060718

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop