Abstract:
An antenna system of any three dimension (3D) geometry and method for constructing said system with an array of basic antenna elements are described. An antenna system beam pattern is specified. The basic antenna element parameters including basic pattern and actual spacing location are measured. The measured parameterized array elements are injected in an exact array frame formula for any 3D array systems to form an array frame. Array calibration is performed by evaluating a dual frame to the array frame and the array control weights are determined based on the dual array frame and the specified system beam pattern. The antenna system and a software tool is then constructed in accordance with the antenna control weights. The present invention enables the high precision beam synthesis with high quality beams for array of any geometry. The present invention is capable of taking account various factors in antenna constructions together in a one-step approach. These factors include, for instance, mutual coupling, element spacing variation, element gain and basic pattern variation, antenna cable and feeds length variation (reflected in phase differences).

Description:
CROSS-REFERENCES TO RELATED APPLICATIONS 
   This application is related to and claims the benefit of U.S. Provisional Patent Application 60/772,909, entitled “Methods and Apparatus for Constructing General Wireless Antenna Systems”, filed on Feb. 13, 2006, which provisional application is incorporated in its entirety by reference into the present application. 

   FIELD OF INVENTION 
   This present invention generally relates to the construction of mobile and fixed link antenna systems with any array geometry. More particularly, this invention allows for the construction of array antenna systems with improved directivity, directive gain, and apertures and improved control over these quantities. 
   BACKGROUND OF THE INVENTION 
   Array antennas of different geometry/shapes are used in all wireless communication systems. An array antenna is typically constructed with a set of basic antenna elements arranged in an array. Such an array can often be linear—elements arranged in a line, can be planar, circular, cylindrical and spherical, etc. Depending on applications, each has their role in wireless communication applications. 
   Array beam pattern synthesis is to combine all elements in an array with complex weights so as to create beam patterns of the desired direction and shape. Effective synthesis and antenna construction has been studied for decades. There are many difficult factors that affect the effectiveness of an array antenna. For instance, mutual coupling among basic elements, element spacing variations in an array (which requires oftentimes high precision mechanical processes), element gain and basic pattern variations, array re-calibration (upon element failure in an array), and high quality beams. 
   Li has described a method and apparatus for constructing linear (including flat panel) wireless array antenna systems (U.S. Pat. No. 6,911,954 B2). The advantages include a beam synthesis method that incorporates all aforementioned factors into a one-step systematic approach. Mutual coupling, element spacing and gain variations, high quality beam can all be accounted for simultaneously. Array re-calibration (upon element failure detection) is also exceedingly easy in Li&#39;s method so that the array can still function in its maximum capacity as allowed by the physics exhibiting in the (remaining) array. 
   While Li&#39;s method is seen to bring in significant improvements to the construction of many useful array antenna systems, array antenna systems of other shapes, for instance, constructed over a cylindrical surface, a spherical surface, or any other three dimensional (3D) geometry (or General Wireless Array Antenna Systems) are still to be constructed in similar ways, taking advantages just as those in Li&#39;s method. 
   The construction of general wireless array antenna systems is considerably more complicated. Precise 3D beam synthesis has not been seen. Existing methods of construction of general array antennas are very limited. None has been able to incorporate aforementioned factors into the design simultaneously. Beam quality is also very restricted. 
   Thus, a better method for general array antenna constructions that provides better beam quality and directivity is needed, one that resolves the many factors in array beam synthesis, provides a precise 3D array synthesis, enables the array re-calibration and render existing construction method obsolete. 
   BRIEF SUMMARY OF THE INVENTION 
   The present invention is directed towards the above need—constructing array antennas with any array geometry. An antenna system in accordance with the present invention has improved directivity and directive gain and improved control over all array element parameters so that the best performance of the antenna system can be achieved. 
   A method in accordance with the present invention is a method for forming a beam for an antenna system that includes a plurality of antenna elements distributed over an arbitrary array geometry. The method includes specifying an antenna system radiation pattern function, determining element radiation pattern functions, determining the value of a set of spacing parameters, forming a frame from the element radiation pattern functions and the array geometry and finding a dual of the frame, and determining the element weight coefficients for the elements. The antenna system radiation pattern function describes the transmission or reception beam of the antenna system. The element radiation pattern functions each include a basic element pattern specification, a frequency of operation and at least one set of spacing parameters that specifies the location of the element in the antenna system. The frame that is formed from the element radiation pattern functions and geometrical relationship among elements arises from a condition, called the frame condition, imposed on the set of element radiation pattern functions and geometrical relationship among elements in an array. The element weight coefficient for each antenna element is based on the elements of a dual frame and the specified antenna system radiation pattern function. More particularly, the element weight coefficients result from the inner product of the dual frame with the specified antenna system radiation pattern function. The inner product is defined because of the frame condition imposed. 
   An apparatus in accordance with the present invention includes an antenna system and a software tool whose beam is formed by means of the method of the present invention. 
   One advantage of the present invention is that it enables the construction of array systems of any geometry precisely. 
   Another advantage of the present invention is that it can precisely include mutual coupling among elements into the description of each antenna element. 
   The other advantage is that non-uniform spacing of the elements is easily and precisely accommodated by the description of the antenna element pattern function. 
   Yet another advantage of the present invention is that real time re-calibration can be carried out if element gain changes or element failures or both are detected. This allows the array antenna to function at its best capacity allowable by the physics principle and allows mobile systems to function without having to replace or repair the antenna immediately. 
   Yet another advantage of the present invention is that computations involved in the method are quick so as to be suitable for re-calibration and reconfiguration of an antenna system after the system has been deployed. 
   Yet another advantage of the present invention is that element functions can include cable length variation, other circuit delays or other irregularities in the currents driving each array element. 

   
     BRIEF DESCRIPTION OF THE FIGURES 
     The mechanism and features, aspects and advantages of the present invention will become better understood with regard to the following description, appended claims, and accompanying drawings where: 
       FIG. 1  shows a flow chart of the steps, in accordance with the present invention, for constructing a general array antenna system of any array geometry. 
       FIG. 2  shows a schematic description of an arbitrary 3D array system in the xyz-coordinate system and in relation to the horizontal and elevation angles. 
       FIG. 3   a  shows a planar array weight control that creates the beam of  FIG. 3   b.    
       FIG. 3   b  shows a broadside beam pattern using the weights as shown in  FIG. 3   a.    
       FIG. 4   a  shows a non-uniformly (with a random variation to the uniform array) distributed planar array system. This is an example where non-uniform distribution can be easily compensated. 
       FIGS. 4   b - 4   c  show, respectively, the real and imaginary components of the weight control that naturally compensates the nonuniform distribution of the elements in  FIG. 4   a . The irregularity of the weight matrix is exactly what needed to compensate the non-uniform placement of the array. 
       FIG. 4   d  shows the result of the beam pattern. The result of the compensation through weight variations preserves almost entirely the original beam pattern of the uniform array (see  FIG. 3   b ), although the far-side sidelobes (below −40 db) are higher. Main beam is identical to the uniform array. 
       FIG. 5   a  shows a weight control matrix that compensates element gain changes. 
       FIG. 5   b  shows the resultant beam pattern of a planar array with element gain changes and compensated with corresponding weight control as shown in  FIG. 5   a.    
       FIG. 6   a  is an example of the weight control matrix a larger planar array. 
       FIG. 6   b  shows the corresponding much narrower beam pattern created by the larger array with weights shown in  FIG. 6   a.    
       FIG. 7   a  depicts a uniformly distributed cylindrical array uniformly distributed on the surface of a cylinder of radius 5.093λ. 
       FIGS. 7   b - 7   c  are associated weight control matrix in real and imaginary part, respectively, using this invention. 
       FIG. 7   d  shows the corresponding beam pattern using the invented method and so generated weight control as shown in  FIGS. 7   b - 7   c.    
       FIG. 8   a  shows an array system non-uniformly distributed on the surface of a cylinder of radius 5.093λ. 
       FIGS. 8   b - 8   c  are the compensated array control weight matrix in real and imaginary part, respectively. The irregular modifications to the weight matrix demonstrate the robustness of this invented method as the resultant beam pattern is kept largely unchanged (see  FIG. 8   d ). 
       FIG. 8   d  shows the corresponding beam pattern compensated using the invented method and weights as given in  FIGS. 8   b - 8   c.    
       FIG. 9  shows the would-be beam pattern if the cylindrical array as shown in  FIG. 7   a  is synthesized using conventional uniform/equal magnitude weights. The quality of the beam is clearly restricted. 
       FIG. 10  shows a sector-spherical array. 
   

   DETAILED DESCRIPTION OF THE INVENTION 
   The present invention relates to a method for constructing an improved array antenna system built over any geometry configuration. The method of the present invention avoids many of the drawbacks and approximations of existing cell and other array antenna systems by a new and rigorous approach, namely a frame theoretical approach. In a typical array antenna construction, the type of element of an antenna array is known, and the basic element pattern is approximated based on a model of the element. However, combination of the elemental radiation pattern to achieve any desired radiation pattern for the antenna array is limited in accuracy and controllability because of many simplifying assumptions must be made to make the problem tractable. The simplification approximation/assumption is particularly necessary in cases where an array geometry is irregular. Some simplifying assumptions include the regular spacing of elements, the allowable spacing of elements, such as half-wavelength, simplified basic element pattern functions, and avoidance of unpredictable time delays among element. 
   The approach of the present invention makes none of these simplifying assumptions to compute from a given basic element pattern the best possible approximation to a given antenna system radiation function for the desired number of elements in a desired array geometry. The present invention allows the synthesized beam to function in its best capacity allowable by the array physics. 
     FIG. 1  is a flow chart of the steps, in accordance with the present invention, for the construction of an antenna system having a desired far-field radiation pattern F(θ, φ)  10  to be transmitted or received by the antenna system, where θ, φ are the elevation and horizontal/azimuth angles, respective. In the first step  15 , the antenna array parameters are identified. The array element parameters include, at least, the element locations in a reference coordinate system {(x mnk , y mnk , z mnk )} m,n,k , where 0≦m≦M−1, 0≦n≦N−1, 0≦k≦K−1, and MNK is the number of total elements in an array, and a basic element pattern ρ mnk (θ, φ) for each element indexed by m,n,k. 
   In step  20 , the array element parameters are collected into a set of functions and identified as an array frame {A mnk } spanning a radiation function space X (a Hilbert space) in which the desired radiation pattern F(θ, φ) is defined or to be generated. The general expression of {A mnk } will be described. 
   A sequence of vectors/functions {a n } in a Hilbert space H is a frame of H if there exist constants 0&lt;C≦D&lt;∞ such that for all vectors/functions ƒ in H, 
             C   ⁢          f        2       ≤       ∑     n   =   0       N   -   1       ⁢            〈     f   ,     a   n       〉          2       ≤     D   ⁢            f        2     .             
Given a frame {a n }, there exists a dual frame {b n } such that for all f in H, we have
 
   
     
       
         
           f 
           = 
           
             
               
                 ∑ 
                 
                   n 
                   = 
                   0 
                 
                 
                   N 
                   - 
                   1 
                 
               
               ⁢ 
               
                 
                   〈 
                   
                     f 
                     , 
                     
                       b 
                       n 
                     
                   
                   〉 
                 
                 ⁢ 
                 
                   a 
                   n 
                 
               
             
             = 
             
               
                 ∑ 
                 
                   n 
                   = 
                   0 
                 
                 
                   N 
                   - 
                   1 
                 
               
               ⁢ 
               
                 
                   〈 
                   
                     f 
                     , 
                     
                       a 
                       n 
                     
                   
                   〉 
                 
                 ⁢ 
                 
                   
                     b 
                     n 
                   
                   . 
                 
               
             
           
         
       
     
   
   Therefore, in step  25 , a dual frame {B mnk } is determined. With a dual frame known, the system radiation pattern function can be expressed in two ways based on the two forms of the frame expansion in the radiation function space X: 
   
     
       
         
           
             
               
                 
                   ∀ 
                   
                     F 
                     ⁡ 
                     
                       ( 
                       
                         θ 
                         , 
                         ϕ 
                       
                       ) 
                     
                   
                 
                 , 
                 
                   
                     F 
                     ⁡ 
                     
                       ( 
                       
                         θ 
                         , 
                         ϕ 
                       
                       ) 
                     
                   
                   = 
                   
                     
                       ∑ 
                       
                         m 
                         , 
                         n 
                         , 
                         k 
                       
                     
                     ⁢ 
                     
                       
                         〈 
                         
                           
                             F 
                             ⁡ 
                             
                               ( 
                               
                                 θ 
                                 , 
                                 ϕ 
                               
                               ) 
                             
                           
                           , 
                           
                             B 
                             mnk 
                           
                         
                         〉 
                       
                       ⁢ 
                       
                         A 
                         mnk 
                       
                     
                   
                 
               
             
             
               
                 ( 
                 1.1 
                 ) 
               
             
           
           
             
               
                 
                   ∀ 
                   
                     F 
                     ⁡ 
                     
                       ( 
                       
                         θ 
                         , 
                         ϕ 
                       
                       ) 
                     
                   
                 
                 , 
                 
                   
                     F 
                     ⁡ 
                     
                       ( 
                       
                         θ 
                         , 
                         ϕ 
                       
                       ) 
                     
                   
                   = 
                   
                     
                       ∑ 
                       
                         m 
                         , 
                         n 
                         , 
                         k 
                       
                     
                     ⁢ 
                     
                       
                         〈 
                         
                           
                             F 
                             ⁡ 
                             
                               ( 
                               
                                 θ 
                                 , 
                                 ϕ 
                               
                               ) 
                             
                           
                           , 
                           
                             A 
                             mnk 
                           
                         
                         〉 
                       
                       ⁢ 
                       
                         
                           B 
                           mnk 
                         
                         . 
                       
                     
                   
                 
               
             
             
               
                 ( 
                 1.2 
                 ) 
               
             
           
         
       
     
   
   The basic method of determine a dual frame {B mnk } is to put the array frame {A mnk } in a matrix A row-by-row for all each and every indices m,n,k. Note for each given set of m,n,k, A mnk  is a matrix in two angles θ, φ, sampled in appropriate manners. One is to put this matrix A mnk  row-by-row first into a row vector, and then put this vector in the matrix A. A has therefore total MNK rows of vectors. Once A is formed, a dual frame can be calculated by finding the pseudo-inverse of the matrix B. B shall now consist of MNK columns of vectors, the i th  column corresponds to the i th  row of matrix A. 
   Thus in step  30 , with a dual frame {B mnk } determined, the array controlling weight coefficients (simply weights) {w mnk } are computed and used to synthesize the desired radiation pattern F(θ, φ) from the given element described by array frame functions {A mnk }, that is
 
 w   mnk =           ( F (θ, φ),  B   mnk             (1.3)
 
and the array antenna system can thereby be constructed is step  35 .

   In accordance with the present invention, the array weights generating a given radiation pattern F(θ, φ) are generally non-unique (when array element spacing is less than the relative half-wavelength, and/or when the number of elements is greater than the number of sampling points in the array beam pattern F(θ, φ). In such cases, there are infinite many dual frame functions {B mnk }, given by the formula the inventor has derived in a previous research article. The selection of a dual {B mnk } can be made to minimize the cost and energy exciting the array system. 
   In the following description, an array frame construction is described in detail for an arbitrary three dimensional (3D) array system. 
   Array System of Any Geometry 
   As shown in  FIG. 2 , assume that the array elements are placed in a Cartesian System with location given by the coordinates {(x mnk , y mnk , z mnk )} m,n,k , where 0≦m≦M−1, 0≦n≦N−1, 0≦k≦K−1, and MNK is the number of total elements in an array. The placement of these elements need not be on regular grids, nor on a flat plane. The array geometry can be of any shape pending on application. 
   Assume that the element at position (x 000 , y 000 , z 000 ) is our reference element to which all other elements are to refer to determine phase differences among antenna elements for a plane wave in the direction of (θ, φ). Here the plane wave directional parameters θ and φ are as indicated in  FIG. 1 . Since the element spacing are all relevant to the reference element, it is customary to assume that x 000 =0, y 000 =0, z 000 =0. That is, the reference element is assumed to locate at the origin of the Cartesian coordinate system. 
   Assume also that the element basic pattern of the element at the location (x mnk , y mnk , z mnk ) is give by ρ mnk (θ, φ). Then the array frame for the 3D array system is given by 
                           {     A   mnk     }     =       {         ρ   mnk     ⁡     (     θ   ,   ϕ     )       ⁢     ⅇ         2   ⁢   π     λ     ⁢     j   ⁡     (         x   mnk     ⁢   cos   ⁢           ⁢   ϕsin   ⁢           ⁢   θ     +       y   mnk     ⁢   sin   ⁢           ⁢   ϕsin   ⁢           ⁢   θ     +       z   mnk     ⁢   cos   ⁢           ⁢   θ       )             }         m   =   0     ,     n   =   0     ,     k   =   0           M   -   1     ,     N   -   1     ,     K   -   1                     =       {         ρ   mnk     ⁡     (     θ   ,   ϕ     )       ⁢     ⅇ         2   ⁢   π     λ     ⁢     j   ⁡     (         x   mnk     ⁢   u     +       y   mnk     ⁢   v     +       z   mnk     ⁢   cos   ⁢           ⁢   θ       )             }         m   =   0     ,     n   =   0     ,     k   =   0           M   -   1     ,     N   -   1     ,     K   -   1                 ,           (   1.4   )               
where λ is the wavelength of the operating frequency, j is the complex symbol, u=cos φ sin θ, v=sin φ sin θ are the direction cosines with which u 2 +v 2 =sin 2 θ≦1 and cos θ=√{square root over (1−(u 2 +v 2 ))}, −π/2≦θ≦π/2. Note that we have assumed that x 000 =0, y 000 =0, z 000 =0. Otherwise, parameters x mnk , y mnk , z mnk  in formula (1.4) are to be replaced by x mnk −x 000 , y mnk −y 000 , z mnk −z 000 , respectively.
 
   The steps involved to construct a 3D antenna system include placing elements in an desired 3D formation/geometry with any non-uniform spacing (roughly around half-wavelength or smaller) as desired; measuring, modeling or specifying the basic element patterns ρ mnk (θ, φ); measuring the element phase differences based on the cables connected to the elements and their lengths and translating the phase differences into spacing parameters x mnk , y mnk , z mnk ; or measuring the phase differences electronically and then translating the phase differences into spacing parameters x mnk , y mnk , z mnk . Next, the frame operator G is formed, inverted and applied to the array frame to compute the dual frame {B mnk } (or a pseudo-inverse of a matrix formed by array frame functions is carried out as specified before step  30 ) and finally, the array controlling weights w mnk =           F(θ, φ), B mnk            are determined from the dual frame {B mnk } and a desired system radiation pattern F(θ, φ).
   A preferred generating function for desired radiation pattern F(θ, φ) at sampling angles is 
   
     
       
         
           
             F 
             ⁡ 
             
               ( 
               
                 θ 
                 , 
                 ϕ 
               
               ) 
             
           
           = 
           
             
               10 
               
                 exp 
                 ⁡ 
                 
                   ( 
                   
                     
                       
                         θ 
                         8 
                       
                       
                         σ 
                         1 
                         2 
                       
                     
                     - 
                     
                       
                         ϕ 
                         8 
                       
                       
                         σ 
                         2 
                         2 
                       
                     
                   
                   ) 
                 
               
             
             . 
           
         
       
     
   
   The followings are some specific applications. 
   Uniform Planar Array 
   In accordance with the present invention, array parameters are measured. Specifically, the basic element patterns ρ mn0 (θ, φ) are measured and specified in a constructed array running at the operating frequency. Measured patterns take the mutual coupling into account. Next, element phase differences are measured that translates into actual element spacing. Array frame {A mn0 } is then formed, and dual frame {B mn0 } computed. Array control weights {w mn0 } are then determined.  FIG. 3   a  shows the weight matrix {w mn0 } of a uniformly distributed planar array with half-wavelength spacing.  FIG. 3   b  is the corresponding beam patterns in direction cosines. The sidelobe level (SLL) of this particular beam is at about −26 db. Traditional uniformly illuminated (controlled) beam has SLL at about −15 db. 
   Non-Uniform Planar Array 
     FIG. 4   a  is a nonuniformly spaced planar array. The spacing variation is clearly visible. One of the advantage of the present invention is that element spacing needs no longer be made mechanically precise. Array control weights will compensate the spacing variations, together with mutual couplings and other factors. As specified in uniform arrays, array parameters are measured first which includes the actual phase differences between and among elements. Actual spacing information is therefore determined. The determination of a dual frame {B mn0 } will then take spacing variations and mutual coupling into consideration. Array control weights {w mnk }are therefore reflecting such spacing variations. 
     FIGS. 4   b - 4   c  are the real and imaginary (the imaginary is practically zero at the magnitude of 10 −16 ) components of the array control weights, respectively. The giggly behavior of the weights ( FIG. 4   b ) reflects exactly the spacing variation, necessary for producing high beam qualities. No existing construction method can handle such issues precisely. 
     FIG. 4   d  is the resulting beam pattern. The main beam is clearly unchanged. Some slight increase of far-side side-lobes can be detected. Beam quality is clearly high. 
   Array with Element Gain Variations 
   In practical antenna constructions, element characteristics can never be identical as we wished for. In accordance with the present invention, element patterns ρ mn0  (θ, φ) is actually measured. Gain variations is therefore reflected. Showing in  FIG. 5   a  is the array control weights of an application where some element gain is notably different. The difference is reflected in the weights, and the resulting beam pattern  FIG. 5   b  is exactly the same as though all elements are identical (compare with  FIG. 3   b ). 
   Larger Arrays for Enhanced Beam Width/Directivity 
   Showing in  FIGS. 6   a  and  6   b  are the array control weights and corresponding beam patterns of a larger array, where the beam width or directivity is clearly much better. 
   Uniform Cylindrical Array 
   Showing in  FIG. 7   a  is a uniformly distributed cylindrical array with half-wavelength spacing. Array parameters are first measured in accordance with the present invention.  FIGS. 7   b - 7   c  are the real and imaginary array control weights, respectively. The weight determination is precise and highly non-trivial, in accordance to the present invention. No such weight matrix has been seen in literature.  FIG. 7   d  is the corresponding beam pattern. 
   Non-Uniform Cylindrical Array 
   To demonstrate the advantage of the present invention, a non-uniform cylindrical array application is showing in  FIGS. 8   a - 8   d.    FIG. 8   a  shows the non-uniform array distributed on the surface of a cylinder. The spacing variation is random. The present invention creates the control weights as shown in  FIGS. 8   b - 8   c.  The giggly behaviors of the weights, both in real and imaginary components, is exactly necessary to compensate the element spacing variation. The resulting beam pattern  FIG. 8   d  is clearly of superb quality because of the compensation provided. 
   If a cylindrical array (uniform) is illuminated/controlled by equal magnitude weights, as is done traditionally, the beam pattern would be seen in  FIG. 9 , which is clearly of limited quality. 
   Spherical and Truncated Conical Array Systems 
   In accordance with the present invention, array systems built on the surface of a sphere or sectional sphere such as  FIG. 10  and truncated conical surface are constructed similarly by first measuring the array parameters, forming the array frame {A mnk } and the calculating a dual array frame {B mnk }, followed by the weight evaluation. 
   Although the present continuation invention has been described in considerable detail with reference to certain preferred versions thereof, other versions are possible. Therefore, the spirit and scope of the appended claims should not be limited to the description of the preferred versions contained herein.