A measurement method of fifth-generation multiple-input multiple-output antenna based on microwave imaging

Abstract

An increase in the quantity and density of antenna elements increases the mismatched failure rate and measurement difficulty of the multiple-input multiple-output. To simplify the measurement method of the S11 parameter utilizing the traditional vector network analyzer, this article proposes a multiple-input multiple-output measurement method based on microwave imaging. The multiple-input multiple-output element was designed, and then the existence of mismatched scattering of the mismatched state through microwave one-dimensional and two-dimensional imaging simulations was verified. A wideband Vivaldi antenna was designed for measurement imaging verification. The research results show that the proposed method is capable of detecting the mismatched scattering of mismatched elements as well as accurately locating the mismatched elements and mismatched position of circuits behind the element, which improves the measurement efficiency.

Keywords

5G IoT MIMO measurement microwave imaging mismatched scattering

Introduction

In recent years, the volume of mobile data services has grown exponentially because the Internet of Things era is developing and expanding. At the same time, with an increasingly large proportion of energy consumption in information technology systems under the fifth-generation (5G) Internet of Things (IoT) scenario, reducing the energy consumption of mobile communication network systems has gradually become a vital goal for the development of mobile communication. The fourth-generation (4G) mobile communication system will be difficult to meet the future mobile communication requirements for spectrum efficiency and energy efficiency, and to solve these problems, the 5G mobile communication system grows up.¹

The use of massive multiple-input multiple-output (MIMO) technology is a basic way to excavate wireless space dimensional resources and improve spectrum efficiency and power efficiency. However, the power consumption produced by massive MIMO will increase significantly. The single station power consumption of a 5G base station is about 2.5–3.5 times that of a single 4G station. The full load power is nearly 3800 W, and the increase in active antenna unit power consumption is the main reason for the increase in 5G power consumption.^2–7 As the quantity of massive MIMO elements gradually increases, the feed circuits also tend to be denser and accumulate heat more readily. Microwave devices, such as connectors and transmission lines, are more sensitive to temperature and are prone to mismatch with antennas. There are many array elements in 5G MIMO, up to 256 × 256,^8,9 which bring tremendous inconvenience to the MIMO measurement.

As a significant part of the entire communication system, it is necessary to separately measure antenna. The research^10–16 used over-the-air (OTA) to measure the MIMO antenna in the near field. The performance parameters of the antenna were obtained.

In addition, some studies have analyzed antennas from their scattering field. Shrestha et al.¹⁷ deduced the received signal formula of the antenna scattering and verified it using a microstrip antenna. Jiang analyzed the single-port antenna scattering theory using the angles in the time and frequency domains and intuitively distinguishes between structural and mode term scattering in the time domain.¹⁸ Moya et al.¹⁹ utilized antenna scattering methods to study the biological crystal and liquid concentration. Microwave one-dimensional (1-D) imaging methods were used to study the Doppler effect of unmanned aerial vehicle wings.^20–24 Microwave two-dimensional (2-D) imaging extraction method was utilized to obtain the antenna scattering in the real state of the fighter.^25,26

In the traditional measurement of MIMO antennas and coaxial line circuits, the vector network analyzer (VNA) is used to detect the S11 parameters of the coaxial line and the antenna one by one, and then the mismatched unit is determined according to the measured S11 parameters. The traditional method is not only inefficient but also unable to determine the position of the mismatched coaxial circuit, and it is often necessary to replace the entire coaxial cable as a whole. The method based on microwave 1-D and 2-D image detection proposed in this article is able to image the MIMO antenna and the circuit as a whole and quickly determine the specific mismatch unit of the MIMO antenna and the position of the mismatched coaxial line.

Based on the aforementioned research background, the rest of this article is organized as follows. The “System model” section introduces the scenario of the proposed measurement method. The “Proposed method” section deduces the basic theory about microwave and imaging and proposes a measurement method. The “Simulation” section presents the simulation results of the proposed method. The “Measurement” section states the microwave chamber measurement results. Finally, the “Conclusion” section concludes the article.

System model

As is shown in Figure 1, 5G signals have the characteristics of severe attenuation and narrow beams and provide the IoT service to distinct users. As a consequence, the deployment of base stations is denser and larger than 4G base stations. For instance, bus stops, office corridors, and street lights are seemingly equipped with 5G base stations. Besides, there are many array elements in 5G massive MIMO, up to 256 × 256. The matching and transmission circuit also tends to be dense, consumes a large amount of energy, and easily accumulates heat. Microwave devices, such as connectors and transmission lines, are sensitive to temperature and are prone to mismatch with antennas. This brings great inconvenience to the measurement. The traditional measurement method is to test the single MIMO element S parameter one by one to ensure whether elements match well with the transmission lines by using VNA, as shown in Figure 2. However, this measurement method is unacceptable under the condition that MIMO contains up to 256 × 256 for the reason that it costs much time. In this article, based on the above scenario, a new measurement method based on microwave imaging is proposed in Figure 3. The MIMO is the device under test (DUT), and it is connected with the matching circuit and transmission lines. The radar is the antenna that transmits the plane wave to illuminate the (MIMO) and detect the matching state of the MIMO array elements by the method of microwave imaging.

Figure 1.

Application scenario of MIMO in 5G IoT.

Figure 2.

Traditional measurement method in MIMO elements.

Figure 3.

Proposed measurement method in MIMO elements.

Proposed method

Microwave imaging theory

The basic imaging method is investigated by the previous scholar and shown.^27–29 However, some specific interpretations about imaging theory are still needed: Figure 4 demonstrates the 1-D imaging schematic, and the single frequency incident wave $u_{T} (t)$ can be defined as follows in equation (1):

u_{T} (t) = e^{j 2 π ft}

(1)

Figure 4.

1-D imaging schematic: (a) the side view of 1-D imaging schematic and (b) the top view of 1-D imaging schematic.

The incident signal that irradiates with a time delay and is weighted by the reflection intensity of the received signal $u_{R} (t)$ in the radar antenna is shown as follows in equation (2)

u_{R} (t) = \int_{0}^{\infty} g (R) e^{j 2 π f (t - \frac{2 R}{c})} dR

(2)

The eventual received signal $G (f)$ in the VNA is obtained after the frequency mixing between the $u_{T} (t)$ and $u_{R} (t)$ as follows in equation (3)

G (f) = \int_{0}^{\infty} g (R) e^{- j 2 π f (\frac{2 R}{c})} dR

(3)

According to equation (3), the target reflectance $g (R)$ and the frequency domain received data $G (f)$ form a pair of Fourier transform pairs. After measuring the different frequency points of the test bandwidth, a series of frequency domain data are obtained. The data are recorded, and the 1-D image of the target is obtained after the inverse Fourier transform.

The 2-D imaging is more intuitive and clear than the 1-D imaging, which is able to display the inverse synthetic aperture radar (ISAR) imaging under the 2-D space with a certain resolution. It is helpful to conduct a multi-dimensional analysis of scattering in engineering. .

The process of microwave 2-D ISAR imaging is shown as follows: 2-D separated processing method is utilized; first, 1-D longitudinal distance resolution of the is obtained, and then the horizontal compression of the aclinic axis is used to obtain the horizontal distance resolution of the aclinic axis. Longitudinal distance resolution of the imaging system is achieved in broadband signal pulse compression, and horizontal distance resolution of the imaging system is achieved in synthetic aperture processing.

The basic geometric relationship of imaging is shown in Figure 5. The transmitting radar antenna is fixed at a certain position and then rotates around the center $O$ . The $u - v$ coordinate system is fixed relative to the position of the transmitting radar antenna, and its coordinate origin is the center O of the device under test (DUT). The distance $R_{0}$ to the rotation center $O$ is a constant. The $x - y$ coordinate system is different from the $u - v$ coordinate system. It is a coordinate system that is fixed on the DUT. The origins of the $x - y$ coordinate system and the $u - v$ coordinate system coincide, and they are both located at point $O$ .

Figure 5.

2-D imaging schematic: (a) the top view of 2-D imaging schematic and (b) the basic geometric relationship of 2-D imaging.

According to the conditions described earlier, the signal transmitted by the transmitting radar antenna is a single frequency signal, where the frequency is $f$ , so the total received signal $G_{t} (f, θ)$ can be demonstrated by equation (4)

\begin{matrix} G_{t} (f, θ) = \int \int_{- \infty}^{\infty} g_{θ} (u, v) e^{- j 4 π fR / c} dudv \\ = \int \int_{- \infty}^{\infty} g (x, y) e^{- j (\frac{4 π f}{c}) (R_{0} + ycos θ - xsin θ)} dxdy \end{matrix}

(4)

The frequency-domain data $G (k, θ)$ is derived as follows after assuming the $k = 2 f / c$

G (k, θ) = \int \int_{- \infty}^{\infty} g (x, y) e^{- j 2 π k (ycos θ - xsin θ)} dxdy

(5)

According to equation (5), the target reflectance $g (x, y)$ and the frequency domain received data $G (k, θ)$ form a pair of 2-D Fourier transform pairs. After measuring the different frequency and different angle points of the test bandwidth, a series of frequency domain data are obtained. The data are recorded, and the 2-D image of the target is obtained after the inverse Fourier transform.

High-frequency signal transmission model

Transmission line is the device to guide and transmit electromagnetic wave. It is the connection between the transmitting antenna and the transmitter or between the receiving antenna and receiver. According to the traditional theory, high-frequency signals are transmitted on the transmission line from the signal source to the load under the condition that they are in impedance matching state. Otherwise, part of the high-frequency signals will be reflected back to the signal source, resulting in energy loss, as shown in Figure 6.

Figure 6.

Transmission line model.

The reflection coefficient $τ$ denotes the ratio of the reflected voltage and the incident voltage on the transmission line. The definition is as follows

τ = \frac{Z_{L} - Z_{0}}{Z_{L} + Z_{0}}

(6)

$Z_{0}$ is the impedance of transmission line, and $Z_{L}$ is the impedance of the load. $τ$ is $0$ when $Z_{0}$ is equal to $Z_{L}$ , which denotes that the transmission line matches well with the load. However, the absolute value of $τ$ is $1$ when the load is in a short state ( $Z_{L}$ is equal to $0$ ) and the load is in an open state ( $Z_{L}$ is equal to $\infty$ ) because of equation (6)

lim_{Z_{L} \to \infty} | τ | = lim_{Z_{L} \to \infty} \frac{Z_{L} - Z_{0}}{Z_{L} + Z_{0}} = lim_{Z_{L} \to 0} | τ | = 1

(7)

Equation (7) also indicates two extreme conditions when the signal is reflected overall from the load to signal source.

MIMO element scattering model

Antenna is a special kind of measured object. Due to its reciprocity, it is able to both radiate and receive high-frequency electromagnetic signals. According to the antenna scattering characteristic, the scattering of the antenna is divided into two parts: the extra and adjoint item scattering and the mismatched item scattering as shown in equation (8). Total antenna scattering field is demonstrated as the $G_{antenna} (f, θ)$

G_{antenna} (f, θ) = \underset{G_{st, i}}{\sum_{i = 0}^{\infty}} (f, θ) + \underset{G_{md, n}}{\sum_{n = 0}^{\infty}} (f, θ) e^{j φ_{n}}

(8)

The first part in the right equation (8) is called the extra and adjoint item scattering or structural item scattering, which has nothing to do with the transmission line circuit but the outer shape. The mismatched scattering is related to the transmission line circuit because the match state between the antenna and the transmission line would determine the reflection coefficient $τ$ , and it will generate the second radiate field as shown in Figure 7. There is still the phase difference $ϕ_{n}$ between two parts.

Figure 7.

Measurement schematic diagram.

As is shown in equation (5), the antenna reflectance $g_{antenna} (x, y)$ , which is also the target image and the frequency domain received data $G_{antenna} (f, θ)$ form a pair of 2-D Fourier transform pairs. After measuring the different frequency and different angle points of the test bandwidth, a series of frequency domain data are obtained. The data are recorded, and the 2-D image of the target is obtained after the inverse Fourier transform as shown in equation (9)

\begin{matrix} g_{antenna} (x, y) = \sum_{i = 0}^{\infty} g_{st, i} (x, y) \\ + \sum_{n = 0}^{\infty} g_{md, n} (x + Δ x, y + Δ y) \end{matrix}

(9)

In equation (9), $g_{st} (x, y)$ represents the antenna reflectance of the received structural item wave and $g_{md} (x + Δ x, y + Δ y)$ represents the antenna reflectance of the received mismatched item wave. However, according to the time shift characteristic of Fourier transform, there is the position shift in the image $Δ x$ and $Δ y$ if the MIMO elements are not capable of matching well with the transmission lines, which is the main method proposed in this article to detect the MIMO elements.

Simulation

According to the antenna scattering theory, when the element antenna matches well with the transmission line, mismatched scattering is 0 because the energy inputs from the antenna are absorbed by the transmission line and circuits. There is only the extra and adjoint item scattering or structural item scattering left. As a consequence, the scattering model of the matched antenna is replaced by the same dimension perfect electric conductor, which owns the same scattering mechanism as the matched element antenna in the equations (8) and (9). In this section, both the mismatched element physical model and equivalent matched element physical model are constructed in the FEKO software and processing the data from the FEKO in MATLAB. Then their microwave 1-D images and 2-D images are obtained and shown.

MIMO element 1-D imaging verification

The FEKO software is used to construct the MIMO microstrip mismatched element (without connection with matched transmission line or load). Its dimension is 50 mm × 80 mm × 3 mm, which is working at the frequency of 2.4, as shown in Figure 8.

Figure 8.

Mismatched element diagram: (a) mismatched element FEKO model and (b) mismatched element dimension diagram.

The equivalent matched element physical model is also constructed in the FEKO software as shown in Figure 9.

Figure 9.

Equivalent matched element model diagram: (a) equivalent matched element model FEKO model and (b) equivalent matched element dimension diagram.

Using the incident plane wave which is vertical to mismatched element and equivalent scattering body plane to irradiate them respectively. The incident direction is as shown in Figure 10.

Figure 10.

The incident direction diagram: (a) Incident wave direction of the mismatched element and (b) Incident wave direction of the equivalent matched element.

The incident wave frequency is from 2 to 6 GHz, and its interval is 100 MHz. After obtaining the electronic field data from the FEKO, the data are processed by microwave imaging method in MATLAB and the 1-D microwave imaging of the mismatched element and the equivalent matched element model as shown in Figure 11.

Figure 11.

Comparison of mismatched and equivalent matched element 1-D images.

These are two phenomena that are obtained after comparing two 1-D imaging images. First, both the scattering intensity of mismatched element and the scattering intensity of equivalent matched element reach the peak of −25.5 dBsm and −23.5 dBsm at the physical position 0 m on the distance axis, which demonstrates that this is the structural term scattering. Second, the scattering intensity of the mismatched element is 7 dBsm higher than the equivalent matched element from 0.1 m to 0.3 m, which denotes that this is the mismatched term scattering and phenomena also correspond to scattering model. Besides, the difference between the highest peak in both two curves can be explained by the energy conservation law.

However, it is incapable of detecting the match state of the horizontal distributed MIMO when there is more than one element in the same longitudinal distance because the 1-D image will alias together in the same longitudinal distance. The second dimension (horizontal dimension) is needed to detect a distinct match state of elements.

MIMO element 2-D imaging verification

To verify the viability of the proposed method under the 2-D imaging, in this section, a 1 × 2 MIMO that contains a mismatched element and an equivalent matched element is constructed in the FEKO software as shown in Figure 12.

Figure 12.

The 1 × 2 MIMO model diagram: (a) top view of the MIMO diagram and (b) side view of the MIMO diagram.

Using the incident plane wave from different angles plane to irradiate the MIMO and the MIMO is always vertical to the incident wave plane. The incident angle $θ$ is 76.4° which is from −38.2° to +38.2°. The incident angle interval is 3.63°. Compared with the 1-D imaging incident wave, the 2-D imaging incident wave frequency is from 2 to 10 GHz, and its interval is 190.4 MHz. The specific incident schematic is shown in Figure 13.

Figure 13.

(a) The 2-D imaging incident schematic of 1 × 2 MIMO and (b) 2-D imaging of MIMO.

The 2-D imaging of the MIMO is obtained after processing the data exported from FEKO software in the MATLAB as shown in Figure 13. As shown in Figure 13, there are two strong scattering areas on the mismatched and equivalent matched MIMO element physical positions, respectively. According to the antenna scattering model derived earlier, these two strong scattering images are the structural item scattering. However, there is still another strong scattering area behind the mismatched (left side) antenna, and there is no strong scattering area behind the equivalent matched (right side) antenna model, which corresponds to the mismatched item scattering derived earlier. As a result, the matched antenna and the mismatched antenna in the MIMO are detected by 2-D imaging.

Measurement

In this section, a specific measurement experiment is conducted. The measurement environment is a microwave anechoic chamber, and the experimental equipment includes a VNA, a signal amplifier, a transmitting/receiving antenna, the measured object is a Vivaldi antenna and a 1 m long coaxial cable, and matched/mismatched loads. The experiment first measured the 2-D imaging of the Vivaldi antenna with a matched and a mismatched load to verify the correctness of the method which is capable of detecting the mismatched element in MIMO and then connected the 1 m long coaxial cable with the matched/mismatched load to the Vivaldi antenna for imaging measurement to verify the viability of determining the mismatched location. The expected imaging in mismatched antenna emerges the strong scattering behind the antenna. On the contrary, the expected imaging in the matched antenna presents the weak or void scattering behind the antenna. Particularly, in the 1-D experiment with 1 m long cable and mismatched load, the strong scattering in the physical mismatched position contributes to verify the method feasibility.

Measurement environment and DUT

As shown in Figure 14, measurement equipment includes a pair of 2.5–18 GHz transmitting and receiving antennas and a cuboid tunnel to restrict the electromagnetic wave plane dimension, which is shown in (a). An electronic 2-D rotation platform, which is shown in (b). And a VNA in (c). An equipment connection diagram and a specific view of the measurement system are shown in Figures 15 and 16.

Figure 14.

Measurement system diagram: (a) transmitting and receiving antenna and cuboid tunnel; (b) 2-D rotation platform; and (c) vector network analyzer (VNA).

Figure 15.

Measurement system connection diagram.

Figure 16.

(a) Top view of the measurement system and (b) side view of the measurement system.

To correspond with the imaging bandwidth from 2 to 10 GHz, a wide band Vivaldi antenna is designed as the MIMO element (DUT), which is shown in Figure 17. Figure 18 shows its E-plane and H-plane normalized pattern after measured by VNA.

Figure 17.

Vivaldi antenna: (a) top view of Vivaldi antenna and (b) backplane of Vivaldi antenna.

Figure 18.

2-D normalized pattern of Vivaldi antenna: (a) E-plane pattern and (b) H-plane pattern.

2-D imaging measurement

To verify the proposed method to detect the MIMO element, a wide band Vivaldi antenna which is placed on the rotation platform aslant is measured under the microwave chamber environment introduced in the part before. Besides, the Vivaldi antenna is connected with the 50 Ω matched load to represent the matched element. Mismatched state element does not connect any load but hang. The measurement scenario diagram is shown in Figure 19, and the specific parameter is shown in Table 1. Note that the 0° sample angle is the circumstance in Figure 19(a), and the anticlockwise direction demonstrates the positive number angle and vice versa.

Figure 19.

2-D imaging diagram: (a) top view of diagram and (b) real state of the element.

Table 1.

Element 2-D imaging parameters.

Parameters	Values
Sample frequency	4.5–7.5 GHz
Frequency points	301
Sample angles	−20° to 20°
Angle points	81
x-axis range	−0.7 to 0.7 m
y-axis range	0–1.4 m

The Vivaldi element under the mismatched and matched conditions on the rotation platform is measured by VNA and the data are processed in the MATLAB software. As a consequence, the 2-D images of the element are shown in Figure 20.

Figure 20.

2-D images of element: (a) 2-D images of mismatched element; (b) 2-D images of matched element.

As shown in Figure 20, the measured element is located on the platform whose location is (0, 0). The images in Figure 20 mainly show the 2-D image behind the antenna port position in the longitudinal direction. Observing the two images above, there is a strong scattering area behind the antenna in (a). The reason is that there is a phase difference between the structural item scattering and the mismatched scattering, and the mismatched scattering will occur behind the physical position of the element, resulting in a delayed effect. Compared with Figure 20(a) and (b) scattering still occurs behind the element, but the scattering effect is lower than Figure 20(a). This is because the impedance of the element is not 50 Ω completely but approximately, which is explained by equation (6). The reflected wave in the load still occurs in the 50 Ω load. According to the aforementioned phenomenon, the 1-D mismatched and matched MIMO elements are capable of being detected by the 2-D imaging and the match state of 2-D MIMO elements is also able to be measured following the steps given below: first, measure the MIMO elements in 2-D imaging method one time to ensure the x-axis coordinate of mismatched elements, and then measure the MIMO elements in 2-D imaging method another time after rotating the elements 90° to ensure the y-axis coordinate of mismatched elements.

1-D imaging measurement

After ensuring the mismatched elements utilizing the aforementioned method, the next step is to ensure the mismatched position between the element and the transmission line. To simulate the real circumstance, in this subsection, the element connected with the transmission line and load is constructed. The 1-D imaging is used to detect the accurate mismatched position between the element and the transmission line.

To verify the proposed method to detect the mismatched position, the element connected with the 1 m transmission line, and load is constructed as shown in Figure 21 and the 1-D imaging parameters are shown in Table 2.

Figure 21.

1-D imaging diagram: (a) side view of the diagram and (b) real state of the element connected with the transmission line and load.

Table 2.

Element 1-D imaging parameters.

Parameters	Values
Sample frequency	4.5–7.5 GHz
Frequency points	301
Sample angle	0°
Angle points	1
y-axis range	0–15 m

Before measuring the element connected with the transmission line and load, standard sphere measurement is needed to ensure the coordinate of the rotation platform along the y-axis because of the existence of the circuit delay. The 1-D imaging of standard sphere schematic is shown in Figure 22, and the 1-D imaging of standard sphere parameters is the same as the element in Table 2. The 1-D image of the standard sphere is shown in Figure 22(b). According to the image, the highest scattering peak coordinate position is 8.424 m where it is also the platform or DUT coordinate position.

Figure 22.

1-D imaging diagram: (a) Side view of diagram and (b) 1-D imaging of standard sphere.

There are two states of element (DUT) in this measurement. The first state is that the element is connected with the 1 m transmission line and open circuit load which means the reflection coefficient $τ$ is 1 on the load position according to equation (7). The second state is that the element is connected with the 1 m transmission line and matched circuit load which means the reflection coefficient $τ$ is 0 on the load position according to equation (6) because the impedance of the load is equal to the transmission line. The 1-D image of the element is shown in Figure 23.

Figure 23.

1-D imaging of element.

As is shown in Figure 23, in addition to the strong scattering near the coordinate position of element, the element connected with the transmission line and mismatched load occurs strong scattering area 1 m behind the element position where there is the mismatched load position approximately. However, there is not a strong scattering area on the curve of element with the matched load, which corresponds to the element scattering theory in Proposed method section and the simulation in Simulation section. In this process, the mismatched position of the single element is detected. The 2-D imaging result and 1-D imaging result of the element are concluded in Tables 3 and 4.

Table 3.

Element without 1 m long coaxial line 2-D imaging result.

Experiment conditions	Results
Mismatched element	Strong scattering behind element
Matched element	Weak scattering behind element

Table 4.

Element with 1 m long coaxial line 1-D imaging result.

Experiment conditions	Results
Mismatched element	Strong scattering behind element
Matched element	Weak scattering behind element

Conclusion

In this research, a measurement method based on the 1-D and 2-D microwave imaging was proposed for MIMO elements. The measurement was demonstrated based on recent work including the theory derivation, FEKO and MATLAB software simulation, and element measurement in the microwave chamber. Compared with the previous studies that primarily measure the total MIMO radiation performance exhibited in the “Introduction” section, this research focused on the detailed scattering mechanism of the element. To solve the mismatching measurement problem, the 1-D and 2-D microwave imaging method is proposed under a dense 5G MIMO elements scenario. The results show that the proposed method is capable of detecting the elements and ensuring the accurate mismatched position of microwave circuits after measuring the elements in 2-D imaging two times and 1-D imaging one time, which improves measurement efficiency. In addition, measurement research under other specific mismatched conditions (such as short circuit mismatched or part mismatched) was not addressed in this research and is a valuable future contribution to this work.

Footnotes

Handling Editor: Bo Rong

Declaration of conflicting interests

The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.

Funding

The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: The work by Tao Hong is supported by National Natural Science Foundation of China under Grant No. 61827901.

ORCID iDs

Hangyu Chen

Tao Hong

References

Wang

Bai

, et al. Recent advances and future challenges for massive MIMO channel measurements and models. Sci China (Inf Sci) 2016; 59(2): 1–16.

Liu

Plets

Goudos

, et al. Multi-objective network planning optimization algorithm: human exposure, power consumption, cost, and capacity. Wirel Netw 2015; 21(3): 841–857.

Mohseni-Ejiyeh

Ashouri-Talouki

. SeVR⁺: secure and privacy-aware cloud-assisted video reporting service for 5G vehicular networks. In: 2017 Iranian conference on electrical engineering (ICEE), Tehran, 2–4 May 2017. New York: IEEE.

Qingfeng

Yindi

. Spectral-energy efficiency tradeoff in mixed-ADC massive MIMO uplink with imperfect CSI. Chin J Electron 2019; 28: 619–624.

Qiu

Chen

Luo

, et al. Downlink precoding with mixed statistical and imperfect instantaneous CSI for massive MIMO systems. IEEE T Veh Technol 2018: 67: 3028–3041.

Huang

Wang

, et al. Spectral and energy efficiency tradeoff for massive MIMO. IEEE T Veh Technol 2018; 67: 6991–7002.

Liu

Wang

Chen

, et al. Distributed energy efficient fair user association in massive MIMO enabled HetNets. IEEE Commun Lett 2015; 19: 1770–1773.

Larsson

Edfors

Tufvesson

, et al. Massive MIMO for next generation wireless systems. IEEE Commun Mag 2013; 52(2): 186–195.

Keusgen

. On limits of wireless communications when using multiple dual-polarized antennas. In: 10th international conference on telecommunications, 2003. ICT 2003, Papeete, 23 February–1 March 2003.

10.

Zhang

, et al. Comparison tests and hand phantom standardization for multi-probe based MIMO OTA. In: 2016 IEEE 5th Asia-Pacific conference on antennas and propagation (APCAP), Kaohsiung, Taiwan, 26–29 July. New York: IEEE.

11.

Wei

Lassi

Kyösti

, et al. Test zone size characterization with measured MIMO throughput for simulated MPAC configurations in conductive setups. IEEE T Veh Technol 2017; 66: 10532–10536.

12.

Kyosti

Fan

Pedersen

, et al. On dimensions of OTA setups for massive MIMO base stations radiated testing. IEEE Access 2016; 4: 5971–5981.

13.

Hong

Jiang

, et al. Multi-beam antenna technologies for 5G wireless communications. IEEE T Antenn Propag 2017; 65: 6231–6249.

14.

Long

Ouyang

Yang

, et al. Calibration method of phased array based on near-field measurement system.

15.

Gao

Wang

Fan

, et al. Beam probability metric for OTA testing of adaptive antenna systems in multi-probe anechoic chamber setups. In: 2019 13th European conference on antennas and propagation (EuCAP), Krakow, 31 March–5 April 2019, vol. 4, pp.1–5. New York: IEEE.

16.

Llorente

Fan

Pedersen

. MIMO OTA testing in small multi-probe anechoic chamber setups. IEEE Antenn Wirel Pr 2016; 15: 1167–1170.

17.

Shrestha

Balachandran

Agarwal

, et al. A method to measure radar cross section parameters of antennas. IEEE T Antenn Propag 2008; 56(11): 3494–3500.

18.

Wang

Zhang

, et al. Radiation and scattering properties of the dual-port fed differential microstrip antennas. In: 2011 4th IEEE international symposium on microwave, antenna, propagation and EMC technologies for wireless communications, Beijing, China, 1–3 November 2011, vol. 1, pp.124–127. New York: IEEE.

19.

Moya

Archilla

Navarro

, et al. Scattering of microwaves by a passive array antenna based on amorphous ferromagnetic microwires for wireless sensors with biomedical applications. Sensors 2019; 19(14): 3060.

20.

Zhao

Yang

, et al. UAV detection and identification in the internet of things. In: 2019 15th international wireless communications & mobile computing conference (IWCMC), Tangier, Morocco, 24–28 June 2019, vol. 1, pp.1499–1503. New York: IEEE.

21.

Giyenko

Cho

. Intelligent UAV in smart cities using IoT. In: 2016 16th international conference on control, automation and systems (ICCAS), Gyeongju, South Korea, 16–19 October 2016. New York: IEEE.

22.

Mozaffari

Saad

Bennis

, et al. Mobile unmanned aerial vehicles (UAVs) for energy-efficient internet of things communications. IEEE T Wirel Commun 2017; 16(11): 7574–7589.

23.

Chen

, et al. Micro-doppler effect in radar: phenomenon, model, and simulation study. IEEE T Aero Elec Sys 2006; 42(1): 2–21.

24.

Morsy

Achtenova

. Rolling bearing fault diagnosis techniques—autocorrelation. In: 2015 23rd Mediterranean conference on control and automation (MED), Torremolinos, 16–19 June 2015, pp.328–334. New York: IEEE.

25.

Chen

, et al. Radar cross section measurement of conformal antenna based on microwave imaging. In: 2017 11th European conference on antennas and propagation (EUCAP), Paris, 19–24 March 2017. New York: IEEE.

26.

Zhao

Chen

. Research on the RCS of serrate gap in real aircraft state. In: 2018 IEEE 18th international conference on communication technology (ICCT), Chongqing, China, 8–11 October 2018, vol. 4, pp.338–341. New York: IEEE.

27.

Rong

Salehian

, et al. Cloud transmission: a new spectrum-reuse friendly digital terrestrial broadcasting transmission system. IEEE T Broadcast 2012; 58(3): 329–337.

28.

Rong

Qian

, et al. Call admission control optimization in wimax networks. IEEE T Veh Technol 2008; 57(4): 2509–2522.

29.

Sun

Kadoch

Gong

, et al. Integrating network function virtualization with SDR and SDN for 4G/5G networks. IEEE Network 2015; 29(3): 54–59.