Abstract
We discuss the mitigation of truncation errors in spherical-scanning measurements by use of a constrained least-squares estimation method. The main emphasis is the spherical harmonic representation of probe transmitting and receiving functions; however, our method is applicable to near-field measurement of electrically small antennas for which full-sphere data are either unreliable or unavailable.
1. Introduction
The transmitting function of an electrically small probe tends to be very broad, so that full-sphere data are needed to compute the spherical-harmonic representation of the pattern by use of standard methods [1]. Unfortunately, backward-hemisphere data are often unavailable or at best unreliable due to support-structure blockage. Simply setting the backward-hemisphere data to zero leads to ringing effects that degrade the accuracy of the probe-pattern representation in the forward hemisphere. In this paper, we approximate the spherical-harmonic expansion of transmit/receive functions using only forward hemisphere data. The solution presented is a least-squares fit with energy constraints (to limit backward-hemisphere radiation) and proper weighting of measurements (to compensate for clustering near the poles). The energy constraint is specified using an estimated value of directivity. Results are given that compare pattern representations obtained using both full-sphere and half-sphere data.
Others have also considered truncation error mitigation in various near-field measurement contexts [2–5].
2. Theory
We wish to approximate a physical quantity by a spherical-harmonic expansion Here, and are vector spherical harmonics [6, chapter 16]. For example, might be the measured transmitting or receiving pattern of a probe, or it might be the “measurement vector” for a spherical near-field scanning measurement [1, 7–9].
When full-sphere data are available, the orthogonality of vector spherical harmonics can be used to determine the coefficients and from discrete measurements of gathered on a uniform grid in and . When so determined, is optimal in the sense that is minimal.
In this paper, we assume that data are available only in the forward hemisphere on the grid with Requirement (5) preserves the standard number of measurement points by halving the maximum sample interval in . Ideally, and are products of small primes to allow efficient application of the fast Fourier transform (FFT) algorithm; however, factorization may not be important when antennas are “electrically small.” Our goal is to minimize the discrepancy where the are positive weights. Equation (6) may be viewed as an adaptation of (3). Because the sampling theorem has been satisfied, the integral can be readily evaluated analytically. The integral, however, has been discretized. (We generally set ).
The discrepancy is minimized subject to the “energy” constraint in order to limit the power radiated into directions where there are no measurements. When is proportional to the transmitting pattern of the test antenna, we have so that may be chosen using an estimate of the directivity in some measurement direction . (Directivity and gain are approximately equal when ohmic losses are small).
Avoiding details for the moment, we write where and are vectors of measurements and predicted values, is a vector of coefficients to be determined, and the matrix represents a known linear relationship. We choose to minimize subject to the constraint . Here, is a diagonal matrix of positive weights, and an asterisk implies Hermitian transpose. It is easy to show that this is determined by the equations The significance of the Lagrange multiplier is seen in the relation that is, is a measure of the residual of the fit. In the absence of the energy constraint (7) of course, we can choose to be the solution of the normal equations .
Since is Hermitian nonnegative definite, we have the singular-value decomposition (SVD) where is unitary, , and . From (13), (10), and (11), we are led to The right side of (14) tends to 0 as , and because of poles at the singular values , there may be numerous solutions for . According to (12), the best choice corresponds to the unique solution for which , where is the smallest singular value. In general, , since the constrained solution is expected to have lower energy than the unconstrained solution. Finally, given , we may compute , Numerical methods used in this paper are discussed in greater detail in [10].
At this point, we have discussed the main ideas of this paper. What follows are some rather unpleasant details that we summarize dutifully and concisely. To begin with, (1) can be written as a Fourier series where the coefficients may be computed from the data with a discrete Fourier transform
Similarly, for (2),
With substitution of (16) and (18), the discrepancy (6) becomes This may be rewritten in the form Here Finally, with the definitions we have Thus, the measurement truncation problem has been reduced, explicitly, to the constrained least-square optimization problem discussed earlier. In particular, (10) and (11) can be written as Because of the block-diagonal structure, the overall SVD can be broken into smaller parts, one at a time, resulting in an algorithm with computational complexity . Comparing unfavorably with the complexity for the full-sphere case, our method may be less useful for larger antennas. The order of is , so there will be more equations than unknowns in every block as long as . Actually, there is always a unique solution to our constrained optimization problem; however, the quality of this solution is observed to degrade dramatically as is decreased below the threshold .
3. Experimental and Simulated Results
To investigate the utility of our least-squares technique for reducing the truncation error that results from zero filling in the rear hemisphere, we examined a number of typical probes for which full-sphere far-field patterns are available. In each case, we computed the far-field pattern using the spherical mode expansion obtained (a) from the standard algorithm with full-sphere data, (b) from the standard algorithm with zero-fill in the backward-hemisphere, and finally, (c) from the constrained least-squares fit to forward-hemisphere measurements described previously. Ideally, the far-field pattern calculated from the spherical-mode expansion should agree with the original pattern. For the full-sphere modal calculation, this is of course the case if the sampling theorem is satisfied. We show E-plane results only, since truncation effects are more important because the patterns tend to be broader. For all cases, the sampling is 2° in and 5° in .
3.1. High-Gain Symmetric Probe
In this case, the input data were simulated using a far-field pattern calculated from a specified set of () modal coefficients. For this antenna, we used , . The directivity of dB was used in (8) to determine the energy constraint. Figures 1 and 2 indicate excellent agreement between the full-sphere and least-squares techniques in the forward hemisphere.
In this example, the constrained least-squares technique also performs well in the backward-hemisphere if the correct energy is specified. This is not necessarily true in the presence of noise. As a test, we added random errors with an RMS value of 0.1% relative to the noise-free RMS pattern level. An unconstrained fit resulted in an on-axis directivity of dB. In other words, the back hemisphere overwhelmed the forward hemisphere. On the other hand, a constrained fit forced the back hemisphere pattern to remain at a reasonable level even though the details were incorrect. Although an unconstrained fit may be somewhat better in the forward directions, a constrained fit is undoubtedly more reasonable given our knowledge of the directivity of the test antenna.
In probe-corrected spherical near-field measurements, for example, the complete probe pattern is required. We often argue that the forward hemisphere pattern is most important, but all bets are off when the back hemisphere pattern dwarfs the forward hemisphere pattern.
3.2. NIST Circularly Cylindrical Waveguide Probe
This cylindrical waveguide probe was designed and built for spherical near-field measurements at 3.3 GHz. The far-field pattern was measured over the entire sphere, although support-structure blockage limits the value of backward-hemisphere information. In this case, we used , . A directivity of dB was calculated using the full-sphere data. Deviations (Figure 3) between the full-sphere and zero-fill methods will lead to significant errors when the zero-fill results are used in spherical scanning measurements. Within the range , the maximum error in the zero-fill result is about dB, while the least-squares technique has a maximum error of dB. Figure 4 shows how the full-sphere, zero-fill, and least-squares results compare in the backward-hemisphere. In this case, the truncation level is as high as dB, which causes a greater discrepancy between the zero-fill and full-sphere patterns than observed in the previous example.
3.3. Rectangular Waveguide Probe
This probe is a section of WR-284 rectangular waveguide. Far-field patterns were obtained over the entire sphere at 3.3 GHz. In this case, we used , . Strictly speaking, rectangular waveguide probes do not have the correct symmetry for spherical near-field scanning applications; however, they often perform satisfactorily when the scan radius is more than about an antenna diameter. Thus, we include results for this type of probe. A directivity of dB, calculated from the full-sphere modal expansion, was used to determine the energy constraint. Measured gain was dB. The results are similar to those for the cylindrical waveguide case. As evident in Figure 5, within the angular range , the zero-fill error is greatest at about 40°, where it is dB. The error is only about dB for the least-squares technique at this angle. In the forward hemisphere, the zero-fill technique exhibits large errors, while the least-squares technique errors remain relatively small. Figure 6 provides a full-sphere comparison of the three techniques.
3.4. Choice of
Several tests indicate that the quality of the least-squares fit in the forward hemisphere is relatively insensitive to the choice of . In a typical case with , the errors observed in the results were generally much less than dB within 20° of the on-axis direction. (Here, is the correct value). When gain is used to approximate directivity, experience indicates that the error in usually will be less than .
4. Summary and Future Work
We have demonstrated a constrained least-squares technique for calculating the spherical mode coefficients of a small antenna from forward-hemisphere far-field data. This technique can significantly reduce truncation errors that arise when the standard algorithm is used with zero-fill in the backward-hemisphere. While the least-squares algorithm is not especially efficient, computational times are still acceptable for small antennas.
Our technique can also be used to process near-field data when reliable full-sphere measurements are not available. The method may be adapted to serve when measurements are made in the range , where is not restricted to . We are considering an iterative algorithm to find both and simultaneously. For well-conditioned systems, this could effectively reduce the computational complexity to . By beginning with a spherical-wave expansion and projecting onto a scanning surface, our method should be useful for mitigation of truncation error in planar and cylindrical scanning measurements, as well.
Acknowledgments
The authors thank Allen Newell, consultant, for bringing this problem to our attention. They also thank Brad Alpert of NIST for insightful discussions concerning the numerical methods employed in this work.