US12480760B2 - Angle width estimation device, angle width estimation method, and target angle measurement device - Google Patents
Angle width estimation device, angle width estimation method, and target angle measurement deviceInfo
- Publication number
- US12480760B2 US12480760B2 US18/207,235 US202318207235A US12480760B2 US 12480760 B2 US12480760 B2 US 12480760B2 US 202318207235 A US202318207235 A US 202318207235A US 12480760 B2 US12480760 B2 US 12480760B2
- Authority
- US
- United States
- Prior art keywords
- null
- width
- angle
- unit
- beams
- 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.)
- Active, expires
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO 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/00—Direction-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/02—Direction-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/14—Systems for determining direction or deviation from predetermined direction
- G01S3/46—Systems for determining direction or deviation from predetermined direction using antennas spaced apart and measuring phase or time difference between signals therefrom, i.e. path-difference systems
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B11/00—Measuring arrangements characterised by the use of optical techniques
- G01B11/26—Measuring arrangements characterised by the use of optical techniques for measuring angles or tapers; for testing the alignment of axes
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO 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/00—Direction-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/02—Direction-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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO 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/00—Direction-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/02—Direction-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/04—Details
- G01S3/043—Receivers
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B17/00—Monitoring; Testing
- H04B17/20—Monitoring; Testing of receivers
- H04B17/27—Monitoring; Testing of receivers for locating or positioning the transmitter
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B17/00—Monitoring; Testing
- H04B17/30—Monitoring; Testing of propagation channels
- H04B17/309—Measuring or estimating channel quality parameters
- H04B17/318—Received signal strength
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B7/00—Radio transmission systems, i.e. using radiation field
- H04B7/02—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
- H04B7/04—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas
- H04B7/08—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the receiving station
- H04B7/0837—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the receiving station using pre-detection combining
- H04B7/0842—Weighted combining
- H04B7/086—Weighted combining using weights depending on external parameters, e.g. direction of arrival [DOA], predetermined weights or beamforming
Definitions
- a target to be observed is a distant target that can be regarded as a point wave source
- a reflected wave from the target does not have an angle width.
- the target to be observed is an adjacent target that cannot be regarded as the point wave source
- the reflected wave from the target has the angle width.
- Non-Patent Literature 1 discloses an angle width estimating technology of estimating an angle width of a reflected wave from an adjacent target.
- a computer estimates the angle width of the reflected wave by executing the Capon method.
- the Capon method is a method of estimating the angle width that can be executed when a plurality of reception array signals at different sampling times is acquired by an array antenna.
- Non-Patent Literature 1 the computer cannot execute the Capon method unless a plurality of reception array signals is acquired. For this reason, there is a problem that the computer cannot estimate the angle width of the reflected wave unless a plurality of reception array signals is acquired.
- the present disclosure has been achieved to solve the above-described problem, and an object thereof is to acquire an angle width estimation device and an angle width estimation method capable of estimating an angle width of a reflected wave from a target from one reception array signal.
- An angle width estimation device includes: beam forming circuitry to acquire one reception array signal of a reflected wave from a target to be observed, and form, from the one reception array signal, a plurality of null beams having nulls in an arrival direction of the reflected wave, and having null widths which are widths of the nulls and different from each other; and angle width estimating circuitry to compare powers of the plurality of null beams formed by the beam forming circuitry with each other, and estimate a null width indicating an angle width of the reflected wave on the basis of a comparison result of the powers.
- FIG. 1 is a configuration diagram illustrating a target angle measurement device including an angle width estimation device 4 according to a first embodiment.
- FIG. 2 is a hardware configuration diagram illustrating hardware of a digital signal processing unit in the target angle measurement device illustrated in FIG. 1 .
- FIG. 3 is a hardware configuration diagram of a computer when the angle width estimation device 4 is implemented by software, firmware or the like.
- FIG. 4 is a flowchart illustrating an angle width estimation method being a processing procedure of the angle width estimation device 4 according to the first embodiment.
- FIG. 5 is an explanatory diagram illustrating each of an elevation angle ⁇ and an azimuth angle ⁇ in a steering vector a( ⁇ , ⁇ ).
- FIG. 6 is an explanatory diagram illustrating a null width in a beam forming direction u d in a null beam.
- FIG. 7 is an explanatory diagram illustrating a relationship between the null width of the null beam and a target.
- FIG. 8 is an explanatory diagram illustrating a relationship between the null width of the null beam and power of the null beam.
- FIG. 9 is a configuration diagram illustrating a target angle measurement device including an angle width estimation device 4 according to a second embodiment.
- FIG. 10 is a hardware configuration diagram illustrating hardware of a digital signal processing unit in the target angle measurement device illustrated in FIG. 9 .
- FIG. 11 is an explanatory diagram illustrating a null beam formed by a DBF unit 12 .
- FIG. 12 is a configuration diagram illustrating a target angle measurement device including an angle width estimation device 4 according to a third embodiment.
- FIG. 13 is a hardware configuration diagram illustrating hardware of a digital signal processing unit in the target angle measurement device illustrated in FIG. 12 .
- FIG. 14 is an explanatory diagram illustrating rotation of a forming direction of the null beam.
- FIG. 1 is a configuration diagram illustrating a target angle measurement device including an angle width estimation device 4 according to a first embodiment.
- FIG. 2 is a hardware configuration diagram illustrating hardware of a digital signal processing unit in the target angle measurement device illustrated in FIG. 1 .
- the target angle measurement device illustrated in FIG. 1 is provided with a reflected wave receiving unit 1 , the angle width estimation device 4 , and a target angle measuring unit 10 .
- the reflected wave receiving unit 1 is provided with reception antenna elements 2 - 1 to 2 -M and analog-to-digital converters (hereinafter referred to as “A/D converters”) 3 - 1 to 3 -M.
- M is an integer equal to or larger than two.
- the reflected wave receiving unit 1 receives a reflected wave from a target to be observed, and outputs a reception array signal r ARE of the reflected wave to the angle width estimation device 4 .
- the reception array signal r ARE of the reflected wave is a signal including a plurality of reception signals r 1 to r M output from the reception antenna elements 2 - 1 to 2 -M, respectively.
- the target angle measurement device illustrated in FIG. 1 is not provided with a radio wave transmitting unit that emits a radio wave toward the target to be observed. Therefore, in the target angle measurement device illustrated in FIG. 1 , it is supposed that an external radio wave transmitting unit not illustrated emits the radio wave toward the target to be observed. However, this is merely an example, and the target angle measurement device illustrated in FIG. 1 may be provided with the radio wave transmitting unit.
- the reception antenna elements 2 - 1 to 2 -M form a reception array antenna.
- the reception antenna element 2 - m receives the reflected wave from the target to be observed, and outputs the reception signal r m included in the reception array signal r ARE of the reflected wave to the A/D converter 3 - m.
- the A/D converter 3 - m converts the reception signal r m output from the reception antenna element 2 - m from an analog signal to a digital signal.
- the A/D converter 3 - m outputs the digital signal to the angle width estimation device 4 , as the reception signal included in the reception array signal r ARE .
- a receiver that detects the reception signal r m output from the reception antenna element 2 - m is not illustrated.
- a receiver may be mounted between the reception antenna element 2 - m and the A/D converter 3 - m , or the receiver may be included in the A/D converter 3 - m.
- the reflected wave receiving unit 1 is provided with the reception array antenna including the reception antenna elements 2 - 1 to 2 -M.
- the reflected wave receiving unit 1 may be provided with an array microphone including a plurality of microphones or an array sensor including a plurality of radio wave sensors, in place of the reception array antenna.
- the reflected wave receiving unit 1 is provided with the array microphone, the reflected wave is a sound wave.
- the reflected wave receiving unit 1 is provided with the array sensor, the reflected wave is a radio wave.
- the angle width estimation device 4 is provided with a beam forming unit 5 and an angle width estimating unit 9 .
- the angle width estimation device 4 estimates an angle width of the reflected wave, by performing signal processing on the reception array signal r ARE output from the reflected wave receiving unit 1 .
- the beam forming unit 5 is provided with coherent integrating units 6 - 1 to 6 -M, pulse compressing units 7 - 1 to 7 -M, and a digital beam forming (DBF) unit 8 .
- coherent integrating units 6 - 1 to 6 -M coherent integrating units 6 - 1 to 6 -M
- pulse compressing units 7 - 1 to 7 -M pulse compressing units 7 - 1 to 7 -M
- DBF digital beam forming
- the beam forming unit 5 acquires one reception array signal r ARE output from the reflected wave receiving unit 1 , and forms, from the one reception array signal r ARE , a plurality of null beams having nulls in an arrival direction of the reflected wave, and having null widths which are widths of the nulls and different from each other.
- the DBF unit 8 is implemented by a DBF circuit 23 illustrated in FIG. 2 , for example.
- the DBF unit 8 is provided with a table 8 a.
- the table 8 a stores a plurality of weights regarding the null beam formed by the DBF unit 8 . Because of the different weights, one or more of a beam forming direction, a null forming direction, and the null width of the null beam changes. Each weight includes M weight elements.
- the DBF unit 8 acquires M signals after the pulse compression from the pulse compressing units 7 - 1 to 7 -M.
- the DBF unit 8 forms the null beam having the null in the arrival direction of the reflected wave, by performing digital beam forming on the M signals after the pulse compression output from the pulse compressing units 7 - 1 to 7 -M.
- the DBF unit 8 multiplies the signals after the pulse compression by the respective weight elements included in the acquired weight, and calculates the sum of a plurality of signals after the weight element multiplication, thereby forming the null beam having the null in the arrival direction of the reflected wave.
- the DBF unit 8 outputs the null beam to the angle width estimating unit 9 .
- the DBF unit 8 After forming the null beam, when the DBF unit 8 acquires the null beam forming command from the angle width estimating unit 9 , the DBF unit 8 acquires one weight different from the already acquired weight from the table 8 a.
- the angle width estimating unit 9 monitors a change in power of the null beam output from the DBF unit 8 , and searches for, as the null width indicating the angle width, the null width of the null beam when the power stops decreasing from a state in which the power decreases.
- the angle width estimating unit 9 outputs the angle width indicated by the estimated null width to the target angle measuring unit 10 .
- Each of the coherent integrating circuit 21 , the pulse compressing circuit 22 , the DBF circuit 23 , and the angle width estimating circuit 24 corresponds to, for example, a single circuit, a composite circuit, a programmed processor, a parallel-programmed processor, an application specific integrated circuit (ASIC), a field-programmable gate array (FPGA), or a combination thereof.
- ASIC application specific integrated circuit
- FPGA field-programmable gate array
- a program for causing the computer to execute each processing procedure in the coherent integrating units 6 - 1 to 6 -M, the pulse compressing units 7 - 1 to 7 -M, the DBF unit 8 , and the angle width estimating unit 9 is stored in a memory 31 . Then, a processor 32 of the computer executes the program stored in the memory 31 .
- the A/D converter 3 - m converts the reception signal r m from the analog signal to the digital signal and outputs the digital signal to the angle width estimation device 4 .
- FIG. 4 is a flowchart illustrating an angle width estimation method being the processing procedure of the angle width estimation device 4 according to the first embodiment.
- a reception array signal r ARE (t) at time t in the reception array antenna with the number of channels M is expressed by the following Expression (1).
- r ARE ( t ) a ( ⁇ , ⁇ ) s ( t )+ n ( t ) (1)
- a( ⁇ , ⁇ ) represents a steering vector when an elevation angle is ⁇ and an azimuth angle is ⁇ .
- s(t) represents a signal complex amplitude at time t
- n(t) represents a thermal noise vector at time t.
- FIG. 5 is an explanatory diagram illustrating each of the elevation angle ⁇ and the azimuth angle ⁇ in the steering vector a( ⁇ , ⁇ ).
- a ⁇ ( ⁇ , ⁇ ) exp ⁇ ( j ⁇ 2 ⁇ ⁇ ⁇ ⁇ R ⁇ L ⁇ ( ⁇ , ⁇ ) ) ( 2 )
- R [ x 1 y 1 ⁇ ⁇ x M y M ] ( 3 )
- L ⁇ ( ⁇ , ⁇ ) [ cos ⁇ ⁇ ⁇ sin ⁇ ⁇ ⁇ sin ⁇ ⁇ ] T ( 4 )
- ⁇ represents a wavelength
- R represents an element coordinate matrix in the reception antenna elements 2 - 1 to 2 -M
- x m represents an x coordinate of a position at which an m-th reception antenna element 2 - m is installed
- y m represents a y coordinate of the position at which the m-th reception antenna element 2 - m is installed.
- L( ⁇ , ⁇ ) represents a line-of-sight direction vector when the elevation angle is ⁇ and the azimuth angle is ⁇ .
- a weight w for forming the null beam having the null with respect to u,v is acquired by the following Expression (5) on the basis of a directionally constrained minimization of power (DCMP) method when beam forming directions u d and v d are determined.
- DCMP directionally constrained minimization of power
- ⁇ represents any normalization coefficient
- R xx ⁇ 1 represents an inverse matrix of a correlation matrix R xx of the reflected wave having the angle width.
- the null width of the null beam in which the null forming directions are u c and v c is defined by the correlation matrix R xx of the reflected wave having the angle width.
- FIG. 6 is an explanatory diagram illustrating the null width in the beam forming direction u d in the null beam.
- the beam forming direction u d in a u-axis direction coincides with the null forming direction u c in the u-axis direction.
- the width in the u-axis direction of the null beam in which a gain of the null beam is equal to or smaller than a threshold Th coincides with the null width in the u-axis direction of the null beam.
- a threshold Th for example, a gain lower, by a signal to noise ratio (SNR), than a gain Gr of the beam having no null in the beam forming direction u d is used, as expressed by the following Expression (6).
- the null width in a v-axis direction in the null beam is handled similarly to the null width in the u-axis direction in the null beam.
- the beam forming direction v d in the v-axis direction coincides with the null forming direction v c in the v-axis direction.
- the width in the v-axis direction of the null beam in which a gain of the null beam is equal to or smaller than a threshold Th coincides with the null width in the v-axis direction of the null beam.
- FIG. 7 is an explanatory diagram illustrating a relationship between the null width of the null beam and the target.
- FIG. 7 illustrates an example in which the null widths in the u-axis direction and the v-axis direction of the null beam are wider than the widths in the u-axis direction and the v-axis direction of the target, so that the gains in the null forming directions u c and v c of the null beam are equal to or smaller than the threshold Th.
- R xx P d ⁇ ⁇ u c + ⁇ ⁇ u m ⁇ i ⁇ n u c + ⁇ ⁇ u max ⁇ v c + ⁇ ⁇ v m ⁇ i ⁇ n v c + ⁇ ⁇ v m ⁇ ax ⁇ ⁇ ( u , v ) ⁇ a ⁇ ( u , v ) ⁇ a ⁇ ( u , v ) H ⁇ dvdu + ⁇ ⁇ I ( 7 )
- P d represents power of the reflected wave received by the reception antenna element 2 - m .
- ⁇ u max represents a constant indicating a maximum angle of the null beam in the u-axis direction
- ⁇ u min represents a constant indicating a minimum angle of the null beam in the u-axis direction.
- ⁇ v max represents a constant indicating a maximum angle of the null beam in the v-axis direction
- ⁇ v min represents a constant indicating a minimum angle of the null beam in the v-axis direction.
- ⁇ (u,v) represents a spatial distribution of desired wave power
- ⁇ represents pseudo noise
- I represents a unit matrix
- a weight w( ⁇ Beam , ⁇ Beam , ⁇ NULL , ⁇ NULL ) capable of forming the null beam having the null in the null forming directions u c and v c can be calculated.
- ⁇ Beam represents a center direction of the beam forming direction u d
- ⁇ Beam represents a center direction of the beam forming direction v d
- ⁇ NULL represents the null forming direction u c
- ⁇ NULL represents the null forming direction v c .
- the null width can be changed by changing each value of ⁇ u max , ⁇ u min , ⁇ v max , and ⁇ v min .
- the coherent integrating unit 6 - m outputs the signal after the coherent integration to the pulse compressing unit 7 - m.
- the pulse compressing unit 7 - m When receiving the signal after the coherent integration from the coherent integrating unit 6 - m , the pulse compressing unit 7 - m performs pulse compression on the signal after the coherent integration in order to suppress an unnecessary wave received by the reception antenna element 2 - m (step ST 2 in FIG. 4 ).
- the pulse compressing unit 7 - m outputs a signal r APC,m (t) after the pulse compression to the DBF unit 8 .
- the table 8 a of the DBF unit 8 stores the weight w( ⁇ Beam , ⁇ Beam , ⁇ NULL , ⁇ NULL ) calculated by substituting the inverse matrix of the correlation matrix R xx in Expression (7) into Expression (5).
- Each of the DBF unit 8 and the angle width estimating unit 9 repeatedly performs the following processing a plurality of times.
- the following processing can be repeated a plurality of times as long as one reception array signal r ARE (t) is acquired.
- null width u c ⁇ u/2 in the u-axis direction is described supposing that the null width v c ⁇ v/2 in the v-axis direction is constant. It is supposed that the null width in the v-axis direction is sufficiently wider than the width in the v-axis direction of the target. Further, it is supposed that the null forming directions u c and v c are constant.
- the DBF unit 8 acquires a weight w 1 ( ⁇ Beam , ⁇ Beam , ⁇ NULL , ⁇ NULL ) for forming the null with the null width of 0 in the u-axis direction or the null with a narrow null width in the u-axis direction from the table 8 a .
- the null with the narrow null width is the null with the null width close to 0.
- the DBF unit 8 multiplies the signal r APC,m (t) after the pulse compression output from the pulse compressing unit 7 - m , by a weight element w 1,m corresponding to the m-th pulse compressing unit 7 - m out of M weight elements w 1,1 to w 1,M included in the weight w 1 ( ⁇ Beam , ⁇ Beam , ⁇ NULL , ⁇ NULL ).
- the DBF unit 8 forms a null beam B NULL,1 by calculating the sum ⁇ r APC,m (t) ⁇ w 1,m of the M multiplication results (step ST 3 in FIG. 4 ).
- the DBF unit 8 outputs the null beam B NULL,1 to the angle width estimating unit 9 .
- the angle width estimating unit 9 calculates power P 1 ( ⁇ Beam , ⁇ Beam , ⁇ NULL , ⁇ NULL ) of the null beam B NULL,1 , as expressed by the following Expression (8) (step ST 4 in FIG. 4 ).
- P 1 ( ⁇ Beam , ⁇ Beam , ⁇ Null , ⁇ null )
- the angle width estimating unit 9 stores the power P 1 ( ⁇ Beam , ⁇ Beam , ⁇ NULL , ⁇ NULL ) of the null beam B NULL,1 in an internal memory.
- the DBF unit 8 acquires a weight w j ( ⁇ Beam , ⁇ Beam , ⁇ NULL , ⁇ NULL ) for forming the null with a wider null width in the u-axis direction than that in (j ⁇ 1)-th processing from the table 8 a.
- the DBF unit 8 multiplies the signal r APC,m (t) after the pulse compression output from the pulse compressing unit 7 - m , by a weight element w j,m corresponding to the m-th pulse compressing unit 7 - m out of the M weight elements w j,1 to w j,M included in the weight w 1 ( ⁇ Beam , ⁇ Beam , ⁇ NULL , ⁇ NULL ).
- the DBF unit 8 forms a null beam B NULL,j by calculating the sum ⁇ r APC,m (t) ⁇ w j,m of M multiplication results (step ST 3 in FIG. 4 ).
- the DBF unit 8 outputs the null beam B NULL,j to the angle width estimating unit 9 .
- the angle width estimating unit 9 calculates power P j ( ⁇ Beam , ⁇ Beam , ⁇ NULL , ⁇ NULL ) of the null beam B NULL,j , as expressed by the following Expression (9) (step ST 4 in FIG. 4 ).
- P j ( ⁇ Beam , ⁇ Beam , ⁇ Null , ⁇ null )
- the angle width estimating unit 9 stores the power P j ( ⁇ Beam , ⁇ Beam , ⁇ NULL , ⁇ NULL ) of the null beam B NULL,j in the internal memory.
- Each of the DBF unit 8 and the angle width estimating unit 9 repeatedly performs the j-th processing until the power P j ( ⁇ Beam , ⁇ Beam , ⁇ NULL , ⁇ NULL ) of the null beam B NULL,j stops decreasing even when the null width in the u-axis direction is widened.
- FIG. 8 is an explanatory diagram illustrating a relationship between the null width of the null beam and the power of the null beam.
- the radio waves reflected by the target increase, so that the power of the null beam increases.
- the radio waves reflected by the target decrease, so that the power of the null beam decreases. Note that, for example, when the null width of the null beam is wider than the width of the target as illustrated in FIG. 7 , even when the null width becomes further wider, the radio wave reflected by the target does not change, so that the power of the null beam does not decrease. That is, the power of the null beam is not smaller than known receiver noise power.
- the processing is performed four times, and the power P j ( ⁇ Beam , ⁇ Beam , ⁇ NULL , ⁇ NULL ) stops decreasing in the third processing.
- the null width in the third processing indicates the angle width of the reflected wave.
- the angle width estimating unit 9 compares J pieces of power P 1 ( ⁇ Beam , ⁇ Beam , ⁇ NULL , ⁇ NULL ) to P J ( ⁇ Beam , ⁇ Beam , ⁇ NULL , ⁇ NULL ) stored in the internal memory with one another.
- the angle width estimating unit 9 searches for a null width B ⁇ in the u-axis direction when the power P j ( ⁇ Beam , ⁇ Beam , ⁇ NULL , ⁇ NULL ) stops decreasing from a state in which the power P j ( ⁇ Beam , ⁇ Beam , ⁇ NULL , ⁇ NULL ) decreases, on the basis of a comparison result of the powers (step ST 5 in FIG. 4 ).
- the DBF unit 8 changes the null width u c ⁇ u/2 in the u-axis direction supposing that the null width v c ⁇ v/2 in the v-axis direction is constant.
- the DBF unit 8 changes the null width v c ⁇ v/2 in the v-axis direction supposing that the null width u c ⁇ u/2 in the u-axis direction is constant. It is supposed that the null width in the u-axis direction is sufficiently wider than the width in the u-axis direction of the target. Further, it is supposed that the null forming directions u c and v c are constant.
- Each of the DBF unit 8 and the angle width estimating unit 9 searches for a null width B ⁇ in the v-axis direction by a method similar to that of the processing of searching for the null width B ⁇ in the u-axis direction while changing the null width v c ⁇ v/2 in the v-axis direction.
- the angle width estimating unit 9 outputs, as an angle width A ⁇ in the u-axis direction of the reflected wave, the searched null width B ⁇ in the u-axis direction to the target angle measuring unit 10 , and outputs, as an angle width A ⁇ in the v-axis direction of the reflected wave, the searched null width B ⁇ in the v-axis direction to the target angle measuring unit 10 .
- the target angle measuring unit 10 measures the angle of the target using the angle widths A ⁇ and A ⁇ output from the angle width estimating unit 9 .
- the target angle measuring unit 10 outputs an angle measurement value of the target to an external display or the like.
- angle measurement processing of the target by the target angle measuring unit 10 is a known technology, the detailed description thereof is omitted. Hereinafter, an example of the angle measurement processing of the target is briefly described.
- the target angle measuring unit 10 acquires a weight w j ( ⁇ Beam , ⁇ Beam , ⁇ NULL , ⁇ NULL ) corresponding to both the angle width A ⁇ in the u-axis direction and the angle width A ⁇ in the v-axis direction from the table 8 a of the DBF unit 8 .
- the target angle measuring unit 10 acquires ⁇ Beam which is a parameter of the acquired weight w j ( ⁇ Beam , ⁇ Beam , ⁇ NULL , ⁇ NULL ).
- the target angle measuring unit 10 also acquires ⁇ Beam which is a parameter of the acquired weight w j ( ⁇ Beam , ⁇ Beam , ⁇ NULL , ⁇ NULL ).
- the angle measurement value in the u-axis direction of the target is ⁇ Beam
- the angle measurement value in the v-axis direction of the target is ⁇ Beam
- the target angle measuring unit 10 provides a latitude in the angle measurement value in the u-axis direction of the target, depending on the angle width A ⁇ in the u-axis direction of the reflected wave. That is, the target angle measuring unit 10 calculates ⁇ Beam ⁇ A ⁇ /2 as the angle measurement value in the u-axis direction of the target.
- the target angle measuring unit 10 Since the reflected wave has the angle width A ⁇ in the v-axis direction, the angle measurement value in the v-axis direction of the target is not ⁇ Beam except for the center direction ⁇ Beam in the beam forming direction v d . Therefore, the target angle measuring unit 10 provides a latitude in the angle measurement value in the v-axis direction of the target, depending on the angle width A ⁇ in the v-axis direction of the reflected wave. That is, the target angle measuring unit 10 calculates ⁇ Beam ⁇ A ⁇ /2 as the angle measurement value in the v-axis direction of the target.
- the angle width estimation device 4 is configured in such a way as to include: the beam forming unit 5 to acquire one reception array signal of the reflected wave from the target to be observed, and form, from the one reception array signal, a plurality of null beams having nulls in the arrival direction of the reflected wave, and having null widths which are widths of the nulls and different from each other; and the angle width estimating unit 9 to compare powers of the plurality of null beams formed by the beam forming unit 5 with each other, and estimate the null width indicating the angle width of the reflected wave on the basis of the comparison result of the powers. Therefore, the angle width estimation device 4 can estimate the angle width of the reflected wave from the target from the one reception array signal.
- the correlation matrix R xx expressed by Expression (10) is calculated, by acquiring an Hadamard product of a correlation matrix R xx ′ of the reflected wave having no angle width expressed by the following Expression (11) and a covariance matrix taper (CMT) matrix T.
- the correlation matrix R xx expressed by Expression (10) is the correlation matrix R xx of the reflected wave having the angle width, and is acquired without numerical integration.
- R xx T ⁇
- R xx ′ T ⁇ ( a ⁇ ( u c , v c ) ⁇ a ⁇ ( u c , v c ) H + ⁇ ⁇ I ) ( 10 )
- R xx ′ ( a ⁇ ( u c , v c ) ⁇ a ⁇ ( u c , v c ) H + ⁇ ⁇ I ( 11 )
- ⁇ represents an Hadamard product.
- FIG. 9 is a configuration diagram illustrating a target angle measurement device including the angle width estimation device 4 according to the second embodiment.
- the same reference sign as that in FIG. 1 represents the same or corresponding portion, so that the description thereof is omitted.
- the beam forming unit 11 is provided with coherent integrating units 6 - 1 to 6 -M, pulse compressing units 7 - 1 to 7 -M, and the DBF unit 12 .
- the DBF unit 12 is implemented by a DBF circuit 26 illustrated in FIG. 10 , for example.
- the table 12 a stores a plurality of weights regarding the null beam formed by the DBF unit 12 as is the case of the table 8 a illustrated in FIG. 1 .
- each of the coherent integrating units 6 - 1 to 6 -M, the pulse compressing units 7 - 1 to 7 -M, the DBF unit 12 , and the angle width estimating unit 9 which are components of the angle width estimation device 4 , is implemented by dedicated hardware as illustrated in FIG. 10 . That is, it is supposed that the angle width estimation device 4 is implemented by a coherent integrating circuit 21 , a pulse compressing circuit 22 , the DBF circuit 26 , and an angle width estimating circuit 24 .
- Each of the coherent integrating circuit 21 , the pulse compressing circuit 22 , the DBF circuit 26 , and the angle width estimating circuit 24 corresponds to, for example, a single circuit, a composite circuit, a programmed processor, a parallel-programmed processor, an ASIC, a FPGA, or a combination thereof.
- the components of the angle width estimation device 4 are not limited to those implemented by the dedicated hardware, and the angle width estimation device 4 may also be implemented by software, firmware, or a combination of the software and firmware.
- angle width estimation device 4 When the angle width estimation device 4 is implemented by the software, firmware or the like, a program for causing a computer to execute each processing procedure in the coherent integrating units 6 - 1 to 6 -M, the pulse compressing units 7 - 1 to 7 -M, the DBF unit 12 , and the angle width estimating unit 9 is stored in a memory 31 illustrated in FIG. 3 . Then, a processor 32 illustrated in FIG. 3 executes the program stored in the memory 31 .
- the null width of the null beam B NULL,j formed by the DBF unit 8 coincides with the angle width of the reflected wave. Therefore, there is a possibility that the null beam B NULL,j having the null with the null width wider than the angle width of the reflected wave is formed by first processing by the DBF unit 8 .
- null widths B ⁇ and B ⁇ cannot be searched for in some cases.
- the DBF unit 12 forms the beam B NULL,j having the beam at the end.
- C represents a constraint matrix
- H represents a constraint response value for the constraint matrix C.
- FIG. 11 is an explanatory diagram illustrating a null beam formed by the DBF unit 12 .
- the null beam B NULL,j is formed in null forming directions u c and v c .
- a hatched region is a null region, the null width in the u-axis direction is (u c + ⁇ u/2) ⁇ (u c ⁇ u/2), and the null width in the v-axis direction is (v c + ⁇ v/2) ⁇ (v c ⁇ v/2).
- the beams a longitudinal direction each of which is a direction parallel to the v-axis are formed at u c ⁇ u/2, and the beams a longitudinal direction each of which are a direction parallel to the u-axis is formed at v c ⁇ v/2.
- the null width of the null beam B NULL,j formed by the DBF unit 12 is likely to be narrower than the null width of the null beam B NULL,j formed by the DBF unit 8 illustrated in FIG. 1 .
- FIG. 12 is a configuration diagram illustrating a target angle measurement device including the angle width estimation device 4 according to the third embodiment.
- the same reference sign as that in FIGS. 1 and 9 represents the same or corresponding part, so that the description thereof is omitted.
- FIG. 13 is a hardware configuration diagram illustrating hardware of a digital signal processing unit in the target angle measurement device illustrated in FIG. 12 .
- the same reference sign as that in FIGS. 2 and 10 represents the same or corresponding part, so that the description thereof is omitted.
- the beam forming unit 13 is provided with coherent integrating units 6 - 1 to 6 -M, pulse compressing units 7 - 1 to 7 -M, and the DBF unit 14 .
- the beam forming unit 13 forms, from one reception array signal r ARE output from a reflected wave receiving unit 1 , a plurality of null beams having nulls in an arrival direction of a reflected wave, and having null widths different from each other, as is the case of the beam forming unit 5 illustrated in FIG. 1 .
- the beam forming unit 13 rotates the forming direction of each null beam.
- the DBF unit 14 is implemented by a DBF circuit 27 illustrated in FIG. 13 , for example.
- the DBF unit 14 is provided with a table 12 a as is the case of the DBF unit 12 illustrated in FIG. 9 .
- the DBF unit 14 is provided with the table 12 a .
- the DBF unit 14 rotates the forming direction of the null beam B NULL,j .
- each of the coherent integrating units 6 - 1 to 6 -M, the pulse compressing units 7 - 1 to 7 -M, the DBF unit 14 , and the angle width estimating unit 9 which are components of the angle width estimation device 4 , is implemented by dedicated hardware as illustrated in FIG. 13 . That is, it is supposed that the angle width estimation device 4 is implemented by a coherent integrating circuit 21 , a pulse compressing circuit 22 , the DBF circuit 27 , and an angle width estimating circuit 24 .
- Each of the coherent integrating circuit 21 , the pulse compressing circuit 22 , the DBF circuit 27 , and the angle width estimating circuit 24 corresponds to, for example, a single circuit, a composite circuit, a programmed processor, a parallel-programmed processor, an ASIC, a FPGA, or a combination thereof.
- the components of the angle width estimation device 4 are not limited to those implemented by the dedicated hardware, and the angle width estimation device 4 may also be implemented by software, firmware, or a combination of the software and firmware.
- angle width estimation device 4 When the angle width estimation device 4 is implemented by the software, firmware or the like, a program for causing a computer to execute each processing procedure in the coherent integrating units 6 - 1 to 6 -M, the pulse compressing units 7 - 1 to 7 -M, the DBF unit 14 , and the angle width estimating unit 9 is stored in a memory 31 illustrated in FIG. 3 . Then, a processor 32 illustrated in FIG. 3 executes the program stored in the memory 31 .
- the DBF unit 14 forms the null beam B NULL,j by a method similar to that of the DBF unit 8 illustrated in FIG. 1 or the DBF unit 12 illustrated in FIG. 9 .
- R rot ( ⁇ r ⁇ o ⁇ t ) [ cos ⁇ ⁇ r ⁇ o ⁇ t - sin ⁇ ⁇ rot sin ⁇ ⁇ rot cos ⁇ ⁇ rot ] ( 16 )
- R r represents the element coordinate matrix after coordinate rotation.
- FIG. 14 is an explanatory diagram illustrating the rotation of the forming direction of the null beam.
- a hatched region is a null region.
- beams are formed at ends in a u-axis direction of the null beam, and beams are formed at ends in a v-axis direction of the null beam.
- the null beam may also be that without the beams formed at the ends in the u-axis direction of the null beam and without the beams formed at the ends in the v-axis direction of the null beam.
- the DBF unit 14 can control null forming directions u c and v c , by calculating a CMT matrix T expressed by the following Expression (18) using the element coordinate matrix R r after coordinate rotation.
- the DBF unit 14 can also control a null expanding direction, by calculating the CMT matrix T expressed by the following Expression (18) using the element coordinate matrix R r after coordinate rotation.
- the angle width estimation device 4 illustrated in FIG. 12 is configured in such a manner that the beam forming unit 13 rotates the forming direction of the null beam. Therefore, as is the case of the angle width estimation device 4 illustrated in FIG. 1 , the angle width estimation device 4 illustrated in FIG. 12 can estimate the angle width of the reflected wave from the target from one reception array signal, and can form the null beam suitable for a shape of the target, a posture of the target or the like.
- the present disclosure is suitable for an angle width estimation device, an angle width estimation method, and a target angle measurement device.
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Electromagnetism (AREA)
- Quality & Reliability (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
Description
-
- Non-Patent Literature 1: Ogawa, Kikuma et al., “Angular Spread Estimation by Derivative Constrained Capon Estimator Using Integrated Mode Vector”, Journal of The Institute of Electronics, Information and Communication Engineers, Vol. J92-B, No. 6, June 2009.
r ARE(t)=a(θ,ϕ)s(t)+n(t) (1)
w=αR xx −1 a(u d ,v d) (5)
Th=Gr−SNR (6)
P 1(θBeam,ϕBeam,θNull,ϕnull)=|w 1(θBeam,ϕBeam,θNull,ϕnull)H r APC(t)|2 (8)
-
- [j-th Processing] (j=2, . . . , J)
P j(θBeam,ϕBeam,θNull,ϕnull)=|w j(θBeam,ϕBeam,θNull,ϕnull)H r APC(t)|2 (9)
In Expression (10), ⊙ represents an Hadamard product.
-
- 1: reflected wave receiving unit, 2-1 to 2-M: reception antenna element, 3-1 to 3-M: A/D converter, 4: angle width estimation device, 5, 11, 13: beam forming unit, 6-1 to 6-M: coherent integrating unit, 7-1 to 7-M: pulse compressing unit, 8, 12, 14: DBF unit, 8 a, 12 a: table, 9: angle width estimating unit, 10: target angle measuring unit, 21: coherent integrating circuit, 22: pulse compressing circuit, 23, 26, 27: DBF circuit, 24: angle width estimating circuit, 25: target angle measuring circuit, 31: memory, 32: processor
Claims (5)
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| PCT/JP2021/007704 WO2022185379A1 (en) | 2021-03-01 | 2021-03-01 | Angle width estimation device, angle width estimation method, and target angle measurement device |
Related Parent Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| PCT/JP2021/007704 Continuation WO2022185379A1 (en) | 2021-03-01 | 2021-03-01 | Angle width estimation device, angle width estimation method, and target angle measurement device |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| US20230314127A1 US20230314127A1 (en) | 2023-10-05 |
| US12480760B2 true US12480760B2 (en) | 2025-11-25 |
Family
ID=83153961
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| US18/207,235 Active 2041-07-23 US12480760B2 (en) | 2021-03-01 | 2023-06-08 | Angle width estimation device, angle width estimation method, and target angle measurement device |
Country Status (4)
| Country | Link |
|---|---|
| US (1) | US12480760B2 (en) |
| EP (1) | EP4283325B1 (en) |
| JP (1) | JP7258264B2 (en) |
| WO (1) | WO2022185379A1 (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US12096242B2 (en) * | 2022-01-19 | 2024-09-17 | Lg Electronics Inc. | Method and apparatus for reducing interference effects in wireless communication systems |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP3738705B2 (en) | 2001-06-05 | 2006-01-25 | 松下電器産業株式会社 | Adaptive antenna device |
| US20070047742A1 (en) * | 2005-08-26 | 2007-03-01 | Step Communications Corporation, A Nevada Corporation | Method and system for enhancing regional sensitivity noise discrimination |
| JP4784976B2 (en) | 2005-09-02 | 2011-10-05 | Kddi株式会社 | Radio wave arrival direction estimation device, radio wave arrival direction estimation program, and recording medium |
| US20170338900A1 (en) * | 2015-02-17 | 2017-11-23 | Mitsubishi Electric Corporation | Receiver apparatus and reception method |
| US20220200665A1 (en) * | 2020-12-18 | 2022-06-23 | Nokia Technologies Oy | Beam steering apparatus and method |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS6122274A (en) * | 1984-07-11 | 1986-01-30 | Natl Space Dev Agency Japan<Nasda> | Method for data processing of synthetic aperture radar |
| US6002983A (en) | 1997-08-27 | 1999-12-14 | Delphi Technologies, Inc. | Angle extent estimation method for a motor vehicle object detection system |
| US5982319A (en) | 1998-03-12 | 1999-11-09 | Northrop Grumman Corporation | UHF synthetic aperture radar |
-
2021
- 2021-03-01 JP JP2023503545A patent/JP7258264B2/en active Active
- 2021-03-01 WO PCT/JP2021/007704 patent/WO2022185379A1/en not_active Ceased
- 2021-03-01 EP EP21928957.6A patent/EP4283325B1/en active Active
-
2023
- 2023-06-08 US US18/207,235 patent/US12480760B2/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP3738705B2 (en) | 2001-06-05 | 2006-01-25 | 松下電器産業株式会社 | Adaptive antenna device |
| US20070047742A1 (en) * | 2005-08-26 | 2007-03-01 | Step Communications Corporation, A Nevada Corporation | Method and system for enhancing regional sensitivity noise discrimination |
| JP4784976B2 (en) | 2005-09-02 | 2011-10-05 | Kddi株式会社 | Radio wave arrival direction estimation device, radio wave arrival direction estimation program, and recording medium |
| US20170338900A1 (en) * | 2015-02-17 | 2017-11-23 | Mitsubishi Electric Corporation | Receiver apparatus and reception method |
| US20220200665A1 (en) * | 2020-12-18 | 2022-06-23 | Nokia Technologies Oy | Beam steering apparatus and method |
Non-Patent Citations (6)
| Title |
|---|
| International Search Report (PCT/ISA/210), issued in PCT/JP2021/007704, dated Apr. 20, 2021. |
| Ogawa et al., "Angular Spread Estimation by Derivative Constrained Capon Estimator Using Integrated Mode Vector", IEICE, B—Abstracts of IEICE Transactions on Communications (Japanese Edition), 2009, vol. J92-B, No. 6, pp. 921-929. |
| Tsukagoshi et al., "A Sparse Signal Processing for Angular Estimation of Spread Wave Sources", IEICE, IEICE Technical Report, 2013, vol. 112, No. 391, pp. 67-71. |
| International Search Report (PCT/ISA/210), issued in PCT/JP2021/007704, dated Apr. 20, 2021. |
| Ogawa et al., "Angular Spread Estimation by Derivative Constrained Capon Estimator Using Integrated Mode Vector", IEICE, B—Abstracts of IEICE Transactions on Communications (Japanese Edition), 2009, vol. J92-B, No. 6, pp. 921-929. |
| Tsukagoshi et al., "A Sparse Signal Processing for Angular Estimation of Spread Wave Sources", IEICE, IEICE Technical Report, 2013, vol. 112, No. 391, pp. 67-71. |
Also Published As
| Publication number | Publication date |
|---|---|
| EP4283325A4 (en) | 2024-03-27 |
| JPWO2022185379A1 (en) | 2022-09-09 |
| EP4283325A1 (en) | 2023-11-29 |
| WO2022185379A1 (en) | 2022-09-09 |
| EP4283325B1 (en) | 2025-02-26 |
| US20230314127A1 (en) | 2023-10-05 |
| JP7258264B2 (en) | 2023-04-14 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US11933909B2 (en) | Method and apparatus for estimating direction of arrival using generation of virtual received signals | |
| US6184830B1 (en) | Compensation of direction finding estimates for polarimetric errors | |
| JP5659472B2 (en) | Direction of arrival estimation apparatus and method | |
| US10523288B2 (en) | Wireless receiver, wireless receiving method, and wireless system | |
| US7961147B1 (en) | Long baseline phase interferometer ambiguity resolution using frequency differences | |
| EP2662699A1 (en) | Detection and ranging apparatus and ranging method | |
| US8130147B2 (en) | Signal wave arrival angle measuring device | |
| KR20160134436A (en) | Method for estimating direction of arrival and apparatus for estimating direction of arrival using the same | |
| US20060208947A1 (en) | Apparatus and method for estimating direction of arrival of radio wave | |
| US10852410B2 (en) | Angle estimating apparatus and angle estimating method | |
| KR101796472B1 (en) | Radar apparatus and DOA estimation method using the same | |
| US10732273B2 (en) | Radar device for vehicle and method for estimating angle of target using same | |
| US20190369223A1 (en) | Radar apparatus | |
| KR102422396B1 (en) | Method of spatial interpolation for linear phased array antenna and appratus thereof | |
| EP3825716B1 (en) | Radar device and target angle measurement method | |
| US12480760B2 (en) | Angle width estimation device, angle width estimation method, and target angle measurement device | |
| US12032053B2 (en) | Angle measuring device, angle measuring method, and in-vehicle device | |
| JP2002243824A (en) | Angle measuring device and angle measuring method | |
| JP7413850B2 (en) | Angle estimation device and method for object position, and radar device | |
| JP7031086B2 (en) | Signal processing equipment, radar equipment and signal processing methods | |
| JPWO2018167952A1 (en) | Adaptive array antenna device | |
| JP6980570B2 (en) | Target detector and signal processing method | |
| JP4977849B2 (en) | Radio wave arrival direction detector | |
| KR101240415B1 (en) | Method for processing sidelobe blanking in adaptive array radar | |
| US12372641B2 (en) | Ego velocity estimator for radar systems |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| AS | Assignment |
Owner name: MITSUBISHI ELECTRIC CORPORATION, JAPAN Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNORS:TAKAHASHI, YOSHIKI;TAKAHASHI, RYUHEI;REEL/FRAME:063895/0649 Effective date: 20230510 |
|
| FEPP | Fee payment procedure |
Free format text: ENTITY STATUS SET TO UNDISCOUNTED (ORIGINAL EVENT CODE: BIG.); ENTITY STATUS OF PATENT OWNER: LARGE ENTITY |
|
| STPP | Information on status: patent application and granting procedure in general |
Free format text: DOCKETED NEW CASE - READY FOR EXAMINATION |
|
| STPP | Information on status: patent application and granting procedure in general |
Free format text: NON FINAL ACTION MAILED |
|
| STPP | Information on status: patent application and granting procedure in general |
Free format text: RESPONSE TO NON-FINAL OFFICE ACTION ENTERED AND FORWARDED TO EXAMINER |
|
| STPP | Information on status: patent application and granting procedure in general |
Free format text: NON FINAL ACTION MAILED |
|
| STPP | Information on status: patent application and granting procedure in general |
Free format text: RESPONSE TO NON-FINAL OFFICE ACTION ENTERED AND FORWARDED TO EXAMINER |
|
| STPP | Information on status: patent application and granting procedure in general |
Free format text: ALLOWED -- NOTICE OF ALLOWANCE NOT YET MAILED Free format text: NOTICE OF ALLOWANCE MAILED -- APPLICATION RECEIVED IN OFFICE OF PUBLICATIONS |
|
| STPP | Information on status: patent application and granting procedure in general |
Free format text: PUBLICATIONS -- ISSUE FEE PAYMENT RECEIVED |
|
| STPP | Information on status: patent application and granting procedure in general |
Free format text: PUBLICATIONS -- ISSUE FEE PAYMENT VERIFIED |
|
| STCF | Information on status: patent grant |
Free format text: PATENTED CASE |