Title of Invention

"A METHOD OF LOCATING ONE OF MORE TRANSMITTERS ON THE BASIS OF AN ARRAY OF SENSORS"

Abstract

A method of locating one or more transmitters on the basis of an array of sensors moving with respect to the transmitters, comprising at least the following steps for an m-th transmitter of the transmitters: determining a direction vector akm, corresponding to the response of the array of sensors to a source with incidence () as a function of the incidence parameters 9, A, and of an unknown parameter p related to the distortion of the phases on the sensors, transforming the direction vector so as to eliminate the unknown parameter p, into a transformed direction vector, using the transformed direction vector to obtain the position of the transmitter using a maximized locating criterion; wherein the transformation step comprises at least the following steps: choosing a reference instant i, and constructing the transformed direction vector c'km whose components correspond to the ratio of the components between the given instant k and the instant akm is the direction vector at the reference instant i.

Full Text

METHOD FOR TRACKING A TRANSMITTER BY MEANS OF A SYNTHETIC SPARSE ANTENNA NETWORK The invention relates to a method of locating one or more transmitters stationary or mobile on the ground on the basis of the running of a carrier and of an onboard sensor producing the associated direction vectors. It is applied for example with an array of synthetic sparse antennas. The prior art describes various procedures for locating one or more transmitters on the basis of running carriers. Figure 1 illustrates an example of airborne locating. The transmitter 1 to be located is at the position (x0,yo,z0); the carrier 2 at the instant & is at the position (x*,y/(,zk) and sees the transmitter at the incidence (0((it,x0,y0,Z0), J(/tox0,yo,zo)). The angles 0(t,Xo,yo,Zo) and A(t,Xo,yo,Zo) evolve over time and depend on the position of the transmitter as well as the trajectory of the carrier. The angles flftx0,yo,Zo) and A(t,Xo,yo,z0) are for example labeled as shown by figure 2 with respect to an array of N antennas that may be fixed under the carrier. There currently exist several families of locating techniques making it possible to determine the position (xm,ym,zm) of a transmitter on the basis of the direction vectors. These locating techniques differ through the geometric constraints or characteristics of the antennal array; they are classed into several categories cited hereafter. Use of unambiguous, paired, standardized and calibrated arrays and of goniometry techniques These techniques are in most cases based on 1D azimuthal goniometry. The azimuths 6U= 6(tkfxm,ym,zm) associated with the mth transmitter are measured for various instants tk. By using the position (xk,yk,zk) of the carrier at this instant k, a position (x.mk,ymk,Zmk) of this transmitter is estimated through a ground intersection. The position (xk,yk,Zk) of the carrier is given by a GPS. Its orientation is given by a compass in the case of a terrestrial carrier, and by a navigation platform in the case of an aircraft. On the basis of all the positions (xm^,y^,zm^) an extraction of data is performed making it possible to determine the dominant position (xm,ym,zm) of the incident transmitter. Locating is done by triangulation or by ground intersection (2D goniometry). The drawback of triangulation techniques is that they require a significant transit. Moreover, in antennal constraint terms, these goniometry techniques must use an unambiguous array of sensors, and require calibration, pairing and standardization of the channels. Use of sparse array of paired antennas The applicant's patent application FR 03/13128 describes a method which makes it possible to effect direct estimation of the position (xm,ym,zm) of the transmitter on the basis of a multichannel parametric analysis of the direction vectors at various instants ^ over a duration At. This method requires a pairing of the channels as well as a correction of the distortions in phase and in amplitude of the receivers. The present invention relies on a different approach which does not require, during normal operation, correction of the distortions of the receivers provided that the receiver exhibits a substantially constant response over the locating durations considered. The invention relates to a method of locating one or more transmitters on the basis of an array of sensors moving with respect to the transmitters characterized in that it comprises at least the following steps: o determining the direction vectors a^ corresponding to the response of the array of sensors to a source with incidence (9, A) as a function of the incidence parameters 9, A, and of the parameter p related to the distortion of the phases on the sensors, o transforming this vector a/c so as to eliminate the unknown parameter p, into a transformed vector c , c' ion km o using the transformed vector to obtain the position of the transmitter using a maximized locating criterion. The method according to the invention exhibits the following advantages in particular: o It allows direct estimation of the positions of each of the transmitters on the basis of a multichannel parametric analysis at various instants t, o ft makes it possible to use of arrays sparse sensor that are unpaired, or even uncalibrated (large-aperture array), o It makes it possible to circumvent the pairing as well as the calibration in amplitude and in phase of the reception channels, o It is possible to take into account a model on the variation in the signal reception level, o it makes it possible to introduce any observation model deduced from the direction vectors at different observation instants, o it is possible according to a variant to consider the whole set of channels, without taking a particular reference channel. Other characteristics and advantages of the invention will be better apparent on reading the description which follows of an example given by way of wholly nonlimiting illustration with appended figures which represent: o figure 1 an exemplary locating of a transmitter by an airplane equipped with an array of suitable sensors, o figure 2 an array of 5 antennas and the angles of incidence of a transmitter. In order to better understand the principle implemented in the present invention, the example which follows is given by way of wholly nonlimiting illustration for a system such as shown diagrammatically in figure 1, comprising an airplane 2 equipped with an array of N sensors and with a processor adapted for executing the steps of the method according to the invention and with a transmitter 1 on the ground to be located. In the presence of M transmitters, the airplane receives at the instant t at the output of the N sensors of the array, the vector x(f). Around the instant tk, the vector x(t+tk) of dimension Wx1 is the mixture of the signals of the M transmitters: M M m=l for \t\ o where b(f) is the noise vector assumed Gaussian, o a(6, A, p) is the response of the array of sensors to a source with incidence (6,A) with complex receiver gains, o pis a parameter representative of the phase distortion on the various reception channels, p=[ pi ... pN]T, o A™=[ a(6«.,, A^p)...,a(QKM, A^p)], a(r)=[cr,(r)...oxX*")]1, 6^^6(7^,^,^) and Axm=J(f/(,xm,ym,zm) and xn(t) is the received signal on the nth sensor. We note that in this model, the mixture matrix A*p depends on the instant tk of observation as well as on the gains of the receivers p. The direction vector a(6,A,p) and the direction vector without phase distortion v(0, A) respectively have expressions: a(6, A, p)= M and v(0, A)= M (2) Where: • bn(0 is the noise vector of channel n, • a(6, A, p) is the response of the array of sensors to a source with incidence (0,A), • v(6, A) is direction vector without distortion for a source with incidence (6,A) • A fc(p)=[a(6/(1, AM, p)... a(efcM, Am p) ], • 6k=6»(^1xm,ym,zm) and A/c=/l(^xm,ym,zm). These vectors a(6fcm, A/ Direction vector akm In the presence of model errors, for example coupling, the measured direction vector akm(p) can be written: a/cm(p) = a(^^,xm,ym,zm), A(ffcxm,ym,zm), p) + e^m with (3) - complex measurement noise dependent on the calibration errors such as coupling. In particular for an array composed of A/=2 sensors spaced apart by a distance of d along the axis of the carrier the direction vector without distortion at the instant k satisfies: 1 ( d ( \ i \exp J2n-~cos(0(tk,xm,ym,zj)cos(&(tk,xm,ym,zm)) \ A ) (4) The incidence (d(tk,xm,ym,zm), A(tk,xm,ym,zm)) can be calculated directly on the basis of the position (xm.ym.Zm) of the carrier at the instant tk and of the position (xm,ym,zm) of the transmitter. The basic idea of the method relies notably on the fact that the parameter p is independent of the incidence (0, A) of the sources, the vector p not being estimated by standardization or by any other procedure before or during normal operation of the system. First variant of implementation of the method using vectors According to a first variant embodiment, the method uses the following direction vector: , A*, (5) It constructs the vector b^ of dimension (N-1)*1 by choosing a reference sensor, for example the sensor of index 1. M )/»*.(!) M where &i The vector b, is a function of p and is expressed in the following manner: bkm = (6) where b(f/c,xm,ym,zm,p)= M The vector bkm is constructed on the basis of the components of the vector a/tf77(p). We choose for this purpose a reference channel and we construct the components of the vector as being the ratio of the components of the vector a/cm and a reference channel which is associated for example with the 1st component of Equation (6) corresponds to the choice of sensor 1 as reference sensor. The vector b^ is thereafter transformed so as to eliminate the components of the complex vector pn. For this purpose, the method constructs the vector £ of dimension (A/-1)x1 by choosing a reference instant, for example k=i: M (7) M Thus in a low noise context || w^ ||«1, C^ (ri) has the following expression: (n) (n) where the non calibrated disturbed vector 8) with: wkm(n)= " " ' We note that the ratio vn(6'(t,Ai,p)/vn(6',,A/,p)is independent of p and equals: A i (9) Under these conditions the vector c^ may be written in the following manner: (10) where c(ffcxm,ym,zm)= M th and wtw (n) is the n component of w^ Locating the transmitter The method having determined the vector c^, it uses it to locate the transmitters, that is to say to obtain the position (xm,ym,zm) of the transmitter. For this purpose the method maximizes the normalized vector correlation criterion l_K(x,y,z) given by expression (11) in the position space (x,y,z) of a transmitter. with -1m M "Km = vK,c(xm,ym,zm) + wK, vK,c(x,y,z) = c(t},x,y,z) M c(tK,x,y,z) w 1m and M (11) W Km In order to refine the estimation of the position (xm,ym,zm) of the transmitters, the method can be implemented in an iterative manner. For this purpose, the method executes together for example the following steps: Step 11 Step 12 at the instant Identification of the vectors Calculation of the criterion LK+I(X,V,Z) in an iterative manner and minimization of U+1(x,y,z) to obtain a new position estimation (Xm,ym,zm) for the transmitter and possible loopback to the preceding step. The noise vector WK has covariance matrix R=E[w/ Under these conditions we obtain the following criterion l_K'(x,y,z): LK'(x,y,z)= with R=E[wKwK ] Other variant of implementation of the method According to another variant embodiment, the method constructs a vector C^of dimension A/x1 on the basis of the direction vector *km- For this purpose, the method applies a transformation which consists, for example, in choosing a reference instant k=i and in forming the ratios of the components of the vector &km with the component corresponding to the instant i. The vector c'te is then expressed in the following manner: km M »*»»/»*,(") M (13) This corresponds to a temporal tracking of the measurements of the vector Sk Thus in a low noise context || ekm ll (14) with: Under these conditions the vector C may be written in the following manner: km V /v iih km 1(15) where c\tk M and wto (n) is the n component of v/ km In the case of large running between the measurements, the method can determine vectors ckm or c\m constructed over a sliding time window (rather than over a constant window defined by a reference instant as is the case in the first variant) by taking /=/f-1 or i=k-L, where L corresponds to the length of the window (number of measurement samples considered) The measurement of the direction vectors &*m is generally obtained to within an undetermined complex factor. For the variant of the method using the vectors c'^, the method can comprise a step which consists in changing the phase reference of the direction vector measured by choosing a virtual channel as reference (and not a real channel as is the case in the first variant) defined, for example, by the phase ban/center (defined to within a constant scalar coefficient that may arbitrarily be fixed at 1). This operation is carried out, for example, by applying to the measured vectors a^ the following transformation: *t 1 T \~krn V / I £ /4 ca*»= lllTM *« ( } V ' "fenvMy The correction coefficient is not fully determined by this expression having regard to the indeterminacy of order N at each instant k of the complex root. A tracking of the phase evolution during the observation period makes it possible to resolve the indeterminacy, as is described hereafter. The complex correction coefficient being defined to within a factor from among the N Nth roots of unity, the phase tracking consists in arbitrarily fixing the first (k=1) correction coefficient (by taking root 1 for example), then in determining at each new iteration k+1, the correction coefficient p, from among the Nth roots of unity, which minimizes the mean phase deviations between the direction vector &(^i)m measured at the instant k+1 and the corrected vector at the instant k a"*m . The minimization criterion, for measurements at the same frequency, can be defined by the following expression: ™re£, ^ ^ v a^(l) j (17) For measurements at different frequencies, it is possible to compare the phases of the components of the two direction vectors by correcting them for a power given by the ratio of these two frequencies. If we consider the vectors C?km, it is then possible to compare them with the theoretical values c'(fycX/TJ.ym,zm) for which the theoretical direction vector a(f/oxm.yAr7.zm) 's calculated with reference to the virtual channel defined by the phase ban/center (theoretical geometric phase barycenter which is the geometric locus for which the theoretical sum of the phase differences vanishes). This locus does not coincide, in general, with the phase center of the array (determined experimentally). Locating the transmitter on the basis of the C! km The locating method thereafter comprises a step which consists in maximizing the following normalized vector correlation criterion L/ i H l_K(x,y,z)= with C'K - • \m M ' (x,y,z) = M (18) Km w 1m and w/ M 1 Km To obtain the position (xm,ym,zm) of the transmitter, the method calculates and maximizes the criterion L/^x.y.z) of equation (18). In order to refine the estimation of the position (xm,ym,ym) of the transmitters the steps of the method can be conducted in an iterative manner, for example in the following manner: o Step 13 Identification of the vectors c'(K+i)m at the instant fK+i- o Step 14 Calculation of the criterion LK+i(x,y,z) in an iterative manner and minimization of L/ The noise vector w/ Under these conditions we obtain the following criterion L-KXx.y.z): (19) LK1(x,y,z)= with R=E[WK WK] Filtering step In the case of a very large number of measurements, the method can for example comprise a step of prior processing of the vectors a* at the K instants k- This processing, executed before the steps for determining the coefficients CIM , makes it possible to reduce the numerical complexity of calculation (which is dependent on the number of measurements) by decreasing K. By way of example, it is possible to perform on the elementary measurements the following processings: - decimation of the instants tk, - filtering (smoothing of the measurements &km) and under sampling, merging of the measurements over a defined duration (extraction by association of direction vector and production of synthesis measurements). In the variant of the method using the vectors c1 lm, it is possible to introduce a model on the variation in level of the received signal based, for example, on a free-space propagation model. The method in its general form, by recourse to the criteria of vector correlation between measurements and models (L/^x.y.z) and L.K'(x,y,z)), can take into account other observation models derived from the direction vectors a*. In the presence of an array of paired antennas, it is notably possible to define an observation model combining the vectors bk and CK. or c' km. According to a variant embodiment, it is possible to use the filtering method described in patent application FR 03/13128. If the decimation ratio is significant (non-negligible transit between two filtered measurements), it is also possible to supplement the decimation with an estimation of the variation in direction vector, and more particularly the rate of variation of the differential phase for each of the channels, and to employ a distance-based model utilizing this variation. Specifically, the phase of each component bk,n of the relative vector b/c is dependent on the incidence (9« ,kk) and on the angle a n formed by the axis defined by the positions of the aerials of the channels n and 1 with the trajectory of the carrier. If y is the angle formed by the direction of incidence and the trajectory of the carrier (pseudo bearing), we have: - bk,n =A,ej^(r) with: Avn(t)=2p-cos(y(t)+an) The instantaneous variation of differential phase may then be written: ~.Bi from which we derive a relation on the distance of the transmitter: (20) This relation then makes it possible, as an adjunct to the vectors b^ or c/f, to employ the estimations of instantaneous variation of differential phase by introducing a distance model into the correlation criterion U(x,y,z). In the case where y =0 (aerials along the axis of the carrier), the relations may be written: WE CLAIM: 1. A method of locating one or more transmitters on the basis of an array of sensors moving with respect to the transmitters, comprising at least the following steps for an m- th transmitter of the transmitters: determining a direction vector akm corresponding to the response of the array of sensors to a source with incidence () as a function of the incidence parameters 0, A, and of an unknown parameter p related to the distortion of the phases on the sensors, transforming the direction vector so as to eliminate the unknown parameter p, into a transformed direction vector, using the transformed direction vector to obtain the position of the transmitter using a maximized locating criterion; wherein the transformation step comprises at least the following steps: choosing a reference instant i, and constructing the transformed direction vector c'km whose components correspond to the ratio of the components between the given instant k and the instant i (Formula Removed) aim is the direction vector at the reference instant i. 2. The method as claimed in claim 1, wherein the locating criterion is equal to: (Formula Removed) noise vector for all the components. 3. The method as claimed in claim 1, further comprising the following iterative steps: identifying vectors c'(k+1)m at an instant tk+1; calculating a criterion Lk+1 (x,y,z) in an iterative manner and minimizing the criterion Lk+1 (x,y,z) to obtain a new position estimation (xm,ym,zm) for the transmitter. 4. The method as claimed in claim 1, wherein a phase tracking is performed for the vectors c'km. 5. The method as claimed in claim 1, further comprising a step of filtering the direction vectors ak. 6. The method as claimed in claim 1 wherein a sliding time window is used by taking i=k-l or i=k-L. 7. The method as claimed in claim 2, further comprising the following iterative steps: identifying vectors c'(k+y1)m at an instant tk+1; calculating a criterion Lk+1(x,y,z) in an iterative manner and minimizing the criterion LK+1(x,y,z) to obtain a new position estimation (xm,ym,zm) for the transmitter. 8. A method of locating one of more transmitters on the basis of an array of sensors moving with respect to the transmitters which comprises at least the following steps: receiving signals from at least one of the transmitters by the sensors and determining a direction vector akm; choosing a reference channel or sensor and constructing a vector b(tK,Xm,Zm,p) whose components correspond to the ratio of 2 channels according to the received signals, where m denotes the m-th transmitter of the transmitters; (Formula Removed) choosing a reference instant and constructing on the basis of the vector b(tk,xm,ym,) a vector whose components are: (Formula Removed) Where akm(n) is the nth component of akm constructing a transformed direction vector ckm (Formula Removed) using the transformed direction vector to obtain the position of the transmitters using a maximized locating criterion; wherein m is the amount of transmitters, and N is the amount of sensors. 9. A method of locating one of more transmitters on the basis of an array of sensors moving with respect to the transmitters, comprising at least the following steps for an m-th transmitter of the transmitters: receiving signals from at least on of the transmitters by the sensors; using the received signals to construct a direction vector; (Formula Removed) where akm(n) is the nth component of akm wherein the vector is a function of p and is expressed in the following manner: (Formula Removed) and wherein the vector is constructed on the basis of the components of the vector akm(); transforming the vector so as to eliminate the components of the complex vector n; constructing a vector ckm of dimension (N-1)xl by choosing a reference instant wherein k=i: (Formula Removed) using the vector ckm to obtain the position of the transmitter using a maximized locating criterion. 10. The method as claimed in claim 9, wherein the locating criterion comprises a normalized vector correlation criterion Lk (x, y, z) in the position space (x, y, z) of a transmitter with: (Formula Removed) 11. The method as claimed in claim 10, further comprising the following iterative steps: identifying the vectors C(K+1)m at an instant tK+1; calculating a criterion LK+i (x,y,z) in an iterative manner and minimizing the criterion LK+1 (x,y,z) to obtain a new position estimation (xm,ym,zm) for the transmitter. 12. The method as claimed in claim 10, further comprising the following iterative steps: identifying the vectors C(K+1)m at an instant tK+1; calculating a criterion LK+i (x,y,z) in an iterative manner and minimizing the criterion LK+1 (x,y,z) to obtain a new position estimation (xm,ym,zm) for the transmitter.

Documents:

6739-DELNP-2006-Abstract-(02-11-2011).pdf

6739-delnp-2006-abstract.pdf

6739-DELNP-2006-Claims-(02-11-2011).pdf

6739-delnp-2006-Claims-(15-04-2014).pdf

6739-delnp-2006-claims.pdf

6739-DELNP-2006-Correspondence Others-(02-11-2011).pdf

6739-delnp-2006-Correspondence Others-(15-04-2014).pdf

6739-delnp-2006-Correspondence Others-(20-02-2014).pdf

6739-delnp-2006-correspondence-others-1.pdf

6739-delnp-2006-correspondence-others.pdf

6739-delnp-2006-description (complete).pdf

6739-DELNP-2006-Drawings-(02-11-2011).pdf

6739-delnp-2006-drawings.pdf

6739-DELNP-2006-Form-1-(02-11-2011).pdf

6739-delnp-2006-form-1.pdf

6739-delnp-2006-form-18.pdf

6739-DELNP-2006-Form-2-(02-11-2011).pdf

6739-delnp-2006-form-2.pdf

6739-DELNP-2006-Form-3-(02-11-2011).pdf

6739-delnp-2006-form-3.pdf

6739-delnp-2006-form-5.pdf

6739-DELNP-2006-GPA-(02-11-2011).pdf

6739-delnp-2006-GPA-(20-02-2014).pdf

6739-delnp-2006-gpa.pdf

6739-DELNP-2006-Petition-137-(02-11-2011).pdf

abstract.jpg