Login| Sign Up| Help| Contact|

Patent Searching and Data


Title:
APPARATUS FOR ESTIMATING A DIRECTION OF ARRIVAL AND CORRESPONDING METHOD
Document Type and Number:
WIPO Patent Application WO/2018/054467
Kind Code:
A1
Abstract:
The invention refers to an apparatus (1) for estimating a direction of arrival. The apparatus (1) comprises an antenna (2), a beamforming network (3), and an e valuator (4). The antenna (2) is configured to receive signals, is circularly polarized, and comprises a plurality of different radiation patterns. The beamforming network (3) is configured to provide based on signals received by the antenna (2) decomposed signals that are received by associated radiation patterns of the plurality of radiation patterns. The evaluator (4) is configured to estimate the direction of arrival based on the decomposed signals and based on information describing signal receiving characteristics of the antenna (2). The invention also refers to a corresponding method.

Inventors:
SAHINBAS BURAK (DE)
WEISGERBER LARS (DE)
SCHÜHLER MARIO (DE)
WANSCH RAINER (DE)
Application Number:
PCT/EP2016/072539
Publication Date:
March 29, 2018
Filing Date:
September 22, 2016
Export Citation:
Click for automatic bibliography generation   Help
Assignee:
FRAUNHOFER GES FORSCHUNG (DE)
International Classes:
G01S3/32; G01S3/40; G01S3/74; H01Q21/06
Foreign References:
US5608411A1997-03-04
Other References:
STEFANO MADDIO ET AL: "Direction of Arrival estimation of an arbitrary oriented targets by a highly directive antenna and coarse RSSI data", MICROWAVE CONFERENCE (EUMC), 2012 42ND EUROPEAN, IEEE, 29 October 2012 (2012-10-29), pages 140 - 143, XP032328231, ISBN: 978-1-4673-2215-7
MADDIO STEFANO ET AL: "A closed-form formula for RSSI-based DoA estimation with Switched Beam Antennas", 2015 EUROPEAN RADAR CONFERENCE (EURAD), EUMA, 9 September 2015 (2015-09-09), pages 341 - 344, XP032824571, DOI: 10.1109/EURAD.2015.7346307
MUEHLDORF E I ET AL: "Polarization Mismatch Errors in Radio Phase Interferometers", IEEE TRANSACTIONS ON AEROSPACE AND ELECTRONIC SYSTEMS, IEEE SERVICE CENTER, PISCATAWAY, NJ, US, vol. AES-10, no. 2, 1 March 1972 (1972-03-01), pages 135 - 141, XP011167139, ISSN: 0018-9251
Attorney, Agent or Firm:
STÖCKELER, Ferdinand et al. (DE)
Download PDF:
Claims:
Claims

Apparatus (1 ) for estimating a direction of arrival,

wherein the apparatus (1 ) comprises an antenna (2), a beamforming network (3), and an e valuator (4),

wherein the antenna (2) is configured to receive signals,

wherein the antenna (2) is circularly polarized,

wherein the antenna (2) comprises a plurality of different radiation patterns, wherein the beamforming network (3) is configured to provide based on signals received by the antenna (2) decomposed signals that are received by associated radiation patterns of the plurality of radiation patterns, and

wherein the e valuator (4) is configured to estimate the direction of arrival based on the decomposed signals and based on information describing signal receiving characteristics of the antenna (2).

Apparatus (1 ) of claim 1 ,

wherein the signal source (10) is linearly polarized. Apparatus (1 ) of claim 1 or 2,

wherein the antenna (2) is either right hand circularly polarized or left hand circularly polarized.

Apparatus (1 ) of any of claims 1 to 3,

wherein the antenna (2) comprises a plurality of antenna elements (5),

wherein the antenna elements (5) of the plurality of antenna elements (5) comprise different radiation patterns, and

wherein the beamforming network (4) is configured to provide the decomposed signals so that the decomposed signals are received with radiation patterns of individual antenna elements (5) or with a combination of radiation patterns of at least two antenna elements (5).

Apparatus (1 ) of claim 4,

wherein the antenna elements (5) of the plurality of antenna elements (5) are all either right hand circularly polarized or left hand circularly polarized.

6. Apparatus (1 ) of claim 4 or 5, wherein the antenna elements (5) of the plurality of antenna elements (5) are located in a plane.

Apparatus (1 ) of any of claims 1 to 6,

wherein the evaluator (4) is configured to estimate an inclination (p) of a signal source (10) emitting the signals received by the antenna (2).

Apparatus (1 ) of any of claims 1 to 7,

wherein the information describing signal receiving characteristics of the antenna (2) refers to a position of a signal source (10) emitting the signals received by the antenna (2) relative to the antenna (2).

Apparatus (1 ) of claim 8,

wherein the information describing signal receiving characteristics of the antenna (2) refers to a co-elevation (Θ) and an azimuth (φ) describing the position of the signal source (10) relative to the antenna (2).

Apparatus (1 ) of claim 8 or 9,

wherein the information describing signal receiving characteristics of the antenna (2) refers to an inclination (p) of the signal source (10) relative to the antenna (2).

Apparatus (1 ) of any of claims 8 to 10,

wherein the information describing signal receiving characteristics of the antenna (2) comprises a set of steering vectors (A) depending on angles (Θ, cp) describing the position of a signal source (10) relative to the antenna (2) and depending on an angle (p) describing an inclination of a polarization of the signal source (10) relative to the antenna (2).

Apparatus (1 ) of claim 1 1 ,

wherein the set of steering vectors (A) is based on dividing a range of possible inclinations (p) of the signal source (10) into partition intervals.

13. Apparatus (1 ) of claim 12,

wherein a width of the partition intervals is a measure for a resolution concerning an estimation of the inclination (p) of the signal source (10).

14. Apparatus (1 ) of claim 12 or 13,

wherein the partition intervals are set based on a projection similarity measure, wherein for chosen reference values of the angles (9ref, cpref) describing the position of a signal source (10) relative to the antenna (2) a reference inclination (prei) is chosen as a reference steering vector,

wherein steering vectors belonging to same reference values of the angles (9ref, c ref) are projected on the reference steering vector, and

wherein in case an obtained projection value lies within a given value interval, then the steering vectors are considered as identical.

15. Apparatus ( 1 ) of claim 14,

wherein steering vectors belonging to same reference values of the angles (9ref, c ref) are projected on the reference steering vector and normalized for different inclination values (p).

16. Apparatus (1 ) of claim 14 or 15,

wherein the given value interval is given by a lower threshold and a value close to one.

17. Method for estimating a direction of arrival,

comprising:

receiving signals emitted by a signal source (10) with a circularly polarized antenna (2) comprising a plurality of different radiation patterns,

providing based on the received signals decomposed signals that are received by associated radiation patterns of the plurality of radiation patterns, and

estimating the direction of arrival based on the decomposed signals and based on information describing signal receiving characteristics of the antenna (2).

Description:
Apparatus for estimating a direction of arrival and corresponding method

Specification The present invention relates to an apparatus for estimating a direction of arrival. Further, the invention relates to a method for estimating a direction of arrival.

Direction-of-arrival (DoA) estimation of a signal emitted by a source or signal source of unknown location requires in the state of art an receive antenna array. Relying on the de- lay the signal is received amongst the antennas, a conclusion on the direction of the source can be drawn.

Not only the location of the signal source, but usually also the source orientation concerning the emission of signals is unknown. Thus, the polarization of the source signal can be arbitrarily inclined with respect to the receive array. Therefore, receive arrays of linearly polarized elements are usually deployed, because the phase difference amongst the elements is preserved for varying source inclination. This, however, can cause polarization mismatches, possibly leading to weak receive signals, whose direction of arrival cannot be accurately estimated.

Circularly polarized receive antennas do not suffer from polarization losses caused by inclined sources, assuming a linearly polarized source signal. Yet because of the limited axial-ratio beam width, which is displayed by a circularly polarized antenna, the phase difference amongst the elements varies with varying source inclination (see Fig. 1 ).

The scientific community has been investigating various antenna arrays and concepts to increase the accuracy of direction finding (DF).

This includes approaches using multibeam antennas (MBA) [2, 3], Typically these anten- nas are linearly polarized and so are the signal sources. Thus, even if the signal source is arbitrarily oriented, the phase of the received signals stays mostly constant for a given impinging angle [4]. This allows for estimating of the DoA without considering source inclination, at the expense of a reduced estimation performance caused by polarization mismatches. To achieve a DoA estimation performance being independent of the source orientation, a circularly polarized direction finding antenna is preferred. However, for such an antenna, the received signal phases vary with changing signal source inclination, due to the limited axial ratio beamwidth.

There is a variety of methods for direction finding such as subspace-based methods, e.g. MUSIC [5]. Amongst all methods, the accuracy strongly depends on the set of steering vectors used and how accurate those represent the antenna characteristics. This includes polarization and the inclination of the source with respect to the antenna. As abovemen- tioned, polarization mismatch causes inaccurate DoA estimates.

To cope with polarization mismatch, receive antenna arrays with dual polarization can be used. These arrays have different steering vectors for different source inclination and are therefore referred as polarization sensitive arrays [6], Each element provides two signals, one for the first polarization and one for the second polarization. A set of two orthogonal polarization is typically used, e.g., linear vertical and linear horizontal polarization or right handed circular and left handed circular polarization (RHCP & LHCP). This, however, doubles the implementation effort, as twice as many signal branches are needed compared to an array supporting one polarization. Hence, an array supporting a single polari- zation only is preferred in practice.

The object of the invention is to provide an apparatus and a corresponding method for estimating the direction of arrival alternative to the state of art. The object is achieved by an apparatus for estimating a direction of arrival according to claim 1 and by a method for estimating a direction of arrival according to claim 17.

The object is achieved by an apparatus for estimating a direction of arrival. The apparatus comprises an antenna, a beamforming network, and an evaluator. The antenna is config- ured to receive - especially electromagnetic - signals, the antenna is circularly polarized, and the antenna comprises a plurality of different radiation patterns. The beamforming network is configured to provide based on signals received by the antenna decomposed signals that are received by associated radiation patterns of the plurality of radiation patterns. Further, the evaluator is configured to estimate the direction of arrival based on the decomposed signals and based on information describing signal receiving characteristics of the antenna. The apparatus comprises an antenna for receiving signals from the signal source. The antenna is circularly polarized and has accordingly just one circular polarization. The antenna comprises a plurality of different radiation patterns that allow to receive a signal emitted by the signal source with different spatial distributions of the sensitivity for receiving electromagnetic signals. Each radiation pattern has a specific spatial distribution for receiving signals. Hence, with the different radiation patterns a signal of the signal source is received with various different spatial sensitivity for - especially electromagnetic - signals. The beamforming network allows to decompose the received signals into decom- posed signals that are associated with the different radiation pattern. Hence, each decomposed signal is received with a radiation pattern of the plurality of available radiations patterns. In an embodiment, the radiation patterns refer to beams of the antenna. Finally, the evaluator estimates the direction of arrival based on the decomposed signals and based on information describing the signal receiving characteristics of the antenna. The information also refers to the different radiation patterns of the antenna. Thus, the apparatus is configured to estimate the direction of arrival based on received signals that are emitted by a signal source.

In an embodiment, the apparatus is configured to estimate additionally the polarization of the signal source. In an embodiment, the estimation of the polarization - or: the inclination of the polarization - is done by using only a single polarization by applying discrete reference inclination values as will be explained in the following.

In an embodiment, the apparatus is especially configured to estimate the direction of arri- val based on signals emitted by a signal source that is linearly polarized.

In an embodiment, the antenna is either right hand circularly polarized or left hand circularly polarized. According to an embodiment, the antenna comprises a plurality of antenna elements. The antenna elements of the plurality of antenna elements comprise different radiation patterns which are the basis for the plurality of radiation patterns of the antenna. Further, the beamforming network is configured to provide the decomposed signals so that the decomposed signals are received with radiation patterns of individual antenna elements or with a combination of radiation patterns of at least two antenna elements. Thus, the plural- ity of radiation patterns of the antenna are either individual radiation patterns of the antenna elements or are combinations of the individual radiation patterns.

Each antenna element - as a subunit of the antenna - comprises its own radiation pat- tern. At least some radiation patterns of the antenna elements differ from each other. In an embodiment, all radiation patterns differ from each other. The beamforming network decomposes the signals received by the antenna in decomposed signals where each decomposed signal is received either by a radiation pattern of one antenna element or by a combination of the radiation patterns of at least two antenna elements.

In an embodiment, the antenna elements are of one design and are, thus, identical to each other.

In one embodiment, the antenna elements of the plurality of antenna elements are all ei- ther right hand circularly polarized or left hand circularly polarized.

In an embodiment, the antenna elements of the plurality of antenna elements are located in a - common or joint - plane. Compare [1] for an example. In an embodiment, the antenna elements are printed antenna elements located and/or produced on the same sub- strate which becomes the common plane.

According to an embodiment, the evaluator is configured to estimate - additionally to the direction of arrival of the received signals - an inclination of a signal source emitting the signals received by the antenna. The inclination is relative e.g. to a normal plane or nor- mal vector of the antenna. In one embodiment, the inclination - called p - of 0" implies that the signal source - with regard to its emitting characteristics - is parallel to a plane of the antenna. A value of 90°, i.e. p = 90°, implies that the signal source is located perpendicular to this plane or that the signal source is parallel to a normal vector of this plane, for example, being a main direction of receiving signals. The inclination of the signal source is here also the inclination of its polarization.

The following embodiments refer to the information that allows the evaluator to estimate the direction of arrival and/or the inclination. In an embodiment, the information describing signal receiving characteristics of the antenna refers to a position of a signal source emitting the signals received by the antenna relative to the antenna. The information used by the e valuator refers to the position of the signal source relative to the antenna.

In an embodiment, the information describing signal receiving characteristics of the an- tenna refers to a co-elevation and an azimuth describing the position of the signal source relative to the antenna. In this embodiment, the position of a signal source relative to the antenna is described by two angles: co-elevation and azimuth.

In one embodiment, the signal source is a satellite associated in an embodiment with an orbit. Such an orbit allows accordingly in one embodiment to set one angle to a fixed value, e.g. for the azimuth: φ = 0°.

In an embodiment, the information describing signal receiving characteristics of the antenna refers to an inclination of the signal source relative to the antenna. In this embodi- ment, the inclination of the signal source is taken into account as it also effects the signals received by the antenna.

According to an embodiment, the information refers to the co-elevation and the azimuth describing the position and refers to the inclination of the signal source.

In an embodiment, the information describing signal receiving characteristics of the antenna comprises a set of steering vectors depending on angles (named Θ for the co- elevation and φ for the azimuth) describing the position of a signal source (or the direction from which the signal received by the antenna stems) relative to the antenna and depend- ing on an angle describing an inclination of a polarization of the signal source relative to the antenna. Hence, the evaluation uses a set of steering vectors which are depending on three angles: two for the direction relative to the antenna (or more precisely to a plane and/or a normal vector of the antenna) and one for the inclination of the polarization of the signal source which is here called the inclination of the signal source.

In an embodiment, the set of steering vectors is based on dividing a range of possible inclinations (named by the angle p) of the signal source into partition intervals. The possible range of inclinations is in one embodiment given by -90° <= p <= 90° and in a different embodiment by -180° <= p <= 180°. The range of inclinations is divided in this embodi- ment into partition intervals. The intervals define in one embodiment the resolution of the estimation of the evaluation as the inclinations angles within the intervals are not discrimi- nated. Thus, each interval is associated with a reference inclination and with a certain width defining an interval of inclinations and a resolution of the estimation regarding the inclination. Hence, in one embodiment, a width of the partition intervals is a measure for a resolution concerning an estimation of the inclination of the signal source. The intervals are associated with one reference inclination p re? for the respective interval. Hence, if an interval has a width of Δρ, then no discrimination is done for inclinations p within the range (p ref - Δρ/2) <= p <= (pref + Δρ/2). Thus, the resolution is given by the width of the respective partition.

In an embodiment, the partition intervals are set based on a projection similarity measure.

For defining the partition intervals, the following steps are performed:

for chosen reference values of the angles (named preferably Q, ef and cp ref ) describing the position of a signal source relative to the antenna a reference inclination (named p ref ) is chosen as a reference steering vector,

steering vectors - named A and being dependent on the three mentioned angles - belonging to same reference values of the angles (the values 9 re < and cp reT ) are projected on the reference steering vector, and

wherein in case an obtained projection value lies within a given value interval, then the steering vectors are considered as identical.

The value interval is in one embodiment given by a lower threshold and by a upper limit of 1 . Hence, if the projection is close to one, the steering vectors are almost identical to the reference steering vector and are, accordingly to an embodiment, considered to belong to the respective partition interval of the reference steering vector.

In an embodiment, steering vectors belonging to same reference values of the angles are projected on the reference steering vector and normalized for different inclination values. Thus, in this embodiment, also a normalization is performed for calculating the projection values. This is also relevant for defining the above mentioned value interval.

In an embodiment, the given value interval is given by a lower threshold and a value close to one. The upper value of 1 is especially based on the normalization of the inclination values. in one embodiment, an overlap of partition values is reduced. Hence, the reference values and the width is set accordingly.

In an embodiment, the set of steering vectors to be used by a following subspace-based estimation such as MUSIC, ESPRIT, Maximum Likelihood, root-MUSIC, CAPON among others [7] for evaluating the decomposed received signals comprises steering vectors (named A) that a dependent on the positions angles (co-elevation Θ and azimuth φ) of the signal source and reference inclinations (p ref ) defining partition intervals. The partitions are used as the multibeam antenna (MBA) allows to estimate a source inclination just within an interval, i.e. its resolution. It was found that this resolution is different for each co-elevation angle Θ (and the azimuth angle φ in generic case). As a result, partitions are used in one embodiment to define the resolution capability of the MBA and they are used later on to form the spectrum indicating the source inclination. They are partitioning only the inclination (p) domain. In one embodiment, the partitions are not required for the estimation of Θ and φ in the case that the antenna has a uniform and very sharp resolution for them throughout their respective domains.

The partitions are constructed by an initial choice of a reference inclination value p Tei after fixing each co-elevation (theta, Θ) and in the general case also the azimuth (phi, φ). Then, for keeping the co-elevation Θ (and the azimuth φ) fixed, the steering vector A(p) is projected on A(p r ef) for all rho values. This is a similarity measure referred and explained in the following. Then, in one embodiment, it is checked for that Θ (and φ), for which p values, the similarity projection produces almost 1 . For that p interval (call it resolution), the steering vector A(p ref ) is used only in the spectrum and the processing of the received signals allow to tell if the received signal is within that rho interval or not, in terms of source inclination. If there are some p intervals that are not representable by that A(p ref ), a different p ref is selected accordingly and same procedure is applied. Other reference inclinations p are selected according to the resolution, after first p ref produces it.

Thus, if required, other reference inclination values p ref are found to yield similarity projec- tion that covers all together above the threshold line for that Θ (and φ). Studies have shown that the number of reference inclination values p ref depends on the characteristics of the antenna.

In short, inclination partitions are used for determining a resolution capability of the anten- na in terms of source inclination. Together with this, spectra are built and used in one embodiment to estimate the direction of arrival and the inclination of the polarization. For this, a multibeam antenna is used, wherein the antenna is unique by having distinct axial ratio patterns in its partitions. The object is also achieved by a method for estimating a direction of arrival.

The method comprises at least the following steps:

• receiving signals emitted by a signal source with a circularly polarized antenna comprising a plurality of different radiation patterns,

· providing based on the received signals decomposed signals that are received by associated radiation patterns of the plurality of radiation patterns, and

• estimating the direction of arrival based on the decomposed signals and based on information describing signal receiving characteristics of the antenna. The above discussed embodiments and features of the apparatus can also be realized via the method and vice versa.

The invention will be explained in the following with regard to the accompanying drawings and the embodiments depicted in the accompanying drawings, in which:

Fig. 1 shows a comparison between phase differences of lobes 2, 3 and 4 with respect to lobe 1 of a multibeam antenna,

Fig. 2 illustrates an embodiment of the apparatus with a 2 x 2 mulitbeam antenna structure and the direction finding scenario including the signal source and its coordinates relation,

Fig. 3 shows example axial ratio values of the respective lobes U to L 5 on the XZ plane, Fig. 4 shows the results of a projection of A(0°, p) on the reference steering vector A(0°, 0°) for -90° <= p <= 90°,

Fig. 5 shows the results of a projection of A(20°, p) on the reference steering vec- tor A(0°, 0°) for -90° <= p <= 90°,

Fig. 6 shows the results of a projection of A(20°, p) on A(0°, p ref 0°) for p ref belonging to the set {-78°, -55°, -32°, -9°, 14°, 37°, 60°, 83°), Fig. 7 shows three spectra for distinct source inclinations: a) p = 0°, b) p = 45°, and c) p = 90° at Θ = -5°,

Fig. 8 shows three spectra for distinct source inclinations: a) p = 0°, b) p = 45°, and c) p = 90° at Θ = 50°,

Fig. 9 shows a different embodiment of the apparatus and

Fig. 10 shows another embodiment of the apparatus. Fig. 1 shows a comparison between phase differences of lobes 2, 3 and 4 with respect to lobe 1 of a multibeam antenna which is similar to the antenna given by [1 ]. Shown is on the x-axis the observation angle Θ in ° and on the y-axis the phase differences δ in °. These phase differences are given for the lobes L 2 , L 3l and L 4 with respect to the first lobe L-] of the four beams of the multibeam antenna taught by [1 ] at an azimuth φ = 0° and to an inclination p = 0°. For each lobe two curves are shown: RP: source inclination p = 0° and MS: source inclination p = 90°, respectively.

An embodiment of the apparatus 1 is shown in Fig. 2. In the apparatus 1 , a multibeam antenna (MBA) 2 is utilized having diverse radiation patterns. The antenna elements 5 support only a single circular polarization: RHCP (in the shown embodiment) or LHCP. If the amplitude and phase response of the patterns are used in a subspace-based method, for instance, a direction of arrival (DoA) estimation of the source signal stemming from the signal source 10 is possible along with an estimation of the source leaning with inclination angle p of the signal source 10 relative to the antenna 2. Fig. 2 shows the geometry by example of a 2-by-2 elements array providing here five beams, based on the antenna provided by [1] with 4 antenna elements 5: four beams are given by the radiation patterns of the four antenna elements and one beam is the result of the combination of the radiation patterns of the antenna elements. That all antenna ele- ments 5 are right handed circularly polarized is indicated by the white arrows.

The array aperture is in the XY plane. The signal source 10 appears at arbitrary co- elevation Θ and azimuth <j> (here: φ = 0°) and is arbitrarily oriented with respect to the receive array denoted by the inclination angle p. A beam-forming network (BFN) 3 provides the excitation vectors necessary for forming the beam patterns, which can be selected via an RF switch (see Fig. 10). The resulting decomposed signals are processed by the eval- uator 4.

The discussion in the following is based on the single source scenario depicted in Fig. 2 for φ = 0°. The method is, however, generally applicable to every multibeam antenna and multiple sources.

The signal model reads

X

where x is the vector of complex received signals, A is the complex responses of the lobes including the round trip phase difference of the signal, N is the vector of complex white Gaussian noise having a variance of σ 2 and s-i is the complex source signal.

A covariance matrix is then formed. It is given by R = xx H = (A.s r + N) (As ' i + N)H = As l Sl H A H + σ 2 1

(2) due to the fact that noise is uncorrelated with the signal.

Equivalently, there is

R a' s 1 AA (3)

The rest of the algorithm depends on the selection of the subspace-based method, which makes use of the covariance matrix defined above, so the following results are legitimate for any subspace-based estimation such as MUSIC, ESPRIT, Maximum Likelihood, root- MUSIC, CAPON among others [7].

For the estimation, the steering vectors Α(θ, p) for each distinct Θ and p are required. This is for the discussed case φ = 0°. Hence, in a general form, the steering vectors Α(θ, φ, p) are required.

The two orthogonal source polarization Α(θ, 0°) and Α(θ, 90°) are measured for -90° < Θ < 90° per each beam output. Using these two steering vectors, the rest follows from

for the each beam i (in our example i = 1. 2, 3, 4, 5). Then, steering vectors for an arbitrary source inclination p is obtained.

The fact is, as long as the steering vector covariance matrix AA H is considerably different for each value of p, the polarization sensitivity is achieved.

Assuming all the beams are perfectly RHCP for all Θ range, then as the linearly polarized source leaning changes, the beam output signals are shifted with exactly the same phase, providing no change on the steering vector covariance matrix. As a result, the polarization sensitivity is not achieved. However, different levels of achieved RHCP provides the sensitivity. This is where the multibeam antenna MBA 2 comes into action. Having an MBA yields an axial ratio pattern as exemplarily portrayed in Fig. 3 in which the axial ratios AR in dB for the five beams of the antenna 2 in Fig. 2 is shown. The AR is given depending on the observation angle Θ in

Each beam or each sub-set of beams - if some beams repeat because of symmetry - has its unique axial ratio characteristics with respect to observation direction. This is peculiar to multibeam antennas and results in polarization sensitivity. Thus, a multibeam antenna as used in the apparatus provides a well resolved direction of arrival estimation characteristics and source inclination estimation.

Before the DoA estimation is done, the steering vector set matching the source inclination is needed. The set is in one embodiment obtained by measurements and is in a different embodiment obtained by measurements and by calculations based on these measurements.

In one embodiment, the range of possible source inclination p is divided into partitions. The partitioning intervals are arranged by a projection similarity measure. For a certain value θ = θ Γβί , a reference p ref is selected, forming the reference steering vector. Then, all the steering vectors having the same d ref value is projected on the reference steering vector and normalized for different p values:

This formula is given for φ = 0°. The general formula is accordingly given for proj(9 re( , cp ref , p).

If the projection value is almost 1 , the steering vectors can be considered identical. In practice, a lower boundary, i.e., a threshold, is selected for differentiating between steering vectors. Above the boundary, steering vectors are considered identical, below steering vectors are different. The threshold is selected with respect to axial ratio beamwidth, signal demodulation and noise; it should be close to 1 , for example, 0.998.

Towards normal direction of the MBA, the projection is close to 1 independent of p as shown in Fig. 4 for p ref = 0° because of the absence of polarization sensitivity (compare Fig. 3). In Figures 4 to 6 the norm of projection is shown in dependency on the inclination angle p in °.

As the co-elevation Θ values get larger, there appears an immediate polarization sensitivi- ty (axial ratio diversity) for some partitions. In Fig.5, the situation is demonstrated by taking 9 r ef = 20° and p ref = 0° for which, there is a polarization sensitivity having a 30° resolution. Fig. 5 also shows that the resolution is about +- 17° and the partition interval will have accordingly a width of about 34°. The resolution is here given by the values for the norm of projection greater than 0.998, i.e. close enough to the value of the reference incli- nation.

Assuming, for the time being, the resolution is identical for all p ref , the range of inclination angles p can be divided into equidistant partitions, as shown in Fig. 6. For Θ = 20° eight partitions are required to cover the area above the threshold line. The resolution interval boundaries are indicated by two circles.

In the example, the following steering vectors were used: A(20°, -78°), A(20°, -55°), A(20°, -32°), A(20°, -9°), A(20°, 14°), A(20°. 37°), A(20°, 60°), and A(20°, 83°). The selection is made in one embodiment regarding to this similarity projection method for the other instances of Θ and p. For the interval having no polarization sensitivity, only a single steering vector, i.e. , a single partition, is used for forming the spectrum. This fact decreases the required number of comparisons. After determining the required steering reference vectors and partitioning the p spectrum for each Θ (and φ), the estimation can be initiated with any subspace based method. The number of partitions, which directly corresponds to the polarization sensitivity of the MBA, increases as Θ increases. In the following table, a sample partitioning for the given MBA of the embodiment of the apparatus of Fig. 2 is demonstrated, which is formed using equidistant approach with some overlapping. Here, for the sake of simplicity, the azimuth is still given by φ = 0°.

° t -28 ° u f28 °, 32 tt j to

-37 °, 33 0= LJ [33 ° , 37% 10

-42 ° s -38 ° : U [38 °, 42 °] 10

-47 °, -43 ° : u [43 ° , 47°] 10

-52 °, -48 ° : u (48 % 52 % 12

-57 ° ; -53 * : u i53 °, 57°j 12

The table gives example numbers of steering vector p partitions for the given range of the co-elevation Θ (here: the azimuth φ is set to 0°). As seen in Fig. 6, for some steering vectors there appears an overlapping above the threshold line as a consequence of the equidistant partitioning. The resolution at the respective steering vectors could therefore be doubled in the worst case, as the estimator favors one of the p intervals in the spectrum while yielding a coarse estimation of the source polarization. By having nonequidistant partitioning, there exists no overlapping as well, hence the MBA achieves a nonoverlapping and a bounded resolution capability. In short, the method yields a guide on how to select the reference partitioning steering vector set, having the constraint of covering all the field above the threshold line. Using the MUSIC algorithm for the verification, some test results are presented to show the capability of muitibeam antenna to resolve the direction of arrival of signals together with the increasing sensitivity to source inclination p as the co-elevation Θ increases. Fig. 7 demonstrates the output of the partitioned spectrum with the source being placed at -5° and for three different inclinations: p = 0°, 45°, and 90°.

Fig. 8 demonstrates the result for the same source inclination set when the source is placed at co-elevation Θ = 50°.

In both figures, Fig. 7 and Fig. 8, it can be seen that the range of possible inclinations p is divided into different partitions. Hence, the spectrum obtained by processing the decomposed signals with regard to their amplitudes and phases and using the set of steering vectors describing the spatial distribution of the signal receiving sensitivity of the different radiation patterns and describing effects of an inclination of the signal source is analyzed based on the different partitions with regard to the inclination of the signal source. The further processing is done by any subspace-based method according to the state of art, e.g. MUSIC. The results supports that the MBA achieves the polarization sensitivity together with its axial ratio pattern. There appears a very good accuracy on direction of arrival estimation, together with the partitioned source inclination estimation. An aspect here is that, the MBA suppresses the high axial ratio angles of the beams with its smaller amplitude responses at those angles, which acts like a weighting, yielding a polarization sensitivity together with a more robust direction of arrival estimation, thus using only RHCP beams without any secondary polarization pattern.

An embodiment of an apparatus 1 for DoA and polarisation estimation is portrayed in Fig. 9.

The "Antenna Array" 2 with different antenna elements 5 connects to an beam forming network (BFN) 3 - e.g. realized by a Butler matrix - which decomposes the received antenna signals into decomposed signals. The estimation is performed in the "Signal Processing " device or evaluator 4. The steering vector set of the respective antenna 2 ob- tained from measurement or simulation is provided via a storage unit 6. The evaluator 4 provides estimations for the direction of arrival and the inclination p of the polarization of the signal source. A controller 7 controls the apparatus 1. The separation of the BFN and the Signal Processing is logical. The BFN can be also part of the Signal Processing unit, representing a digital signal decomposition. Another embodiment of the apparatus 1 for DoA and polarization estimation is portrayed in Fig. 10.

The "Antenna Array" 2 with antenna elements 5 (having all the same circular polarization) connects to an BFN 3, which decomposes the antenna signals into decomposed signals. In contrast to the embodiment in Fig. 9, only one decomposed signal is selected at an instant of time via a switch unit 8. Supposing an analogue implementation of the BFN 3 and the switch unit 8, this allows for reducing the number of signal branches to one. The selection of the signal is controlled by the controller 7. It ensures a time-synchronous switching. The estimation is performed in the evaluator 4. The steering vector set of the respective antenna 2 is provided by the storage unit 6. The separation of the beam forming network 3, the switch unit 8, and the evaluator 4 is logical. The BFN 3 and the switch unit 8 are in a different embodiment part of the evaluator 4, representing a digital signal decomposition and selection. The provided apparatus and method provide at least the following benefits:

Reduced costs for implementation as single polarization requires only a single signal branch for each element, halving the effort compared to dual polarized solutions. Because of the diverse patterns, the DoA and the polarization (inclination, orientation) of the source signal can be concurrently estimated using an array with one polarization.

Applying partitioning of the range of inclination angles reduces the signal processing effort for estimation, allowing for usage of cheaper processing platforms such as microcontrol- lers (e.g. ARM family) or digital signal processors.

Possible application areas are:

Direction and orientation estimation of radio transmitters/transponders such as RFID transponders, wireless sensors, mobile devices (e.g. mobile phones, laptops, tablet computers), vehicles, aircrafts. Direction finding for military purpose and security (radio reconnaissance).

Sensor applications: determination of inclination / positioning / orientation

Although some aspects have been described in the context of an apparatus, it is clear that these aspects also represent a description of the corresponding method. Analogously, aspects described in the context of a method step also represent a description of a corresponding block or item or feature of a corresponding apparatus.

Some embodiments according to the invention comprise a data carrier having electronically readable control signals, which are capable of cooperating with a programmable computer system, such that one of the methods described herein is performed. Generally, embodiments of the present invention can be implemented as a computer program product with a program code, the program code being operative for performing one of the methods when the computer program product runs on a computer. The program code may for example be stored on a machine readable carrier. Other embodiments comprise the computer program for performing one of the methods described herein, stored on a machine readable carrier.

In other words, an embodiment of the inventive method is, therefore, a computer program having a program code for performing one of the methods described herein, when the computer program runs on a computer.

A further embodiment of the inventive methods is, therefore, a data carrier (or a digital storage medium, or a computer-readable medium) comprising, recorded thereon, the computer program for performing one of the methods described herein. The data carrier, the digital storage medium or the recorded medium are typically tangible and/or non-transitory.

A further embodiment of the inventive method is, therefore, a data stream or a sequence of signals representing the computer program for performing one of the methods de- scribed herein. The data stream or the sequence of signals may for example be configured to be transferred via a data communication connection, for example via the Internet. A further embodiment comprises a processing means, for example a computer, or a programmable logic device, configured to or adapted to perform one of the methods described herein.

A further embodiment comprises a computer having installed thereon the computer program for performing one of the methods described herein.

A further embodiment according to the invention comprises an apparatus or a system con- figured to transfer (for example, electronically or optically) a computer program for performing one of the methods described herein to a receiver. The receiver may, for example, be a computer, a mobile device, a memory device or the like. The apparatus or system may, for example, comprise a file server for transferring the computer program to the receiver.

In some embodiments, a programmable logic device (for example a field programmable gate array) may be used to perform some or all of the functionalities of the methods described herein. In some embodiments, a field programmable gate array may cooperate with a microprocessor in order to perform one of the methods described herein. Generally, the methods are preferably performed by any hardware apparatus.

The apparatus described herein may be implemented using a hardware apparatus, or using a computer, or using a combination of a hardware apparatus and a computer. The methods described herein may be performed using a hardware apparatus, or using a computer, or using a combination of a hardware apparatus and a computer.

References [1 ] L. Weisgerber and A. E. Popugaev, " Multibeam antenna array for RFID applications, ' ' in 2013 European Microwave Conference (Eu C), Nuremberg, Germany, Oct. 10-16, 2013, pp. 84-87.

[2] K. A. Gotsis, K. Siakavara, and J. N. Sahalos, "On the direction of arrival (doa) estimation for a switched-beam antenna system using neural networks, " IEEE Transactions on Antennas and Propagation, vol. 57, no. 5, pp. 1399-141 1 , May 2009.

[3] W. X. Sheng, J. Zhou, D. G. Fang, and Y. C. Gu, "Super-resolution DOA estimation in switch beam smart antenna," in Proc. 5th Int. Symp. Antennas, Propagation and EM Theory ISAPE 2000, Aug. 2000, pp. 603-606.

[4] L. Weisgerber and A. E. Popugaev, "Antenna concepts for DoA estimation in traffic scenarios," in Proc. IEEE-APS Topical Conf. Antennas and Propagation in Wireless Communications (APWC), Aug. 2014, pp. 423-426.

[5] R. Schmidt, "Multiple emitter location and signal parameter estimation," IEEE Transactions on Antennas and Propagation, vol. 34, no. 3, pp. 276-280, Mar. 1986.

[6] Y. Han, Q. Fang, F. Yan, M. Jin, and X. Qiao, "Joint DOA and polarization estimation for unequal power sources based on reconstructed noise subspace," Journal of Systems Engineering and Electronics, vol. 27, no. 3, pp. 501-513, Jun. 2016.

[7] J. Foutz, A. Spanias, and M. K. Banavar, Narrowband Direction of Arrival Estimation for Antenna Arrays, ser. Synthesis Lectures on Antennas #8. Morgan & Clay- pool, 2008.