Method for estimating the precise orientation of a satellite-borne phased array antenna and bearing of a remote receiver

Abstract: The precise three-axis attitude of a space-borne phased-array antenna is estimated based on the assumption that the array geometry, consisting of the number of radiating elements and their relative spacing in three dimensions, is known and that the array position and coarse knowledge of the array attitude are available a priori. An estimate is first made (21 1 , 21 2 ... 21 M ) that define each element's straight-through contribution to the signals received at each of two or more remote calibration sites, where a "straight-through" antenna configuration is defined as the condition in which all elements are made to radiate with the same amplitude and phase. An optimization strategy is then used (22,23) to determine which array attitude lying in the neighbourhood of the coarsely known attitude is most consistent with the full set of straight-through gain values. Another technique for estimating the precise angular location of a receiver with respect to the coordinates of the space-borne phased-array antenna is based on the assumptions that the array geometry is known, and that the receiver bearing is coarsely known or available. After an estimate is made of the set of complex-valued gains that define each element's straight-through contribution to a composite signal measured at the receiver site, an optimization strategy is used to determine which receiver direction lying in the neighbourhood of the coarsely known direction is most consistent with the latter set of straight-through gain values.

Description: This invention relates to satellite communications and, more particularly, to a method for estimating the precise three-axis attitude of a space-borne phased-array antenna and the precise angular location of a receiver with respect to the coordinates of the space-borne phased-array antenna. Precise attitude knowledge of the orientation of a satellite-borne phased-array antenna is critical when the antenna pattern is highly directed, especially if the satellite serves multiple ground-based transmitter/receiver sites with a high degree of geographic selectivity. Attitude control systems employed on current state-of-the-art commercial communication satellites are capable of sensing and maintaining attitude to within approximately 0.1° in each of three rotational coordinates. For a satellite orbiting the earth at geosynchronous altitude, this corresponds to an uncertainty of approximately 60 km on the ground. However, the orientation of a space-borne phased-array antenna needs to be measured with significantly greater precision than the levels just cited for the next generation of geostationary communication satellites. In addition, calibration of a satellite-borne phased-array antenna from the ground (or from any remote site) requires precise knowledge of the bearing of the calibration site with respect to the radiation pattern of the array. This is because one needs to distinguish the effects of attitude disturbances from drifts in the phasing circuits of the array elements, both of which are observed as phase shifts at the receiver. Station-keeping maneuvers employed on current state-of-the art commercial communication satellites maintain positional stability to within approximately 75 km. For geostationary satellites, this implies that fixed locations on the earth's surface have a directional uncertainty of approximately 0.1 to 0.2° with respect to a coordinate system local to both the satellite and the array. This level of uncertainty significantly limits the precision with which the array can be calibrated. As a case in point, the phase shifters located at the comers of a 16x16 array with a three wavelength element spacing can drift up to approximately 0.04 cycles in phase before the effect seen at a receiver on the ground begins to exceed that of attitude and position uncertainty. This implies that the maximum phase resolution achievable through ground-based calibration is between four and five bits. Phased-array payloads being designed for deployment in the next generation of geostationary communication satellites will employ up to 256 levels (i.e., eight bits or 2 8 ) of phase resolution. To calibrate such systems from the ground will require at least an order of magnitude improvement either in position and attitude sensing capability or in other means for ascertaining the precise angular coordinates of the calibration site. It is therefore an object of the present invention to provide a computer implemented method for estimating the precise orientation of a satellite-borne phased-array antenna during calibration of the array from two more remote sites. It is another object of the invention to provide a computer implemented method for estimating the precise bearing of a remote receiver with respect to the radiation coverage of a satellite-borne phased-array antenna. According to one aspect of the invention, a computer implemented technique is provided for estimating the precise three-axis attitude of a space-borne phased-array antenna. The technique assumes that the array geometry, consisting of the number of radiating elements and their relative spacing in three dimensions, is known, and that the array position and coarse knowledge of the array attitude are available a priori . A hypothetical "straight-through" antenna configuration is defined as the condition in which all elements are made to radiate with the same amplitude and phase. The technique according to this aspect of the invention consists of two steps. First, an estimate is made of the set of complex-valued gains that define each element's straight-through contribution to the signals received at each of two or more remote calibration sites. Second, a determination is made by means of a mathematical optimization strategy as to which array attitude lying in the neighborhood of the coarsely known attitude is most consistent with the full set of straight-through gain values determined in the first step. According to another aspect of the invention, a computer implemented technique is provided for estimating the precise angular location of a receiver with respect to the coordinates of a space-borne phased-array antenna. This technique is based not on any assumption that the array position and attitude are known or available, but instead on the assumptions that the array geometry is known, as in the first-described technique, and that the receiver bearing is coarsely known or available. This technique, like the first-described technique, consists of two steps. First, an estimate is made of the set of complex-valued gains that define each element's straight-through contribution to a composite signal measured at the receiver site. Second, a determination is made by means of a mathematical optimization strategy as to which receiver direction lying in the neighborhood of the coarsely known direction is most consistent with the straight-through gain values determined in the first step. The features of the invention believed to be novel are set forth in the appended claims. The invention, however, together with further objects and advantages thereof, may best be understood by reference to the following description taken in conjunction with the accompanying drawings, in which: Figure 1 is a pictorial diagram illustrating a satellite-borne phased-array antenna and a plurality of remote ground-based receivers; Figure 2 is a block diagram illustrating the flow of the satellite-bome phased-array attitude estimation technique according to one aspect of the invention; and Figure 3 is a block diagram illustrating the flow of the receiver bearing estimation technique according to a second aspect of the invention. Figure 1 illustrates a satellite-borne phased-array antenna 10 made up of a plurality of radiating elements, and a plurality of remote ground-based receivers 11 and 12, here referred to as Receiver #1 and Receiver #2, respectively. Orientation of space-borne phased-array antenna 10 according to a first aspect of the invention requires use of two or more earth-based receivers 11 and 12 whose precise geographical coordinates are known. The technique itself is a two-step procedure which is schematically represented in the block diagram of Figure 2, to which reference is now made. The first step requires measurement at each receiver site of the so-called "straight-through" signal path gains, as generally indicated at function blocks 21 1 to 21 M . These straight-through gains, which are complex-valued, represent the magnitude and phase that a unit signal attains as it flows through the amplifier chain and propagation path associated with each element in an unsteered array. An unsteered array is defined as one whose elements are made to radiate with a uniform amplitude and phase, represented by a single complex gain value k . In the description that follows, it is assumed that the receiver lies within a region over which the array elements radiate isotropically and that the propagation path is free of atmospheric disturbances. Let G m / n denote the gains measured at receiver site m , where m =1,2,..., M , and M is the number of receiver sites used in the procedure. As seen from the m th receiver site, the straight-through gain for the n th element is given by G m n = k e j 2π λ R m - r n R m - r n , where R m is the receiver position, r n are the element positions expressed in the local coordinate frame, and λ is wavelength. In the far field, i.e., where | r n | 2 << λ| R m |, G m / n can be rewritten as G m n = k e j 2π λ R m R m e - j 2π λ ( û m - r n ) = κ m e - j 2π λ ( û m - r n ) , where R m = | R m |, and û m is a unit vector directed toward the receiver from the local origin. In a steered array, the total gain imposed by each element is the product of G m / n and a selectable gain A n , which, in combination, fully characterize the signal response of the array at the given receiver site. The attitude estimation method described here makes use of the straight-through gains G m / n measured at two or more receiver sites, but requires no knowledge of the selected gains A n . Any method deemed suitable for measuring these straight-through gains can be successfully used in the attitude estimation procedure. One such procedure encodes coherent signals from the phased array elements using controlled switching of the gain and phase shifter delay circuits. Such procedure is set forth in Silverstein et al., U.S. patent 5,572,219, issued November 5, 1996. For N elements, the control circuit switching is dictated by matrix elements of an N x N Hadamard matrix. The encoded signal vectors are decoded with the inverse of the same Hadamard matrix used in the control circuit encoding. Other methods can be used in the attitude estimation procedure, and the invention is not dependent on the particular method used. To implement the second step in the attitude estimation procedure, a model is constructed for the full set of straight-through gains: In this expression, α m is a site-dependent, unknown complex amplitude, and Θ represents a set of angles that define the attitude of the array. As the array position and all receiver positions are assumed known, the array attitude determines all receiver directions û m . It is convenient to think of Θ as consisting of three orthogonal component angles which specify the rotation that the nominal known attitude must undergo to give the true array attitude. The attitude estimation problem thus reduces to finding that set of rotational angles (i.e., roll, pitch and yaw) and complex amplitudes α m for which G m / n best "matches" G m / n . To do this, the measurement vectors g m = [ G m / 1,..., G m / N ]' and signal model vectors e m (Θ) = [ 1 m (Θ), 2 m (Θ),..., Nm (Θ)]' are first defined, where N is the total number of elements and (') denotes the matrix transpose operation. Therefore, g m = α m e m ( Θ ) + n m , where n m is a complex random vector of noise values representing the errors in the measurements G m / n . Next, vectors g, a , and n and matrix E are constructed as follows: Therefore, g = E ( Θ ) a + n . It is assumed that the components of n are zero-mean complex Gaussian variables with E { Re ( n ) Im ( n H )}=0 and E { Re ( n ) Re ( n H )} = E { Im ( n ) Im ( n H )}, where the H denotes Hermitian transpose and E () denotes the expectation operation. A further definition is Σ=E { nn H }. With these definitions in place, it is then possible to write an expression that specifies the maximum likelihood (ML) solution to the attitude estimation problem. Denoting by the corresponding ML estimates of ( a ,Θ), then with F (Θ) defined as F (Θ)= g H Σ -1 E ( E H Σ -1 E ) -1 E H Σ -1 g , where the explicit dependence of E on Θ has been suppressed for clarity of notation. The amplitude estimate, though not explicitly required for attitude estimation, is given by â = ( E H Σ -1 E) -1 E H Σ -1 g, The expressions above simplify greatly for the degenerate case in which the measurement errors are identically distributed; i.e., where Σ=σ 2 I . In this case, the ML estimate for the angle vector specifying the array attitude is given by where ∥ν∥ 2 = ν H ν. The corresponding amplitude estimate is In the process illustrated in Figure 2, the gains G m / n are fit to the model by evaluating F (Θ) and choosing Θ that maximizes F , as indicated at step 22. Maximization of the function F (Θ) can be carried out efficiently in practice by making use of any standard gradient search method 23. As shown in Figure 2, the search begins at Θ = (0,0,0), which implies no rotation at all, and thus represents the initial coarse knowledge of the array attitude. The solution obtained in this manner will be unique if the initial attitude uncertainty is commensurate with the level noted earlier. Simulations based on a hypothetical 16x16 array in a geostationary position above a pair of receiver sites displaced ±3° from the boresight axis of the array demonstrate that approximately 0.001 to 0.01° of attitude precision can be obtained with the method just described. The experiments assume operation at 12 GHz with an element spacing of three wavelengths and a receiver signal-to-noise ratio (SNR) of 20 dB. This represents an improvement of one to two orders of magnitude with respect to the initial three-axis attitude uncertainty of 0.1°. The method for estimating the precise bearing of a remote receiver with respect to the radiation coverage of a satellite-borne phased-array antenna 10 (as shown in Figure 1) is a similar two-step process. As shown in Figure 3, the first step 31 of this process requires measurement of the so-called "straight-through" signal path gains, as above. The straight-through gain for the nth array element, as seen from the receiver, is given by G n = k e j 2π λ R - r n R - r n , where R is the receiver position, r n are the element positions expressed in the local coordinate frame, λ is wavelength, and k again represents the magnitude and phase of the radiation from the array in its "unsteered" state. In the far field, i.e., where | r n | 2 << λ| R |, G n can be rewritten as G n = k e j 2π λ R R e - j 2π λ ( û - r n ) = κ e - j 2π λ ( û - r n ) , where is another complex constant, R = | R |, and û is a unit vector directed toward the receiver from the local origin. In a steered array, the total gain imposed by each element is the product of G n and a selectable gain A n , the values of which are chosen to achieve a desired antenna beam orientation and shape. The two quantities, G n and A n , fully characterize the signal response of the array. However, only the straight-through gains G n are required for implementing the method according to this aspect of the invention, namely, estimation of the receiver bearing û . Any method deemed suitable for measuring these straight-through gains can be successfully used in the bearing estimation procedure. The second step in the bearing estimation procedure is to construct a model for the straight-through gains, as follows: G n (α,  1 , 2 ) = α e - j 2π λ [ û ( 1 , 2 )· r n ] = α  n ( 1 , 2 ). In this expression, α is an unknown complex amplitude, and  1 and  2 are angles that define the receiver direction û . The bearing estimation problem then reduces to finding that set of angles ( 1 ,  2 ), along with the corresponding α for which G n best "matches" G n . This is done by defining a measurement vector g =[ G 1 , G 2 ,..., G N ]' and a signal model vector e ( 1 , 2 )=[ 1 ( 1 ,  2 ),  2 ( 1 ,  2 ),...,  N ( 1 ,  2 )]', where N is the total number of elements and (') denotes the matrix transpose operation. Therefore g = α e (1, 2) + n , where n is a complex random vector of noise values representing the errors in the measurements G n . By assuming that the components of n are zero-mean complex Gaussian variables with E { Re ( n ) Im ( n H )=0 and E { Re ( n ) Re ( n H )] =E { Im ( n ) lm ( n H )}, where the H denotes Hermitian transpose and E () denotes the expectation operation, and by defining Σ= E { nn H }, it is then possible to write an expression that specifies the maximum likelihood (ML) solution to the bearing estimation problem. Denoting by the corresponding ML estimates of (α, 1 , 2 ), then with F ( 1 ,  2 ) defined as F (  1 ,  2 ) = g H Σ -1 e ( e H Σ -1 e ) -1 e H Σ -1 g , where the explicit dependence of e on ( 1 ,  2 ) has been suppressed for clarity of notation. The amplitude estimate, though not explicitly required for bearing estimation, is given by As before, the expressions above simplify greatly for the degenerate case in which the measurement errors are identically distributed; i.e., where Σ=σ 2 1 . In this case, the ML estimates for the angles specifying the receiver direction are given by and the corresponding amplitude estimate is α = 1 N e H g . Maximization of the function F ( 1 ,  2 ) at step 32 of Figure 3 can be carried out efficiently in practice by making use of any standard gradient search method, as indicated at step 33. As shown in Figure 3, the search begins at the values for ( 1 ,  2 ) that correspond to the initial coarse knowledge of the receiver direction with respect to the array. The solution obtained in this manner will be unique if the initial direction uncertainty is commensurate with the level noted above. Simulations based on a hypothetical 16x16 array in a geostationary position above a receiver site displaced 5° from the boresight axis of the array demonstrate that approximately 0.001 to 0.004° of directional precision can be obtained with the method just described. The experiments assume operation at a frequency of 12 GHz with an element spacing of three wavelengths and a receiver signal-to-noise ratio (SNR) of 20 dB. This represents an improvement of one to two orders of magnitude with respect to the initial uncertainty of 0.1 to 0.2°.

Claims: A method for estimating in a computer the precise three-axis attitude of a space-borne phased-array antenna made up of a plurality of radiating elements, comprising the steps of: inputting to the computer the array geometry, including the number of radiating elements and their relative spacing in three dimensions, and the array position and coarse knowledge of the array attitude; defining a hypothetical "straight-through" antenna configuration as a condition in which all of the radiating elements are made to radiate with the same amplitude and phase; estimating in the computer a set of complex-valued gains that define a straight-through contribution by each of the radiating elements to the signals received at each of two or more remote receiver calibration sites; and employing an optimization strategy in the computer to determine which array attitude lying in the neighborhood of the coarsely known attitude is most consistent with the set of straight-through gain values determined in the estimating step. The method for estimating in a computer the precise three-axis attitude of a space-borne phased-array antenna of claim 1 wherein the step of estimating in the computer a set of complex-valued gains comprises the steps of: measuring at each of said two or more remote receiver calibration sites straight-through signal path gains; and constructing a model for a full set of straight-through gains based on the measured straight-through signal path gains. The method for estimating in a computer the precise three-axis attitude of a space-borne phased-array antenna of claim 2 wherein G m / n denotes the gains measured at a receiver calibration site m , where m =1,2,..., M , and M is the number of receiver sites and, as seen from the m th receiver site, the straight-through gain for the n th element of the phased-array antenna is given by G m n = k e j 2π λ R m - r n R m - r n , where R m is the receiver position, r m are the element positions expressed in a local coordinate frame, and λ is wavelength, and in the far field where | r n | 2 << λ| R m |, G m n = k e j 2π λ R R m e - j 2π λ ( û m - r n ) = κ m e - j 2π λ ( û m - r n ) , where R m = | R m | , and û m is a unit vector directed toward the receiver calibration site from the local origin, and wherein the model constructed for the full set of straight-through gains is expressed as where α m is a site-dependent, unknown complex amplitude, and Θ represents a set of angles that define the attitude of the array, and wherein the step of employing an optimization strategy in the computer to detemine which array attitude lying in the neighborhood of the coarsely known attitude is most consistent with the set of straight-through gain values comprises finding a set of rotational angles Θ and complex amplitudes α m for which G m / n best matches G m / n . A method for estimating in a computer the precise angular location of a receiver with respect to the coordinates of a space-borne phased-array antenna made up of a plurality of radiating elements, comprising the steps of: inputting to the computer the array geometry, including the number of radiating elements and their relative spacing in three dimensions, and coarse knowledge of the receiver bearing; defining a hypothetical "straight-through" antenna configuration as a condition in which all of the radiating elements are made to radiate with the same amplitude and phase; estimating in the computer a set of complex-valued gains that define a straight-through contribution by each of the radiating elements to a composite signal measured at the receiver site; and employing an optimization strategy in the computer to determine which receiver direction lying in the neighborhood of the coarsely known bearing is most consistent with the set of straight-through gain values determined in the estimating step. The method for estimating in a computer the precise angular location of a receiver with respect to the coordinates of a space-borne phased-array antenna of claim 4 wherein the step of estimating in the computer a set of complex-valued gains comprises the steps of: measuring at said remote receiver site straight-through signal path gains; and constructing a computer model for a full set of straight-through gains based on the measured straight-through signal path gains. The method for estimating in a computer the precise angular location of a receiver with respect to the coordinates of a space-borne phased-array antenna of claim 5 wherein G n denotes the straight-through gain for the n th array element as seen from the receiver, and is given by G n = k e j 2π λ R - r n R - r n , where R is the receiver position, r n are the element positions expressed in a local coordinate frame, λ is wavelength and k represents the magnitude and phase of the radiation from the array in an unsteered state and, in the far field where r n 2 << λ| R , G n = k e j 2π λ R R e - j 2π λ ( û - r n ) = κ e - j 2π λ ( û - r n ) , where R =| R | , and û m is a unit vector directed toward the receiver from the local origin, and wherein the model constructed for the set of straight-through gains is expressed as where α is an unknown complex amplitude and  1 and  2 are angles that define the receiver direction û , and wherein the steps of employing an optimization strategy in the computer to determine which receiver direction lying in the neighborhood of the coarsely known bearing is most consistent with the set of straight-through gain values determined in the estimating step comprises finding a set of angles ( 1 ,  2 ), along with the corresponding α for which G n best matches G n .

Inventors: ASHE JEFFREY MICHAEL, BONANNI PIERINO GIANNI, SILVERSTEIN SETH DAVID, ASHE, JEFFREY MICHAEL, BONANNI, PIERINO GIANNI, SILVERSTEIN, SETH DAVID
Applicants: GEN ELECTRIC, GENERAL ELECTRIC COMPANY
Application year: 1997
Main classification: H01Q1
IPC codes: H01Q21, H01Q3, H01Q1