Abstract

This paper discusses over the air (OTA) testing for multiple input multiple output (MIMO) capable terminals with emphasis on modeling bidirectional spatial channel models in multiprobe anechoic chamber (MPAC) setups. In the literature, work on this topic has been mainly focused on how to emulate downlink channel models, whereas uplink channel is often modeled as free space line-of-sight channel without fading. Modeling realistic bidirectional (i.e., both uplink and downlink) propagation environments is essential to evaluate any bidirectional communication systems. There have been works stressing the importance of emulating full bidirectional channel and proposing possible directions to implement uplink channels in the literature. Nevertheless, there is no currently published work reporting an experimental validation of such concepts. In this paper, a general framework to emulate bidirectional channels for time division duplexing (TDD) and frequency division duplexing (FDD) communication systems is proposed. The proposed technique works for MPAC setups with arbitrary uplink and downlink probe configurations, that is, possibly different probe configurations (e.g., number of probes or their configurations) in the uplink and downlink. The simulation results are further supported by measurements in a practical MPAC setup. The proposed algorithm is shown to be a valid method to emulate bidirectional spatial channel models.

1. Introduction

Multiple-input multiple-output (MIMO) systems have demonstrated their capabilities to improve system capacity and link reliability in rich multipath environments [1]. The propagation channel and antenna design are the two key parameters that together impact MIMO system performance [2]. Over the air (OTA) testing of MIMO and diversity capable terminals has attracted great research attention, since OTA testing is the only way to evaluate the radio in a realistic way [3]. Several different OTA systems have been proposed to evaluate MIMO capable terminals, for example, the two stage systems (with the conducted/radiated option) [4, 5], reverberation chamber based systems (with/without a radio channel emulator) [6, 7], and multiprobe anechoic chamber (MPAC) systems (with two-dimensional (2D)/3D probe configurations) [8, 9]. The MPAC system is a promising candidate, due to its capability to physically reproduce arbitrary multipath environments in the laboratory.

Research work in MPAC setups has been focused on emulating downlink (i.e., communication from base station (BS) to mobile station (MS)) channel models, where the uplink (i.e., communication from MS to BS) channel is often realized through a single communication antenna and a cable [3]. As shown in Figure 1, the uplink signal is radiated by the device under test (DUT) and received by an uplink communication antenna inside the anechoic chamber. Then the uplink antenna is directly connected to the input of the BS emulator with a cable, and therefore the uplink channel is constant and does not experience any fading. However, modeling uplink channel as a static free space line-of-sight (LOS) is not realistic.

Emulating realistic bidirectional channels (i.e., both uplink and downlink) is gaining significant research interest, as it is useful for evaluating any bidirectional radio communication system [10]. Bidirectional channel emulation is mandatory, especially for device testing in multiuser interference scenarios and cognitive radio terminals testing [10, 11]. The uplink channel state information can be critical for downlink performance in some closed loop communication systems with adaptive modulation, coding, and MIMO transmission modes. In the time-division duplexing (TDD) mode, as the downlink and uplink propagation channels are reciprocal, the BS may estimate the radio channel state information from the uplink channel. Then the information is utilized in selecting downlink transmission mode or even in weighting Tx antenna signals. Also in the frequency-division duplex (FDD) case, it may be important to perform testing with realistic uplink fading, for example, in order to get typical feedback channel errors. Furthermore, to reduce the difference in peak data rate between uplink and downlink, uplink with up to 4 transmission antennas (i.e., 4 transmission antennas on MS) has been introduced in LTE-Advanced (Release 10 and beyond) [12]. Furthermore, as discussed in [13], the terms downlink and uplink are not relevant to communication scenarios utilizing ad hoc and mesh networks. It would be desirable to model realistic bidirectional propagation environments to evaluate the MIMO terminal performance.

Very few papers have addressed emulation of bidirectional channel models in MPAC setups. The concept of emulating bidirectional channel models for TDD systems in MPAC setups was briefly mentioned in [14], though no details were given. In [10, 11, 13], the importance of bidirectional channel emulation is well described for closed-loop OTA testing, OTA testing of cognitive radio terminals, and OTA testing in car to car/infrastructure scenarios. In [10], it was proposed to apply wave field synthesis principles to generate uplink channel models, though no results were repoted to validate the proposal. The main contributions of this paper lie in the following aspects:(i)It is often assumed that frequency separation of two different carrier frequencies is sufficiently large compared with the channel coherence bandwidth, and therefore fading characteristics for the frequency-separated channels are uncorrelated. This is not always correct, as demonstrated in this work using data from a channel measurement campaign. The knowledge gained in the measurements is then adopted to emulate bidirectional channels in MPAC setups for FDD systems.(ii)A general framework to emulate bidirectional channels both for TDD and FDD communication systems is proposed. The proposed technique works for MPAC setups with different uplink and downlink probe configurations.(iii)Several multiprobe anechoic chamber setups, capable of emulating realistic bidirectional channels, are proposed and evaluated. The proposed technique is demonstrated by simulation results for some representative channel models.(iv)To validate the proposed algorithm, the algorithm has been implemented in a commercial channel emulator and validated in a practical MPAC setup. A measurement setup for validation purpose is proposed in the paper as well. To the best of the authors knowledge, it is the first time that bidirectional channel models have been emulated in a practical MPAC setup for MIMO OTA testing.

The target channel models adopted for the downlink in MIMO OTA testing are the well-accepted and standardized geometry-based stochastic channel (GBSC) models [3, 15]. MIMO channel models like SCME, WINNER, and IMT-Advanced models belong to this family [16, 17]. Geometry-based modeling enables separation of channel propagation and antennas at the transmitter (Tx) and receiver (Rx) side. A GBSC model consists of multiple clusters (i.e., groups of multipath components), each of which is modeled by its cluster power, delay, nominal angle of arrival (AoA), nominal angle of departure (AoD), angle spread of arrival (ASA), angle spread of departure (ASD), and cross polarization power ratio (XPR). The GBSC models are used to model the uplink channels in the paper as well, to be aligned with the downlink channels.

In TDD communications systems, as the same frequencies are utilized for the downlink and uplink, the uplink and downlink channels are considered as reciprocal. In the FDD communication systems, the uplink and downlink channels are frequency separated. In this section, a literature review on cross-correlation of channel parameters between frequency separated channels is firstly given.

2.1. Cross-Correlation between Frequency Separated Channels

The correlation properties between frequency separated channels have attracted research attention, as they are important to the design and performance prediction of multicarrier communications [18].

2.1.1. Large Scale and Small Scale Channel Parameters

The same terminologies as in GBSC models are adopted in this paper. Large and small scale parameters are specified to generate channel coefficients [19]. Large scale channel parameters, which include shadow fading, angle spread, delay spread, and Rician factor, describe GBSC models at a macroscopic level [19]. Small scale parameters, that is, cluster parameters such as powers, delays, and angles, describe the channel at a microscopic level.

With a relatively small frequency separation with respect to the carrier frequency in LTE systems, large and small scale parameters are often assumed to be the same in frequency separated channels [20, 21]. In [22], simulation results demonstrated that power delay profiles (PDPs) on different frequency bands (with a separation up to 400 MHz) were highly correlated for both line-of-sight (LOS) and non-LOS scenarios. In [23], experimental results showed that angles of arrival (AoAs) of strong specular signals did not change over the 25 MHz frequency range. Also in [24], it was concluded that power angle spectra (PAS) are approximately identical for the uplink and downlink with a frequency separation of 68 MHz. Results showed that both PDP and PAS did not depend on carrier frequency for the LOS and non-LOS scenarios based on extensive measurement campaigns, if the carrier frequency is below 6 GHz [20, 21]. Measurement results in [25] further supported the spatial reciprocity of uplink and downlink channel with a frequency separation up to 200 MHz.

2.1.2. Fading Coefficients

It is often assumed that the frequency separation in mobile cellular systems is significantly larger than the coherence bandwidth of the propagation channel, and hence the instantaneous fast fading characteristics of frequency separated channels are uncorrelated. However, in practical frequency-diversity systems, the frequency separation of different channels could be smaller than the coherence bandwidth. For example, the minimum band gap between downlink and uplink can be as low as 5 MHz in LTE FDD [26]. Carrier aggregation technology has been introduced in LTE-Advanced, where frequency separation between two adjacent subbands (intra- or interband) can be very small [27]. Also, in some frequency-diversity systems, small frequency separation exists due to spectrum and frequency planning constraints [28]. Consequently, fast fading characteristics of frequency separated channels might be correlated, depending on the frequency separation and the channel frequency correlation function (FCF) [29]. In the literature, several techniques have been proposed to generate two or any number of frequency correlated Rayleigh fading sequences [30].

2.2. Measurement on Frequency Separated Channels

Very few measurement results were reported in the literature, with a focus on investigating fading correlation of frequency separated channels. In [29], measurement results showed that the correlation decreased rapidly in an example of Rayleigh channel, as the frequency separation increases. As discussed in Section 2.1, fast fading characteristics of frequency separated channels might be correlated, depending on the frequency separation and the propagation channel. This is demonstrated in the following using experimental data.

2.2.1. Measurement Campaign Description

A detailed description of the measurement campaign was given in [31] and hence is only outlined here. The measurements were performed inside a building in a typical microcell non-LOS urban scenario with a wideband MIMO channel sounder [32]. The channel was sampled at 60 Hz, to cope with channel variations due to the movements of the Rx antennas and other changes in the channel. During the measurement, the Rx antennas were moved along a 1 m 1 m square marked on the floor with 20 s measurement time. The measurements were performed at MHz with an effective bandwidth of MHz. Note that the Rx antennas were kept at the same orientation with respect to the environment during the measurement, to ensure that the signals received by the Rx antennas have similar impinging spatial angles, and hence the recorded channels were stationary. Note that measurement data for two representative cases, selected from a huge pool of available measurements, are presented as examples here.

2.2.2. Analysis and Results

The time-variant transfer function of the measured channel is calculated by taking the Fourier transform of the measured time-variant channel in the delay domain. An illustration of the time-variant transfer function of the measured channel is shown in Figure 2. The frequency spacing between subchannels is MHz. The measured channel is Rayleigh faded with Rician factor close to . The fading correlation among the frequency separated subchannels is shown in Figure 3, analyzed with the method suggested in [29]. The correlation generally decreases as the frequency separation between subchannels increases, as expected. The measured correlations are up to 0.8 and 0.6 for a frequency separation of 5 MHz and 10 MHz, respectively.

The cross-correlation among frequency separated subchannels for a different measurement case is shown in Figure 4. The correlation drops much faster as the frequency separation increases, compared with the results shown in Figure 3. The measured correlations are close to 0.4 and 0 for a frequency separation of 5 MHz and 10 MHz, respectively.

As demonstrated by the experimental results, fast fading characteristics of frequency separated channels might be correlated, depending on the frequency separation and the propagation environment.

2.3. Summary of Requirements

Based on the measurement results and the results reported in the literature, it can be summarized that(i)large and small scale parameters in uplink and downlink channels are highly correlated both for the TDD and FDD systems, and hence large and small scale parameters used in the downlink channels can be adopted for the uplink channels directly;(ii)fading coefficients for the uplink and downlink are fully correlated in TDD systems, while fading coefficients of downlink and uplink might be correlated in FDD systems, depending on the frequency separation, and the channel FCF. Realistic uplink channel models should be reciprocal to downlink channel models in TDD systems. More specifically, we have (1)identical temporal characteristics, for example, Doppler spectrum and temporal correlation.(2)spatial reciprocity, that is, the same continuous PAS is selected as the target for the uplink and downlink;(3)fully correlated uplink and downlink fading CIRs experienced by the DUT, that is, . We havewhere and are the correlation operator and covariance operator, respectively.

As for the FDD case, the difference is that the fading CIRs in uplink channel might or might not be correlated with those in the downlink channel, depending on the downlink and uplink frequency separation and the channel coherence bandwidth. That is, it would be desirable to create the uplink and downlink fading CIRs with an arbitrary target cross-correlation :

3. Emulating Bidirectional Channels in MPAC Setups

In this section, the techniques for emulating bidirectional channels in MPAC setups are presented. First, the technique to emulate downlink channel models is revisited. After that the techniques to emulate uplink channel models are proposed. Since multicluster channel models consist of independently generated single clusters, the discussion in this section is limited to single cluster channel models. The proposed technique can be applied for multicluster channel models, simply by repeating the same procedure for each cluster [8, 33].

3.1. Downlink Channel Modeling

The prefaded signal synthesis (PFS) technique, which was proposed in [8], has been widely used for emulating downlink channel models in MPAC setups. As shown in Figure 1, each probe is connected to a channel fading module in the channel emulator, where an arbitrary fading channel can be generated. The objective is to reproduce channels with the desired temporal and spatial characteristics, with the help of the channel emulator and multiple probes.

Assume the number of probes utilized for downlink is and the channel impulse response (CIR) for the th probe for downlink is denoted by with , where is the power weight vector for the downlink probes. are the unweighted complex CIRs. is the length of the CIR. Then the CIR experienced by a DUT in the downlink, assuming an isotropic antenna pattern, isNote that the propagation coefficients from the OTA probes to the test area center are neglected, as they are the same after calibration. The unweighted fading CIR for th probe, for a specified Tx antenna and specified polarization, can be simplified according to [8], as follows: where , denote the complex constant amplitude and the Doppler frequency for the th subpath of the cluster, respectively. denotes the random initial phase for the th subpath for th probe. is the subpath index. The unweighted fading CIRs transmitted from probes belong to Class II stochastic channel simulator as detailed in [34]. That is, each OTA antenna, contributing to the same cluster, transmits independent fading CIRs with identical statistics [8]. As is formed as a linear combination of contributions from multiple probes, it has identical statistics as well.

With the PFS technique, the desired channel temporal characteristics, for example, fading distribution, temporal correlation, and Doppler spectrum, can be reproduced even with a single probe. Multiple probes around the DUT are utilized only to reproduce the spatial characteristics of the channel models at the Rx side, which are crucial for MIMO capable terminals. To reproduce the channel spatial characteristics in the test zone at the Rx, appropriate power weights are allocated to the OTA probes. The goal is to minimize the deviation between the theoretical spatial correlation resulting from the target continuous PAS and the emulated spatial correlation resulting from the discrete PAS, with its shape characterized by the discrete angular positions of the probes and the power weights [33]. The weighted CIRs for a Laplacian shaped spatial cluster with azimuth angle of arrival (AoA) and azimuth spread (AS) are shown in Figure 5. To approximate the target continuous PAS, a dominant power weight is allocated to probe 1 (with angular location ), while a very small (close to 0) power weight is allocated to probe 5 (with angular location ). The discrete PAS follows the target PAS well. The optimization technique to obtain the power weights is detailed in [33] and not repeated here.

3.2. Uplink Channel Modeling
3.2.1. Probe Configuration

An example of the multiprobe anechoic chamber setup that is capable of emulating both uplink and downlink channels is shown in Figure 6. The goal here is to emulate the realistic uplink channels experienced by the signals transmitted from the DUT, with the help of probes and the channel emulator in the uplink. A mapping of the uplink channel to each of the uplink probes is, therefore, required for adding realistic uplink fading.

In this paper, symmetric probe configuration indicates that the same set of probes are used in the uplink and downlink, that is, colocated probes for the uplink and downlink. With a symmetric probe configuration, one practical constraint is that the circulators to separate downlink and uplink signals (in Figure 6) are quite limited in frequency range. Moreover, circulators are practically limited in isolating the uplink and downlink (e.g., with an isolation loss of 20 dB).

An alternative is to use different (spatially separated) probes to separate the downlink and uplink. Asymmetric probe configurations, where different number of probes or spatially separated probes are adopted for the uplink and downlink, might be more popular in practical setups. The test zone is a geometrical volume inside which the DUT is located and it is required that the target channel models are reproduced accurately inside the test zone. It has been demonstrated that the test area size depends directly on the number of probes [33]. There might be fewer uplink probes by design if a subscriber terminal has fewer transmitting antennas than receiving antennas and the smaller transmitting antenna array size requires fewer uplink probes for the same accuracy of OTA testing. Furthermore, downlink probe configuration has been specified in standardization (e.g., 3GPP), while uplink probe configuration has not been defined.

3.2.2. Objective

Assume that we have probes for the uplink and the CIR for the th probe is denoted by with , where is the power weight vector for uplink probes. are unweighted complex CIRs. Then the CIRs experienced by the DUT in the uplink, assuming isotropic antenna pattern and neglecting propagation coefficients, are

In the following sections, techniques to emulate the uplink channel models for different probe configurations are detailed for both TDD and FDD systems. The goal is to find appropriate fading CIRs and power weights for the uplink, so that the requirements discussed in Section 3.2.2 are fulfilled for FDD and TDD systems.

3.2.3. Symmetric Probe Configurations

(a) TDD Systems. Channel reciprocity for TDD systems in symmetric probe configurations can be achieved simply by reusing the downlink power weights and fading CIRs in the uplink, that is, and with .

(b) FDD Systems. As explained earlier, in typical systems the frequency separation is high with respect to the channel coherence bandwidth, and it is a valid assumption that the uplink and downlink fading are independent and channel realizations should be uncorrelated. Anyhow there might be cases, where frequency separation is small and/or the FCF decreases very slowly, resulting in correlated uplink and downlink channels.

The procedure to create correlated uplink and downlink CIRs for symmetric probe configurations is detailed below.

Firstly, two sets of unweighted CIRs and , each CIR with independent and identical distribution, are generated according to (4). Note that the created CIRs are with identical temporal characteristics.

Secondly, the CIRs for the downlink and uplink with desired correlation can be created by multiplying the uncorrelated fading CIRs and with Cholesky factorization of the correlation matrix, according to [35], aswhere with being the complex conjugate operator and . is the Cholesky factorization operator [35].

As the probes are colocated, the spatial reciprocity can be achieved simply by setting with . The generated fading CIRs for the uplink are with desirable temporal characteristics and correlated with the downlink CIRs with a desired correlation value .

3.2.4. Asymmetric Probe Configurations

The methods proposed in Section 3.2.3, however, only work for the cases when the uplink and downlink probes are colocated (i.e., symmetric probe configurations). In this section, a novel technique is proposed for asymmetric probe configurations. In the following, is assumed for notational simplicity, although the proposed technique is not limited to this assumption.

As explained in Section 3.2.2, to obtain the spatial reciprocity, the same target PAS should be used when calculating the power weights in the uplink and downlink. The optimization techniques to obtain probe power weights for a specific probe configuration and target PAS were detailed in [33]. An illustration of the power weights for the downlink and uplink channel model for an asymmetric probe configuration is shown in Figure 7. It is expected that the emulation accuracy in terms of spatial correlation in the uplink is worse, as few probes are utilized. However, the required test area size for the uplink is also smaller, as fewer Tx antennas than Rx antennas at the MS side are expected.

The procedure to generate correlated uplink and downlink CIRs with the same temporal characteristics for asymmetric probe configurations is detailed as follows.(1)Divide downlink fading CIRs into subgroups so that each uplink CIR has an associated subset of downlink CIRs. Assume that the number of downlink fading CIRs in th subgroup is , with and . Note that there is no restriction on how downlink CIRs are divided into subgroups. th () CIR in th subgroup corresponds to the downlink CIR with in downlink CIRs.(2)Generate a set of unweighted CIRs with independent and identical distributions, according to (4).(3)Obtain uplink CIRs by expressing the uplink CIR in terms of the downlink CIRs belonging to th subgroup and th unweighted CIR . This involves expressing in the formwhere and are the complex coefficients to be calculated. th uplink CIR may then be defined in terms of the downlink CIRs belonging to th subgroup and th complex CIR . The objective is to obtain and , such that (2) is fulfilled. As explained earlier, can be obtained based on the uplink probe configuration and the same target PAS used for the downlink channel. The target correlation between the uplink and downlink fading CIRs can be calculated by methods known in the art. An example is through the FCF of radio channels, which can be obtained via Fourier transform of the channel average PDP [29]. As indicated above, for TDD systems and a predefined for FDD systems are targeted.(4)Based on the properties of and , we have and is calculated by solving the equation,where is a set of coefficients that satisfy . Note that we have coefficients to solve with only one equation. Any set of satisfying the above equation can be used. It can be easily checked that the fading correlation between uplink channel and downlink channel can reach the target correlation . A simple illustration of the procedure to generate the desired uplink CIRs (the downlink and uplink probes are selected to be and , resp.) with an example case is shown in the following.

Split downlink Irs into two subsets:

Select the appropriate coefficients :

Calculate the appropriate coefficients :

Obtain uplink Irs:

3.2.5. Discussions

The maximal achievable correlation between the uplink and downlink depends on the coefficients , uplink power weights , and downlink power weights and how the downlink CIRs are divided into subgroups. From (9), we have , and hence the maximal achievable fading correlation is

Note that symmetric probe configurations are a special class of assymetric case, where we have , , , and for . is reached when is selected. That is, the same set of CIRs for the downlink is reused for the uplink. Arbitrary correlation can be achieved by selecting appropriate values for . For example, uplink and downlink fading CIRs can be uncorrelated by setting for

From (14), we can conclude that, for an asymmetric probe configuration, can only be achieved when the PAS is uniform (i.e., for and for ), subgroups are equal (i.e., for ), and is set. For nonuniform PAS cases, the achievable correlation is upbound according to (14) and will be slightly smaller than 1.

4. Simulation Results

Three representative scenarios are selected in the simulations to demonstrate the proposed technique, as detailed in Table 1. Note that downlink CIRs are evenly divided into subgroups in all three scenarios.

4.1. Spatial Characteristics Reciprocity

As explained earlier, to ensure spatial reciprocity between uplink and downlink channels, the same target continuous PAS is selected for the downlink and uplink. The emulation accuracy in terms of spatial correlation depends directly on number of probes, as explained in [33, 36]. The spatial correlation is used as a measure to determine how well the spatial characteristic is reproduced. The spatial correlation emulation accuracy for the downlink and uplink for the scenario A is shown in Figure 8. The target curve follows zeroth-order Bessel function of the first kind, as the target PAS is uniform. With this asymmetric probe configuration, the maximal antenna separation in transmitting array on the MS is less than for the uplink, while the antenna separation in the receiving array on the MS can be as large as . Similar conclusions can be drawn for the other two scenarios. The discrete PAS for the scenario B is shown in Figure 7.

4.2. Temporal Fading Characteristics Reciprocity

As explained earlier, is formed as a linear combination of contributions from multiple probes, and hence it has identical fading statistics to that of the CIRs from individual probes. As shown in Figure 9, the temporal autocorrelation of the generated downlink CIR and the unweighted fading CIR from 1st probe match very well. The same results applies for the uplink CIR and , as shown in Figure 9. As both and are generated according to (4), the temporal autocorrelation of matches that of well. The same conclusions can be drawn for the other two scenarios. As a conclusion, the temporal fading characteristic reciprocity in the uplink and downlink can be achieved with the proposed technique.

4.3. Fading Correlation between the Downlink and Uplink Channels

To achieve the maximal correlation between the downlink and uplink fading channel, for is selected for three representative scenarios (i.e, for scenario A, for scenario B, and for scenario C). The achieved correlations between the downlink and uplink fading CIRs and , calculated based on correlation definition, are , , and for scenarios A, B, and C, respectively. The values match perfectly the maximal achievable correlation specified in (14). As we can see, though cannot be reached for nonuniform PAS channel models (e.g., scenario B and scenario C) in asymmetric probe configurations, the achieved correlations are very high (close to 1).

Cross-correlation among the fading CIRs in the downlink and in the uplink for the scenario B is shown in Figure 10. We can see that(1)correlation between the th CIR in the uplink and the downlink CIRs in th subset is , which is for the scenario B;(2)correlation between and the downlink CIRs that are not in th subset is 0;(3)correlation among uplink CIRs is 0;(4)correlation among downlink CIRs is 0. The results are expected from the equations shown in Section 3.2.4.

Note that the generated uplink CIR can be made correlated with the downlink CIR with an arbitrary value with its magnitude ranging from to the maximal achievable limit, for example, 0.89 for scenario C.

5. Measurement Validation in the Anechoic Chamber

A measurement campaign was carried out in a practical setup at Anite Telecoms Oy, Finland, to verify the algorithms presented in Section 3. Figure 11 illustrates the test setup. Eight uniformly located probes on the azimuth ring were connected to an Anite Propsim F32 to synthesize the downlink channel, while four probes were connected to the channel emulator to generate the uplink channel. Note that a channel emulator that is capable of separating uplink and downlink signal (i.e., bidirectional fading modes) is required. Phase and amplitude calibrations are performed for each probe before the measurements with a calibration antenna. The measured spatial correlation and temporal correlation of downlink channels have been well investigated in the literature [3] and not detailed here. The focus of the measurement section is to investigate whether the target fading correlation can be achieved between the downlink and uplink channel. The basic idea is to record the channel impulse responses using the calibration dipole and spectrum analyzer (SA) in the downlink and uplink separately and then calculate the measured fading correlation.

A uniform PAS is selected in the validation measurements for the sake of simplicity. Both TDD and FDD systems with target fading correlation were considered in the measurement campaign. The measured fading correlations match the target values quite well, 0.99 and 1 for the TDD system and FDD system, respectively. The deviation in the TDD system might be due to the system uncertainties in the measurement system. The measured fading profiles for the uplink and downlink for the FDD system are shown in Figure 12. Note that an arbitrary fading correlation ranging from 0 to 1 can be achieved for the FDD system in this case, as explained earlier.

6. Conclusion

This paper addresses techniques to emulate bidirectional channel models in MPAC setups. A general framework for emulating bidirectional spatial channels for both TDD and FDD communication systems is proposed for arbitrary uplink and downlink probe configurations. Large scale and small scale parameters are modeled as identical in uplink and downlink channels, based on an extensive literature review. Fading coefficients are modeled as identical as well for TDD systems to achieve channel reciprocity. For FDD systems, there might be cases, where frequency separation is small and/or the FCF decreases very slowly, resulting in correlated uplink and downlink channels. The proposed algorithm, which basically attempts to generate uplink channels based on downlink channels in the MPAC, works well in both simulations and measurements. Results have demonstrated that, for a symmetric downlink and uplink probe configuration, channel reciprocity (i.e., fully correlated uplink and downlink channel) for TDD systems and desired arbitrary fading correlation between uplink and downlink channel can be achieved for FDD systems, whereas, for asymmetric probe configurations, the maximal achievable fading correlations between uplink and downlink channels, though generally high, are limited.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

This work has been supported by the Danish National Advanced Technology Foundation via the Virtuoso project. The authors also appreciate the assistance from colleagues in Anite Telecoms Oy, Finland, especially Lic.Tech. Jukka-Pekka Nuutinen, now with Intel Mobile Communications.