#### Abstract

User effects on portable antennas can decrease the global terminal performance, such as reducing the battery life or the wireless coverage. People presence will modify the input impedance, the far field radiation pattern, the radiation efficiency, and the near field. In this paper, we investigate the influence of the user on the input impedance by means of a new statistical approach which relies on uncertainty, where a handheld antenna works in two modes (data and talk) while people couple to the antenna. Studying some reflection coefficient measurements performed with a vector network analyzer (VNA), a standard framework based on uncertainty ellipses is developed. Finally, the method is applied to some antenna prototypes, showing its advantages and reliability to quantify and compare proximity effects.

#### 1. Introduction

Mobility and user proximity impact on handheld terminals and radiating structures has become a key issue in order to improve standard portable antenna designs. During the last 20 years, some researchers have focused on the evaluation of the human body impact [1–5], which has been mostly carried out by means of electromagnetic simulations and verified sometimes by measurements.

On the one hand, in spite of the clear general understanding of the head, hands, and even other body parts and surrounding objects on antenna circuit and field parameters, the investigations have been carried out presenting averaged ratios of usual quantities, such as the input impedance [6, 7] or the radiation pattern [8]. However, we should not forget that neither the hands nor the user’s head keeps the same position relatively to the antenna, and consequently, there is some variability around those averages [9]. It is always possible to study a few feasible configurations, but some magnitudes like the antenna input impedance will be no longer static in real scenarios. Therefore, the actual “instantaneous” input impedance will be near the obtained values for a finite set of suitable states. Even though modern computational electromagnetic techniques can deal efficiently with each state separately, calculating exactly the actual antenna radiation pattern and the input impedance might be difficult.

On the other hand, there is also a growing concern about the final effects on the system where the antenna is integrated, working in transmitting or receiving mode (e.g., the signal to noise ratio (SNR) variations due to antenna and receiver mismatching [10, 11] or the degradation on power amplifier (PA) performance). According to those works, the impact of mismatching is not negligible and must be considered so that the communication system can be properly designed. Mismatch may be even more relevant in the case of wireless technologies since there are additional sources of variability introduced by mobility and user interaction with the antenna. In fact, some solutions have been suggested to reduce mismatch, such as impedance tuning networks or reconfigurable antennas [12–14]. They are supposed to reduce the reflection coefficient dynamically by means of variable lumped or adaptive transmission line elements.

In contrast to other radiating structures, portable antennas should be characterized not only in terms of the classical static magnitudes (input impedance, radiation pattern, axial ratio, etc.) but also from the dynamical point of view. There is a lack of dynamic models in the literature which would be very helpful to quantify the variability and to develop improved portable antennas. In addition, these models could be relevant to design receivers and transmitters because suitable mitigation policies could avoid a bad reception and improve the battery life or the handheld coverage.

The proposed novel approach analyzes the impedance variability. It is not exactly a dynamical model but aims at estimating the impedance variability region on the Smith chart and is based on statistical measurements. The advantage of the method relies on its generality since it can be applied to any antenna geometry, object, and frequency band. It allows to compare portable antenna performances and, moreover, provides to the antenna and microwave engineer a procedure to build the variability region with a few parameters in order to evaluate the performance degradation or to simplify and improve impedance tuning networks or reconfigurable circuit designs.

The outline of the paper is as follows. First, we establish the connection between proximity effects and variability by means of the support of the probability density function (PDF) corresponding to the reflection coefficient. Second, we estimate the variability with the aid of the uncertainty ellipses which come naturally from the covariance matrix corresponding to the experimental measurements accomplished to study the effect of certain objects. We focus on a typical and practical mobile scenario with real antennas in the UTRA/FDD LTE B8 and UTRA/FDD LTE B3 frequency bands. The ellipse framework requires five geometrical parameters which may be estimated with certain confidence. We present the evaluation of this new method performed on three prototypes and give the meaning of those parameters in terms of the variability. Finally, we propose the 95% ellipse as its measure and some applications on reconfigurable tuning networks.

#### 2. Statistical Approach to the Impedance Variation Problem

Figure 1 shows a typical wireless terminal where the antenna is characterized as a complex input impedance from the circuit point of view and the interaction with a user in two possible modes (talk and data transmission configurations). The input impedance (at the frequency ) is denoted as , where is the resistance and is the reactance.

Antenna proximity effects can be geometrically identified on the Smith chart. Dynamic coupling transforms the point which represents the static impedance at a fixed frequency onto a set of probable impedance values which constitute a variability area. Thus, some dynamic coupling features can be analyzed by studying and quantifying this Smith chart subdomain. Identically, the input impedance curve that models the antenna frequency band behavior is also transformed into a certain area, and thus, similar techniques can be applied to broadband or narrowband problems.

##### 2.1. Variability Impedance Region through the Reflection Coefficient in the Matching Domain

Let a radiating structure be placed close to the user in a wireless environment working at the frequency . The antenna will have an input impedance at the time . Then, its instantaneous reflection coefficient is defined by where is the reference impedance and is the normalized antenna impedance. The modeling process can be carried out either on the antenna input impedance or on the instantaneous reflection coefficient, because both reflect the same concept, but in different domains. They are complex random variables in discrete time ().

In order to keep the input impedance study as simple as possible, it is important to find out which is the most convenient representation. There are four different possibilities, but the Cartesian form where the reflection coefficient is mathematically represented by its real and imaginary parts , being that and are real numbers, is preferred because the components of the random vector have identical compact support and dimensional units. In addition, it is easy to estimate the input impedance from where denotes the normalized impedance.

Rigorously, the reflection random vector which models the impedance variations is a stochastic process in discrete time. The cross-correlation, covariance, or higher order moments are typically required to describe its properties from a mathematical point of view. However, if antenna designers only need to compute the main magnitudes and do not deal with prediction issues in time, the probability density of can be accurate enough.

We define mathematically the impedance variability region (IVR), , as the closure of the set of possible values of having the density on the Smith chart space such that

The impedance variability region (whose limit is the Smith chart circle ) contains the proximity effects and depends on the antenna geometry and the nature of the interaction, hence the interest of estimating such a region with certain confidence. As mentioned, it can be applied to a single frequency or to a complete frequency band.

##### 2.2. : Statistical Estimation

A VNA can collect a set of reflection coefficient samples from an antenna every second. The measurement process is performed while the antenna is moving on a certain path or, for example, simulating a phone call with the radiating structure close to the user’s head or pocket. The coupling phenomena appear in these data and will be very different from the static values.

The sampling time (or sampling frequency) depends on the VNA capabilities and the frequency of interest. Once the vector samples are measured , the estimated impedance variability region may be calculated from the approximated probability density of .

The samples obtained after the data registration can be arranged in a complex vector which can be partitioned into a number of equal-sized cells of the area because the support of the random components is compact (). The density estimation at a point [15] is where is the number of samples in the cell, called bin. According to the approximated numerical estimation, is

For example, the impedance variations are very small when the antenna is isolated. As a result, the impedance density should tend to a Dirac distribution and the estimated impedance variability region should approximate to

#### 3. Reflection Coefficient Statistical Measurements under User Influence

##### 3.1. Method

The computation of the approximated impedance variability region begins with the experimental antenna measurement sampling (Figure 2(a)). There are many feasible situations where an impedance fluctuation analysis is important. Since the complex environment includes so many objects interacting with the antenna, we reduce the research in the present work to a set of the most probable scenarios and study the radiating structure proximity effect in two situations: the so-called “data” and “talk” modes. When the terminal is operating in a “data mode,” the antenna is held by a user close to the body to simulate the access to internet applications. When the user is performing a phone call, the geometry is held with the hand near the head and this is the “talk mode.” Besides, the user has certain mobility around to include the dynamic effect. The VNA (Anritsu 37000) captures the reflection coefficient while the person performs his tasks and registers the interaction.

**(a)**

**(b)**

**(c)**

**(d)**

As far as the frequency band is concerned, the study is carried out in the UTRA/FDD LTE B8 and UTRA/FDD LTE B3 frequency bands. This aims at showing the proposed method advantages in real applications, but the approach is applicable to other frequency bands as well. The measurement process is made by means of in-house software which samples the instantaneous reflection coefficient through a GPIB device. The maximum sampling frequency depends on the VNA and the equipment calibration. In this work, data have been collected using a 1 ms sampling time. The user dynamics is very slow compared to the signal frequency, and this sampling frequency can record most of the proximity effects.

##### 3.2. Antennas under Test, Experimental Setup, and Results

Figures 2(b)–2(d) show the test configurations to capture the reflection coefficient variability. Three different antennas [16–18] suitable for medium-size mobile devices have been designed on FR4 substrate: a single-sided fractal antenna (Cantor model, Figure 2(b)), a CPW-feed printed slot antenna (Descartes model, Figure 2(c)), and a double-sided printed antenna (T-monopole model, Figure 2(d)). The prototypes are intended to work according to the UTRA/FDD LTE B8 and UTRA/FDD LTE B3 standard. The requirement concerning the return losses is 6 dB, but they have different scattering parameters () in those frequency bands. The antennas differ in field and directivity performance because of the current paths on the conductors, in spite of the underlying radiation physics similarity (the structures are based on printed technology). Figure 3 presents the experimental reflection coefficient and points out the differences between the prototypes.

Three experiments are carried out to measure the reflection coefficient magnitude and its phase in independent realizations. The dynamical reflection coefficient is captured for every frequency of interest in three situations. The antenna is placed inside the anechoic chamber and works under no proximity effects in the static case. Each recording contains 5000 points because the variability should be very low. Data and talk modes are measured when the user interacts with the radiating structure during 4 minutes approximately, which turns out in 25000 samples for each frequency. The UTRA/FDD LTE B8 and UTRA/FDD LTE B3 bands [19] are tested in the up- and downlink operative frequencies: (i)UTRA/FDD LTE B8 uplink (UL) test frequencies: 880, 897.5, and 915 MHz(ii)UTRA/FDD LTE B8 downlink (DL) test frequencies: 925, 942.5, and 960 MHz(iii)UTRA/FDD LTE B3 uplink (UL) test frequencies: 1710, 1747.5, and 1785 MHz(iv)UTRA/FDD LTE B3 downlink (DL) test frequencies: 1805, 1842.5, and 1880 MHz

Finally, the samples are mixed to estimate the total in-band variability region, which is of interest, for example, to determine the required coverage area for an adaptive antenna matching network or to compute the element values in reconfigurable antenna structures. Unlike the isolated antenna, with very small variations, the matching domain representation states the impedance range as a result of the user interaction, as shown in Figure 4 in the particular case of the Descartes model.

**(a)**

**(b)**

In fact, Figure 4 points out the difference in the measured scenarios. Both UTRA/FDD LTE B8 and UTRA/FDD LTE B3 impedances are clearly observed at the six corresponding frequencies. These curves spread to the variability regions in data or talk modes, and the final dynamical impedance region is their ensemble. Furthermore, the data distributions are not uniformly spread since not all the measurements have the same probability. These PDFs depend on the frequency and the antenna geometry. According to the measurements, the variability decreases when the frequency increases due to a stronger coupling in lower bands (the user is closer in terms of wavelength in the UTRA/FDD LTE B8 standard).

Even though the impedance variability region can be estimated from the measured histogram, the impedance region is strongly dependent on the measured reflection coefficient PDF which is a disadvantage because the comparison between different antennas can be unavailable. Thus, it may be more reliable for engineers to seek a statistical alternative with few parameters which is applicable to any radiating structure.

#### 4. Theoretical Analysis Framework: Uncertainty Impedance Ellipses

The proposed general framework is based on the analysis of the uncertainty [20–22] in the reflection coefficient. The study of the impedance PDF support has defined a nonuniform area (), which is different even for the same antenna operating as a multiband structure. We propose to approximate by means of an impedance ellipse, which may be a convenient method to express the uncertainty of the antenna input impedance in a graphical format.

If the observed reflection coefficients are referred to the Cartesian coordinate system on the space, they can be written in a vector notation as

Then, we could calculate the estimated mean corresponding to samples as

The estimated covariance matrix should be computed as where and denote the standard deviations corresponding to the real part and the imaginary part of the reflection coefficient and is the correlation index for and . The quadratic form is the probability impedance ellipse and an approximation of . A certain percentage of the actual impedance values (depending on the statistical distribution) is within the ellipse. Typically, is chosen in terms of the standard deviation . If , or , these are called , and impedance ellipses. If this parameter increases, the ellipses have larger axes and more impedance values are contained inside. We will discuss a method to choose this value in the following section, but the idea is that it suggests a confidence bound which depends on the statistical distribution.

The impedance ellipse can cast as containing the information on the dynamic coupling and proximity effects caused by the user. The main advantage of the proposed method is that coupling effects can be viewed with a novel geometric perspective. Furthermore, only five geometrical parameters are required for modeling: the ellipse center, the major and minor axis, and the rotation angle (related to the space axes [23], as Figure 5 shows). Therefore, the antenna engineer, with the knowledge of these five magnitudes, is able to build an approximate variability region with a certain confidence level, which can be used for developing proper matching networks or evaluating antenna performance.

The algorithm to compute those five geometrical parameters relies on the analysis of the statistical data.

##### 4.1. Ellipse Center

Compute the ellipse center and the first and second values, by using (11).

##### 4.2. Principal Axes

Calculate the principal axes with the aid of the matrix. Notice that we can perform a single-value decomposition (SVD) on to obtain two eigenvalues, and , and two eigenvectors, and , so that where and are the diagonal and orthogonal matrices whose columns are the eigenvectors such that

The principal axes are given by where is the rotation angle of the ellipse (third parameter).

##### 4.3. : Major and Minor Axes

The fourth and fifth parameters are computed which provided a specified value of . Assuming that is the maximum eigenvalue,

To sum up, our approach evaluates the proximity user effect through the uncertainty impedance ellipse with parametric boundary where is a variable which generates the ellipse contour

#### 5. Antenna Evaluation and Results

The proposed evaluation method is applied to the UTRA/FDD LTE B8 and UTRA/FDD LTE B3 antenna prototypes in order to show how the impedance ellipse fits the variability region. The algorithm estimates the geometrical parameters in the multiband dynamic situations (using the ensemble of individual frequencies). Table 1 shows the results for those parameters.

The uncertainty impedance ellipse contains more points when the factor increases according to Figures 6 and 7. In fact, almost all the observations are included at .

The geometric values allow obtaining different comparisons and effects on the antennas as well as the approximate IVR. As far as the meaning of each parameter is concerned, we realize that the larger the impedance ellipse area, the higher the covariance matrix eigenvalues. This comes out as a result of more influence of the surrounding objects on the radiating structures. Therefore, the eigenvalue product, proportional to the ellipse surface, is a measure of the impedance spreading.

The ellipse angle is related to the correlation between the antenna resistance and the reactance because it addresses the direction of the main axis, being close to the most usual impedance values. Thus, a positive angle suggests that random changes increasing the antenna resistance tend to result in a reactance displacement towards the inductive Smith chart region (i.e., either reducing the capacitance or increasing the inductance) whereas negative values state the opposite. As the deep connection between the reactance and the near radiated field is well established, we believe the ellipse rotation is strongly influenced by the geometry because it is measuring the reactance sensitivity which provided a change in the resistance (e.g., if , only the resistance is affected by the user). However, we should also notice if the eigenvalues are very close, the ellipse is approximately a circumference and, thus, the rotation angle is meaningless.

Finally, the ellipse center is the averaged reflection coefficient which locates the variability region inside the Smith chart. This value may not be the statistical mode (the most frequent reflection coefficient value) because the impedance distribution might not be symmetrical or could be bimodal like in the case of the Descartes antenna in the UTRA/FDD LTE B8 frequency band. The ellipse center can be cast as the averaged dynamic impedance. Its nature (capacitive or inductive) takes into account the fact that the dynamic antenna impedance is more spread in the capacitive or inductive Smith chart subspace.

The proximity effects are compared in terms of three aspects according to the interpretation of the ellipse geometry parameters: (i)Influence of the frequency band: the eigenvalues are higher in the UTRA/FDD LTE B8 frequency band, and proximity effects are more severe, as expected, at lower frequencies. Since the ellipse surface is a measure of the variability, we define the following factor of merit: The evaluation results are 2.84, 1.92, and 1.11 corresponding to the Cantor, Descartes, and T-monopole structures, respectively.(ii)Antenna geometry effects: instead of taking into account the multiband aspects, the same idea may be applied to evaluate the structure influence, defining a geometric figure of merit as For example, according to the numerical results regarding the UTRA/FDD LTE B8 frequency band (, , and ), the T-monopole is more affected by the user proximity than the Cantor structure. This conclusion remains also valid in the UTRA/FDD LTE B3 frequency band.(iii)Scenario dependence: data and talk modes have strong differences in their impedance ellipses. We have found that talk mode variations are usually smaller than data mode variations. This result is reasonable due to the user mobility under a phone call regime, which is less than in the data mode (the head dynamics is not as high as around the body). The final uncertainty impedance ellipse (ensemble of the two scenarios, named “data-talk” or “D + T”) may be inside, intersecting, or completely outside the data ellipse. In the first case, the trend is that it is slightly smaller than the data ellipse. Unlike the inside case, if the data-talk ellipse is outside; its size is bigger than in data scenarios. The talk mode rotation angle is always positive, suggesting displacements towards the inductive Smith chart region when the user effects increase the antenna resistance.

#### 6. Confidence Level Analysis: 95% Uncertainty Ellipse

The uncertainty impedance ellipse is a conic curve which depends on the parameter. The proximity effects may be analyzed by studying the covariance matrix, but some final variability region is required when designers deal with matching issues. We propose to estimate the final IVR by a confidence analysis, that is, the number of observed values which is contained by .

The confidence is estimated by counting the number of observations inside the ellipse. This number depends on the reflection coefficient distribution. If the data were Gaussian, the value could be found relying on the distribution. However, the impedance PDF can be clearly far from Gaussian and a method to approximate the confidence must be developed.

The following function counts the number of observations: where the condition if is equivalent to being that and are the ellipse foci and is the major ellipse axis. Finally, the confidence is then estimated as

Table 1 shows (columns from to ) the estimated confidence for every evaluated case and the value corresponding to a 95% confidence (the column). A major confidence is achieved when , but it is not constant due to the PDF shape. If the distribution was one-dimensional, at least of the data would be accurately modeled, despite the density of the measurements. Although this result is not directly applicable to multivariate distributions, it includes a representative set of observed values [24]. In addition, more than 95% data on average is inside the uncertainty impedance as a result of the proposed counting approach.

The usual suitable trade-off sets a confidence level. This compromise can be used for estimating the value by solving . Nevertheless, the approximated parameter should be near the real number because the observed sample number *N* is very large.

Figures 8 and 9 show the T-monopole antenna case ellipse and the investigated scenarios (data + talk, data, and talk modes). The semiaxis magnitudes ( and ) presented in Table 1 are computed experimentally with an estimated confidence level.

One of the main applications of such ellipses in matching issues is the computation of real reconfigurable tuning unit coverage areas. The criteria to design adaptive impedance networks cover the whole Smith chart space according to the classical point of view. In fact, the non-idealities introduced by circuit elements (switches, capacitors, and PIN-diodes) reduce the coverage and make the development to be a difficult task. The impedance ellipse can be taken as a starting point to focus the tuning unit design onto that area, which may result in more simple matching networks.

Indeed, the coverage factor reduction can be evaluated with the aid of the proposed method. Let us think in the ideal tuning unit which covers all the Smith chart and an optimized network which fulfils the actual requirement. Then, the coverage factor () is the ratio of the areas of the total smith chart, and the uncertainty impedance ellipse

The numerical evaluation (included in Table 2) shows that only 24% of the coverage is actually required in the worst case to design the adaptive impedance matching network.

#### 7. Conclusion

User interaction effects on portable antenna input impedance can be quantified in geometrical terms using a novel framework which is applicable to any structure, coupling objects, and frequency band. This technique relies upon the mathematical support of the experimental statistical probability density functions for the reflection coefficient, which is estimated by means of stochastic magnitudes (the covariance matrix) and the so-called uncertainty impedance ellipse. Its main parameters (center, major and minor axes, and rotation angle) allow engineers to know with certain confidence level the set of feasible impedances which could have an antenna and, thus, to evaluate dynamic mismatching. Furthermore, the method has been applied to three structures predicting some trends in the antenna impedance. This algorithm is a powerful tool for making an accurate choice among possible portable antenna candidates and for reducing coverage requirements on impedance tuning networks.

#### Conflicts of Interest

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