COMPUTATION OF ANTENNA SIDE-LOBE COUPLING IN THE NEAR FIELD USING APPROXIMATE FAR-FIELD DATA
Michael H. Francis Arthur D. Yaghjian
National Bureau of Standards U.S. Department of Commerce Boulder, CO 80303
August 1982
"*QC™ 100
, U56
I; 82-1674 1982 C.2
NATIONAL BUREAU OF STANDARDS! LIBRARY
NQV 2 1932 NBSSR 82-1674 nc*< GlQ, \0O COMPUTATION OF ANTENNA SIDE-LOBE COUPLING IN THE NEAR FIELD USING APPROXIMATE FAR-FIELD DATA
Michael H. Francis Arthur D. Yaghjian
Electromagnetic Fields Division Center for Electronics and Electrical Engineering National Engineering Laboratory National Bureau of Standards U.S. Department of Commerce Boulder, CO 80303
August 1982
U.S. DEPARTMENT OF COMMERCE, Malcolm Baldrige, Secretary
NATIONAL BUREAU OF STANDARDS, Ernest Ambler, Director I
I CONTENTS Page Glossary of Symbols for the Text..... iv
Abstract. 1
1. Introduction. 2
2. Side-Lobe Coupling of Two Antennas Using Their Exact and Approximate Far Fields...... 4 2.1 Statement and Formulation of the Problem.... 4 2.2 The Hypothetical Circular Antenna. 5 2.2.1 The Exact Far Field..... 5 2.2.2 The Approximate Far Field for Approximation-1...... 6 2.2.3 The Approximate Far Field for Approximation-2...... 6 2.2.4 Limitations on the Approximations...... 7 2.3 Results and Comparisons... 8 2.4 Conclusions..... 10
3. Mathematical Upper Bounds Derived from Equation (1 11 3.1 Assumptions and Integration. 11 3.2 Comparison of Upper Bound from Formulas to the Exact Upper Bound...... 15 3.3 Conclusion...... 15
4. Summary and Concluding Remarks 16
5 . References 18
6. List of Figures 19
Appendix A. Derivation of Ay and A^y in Terms of MA y A Antenna Gains and Rel ati ve Side-Lobe Levels...... 40
Appendix B. Documentation of ENVLP, the Computer Program to Estimate Coupling Between Two Antennas... 42 B.l General Overview of Computer Program ENVLP. 42 B.2 Computer Code and Sample Output...... 48
Appendix C. Subroutine FOURT...... 70
Appendix D. Subroutine PLT120R..... 78 : .
Glossary of Symbols for The Text
a radius of the antenna aperture.
a ,a : radius of the transmitting and receiving antennas, respectively, T R a°: incident wave-amplitude in the transmitting antenna feed,
a : the of the far field in direction MAX maximum amplitude the of the axis of separation.
A A : the maximum amplitude of the far field in the direction of TMAX’ RMAX the axis of separation for the transmitting and receiving antennas, respectively. the emergent wave-amplitude in the receiving antenna feed.
1 C' * n Z (1 - r r ) oo o L d: the separation distance between antennas.
D ,D : the diameter of the smallest sphere which circumscri bes the t r transmitting and receiving antennas, respectively,
E_: the electric field. the far field of the receiving and transmitting antennas,
respectively, with components f f', f' and f f f . x , y z x , y , z G: the antenna gain.
G ,G : the gain of the transmitting and receiving antennas, t r respecti vely V the Bessel function of first order. k_: the propagation vector with components k k y. x , y ,
k: /k « k = 2tt/A .
k : D )/2d. o MD t + R
K_: k e + k e . xx y y K: /irnr .
the input power to the antenna, ^i nput L position vector with components x, y, z. r magnitude of _r« A A R xe + Jye . x y S the relative side-lobe level in the direction of the separation axis. S ,S the value of S for the transmitting and receiving antennas, T R respectively. Z the wave impedance of free space, o Y the z-component of k_. r the reflection coefficient of the receiving load. L r the reflection coefficient of the receiving antenna. o the characteristic for the propagated mode in the n 0 admittance waveguide feed of the receiving antenna. A the wavelength. the Euleri an angles of the transmitting antenna, e antenna. ^R’ R’^R the Euleri an angles of the receiving COMPUTATION OF ANTENNA SIDE-LOBE COUPLING IN
THE NEAR FIELD USING APPROXIMATE FAR-FIELD DATA
M. H. Francis and A. D. Yaghjian
National Bureau of Standards Boulder, Colorado 80303
Computer programs, in particular CUPLNF and CUPLZ, are presently in existence to calculate the coupling loss between two antennas provided that the amplitude and phase of the far field are available. However, for many antennas the complex far field is not known accurately. In such cases it is nevertheless possible to specify approximate far fields from a knowledge of the side-lobe level of each antenna along the axis of separation, and the electrical size of each antenna. To determine the effectiveness of using approximate side-lobe level data instead of the detailed far fields, we chose as our test antennas two hypothetical, linearly polarized, uniformly illuminated circular antennas for which the exact far fields are given by a simple analytic expression. The exact far fields are supplied to the program CUPLNF to compute the exact near-field coupling loss. Approximate fields are supplied to a new program ENVLP developed for the purpose of computing the approximate near-field coupling loss. The comparison of the results from ENVLP to those of CUPLNF indicates that the use of approximate far fields gives an estimate of the coupling loss which is good to about +5 dB. In addition, the plane-wave transmission formula for coupling between two antennas is used to estimate upper-bound values of coupling loss. These upper bounds are compared with the maximum coupling losses obtained from programs CUPLNF and ENVLP.
Key words: antenna coupling; antenna theory; coupling loss; near-field measurements.
1 1. Introduction
Three years ago, the Antenna Systems Metrology Group, Electromagnetic Fields
Division of the National Bureau of Standards, under the sponsorship of the
Electromagnetic Compatibility Analysis Center, completed the development of a highly efficient computer program, CUPINF, for calculating the coupling loss between two antennas given the far fields of each antenna [1,2]. For antennas arbitrarily oriented and separated in free space, CUPLNF computes the coupling loss at a single frequency as a function of antenna displacement in a plane transverse to the axis of separation of the antennas. Multiple reflections between antennas are neglected but no other restrictive assumptions are involved. CUPLNF automatically computes electric fields from the transmitting antenna when a "virtual antenna" of uniform far field replaces the receiving antenna. Coupling loss for antennas hundreds of wavelengths in diameter is computed in a few minutes of CPU time within the central memory core,
1 e.g. , of a CDC 6600.
The major limitation of CUPLNF is its requirement for the magnitude and phase of the electric far-field components of each antenna within the solid angle mutually subtended by the antennas. For example, if two antennas are oriented so that coupling occurs mainly through their side-lobe region, the complex vector far field of each antenna covering this angular side-lobe region must be supplied to the program CUPLNF.
In practice, one may not have such detailed information on the side-lobe far fields to estimate the coupling loss between two antennas co-sited in their near field. Often, one has some knowledge of the side-lobe levels of one's antennas, even if detailed phase and amplitude information of the field components is not avail- able. Thus a natural and important question to answer is whether this limited information, specifically side-lobe level of each antenna near their axis of separation, can be used to estimate the antenna coupling in the near field using a new program, ENVLP. Of course, as the separation distance between the antennas approaches the mutual Rayleigh distance, i.e., the antennas lie in each others far field, the ordinary far-field formula for coupling between two antennas can be used to estimate coupling loss. For separation distances much less than the mutual Rayleigh distance, the far-field formula would not be expected to give realistic estimates of coupling.
It is in this near-field region that we set out to decide if ENVLP will give reason- able values of side-lobe coupling (±6 dB accuracy) by utilizing only side-lobe levels
i The specific computer is identified in this paper to adequately describe the computer program. Such identification does not imply recommendation or en- dorsement by the National Bureau of Standards nor does it imply that the computer identified is necessarily the best available for the purpose.
2. .
along the axis of separation of the antennas (in addition to their size, separation distance, and feed characteristics) To determine the efficacy of ENVLP using approximate side-lobe level information instead of the detailed far fields, we chose as our test antennas two hypothetical, linearly polarized, uniformly illuminated circular aperture antennas for which the exact far fields are given by a fairly simple analytic expression involving functions no more complicated than the first order Bessel function* With these hypothetical antennas, the exact far fields could be supplied to the program CUPLNF to compute near-field coupling loss without introducing the errors incumbent with measured far- field data. An approximate side-lobe far field can be substituted for the exact far field and a comparison of the near-field coupling loss computed with the exact and approximate far fields can be made simply and straightforwardly.
At first thought it may seem overly optimistic to want to obtain reasonable estimates of near-field coupling having only side-lobe level information, i.e., having only the envelope of the magnitude of the far fields in the direction of separation.
Fortunately, however, we can also estimate the phase variation of the side lobes of most microwave antennas merely from the frequency of operation and the radii of the antennas, more precisely, from the radius in wavelengths of the smallest sphere which circumscribes all significantly radiating parts of each antenna. Also, even though polarization characteristics of the side-lobe far fields may be unknown, we can assume polarization match of the far-fields of the antennas over the relevant range of integration. This assumption of polarization match tends to produce an upper-bound estimate of coupling loss, which may be of appreciable value to the user. We also assume the user has a knowledge of (1) the separation distance between the antennas, and (2) the reflection coefficients in the waveguides feeding the antennas--! n all a rather minimal amount of information about one's antennas.
In Appendix B, we document program ENVLP, which is a modification of the program
CUPLNF requiring as input this minimal amount of input data, i.e., the antennas 1 radii, separation distance, reflection coefficients, and relative side-lobe levels in the direction of separation.
Finally, we work directly with the plane-wave transmission formula for coupling between two antennas to estimate upper-bound values of coupling loss, comparing these upper bounds with the maximum coupling losses obtained from program CUPLNF and program ENVLP using exact and approximate far-field input. We anticipate these simple upper- bound formulas being especially valuable when the major requirement for co-sited antennas is that the coupling interference or fields lie below a certain threshold.
3 —
2. Side-Lobe Coupling of Two Antennas Using Their Exact and Approximate Far Fields
2.1 Statement and Formulation of the Problem
Yaghjian [1] showed that the coupling between two antennas in free space, neglecting multiple reflections, is given by
iyd i K *R
b' - f(k)e e o
II dK , ( 1 ) a o K = = =~ where C -, and 1/(1 - r r. )is the mismatch factor of the receiving n r antenna, is the characteristic admittance of the propagated mode in the n 0 waveguide feed of the receiving antenna, Z is the wave impedance of free space, k is the Q magnitude of the propagation vector, K_ = k e + k e is the transverse part of the x x y y propagation vector, = (k 2 - K 2 1 ^ 2 dK_ is an abbreviation for dk dk and the y ) , x y , coordinates (R_, z = d) give the position of the origin O' fixed in the receiving antenna with respect to the (x,y,z) coordinate system fixed at 0 in the transmitting antenna (see fig. 1). The functions _f and _f‘ are the electric far fields of the transmitting and receiving antennas, respectively, without the presence of the other, r and are the reflection coefficients of the receiving antenna and receiving load Q respectively, a is the amplitude of the input to the transmitting antenna, and is Q the amplitude of the output of the receiving antenna. Formula (1) assumes the receiving antenna is reciprocal, although the formula easily generalizes to include non-reci procal antennas. Yaghjian [1] also showed that for most practical purposes the integration range in (1) could be limited to K/k < (Dy + D^)/2d for R = 0 provided ( D-j- + D^)/2 < d < (Dj + Dr)^/A. Dy is the diameter of the smallest sphere which circumscribes the transmitting antenna and D is the diameter of the smallest sphere R ci rcumscribi ng the receiving antenna. (If Dy or D^ is less than twice the wavelength, y, Dy, or is replaced by 2A.) For coupling in the transverse plane over the range - (Dy + D^) < R < Dy + D the integration range should be doubled [1]. r 1 There are some situations where detailed information of the far fields (_f , f) is lacking and it is desirable to get an estimate of the coupling between two antennas. For example, both amplitude and phase information may not be available. We would like, therefore, to consider some possible methods of approximating the coupling and in particular of finding a good estimate for the maximum coupling between two antennas in the general direction of the separation axis (the separation axis is drawn from a 4 point centrally located on the transmitting antenna to a point centrally located on the receiving antenna for the antennas in their initial position, R_= 0). To find the maximum coupling in the general direction of the separation axis we shall calculate the coupling over a range of values of R_ (R_ is perpendicul ar to the separation axis). Specifically, we shall move the receiving antenna in the x and y directions, over a range from - (Dy + D ) to + (Dy + D ) for both the x and direc- R R y tions. In order to save computer time and cost we will hold y = 0 while varying x (the XO cut) and hold x = 0 while varying y (the YO cut). If we are to determine how accurate our approximate results for coupling loss are, we need to compare them to the results that we would have obtained if we knew and used the exact far fields of the antennas for finding the coupling quotients. For this purpose we will use two hypothetical circular antennas and compare results obtained using the exact far fields to those obtained using approximate far fields. These results will be compared for different values of the separation distance, diameters of the antennas, frequency, and orientation of the antennas. 2.2 The Hypothetical Circular Antenna 2.2.1 The Exact Far Field We consider the case of a uni formally illuminated circular aperture with the transverse electric field polarized in the x-di recti on. Then the far field for a point in the (9, Ji(ka sine ) f — cose (2a) x "ka sine /tt (2b) J ( ka sine) ka ! f s1ne “* " • (2c) z Tea We- /tt For this formula, the circular aperture lies in the xy plane with the center of the aperture at the origin. The symbols, e,<|> denote the usual angles in spherical coordinates, k is the magnitude of the propagation vector, a is the radius of the aperture, and J is the Bessel function of first order. An example of this pattern L can be seen in figure 2, where ka = 20.94. D ) 2.2.2 The Approximate Far Held for Approximation-1 Approximation-1 is specifically geared to the uniformly illuminated circular aperture antenna. For this approximation we replace cos(ka(e - e )) J , (ka sine) by ka sine Hka sine ka sine o o where e is the direction of the separation axis relative to the preferred antenna Q coordinate system. 0 equals for the transmitting antenna and e for the receiving Q 0^ R antenna (see the sketches on figures 6 through 20). This electric field pattern can be found for e = 60° in figure 3. The values of the field are given relative to the main o beam at e 0° in figure i.e. the field is normalized to a main beam equal to 1. o = 2, The above approximation may be justified as follows. J^ka sine)/ka sine may be replaced using the asymptotic expansion for J . This yields: J ( ka sine) -j /ZlZ cos(ka sine - 3 tt/ 4 ka sine /irka sine ka sine The cosine term represents the' variation of the sidelobe while the rest of the term is the amplitude of the envelope. We replace the terms involving the amplitude of the envelope by their values at 0 and expand the cosine term about e . Thus, O Q cos(ka sine - 3u/4) ~ cos (ka sine_ - 3 tt/ 4 + ka(0 - eHcoseH. The ka sine,, can be lumped with -3 tt/ 4 into a phase factor which we will ignore since we are finding only an average coupling over a solid angle. The cos0 term has been replaced by 1 because Q actual antennas do not display this cose dependence of sidelobe width peculiar to the o hypothetical finite planar aperture distribution. 2.2.3 The Approximate Far Field for Approximation-2 The approximation-2 far field is a general approximation that can be made for a 1 large class of coupling cases. For this approximation we replace ( -k ) *f JO by JF ( A cosk a t o sk a A cosk a cosk a TMAX x T y T RMAX x R y R . ' A anc* are rnax1 mum m n i Guides of the far fields in the general direction TMAX Armax a§ of the separation axis (i.e., the side-lobe level) for the transmitting and receiving antennas, respectively, k = ksin0 cos 6 . separation axis and are given in absolute SI (mksA) units. In Appendix A we derive in terms the side-lobe levels and the gains of the aTMAX and a RMAX transmitting and receiving antennas, ay and a^ are the radii of the smallest spheres circumscribing the transmitting antenna and the receiving antenna, respecti vely /v- We have chosen the above approximation because for many antennas the components of the far-field patterns vary roughly as cosk a where a is the radius of the x coskxa, smallest sphere circumscribing all the significantly radiating parts of the antenna.^ Note that for approximation-2 we assume that the polarization of the receiving and transmitting antennas are matched and thus approximation-2 tends to yield an upper bound. The approximation-2 pattern for = 0° or 90° and the separ- ation axis in the direction of 9 = 60° can be found in figure 4. The values of the field are given relative to the main beam at e = 0° in figure 2. 2.2.4 Limitations on the Approximations Both approximation-1 and approximation-2 are valid over a relatively narrow range of angles. In particular, the approximations should be valid (under the above stated limitations) if the amplitude of the envelope of the far field of each antenna varies by not more than about 3 dB over the range of angles mutually subtended by the transmitting and receiving antennas. Thus, for any case in which the integration of equation (1) must be performed over a large range of angles, approximation-1 and approximation-2 may be poor. This will, in general, include those cases involving coupling with the main beam. It will also include electrically small broad beam antennas and the unusual cases where both the receiving and transmitting antennas are identical and also have identical Eulerian angles [1] (see figure 5 for a definition of these angles). In addition, for 9 = 0 the approximation of the far field for q p ‘•The far side-lobe region of the far-field pattern is predominantly caused by diffraction from edge points of the antenna. These edge points are usually separated by a distance of approximately an antenna diameter and this leads to a side-lobe pattern which has a null approximately every x/2a radians, as numerous hypothetical and measured far-field patterns confirm (see for example Johnson et al . [4], and Newell and Crawford [5]). We assume a smooth varia- tion in the form of a cosine, and that the k and k variation is separable. This leads to the cosk a cosk a dependence aoout the axis of separation - a x y fairly reasonable approximation provided the antenna is not highly elongated (i.e., the length in one direction is not much greater than the length in the other direction). In the case of a highly elongated antenna the value of a is approximately half the longer side so that we obtain a variation which is good in the direction of the long side but poor in the direction of the short side. For such elongated antennas, a better approximation can be obtained by using a variation of the form cosk a cosk a where a is half the length of the long x x y x side and a is half the length of the short side, assuming the x and y axes are aligned with the long and short sides of the antenna, respectively. 7 approximation-1 goes to infinity. Thus, approximation-1 is never good for 2 Dr + D (D + D ) r the main beam. Finally, d must be in the range ^ < d < ^ . 2.3 Results and Comparisons The results using the exact far fields and the approximation-1 and approximation-2 far fields in equation (1) were obtained using computer programs based on CUPLNF [2]. The program for approximation-2 (titled ENVLP) was written for general use and can be found documented in Appendix B. Results were obtained for various values of antenna orientation, antenna diameter, separation distance and frequency, and are summarized in table 2.1 and in figures 6 through 20. Column 1 of table 2.1 is an identifier which corresponds to the identifier in the caption of figures 6 through 20. Columns 2 to 4 give the Eulerian angles of the transmitting antenna; columns 5 through 7 give the Eulerian angles of the receiving antenna; columns 8 and 9 give the diameters of the transmitting and receiving anten- nas, respectively, in units of wavelength. Columns 10 and 11 give the side-lobe levels relative to the main beam and within the integration range (as specified above) for the transmitting and receiving antennas, respectively. Column 12 gives the separation distance in terms of the mutual Rayleigh distance; column 13 gives the frequency in hertz. Columns 14 through 16 give the maximum coupling amplitude as given by the use of the exact far field, the approximation-1 far field, and the approximation-2 far field, respectively. The maximum coupling magnitude implied by the far-field formula is given in column 17. The far-field formula which is used to determine the maximum coupling is b' A A 0 TMAX RMAX X a n z (1 - r r ) d 0 oo o L Column 18 gives the maximum coupling amplitude as implied by the upper-bound equations (9), which will be discussed in section 3. Finally, columns 19, 20, and 21 contain the RMS mean coupling for the exact, approximation-1, and approximation-2 far fields, respectively. Figures 6 through 20 are plots of the magnitude of the coupling quotients in the x-di recti on at y = 0 (the X0 cut) for the exact far fields (_), approximation-1 (•...), and approximation-2 ( ). The Y0 cut is not shown because it is similar to the X0 cut. Each figure corresponds to one of the cases in table 2.1. 8 4444 4 1144 —4441 44 4 4 — H H H H 8 3 1 OJ hi »— CO o LO CO CD LO p^. oo O OJ LO o t— OJ S.B.5 CO d d i—4 i—4 £B CL 4-> e r-. o- rd k p^ p*. 8 8 » i i—4 i l 1 i 4 i 'f 4 1 1 «M $ 1 1* t—4 r-H i—4 - | rH OJ LO OJ LO OJ CD CO o LO OJ CD OJ LO i-H p^. d 1—4 r~- o- fS P*j» r^- ffl <\l o- i—4 III I i i 1 ? 4 4 i t CL 4-> 1 §5-1 s O LO LO P^ o o. OJ OJ CO LO LO t—4 co OJ r— CO i-H P^ CD i—4 d p*. § “o “o “S’ “S “Id “o “u “u “S’ “o “u “3> “qj ’£ 2 2 « o i—4 LO CD CO LO LO CO o 00 LO CD LO OJ C'- LO CO t—4 d i—4 SR a 8 5 t—4 O- r^. r—4 i I 4 4 i t 1 i 4 t 1 t 1 CL 1 OJ -69.0 -66.3 -71.2 -66.0 -63.8 -60.3 -64.4 -56.7 -78.6 129.3 -91.8 -50.3 -75.4 -80.6 -41.9 LO e i— E qlV> 15-i 1 1 «—4 ’x s LO o LO CD i—1 CO LO CO LO CO i—4 LO co LO e.i e LO rH r—4 CO CO LO CO LO t—4 CO CO iH CL4-> 5? r*» p*. O- LO o. r—4 r^- 8B 'T i i 1 I t i—4 8 i i 5-i 1 -»-> LO r-. CD 00 LO OJ 00 i—4 OJ 00 o LO i—4 "3" e LO CD d i— LO LO t—4 iH 3 >? p^ p*. LO i—4 s i t—4 1 UJ 4 'f i 4 4 i 1 3 o o o o 3 o o O o o 3 3 O O O d d d d o d d d d d N rH r— i-h rH i—4 t—4 i— r—4 2 2 r—4 r—4 3 co 1 S r— X X X X X X X X X X X X X X X 4O- o o o o o o o o o o o o o o J: o I— r— r-H r— i—4 i—4 LO OJ i—4 i-H i-H OJ 'te "re -r- ^ O'* cr* O'* i—4 LO OJ t—4 i— i— SB r— rH 8 8 8 8 8 o a s 8 s S >> > d d d d d d d d d d d d d d d E3 (D 00 00 00 oo 5 t—4 C£ 31.4 31.4 31.4 39.5 43.1 43.6 37.8 60.6 42.4 31.8 14.3 * 4 to « a| a s 88 8 > c ii o LO r^» 00 LO 00 CO CO LO 00 CO co CO i— u*“ r-H p*^ CO p^- i—4 CO 1—4 CO OJ CO 8 S3 » co a s LO CO 88 CO w t—4 r- p^. ^C to LO Co 8 8 to 8 8 8 8 8 Co to LO 8 LO LO LO CO CO CO CO CO LO LO LO 1—1 a a i—4 a r—4 P- oh— Co to to to LO to 8 8 8 Co to Co to to 8 LO LO LO LO LO LO co CO CO LO LO LO LO i—4 r-H t—4 a 8 8 8 8 8 8 8 8 8 8 8 s 8 8 8 1 1 T— r— i—4 i—4 i—4 t—4 i—4 i—4 r—4 i H i 4 i—4 i—4 1 oc 1 LO CD 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 LO O O O O o O O O O O O O O O O I c O O O O o O O o O O O O O O o — 1 CD1 09 60 09 CO 1 8 R 8 8 8 8 8 8 8 8 8 8 OJ — o o o O o O o O O O O O O O O i—4 OJ CO i— t— p^ i—4 OJ i—4 OJ 1 T— OJ CO 1 1 s 9 An inspection of table 2.1 will show that the maximum coupling predicted by approximation-1 compares favorably with the result given by the exact far fields. With the exception of case Omega, the difference in the maximum coupling loss between the exact far-field result and the result from approximation-1 is less than 5 dB while for case Omega, the difference is about 10 dB. The explanation for the large differ- ence for case Omega is that the wavelength for this case is so large that the electric field pattern lies near the first null (see, e.g., figure 2); thus, approximation-1 which is based on the envelope of the Bessel function is very poor in this case. We find that with the exception of two cases (Beta-3 and Delta-1) the maximum coupling loss predicted by approximation-2 is within 10 dB of that given by the use of the exact far fields. For Beta-3, Dy = 0.2m, D^ = 0.8m, and d = 1.0 m; thus, for this case the coupling integration covers a broad angular region and the polarization does not match at the wider angles as it does at the center. This means that approxima- tion-2, which assumes perfect polarization match over the entire range of integration, will appreciably overestimate the coupling for this case. In the case of Delta-1, even though the approximation-2 estimate of the maximum coupling is about 13 d3 low, this is a relatively good estimate of coupling since the exact coupling lies below -116 dB. If we now examine figures 6 through 20 and columns 19 through 21, we find that, in general, the RMS mean level of coupling generally agrees to within about ±5 dB of the exact for approximation-2 and considerably closer for approximation-1. On the other hand, approximate and exact coupling loss at individual points can differ by a substantially greater amount (occasionally more than 20 dB). It is of interest to compare our results to those obtained from the far-field formula, equation (3). It can be seen that, in general, if the separation distance between the antennas is greater than about one-quarter of a mutual Rayleigh distance, 2 (D-j. + D ) /A, the far-field formula gives results closer to the exact results than do R approximation-1 and approximation-2. On the other hand, if the separation distance is less than one-tenth of a mutual Rayleigh distance, the results from approximation-1 and approximation-2 are in general substantially better than those given by the far- field formula. For separation distances between one-tenth and one-quarter of a mutual Rayleigh distance both the far-field formula and approximation-1 or approximation-2 yield values of coupling loss of about the same degree of accuracy. 2.4 Conclusions We conclude that approximation-1 and approximation-2 can be used in the computer programs to obtain a reasonable estimate of the maximum coupling in the general direction of the separation axis. It should be emphasized that while we can use 10 approximation-1 and approximation-2 to find the maximum coupling or the mean coupling over a narrow range of directions in the general direction of the separation axis we do not necessarily get a good estimate in the exact direction of the separation axis. Approximation-1 is limited to uniformly illuminated circular aperture antennas while approximation-2 can be used for general antennas. In using either approximation-1 or approximation-2 care must be used so as to stay within the limitations stated in section 2.2.4. In the event that the separation distance is greater than about one-quarter of a Rayleigh distance the far-field formula (equation (3)) gives better results than either approximation-1 or approximation-2 and the far-field formula should be used in those cases to estimate the maximum coupling. The real advantage of the computer program for approximation-2 (Program ENVLP) is that it estimates the coupling loss between arbitrary antennas arbitrarily oriented in the near field of each other from a mere knowledge of (1) the separation distance between antennas, (2) the side-lobe level of each antenna in the direction of the axis of separation, and (3) the radius of each antenna. 3. Mathematical Upper Bounds Derived from Equation (1) 3.1 Assumptions and Integration It is possible to derive an upper bound to the magnitude of b' /a if we make some Q Q simplifying assumptions that will allow us to integrate equation (1); initially, let us make the following extremely crude assumptions which will allow us to immediately derive an upper bound for coupling quotient by performing the integration in equation (1). Afterwards a smaller, more realistic upper bound will be derived. -'(“-* * A A (4a) -- TMAX RMAX (d + d )2\1/2 t r Y Y (4b) min M ^ 4d 2 iyd i K*R , ? = e = 1 . (4c) ’ a and tlie _f f^_, respectively, near the axis of TMAX aRMAX are max1 muP1 values of and separation. With the assumptions of equations (4) we find after integrating equation -k(D + D ) +k(D + D ) (1) from vyr to «-r—— that 11 2 b* a a ( D ) 0 tmax rmax t+Dr < ( 5 ) - )2 2 a n Z (1 r r ) (d + d d 0 oo o L t r 1 - By using this integration range we limit the validity of equation (5) to 2 (D D )/2 < d < (D D ) /A. t + r T + R Equation (5) gives an absolute upper bound to the magnitude of the coupling quotient, neglecting multiple reflections, provided that coupling does not occur 2 through the main beam of either antenna, that (Dy + D )/2 < d < (Dy + D ) /A and that R R the integration range we have used [i.e., -k ( Dy + D )/2d to +k(Dy + D /2d)] is R R adequate (as explained in section 2.2.4). As d approaches the mutual Rayleigh distance we expect equation (5) to approach the far-field formula, equation (3). We can see that this is indeed the case if we (d D ry/2 t R - 2 al low \ 1 to approach 1 (note that (Dy + D ) will be much less R than 4d2 at the mutual Rayleigh distance) and replace one of the d's in d 2 2 by (Dy + D ) / A . We thus obtain R b' A A 0 TMAX RMAX A < a n Z (1 - r r ) d 0 oo o L which is just the far-field formula. If we compare equation (5) to the far-field formula for separation distances less than a mutual Rayleigh distance, equation (5) always gives a larger value. Instead of the assumptions of equations (4), we can derive more realistic upper bounds by making the more realistic approximations Dy Dy D R • f'(k) f(k) = Ay cos (k + ( + < ) A cos(k + °R X + ) COSUy-yy (f) 4 ( 0 + D )^/2 T R = = - !6b) Y Y . Ml 2 mi n 4d / iYd ia = e e , (6c) 12 k 2 where a is a constant and the lycl 1 — (6c), is made because e varies more slowly than the cosine terms or the e *— term for the integration range being considered. Under assumptions (6), equation (1) becomes 1 a k k A A e o o TMAX RMAX cos(k + *i>cos(k- + - 2 n Z (1 r r )k /I - (dt ro / xT oo o L I T R ml -k -k V 4 Df Dn + + X cos(k^-^-+ ) cos ( — i a b' A A G o TMAX RMAX Z r r )k 2 (D_ + D_) Z n (1 \V?' ~ oo oL / T R 1 " 4d^~ ( 8 ) We wish to find the maximum value of equation (8). First, we notice that ~ ^ each term in the two pairs of brackets, is of the form e”^i ^(sinBK )/B. [ ], Q We might be tempted to say that the maximum value of a term of this form is 13 ) , (1/B) but this would be true only if Bk > 1. For Bk. < 1 this term has a maximum value of at B = 0. We note that there are four values of x and LU four values of y for which B+0. They are + 0 ) (d - D ) {d - (0 (°T R t r t T = - - + V + V x, y s and In general 2 2 2 2 only one term of the x bracket and one term of the y bracket at a time have Bk < 1. However, if either Dy or D or both are sufficiently small R o | more than one term in the x bracket and more than one term in the y bracket can have Bk. < 1. The maximum value of equation (8) as a function of phase will occur when each of the phases is an integral number of 2-rr radians. Keeping this in mind and substituting for k we have the following five Q , expressions for the maximum coupling: " b' a a k(D + 0 < tmax rmax i v 1 1 1 + + (9a) a (d + h 2d 6 6 - - - 2 0 n Z (1 r r, )( 1 t i 4k 7 T ivy oo o L ^ /i a 2 ) 4d for DU 2 D 2 > — DU 2 > — T °R i k R k lb* a a 2k(D + D 0 < tmax rmax T r’ 1 i 1 [ . + _ _ T (9b) 3 Z ( 1 - r r 2d D - Max(D ,D ) i n l o oo o L | r °R T R 4d 2 - 2 2 2 2 for |D D|| max (D ,D )> min (D ,D ) ~ k ^ b' a a 2k(D + 2 0 tmax rmax T V 1 1 1 < + -f (9c) - r 2 a n Z (1 rV? (D + D ) \V2 2d 0 T d °T °R (•- «* ,/ ) 4d 4d for D 2 - nJ 2 D 2 > D| > l R T» ~T 14 ) b' a a 3k( D D 2 ) 0 tmax rmax r x ] / (9d) 2 - r r)/ (D 2d Max(D ,D a n Z (1 T + D) \V2 ) 0 T p °° (l - 4k2 4 / ) 2 - D 2 < ijjj-, min ,Dg) < max (D 2 ,D for |D| | ( ) > ^ 2 b' a a (d + o tmax rmax t y (9e) 2 Z(l-rr)/ (D_ + D ) AV2 d oo o L|. T R 1 1 - 4d 2 for IDf-D|l <« D|<« D|<« Of course, in equation (9e), the latter two conditions imply the first. Notice that equation (9e) gives the same result as does equation (5). 3.2 Comparison of Upper Bound from Formulas to the Exact Upper Bound The results obtained from equations (9) for each case in section 2 is given in the last column of table 2.1. The small letter in parentheses in column 18 of table 2.1 indicates which set of conditions of equations (9) is applicable for each case. There is a minimum of one case for each set of conditions in equations (9). We notice that equations (9) give results which differ from the exact results by more than 20 dB in some cases. However, we further notice that equations (9) always give results which are greater than the results found using the exact far fields. 3.3 Conclusion We conclude that equations (9) can be used to obtain an upper bound to the coupling between two antennas neglecting multiple reflections, provided that coupling does not take place through the main beam of either antenna and provided that the to integrate k and ky integration range we have chosen is valid, i.e., it is adequate x only over the range -k ( - D )/2d to +k(Dy + Dp /2d and we limit d to the range D j p ) 2 (Dj + D )/2 < d < (Dy + Dp / A . It has been shown that this range of integration is p adequate for determining side-lobe coupling for nearly all realistic microwave antennas [1]. 15 4. Summary and Concluding Remarks In this report we have performed a feasibility study to see if it is possible to get an estimate of the maximum coupling between two antennas when the details of the far-field amplitude and phase are unknown. We conclude from a study using hypothetical uniformly illuminated circular antennas that approximation-1 and approximation-2 are especially useful methods of estimating the maximum coupling of two antennas when the separation distance is less than one- tenth of a Rayleigh dis- tance. Approximation-1 is limited to uniformly illuminated circular aperture antennas. Approximation-2 replaces the dot product of the fields in equation (1) by the product of the maximum field magnitudes in the general direction of the separation axis multiplied by cosine functions of k and k^ and is an approximation applicable to x general antennas. Approximation-2 is used in computer program ENVLP, which is documented in Appendix B. For distances greater than about one-quarter of a mutual Rayleigh distance, the far-field formula [equation used with (3)] Ay^ and A RMAX gives an increasingly accurate estimate of the maximum coupling between two antennas that is generally more accurate than the estimates provided by approximation-1 and approximation-2. For distances between one-tenth and one-quarter of a mutual Rayleigh distance approximation-2 and the far-field formula give equally reasonable estimates. For separation distances less than about one-tenth of a mutual Rayleigh distance approximation-2, i.e., program ENVLP, gives substantially more accurate values of coup! i ng. In section 3 we derived a set of expressions [equations (9)] that allow one to determine an upper bound to the coupling between two antennas at arbitrary separation distance. For many applications the upper-bound expression may be adequate, especially when the requirement is simply that coupling lies below a given value. The appeal of approximation-2 and the correspond!’ ng computer program ENVLP, is also their simplicity and the few input parameters that they require. In particular, the near-field coupling is computed between arbitrary antennas from a mere knowledge of (1) the separation distance of the antennas, (2) the side-lobe level of each antenna along the axis of separation, and (3) the radius of each antenna. Having obtained encouraging results for approximation-2 and the upper-bound equations (9) for hypothetical antennas we suggest that these results be tested experimentally using real antennas, and that a similar approximation be formulated and tested for the computer program CUPLZ. 16 The authors wish to thank Carl Stubenrauch of MBS for the original suggestion of using approximate far fields as well as for helpful discussions and suggestions throughout the work. Helpful discussions were also held with Allen Newell and Andrew Repjar of MBS. In addition, we wish to gratefully acknowledge the Electromagnetic Compatabil ity Analysis Center of the Defense Department for its financial support of this work. 17 5. References [1] Yaghjian, A. D. Efficient computation of antenna coupling and fields within the near-field region, IEEE Trans. Antennas Rropag. AP-30, pp. 113-128; January 1982. [2] Stubenrauch, C. F.; Yaghjian, A. D. Determination of mutual coupling between co- sited microwave antennas and calculation of near-zone electric field, Nat. Bur. Stand. (U.S.) NBSIR 80-1620, pp. 97-104, 113-114, 126-150; June 1981. [3] Johnson, C. C. Field and wave electrodynamics. New York: McGraw-Hill, Sec. 10.5; 1965. [4] Johnson, R. C. ; Ecker, H. A.; Hollis, J. S. Determination of far-field patterns from near-field measurements, Proc. IEEE 61(2), pp. 1668-1694; December 1973. [5] Newell, A. C.; Crawford, M. L. Planar near-field measurements on high performance array antennas, Nat. Bur. Stand. (U.S.) NBSIR 74-380; July 1974. [6] Kerns, D. M. Plane-wave scatteri ng-matrix theory of antennas and antenna-antenna interactions: Formulation and applications, J. Res. Nat. Bur. Stand. (U.S.) 80B, p. 24; January-March 1976. 18 . 6. List of Figures Figure 1. Schematic of two antennas arbitrarily oriented and separated in free space. Figure 2. The exact far-field pattern of the x-component of electric field as a function of 0 with X = 3 cm and a = 10 cm for a uniformly illuminated circular antenna. Figure 3. The approximation-1 far-field pattern of the x-component of electric field as a function of 0 with A = 3 cm and a = 10 cm for a uni formly illuminated :ircular antenna. = 0° Figure 4. Ihe approximation-2 far-field pattern as a function of 0, with Figure 13. 2 coupling quotient for the X0 cut using the exact, approximation-1 and Droximation-2 far fields: Case Gamma-1. Figure 14. 2 coupling quotient for the X0 cut using the exact, approximati on-1 and Droximation-2 far fields: Case Gamma-2. Figure 15. 2 coupling quotient for the X0 cut using the exact. approximation-1 and Droximation-2 far fields: Case Delta-1. Figure 16. ? coupling quotient for the X0 cut using the exact, approximation-1 and Droximation-2 far fields: Case Delta-2. Figure 17. d coupling quotient for the X0 cut using the. exact. approximation-1 and Droximation-2 far fields: Case Epsilon-1. Figure 18. d coupling quotient for the X0 cut using the exact. approximati on-1 and Droximation-2 far fields: Case Epsilon-2. Figure 19. 2 coupling quotient for the X0 cut using the exact, approximation-1 and Droximation-2 far fields: Case Epsilon-3. Figure 20. 5 coupling quotient for the X0 cut using the exact, approximation-1 and Droximation-2 far fields: Case Omega. 19 free in separated and oriented arbitrarily antennas two of Schematic space. 1. Figure 20 Figure 2. The exact far-field pattern of the x-component of electric field as a function of 0 with X = 3 cm and a = 10 cm for a uniformly illuminated circular antenna. 21 kass 20.94 Figure 3. The approximation-1 far-field pattern of the x-component of electric field as a function of 6 with X = 3 cm and a = 10 cm for a uniformly illuminated circular antenna. 22 ka = 20.94 = 0° Figure 4. The approximation-2 far-field pattern as a function of 6, with 23 ,z to the Figure 5. Eulerian angles ( 24 and approximation-1 a(E o oo exact, co the using CM OC O Alpha-1. O cut d XO Case CM o r X the for fields: r CO r quotient far CO N j 33.3A i GH; 10 © FREQ The 6. o o .m m Figure 25 and ooc approximation-1 exact, the using Alpha-4. cut XO Case the for fields: quotient far QH; 33.3A 10 - = coupling d approximation-2 FREQ The 7. „ ©i © n j« Figure 26 and approximation-1 exact, the using Alpha-7. cut XO Case the for fields: A • o quotient far 6.67 O •© © VCO = II II II r esc GHz X D 33.3 10 0.67A coupling = = © © © © approximation-2 & © © © © q7 d ® o d 80 © in d TO o TO © © N- © to CO K ® © CM TO TO FREQ = II II IS 8 I 8 8 8 ! 5 l The ' H K H ! CD C 8. _ O o XI & Figure 27 and ooc approximation-1 V- cm o 8 1 T~ e ©S o — cq <3=3 « 0 E E exact, K X <0 *5u © o & a a x a a the UJ < < ’ 1 : i CM using 9 • B tr Q Beta-1. © + cut © QK W XO Case CM o 1 X the — for r fields: CO i‘ quotient far is n si si CO c tr s oc 1° QM; O ^ 10 1 1 q = ° © O © © © © © © © ® coupling ® 9 approximation-2 ad «> id © « *o S ° © m 6 6 6 © 3*® 5 ** c© I* ® © CM CO m FREQ ‘J®- si a ii si V 8 i 5 I i t- H- K »- 8 1 1 i © ©* «> •*• The ffl “D C 9. _ o © jQ 0 Figure 28 X3 c «T3 ci o •r— +> fD E x o <*» s- £M a. Cl i I -4se OS oo "4 II 4-3U « © fO E E X X X <0 a> *s © o O) a a a 4-> x a a m < < eCD CM i to CM ^ 3 BE M Beta- Q 3 © + U “ o d X Case V • CM o £ 10 = 6.67A CL = © © © © © © © © © o ©7 3 approximation-2 d o . © • d d d d d o’ d o’ d o o ® o o <_> 'F* CO © K © CM CO tf) FREQ = F3* ” II II II i e i i s 1 T» a; r- H- » I- i I 1 a O •$. <©. <*. CD 29 o c fO ci o +j a fO Q e X o i~ CL CL to +> u X cn c • •f- oo tn 3 Beta- 4->3 o£E o o + X Case £ +-> cn c CL 3 approximation-2 o u O) * O © 3l. n a cn 30 -o e m Beta- Case fields: far Ql 10 approximation-2 = FREQ 31 r *0 c ra C! o 4-> ra E X tz o u Cl CL (O u XfO O) a> .c 4-> CD • i/) Gamma- +J 3 O o X Case Ga; +J O1. 4- fields: +J c A a> far +j 3.33 o 3 = c X Do o~> GHz 7 c . 5 18 0. = = A 3 approximation-2 33 o d . 3 u FREQ aj .= .c ro C cu _ © o 3 m CD 32 o c <13 cI o 4-> O xre QJ Q> CD • C CM ooc + Gamma- *->3 O o o X Case X CQJ +J o1. fields: 4->c QJ r— far +J o GHz 66.7A 3 Cr 20 CCD = = d Cl 3 approximation-2 FREQ O a Q> •o c - o o .O CO Q) 3i~ CD 33 A — o c fO c o M «3 E •r-X o s- CL CL <13 -MO XfC O) QJ er> c . to3 Delta- +J 3 o X Case o> O JC X +-> o fields: +J e o> •r far o 3 1688./ CT GHz ccn 10 CL = ~ 3 approximation-2 o d o FREQ QJ tr> O) . © © s- 83 3 I3> 34 o c ci o 4-> u xre OJ 0) M CT> • sr c\j to ac 3 Delta- Q +-> 4- U=3 o X Case •s— far 4-> = o A 3 b er GHz D 233.3 cn 10 c: = = 0.87A D. d 3 approximation-2 O FREQ = O a; o c CD . © o jQ 03 a> 3s~ CD 35 TJ c: to ci o +-> TO a. X o s_ CL Cl TO O XTO 0> t eCD £ */> Q Epsilon- +> 3 a o Case O X X a> sz 4-> oS_ fields: ?A 6 •Mc ®. QJ far = O 3 0Hz 7A Do cr 10 16. c = - ®J?A CL d 3 approximation-2 O FREQ S o 0J Jc •a „ © © jq a a> 3i- cn 36 TT E fO cl o +J fD cc E o X o a. Cl to 4-> O Xra QJ QJ JE +-> . CM en CO Epsilon- o fields: +-> 6.67A c D. approximation-2 6.87 3 d o o FREQ = a> c co _ o o 37 and QJE approximation-1 exact, the using Epsilon-3. cut XO Case the ®.®/A for fields: = A 0o quotient far GHz 133.3 10 ®.@?A = ss d FREQ coupling .= approximation-2 T3 The - © o M & 19. Figure 38 and approximation-1 exact, the using o£E Omega. cut X0 Case the for fields: quotient far coupling approximation-2 oGO The 20. _ o o JO CO Figure 39 and A in of Appendix A. Derivation of AyMAX RMAX Terms Antenna Gains and Relative Side-Lobe Levels In this section we derive the amplitudes, Ay and A of the approximate MAX RMAX , pattern for the transmitting and receiving antennas, respectively (see section 2.2.3) in terms of the antenna gains, Gy and G the side-lobe levels in dB, Sy and S the R , R , reflection coefficients of the antennas, r _ and r the characteristic admittances oT oRn , for the propagated mode in the waveguide feeds of the antennas, n and n and the Qy QR , wave impedance of free space, Z . Q The magnitude of the vector _f (which is just A^) in equation (1) of section 2.1 in the direction r_ is given by [1]: |E(r)|r * 1 * I fill I y^ (AD where E_(_r) is the electric field in the direction _r, and a is the amplitude of the Q mode of the waveguide feed antenna. is incident to the AMAX either AyMAX or ARMAX depending on whether the transmitting or receiving antenna is being considered. The relative side-lobe level S in the direction (Sy for the transmitting antenna and S £ R for the receiving direction), is: S = -20 log where £0 is the direction of the main beam. It is well known that the gain for an antenna, G (Gy for the transmitting antenna and G for the receiving antenna), is R given in dB by 2 Z 4tt r I— (^o^ | G = 10 log * (A3) 2 7 p. 0 input Assuming a single propagating mode in the waveguide feeding the antenna, the input an antenna, P-j can be expressed as power to npu y, [6]: 1 ( A4) input 7 where n is n for the transmitting antenna and n for the receiving antenna and Q Qy 0R r is r for the transmitting antenna and t for the receiving antenna. Q oT qR Substituting (A4) into (A3) for Pjp we find that the gain is t 2 2 4tt E (r ) r — -0 = 10 log ( A5) G 2 2 Z n a (1 r ) 0 0 o 0 40 — 1 ‘ ) |E_ ( r ) We solve for - r- r in (A5) to get a p l ol ( A6) Substitution of from (A2) into (A6) gives |£(l0 )j £(r^)jr V2 s,2° G/20 10 a 10 10 ( A7 -p^j- max I>.w) turns out to be simply Thus, A^ax /2 - = iqC CD 1 I GO ro \V2 —1 0 = 2 — a r t— O (A9a) tmax ( >) l ol| ((G - S )/20) \V2 r r = 2 K K A ^ r 10 ( A9b) RMAX oR 1 7 In summary, equations (A9) give the magnitude of the approximate electric far- field pattern of the transmitting and receiving antenna along their axis of separation in terms of the gain (in dB), side-lobe level (in dB) along the axis of separation, the antenna input reflection coefficient, the characteristic admittance of the waveguide feeding the antenna, and, of course, the impedance of free space. 41 . Appendix B. Documentation of ENVLP, the Computer Program to Estimate Coupling Between Two Antennas This appendix documents the program ENVLP which estimates the coupling loss between two antennas given their radii, separation distance, and side- lobe levels along their axis of separation. Subroutines used by this program have also been used by CUPLNF (except for CRTPLT3) and are documented in NBSIR 80-1630[2] B.l General Overview of Computer Program ENVLP The techniques used by ENVLP for evaluating the coupling loss are basically the same as those used by the computer program CUPLNF [2]. The flow chart for ENVLP is presented below to give the reader a general understanding of the program package. Purpose : To calculate an estimate of the coupling loss between two antennas given the radii of the antennas, their separation distance, and their side-lobe levels along their separation axis. Method : Evaluate equation (1) of the main text using approximation-2 (see section 2.2.3) to estimate the dot product of the far fields. General Discussion : The main program ENVLP divides into the following sections: 1. Genera] information, 2. Specification statements, 3. Definition and reading of input data, 4. Far-field coupling computation, 5. Limits of integration and number of integration points, 6. Filling the input matrices to the FFT subroutine FOURT, and 7. Calculation of maximum, minimum, mean coupling; printout and plotting of the XO and YO cuts. General Information : This section provides a general description to the program user of what the program does and also defines the more important parameters of the program. A reading of this section will be sufficient for most users to begin using the program. Specification Statements : This section dimensions arrays, places arrays in correnon, and declares complex variables. 42 FLOW CHART FOR PROGRAM ENVLP . . . Definition and Reading of Input Data : This section defines basic constants such as the speed of light and reads from data cards the antenna parameters. The required data cards are: Card 1 Col. 1-10 An alphanumeric identifier of the user's choice used to identify the case being computed; the identifier appears at the top of the printout, (HEAD (2)). Card 2 Col. 5 The maximum value that the user wishes XLIM to assume; it should be 1 or 2. XLIM adjusts the integration range (see below), ( ILIMAX) Col. 10 The maximum value that the user wishes BFAC to assume; it should be 1 or 2. BFAC adjusts the integration step size (see below) , (IBFMAX). Card 3 Col. 1- 10 The separation distance in meters, (ZO). Col 11 -20 The frequency in Hz. Card 4 Col. 1- 10 The radius in meters of the smallest sphere circumscribing the effective transmitting antenna, (RADT). Col. 11 -20 The gain of the transmitting antenna in dB, ( GAINT ) Col 21 -30 The relative side-lobe level in dB of the transmitting antenna along the axis of separation (ST). Col 31 -40 The real part of the input reflection coefficient, GAMMAOT, for the transmitting antenna in free space. Col 41 -50 The imaginary part of the input reflection coefficient, GAMMAOT, for the transmitting antenna in free space. Card 5 Col 1- 10 The radius in meters of the smallest sphere circumscribing the effective receiving antenna, (RADR). Col 11 1 CM o The gain of the receiving antenna in dB, (GAINR). Col. 21 -30 The relative side-lobe level in dB of the receiving antenna along the separation axis (SR). Col 31- 40 The real part of the input reflection coefficient, GAMMAOR, for the receiving antenna in free space. Col. 41 -50 The imaginary part of the input reflection coefficient, GAMMAOR, for the receiving antenna in free space. Col. 51 -60 The real part of the reflection coefficient, GAMMALR, for the passive termination on the receiving antenna. Col. 61 -70 The imaginary part of the reflection coefficient, GAMMALR, for the passive termination on the receiving antenna. Card 6 Col. 1- 10 The characteri stic admittance of the propagating mode in the waveguide feed of the transmitting antenna, (ETAT). If 44 . ETAT=0, the program will set ETAT=1/CAPZ0. Col. 11-20 The characteristic admittance of the propagating mode in the waveguide feed of the receiving antenna, (ETAR). If ETAR=0, the program will set ETAR=1/CAPZ0. Far-Field Coupling Computation : This section computes the coupling between two antennas directly from their far fields if their separation distance is greater than a mutual Rayleigh distance (equal to the square of the sum of the effective antenna diameters, divided by the wavelength; see [1]). Limits of Integration and Number of Integration Points : For XLIM=1, the limits of integration are computed using the specification of section 2.1 of the main text. These limits can be doubled to see if a wide enough integration range has been included by setting XLIM=2. Strictly speaking, approximation-2 (see section 2.2.3) is only good for the XLXM=1 integration range. However, if the results of the XLIM=2 integration range for the mean coupling differs by a huge amount (e.g., 10 dB) from the mean coupling result for XLIM=1, it is ques- tionable whether ENVLP gives a good estimate of the coupling loss between the two antennas being considered. This section also chooses the integration increment size small enough to prevent aliasing. The increment size can be reduced by increasing BFAC. Normally, BFAC is set equal to 1 and convergence is tested by halving the integration increment size, i.e., letting BFAC=2. Filling the Input Matrices to the FFT Subroutine FOURT : This section calculates the dot product of the far fields in a square array using approximation-2. The dot products are then summed in the k (k ) direction for each value of k (k ) y x x y and placed in the array AX(AY). FOURT then performs a one dimensional FFT on the array AX (AY) to obtain the values of the coupling quotient along the XO(YO) cut. Calculation of Maximum, Minimum, Mean Coupling; Printout and Plotting of the XO and YO Cuts : This section computes the maximum, minimum, mean, and RMS mean coupling for the XO and YO cuts. It further computes the standard deviation in the coupling for each cut. It prints out the values of the coupling quotient for each cut and then plots the values of the coupling quotient for each cut from -(DIAMR+DIAMT) to +(DIAMR+DIAMT) 45 Symbol Dictionary (in alphabetical order): ABL Intermediate value for defining the range of k /k and k /k. The x y range of ABL beyond XKLIM is zero filled. ACL CUT A real array used to store the magnitude of the coupling quotient along XO and YO cuts. AMAX The maximum value of the coupling quotient for either the XO or YO cut. AMEAN The mean value of the coupling quotient for either the XO or YO cut. AM!NX,(AMINY) The minimum value of the coupling quotient for the XO(YO) cut. ARMAX, (ATMAX) The side-lobe level for the receiving (transmitting) antenna. AX, (AY) Complex arrays used to store first the far-field product then the coupling quotient along the XO(YO) cut. (A1,A2) ,(B1,B2) - The limits of integration of k /k and k /k respectively. x y BFAC Variable which adjusts the integration increments and should be about 1 or 2; making BFAC larger tests for convergence by making the integration increments proportionally smaller. C The increment size for the XO and YO cuts when the far-field formula is used. CAPZO The wave impedance of free space. CEE The speed of light in a vacuum in nmeters/second. CLCUTXA( CLCUTYA) = The values of the coupling quotients in the XO(YO) cuts stored for plotting. COEF 2.*PI*FMM/ETA/CAPZ0/XK for the far-field formula, -FMM*C1*C2/ETA/ CAPZO otherwise. Cl C2 The k /k and k^/k increments respectively. s x DIAMR, (DIAMT) Twice the larger of RADR (RADT) or WAVLGTH. DIAMSUM DIAMR plus DIAMT. DKOK The approximate k /k and k /k increments, x y DX,(DY) The increments in XO(YO) over which the coupling quotient is computed by the FFT. etat,(etar) The characteristic admittance for the propagated mode for the waveguide feed of the transmitting (receiving) antenna. FDOTFP The dot product of the electric far-field pattern of the two antennas. FMM The mismatch factor, l/(l-r r£, for the receiving antenna. 0 FREQ Frequency in Hz. 46 M ) . FX,FY Intermediate complex variables. GAINT, (GAINR) The gain in dB of the transmitting (receiving) antenna. GAMMALR,GAMMAOR,GAMMAOT = The reflection coefficients of the receiving load. receiving antenna, and transmitting antenna, respectively. HEAD Integer array identifier. IBFAC Loop index for varying BFAC. IBFMAX The maximum value assumed by BFAC; it should equal to 1 or 2. ILIMAX The maximum value assumed by XLIM; it should be 1 or 2. I XL I Loop index for varying XLIM. J1,J2 Loop indices used in the filling of AX (AY) before the FFT. Ml, M2 Loop indices used in completing the coupling quotient computation after the FFT. N The number of points in an XO(YO) cut when the far-field formula is used. NN1,NN2 Integer arrays of dimensions one used in the call to the FFT subroutine FOURT. NSX , ( NSY The number of values of the coupling quotients to be plotted in the XO(YO) cut. N1,N2 The number of k and k integration points respectively, x y N1MAX , (N2MAX) Integers determining the maximum of the XO(YO) range over which the coupling is printed. M1MIM, (M2MIN) = Integers determining the minimum of the XO(YO) range over which the coupling is printed. N10,N20 Intermediate integers used to determine ( N1MIN ,N1MAX) and ( N2MI N , N2MAX ) PI it = 3.14159... RADR, (RADT) The radius of the smallest sphere circumscribing the receiving (transmitting) antenna in meters. RMS The RMS mean coupling quotient for the XO or YO cut. RO The total distance between the two antennas = (X 2 + Y 2 + l z fa. SDEV The standard deviation in the coupling quotient for the XO or YO cut. SR, ST The relative side-lobe levels, in dB, of the receiving and transmitting antennas along their axis of separation. = SUM , SUMRMS , SUMSQ Intermediate variables used to calculate AMEAN, RMS and SDEV, respectively. SUM2 Summation variable used in the filling of the AX and AY matrices. TEMP An intermediate variable used to calculate ATMAX and ARMAX. 47 ) . . TSUM21 Summation variable used to compute the coupling quotient at XO = 0 by summing directly without the use of the FFT. WAVLGTH = Wavelength in meters. WORK = Complex array required by the FFT subroutine FOURT. = X , ( Y Array containing the abscissa value for the plots of the XO(YO) cuts. XK = 2tt / X XKLIM - Variable which limits the range of k /k and k /k integration when x y its value is less than XKMAX. XKMAX — An upper limit (less than 1.0 and usually chosen at .9) on XKLIM; except for close antennas XKLIM will usually be less than XKMAX. XKMIN = Sum of the diameters of the two antennas divided by their separation distance. XKOK - The square root of the sum of the squares of XKXOK and XKYOK XKXQK ,XKY0K k /k and k /k, respectively, x y XLIM Variable used for adjusting XKLIM; making XLIM larger tests to see if a large enough integration range has been included; XLIM should equal 1 or 2. XNX ,XNY ,XNZ = Variables used for incrementing k /k, k /k, and y/k, respectively. ” y XYO = Intermediate variable used for calculating RO. antenna in the XO, Y0 S ZQ X S Y,Z coordinates of the origin of the receiving mutual coupling coordinate system of the transmitting antenna; specifically ZO is the separation distance. ID, IFL, IU, NOFRAME = Variables used by the CRT plotting routine CRTPLT3, which is a specialized routine for the NOAA/NBS CYBER 170/750. Subroutines Not Available in the FORTRAN Library : FOURT (a standard FFT subroutine documented with program CUPLNF in NBSIR 80-1630C2].) PLT120R (printout plotting routine documented in NBSIR 8Q-1630[2].) CRTPLT3 (a CRT plotting routine specifically written for the NOAA/NBS Cyber 170/750; we suggest you substitute your own subroutine). Note : If the electric far field, including phase as well as amplitude, is available use the program CUPLNF documented in NBSIR 80-1630[2]. B.2 Computer Code and Sample Output A copy of the computer code for program ENVLP and a sample output for the case BETA-1 are found below. 48 • J ! « o flu ©O' © o 1 1 © mJ 1 o mJ ft 1 • ft © ft ft •> © & I © © z z z ft • I z o © © Z X •< © <4 COMMON SEPARATION © © z z © © i «k ft ft ft © © o z © © © © © © X ft’ mJ © <=£ © C— © > c=3 © © CJ z © O U © * O © X (S=9 « z «« &L © ft €3 © ft X a. ZD © © ft <» © 40 O 1 tag? ft ft ft ft © ft © to OFTFPMfNEO CNJ © ft © © Z iii z ft © ft ft © ft z ft z © © z z © X ft ft ft X ft X in ft ft 40 © ft THF FOU»T. © © o © M. <©> <3 < © © © © ft ft © > ft •S3 • z © z c X ft c ft © y. z © a © THE ft © © © c a. © © © 40 I ft z z - ft ft o ft ft mJ ft X ft © © c © ft X ft «s> ft <80 X z © ft ft ft APF a ft < =« UCt-tC^LUCCUUCL’CULLCLJUL-CUCCtLUt'UUUUCLJUL'UU^UUCCLCULCCOC’ULO © m o © jp o « • it iCVi cu ^ * IT ; — » • ; M ©ft G. & © 4P «a ® 2 o ft 2. ft I Ui 3- O 40 •k S3 ft 2 ! ft t=* c=*> «» : a> SO ft 2 ft X •i m o 1 © > 1 z° © c- ft 4\> sc © 15 z 2 OC £=3 2 1 1S> 40 ft 2 flD © &=j ft © £ yj ft 15 1 ©3. 40) ec ag ® ft- 1 40 ft &=* tss ft® ft - o «*S ft O r ft «=*. *9. • ft <*a © •' . ft <*• ft O. ft ft z LL © <* c ft 40 e ft «a <39 c=- ttes. ft ft ft ft Rs* ft ft © <$> e «S!> 8=3 : <^B. <1 «a as »- ft 5=^ 8=0 ! a**. f. ft ft ft ft- X z > “ aJ ft ft Z CD k « © z © 3 c c © © © >» <33 “Sis, 0=* : s~ >pv Ok 4f X 40« 2 40 40 ^* 40. ft=« w 80 0 X <=5 © «e c* X ft ft 2 a 2 40 ft ft ft ft 2 z &® <« ft Ift ft ft « Ol U, ft- Ux y. u. ft- X ft- © o ft N 2 ft- Jr © 3 2 ft ft © «® « u. ft 2 2 © © 2 ft O. z X — y 8- ^°» @> 2 Qa. 40 ft— lx ^ e=> <30 z — ft ^ 2 > a 2 ft ft St © a ft © X ft X © *«. ft ft © 2 « *« 53. ** u- « x o=> <* x <5® z - X £=5 ©' t tut © © © © © 4=, © © © © ©©©©©©©©©>©©©© I 1 i 1 1 i ; 1 i ! ! ! 1 : j 1 1 j ! i ! 1 © m © cn <_> 40. © & °4j (ft CL ! f** Ov & © j © i j j j sq l » m UJ© C*I o. o M CL FPP JO c PHONF © AN?) CVJ © fV © f\j «\- N4MF h ft & Q (T C IT «s o © O CL I I © <* © X t- X fife et <*a x «a X Z X X UJ XXX o e®» « © ^ X « X © «a © <1 CnWIMMS © r C C r r © c © © © • l j . 1 ! ! ?5 u m !© ir © ir r m ,* i 1 - 4 I I*- "i i I i i 1 51 i ©4b a. X fn © (SJ o *** x » . 9 i o 40 ^ X X* © © X* X i 2 .=3 • © 2 O NUC «S» ^ X JMS 1 40 X £ _J i « c X ®~» © X u © 9 X C\J $ 4b ab CO 3i Lb • CL. 1 ®*fc. 4£3 «« <*» «*=> • -^ 4- fiC X • x x e *l & X X r x 0 © x « a u <^ 1 *—» a x 40 9=» X ©- CO © © x ^ * <« £H fe • X u. • X © CO X f* e ^ x < X a X <1 x & 2 x g CL X X ^ X 4r X 1 & L C © u. c^ s *- X Lb X x XU X a a 2 © © X «a X < © X X c c =j a x <® — C X & c z © > X « z « <* u u « a *9 — © 9- Z «. o © ^ to x ^ w « < ft «k = «9 w * 4b © * <* © 4^ 40• X * © X 9- X © X © ^ f-» 4/ » «*=* © «=i X >*. X z x a x X X © © © *. • *** X- « X ^ x a x , • ©*’-'© 9 P* & r* ** e— 40 X © X V «& * 9 © © U 9* X x c u x *-• x *-* U & X* • • C fL CL r" & 9* W X 9* *=* 0~ Q. u u © cr 0 “9 w « © © £ 9 x X co "> Co IT X O £k © C © © © © ©I © ir> if ©* 52 . i . © a. o IVI Ca.> © © 2& fto > ft. © 1 o © © ft <** IS Z © «es irv z c> ft v> ®L X © © *“ e& a €S ©> © © © © © CM-) C L O © © © © «r r. *r '© OL a . K CSAJ fSJ a* fV. is 4SJ N i is 53 T ; o<*- oot © z 4^ © « Z • z X @2 0“ CA4 I to ft NP A ft C (VJ <**=• X X 40 to © c i © en ft © x= © to to © to © a- Z m ft 40 o er <<1 to* z J to ft to to ft z 40 (to ft Z © 40 as ft i z » to ft @2 0=s at ; © ft X 3 ® z ' & to tv o at IO V t!*5> 1 *a <& ft at © o * i Z G © © «0 © # at z H* z to» ft> o to © ft ft » ft to z ft 9 » 1 3 > a (f© at to H- & SL « to to a ft X m. C=i * 0 0 e=a 11N^»-N1- tS <£3 to at <» z at 4^ ft 8 «r c ft © ® to ft h fvj * X <« 4T> © © £*•*==> «** 1 o. ®, 5=0 I 5=9 to ft ^ O' 6 0 » • e= z o» C >- ft K a ft * ft z © to to to a © F0 «s» y> C=:' ft • a. G • 0 © X © © 40 © © . © X to c> to 4/ © © ft <1 # X 9 <» (VJ tv. cv. } a* to- © 0= © y. c. © *- e © a •*0 «o to to OB* ft ft C x c « © © © © ir» a e o I© to «N- fr> fT) fw m |W «n ;m HI !«*> 54 , ; ' © m © o«r ft • ID f^J «VJ sj ID u. X u . ID ft ft # * ft ID 1 //» > fiD ID a o ft ft ! ** «* # ft *- ft 1 O Ol ft y» ID * © 1 Uj C ft> © © 01 1 in «» 1 ID © 1 © ft © EVEN ft in en e\ (Vj ID © ft ID ID ®n ft ft • 9 ft z | Z ID ft Z uj z Z S\J ru it z j < CC z ID) ft ID ID ft ft * BE © I ID ft s m y» ID in ! iS © X ft ft in ID ID * « X tor 8 ft If»'5*»*THFV c D Z ID 5D U ft X C ft X © ft X ft > ft ft. id ft X © 9 ft «a <**. •* 1 ft < «<*»• ft * < Z U CV Z it' © -j 2 fr- HX 2 g X Z ' 1 X • • © I IT W r+ «»> «* — IT to* ft © If IT © ft cr ifi ^ f. oc © & CV fv CV i 1 i © 1 © © g C g g g u g 1 i i i 1 i 1 j j i i j ' ! | j i i «v IfN r' © o i . © - o s- ! If. if* *> D. ft i I oc cc, a o i cn .fr> n «*i ^ icn mi *" m m i l ! j 1 1 ! 55 i 1 CD © ft W ft AC BF a * IN* 2 O • © XL j © © 8 © © © FOP © in © <3 a i — 8 «t © s — 1 ft a © ^4 , ©) 0 ft ft i* ft 0 ft ft ! ANALYSIS 1 0 «* 1 £3=B ft-. 49 X © © 0 8 0 © ft x # ^ ft ft © o 8 “8 U. 1 to «» «^ ^ a ^ If © 8 ^ at ! «=»<=) St i in ft «»; to © © a a : X 2 t I «3 • to ft © a a «=* x «-> ® : - X $ to # X **> «=* X 2 © > * : x m ft © © «k K « .to «=« X X <*“> ^=« <** »Nf> <£ 9 x C X * # © ft 8 c © DH dj in x © :« l ^ j x in «8 “5 x X } to 4t FflirNS. a © e «£ © 2 8- 8- © 8 to © r> © •-« <® «** to 4n © © © © ft © in MUTNtPFaru.rooooi «n © Q X X © * © « ©=* ^ ^ y, ® X 8>- s « *=t © X 3 O X 2 a a a © S> o m> y 2 PPTMTOUT <33 ©' 1 X • I ft ^ in © «w ft • 9 rsj « » u m a C fsj »« to e § © m ^ 8ft© £ C © © ft in # 2 0 0 W> ^ ft ft 2 H O «zJ a ft e> is> «=j ft in ir ext; X 2 ft 4> *=• © © r © c © S 2 •> fc cvj y u ^<2 a «» © r* x c c <* a a y> ft © 0-* © «a m ^ «a @i (V. ®* ®» ft «w ft to* is, o ^ ^ tr> © « ft THAN «*- «* «k <«. 1 ^ © a. © 0 in | « ^ * 8 ^ a © ^ t b m i • 4/ • » * O ~ •-» a~ e=» «r X • © *» © • ft S « x a ft a ft IT © © 9* © 0 a i/ ft ©> •» ^ X ft ft © St ft x © “> ^ s in v~} ft ft ft ^ © ft ?! ** o a z> %- in in ^ r>© © © © SV ^) i. t 0 ©©• a D-in in to ^ ^ «O ft © c 8 < <£ r* © «« e 0 St ft s : in if- * « in. fv fw m 2 2 x © X 2 m X © c — 2 x it x 7 o > 2 -Ni . >- ft 8 a m «e z © S5 X X 8 8 60 XXXI « 8 a © © © a m » '© I © *v © © om lf> © © s=a © f-i ft © IT ^ Sfv 4C r® *. u © . © © © © fy. ,«t «^> e 56 , 1 J to o (Vi I X © to ro ; OL to ** • y- x <**> «s*» ft Ui > X X 9=4 2 Z C ft ft 9k VMh h © © 4 X Z iti ft * ft 4 © ^ x > O © fo © © U to «r. 4 4 tfVJ to U. _i 2 -J u- X y=* x Z Z 4 A » R 3 X to x to ft ft I z t 2 O >-3 3 % % IV » LU - t= O >- X X X X -3 X •vi 3 C O l/* ^ ^ U'i o» »SJ (H> C cf U. 2 Z Z Z x z *- Z ft ft • • © I 3 I u to u «si © © © © -j -» «j at 33 rv tr* IT IT u — (tan to to »* 8 8 0 I El © a © 4/1 4 to • Z ft ft ft X © X © © to 4 o • a a ** 4 X to 4 1 to CO C ^ & o © © o •Si a © ft z © Z to X X © ft ft ft ft ^ 4 44 o to CD 4 4 © c— ft © CQ © © c- X to 4 to to to X to x X to • • a • to to X © «» ft to X o ft © © to to X to to to -* • 4 H ar 4 to © ft to x 8 • 0 1 ^ ** * © X to 4a to 4 c 4 X © Z X © © o © 1 4 I to • X Z **» © 4 X to a to z C30 © to V to X a a a • ^ h <4 to 0=1 © to> to 4 to 4 to to X to 4 to to to to «Sb to to c. ft ft ft ^ C 2 C # X 4 o a Z © # Z 4 a i * X © ft «j ^ • • to Z 2 © fit to to • ® ft Z z 4 X • c 4 S* 4 4 4 4 X 4 > ft — O z •— <3 © X 4 > z to O z to 4 to 1 4 X U. X 4 2x X X X CNj c trvi © '© 4) Xi V I cr © 4* •4* ,«r S| 57 ' o mK 9 © o «SJ OL © «0 S «J X 1 K a. Ok a K l a c Ci a i 1 fir © cn fr ,*r ir% 58 i I o o o o• PAGE o o 55. Of. 0 12. fit o © o' Pl/11/12. © © o o © o oo o• 1.5-52R © o © © © © © © • • lOSOFR © o Q • rYPFO N-© ©& <31 s. 2. S SL © © © © © w (f*> IT. • # © Q- m rri 0 0 K- tt © , © © © ! © © CNJ 4 (VJ Oi <# «-• © I • KN © | • INJ | rsj c © ** © ©, « a , 0 1 • z Z ft fNVto o « < « h 4 _ u. tt © M«p © e © © © is. r,( © u. I & n u. u. O r — . © © • IN* L ^ 4U C © f^ U- O © C AM a. © F-i * • # °i* | >= i .1 OOOOOOOOOOOOOOOOOOOOOOOOOOOOO I Ui I I I | | | | I | • I I I I | • | I | | | I I I I I I I CIUIt •Mssortfn^-«*'tr>i>4 i©©s«4.©.e.-»(rs(vi©«<»'S,®tr lF,«>f-ir.#©or»©r»,-«.c«rfn 8 « • 1 8 8 0 8 s 8 Vj 1 o ^ f* # (« W mw (n # (11 fn fn #n ^ fn ^ <»> < fn m ft fn ft P P P « © © © © OOOOOGOO © o © © o © © o o © © © © O O G © 1 yy 0 9 9 8 8 9 8 i 8 9 8 0 9 II 0 0 9 0 0 9 1 9 e ! » i 1 9 9 1 i Ul ut u UJ fey UJ yy uw u> UJ UJ Uj UJ ft ft ft UJ UJ UJ UJ ft ft ft ft ft ft ft ft UJ © © . © , • If* IP i 8 * 9 8 0 8 i ! 9 «sj 1 fV8 © & © © c=« © + # ^ in m fn fn fn uj ft Ut 9 9 9 s 0 8 0 8 e 0 0 0 0 9 0 0 1) 9 » 0 1 9 3 9 « 8 1 1 i ! ss SL s© 9 feW few ft few m m : yj U u U UJ UJ UJ u u ft ft ft ft ft ft ft UJ ft- ft ft ft ft ft uj METERS (fvj 5=0 (Ph 0=4 gs aa m , © r> ® © IT © ©«-<» <# «# (Vi & (TS © fv © p^ fn © ft O IP P ® ( yy i ft at tfN ® © m IPs © fw ® «* © O © IP * m ® r* 0-4 « m o IP ft ff © r no P »i o © rv © BT o 0=0 s •*• © & ® a*, s sf «# c c ft CD C © ® P P IP fv w Mj £ z 3 m m fn P* ^ P=4 ^ *=4 © © * ^ fy p- P» ft © © ft no (VJ •c © p » & ft ft © 9 «V r-> fSJ «-o ^ (Si ru (Vi «n csi gv 1 2 1 1 » 8 0 9 1 • i er. an ft ft «fc 1,00000 «sj no © ^ ^0 ^ fn fn IPS fn fr . «?00O0 © 8 9 6 • 9 9 1 1 i l 1 » 1 | ft U- £ (n o & fc=> fn fn uj UJ a. c=3 «D © fw © • s. E. © © ft >- ® IfV & ft y* 5^ c~ 4T «• Cfc fw fVi P» G © fV =4, ®> G *-4 0=^ P-4 IS p* 0^ ft p* G o © ft =4 0-4 0=00 0=4 m OL © © & & 2 OC z «g 40 ft cu y > ( © ' (o=* «=« S3 ft ft «> ft c #- ft ft®* O 6- c U* v o z c. 5V e— © u. SC 4«j GVj 2 ft © 2 © £ © ft c e &=* 1 ft O «- ft- ea ft Ux U bt X Z Si —j e- ft ft ft * # © 2 @ ft ft X Z ffi ft & ft ft < ft •tf >* ft © c ft ter C ^ u «J ft u X ^-j 4T XX© > 3 ft 5=5 X ft C (M c ft ft o X »-• ^ i c ft ft ft c te. ft C c © ft © »»>K K D IS MTCPHFUH I i ff> 4CCLZt^C c z. o z <* K I CVj © X- St U 40 <38 QD 2 9- X 0 2 CSj 3B O © 3 O Ql «• CO x« 2 ^ *- ! © © X» e- U. © ec yj <3fi © 2 u» ® Ql aJ ©w 2 X ut X X 4b (x x a X © O X» X> 8^> 2 © mJ © 2 e ft- to © & « i © u. X X © > £ X aJ X & CO « & <33 u. © © ^=6 X ®5> «b W ® £ U K 2 © 3 # i © © Cj ^ © X 2 © © © ft- x «J 40 X © 4b ^ — © X © — “9 ® a « © sr <** <& c. C X to e © <23 © © « -a*- CO •> a. &t @s U- 9- U. C 0 ll X X * * — 0 X X- © x © © © to 2 * CO u a x a © © © «J ^ ^ © x- X 0 CO o c — 2 -> • i >• «»• © V) *- 3 X tv ta z © c ft ft © K K ft ~7 ft- tar X © ft © © ft ft © © X © tar 3 © X ft! c 3 S ft «2 ft- © IS ^ © © ft taJ © 3 O 3 ^ *** X oc S> ft £S, ft © <33. © © IS/ rsj © ft- © <2 © © © X z ft <** «9 •» e» © m. ft ft a & ft z Z ft taw © P 2 ft © © © v. X z X &&» v-» @» ft ta» X © © X ft- e=» * V* © ft 3 X a =3 a b- _» taJ (NJ X > © z ft X z X c ft X 31 X © ft -7 <3 tar tar tar tar ft H ft. Ol o-d ft-" z ta O 0 © « »= K 40 *- “a &= SL to e*= -a to ft y. « ft to to- to to * ^ U I 0= «* — © <*»> G ft z. a* a* a. a- a- © ® 2 & li 8 <$> O 3 -T © © ® S3 © cc © ki F=S a Xr &> 9 ®=3 © y y ® s» H «x & a Ok «S y* «• <3x c- c 2 «r a: «» 90 © y K © eg, K & © y © «s V) v) © K K c=H <& «> «a, »=3 2 2, y e V* • a 5K y y y a 2 X: x o e Vi « © © «= 3 9“ 5* 3 • e O ^ © e=j © © e=!l K <3=b (S* 5a M_ O- <2> a a © ^ © te» ^ X y 5a 3. K a. © <=j & c® K K X © a P- Sf' K <3 K © 6*v o- ©r ® © © y « k 2 «^* C« K i