[go: up one dir, main page]
More Web Proxy on the site http://driver.im/
Next Article in Journal
A Consistent Theory of Tachyons with Interesting Physics for Neutrinos
Next Article in Special Issue
Cross-Layer Optimization-Based Asymmetric Medical Video Transmission in IoT Systems
Previous Article in Journal
Research on Real-Time Communication Algorithm of Substation Based on Time-Sensitive Network
Previous Article in Special Issue
Joint Offloading Decision and Resource Allocation in Mobile Edge Computing-Enabled Satellite-Terrestrial Network
You seem to have javascript disabled. Please note that many of the page functionalities won't work as expected without javascript enabled.
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Location Information-Assisted Robust Beamforming Design for Ultra-Wideband Communication Systems

1
School of Information and Control Engineering, China University of Mining and Technology, Xuzhou 221116, China
2
Autonomous Systems and Robotics Lab, ENSTA Paris, Institut Polytechnique de Paris, 91120 Palaiseau, France
*
Author to whom correspondence should be addressed.
Symmetry 2022, 14(6), 1171; https://doi.org/10.3390/sym14061171
Submission received: 9 May 2022 / Revised: 31 May 2022 / Accepted: 3 June 2022 / Published: 7 June 2022

Abstract

:
The future of mobile communication systems is evolving rapidly toward being more intelligent, while also having the ability to interconnect with everything and be aware of the current wireless environment. For complex scenarios, the asymmetry features of channel state information (CSI) will seriously restrict the performance of the communication system. With its accurate positioning technology and high-speed communication rate, ultra-wideband (UWB) has a promising future as a solution to integrate communication and positioning functions. The traditional CSI channel estimation usually requires channel training, which will greatly increase the overhead of the system. This paper proposes a location information-assisted beamforming to replace the traditional channel training process. First, we use the user location parameter information to reconstruct the channel model and derive the CSI error distribution based on the location distribution. Second, considering the uncertainty of the user positioning error, we model the robust beamforming optimization problem that minimizes the total transmit power. To solve this non-convex problem effectively, we design a new beamforming algorithm by using semidefinite relaxation (SDR) and Bernstein-type inequalities. Finally, simulation results verify the robustness of the proposed robust beamforming compared to the worst-case robust beamforming.

1. Introduction

With the development of mobile connectivity and the Internet of Things (IoT), short-range communication has attracted a lot of attention. For financial and practical reasons, it is more attractive to integrate positioning functions into existing communication technologies than to have two separate communication and positioning systems [1]. To provide better communication services for users with different channel state information, the base station (BS) needs to allocate resources asymmetrically according to different CSI [2]. However, due to information asymmetry, it is difficult to obtain the pecfect status information at the base station. Therefore, this paper will obtain the user’s location information through the positioning algorithm at the base station to assist beamforming. However, positioning performance in indoor environments where a global positioning system (GPS) is unavailable is necessarily affected by multipath propagation and wall obstruction. The current commonly used positioning methods, such as Bluetooth low Energy (BLE), Wi-Fi, and Zig-Bee. It is difficult to meet the communication and location needs of the next generation Internet of things mobile system. Among them, BLE and Zig-Bee technology can greatly reduce the complexity and energy consumption of system design, but their main drawback is the limited data rate [3]. The latest Wi-Fi 6 is capable of reaching a maximum data rate of 9.6 Gb/s, but it suffers from high system complexity and high power consumption, and it may not be used in scenarios with high-precision positioning requirements [4].
Ultra-wideband (UWB) technology, with its significant characteristics [5] such as high transmission rate, low power consumption, and cost-efficiency, has become an indispensable technology in the field of the Internet of Things. Based on the Shannon channel capacity theory, conventional UWB systems can achieve extremely high data transmission rates under a huge system bandwidth by transmitting very narrow Gaussian pulses [6]. For UWB positioning performance, the wide bandwidth of the UWB signal can support centimeter-level ranging accuracy and can resolve multipath delay at nanosecond resolution, so it is very suitable for short-range accurate positioning, communication, and sensing [7,8]. In addition, UWB technology possesses high robustness to multipath fading and good coexistence with other narrowband systems. UWB communication technology is a promising candidate for improving the data rate. Therefore, UWB is suitable for short-range indoor wireless communication and can be combined with cognitive radio technology, beamforming, and multiple-input multiple-output (MIMO) technology to achieve better system performance [2,9,10].
The main problem of UWB is the limited communication transmission range due to low power transmission, the multipath effect, and interference. The combination of MIMO and UWB technologies offers a feasible solution to the power limitation problem of UWB. The integration of UWB and MIMO enables wireless communication systems to increase data rates by transmitting signals over multiple channels without increasing the transmit power [11]. In indoor, dense multipath propagation environments, MIMO technology makes full use of array gain and spatial multiplexing gain to significantly improve spectral efficiency [12]. Beamforming, an important technology for MIMO, can increase the system’s anti-interference capability in order to suppress interference, improve data rates, etc.
UWB beamforming is different from narrowband signal beamforming. Compared with narrowband beamforming, UWB beamforming needs fewer sensors to obtain the same angle beam, and the uniform power allocation can provide a sidelobe-free beam mode [13,14]. Through the design of a sparse array, the intercoupling between UWB antennas is effectively reduced. In [14], Wang, M. et al. proposed a time-domain beamforming algorithm that adaptively guides the array beam by estimating the direction of arrival of the signal and adjusting the corresponding time delay. Considering the nanosecond duration of the UWB pulse signal in the time domain, a least-squares-based beamforming achieved by taking full advantage of frequency invariance was proposed in [15]. Based on the geometry of the UWB antenna array, a new pointing beam was proposed to eliminate other interferences in the same frequency [16]. However, for multi-user communication systems, the performance of the proposed directional beamforming (OBF) technology was seriously degraded. A nonlinear OBF technique based on an improved radial basis function neural network was proposed in [17], which could effectively eliminate the interference caused by directional UWB. However, the array elements of the two antennas needed to be matched one by one; consequently, the implementation complexity of the system was high.
It is still worth pointing out that most of the systems designed in the literature implement communication and location functions independently, which leads to a waste of spectrum resources. At present, some researchers have paid more and more attention to location information-assisted communication, namely the use of location estimation information to improve the system capacity, thus improving the Quality of Service (QoS) of communication services [18]. The joint optimization problem of maximizing the total rate under the localization constraint and minimizing the localization error under the communication constraint was presented in a large-scale MIMO system, and the results showed that high data rate communication and mobile 3D localization can be achieved simultaneously [19]. A joint location and beamforming scheme was proposed to jointly estimate the user’s position and instantaneous CSI, and then to optimize the beamforming vector according to the estimated results [20]. The existing research mainly use the estimated user location information to assist the communication process and then improve the communication quality. The common method is to take the location estimation as a constraint in order to maximize the communication rate. However, the coupling relationship between communication and location has not been considered in the existing research, and the impact of the location estimation error on CSI estimation and the communication rate has not been taken into account.
In this paper, we will make use of the inherent coupling relationship between communication rate and positioning accuracy, robust beamforming is designed by introducing CSI estimation error in the presence of multipath interference. The mechanism mainly consists of two key components. For the location, a weighted least squares (WLS) localization algorithm using reference point coordinate information is proposed. Combined with Taylor expansion and multipath information, an approximate CSI estimation error distribution is derived from the user position error distribution constructed by the Cramer–Rao Lower Bound (CRLB). On this basis, a robust beamforming optimization problem to minimize the total transmit power is formulated. Since the non-convex optimization problem is difficult to solve, a suboptimal method based on the Bernstein-type inequality and SDR is proposed. More specifically, robust beamforming is accomplished by converting the SINR probability constraint into a deterministic constraint and by transforming the unsolvable nonconvex problem into a solvable convex problem.

2. System Model

The paper considers the downlink transmission in a cellular network consisting of J cells, each having an N multi-antenna BS to transmit independent messages to one active single-antenna MS, as shown in Figure 1. It is assumed that all BSs share the same narrowband spectrum for downlink transmission. This is assuming that the exact locations of the base station (BS) and the reference node (RN) are fixed and known after their deployments. Let us denote J = Δ 1 , 2 , , J as the set of the BS, the location of the jth BS is u j = [ u j , x , u j , y ] T , and the location of the RN is p f = [ p f , x , p f , y ] T . Each user equipment (UE) is equipped with a transmitting module for broadcasting and transmitting UWB pulse signals, and the location of the kth UE is expressed as p k = [ p k , x , p k , y ] T .
Each symbol is conveyed by repeating one pulse per frame (frame duration T f T w ) over N f frames, resulting in a low duty cycle transmission form that is beneficial for distinguishing the multipath signal in terms of time resolution; T w denotes each pulse duration. The transmitted signal of this system adopts the DS-UWB modulation mode, and the transmitted symbol from the j t h BS to the kth UE is expressed as follows [21]:
s j , k t = m = 0 d j , m ( k ) i = 0 N f 1 A m ( k ) p t m T s i T f .
where T s = N f T f indicates the symbol duration, d j , m ( k ) is the mth message sent to kth UE, and A m ( k ) ± 1 is the spread spectrum sequence corresponding to the kth UE. The second-order derivative of the Gaussian function has been widely used in practice and is used as the model of UWB signal p ( t ) . Subsequently, we will explain in detail the geometrical model, the TDOA correction method, and the CSI estimation error model.

2.1. Geometrical Model

The complex channel vector between the n t h antenna of the jth BS and the kth user is expressed as follows:
h j , k n ( t ) = m = 1 M 1 β j k , m n δ ( t τ j k , m n ) ,
where τ j k , m denotes the TOA of the mth path between the kth UE and the jth BS. The complex channel fading amplitude β j k , m n of the mth path is modeled as [22]:
β j k , m n = λ 4 π d j k , 1 n e j 2 π d j k , 1 n λ , LoS , λ Γ R 4 π d j k , m n e j 2 π d j k , m n λ , reflector ,
where Γ R denotes the reflection coefficient, d j k , 1 n = u j p k , n denotes the direct shooting distance, d j k , m = a j k , m p k , n denotes the distance from the jth BS and the virtual point. According to the application of the image-source model for a first-order MPC, a j k , m denotes the mirror position of the jth BS on the m t h path, which is described in detail below.
We consider that multipath signals are reflected from the plane when the prior information of the position is known, such as walls, doors, and tables, which are represented as wall segments below [23,24]. Each wall segment s S = { 1 , . . . , S } is described by its location p s R 2 (an endpoint of the wall segment) and orientation l s e s , with l s as length and the unit-vector e s R 2 as the direction of the wall segment. To relate the estimated path delays to the environment, we employ a geometric model, as shown in Figure 2. The application of the image-source model for a first-order MPC is equivalent to mirroring the BS position u j at segment s S , expressed as:
a j , s = u j 2 T s u j p s .
where T s = U π 2 e s e s U π 2 is the rotation matrix, and U π 2 = 0 1 1 0 .

2.2. TDOA Correction Method

In this paper, we will transmit communication and positioning signals simultaneously in the TDMA mode. Considering that the distance between the antennas is much smaller than the distance between the BS and the UE, it is reasonable to assume that the delay in arriving at each antenna of the base station is the same τ j k , m n = τ j k , m . For the kth UE, the received signal vector y k ( p ) ( t ) = [ y 1 , k ( p ) ( t ) , , y J , k ( p ) ( t ) ] is expressed as follows:
y j , k ( p ) ( t ) = m = 1 M h j k , m w j s j ( p ) ( t τ j k , m ) + n k ( t ) , t 0 , T p ,
where the positioning signal of the jth BS is s j ( p ) ( t ) , t [ 0 , T p ] , w j C N represents the beamforming vector of the jth BS, and w j = [ w j , 1 , , w j , N ] T . h j k , m C 1 × N = β j k , m 1 , , β j k , m N represents the channel model between the jth BS and the kth UE, and τ j k , m represents the TOA of the mth path. In addition, T p = N p T s indicates the duration of the positioning signal, and N p denotes the number of positioning frames.
For the receiver of the UWB, the first path component (FPC) plays a very important role in finding out the receiving time of the message [25]. The received signal power is likely to be more concentrated in the FPC than in the other multipath components (MPC). In reality, UE localization is subject to equipment installation error of BS, continuous start-up errors of the UWB, and clock synchronization errors related to the external environment such as the propagation distance of pulse signals and temperature [26]. To obtain accurate real-time localization of the UE, we design a distance-corrected localization algorithm that uses RN to correct the frequency and time differences between base stations. In addition, this paper will adopt the optional bilateral two-way (AltDS–TWR) ranging mode [27]. The jth BS measures the distance from the kth UE as [28]:
d ˜ j , k ( t ) = 1 + e j d j , k ( t ) + B j + ω j , k ,
where d j , k ( t ) denotes the real distance between the j t h BS and k t h UE, e j is the proportional error of the j t h BS, B j is the measurement deviation of the j t h BS, and ω j , k is a zero-mean white Gaussian noise with a power spectral density N 0 . The TOA value measured by the j t h BS from the RN is as follows:
d ˜ j , f ( t ) = 1 + e j d j , f ( t ) + B j + ω j , f ,
where d j , f ( t ) denotes the distance between the j t h BS and RN. The unknown error vectors δ j = [ e j B j ] T of the j t h BS are estimated by the classical LS algorithm using the known coordinates of the BS and RN combined with the measured TOA data, which as follow:
δ ^ j = F j T F j 1 F j T f j ,
with
F j = d j , f ( t ) 1 d j , f ( t Δ T ) 1 d j , f ( t 2 Δ T ) 1 , f j = d ˜ j , f ( t ) d j , f ( t ) d ˜ j , f ( t Δ T ) d j , f ( t Δ T ) d ˜ j , f ( t 2 Δ T ) d j , f ( t 2 Δ T ) ,
where Δ T denotes the sample period. After obtaining the error vector of the j t h BS, this paper will correct the measured distance data between the k t h UE and the j t h BS in real time as follow:
d ˜ j , k t B ^ j c 1 + e ^ j = u j p k t + ω j , k ,
Square the two sides of (10) and subtract the equation of a common RN from the corresponding equations of the other BS, then a set of linear equations of the UE location as follow:
A k p k ( t ) = B k ( t ) ,
where
A k = 2 u 2 T u 1 T u J T u 1 T , B k ( t ) = d ˜ 2 , k ( t ) B ^ 2 1 + e ^ 2 2 d ˜ 1 , k ( t ) B ^ 1 1 + e ^ 1 2 u 2 T u 2 + u 1 T u 1 d ˜ J , k ( t ) B ^ J 1 + e ^ J 2 d ˜ 1 , k ( t ) B ^ 1 1 + e ^ 1 2 u J T u J + u 1 T u 1 ,
The classical LS localization method does not consider the effect of noise characteristics. To improve the localization accuracy, the WLS method using RN distance information will be designed, which is given by:
p ^ k ( t ) = A k T Q k A k 1 A k T Q k B k ( t ) .
where Q k denotes the weighting coefficient matrix.

3. CSI Estimation Error Model

The goal of this paper is to derive an error model for CSI estimation based on the UE location and known environmental parameters. We intend to achieve this in two steps: firstly, the position error bound of UE is derived from the observed signal. Secondly, the channel error model is reconstructed under the assumption that the localization error distribution follows a Gaussian distribution [29]. After the location information of UE is obtained by the location algorithm, in order to analyze the location accuracy of UE, we will derive the CRLB of UEs to construct the location error distribution. Finally, we define the following vectors:
τ k T = [ τ 1 , k T , τ 2 , k T , . . . , τ J , k T ] , β k T = [ β 1 , k T , β 2 , k T , . . . , β J , k T ] ,
where τ j , k = [ τ j k , 1 , τ j k , 2 , . . . , τ j k , M ] T and β j , k = [ β j k , 1 , β j k , 2 , . . . , β j k , M ] T . Then, define the parameter vector of the k t h UE as follows:
φ k = [ τ k T , R β k T , I β k T ] T ,
where R β k = Δ β k and I β k = Δ β k represent the real and imaginary parts of the amplitude between the UE and BS, respectively. The mean square error matrix of unbiased estimation φ ^ k satisfies the information inequality as follows [30]:
E ( φ ^ k φ k ) ( φ ^ k φ k ) T J k 1 ,
The element of the FIM matrix J k C 3 J M × 3 J M can be expressed as:
[ J k ] u , v = Δ 1 N 0 0 T p { μ k H ( t ) φ k , u μ k ( t ) φ k , v } d t ,
where μ k ( t ) is the noiseless part of the received signal in (5), given by the following equation:
μ k ( t ) = Δ n = 1 N m = 1 M β 1 k , m w 1 , n s 1 ( p ) ( t τ 1 k , m ) n = 1 N m = 1 M β J k , m w J , n s J ( p ) ( t τ J k , m ) ,
Therefore, the corresponding FIM matrix J k C 3 J M × 3 J M , partitioned into J M × J M submatrices, is structured as:
J k = J τ k , τ k J τ k , R β k J τ k , I β k J τ k , R β k T J R β k , R β k J R β k , I β k J τ k , I β k T J I β k , R β k J I β k , I β k ,
Each submatrix in (19) is expressed as follows:
J τ k , τ k = diag { J τ 1 , k , τ 1 , k , . . . , J τ J , k , τ J , k } , J τ k , R β k = diag { J τ 1 , k , R β 1 , k , . . . , J τ J , k , R β J , k } , J τ k , I β k = diag { J τ 1 , k , I β 1 , k , . . . , J τ J , k , I β J , k } , J R β k , R β k = diag { J R β 1 , k , R β 1 , k , . . . , J R β J , k , R β J , k } , J I β k , I β k = diag { J I β 1 , k , I β 1 , k , . . . , J I β J , k , I β J , k } ,
where J τ j , k , τ j , k C M × M , J τ j , k , R β j , k C M × M , J τ j , k , I β j , k C M × M , J I β 1 , k , I β 1 , k C M × M , and J I β 1 , k , I β 1 , k C M × M are denoted as:
J R β j , k , R β j , k = J I β j , k , I β j , k = N p N 0 n = 1 N w j , n 2 { R 0 } , J R β j , k , I β j , k = j N p N 0 n = 1 N w j , n 2 { R 0 } , J τ j , k , τ j , k = N p N 0 n = 1 N w j , n 2 { B k H R 2 B k } , J τ j , k , R β j , k = N p N 0 n = 1 N w j , n 2 { R 1 B k } ,
where B k = diag β k ; for details, see Appendix A. Assuming that the symbols are independent and identically distributed, the matrix R i elements are expressed as follows:
[ R i ] u , v = Δ W W 2 2 W W 2 2 ( 2 π f ) i P ( f ) 2 e j 2 π f Δ τ u , v d f ,
in which Δ τ u , v = τ j k , u τ j k , v , i = { 0 , 1 , 2 } , and W is the signal bandwidth. As can be seen from the formula above, when i is 0 or 2, the diagonal element of R 0 and R 2 is the largest 1 and 4 π 2 W eff 2 , respectively, and W eff 2 = W W 2 2 W W 2 2 f 2 P ( f ) 2 d f is the effective bandwidth. In addition, for any bandwidth W, diag ( R 1 ) = 0 M , lim W R 1 = 0 M × M . In order to directly calculate the CRLB of the UE position, it is necessary to redefine the unknown parameter vector according to the chain rule, and the equivalent Fisher matrix of the position can be expressed as follows [28]:
J p k = b k H J τ k τ k b k ,
where the elements of the vector b k R J M × 1 are expressed as:
[ b k ] n = p k a j n , k , m n c p k a j n , k , m n 2 .
where a j n , k , m n denotes the reflection point coordinates of the m t h from the j t h BS for the kth UE, where m n and j n are related to n, respectively, by j n = n n M M , m n = n mod M .
According to the reference point auxiliary location algorithm proposed in the previous section, the location p ^ k ( t ) of the UE can be estimated in real time, but in practice, due to the BS coordinate error and measurement error, it is impossible for the estimation algorithm to obtain accurate UE coordinates. The estimation error at this time is expressed as e k ( t ) = p k ( t ) p ^ k ( t ) . Considering the assumption of Gaussian white noise, the least square equivalent and maximum likelihood estimation method can achieve the estimation result under CRLB; that is, the positioning error e k obeys the mean value 0, while the variance is J p k 1 .
e p k ( t ) CN ( 0 , J p k 1 ) ,
with the location of the kth UE and the known environmental parameters, the partial channel model β ^ j k , m ( t ) between the kth UE and the j t h BS can be constructed, and its elements are expressed as β ^ j k , m ( t ) = β ^ j k , m 1 ( t ) , , β ^ j k , m N ( t ) . The channel estimation error is expressed as follows:
Δ β j k , m n ( t ) λ 4 π p ^ k ( t ) u j , n T e k ( t ) p ^ k ( t ) u j 2 3 , m = 1 λ Γ R p ^ k ( t ) u j , n T e k ( t ) 4 π ( p ^ k ( t ) a k , m 2 + a k , m u j 2 ) 3 , m 2
Furthermore, this paper uses the first-order Taylor expansion formula, given by the following equation:
p ^ k ( t ) + e k ( t ) u j 2 p ^ k ( t ) u j 2 p ^ k ( t ) u j T p ^ k ( t ) u j 2 e k ( t ) .

4. Robust Beamforming

This section is based on the estimated position of kth UE for beamforming. The distance between the UE and all BSs is first calculated, and then the closest base station is selected for communication services. For the convenience of representation, we consider the jth UE and the jth base station to be the closest. The received signal of the j t h UE is represented by the following equation:
y j ( c ) ( t ) = m = 1 M β ^ j j , m + Δ β j j , m w j s j ( c ) ( t τ j j , m ) + i j J m = 1 M ( β ^ i j , m + Δ β i j , m ) w i s i ( c ) ( t τ i j , m ) + n j ( t ) , t [ T p , T p + T c ] ,
where s j ( c ) ( t ) denotes the communication signal of the jth BS, and T c is the communication duration. The rate of the jth UE served by the jth BS is expressed as follows:
R j = T c T p + T c log 1 + ( g ^ j j + Δ g j j ) w j 2 i j J ( g ^ i j + Δ g i j ) w i 2 + σ j 2 ,
where g ^ j j = m = 1 M β ^ j j , m , Δ g j j = m = 1 M Δ β j j , m . The nth element of Δ g j j is represented by Δ g j j , n = b j j , n T e j , with:
b j j , n = λ 4 π p ^ j u j , n p ^ j u j , n 2 3 + m = 2 M λ Γ R p ^ j u j , n 4 π ( p ^ j a j 2 + a j u j , n 2 ) 3 ,
Each BS is scattered throughout the coverage area according to certain rules. When the BS receives the positioning data signal from the UE, it is connected to the central processing unit through high-speed transmission media such as optical fiber. The robust beamforming optimization problem to minimize the total transmit power is formulated as follows:
min { w j } j = 1 J j = 1 J w j 2 2 s . t . w j 2 2 p UWB max , j = 1 , . . . , J , Pr { SINR j ( Δ g j j H ) γ j } ρ j ,
where p UWB max represents the maximum power, γ j = 2 T p + T c R ¯ T p + T c R ¯ T c T c 1 denotes the minimum SINR constraint of the jth UE and the maximum outage probability ρ j of the jth UE. We define W j = Δ w j w j H . To simplify the calculation, the SINR of the jth BS is expressed as follows:
f ( Δ g j j ) = g ^ j j H W j g ^ j j + 2 Re Δ g j j H W j g j j + Δ g j j H W j Δ g j j ,
Therefore, the constraint (31) of the optimization problem is represented by the following equation:
Pr f ( Δ g j j ) γ j i = 1 , l j J f ( Δ g i j ) σ j 2 p out ,
According to formula (25), the signal amplitude error is expressed as:
Δ g j j C N ( 0 , Δ B j j T Δ B j j J p k 1 ) ,
where Δ B j j R N × 2 , Δ B j j = b j j , 1 T , , b j j , N T T . The signal amplitude error is rewritten as:
Δ g j j = E j j v j ,
where v j N ( 0 , I ) , E j j = Δ B j j H Δ B j j J p j 1 . In order to simplify the equation, we defined E ¯ j = diag { E 1 , j , . . . , E J , j } and g ¯ j = g 1 , j T , . . . , g J , j T T . Then, we can obtain the following equation:
Δ g ¯ j = E ¯ j g ¯ j ,
where
Δ g ¯ j = [ Δ g 1 , j T , . . . , Δ g J , j T ] T , W ¯ j = diag { W 1 , . . . , W j 1 , W j γ j , W j + 1 , . . . , W J } ,
then
f ( Δ β j j ) γ j i = 1 , i j J f ( Δ β i j ) = v ¯ j T A ¯ j v ¯ j * + 2 Re v ¯ j T a ¯ j + g ¯ j T W ¯ j g ¯ j * ,
where
A ¯ j = E ¯ j W ¯ j E ¯ j , a ¯ j = E ¯ j W ¯ j g ¯ j ,
where ⊙ denotes the Hadamard product. Therefore, the constraint item (31) can be expressed as:
Pr { v ¯ j T A ¯ j v ¯ j * + 2 Re { v ¯ j T a ¯ j } c j } ρ j ,
where c j = σ j 2 g ¯ j T W ¯ j g ¯ j * . In order to transform probabilistic constraints into deterministic constraints, this paper uses the following Bernstein-type inequalities: G = v H Av + 2 v H a , A H N , a C N , v C N ( 0 , I ) . For any ς > 0 , the following inequality holds:
Pr { G Tr ( A ) 2 ς vec ( A ) 2 + 2 a 2 ς s ( A ) } e x p ( ς ) ,
where s ( A ) = max { λ m a x ( A ) , 0 } , s + ( A ) = max { λ m a x ( A ) , 0 } , λ max ( A ) is the largest eigenvalue of matrix A. According to the Bernstein-type inequality, the constraint term (33) is the following deterministic constraint:
Tr ( A ¯ j ) 2 ς j vec ( A ¯ j ) 2 + 2 a ¯ j 2 ς j s ( A ¯ j ) c j ,
where ς j = ln ( ρ j ) . Furthermore, the constraints of (42) can be reformulated as:
(43a) Tr ( A ¯ j ) 2 ς j μ j ς j v j c j 0 , j J , (43b) vec ( A ¯ j ) , 2 a ¯ j μ j , j J , (43c) v j I j × j + A ¯ j 0 , v j 0 , j J ,
where μ j and v j represent a relaxation variable. We have transformed the second non-convex constraint of formula (31) into the convex constraints 43a, 43b, 43c, so that problem (31) can be effectively solved by the convex optimization tool: CVX.

5. Simulation Results

In this section, the performance of the proposed location information-assisted method is verified by simulation, and the effect of the number of antennas, SINR threshold, and room size on the proposed robust power allocation is analyzed. Assuming that J = 6 , N = 5 , K = 4 , the reflection coefficient Γ R = 0.7 . The probability of the constraints are set to 0.05 , i.e., ρ j = 0.05 , j .
The SINR generated by 5000 random channels is summarized as a cumulative distribution function (CDF) in Figure 3. It can be seen in Figure 3 that the CDF rate of the non-robust power distribution cannot guarantee the minimum rate requirement, while the proposed robust power distribution satisfies the minimum rate constraint. In the non-robust beamforming method, there is no way to always assign extra power to users to suppress interference from other users; therefore, it is difficult for SINR to maintain at a certain level. However, the robust method proposed can allocate a reasonable amount of power to the users in order to improve the signal-to-noise ratio, thus the SINR value is always greater than the threshold value.
In order to demonstrate the performance of the proposed robust beamforming method, we compare the performance with the non-robust beamforming method (which ignores the estimation errors) and the worst-case scheme (which is presented in [31]). The variation curve of the average transmit power with the SINR threshold in the robust, non-robust, and worst-case is shown in Figure 4. Compared with the non-robust power allocation, robust power allocation requires more power to maintain the SINR at a stable level. In addition, compared with the worst-case robust beamforming, the proposed method yields less conservative results.
Figure 5 analyzes the situation in different scene sizes. It can be seen that as the room size becomes larger, the transmit power for both the robust and the non-robust beamforming gradually becomes larger. Since this paper assumes that the BSs are mainly arranged in the four corners of the scene, the measurement distance error between the BSs and the UE inevitably increases as the scene size increases, which leads to both robust beamforming and non-robust beamforming that require more power to cope with the channel error.
Figure 6 gives the average transmit power versus the transmit antenna number. As the number of antennas increases, the power required for the non-robust, robust, and worst-case scenarios gradually decreases. We can see that the method is more effective than the worst-case method, and therefore the method is less conservative. In addition, form the Equation (21), the CRLB of the user location is inversely proportional to the antenna number. Increasing the antennas number not only improves the accuracy of the user location but also improves the reliability of the reconstructed channel.

6. Conclusions

In this paper, we derived the CSI error bounded in the case of UE location uncertainty, with the help of the Taylor approximation, when using the estimated UE position to obtain the CSI in UWB wireless communication systems. Then, to resist the impact of the position error on the beamforming design, we proposed a robust beamforming optimization problem to minimize the total transmit power. To efficiently solve this nonconvex problem, we proposed an SDR-based algorithm and evaluated the performance of the UE in terms of SINR, scene size, and the antenna number. The simulation results verify the theoretical derivation of the CSI error bound, and the robustness of the method was verified by its comparison with the worst-case and non-robust methods. Therefore, the beamforming scheme proposed in this paper can have a good application prospect in promoting the development of location information-assisted UWB communication.

Author Contributions

Conceptualization, L.Y. and S.L.; software, Z.Z., X.F., and Z.S.; writing—review and editing, L.Y. and S.C.; formal analysis, S.L.; funding acquisition, S.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China under Grant 61771474.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors would like to thank the editors and the reviewers for their helpful suggestions and constructive comments.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Taking the derivative of (18) with respect to the unknown parameters, we can write the following equation:
μ f τ j k , m j = n = 1 N β j k , m w 1 , n s k U ( t τ j k , m ) τ j k , m , μ f R β j k , m j = j μ f R β j k , m j = n = 1 N w j , n s k U ( t τ j k , m ) ,
Therefore, the submatrix of FIM can be expressed as follows:
[ J τ j k τ j k ] u , v = 1 N 0 N 2 N s β j k , u β j k , v [ R 2 ] u , v ,
where
[ R 2 ] u , v = Δ 0 T o s k * ( t τ j k , u ) τ j k , u s k ( t τ j k , v ) τ j k , v d t = W W 2 2 W W 2 2 ( 2 π f ) 2 P ( f ) 2 e j 2 π f Δ τ u v d f ,
We can rewrite J τ j k τ j k , u as:
J τ j , k τ j , k = N s N 0 n = 1 N w j , n 2 { B H R 2 B } .
The other submatrices of (21) can be obtained similarly.

References

  1. Cui, Y.; Liu, F.; Jing, X.; Mu, J. Integrating Sensing and Communications for Ubiquitous IoT: Applications, Trends, and Challenges. IEEE Netw. 2021, 35, 158–167. [Google Scholar] [CrossRef]
  2. Fanti, A.; Schirru, L.; Casu, S.; Lodi, M.B.; Riccio, G.; Mazzarella, G. Improvement and Testing of Models for Field Level Evaluation in Urban Environment. IEEE Trans. Antennas Propag. 2021, 68, 4038–4047. [Google Scholar] [CrossRef]
  3. Decuir, J. Bluetooth Smart Support for 6LoBTLE: Applications and connection questions. IEEE Consum. Electron. Mag. 2015, 4, 67–70. [Google Scholar] [CrossRef]
  4. Naik, G.; Park, J.; Ashdown, J.; Lehr, W. Next Generation Wi-Fi and 5G NR-U in the 6 GHz Bands: Opportunities and Challenges. IEEE Access 2020, 8, 153027–153056. [Google Scholar] [CrossRef]
  5. Gezici, S.; Poor, H.V. Position estimation via ultra-wide-band signals. Proc. IEEE 2009, 97, 386–403. [Google Scholar] [CrossRef] [Green Version]
  6. Win, M.; Scholtz, R. Ultra-wide bandwidth time-hopping spread spectrum impulse radio for wireless multiple-access communications. IEEE Trans. Commun. 2000, 48, 679–691. [Google Scholar] [CrossRef] [Green Version]
  7. Zhu, X.; Yi, J.; Cheng, J.; He, L. Adapted Error Map Based Mobile Robot UWB Indoor Positioning. IEEE Trans. Instrum. Meas. 2020, 69, 6336–6350. [Google Scholar] [CrossRef]
  8. Fang, Z.; Wang, W.; Wang, J.; Liu, B.; Tang, K.; Lou, L.; Heng, C.; Wang, C.; Zheng, Y. Integrated Wideband Chip-Scale RF Transceivers for Radar Sensing and UWB Communications: A Survey. IEEE Circuits Syst. Mag. 2022, 22, 40–76. [Google Scholar] [CrossRef]
  9. Zhang, S.; Pedersen, G.F. Mutual Coupling Reduction for UWB MIMO Antennas With a Wideband Neutralization Line. IEEE Antennas Wirel. Propag. Lett. 2016, 15, 166–169. [Google Scholar] [CrossRef]
  10. Kaiser, T.; Zheng, F.; Dimitrov, E. An Overview of Ultra-Wide-Band Systems With MIMO. Proc. IEEE 2009, 97, 285–312. [Google Scholar] [CrossRef]
  11. Drobczyk, M.; Lehmann, M.; Strowik, C. EMC characterization of the UWB-based wireless positioning and communication experiment (wireless Compose) for the ISS. In Proceedings of the 2017 IEEE International Conference on Wireless for Space and Extreme Environments (WiSEE), Montreal, QC, Canada, 10–12 October 2017. [Google Scholar]
  12. Schirru, L.; Ledda, F.; Lodi, M.B.; Fanti, A.; Mannaro, K.; Ortu, M.; Mazzarella, G. Electromagnetic Field Levels in Built-up Areas with an Irregular Grid of Buildings: Modeling and Integrated Software. Electronics 2020, 9, 765. [Google Scholar] [CrossRef]
  13. Hussain, M.G.M. Principles of space-time array processing for ultrawide-band impulse radar and radio communications. IEEE Trans. Veh. Technol. 2002, 51, 393–403. [Google Scholar] [CrossRef]
  14. Wang, M.; Yang, S.; Wu, S. Beamforming of UWB pulse array and its implementation. Digit. Signal Process. 2006, 16, 333–342. [Google Scholar] [CrossRef]
  15. Zebiri, C.; Benabdelaziz, F.; Lashab, M. Beamforming of ultra wideband signals in an IEEE 802.15.3a channel environment. In Proceedings of the 2014 IEEE International Conference on Ultra-WideBand (ICUWB), Paris, France, 1–3 September 2014. [Google Scholar]
  16. Aye, S.; Ng, B.P. Introduction to a New Array Processing Concept: Orientational Beamforming. IEEE Commun. Mag. 2019, 57, 51–57. [Google Scholar] [CrossRef]
  17. Han, J.; Ng, B.P.; Er, M.H. Orientational Beamforming via a Modified RBF Neural Network for Orientational UWB Interference Rejection. IEEE Trans. Veh. Technol. 2022, 71, 2900–2913. [Google Scholar] [CrossRef]
  18. Hu, X.; Zhong, C.; Zhang, Y.; Chen, X.; Zhang, Z. Location Information Aided Multiple Intelligent Reflecting Surface Systems. IEEE Trans. Commun. 2020, 68, 7948–7962. [Google Scholar] [CrossRef]
  19. Jeong, S.; Simeone, O.; Kang, J. Optimization of Massive Full-Dimensional MIMO for Positioning and Communication. IEEE Trans. Wirel. Commun. 2018, 17, 6205–6217. [Google Scholar] [CrossRef] [Green Version]
  20. Zhou, B.; Liu, A.; Lau, V. Successive localization and beamforming in 5G mmWave MIMO communication systems. IEEE Trans. Signal Process 2019, 67, 1620–1635. [Google Scholar] [CrossRef]
  21. Parihar, A.; Lampe, L.; Schober, R.; Leung, C. Equalization for DS-UWB Systems—Part I: BPSK Modulation. IEEE Trans. Commun. 2007, 55, 1164–1173. [Google Scholar] [CrossRef]
  22. Abu-Shaban, Z.; Zhou, X.; Abhayapala, T.; Seco-Granados, G.; Wymeersch, H. Error Bounds for Uplink and Downlink 3D Localization in 5G Millimeter Wave Systems. IEEE Trans. Wirel. Commun. 2018, 17, 4939–4954. [Google Scholar] [CrossRef] [Green Version]
  23. Kulmer, J.; Leitinger, E.; Meissner, P.; Hinteregger, S.; Witrisal, K. Cooperative localization and tracking using multipath channel information. In Proceedings of the 2016 International Conference on Localization and GNSS (ICL-GNSS), Barcelona, Spain, 28–30 June 2016. [Google Scholar]
  24. Kulmer, J.; Leitinger, E.; Grebien, S.; Witrisal, k. Anchorless Cooperative Tracking Using Multipath Channel Informatio. IEEE Trans. Wirel. Commun. 2018, 17, 2262–2275. [Google Scholar] [CrossRef]
  25. DW1000 User Manual. Available online: https://www.decawave.com/sites/default/files/resources/dw1000_user_manual_2.11.pdf (accessed on 10 November 2019).
  26. Zhang, R.; Hoflinger, F.; Reindl, L. Inertial Sensor Based Indoor Localization and Monitoring System for Emergency Responders. IEEE Sens. J. 2013, 13, 838–848. [Google Scholar] [CrossRef]
  27. Kulmer, J.; Leitinger, E.; Meissner, P.; Hinteregger, S.; Witrisal, K. An Alternative Double-Sided Two-Way Ranging Method. In Proceedings of the 2016 13th Workshop on Positioning, Navigation and Communications (WPNC), Bremen, Germany, 19–20 October 2016. [Google Scholar]
  28. Bottigliero, S.; Milanesio, D.; Saccani, M.; Maggiora, R. A Low-Cost Indoor Real-Time Locating System Based on TDOA Estimation of UWB Pulse Sequences. IEEE Trans. Instrum. Meas. 2021, 70, 1–11. [Google Scholar] [CrossRef]
  29. Kay, S.M. Fundamentals of Statistical Signal Processing: Estimation Theory; Oppenheim, A.V., Ed.; Technometrics: Upper Saddle River, NJ, USA, 1993; p. 50. [Google Scholar]
  30. Shen, Y.; Win, M.Z. Fundamental Limits of Wideband Localization— Part I: A General Framework. IEEE Trans. Inf. Theory 2010, 56, 4956–4980. [Google Scholar] [CrossRef]
  31. Huang, Y.; Li, Q.; Ma, W.; Zhang, S. Robust Multicast Beamforming for Spectrum Sharing-Based Cognitive Radios. IEEE Trans. Signal Process. 2012, 60, 527–533. [Google Scholar] [CrossRef] [Green Version]
Figure 1. The considered UWB communication system model.
Figure 1. The considered UWB communication system model.
Symmetry 14 01171 g001
Figure 2. Illustration of multipath propagation using an image-source model.
Figure 2. Illustration of multipath propagation using an image-source model.
Symmetry 14 01171 g002
Figure 3. CDF rate of over 5000 random channel realizations.
Figure 3. CDF rate of over 5000 random channel realizations.
Symmetry 14 01171 g003
Figure 4. Average transmit power versus the SINR γ (dB).
Figure 4. Average transmit power versus the SINR γ (dB).
Symmetry 14 01171 g004
Figure 5. Average transmit power versus the scene size.
Figure 5. Average transmit power versus the scene size.
Symmetry 14 01171 g005
Figure 6. Average transmit power versus the antenna number.
Figure 6. Average transmit power versus the antenna number.
Symmetry 14 01171 g006
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Yang, L.; Zhang, Z.; Fang, X.; Cao, S.; Shangguan, Z.; Li, S. Location Information-Assisted Robust Beamforming Design for Ultra-Wideband Communication Systems. Symmetry 2022, 14, 1171. https://doi.org/10.3390/sym14061171

AMA Style

Yang L, Zhang Z, Fang X, Cao S, Shangguan Z, Li S. Location Information-Assisted Robust Beamforming Design for Ultra-Wideband Communication Systems. Symmetry. 2022; 14(6):1171. https://doi.org/10.3390/sym14061171

Chicago/Turabian Style

Yang, Lei, Zhi Zhang, Xiao Fang, Shiyu Cao, Zhegong Shangguan, and Shiyin Li. 2022. "Location Information-Assisted Robust Beamforming Design for Ultra-Wideband Communication Systems" Symmetry 14, no. 6: 1171. https://doi.org/10.3390/sym14061171

APA Style

Yang, L., Zhang, Z., Fang, X., Cao, S., Shangguan, Z., & Li, S. (2022). Location Information-Assisted Robust Beamforming Design for Ultra-Wideband Communication Systems. Symmetry, 14(6), 1171. https://doi.org/10.3390/sym14061171

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