Next Article in Journal
Ultrasound Standing Wave-Based Cell-to-liquid Separation for Measuring Viscosity and Aggregation of Blood Sample
Previous Article in Journal
Vibration Control of Scanning Electron Microscopes with Experimental Approaches for Performance Enhancement
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Orbit Angular Momentum MIMO with Mode Selection for UAV-Assisted A2G Networks

1
School of Communication and Information Engineering, Chongqing University of Posts and Telecommunications, Chongqing 400065, China
2
Department of Electronic and Electrical Engineering, The University of Sheffield, Sheffield S10 2TN, UK
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Sensors 2020, 20(8), 2289; https://doi.org/10.3390/s20082289
Submission received: 24 February 2020 / Revised: 15 April 2020 / Accepted: 16 April 2020 / Published: 17 April 2020
(This article belongs to the Section Communications)

Abstract

:
As an emerging solution for line-of-sight (LOS) wireless communications, in air-to-ground (A2G) channels, the unmanned aerial vehicle (UAV), and allowing the dynamic and flexible network deployments enables the supplement or/and replacement of the terrestrial base stations (BSs). However, in conventional multiple-input-multiple-output (MIMO) systems, high-speed communications are significantly limited by channel crosstalks and spectrum scarcities. An orbit angular momentum (OAM) wireless network, allowing co-existence of multiple physical channels within the same frequency band, offers new degrees of freedom to address this dilemma. In this paper, we investigate the UAV-based A2G radio vortex wireless networks and study its channel model. Then we propose a branch and bound search-based mode selection (BBS-MS) scheme, which uses the spatial distribution characteristics of vortex beams to optimize the spectrum efficiency (SE). Theoretical derivations and numerical results demonstrate that our developed BBS-MS scheme can obtain the optimal performance, which outperforms conventional OAM-based MIMO systems. Also, it possesses a lower complexity compared with exhaustive searches.

Graphical Abstract

1. Introduction

Orbit angular momentum (OAM) is a physical property of electromagnetic waves, which can be used to generate multiple orthogonal vortex beams [1]. Recently, OAM multiplexing has attracted much attention as a novel approach for improving the spectral efficiency (SE) of multiple-input-multiple-output (MIMO) wireless communication systems [2,3,4]. The first experimental test has demonstrated that OAM can be applied to wireless communications [2]. Then, a radio vortex-multiple-input multiple-output (RV-MIMO) system with high SE is introduced in [3], where multiple spiral phase plate (SPP) antennas are used for OAM generations. Also, the feasibility that the uniform circular array (UCA) antenna can be used to creat twist beams was theoretically and practically verified [4]. Furthermore, OAM techniques with the characteristic of orthogonal transmissions are widely used in optical communications [5,6], line-of-sight (LOS) communications [7,8,9], and radar imaging [10,11].
In the case of natural or human-made disasters, the traditional terrestrial communication networks will be damaged and communication links can be destroyed. The flexible and low complexity air base stations (BSs) can be rapidly established to provide in-time communication supports in this situation [12,13,14]. The unmanned aerial vehicle with its flexibility and low complexity was considered to be a candidate of the air BS. The unmanned aerial vehicle (UAV) with its flexibility and low complexity is viewed as a candidate of the air BS. Besides, the UAV techniques are widely used in atmospheric sciences [15,16,17,18,19,20,21,22,23]. The authors in [15] presented the information about the use of rotary-wing droned to sample near-surface fields. In [16], a UAV-based system was proposed to measure the profiles and turbulence in the boundary layer. Moreover, some measurements relative to the wind and turbulence in the boundary layer were carried out in [17,18,19]. The researchers in [20] proposed a method for detecting atmospheric Lagrangian coherent structures using a fixed-wing drone. Then, a coordinated drone and ground-based weather measurement system was introduced to predict Lagrangian structures [21]. Besides, the UAV can be used to validate the weather modeling [22] and monitoring trace tropospheric gases [23].
Meahwhile, in air-to-ground (A2G) wireless networks, UAV-based BS (UAV-BS) plays a fundamental role in serving temporal hotspots to improve the quality of service [24,25]. In addition, UAV-BS also plays a crucial role in enabling reliable and high-speed wireless backhaul connectivity for the traditional terrestrial networks [26]. Temporal hotspots such as big tournaments, open-air concerts, and huge gatherings, where users tend to be fixed or less mobile, often yield severe network congestion. Fortunately, in such scenarios, reliable and ideal LOS communication links can be established to improve the quality and capacity of communications by optimally placing drones to avoid obstacles. Also, OAM-based MIMO can significantly improve the SE of A2G network due to the dominant LOS channels in A2G communication links. Generally, both of the scenarios need timely nodes for assistance but still suffer from accompanying interference and lack of frequency bands. According to the physical property of Laguerre Gauss (LG) beams, the energy distribution of OAM waves is in the shape of a doughnut [1,9] as shown in Figure 1. As the mode absolute value increases, the angle between maximum radiation direction and propagation direction becomes extensive. It is notable that the circles shown in Figure 1 denote the range of the power distribution which coresponding to the OAM main lobes. Thus, any terrestrial users equipped with a general OAM receiving antenna cannot receive the whole OAM signal ultimately.
To improve the SE of the A2G wireless communication links, we first model the UAV-based LOS A2G channel that makes mode division multiplexing (MDM) transfer possible. Then we continue with a design of the mode selection scheme, named by branch and bound search-based mode selection (BBS-MS), touching upon the improvements in performance and the reductions of computational complexities. Also, we derive capacity expressions for the mode selection system, and the numerical simulations demonstrate the superiority of BBS-MS scheme over conventional MDM systems is obtained.

2. System Model and Problem Formulation

In this section, we consider an air-to-ground-radio vortex (A2G-RV) wireless communication system, as shown in Figure 1. At the transmitter, the drone is equipped with an N t -element UCA antenna, where T X n ( 1 n N t ) is n th transmit antenna element, the parameter α is the transmitting rotation angle, and 2 π / N t is the angle between adjacent antenna elements. However, at the receiver, the ground user is equipped with a 2-element OAM receiver ( R X 1 and R X 2 ), where φ denotes the receiving azimuth angle, and the angle between two receiving antenna elements is denoted by δ . Please note that the OAM signal can be precisely captured once the constant B < π | l max | , where l max is the maximum value of the OAM topological charge [27]. The parameter d n 1 ( d n 2 ) denotes the distances between n th ( n [ 1 , N ] ) transmitting element and the receiving antenna R X 1 ( R X 2 ) . The altitudes of the UAV-based OAM transmitter with the radius of R and the ground user are h t and h r , respectively.
Also, A2G-RV wireless channels experience the dominant LOS channels even in dense urban environments, which is different from the conventional terrestrial communications [14,24,25]. Therefore, the LOS A2G channel matrix with multi-path, denoted by H C N r × N t , which can be expressed as
H = H LOS
where H LOS is the LOS channel matrix from the transmitter to the receiver.

2.1. The LOS A2G Channel

To facilitate channel modeling, as depicted in Figure 2, we employ the cylindrical coordinate system ( r , θ , ξ ) to express the position of the antenna element, in which r, θ , and ξ denote the polar axis, angular coordinate and longitudinal axis, respectively. Furthermore, the position of the centre of the transmit UCA T X n is given by ( 0 , 0 , 0 ) , and the positions of the receive antenna R X 1 and R X 2 are expressed by ( r , φ , h t h r ) and ( r , φ + δ , h t h r ) , respectively, where r is the radial distance of the receive antenna pair corresponding to the propagation distance as h t h r in the longitudinal axis. Thus, the direct distances from n th transmitting element to the receiving elements R X 1 and R X 2 can be calculated as follows
d LOS ( n , 1 ) = a n 1 2 + b n 1 2 + d ˜ n 1 2
d LOS ( n , 2 ) = a n 2 2 + b n 2 2 + d ˜ n 2 2
with d ˜ n 1 = d ˜ n 2 = h t h r and
a n m = { r cos φ R cos ( ϕ n + α ) , m = 1 r cos ( φ + δ ) R cos ( ϕ n + α ) , m = 2
b n m = { r sin φ R sin ( ϕ n + α ) , m = 1 r sin ( φ + δ ) R sin ( ϕ n + α ) , m = 2
where ϕ n = 2 π n / N t denotes the basic azimuth angle for the transmitting UCA. Thus, the d LOS ( n , 1 ) and d LOS ( n , 2 ) can be rewritten as follows
d LOS ( n , 1 ) ( a ) h t 2 + h r 2 + r 2 + R 2 h t h r h t 2 + h r 2 + r 2 + R 2 r R cos ( ϕ n + α φ ) h t 2 + h r 2 + r 2 + R 2
d LOS ( n , 2 ) ( a ) h t 2 + h r 2 + r 2 + R 2 h t h r h t 2 + h r 2 + r 2 + R 2 r R cos ( ϕ n + α φ δ ) h t 2 + h r 2 + r 2 + R 2
where (a) denotes the total sum of squares approximation with the example of a 2 2 b a b a .
In free-space LOS scenarios, there is an only direct path at both ends of the transceiver, and the propagation of electromagnetic waves leads to amplitude corruption and periodic phase change of the transmitted signal. Therefore, the channel responses between n th antenna element of the transmitting UCA and m th antenna element of the receiving UCA is denoted by [3]
h LOS ( n , m ) = β λ 4 π d LOS ( n , m ) e j k d LOS ( n , m )
where β is the attenuation corresponding to the pattern of the antenna element, d LOS ( n , m ) is the direct distance from n th transmitting element to m th receiving element, λ denotes the wavelength of the carrier wave, and k = 2 π / λ denotes the wavenumber. Consequently, in the A2G radio vortex communication system, the channel responses between n th antenna element in the transmitting UCA and the receiving antenna R X 1 and R X 2 can be rewritten as
h LOS ( n , 1 ) ( b ) β λ 4 π d ˜ LOS e j k d ˜ LOS e j k [ h t h r + r R cos ( ϕ n + α φ ) ] d ˜ LOS
h LOS ( n , 2 ) ( b ) β λ 4 π d ˜ LOS e j k d ˜ LOS e j k [ h t h r + r R cos ( ϕ n + α φ δ ) ] d ˜ LOS
with d ˜ LOS = h t 2 + h r 2 + r 2 + R 2 , where (b) denotes the approximate condition that d LOS ( n , 1 ) = d LOS ( n , 2 ) h t 2 + h r 2 + r 2 + R 2 for the amplitude term.

2.2. A2G-RV Channel

In the UCA-based multiple mode transmitter, the vortex beams can be generated by feeding the same input symbol s l to each antenna element in transitting UCA with a successive phase different. Consequently, the transmitted signal at n th antenna element, denoted by x n , is given as
x n = 1 N t l L ρ l s l e j [ ϕ n + α ] l
where L = l 0 , l 1 , , l i , , l N t 1 is the set of the OAM modes with the value of l i [ N t / 2 , N t / 2 ) , ρ l is the transmitted power allocated to l th OAM mode, and ϕ n + α is the transmitted azimuthal angle. Therefore, the matrix expression of the complex baseband signal mode of the A2G-RV scheme would be described as
y = H x + w
where y C M × 1 is the received vector, H C M × N t is the A2G channel matrix, w C M × 1 is the vector of the received additive complex white Gaussian noise (AWGN) with zero-mean and variance σ 2 , and the x C N t × 1 = F Ps , where F C N t × N t is the discrete Fourier transform (DFT) matrix, P = diag [ ρ 0 ρ 1 ρ N t 1 ] is the transmitting power allocation matrix, and s = [ s 0 s 1 s N t 1 ] T is the data symbols vector.
In the proposed A2G-RV system, the channel responses between n th transmit antenna, and the receive antenna R X 1 and R X 2 have the equations as
h LOS l ( n , 1 ) = β ˜ l , 1 λ 4 π d ˜ LOS e j k d ˜ LOS e j ( φ ϕ n α ) l e j k [ h t h r + r R cos ( ϕ n + α φ ) ] d ˜ LOS
h LOS l ( n , 2 ) = β ˜ l , 2 λ 4 π d ˜ LOS e j k d ˜ LOS e j ( φ + δ ϕ n α ) l e j k [ h t h r + r R cos ( ϕ n + α φ δ ) ] d ˜ LOS
where β ˜ l , 1 and β ˜ l , 2 denote the attenuation outside of the circle regions of OAM beams. The attenuation coefficients satisfy β ˜ 1 = β ˜ 2 = β when the two-element receive antenna locates in the circle regions of OAM beams. In this case, the fading of OAM signals obeys the Friis Law. However, in a general A2G-RV communication system, the receiver (ground user) could be outside of the OAM circle regions, where the attenuation coefficients β ˜ l , 1 and β ˜ l , 2 depend on the location of the receiver. Furthermore, in the proposed system, the attenuation coefficients satisfy β ˜ l , 1 = β ˜ l , 2 = β ˜ l due to the equal radial distance r of each receive element. The following will present the mathematical derivation of β ˜ l .
The Laguerre-Gaussian (LG) beams, the most widely used example of OAM research, were considered to describe the radio OAM beams [27]. By solving the partial differential equation u l p I l r r = 0 , where u l p is the expression of the LG beams [27], the radius of the maximum beam intensity can be given by
r max , l z = l 2 w l z = w l | l | 2 ( 1 + ( z z R ) 2 )
Given the A2G-RV mode selection wireless communication system, the response between n th transmit antenna and the receive antenna R X m , and the response between n th transmit antenna and the receive antenna R X m can be denoted by [27]
u l ( n , m ) = h LOS l ( n , m ) x ^
u l ( n , m ) = h LOS l ( n , m ) x ^
where u l is the expression of LG beams [8,27], x ^ is the unit pulse input, R X m and R X m donate the receiver are inside and outside of the OAM regions, respectively. Based on this assumption, the positions of the receive antenna R X m and R X m are expressed by ( r max , l ( z ) , φ , h t h r ) and ( r n m , φ , h t h r ) . Furthermore, and the responses can be derived by using Equations (13) and (14) as follows
u l ( n , m ) = β ˜ l ( 2 r max , l ( z ) ) | l | w l ( z ) | l | + 1 e r max , l 2 ( z ) w l 2 ( z ) e j k r max , l 2 ( z ) 2 R l ( z ) e j ( | l | + 1 ) e j l φ m
u l ( n , m ) = β ˜ l ( 2 r n m ) | l | w l ( z ) | l | + 1 e r n m 2 w l 2 ( z ) e j k r n m 2 2 R l ( z ) e j ( | l | + 1 ) e j l φ m
Combining Equations (13)–(19), the attenuation β ˜ l , 1 can be calculated as
β ˜ l , m = β d LOS ( n , m ) d LOS ( n , m ) ( r n m r max , l ( z ) ) | l | e r n m 2 r max , l 2 ( z ) w l 2 ( z ) e j k ( r n m 2 r max , l 2 ( z ) ) 2 R l ( z ) e j k [ d LOS ( n , m ) d LOS ( n , m ) ]
where r n m is the radial distance between n th transmit antenna and the receive antenna R x m . When the receiver locates inside of the circle region corresponding to OAM mode l, the attenuation coefficientsatisfies β ˜ l , m = β . Substitute Equation (20) into Equation (13), the channel responses can be rewritten as
h LOS l ( n , m ) = β λ 4 π d ˜ LOS ( r n m r max , l ( z ) ) | l | e r n m 2 r max , l 2 ( z ) w l 2 ( z ) e j k ( r n m 2 r max , l 2 ( z ) ) 2 R l ( z ) e j k [ h t h r + r R cos ( ϕ n + α φ m ) ] d ˜ LOS e j k d ˜ LOS e j ( φ m ϕ n α ) l
Combining Equations (11), (12) and (21), the A2G-RV channel matrix, denoted by H RV which can be expressed as follows
H RV = h 11 h 1 N t h 21 h 2 N t h M 1 h M N t 1 1 1 e j Φ 2 N t l N t 1 1 e j Φ M N t l N t 1
where the entry h m n can be replaced by h LOS ( n , m ) , and Φ m n = ϕ n + α φ m is the azimuth angle.

3. Capacity and Mode Selection

3.1. Exact Upper Bounds

In this paper, we suppose that the channel gains related to the transmit modes is known at the transmitter, and the independent Gaussian information symbols are transmitted. Furthermore, we consider a time division duplexing (TDD) MIMO system, where only uplink channels need to be estimated due to channel reciprocity. The capacity can be written as [28]
C = log 2 [ det ( I N t + P T N t H RV H RV H ) ]
where P T is the total transmitted power, and H RV = [ h 1 , , h n , , h N t ] is the A2G-RV channel matrix. Please note that the uniform power allocation policy is considered in this paper for simplicity, where the transmitting power allocated to each effective OAM mode is ρ = P T L with 1 L N t . With the channel gains at the transmit side, we can sort the channel gains | | h n | | 2 in descending order, i.e., | | h 1 | | 2 | | h 2 | | 2 | | h N t | | 2 . Then, the transmitter can obtain the best index subset S L = 1 , , L , which determines the best L out of N t transmit OAM modes from the ordered index set S = 1 , , N t by using a rate-limited return channel [28]. Finally, the optimal transmit channel matrix relative to the optimal subset S L , denoted by H ˜ RV = [ h ˜ 1 , , h ˜ L ] , can be obtained.
Based on the description above, the mode selections for maximizing the channel capacity can be built into an optimization problem as follows
H ˜ opt = arg max H ˜ RV H { log 2 [ det ( I L + ρ H ˜ RV H H ˜ RV ) ] }
Subject to ρ L P T ρ 0 1 L N t
where H is the sets of total C ( N t , L ) possible submatrices for the mode selection system. Thus, the capacity of the mode selection system is derived as
C opt = log 2 [ det ( I L + ρ H ˜ opt H H ˜ opt ) ]
Benefit from the orthogonality between different OAM subchannels, the upper bound on the capacity for the A2G-RV MIMO system with mode selection can be defined as [28]
C A 2 G RV Upper = l = 1 L log 2 ( 1 + ρ γ l )
with γ l = | h l | 2 σ l 2 , where h l is the A2G-RV channel amplitude gain, which was derived in Equation (21) for LOS path, and σ l 2 is the variance of the received complex AWGN which corresponding to the OAM mode l.
Given an ideal A2G-RV communication system, we can accurately derive the upper bound on the mode selection channel capacity for the LOS path. In this case, the expression of the capacity equation is a closed-form. However, in practical applications, the channel capacity may be variable as well as the random channel matrix H RV . To the best of our knowledge, there is no statistics model of RV-MIMO channels so far. In this paper, we focus on the superiority of the mode selection scheme over the conventional MDM scheme in the ideal OAM channel rather than in the statistical model.

3.2. The Branch and Bound Search for Mode Selection

In the proposed mode selection system, the number of the selected transmit OAM modes should satisfy the constraint that L min ( N t , M ) . The OAM mode selection problem given in Equation (24), where several best OAM subchannels are chosen for transmission, which is similar to the antenna selection problem [29]. Both the mode selection and antenna selection problem can be classified as a feature subset selection problem. Based on this similarity, we can directly apply many optimized search algorithms used by the antenna selection to the mode selection. In terms of algorithm selection, exhaustive search and greed search are the two most popular algorithms for antenna selection and beam selection [29,30]. On the one hand, exhaustive search provides the optimal performance as well as tremendous computational cost, especially with a massive MIMO system [30]. On the other hand, a greedy search that relies on its ultra-low complexity has been widely applied in beam/antenna selection systems. Therefore, the choice of algorithm mainly depends on the system configurations. To address this problem, we proposed a fast and optimal searching algorithm, named by branch and bound search-based mode selection (BBS-MS), which can also be used in a large-scale system. Based on the optimization problem in Equation (24), furthermore, we employ the minimum loss channel capacity to construct the criterion function for the BBS-MS algorithm, and the definition of the minimum loss channel capacity will be given in the following content.
It is well known that the global optimal BBS scheme is very efficient without exhaustive search once the corresponding monotonic feature function can be constructed [31]. Based on the searching direction, the BBS scheme can be classified into upward and downward searching algorithms [32]. In a UCA-based short-distance OAM wireless communication system, we certainly expect to choose as many as possible OAM subchannels for transmission. For mode selection in this paper, the downward BBS algorithm is selected, where the maximum value of a monotonically decreasing criterion function needs to be found out. Specifically, the procedure of designing a criterion function can be divided into three steps as follows

3.2.1. Selection of the Sub-Criterion Function

We introduce a loss capacity function, which composed of the removed N t L subset to be the sub-criterion function due to the search principle of the downwards BBS algorithm. Please note that the downwards BBS algorithm starts with a full set and then discards one OAM mode resource per step from this full set until there is only an L feature subset left.
Let H ˜ n C M × n ¯ be the sub-channel matrix of H RV relative to n ¯ = N t n best OAM modes by discarding n worst OAM modes. Thus, the loss capacity function is given by C loss ( n ) = C log 2 [ det ( I M + ρ ˜ H ˜ n H ˜ n H ) ] with ρ ˜ = P T / N t . Suppose that in the n + 1 step of OAM mode discard, the sub-channel h ˜ j corresponding to j th column of H RV is discarded, which yields to a minimum contribution to the capacity function. The new sub-channel matrix H ˜ n + 1 C M × n + 1 ¯ can be obtained by dropping h ˜ j from H ˜ n . Then, it has
C loss ( n + 1 ) = log 2 [ C det ( I M + ρ ˜ H ˜ n + 1 H ˜ n + 1 H ) ]
with H ˜ n + 1 H ˜ n + 1 H = H ˜ n H ˜ n H h ˜ j h ˜ j H and C = det ( I N t + ρ ˜ H RV H RV H ) . By applying the Sherman-Morrison formula [28] to Equation (28), we obtain the transformation
C loss ( n + 1 ) = log 2 [ C det ( I M + ρ ˜ H ˜ n H ˜ n H ) ] log 2 [ det ( I M + ρ ˜ H ˜ n + 1 H ˜ n + 1 H ) det ( I M + ρ ˜ H ˜ n H ˜ n H ) ] = C loss ( n ) log 2 [ det ( I M ρ ˜ h ˜ j h ˜ j H I M + ρ ˜ H ˜ n H ˜ n H ) ] = C loss ( n ) log 2 ( 1 ρ ˜ h RV , j H B n h RV , j ) = C loss ( n ) log 2 ( 1 ρ ˜ ϵ j , n ) = C loss ( n ) j , n
where B n = ( I M + H ˜ n H ˜ n H ) 1 and ϵ j , n = h ˜ j B n h ˜ j H are notations, suppose that the mode index J n , which can minimize C loss n + 1 in Equation (29), was found in n + 1 step. Then, the updated matrix inverse B n + 1 can be expressed as
B n + 1 = ( I M + H ˜ n + 1 H ˜ n + 1 H ) 1 = B n + b n + 1 H b n + 1
where b n + 1 = B n h ˜ J n ρ ˜ 1 + ϵ J n .

3.2.2. Construction of the Criterion Function

The loss capacity function in Equation (29) is monotonically increasing due to j , n < 0 . In this case, we introduce an offset to reverse the monotonicity of loss capacity function C loss ( n ) which is similar to the scheme in [32], and then the criterion function can be defined as
C ˜ loss ( n ) = s = 0 n 1 Λ s C loss ( n )
where Λ s = log 2 ( 1 ρ ˜ γ s 2 I M + ρ ˜ H ˜ s + 1 H ˜ s + 1 H ) is the offset, and γ s = min j I s | | h ˜ j | | F for the mode selections, where I s is the mode set of all of the candidate mode sources in s th step for discarding the worst mode. In the following, we need to prove that C ˜ loss ( n ) is a monotonically decreasing function of n.

3.2.3. Proof of Monotonicity

Combining Equations (29) and (31), the recursive form of C ˜ loss ( n ) is given by
C ˜ loss ( n + 1 ) = C ˜ loss ( n ) + Λ n + j , n
with n = 0 , 1 , , N t L 1 . From Equation (32), we need to prove that Λ n + j , n 0 , which ensures C ˜ loss ( n + 1 ) is a monotonically decreasing function of n. In other words, we need to verify that Λ n | j , n | .
Lemma 1.
Given j , n = log 2 ( 1 + ρ ˜ ϵ j , n ) in Equation (29) and Λ n in Equation (31). Then we have Λ n | j , n | , n = 0 , 1 , , N t L 1 .
Proof. 
Thanks to a similar justification in [32]. From Equations (29) and (30), we have
ϵ j , n = h ˜ j h ˜ j H I M + ρ ˜ H RV H RV H h ˜ j ( 1 I M + ρ ˜ H RV H RV H + s = 1 n b s H b s ) h ˜ j H
where n = 0 and n = 1 , 2 , , N t L 1 are the conditions for equality in Equation (33).
When n = 0 , we have
| j , 0 | = | log 2 ( 1 ρ ˜ ϵ j , 0 ) | = | log 2 ( 1 ρ ˜ h ˜ j h ˜ j H I M + ρ ˜ H RV H RV H ) | ( a ) | log 2 ( 1 ρ ˜ γ s 2 I M + ρ ˜ H ˜ n + 1 H ˜ n + 1 H ) | = Λ 0
where the “=” in the condition (a) holds when the candidate mode index j I s in 0th discard step.
When n = 1 , 2 , , N t L 1 . The matrix s = 1 n b s H b s is a positive semi-definite matrix since the gram matrix b s H b s is positive semi-definite. Then
B n 1 I M + ρ ˜ H RV H RV H = s = 1 n b s H b s 0
Thus, we have
| j , n | = | log 2 ( 1 ρ ˜ ϵ j , 0 ) | = | log 2 [ 1 ρ ˜ h ˜ j ( 1 I M + ρ ˜ H RV H RV H + s = 1 n b s H b s ) h ˜ j H ] | ( b ) | log 2 ( 1 ρ ˜ γ s 2 I M + ρ ˜ H ˜ n + 1 H ˜ n + 1 H ) | = Λ n
where the “=” in the condition (b) holds when the candidate mode index j I s in n th discard step. □
Please note that the criterion function C ˜ loss ( n ) is a monotonically decreasing function of n, which equals to the offset s = 0 n 1 Λ s minus the sub-criterion function C loss ( n ) according to Equation (31). Moreover, the selected mode subset yielding a minimum value of C loss ( n ) will produce the globally best value of C ˜ loss ( n ) since the offset s = 0 n 1 Λ s is a constant value in each channel realization and independent of the selected mode subset. The optimal mode subset leading to the maximum value of the criterion function C ˜ loss ( n ) can be found by using the proposed branch and bound algorithm.

3.2.4. Searching Procedure

Thanks to the pseudocode of BAB-search-based beam selection in [24], the proposed BBS-MS scheme is shown in Algorithm 1. It is notable that if the initial value of B is close to the optimal C ˜ loss opt ( N t L ) , lesser nodes need be evaluated than that with a smaller B, which leads to a dependable B. For simplicity, we choose the initial value as B = in this paper.
To better understand the downward BBS-MS algorithm, Figure 3 shows an example of searching trees by discarding 3 out of 8 OAM modes. In addition, the main steps of the BBS-based searching can be summarized as follows
Step 0: Initialize the parameters. Set S = 1 , , N t , S = 1 , 2 , , N t L , S loss = , N 0 , 1 = L + 1 , J 0 = 1 , B = , and n = 1 ;
Step 1: Initialize S . Compute C ˜ loss ( n ) for all OAM subchannels corresponding to the OAM mode set S and record the values of all C ˜ loss ( n ) into the set j in descending order.
Step 3: Check the lower-bound. If C ˜ loss ( n ) < B , return the entry J to the set S and shift to Step 4. In addition if n = N t L 1 , shift to Step 5. Or else, let n = n + 1 and then shift to Step 1.
Step 4: Backtrack. Let n = n + 1 . If n = N t L , terminate the BBS-MS algorithm. Or else, return the entry J to the set S and shift to Step 2.
Algorithm 1 The BBS-MS scheme
Require:
1:Initialization parameters: B = 1 I M + ρ ˜ H RV H RV H , B = , C ˜ loss ( i ) = M , n = 0 , J 0 = 1 , S = 1 , 2 , , N t , S = 1 , 2 , , N t L , and the discared mode index set S loss =
2: ν j = | | h ˜ j | | 2 , j S
3: ϵ j = h ˜ j H B h ˜ j , j S
4: γ s = min j I s ν j , s S
5: Λ s = log 2 ( 1 ρ ˜ γ s 2 I M + ρ ˜ H ˜ s + 1 H ˜ s + 1 H ) , s S
6:if n = N t L 1  then
7: C ˜ loss ( n + 1 ) j = C ˜ loss ( n ) + j , n + Λ N t L 1 , j I n
8:if min s I N t L 1 C ˜ loss ( n + 1 ) s < B then
9:   J = arg max s I N t L 1 C ˜ loss ( n + 1 ) s
10:  update S loss = S loss J
11:   B = max s I N t L 1 C ˜ loss ( n + 1 ) s
12:end if
13:else
14: C ˜ loss ( n + 1 ) j = C ˜ loss ( n ) + j , n + Λ n , j I n
15: sort C ˜ loss ( n + 1 ) j , j I n in a desend order to get an ordered index set j
16: B ¯ = B , ν ¯ j = ν j , j S
17:for i = 1 : | I n | do
18:   J = [ j ] i
19:  if C ˜ loss ( n + 1 ) j > B then
20:    K = [ j ] i + 1 , [ j ] i + 2 , , [ j ] N t n + 1
21:   update the index set S loss = S loss J
22:    b n + 1 = B n h ˜ J n ρ ˜ 1 + ϵ J , and B = B ¯ + b n + 1 H b n + 1
23:    C ˜ loss ( n ) = C ˜ loss ( n + 1 ) J
24:    ϵ s = h ˜ s H ( ϵ ˜ s + b n + 1 H b n + 1 ) h ˜ s , s K
25:    s = log 2 ( 1 ρ ˜ ϵ s ) , s K
26:    n = n + 1 , go to line 6
27:  else
28:   break the loop
29:  end if
30:end for
31:end if
32:Output: the final mode index set S selected = S \ S loss
Step 5: Update the lower-bound. Let B = max s I N t L 1 C ˜ loss ( n + 1 ) s , where C ˜ loss ( n + 1 ) is obtained by discarding OAM modes 1 , , N t L and then save it as the temporary optimal solution. Return the entry J to the set S loss and shift to Step 4.

4. Numerical Results

In this section, we evaluate the performance of the proposed BBS-MS, the greedy search, the exhaustive search, and the conventional MDM schemes. First, we analyze the ability of OAM waves corresponding to different mode number for information transmission when considering different locations of receivers. Then, we evaluate the spectrum efficiency of the proposed BBS-MS scheme considering the different numbers of selected transmitting modes. The default simulation parameters are summarized in Table 1. Notably, the standard value of the available L is default as 2 in this paper due to the same performance between the two modes that are negative to each other.
Figure 4 delineates the amplitude gains relative to the different locations of receivers. As illustrated in Figure 4a, within a relatively short transmission range, all OAM waves almost possess the same amplitude gain with the horizontal distance increasing. However, in Figure 4b, we can observe that the OAM waves with different mode numbers consecutively obtain the optimal amplitude gain in the process of the receiver away from the centre point. With the altitude of UAV increasing, both in Figure 4c,d, this phenomenon is more visible, and the amplitude gains of the high-order OAM waves are close to zero. The result also validates that the fundamental field distribution property of vortex electromagnetic waves, i.e., concentric zone pattern.
In Figure 5, we evaluate the performance of the traditional MDM system and other mode selected-based schemes, including the proposed BBS-MS, the greedy search, and the exhaustive search algorithms. Notably, both the MDM and the conventional MIMO systems possess the same performance in the LOS channel. Since the sub-optimal greedy search can achieve the most optimal performance in such a small-scale mode selection system, we can observe that the three mode selected algorithms obtain the same performance, for both N t = 8 and N t = 4 cases. Also, in the similar conventional antenna selection systems, the same conclusion was drawn in [12]. Since the properties of divergence and circular radiant field of OAM waves, the received intensities emitted by different OAM signals depend on the location of the receiver. In other words, In a medium or long distance OAM wireless network, only a part of OAM modes are suitable for transmission. Thus, the full mode MDM system obtains the worst performance.
The evaluation of the different numbers of the selected OAM modes is presented in Figure 6, where we set h t = 20 m, SNR = 20 dB, and R = 10 λ and employ the uniform power allocation scheme for comparison. It can be observed that the smaller the number of selected modes L, the better SE enhancement the mode selection system can obtain over the MDM systems. As is well know, the optimal system performance can be achieved by adopting a water-filling power allocation strategy. Due to different transmission gains of OAM waves in wireless channels, there is a most energetic OAM beam, which enables the mode selection scheme to obtain the optimal system performance. That is to say, a single-mode with the best transmission gain can be selected to maximize the SE.

5. Conclusions

In this paper, by reviewing the major field distribution of electromagnetic vortex waves, we explored the mode selected-based radio vortex communication system in the A2G scenario. To maximize the SE of the system, the BBS-MS was proposed for OAM mode selections, which can obtain the optimal performance with low computational complexity and is suitable for large-scale OAM systems. By analyzing the simulation results, we found that the mode reuse technique does not offer additional SE for single user OAM wireless communications. Thus, we made a conclusion that this orthogonal transfer technology only really delivers performance gains when multiple data transmissions simultaneously and multi-user environments are considered. Besides, in our future work, more practical conditions or even experience tests need to be considered for the theoretical verification.

Author Contributions

Conceptualization, T.H. and Y.W.; methodology, T.H.; software, T.H.; validation, T.H., Y.W. and J.Z.; formal analysis, T.H.; investigation, T.H.; resources, T.H.; data curation, T.H.; writing–original draft preparation, T.H.; writing–review and editing, T.H., and Y.W.; visualization, T.H. and B.M.; supervision, Y.W. and J.Z.; project administration, Y.W. and J.Z.; funding acquisition, Y.W. and J.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the International Science and Technology Cooperation Program of China under Grant 2017YFE0118900, the National Natural Science Foundation of China under Grant 61801062,61831002, the Chongqing Natural Science Foundation under Grant cstc2019jcyj-msxmx0288, and in part by the High-end Project for Doctoral Students at Chongqing University of Posts and Telecommunications under Grant BYJS201801.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Thidé, B.; Tamburini, F.; Then, H.; Someda, C.G.; Ravanelli, R.A. The physics of angular momentum radio. In Proceedings of the 2015 1st URSI Atlantic Radio Science Conference (URSI AT-RASC), Las Palmas, Gran Canaria, Spain, 16–24 May 2015. [Google Scholar]
  2. Tamburini, F.; Mari, E.; Sponselli, A.; Thidé, B.; Bianchini, A.; Romanato, F. Encoding many channels in the same frequency through radio vorticity: First experimental test. New J. Phys. 2011, 14, 78001–78004. [Google Scholar]
  3. Zhu, Q.; Jiang, T.; Qu, D.; Chen, D.; Zhou, N. Radio Vortex–Multiple-Input Multiple-Output Communication Systems with High Capacity. IEEE Access 2015, 3, 2456–2464. [Google Scholar] [CrossRef]
  4. Gaffoglio, R.; Cagliero, A.; Vita, A.D.; Sacco, B. OAM multiple transmission using uniform circular arrays: Numerical modeling and ex- perimental verification with two digital television signals. Radio Sci. 2016, 51, 645–658. [Google Scholar] [CrossRef] [Green Version]
  5. Liu, S.; Zhang, Y.; Fu, C.L.; Bai, Z.; Li, Z.; Liao, C.; Wang, Y.; He, J.; Liu, Y.; Wang, Y. Temperature Insensitivity Polarization-Controlled Orbital Angular Momentum Mode Converter Based on an LPFG Induced in Four-Mode Fiber. Sensors 2018, 18, 1766. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Wang, J.; Yang, J.Y.; Fazal, I.M.; Ahmed, N.; Yan, Y.; Huang, H.; Ren, Y.; Yue, Y.; Dolinar, S.; Tur, M.; et al. Terabit free-space data transmission employing orbital angular momentum multiplexing. Nat. Photonics 2012, 6, 488–496. [Google Scholar] [CrossRef]
  7. Zhang, W.T.; Zheng, S.L.; Hui, X.N.; Dong, R.; Jin, X.; Chi, H.; Zhang, X. Mode Division Multiplexing Communication Using Microwave Orbital Angular Momentum: An Experimental Study. IEEE Trans. Wirel. Commun. 2017, 16, 1308–1318. [Google Scholar] [CrossRef]
  8. Yan, Y.; Xie, G.; Lavery, M.; Huang, H.; Ahmed, N.; Bao, C.; Ren, Y.; Cao, Y.; Li, L.; Zhao, Z.; et al. High-capacity millimetre-wave communications with orbital angular momentum multiplexing. Nat. Commun. 2014, 5, 1–9. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Hu, T.; Wang, Y.; Liao, X.; Zhang, J.; Song, Q.L. OFDM-OAM modulation for future wireless communications. IEEE Access 2019, 7, 59114–59125. [Google Scholar] [CrossRef]
  10. Bu, X.; Zhang, Z.; Liang, X.D.; Chen, L.; Tang, H.; Zeng, Z.; Wang, J. A Novel Scheme for MIMO-SAR Systems Using Rotational Orbital Angular Momentum. Sensors 2018, 18, 3511. [Google Scholar] [CrossRef] [Green Version]
  11. Liu, K.; Cheng, Y.; Li, X.D.; Qin, Y.; Wang, H.; Jiang, Y. Generation of Orbital Angular Momentum Beams for Electromagnetic Vortex Imaging. IEEE Antennas Wirel. Propag. Lett. 2016, 15, 1873–1876. [Google Scholar] [CrossRef]
  12. Mozaffari, M.; Saad, W.; Bennis, M.; Nam, Y.H.; Debbah, M. A Tutorial on UAVs for Wireless Networks: Applications, Challenges, and Open Problems. IEEE Commun. Surv. Tutor. 2018, 21, 2334–2360. [Google Scholar] [CrossRef] [Green Version]
  13. Jiang, H.; Zhang, Z.C.; Wu, L.; Dang, J. Three-dimensional geometry-based UAV-MIMO channel modeling for A2G communication environments. IEEE Commun. Lett. 2018, 7, 1438–1441. [Google Scholar] [CrossRef]
  14. Liu, Y.W.; Qin, Z.J.; Cai, Y.L. UAV communications based on non-orthogonal multiple access. IEEE Wirel. Commun. 2019, 26, 52–57. [Google Scholar] [CrossRef] [Green Version]
  15. Lee, T.R.; Buban, M.; Dumas, E.; Baker, C.B. On the Use of Rotary-Wing Aircraft to Sample Near-Surface Thermodynamic Fields: Results from Recent Field Campaigns. Sensors 2019, 19, 10. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Alaoui-Sosse, S.; Durand, P.; Medina, P.; Pastor, P.; Lothon, M.; Cernov, I. OVLI-TA: An Unmanned Aerial System for Measuring Profiles and Turbulence in the Atmospheric Boundary Layer. Sensors 2019, 19, 581. [Google Scholar] [CrossRef] [Green Version]
  17. Rautenberg, A.; Schön, M.; zum Berge, K.; Mauz, M.; Manz, P.; Platis, A.; van Kesteren, B.; Suomi, I.; Kral, S.T.; Bange, J. The Multi-Purpose Airborne Sensor Carrier MASC-3 for Wind and Turbulence Measurements in the Atmospheric Boundary Layer. Sensors 2019, 19, 2292. [Google Scholar] [CrossRef] [Green Version]
  18. Rautenberg, A.; Graf, M.S.; Wildmann, N.; Platis, A.; Bange, J. Reviewing Wind Measurement Approaches for Fixed-Wing Unmanned Aircraft. Atmosphere 2018, 9, 422. [Google Scholar] [CrossRef] [Green Version]
  19. Witte, B.M.; Singler, R.F.; Bailey, S.C.C. Development of an Unmanned Aerial Vehicle for the Measurement of Turbulence in the Atmospheric Boundary Layer. Atmosphere 2017, 8, 195. [Google Scholar] [CrossRef] [Green Version]
  20. Nolan, P.J.; McClelland, H.G.; Woolsey, C.A.; Ross, S.D. A Method for Detecting Atmospheric Lagrangian Coherent Structures Using a Single Fixed-Wing Unmanned Aircraft System. Sensors 2019, 19, 1607. [Google Scholar] [CrossRef] [Green Version]
  21. Nolan, P.J.; Pinto, J.; González-Rocha, J.; Jensen, A.; Vezzi, C.N.; Bailey, S.C.C.; De Boer, G.; Diehl, C.; Laurence, R., III; Powers, C.W.; et al. Coordinated Unmanned Aircraft System (UAS) and Ground-Based Weather Measurements to Predict Lagrangian Coherent Structures (LCSs). Sensors 2018, 18, 4448. [Google Scholar] [CrossRef] [Green Version]
  22. Schuyler, T.J.; Gohari, S.M.I.; Pundsack, G.; Berchoff, D.; Guzman, M.I. Using a Balloon-Launched Unmanned Glider to Validate Real-Time WRF Modeling. Sensors 2019, 19, 1914. [Google Scholar] [CrossRef] [PubMed]
  23. Schuyler, T.J.; Bailey, S.C.C.; Guzman, M.I. Monitoring Tropospheric Gases with Small Unmanned Aerial Systems (sUAS) during the Second CLOUDMAP Flight Campaign. Atmosphere 2019, 10, 434. [Google Scholar] [CrossRef] [Green Version]
  24. Madridano, Á.; Al-Kaff, A.; Martín, D.; de la Escalera, A.A. 3D Trajectory Planning Method for UAVs Swarm in Building Emergencies. IEEE Access 2020, 20, 642. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Da Rosa, R.; Aurelio Wehrmeister, M.; Brito, T.; Lima, J.L.; Pereira, A.I.P.N. Honeycomb Map: A Bioinspired Topological Map for Indoor Search and Rescue Unmanned Aerial Vehicles. Sensors 2020, 20, 907. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Challita, U.; Saad, W. Network formation in the Sky: Unmanned aerial vehicles for multi-hop wireless backhauling. In Proceedings of the IEEE Global Telecommunications Conference (GLOBECOM), Singapore, 4–8 December 2017. [Google Scholar]
  27. Ge, X.H.; Zi, R.; Xiong, X.S.; Li, Q.; Wang, L. Millimeter wave communications with OAM-SM scheme for future mobile networks. IEEE J. Sel. Areas Commun. 2017, 35, 2163–2177. [Google Scholar] [CrossRef] [Green Version]
  28. Gharavi-Alkhansari, M.; Gershman, A.B. Fast antenna subset selection in MIMO systems. IEEE Trans. Signal Process. 2004, 52, 339–347. [Google Scholar] [CrossRef]
  29. Molisch, A.F.; Win, M.Z.; Winters, J.H. Capacity of MIMO systems with antenna selection. IEEE Trans. Wirel. Commun. 2005, 4, 1759–1772. [Google Scholar] [CrossRef] [Green Version]
  30. Pal, R.; Srinivas, K.V.; Chaitanya, A.K. A Beam Selection Algorithm for Millimeter-Wave Multi-user MIMO Systems. IEEE Commun. Let. 2018, 22, 852–855. [Google Scholar] [CrossRef]
  31. Narendra, P.M.; Fukunaga, K. A Branch and Bound Algorithm for Feature Subset Selection. IEEE Trans. Comput. 1977, 26, 917–922. [Google Scholar] [CrossRef]
  32. Gao, Y.; Khaliel, M.; Zheng, F.; Kaiser, T. Rotman Lens Based Hybrid Analog-Digital Beamforming in Massive MIMO Systems: Array Architectures, Beam Selection Algorithms and Experiments. IEEE Trans. Veh. Technol. 2017, 66, 9134–9148. [Google Scholar] [CrossRef]
Figure 1. The scenario of the UAV-based radio vortex wireless network.
Figure 1. The scenario of the UAV-based radio vortex wireless network.
Sensors 20 02289 g001
Figure 2. The LOS geometric channel model.
Figure 2. The LOS geometric channel model.
Sensors 20 02289 g002
Figure 3. Illustration of the BBS-MS searching tree with N t = 8 and L = 3 .
Figure 3. Illustration of the BBS-MS searching tree with N t = 8 and L = 3 .
Sensors 20 02289 g003
Figure 4. Amptitude gain | h l | , SNR = 20 dB and R = 20 λ : (a) h t = 10 m ; (b) h t = 20 m ; (c) h t = 50 m ; (d) h t = 100 m .
Figure 4. Amptitude gain | h l | , SNR = 20 dB and R = 20 λ : (a) h t = 10 m ; (b) h t = 20 m ; (c) h t = 50 m ; (d) h t = 100 m .
Sensors 20 02289 g004
Figure 5. Spectrum efficiency, R = 20 λ , SNR = 20 dB, and L = 2 .
Figure 5. Spectrum efficiency, R = 20 λ , SNR = 20 dB, and L = 2 .
Sensors 20 02289 g005
Figure 6. Spectrum efficiency, R = 20 λ and SNR = 20 dB.
Figure 6. Spectrum efficiency, R = 20 λ and SNR = 20 dB.
Sensors 20 02289 g006
Table 1. Simulation Parameters.
Table 1. Simulation Parameters.
ParametersValues
The number of transmitting antennas N t = 8
The number of receiving antennas M = 2
The number of the available modes L = 2
The relative rotation angle α = 0
The radius of the transmitter array R = 20 λ
The transmission SNRSNR = 20 dB
The altitude of the ground user h r = 1.5 m
The altitude of the drone h t = 20 m
The centre frequency of carrier wave f = 70 GHz

Share and Cite

MDPI and ACS Style

Hu, T.; Wang, Y.; Ma, B.; Zhang, J. Orbit Angular Momentum MIMO with Mode Selection for UAV-Assisted A2G Networks. Sensors 2020, 20, 2289. https://doi.org/10.3390/s20082289

AMA Style

Hu T, Wang Y, Ma B, Zhang J. Orbit Angular Momentum MIMO with Mode Selection for UAV-Assisted A2G Networks. Sensors. 2020; 20(8):2289. https://doi.org/10.3390/s20082289

Chicago/Turabian Style

Hu, Tao, Yang Wang, Bo Ma, and Jie Zhang. 2020. "Orbit Angular Momentum MIMO with Mode Selection for UAV-Assisted A2G Networks" Sensors 20, no. 8: 2289. https://doi.org/10.3390/s20082289

APA Style

Hu, T., Wang, Y., Ma, B., & Zhang, J. (2020). Orbit Angular Momentum MIMO with Mode Selection for UAV-Assisted A2G Networks. Sensors, 20(8), 2289. https://doi.org/10.3390/s20082289

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