WO2000067042A1 - Robust estimation of doa for antenna arrays - Google Patents

Robust estimation of doa for antenna arrays Download PDF

Info

Publication number
WO2000067042A1
WO2000067042A1 PCT/US2000/007472 US0007472W WO0067042A1 WO 2000067042 A1 WO2000067042 A1 WO 2000067042A1 US 0007472 W US0007472 W US 0007472W WO 0067042 A1 WO0067042 A1 WO 0067042A1
Authority
WO
WIPO (PCT)
Prior art keywords
doa
dbf
complex
signal
subarray
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Ceased
Application number
PCT/US2000/007472
Other languages
English (en)
French (fr)
Inventor
Sam Mordochai Daniel
Stephen Chihhung Ma
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Motorola Solutions Inc
Original Assignee
Motorola Inc
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Motorola Inc filed Critical Motorola Inc
Priority to DE60020693T priority Critical patent/DE60020693T2/de
Priority to EP00919499A priority patent/EP1177456B1/en
Priority to AU40178/00A priority patent/AU4017800A/en
Priority to JP2000615827A priority patent/JP2002543436A/ja
Publication of WO2000067042A1 publication Critical patent/WO2000067042A1/en
Anticipated expiration legal-status Critical
Ceased legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S3/00Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received
    • G01S3/02Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received using radio waves
    • G01S3/74Multi-channel systems specially adapted for direction-finding, i.e. having a single antenna system capable of giving simultaneous indications of the directions of different signals
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S3/00Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received
    • G01S3/02Direction-finders for determining the direction from which infrasonic, sonic, ultrasonic, or electromagnetic waves, or particle emission, not having a directional significance, are being received using radio waves
    • G01S3/04Details
    • G01S3/06Means for increasing effective directivity, e.g. by combining signals having differently oriented directivity characteristics or by sharpening the envelope waveform of the signal derived from a rotating or oscillating beam antenna
    • HELECTRICITY
    • H01ELECTRIC ELEMENTS
    • H01QANTENNAS, i.e. RADIO AERIALS
    • H01Q3/00Arrangements for changing or varying the orientation or the shape of the directional pattern of the waves radiated from an antenna or antenna system
    • H01Q3/26Arrangements for changing or varying the orientation or the shape of the directional pattern of the waves radiated from an antenna or antenna system varying the relative phase or relative amplitude of energisation between two or more active radiating elements; varying the distribution of energy across a radiating aperture

Definitions

  • This invention relates generally to phased array antennas and, more particularly, to a method and apparatus for robust estimation of directions of arrival for beams associated with antenna arrays.
  • Satellite communication systems have used phased array antennas to communicate with multiple users through multiple antenna beams.
  • efficient bandwidth modulation techniques are combined with multiple access techniques, and frequency separation methods are employed to increase the number of users.
  • frequency separation methods are employed to increase the number of users.
  • electronic environments are becoming increasingly dense, more sophistication is required for wireless communication systems.
  • DBF systems have been developed for use in communications systems and radar systems.
  • DBF systems require accurate direction of arrival (DOA) information to efficiently position beams and nulls in their antenna radiation/reception patterns. In DBF systems in which transmitters and receivers move relative to each other, DOA information is continually updated to maintain accuracy.
  • DOA direction of arrival
  • a DBF system can improve its allocation of resources by placing beams on active and high-traffic areas, while avoiding unnecessary coverage of large inactive regions. This is particularly important when satellites are being used because the efficient determination of directions of arrival can decrease the required processing load and decrease the amount of on-board power required.
  • FIG. 1 shows a simplified block diagram of a satellite communication system within which the methods and apparatus of the invention can be practiced
  • FIG. 2 illustrates a perspective drawing of a planar array of antenna elements and an incident signal from a generalized point source
  • FIG. 3 shows a simplified block diagram of a direction of arrival (DOA)-aided DBF subsystem that includes a digital beamformer and direction of arrival estimator (DOAE) in accordance with a preferred embodiment of the invention
  • FIG. 4 illustrates a flow diagram of a procedure for determining directions of arrival in a DOA-aided DBF subsystem in accordance with a preferred embodiment of the invention.
  • FIG. 5 shows a graph that illustrates the directional signal activity estimated by the DOA algorithm.
  • the present invention provides a method and apparatus that increase the frequency and code reuse factor in communication systems by accurately and efficiently determining directions of arrival, thereby, allowing beams to be more accurately directed and closely spaced.
  • the method and apparatus of the invention also provide more efficient processing of antenna beam patterns in communication systems.
  • the method and apparatus of the invention are especially significant for use with non-geostationary satellites in satellite communication systems.
  • Digital beamforming is essentially an open-loop concept.
  • a necessary input to a digital beamformer is the DOA information necessary to point the beams and nulls in the desired directions. When the DOA information is not precise, then these directions are not exactly correct, and the beams point away from the actual sources accordingly. In addition, their mutual nulls are shifted from their ideal locations, degrading the overall performance.
  • the invention combines an enhanced DOA estimation algorithm with a DBF based system to significantly improve the capacity of current and future communication systems, while remaining compatible with existing modulation techniques.
  • digital beamforming techniques are enhanced by, among other things, determining directions of arrival for desired and undesired incident signals.
  • FIG. 1 shows a simplified block diagram of a satellite communication system within which the method and apparatus of the invention can be practiced.
  • FIG. 1 illustrates a single satellite 110 with digital beamformers in a typical spectrum-sharing scenario. As illustrated, there are several communication paths between satellite 110 and terrestrial-based communication devices 120, 130, and 140.
  • a first beam 125 can be used to establish a link between satellite 110 and terrestrial-based communication device 120.
  • a second beam 135 can be used to establish a link between satellite 110 and terrestrial-based communication device 130.
  • a third beam 145 can be used to establish a link between satellite 110 and terrestrial-based communication device 140.
  • DOA estimates are determined at satellite 110 for first beam 125, second beam 135, and third beam 145.
  • a number of communications satellites such as illustrated by satellite 110, reside in non-geostationary orbits and are interconnected using crosslinks (not shown). In an alternate embodiment, the communications satellites are not all interconnected. For example, some communication satellites may be in geostationary orbits.
  • satellites In non-geostationary orbits, satellites can move at high speed relative to any given point on the surface of the earth. This means that these satellites can come into view at various times with respect to a point on the surface of the earth.
  • a first interference source 170 and second interference source 180 are also illustrated.
  • a first interference signal path 175 exists between satellite 110 and first interference source 170.
  • a second interference signal path 185 exists between satellite 110 and second interference source 180.
  • DOA estimates are determined at satellite 110 for interference signal paths 175 and 185.
  • Digital beamforming methods require that users, illustrated by 120, 130, and 140, are detected and that DOA information be obtained or estimated for them. In some cases, beamforming calculations can be significantly degraded by errors in the DOA data.
  • Several algorithms for estimating directions of arrival are known to those skilled in the art, including decoupled maximum likelihood (DML) and multiple signal classification (MUSIC), but these algorithms have limitations which restrict their performance in a dynamically changing environment.
  • DML decoupled maximum likelihood
  • MUSIC multiple signal classification
  • satellite 110 employs a digital beamformer (not shown) and a DOAE (not shown) as described below. Satellites 110 and terrestrial- based communication devices 120, 130, and 140 can be viewed as nodes in satellite communication system 100. Those skilled in the art will recognize that the below- discussed features of a preferred embodiment of the invention can be practiced at any node of satellite communication system 100 or any node of other radio frequency (RF) communications systems.
  • RF radio frequency
  • the directions of arrival for the desired signals and undesired signals in the receive mode and transmit mode are estimated.
  • DOA information is used for the positioning of beams and nulls in both the receive and the transmit modes.
  • DOA estimates are periodically updated, and control matrices for the digital beamformers are periodically adjusted to maintain the positions for the beam and nulls.
  • DOA estimation algorithms can also be used to produce combinations of narrow and wide nulls or, more generally, nulls having desired widths. This flexibility also allows concurrent beams with different beamwidths and different null widths.
  • the satellite In a receive mode, the satellite, at a particular point in time, desirably points a receive beam at a particular terrestrial-based communication device while preferably providing nulls in the antenna's receive pattern in the direction of any interfering signal transmitters. Accordingly, any interference received on an undesired signal path is significantly reduced.
  • interfering signal transmitters can be other users in this communication system, signal sources in other systems, or jamming signals.
  • at least one null in the receive antenna pattern of a satellite is directed toward and tracks each undesired signal, which is transmitted within the field of view of that satellite.
  • the field of view can be defined by a current operational field of view or by the entire field of view of the satellite.
  • the satellite In a transmit mode, the satellite desirably points at least one communications beam at a particular terrestrial-based communication device while preferably providing nulls in the antenna's transmit pattern in the direction of any known desired or undesired signal receivers. Again, directions of arrival are used to determine the angular positions for the beam and nulls.
  • the DOA can be determined by a DOAE, which can use, among other things, information associated with the location of the device. For example, GPS information can be used to determine the location of a node in the system.
  • a digital beamformer when employed in terrestrial-based communication devices 120, 130, and 140, desirably adjusts its transmit and receive antenna beam characteristics to point at least one beam at the desired satellite while directing at least one null in the direction of an interfering (undesired) signal.
  • an interfering signal can be associated with another satellite in this communications system or in another system.
  • the accuracy of the DOA estimates is controlled for the most part by the DOA estimation algorithm and its implementation.
  • the accuracy required is governed by the angular separation between desired users, the angular separation between interference sites, and the angular separation between desired users and interference sites.
  • FIG. 2 illustrates a perspective drawing of a planar array of antenna elements and an incident signal from a generalized point source.
  • a three-dimensional coordinate system is illustrated having an x-axis, a y-axis, and a z-axis.
  • Planar array 200 comprises a number of elements with mutual separation d x and d y in the x and y directions, respectively.
  • a receive mode a first set of elements is used
  • a second set of elements is used.
  • the center of the coordinate system is positioned at the array's geometrical center.
  • Point source 210 corresponds to the j-th co-channel signal, where r represents the distance from the origin of the coordinate system and 6j and ⁇ correspond to the elevation and azimuth angles, respectively. Desirably, (6j, ) are used to designate a DOA estimate for point source 210.
  • At least one signal path 215 exists between point source 210 and planar array 200.
  • signal path 215 has elevation and azimuth angles (6j and ⁇ ) associated with it.
  • Directions of arrival can be determined relative to planar array 200 or to point source 210.
  • Point source 210 can represent a receiver in one case, a transmitter in a second case, or a combined receiver/transmitter in another case.
  • Transceivers and transmitters can be desired signal sources and/or undesired signal sources.
  • Subarray 220 comprises a first subset of the elements in planar array 200.
  • Subarray 220 comprises a large number of antenna elements with individual complex weights, usually at baseband, represented by the vector w.
  • Reference antenna 230 comprises a second subset of the elements in planar array 200. For simplicity of illustration, a single frequency of operation / and uniform element spacing is assumed, although these restrictions are not required for the invention.
  • Reference antenna 230 has a substantially uniform gain over a desired field of view within a prescribed conical region.
  • Field of view (FOV) 235 is the operational FOV associated with reference antenna 230.
  • FOV Field of view
  • the look direction can be altered to a limited extent. For example, the look direction can be altered to achieve and maintain improved coverage of existing communication traffic as the satellite moves or to compensate for lost or heavy service demands of neighboring satellites in a constellation.
  • Footprint 240 is associated with the FOV of reference antenna 230.
  • Footprint 250 illustrates the potential coverage region for a reference beam associated with planar array 200.
  • the received baseband signal from reference antenna 230 is combined with a weighted sum of the individual baseband signals from the antenna elements in subarray 220.
  • Minimizing the composite received power by adjusting the subarray weights gives rise to a composite antenna power pattern that exhibits nulls and low-gain regions corresponding to incident directional signal activity.
  • Standard digital beamforming is effective in increasing the isolation among the beams beyond that of the sidelobe structure, because beams are formed with mutual narrow nulls among them. Typically, each formed beam can exhibit a single null toward the prescribed direction of an interfering beam.
  • standard digital beamforming is limited when considered with respect to its sensitivity to pointing errors, bandwidth, and orbit dynamics, particularly in the context of sticky beams.
  • Sticky beams are beams which have an endpoint that remains substantially fixed on the earth's surface as the satellite moves overhead.
  • the sharp nulls do not coincide with the actual directions. This occurs because the assigned beam directions are no longer centered over the user locations. Consequently, standard digital beamforming provides marginal benefit.
  • the limited effectiveness of the narrow nulls of standard digital beamforming is substantially eliminated by enhancing the DOA estimation algorithms to more accurately position beams and nulls. Further, wide mutual nulls add to the isolation among the formed beams even with the motion of the satellite relative to the ground sources.
  • FIG. 3 shows a simplified block diagram of a direction of arrival (DOA)-aided DBF subsystem that includes a digital beamformer and direction of arrival estimator (DOAE) in accordance with a preferred embodiment of the invention.
  • DOA-aided DBF subsystem 300 includes array-antenna 310, which comprises a plurality of receive (Rx) elements 312, a plurality of transmit (Tx) elements 314, a plurality of receiver/transmitter (Rx/Tx) modules 320, digital beamformer 330, DOAE 340, and controller 350.
  • array-antenna 310 which comprises a plurality of receive (Rx) elements 312, a plurality of transmit (Tx) elements 314, a plurality of receiver/transmitter (Rx/Tx) modules 320, digital beamformer 330, DOAE 340, and controller 350.
  • Array-antenna 310 includes elements that are preferably arranged in a linear or planar two-dimensional array; however, other array configurations are suitable. Received radio frequency (RF) signals are processed at the element level.
  • RF radio frequency
  • antenna elements a ! through a ⁇ are illustrated as receiving elements 312 for incident directional signals b 1 through b,. Desirably, there are J directional signals incident on a K-element receive array.
  • antenna elements a/ through a k -' are illustrated as transmitting elements 314 for transmitted directional signals b through b . Desirably, there are J' directional signals transmitted from a K'-element transmit array.
  • At least one separate transmit (Tx) array antenna is used, and at least one receive (Rx) array antenna is used.
  • Rx/Tx modules 320 comprise some separate receive functions and some separate transmit functions.
  • the plurality of receiving elements 312 and the plurality of transmitting elements 314 are controlled using digital beamforming techniques.
  • the antenna pattern from an array of receive elements 312 or an array of transmit elements 314 can be steered by applying linear phase weighting across the array. For example, an array pattern can be shaped by the amplitude and phase weighting of the outputs of the individual receive elements 312, and another array pattern can be shaped by the amplitude and phase weighting of the inputs of the individual transmit elements 314.
  • Rx/Tx modules 320 are coupled to Rx elements 312 and Tx elements 314. Preferably, at least one antenna element is coupled to a Rx/Tx module.
  • Rx/Tx modules 320 When operating in the receive mode, Rx/Tx modules 320 perform, among other things, the Rx functions of frequency down-conversion, filtering, amplification, and A/D conversion.
  • Rx/Tx modules 320 In response to received signals, Rx/Tx modules 320 generate digital data using in-phase (I) and quadrature (Q) A/D converters. I and Q digital data respectively represent real and imaginary parts of a complex analog signal envelope and are processed by DBF 330.
  • Rx/Tx modules provide K digital received signals (r, through r ⁇ ) to DBF 330.
  • Rx/Tx modules 320 desirably perform, among other things, the Tx functions of frequency up-conversion, filtering, amplification, and D/A conversion. D/A converters convert digital data into corresponding analog signals for each Tx array element.
  • Rx/Tx modules 320 generate signals suitable for transmission by Tx array elements from digital data received from DBF 330.
  • Rx/Tx modules are provided K' digital transmit signals (ry through r K -') by DBF 330.
  • Digital beamformer 330 is coupled to Rx/Tx modules 320. Digital data is exchanged between DBF 330 and Rx/Tx modules 320. DBF 330 implements beam forming and beam steering functions necessary to form antenna beam patterns with the desired characteristics. DBF 330 forms receive and transmit beams for the reception and transmission of directional signals with minimal inter-beam interference.
  • DBF 330 is coupled via 335 to controller 350.
  • Digital data is exchanged between DBF 330 and controller 350.
  • Digital data includes data for operating in the receive mode and the transmit mode as well as other data for control.
  • DBF 330 provides the received beam port signals (b. through bi ) to controller 350. These signals are optimal estimates of the incident signals (b. through ).
  • DBF 330 obtains the transmit beam port signals (b'> through b'j ) from controller 350.
  • the actual transmitted signals (b ⁇ through b ⁇ ) are optimal representations of the original versions (b' ⁇ through b'j).
  • One of the functions of the digital beamformer is to derive T and T' transformation matrices which, operating on received and transmitted signals, respectively, produce optimal estimates of incident and actually transmitted signals.
  • DBF 330 uses a best linear unbiased estimator
  • DOAE 340 is coupled 345 to DBF 330 and is coupled 355 to controller 350. Digital data is exchanged between DOAE 340, DBF 330, and controller 350. In a preferred embodiment, DOAE 340 obtains antenna configuration data from controller 350.
  • controller 350 can provide initial values for directions of arrival, calculated values for directions of arrival, and stored values of directions of arrival.
  • DOAE 340 provides, among other things, reference antenna weights, and subarray weights to DBF 330.
  • DOAE 340 uses knowledge of existing beams to enhance its search strategy for new directional incident signal activity. Desirably, DOAE 340 obtains existing beam information from DBF 330 and controller 350. Also, DOAE 340 uses, among other things, beam width and null width information to determine its beam assignment strategy.
  • DOAE 340 is used to compute, among other things, covariance matrices, and cross-correlation vectors as discussed below.
  • DOAE 340 employs a relaxation algorithm that is numerically stable and does not impose a requirement that the underlying covariance matrix have an inverse.
  • DOAE 340 comprises one or more parallel processors.
  • DOAE 340 stores data that serve as its instructions and that, when executed, cause DOAE 340 to carry out procedures that are discussed below.
  • DOAE 340 can be implemented using digital signal processors.
  • DOAE 340 can be implemented using special processors, which can include logarithm converters, inverse logarithm converters, and parallel processors.
  • the processors used in DBF 330, DOAE 340, and controller 350, can use logarithmic number system (LNS) arithmetic.
  • LNS-based arithmetic provides an advantage because multiplication operations can be accomplished with adders instead of multipliers.
  • a LNS-based processor can include log converters, summing circuits, weighting circuits, and inverse-log converters 1.
  • DBF 330 and DOAE 340 can share computational resources.
  • DOAE 340 uses an N-element receive subarray to determine the distribution of directions of arrival over a desired field of view.
  • N is a positive integer that is less than the total number of elements in the receive antenna array.
  • DOAE 340 is not restricted to use an N-dimensional subarray but can actually employ the full array to achieve a finer definition of the directional distribution of incident signal power. As illustrated in FIG. 3, DOAE 340 receives sampled signals s 0 , s terme. . . s N , where s 0 is the reference port signal, and s ([s ⁇ . . . s N ] T is the N-component subarray signal vector.
  • C w + d 0 in some appropriate sense (e.g., least-mean-square sense, etc.).
  • C is the N x N complex subarray signal covariance matrix
  • d is the N-component subarray/reference signal cross-correlation vector, given respectively by
  • BCR technique is that it is numerically stable, even with limited arithmetic resolution, and the BCR technique can be applied iteratively to improve a current estimate for w, thereby allowing incremental refinements in a varying environment.
  • (x n , y n ) is the half-wavelength normalized xy-location of the n-th antenna element.
  • E c ( ⁇ , ⁇ ) is a composite pattern
  • E 0 ( ⁇ , ⁇ ) is a reference pattern
  • E n ( ⁇ , ⁇ ) is an embedded element pattern.
  • the reciprocal normalized composite power pattern exhibits an incident signal directional profile in the form of peaks and high regions relative to the reference antenna power pattern, as a function of the direction angles ⁇ , and ⁇ .
  • P( ⁇ , ⁇ ) is an estimate of incident directional power referenced to unity.
  • Knowledge of the normalized combined power pattern allows peak-finding techniques to be applied to identify the directions of incident signal activity within the field of view of a DOA estimation subsystem. Then, using knowledge of achievable beamwidths for DBF beams, it is possible to specify beam directions that provide optimal coverage, in some practical sense. Furthermore, with the knowledge of the current beam directions, the DOA estimation subsystem needs to only search for new beam directions and this technique eliminates unnecessary computations.
  • the computational complexity involved with this search technique depends on a number of factors including the step sizes used for ⁇ an ⁇ ⁇ . For example, the complexity is also dependent on the accuracy required for the directions of arrival and the size of the search range in the desired FOV.
  • incident signal energy having certain spectral characteristics can be used to determine where to induce nulls into the composite power pattern or a peak directional power estimate.
  • incident signal energy having certain spectral characteristics can be used to determine where to induce nulls into the composite power pattern or a peak directional power estimate.
  • incident signal is identified as an undesired signal (interference)
  • its DOA information is important to a DBF system.
  • the DBF system utilizes this knowledge to place a null at this point in each of the beams the DBF system forms.
  • a DBF system can form a beam in this direction in an attempt to further identify the undesired interference. By intelligently positioning nulls, undesired signals or jamming signals are suppressed using the DOA information.
  • digital beamforming refers to the creation of multiple, simultaneous, independently controlled beams at baseband, which are controlled through digital signal processing.
  • a DOA-aided DBF subsystem with an enhanced DOA estimation capability, allows a communication system to employ a resource allocation policy that assigns beams to sources and nulls to interferers to optimize total capacity and quality of service.
  • antenna patterns can be established to have beam and null positions that vary according to traffic density.
  • DOA-aided DBF subsystem 300 shown in FIG. 3 has advantages over conventional systems with fixed beam antennas because it, among other things, adaptively adjusts antenna patterns and produces accurately positioned beams and nulls in response to received data. Desirably, the received data is used to accurately estimate the DOA for desired and undesired incident signals. In other words, DOA- aided DBF subsystem 300 provides accurate antenna beam pointing in response to demand for communication sen/ices, and it provides improved null positioning to lessen the impact of unwanted RF signals. These features are implemented through appropriate software embedded in DBF 330, DOAE 340, and controller 350.
  • the DOA information must be updated as the satellites and users move relative to each other.
  • the DOA is estimated for signals incident to at least one antenna array on a satellite.
  • This information can be used, for example, to improve the allocation of resources by placing beams on active and high-traffic areas and by not placing beams on inactive regions. This improvement in resource allocation can enhance overall performance and reduce on-board power consumption.
  • FIG. 4 illustrates a flow diagram of a procedure for determining directions of arrival in a DOA-aided DBF subsystem in accordance with a preferred embodiment of the invention.
  • Procedure 400 starts in step 402.
  • a reference antenna is defined.
  • the reference antenna forms a reference beam that covers a desired field of view while rejecting the remaining space with a sufficiently low sidelobe structure.
  • the elements that are used to define the reference antenna are not hardware dependent.
  • the reference antenna is dynamically controlled by the DOAE.
  • the reference antenna varies as the size of the field of view changes and as the direction for the field of view changes. For example, in a non-geostationary satellite system fields of view change as the satellites move with respect to the surface of the earth. In addition, fields of view can be changed when beams are moved from light traffic areas to heavy traffic areas.
  • an N-element subarray is selected.
  • the N-element array is a subarray in the receive antenna array.
  • the N-element array is the entire receive array.
  • N changes as requirements change.
  • the invention uses a weighted subarray in combination with a reference antenna to determine the distribution of directions of arrival over a desired field of view.
  • the relaxation algorithm employed is numerically stable and robust imposing no requirement that the underlying covariance matrix have an inverse. Consequently, the system is not restricted to use a small subarray but can actually employ the full array to achieve a finer definition of the directional distribution of incident power.
  • the field of view can be changed to some extent when DOA estimates are required for a specific area.
  • step 410 the covariance matrix is computed.
  • a cross-correlation vector is computed.
  • step 412 the linear system of step 408 is solved for w, which is an N- dimensional complex weight vector over the N-element subarray.
  • w is an N- dimensional complex weight vector over the N-element subarray.
  • BCR or related eigenspace techniques are used.
  • a composite antenna pattern is computed using the previously obtained value for w.
  • the composite antenna pattern, E c is computed using the previously obtained value for w.
  • step 416 the reciprocal normalized composite power pattern P( ⁇ , ⁇ ) is computed for a plurality of points ⁇ .- with in a desired field of view.
  • the DOA estimates are derived for the plurality of incident signals from the reciprocal normalized composite power.pattern P( ⁇ , ⁇ computed at the plurality of points ( ⁇ , ⁇ ).
  • peak finding techniques are used to locate the peaks of P( ⁇ , ⁇ ) in ( ⁇ , ⁇ ) space corresponding to the DOA estimates.
  • procedure 400 ends.
  • Procedure 400 is used to calculate a composite power pattern for a plurality of points within a desired field of view, and to derive DOA estimates for at least one incident signal by applying at least one peak- finding technique to the composite power pattern.
  • the DOA-aided DBF subsystem establishes a beam in a direction substantially equal to the DOA estimate, when said incident signal is a desired signal; and the DOA-aided DBF subsystem directs a null in a direction substantially equal to the DOA estimate, when the incident signal is from an undesired source.
  • FIG. 5 shows a graph that illustrates the directional signal activity estimated by the DOA algorithm.
  • peaks 510 indicative of signal activity at three different directions of arrival, were produced in this example using DOA estimation techniques. Calculations were performed using BCR and pseudo-inverse techniques to solve for w for a multi-element linear array. In this example, the directions of arrival for the incident signals occurred at 49, 50, and 51 degrees.
  • a DOA-aided DBF subsystem includes both beam steering and null steering, the computing load increases and becomes more complex. Suitable techniques are used to accomplish this. For example, using the BLUE DBF algorithm in combination with DOA estimation techniques provides a good tradeoff between performance and complexity. This combination allows the simultaneous steering of beams in the directions of the desired users and the steering of nulls in the direction of interfering signals.
  • Computational complexity is also managed by controlling when calculations are performed. Desirably, the DOA estimation calculations are performed as required. For example, the update rate is dependent on a number of factors including, among other things, changes in satellite position relative to the desired and undesired transmitters and receivers. If the rate of change of the elevation angle is 0.0569°/ sec, then this corresponds to 0.005 frame, assuming a frame length of 90 msec. When a correction of the beamformer is required every 0.5°, then DOA information needs to be updated every 10Oth frame. This can be done using appropriate digital signal processors. In a DOA-aided DBF system, the computational complexity can be reduced by making iterative updates to the DOA estimates.
  • the BCR technique allows iterative w updates to be made, thereby providing an essentially continuous refinement of DOA estimates.
  • DOA computations can be distributed and time-shared with DBF computations. Desirably, these options allow the economic use of available computational resources.
  • peak detection can be made more computationally efficient by using iterative updating and searching near previously identified peaks.
  • the computational complexity can also be reduced by using approximate estimates of the DOA for at least some of the incident signals.
  • approximate estimates for directions of arrival are obtained for users during an acquisition phase. For example, wider FOVs and wider beams can be used during the acquisition phase.
  • DOA estimation techniques can be formulated for a single frequency of operation, presumably the center frequency of a communication system's operational bandwidth.
  • the bandwidth B is a small fraction of the carrier frequency f c
  • DOA estimates may be derived for individual subbands in a time-shared fashion.
  • a GPS-aided location system can also be used to obtain the directions of arrival.
  • location system data can be used to provide limits and initial conditions for BCR.
  • DOA estimation algorithms could be used to identify unknown sources of interference.
  • coarse location data can be provided to terrestrial based devices by satellites in a communications system.
  • the method and apparatus of the present invention enable the capabilities of a satellite communication system to be greatly enhanced by using antenna patterns with accurately positioned beams and nulls. Using accurately positioned nulls to minimize the effect of interfering signals can be optimized for various missions, and additional cost benefits can be accrued by the system using the method and apparatus of the invention. Furthermore, using the method and apparatus of the present invention with wide nulls adds to the overall system robustness.
  • the invention has been described above with reference to a preferred embodiment. However, those skilled in the art will recognize that changes and modifications can be made in this embodiment without departing from the scope of the invention. For example, while a preferred embodiment has been described in terms of using a specific block diagram for the transceiver, other systems can be envisioned which use different block diagrams. Also, although BCR is the preferred algorithm for updating the adaptive weight vector, other equivalent techniques can be applied. Accordingly, these and other changes and modifications which are obvious to those skilled in the art are intended to be included within the scope of the invention.

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Variable-Direction Aerials And Aerial Arrays (AREA)
  • Radio Transmission System (AREA)
PCT/US2000/007472 1999-05-03 2000-03-21 Robust estimation of doa for antenna arrays Ceased WO2000067042A1 (en)

Priority Applications (4)

Application Number Priority Date Filing Date Title
DE60020693T DE60020693T2 (de) 1999-05-03 2000-03-21 Robuste schätzung der empfangsrichtung für gruppenantennen
EP00919499A EP1177456B1 (en) 1999-05-03 2000-03-21 Robust estimation of doa for antenna arrays
AU40178/00A AU4017800A (en) 1999-05-03 2000-03-21 Robust estimation of doa for antenna arrays
JP2000615827A JP2002543436A (ja) 1999-05-03 2000-03-21 アンテナ・アレイ用のdoaの堅実な推定

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US09/305,347 1999-05-03
US09/305,347 US6075484A (en) 1999-05-03 1999-05-03 Method and apparatus for robust estimation of directions of arrival for antenna arrays

Publications (1)

Publication Number Publication Date
WO2000067042A1 true WO2000067042A1 (en) 2000-11-09

Family

ID=23180428

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/US2000/007472 Ceased WO2000067042A1 (en) 1999-05-03 2000-03-21 Robust estimation of doa for antenna arrays

Country Status (6)

Country Link
US (1) US6075484A (enExample)
EP (1) EP1177456B1 (enExample)
JP (1) JP2002543436A (enExample)
AU (1) AU4017800A (enExample)
DE (1) DE60020693T2 (enExample)
WO (1) WO2000067042A1 (enExample)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7427953B2 (en) 2004-10-15 2008-09-23 Interdigital Technology Corporation Wireless communication apparatus for determining direction of arrival information to form a three-dimensional beam used by a transceiver
KR102516903B1 (ko) * 2021-11-11 2023-04-03 한국전자통신연구원 통신 시스템에서 다수의 신호원들의 동시 전력 추정을 위한 방법 및 장치

Families Citing this family (74)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
SE509776C2 (sv) * 1997-07-04 1999-03-08 Ericsson Telefon Ab L M Anordning och förfarande vid antennlobsstyrning i radiokommunikationssystem
GB2337171A (en) * 1998-05-06 1999-11-10 Motorola Ltd Direction finder
JP3662772B2 (ja) * 1999-05-24 2005-06-22 東芝テック株式会社 無線通信システム
CN1107424C (zh) * 2000-06-12 2003-04-30 信息产业部电信科学技术研究院 在频分双工无线通信系统中使用智能天线的方法与装置
US7245880B1 (en) * 2000-08-31 2007-07-17 Intel Corporation Transmit power control within a wireless transmitter
JP2002084217A (ja) * 2000-09-08 2002-03-22 Matsushita Electric Ind Co Ltd 基地局装置および到来方向推定方法
US6606055B2 (en) * 2000-12-06 2003-08-12 Harris Corporation Phased array communication system providing airborne crosslink and satellite communication receive capability
CN100356711C (zh) * 2000-12-29 2007-12-19 诺基亚公司 基站、基站模块和估计上行链路信号的参数的方法
JP4569015B2 (ja) * 2001-02-28 2010-10-27 ソニー株式会社 広帯域アレイアンテナ
JP3767799B2 (ja) * 2001-04-09 2006-04-19 日本電気株式会社 アレーアンテナのヌル方向制御方法及び装置
DE60141741D1 (de) * 2001-04-27 2010-05-20 Mitsubishi Elec R&D Ct Europe Verfahren zur Einfallsrichtungsschätzung
US7327798B2 (en) 2001-10-19 2008-02-05 Lg Electronics Inc. Method and apparatus for transmitting/receiving signals in multiple-input multiple-output communication system provided with plurality of antenna elements
KR100463526B1 (ko) * 2002-01-04 2004-12-29 엘지전자 주식회사 다중 입력 다중 출력 시스템에서의 전력 할당 방법
US6738018B2 (en) * 2002-05-01 2004-05-18 Harris Corporation All digital phased array using space/time cascaded processing
US7084801B2 (en) * 2002-06-05 2006-08-01 Siemens Corporate Research, Inc. Apparatus and method for estimating the direction of arrival of a source signal using a microphone array
US7151951B2 (en) * 2002-12-23 2006-12-19 Telefonktiebolaget Lm Ericsson (Publ) Using beamforming and closed loop transmit diversity in a multi-beam antenna system
US7277730B2 (en) * 2002-12-26 2007-10-02 Nokia Corporation Method of allocating radio resources in telecommunication system, and telecommunication system
KR100770875B1 (ko) 2004-05-24 2007-10-26 삼성전자주식회사 배열 안테나 시스템에서 간섭전력 추정을 이용한 빔 형성장치 및 방법
US7382743B1 (en) * 2004-08-06 2008-06-03 Lockheed Martin Corporation Multiple-beam antenna system using hybrid frequency-reuse scheme
US8164533B1 (en) 2004-10-29 2012-04-24 Lockhead Martin Corporation Horn antenna and system for transmitting and/or receiving radio frequency signals in multiple frequency bands
JP2006258529A (ja) * 2005-03-16 2006-09-28 Fujitsu Ten Ltd 電波到来方向推定装置及び方法
US7239270B2 (en) 2005-05-31 2007-07-03 Research In Motion Limited Mobile wireless communications device comprising a satellite positioning system antenna and electrically conductive director element therefor
US7436370B2 (en) * 2005-10-14 2008-10-14 L-3 Communications Titan Corporation Device and method for polarization control for a phased array antenna
KR100749451B1 (ko) * 2005-12-02 2007-08-14 한국전자통신연구원 Ofdm 기지국 시스템에서의 스마트 안테나 빔 형성 방법및 장치
US7414578B1 (en) * 2006-09-27 2008-08-19 Rockwell Collins, Inc. Method for efficiently computing the beamforming weights for a large antenna array
US9733349B1 (en) 2007-09-06 2017-08-15 Rockwell Collins, Inc. System for and method of radar data processing for low visibility landing applications
US9354633B1 (en) 2008-10-31 2016-05-31 Rockwell Collins, Inc. System and method for ground navigation
US9939526B2 (en) 2007-09-06 2018-04-10 Rockwell Collins, Inc. Display system and method using weather radar sensing
FR2931302B1 (fr) * 2008-05-16 2010-07-30 Centre Nat Etd Spatiales Systeme antennaire multi-faisceaux pour couverture multispots et satellite comprenant un tel systeme
US7737904B2 (en) * 2008-06-11 2010-06-15 Lockheed Martin Corporation Antenna systems for multiple frequency bands
US8558731B1 (en) 2008-07-02 2013-10-15 Rockwell Collins, Inc. System for and method of sequential lobing using less than full aperture antenna techniques
US8077078B1 (en) 2008-07-25 2011-12-13 Rockwell Collins, Inc. System and method for aircraft altitude measurement using radar and known runway position
WO2010050269A1 (ja) * 2008-10-30 2010-05-06 三菱電機株式会社 通信装置および通信システム
US8170088B2 (en) * 2008-11-19 2012-05-01 Harris Corporation Methods for determining a reference signal at any location along a transmission media
US20100123618A1 (en) * 2008-11-19 2010-05-20 Harris Corporation Closed loop phase control between distant points
US7855681B2 (en) * 2008-11-19 2010-12-21 Harris Corporation Systems and methods for determining element phase center locations for an array of antenna elements
US7970365B2 (en) * 2008-11-19 2011-06-28 Harris Corporation Systems and methods for compensating for transmission phasing errors in a communications system using a receive signal
US7969358B2 (en) * 2008-11-19 2011-06-28 Harris Corporation Compensation of beamforming errors in a communications system having widely spaced antenna elements
US20100125347A1 (en) * 2008-11-19 2010-05-20 Harris Corporation Model-based system calibration for control systems
US20100124263A1 (en) * 2008-11-19 2010-05-20 Harris Corporation Systems for determining a reference signal at any location along a transmission media
US8804612B1 (en) * 2009-02-06 2014-08-12 Qualcomm Incorporated Triggering and transmitting sounding packets for wireless communications
KR101097843B1 (ko) 2010-12-08 2011-12-23 엘아이지넥스원 주식회사 위상각 검출 장치 및 그 방법
US9019145B1 (en) 2011-07-14 2015-04-28 Rockwell Collins, Inc. Ground clutter rejection for weather radar
CN102385048A (zh) * 2011-08-10 2012-03-21 西安交通大学 基于均匀线阵的混合信号方向估计方法
US8947294B1 (en) * 2011-11-02 2015-02-03 Lockheed Martin Corporation Method and system for adaptively cancelling clutter from the sidelobes of a ground-based radar
ES2596367T3 (es) * 2012-06-01 2017-01-09 Airbus Ds Gmbh Método de detección de una dirección de llegada de al menos una señal de interferencia y sistema para realizar dicho método
US9262932B1 (en) 2013-04-05 2016-02-16 Rockwell Collins, Inc. Extended runway centerline systems and methods
JP2015012568A (ja) * 2013-07-02 2015-01-19 三星電子株式会社Samsung Electronics Co.,Ltd. 指向性制御装置、および指向性制御方法
KR102177553B1 (ko) * 2014-03-27 2020-11-11 삼성전자주식회사 다중 사용자 지원을 위한 빔포밍 방법 및 장치
US10270524B2 (en) 2014-04-15 2019-04-23 Space Systems/Loral, Llc Broadband satellite payload architecture
CN104298850A (zh) * 2014-07-18 2015-01-21 哈尔滨工业大学深圳研究生院 信源数未知的相干信号测向方法及系统
CN104181529A (zh) * 2014-07-21 2014-12-03 中国科学院电子学研究所 一种Ka波段合成孔径雷达SAR信号处理方法和设备
US10928510B1 (en) 2014-09-10 2021-02-23 Rockwell Collins, Inc. System for and method of image processing for low visibility landing applications
US20160301463A1 (en) * 2015-04-10 2016-10-13 Space Systems/Loral, Llc Broadband satellite payload architecture
JP6275330B2 (ja) * 2015-04-30 2018-02-07 三菱電機株式会社 送信局、制御局、受信局、データ伝送システムおよびデータ伝送方法
US10705201B1 (en) 2015-08-31 2020-07-07 Rockwell Collins, Inc. Radar beam sharpening system and method
US10700762B2 (en) * 2016-05-04 2020-06-30 Telefonaktiebolaget Lm Ericsson (Publ) Beam forming using an antenna arrangement
US10228460B1 (en) 2016-05-26 2019-03-12 Rockwell Collins, Inc. Weather radar enabled low visibility operation system and method
US10353068B1 (en) 2016-07-28 2019-07-16 Rockwell Collins, Inc. Weather radar enabled offshore operation system and method
US10288715B2 (en) * 2016-09-09 2019-05-14 Raytheon Company Systems and methods for direction finding using augmented spatial sample covariance matrices
CN106788648B (zh) * 2016-11-30 2020-04-24 西南交通大学 一种智能天线系统的自适应波束形成方法
EP3616332A1 (en) 2017-04-28 2020-03-04 Fraunhofer Gesellschaft zur Förderung der Angewand Specular component estimation in a wireless communication network
JP6926702B2 (ja) * 2017-06-09 2021-08-25 富士通株式会社 無線通信装置及びビーム制御方法
CN110426670B (zh) * 2018-12-26 2022-09-23 西安电子科技大学 基于tls-cs的外辐射源雷达超分辨doa估计方法
CN110380770B (zh) * 2019-06-10 2021-07-06 浙江大学 一种低轨移动卫星通信网络的自适应对星方法
US11316268B2 (en) * 2019-06-24 2022-04-26 Rockwell Collins, Inc. Method for antenna beam and null steering under high platform dynamics
KR102812754B1 (ko) 2019-10-15 2025-05-23 삼성전자주식회사 통신 장치 및 통신 장치의 데이터 수신 방법
WO2022011510A1 (zh) * 2020-07-13 2022-01-20 北京小米移动软件有限公司 一种通信方法、通信装置及存储介质
KR102499336B1 (ko) * 2021-07-08 2023-02-14 국방과학연구소 적응형 널 폭 제어 장치와 방법 및 이를 이용한 기지국
US11929798B2 (en) * 2021-08-30 2024-03-12 Rockwell Collins, Inc. Technique for post-correlation beamforming
SE545126C2 (en) * 2021-08-31 2023-04-11 Saab Ab Multi-channel active array system and Method for obtaining positional information of an object
US20230106766A1 (en) * 2021-10-05 2023-04-06 Qualcomm Incorporated Nulling for inter-user equipment interference cancellation
US20240393427A1 (en) * 2023-05-23 2024-11-28 Rockwell Collins, Inc. Aesa-based synthetic nulling for enhanced radar ground clutter suppression
CN119619985A (zh) * 2024-12-03 2025-03-14 中国空间技术研究院 通信卫星干扰源定位方法及系统

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4353119A (en) * 1980-06-13 1982-10-05 Motorola Inc. Adaptive antenna array including batch covariance relaxation apparatus and method
US4720712A (en) * 1985-08-12 1988-01-19 Raytheon Company Adaptive beam forming apparatus
US5754138A (en) * 1996-10-30 1998-05-19 Motorola, Inc. Method and intelligent digital beam forming system for interference mitigation

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4771289A (en) * 1982-05-28 1988-09-13 Hazeltine Corporation Beamforming/null-steering adaptive array
US5144322A (en) * 1988-11-25 1992-09-01 The United States Of America As Represented By The Secretary Of The Navy Large-aperture sparse array detector system for multiple emitter location

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4353119A (en) * 1980-06-13 1982-10-05 Motorola Inc. Adaptive antenna array including batch covariance relaxation apparatus and method
US4720712A (en) * 1985-08-12 1988-01-19 Raytheon Company Adaptive beam forming apparatus
US5754138A (en) * 1996-10-30 1998-05-19 Motorola, Inc. Method and intelligent digital beam forming system for interference mitigation

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
BARRET M ET AL: "ADAPTIVE ANTENNAS FOR MOBILE COMMUNICATIONS", ELECTRONICS AND COMMUNICATION ENGINEERING JOURNAL,GB,INSTITUTION OF ELECTRICAL ENGINEERS, LONDON, vol. 6, no. 4, 1 August 1994 (1994-08-01), pages 203 - 214, XP000469556, ISSN: 0954-0695 *
VEEN VAN DER A -J ET AL: "SUBSPACE-BASED SIGNAL ANALYSIS USING SINGULAR VALUE DECOMPOSITION", PROCEEDINGS OF THE IEEE,US,IEEE. NEW YORK, vol. 81, no. 9, 1 September 1993 (1993-09-01), pages 1277 - 1308, XP000408616, ISSN: 0018-9219 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7427953B2 (en) 2004-10-15 2008-09-23 Interdigital Technology Corporation Wireless communication apparatus for determining direction of arrival information to form a three-dimensional beam used by a transceiver
US7705779B2 (en) 2004-10-15 2010-04-27 Interdigital Technology Corporation Wireless communication apparatus for determining direction of arrival information to form a three-dimensional beam used by a transceiver
KR102516903B1 (ko) * 2021-11-11 2023-04-03 한국전자통신연구원 통신 시스템에서 다수의 신호원들의 동시 전력 추정을 위한 방법 및 장치

Also Published As

Publication number Publication date
US6075484A (en) 2000-06-13
AU4017800A (en) 2000-11-17
EP1177456B1 (en) 2005-06-08
DE60020693T2 (de) 2005-10-20
EP1177456A1 (en) 2002-02-06
JP2002543436A (ja) 2002-12-17
DE60020693D1 (de) 2005-07-14

Similar Documents

Publication Publication Date Title
US6075484A (en) Method and apparatus for robust estimation of directions of arrival for antenna arrays
EP0956612B1 (en) Antenna beam patterns having wide nulls
US6697633B1 (en) Method permitting increased frequency re-use in a communication network, by recovery of transmitted information from multiple cochannel signals
US6215983B1 (en) Method and apparatus for complex phase equalization for use in a communication system
US6018317A (en) Cochannel signal processing system
US6208295B1 (en) Method for processing radio signals that are subject to unwanted change during propagation
US6535666B1 (en) Method and apparatus for separating signals transmitted over a waveguide
US6658234B1 (en) Method for extending the effective dynamic range of a radio receiver system
US6392596B1 (en) Single-port weighting systems for GPS reception in multiple-interference environments
CA2272930C (en) Cochannel signal processing system
CN115396005B (zh) 多波束卫星的波束间干扰及用户信道向量确定方法及装置
JPH10145275A (ja) トラフィック要求に応答するインテリジェント・ビーム形成方法およびシステム
Huang et al. Frequency-domain AoA estimation and beamforming with wideband hybrid arrays
US7068219B2 (en) Communications system including phased array antenna providing nulling and related methods
US5977907A (en) Method and system for antenna pattern synthesis based on geographical distribution of subscribers
US7414578B1 (en) Method for efficiently computing the beamforming weights for a large antenna array
EP1894320B1 (en) Communications systems including adaptive antenna systems and methods for inter-system and intra-system interference reduction
AU727787B2 (en) Cochannel signal processing system
AU763322B2 (en) A method for increasing the capacity of a two-way communication system by permitting cochannel signals
Alio et al. Adaptive array antennas for mobile earth stations: A review
AU6654200A (en) Cochannel signal processing system
Senapati A thesis submitted To
Xu Array signal processing for communications from high altitude platforms and other applications

Legal Events

Date Code Title Description
AK Designated states

Kind code of ref document: A1

Designated state(s): AE AL AM AT AU AZ BA BB BG BR BY CA CH CN CU CZ DE DK EE ES FI GB GD GE GH GM HR HU ID IL IN IS JP KE KG KR KZ LC LR LS LT LU LV MD MG MK MN MW MX NO NZ PL PT RO RU SD SE SG SI SK SL TJ TM TR TT UA UG UZ VN YU ZA ZW

AL Designated countries for regional patents

Kind code of ref document: A1

Designated state(s): GH GM KE LS MW SD SL SZ TZ UG ZW AM AZ BY KG KZ MD RU TJ TM AT BE CH CY DE DK ES FI FR GB GR IE IT LU MC NL PT SE BF BJ CF CG CI CM GA GN GW ML MR NE SN TD TG

121 Ep: the epo has been informed by wipo that ep was designated in this application
DFPE Request for preliminary examination filed prior to expiration of 19th month from priority date (pct application filed before 20040101)
WWE Wipo information: entry into national phase

Ref document number: 2000919499

Country of ref document: EP

ENP Entry into the national phase

Ref country code: JP

Ref document number: 2000 615827

Kind code of ref document: A

Format of ref document f/p: F

WWP Wipo information: published in national office

Ref document number: 2000919499

Country of ref document: EP

REG Reference to national code

Ref country code: DE

Ref legal event code: 8642

WWG Wipo information: grant in national office

Ref document number: 2000919499

Country of ref document: EP