CN109613475A - A kind of arrangement method and target DOA evaluation method of the relatively prime battle array of adaptive optimization - Google Patents
A kind of arrangement method and target DOA evaluation method of the relatively prime battle array of adaptive optimization Download PDFInfo
- Publication number
- CN109613475A CN109613475A CN201910122620.5A CN201910122620A CN109613475A CN 109613475 A CN109613475 A CN 109613475A CN 201910122620 A CN201910122620 A CN 201910122620A CN 109613475 A CN109613475 A CN 109613475A
- Authority
- CN
- China
- Prior art keywords
- battle array
- relatively prime
- prime battle
- difference
- optimization
- 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.)
- Granted
Links
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
Abstract
The invention discloses a kind of arrangement method of the relatively prime battle array of adaptive optimization and target DOA evaluation methods, include the following steps, step 1: any one submatrix position in relatively prime battle array being inverted, i.e., coordinate position is inverted to negative axis from positive axis；Step 2: the submatrix after inverting in step 1 being translated, shift factor is obtained by optimization method as FA；The arrangement algorithm of the relatively prime battle array of adaptive optimization of the present invention is easy, easy, has stronger universality；It does not need to recalculate optimal arrangement mode for the relatively prime battle array of different units number；There are more continuous unit numbers with the relatively prime battle array of conventional differential compared with using the obtained relatively prime battle array of difference of the present invention, aperture number is bigger, can be under the premise of not adding unit number for the fields such as communicating, being imaged, position, acquisition better information treatment effect.
Description
Technical field
The present invention relates to field of antenna, and in particular to the arrangement method of a kind of relatively prime battle array of adaptive optimization and relatively prime based on this
The target DOA evaluation method of the battle array relatively prime battle array of difference generated.
Background technique
Antenna array as communication, imaging, positioning etc. in numerous areas the transceiver of signal and be widely studied.How
Realize that aperture as big as possible is the technical problem in field of antenna in the case where Finite Array unit.Traditional N member antenna array
N-1 target of detection is at best able to since aperture is limited；In order to solve this problem, relatively prime battle array is suggested.Although by relatively prime
The relatively prime battle array of difference that battle array generates has the aperture bigger than traditional relatively prime battle array or linear battle array, i.e. freedom degree wants high compared with traditional array,
To further significantly improve the reception ability and detectability of signal.But sequential cells number is again without sufficiently being dug
The problem of pick, i.e. freedom degree do not reach achievable maximum value, hole is not reduced fundamentally.In addition, existing reduction hole is
Improve freedom degree method it is complicated for operation, it is computationally intensive, for different mutual prime rwmbers composition relatively prime battle array needs recalculate optimization
Arrangement mode, operability wretched insufficiency.
Summary of the invention
The present invention provide it is a kind of it is easy to operate, freedom degree is high, do not need to recalculate for the relatively prime battle array of different units number
The arrangement method of the relatively prime battle array of the adaptive optimization of optimal arrangement, and the relatively prime battle array of the relatively prime obtained difference of battle array is used for target
The method of DOA estimation.
The technical solution adopted by the present invention is that: a kind of arrangement method of the relatively prime battle array of adaptive optimization, comprising the following steps:
Step 1: any one submatrix position in relatively prime battle array being inverted, i.e., coordinate position is inverted to negative axis from positive axis；
Step 2: the submatrix after inverting in step 1 is translated into the relatively prime battle array after being optimized, shift factor FA.
Further, the shift factor FA is as follows:
In formula: AdisFor the location matrix of subarray A, BdisFor the location matrix of subarray B, function g (a, b) is matrix a
With maximum value continuous in the difference of b.
One kind being based on the relatively prime battle array target DOA evaluation method of difference, comprising the following steps:
S1: the relatively prime battle array of difference that the relatively prime battle array after the optimization obtained according to step 2 is optimized makees the relatively prime battle array of the difference
For dual-mode antenna battle array；
S2: penetrating direct impulse using the dual-mode antenna paroxysm in step S1, and returning for target is received after extraterrestrial target reflects
Wave signal；
S3: the echo-signal that step S2 is obtained carries out matrixing processing；
S4: the DOA from different path targets is calculated using Capon method, completes the estimation of target DOA.
The method of the present invention the utility model has the advantages that
(1) the arrangement algorithm of the relatively prime battle array of adaptive optimization of the present invention is easy, easy, has stronger universality；It does not need
Optimal arrangement mode is recalculated for the relatively prime battle array of different units number；
(2) the present invention relatively prime battle array of difference generated has compared with the relatively prime battle array of difference that the relatively prime battle array of tradition obtains and more connects
Continuous unit number, aperture number is bigger, for the fields such as communicating, being imaged, position, can obtain under the premise of not adding unit number
Obtain better information treatment effect；
(3) the target DOA evaluation method in the present invention can significantly improve the positioning accuracy and resolution ratio and reduction of DOA
Secondary lobe.
Detailed description of the invention
Fig. 1 is the schematic diagram of the relatively prime battle array of tradition and submatrix.
Fig. 2 is to work as Mc=3, NcWhen=8, by the arrangement schematic diagram for the relatively prime battle array of conventional differential that the relatively prime battle array of tradition generates.
Fig. 3 is the flow chart that the present invention optimizes relatively prime battle array arrangement method.
Fig. 4 is the schematic diagram of the relatively prime battle array and its submatrix after optimizing in the present invention.
Fig. 5 is to work as Mc=3, NcWhen=8, the arrangement of the relatively prime battle array of optimization difference generated by the relatively prime battle array of optimization of the invention is shown
It is intended to.
Fig. 6 is the power spectrum comprising DOA information based on the relatively prime battle array of conventional differential and the optimization relatively prime battle array of difference.
Fig. 7 is the enlarged diagram of Fig. 6.
Specific embodiment
The present invention will be further described in the following with reference to the drawings and specific embodiments.
A kind of arrangement method of the relatively prime battle array of adaptive optimization, comprising the following steps:
Step 1: any one submatrix position in relatively prime battle array being inverted, i.e., coordinate position is inverted to negative axis from positive axis；
Step 2: the submatrix after inverting in step 1 is translated into the relatively prime battle array after being optimized, shift factor FA.
Assuming that any one submatrix location A in relatively prime battle array inverts, the battle array A of relatively prime a period of time after being optimized；After optimizing
Relatively prime submatrix A and former relatively prime battle array in the new relatively prime battle array of submatrix B composition, and it is mutual using the difference that new relatively prime battle array is optimized
Matter battle array.
Shift factor FA is as follows:
In formula: AdisFor the location matrix of subarray A, BdisFor the location matrix of subarray B, function g (a, b) is matrix a
With maximum value continuous in the difference of b.
Choose any two relatively prime natural number McAnd Nc, it is N=according to traditional relatively prime battle array arrangement array element number
2Mc+Nc-1；If not carrying out any its effective aperture number of operation to the battle array is N, the number of targets that can be detected is up to N-1.
Bigger aperture in order to obtain then needs to deform the relatively prime battle array of tradition to realize virtual uniform linear array, which has more
More sequential cells.Therefore, available bigger aperture.The deformation operation is formed by array and is known as the relatively prime battle array of difference.Phase
The array element position answered is SDCA={ ± (ncMcd-mcNcD) }, wherein (1≤nc≤Nc- 1), (1≤mc≤2Mc-1).This kind of side
The sequential cells number that the relatively prime battle array of the difference that formula generates generates is limited, i.e., freedom degree is limited or more hole occurs.
Precisely in order to the above-mentioned problems in the prior art is solved, any one submatrix position of the relatively prime battle array of tradition is anti-
Turn, i.e., coordinate position is inverted to negative axis from positive axis；Submatrix after reversion is translated, shift factor FA.
FA is obtained by solving the maximum successive value of the range difference matrix of two submatrixs, the most Dalian of the range difference matrix
The physical significance of continuous value is the maximum sequential cells number of the relatively prime battle array of difference.
Detailed process is as follows:
Two sub- array element numbers for taking relatively prime battle array are respectively McAnd Nc, it is assumed that Mc< Nc.At this point, obtained by conventional method
Relatively prime array element number is N=2Mc+Nc-1.If not carrying out any its effective aperture number of operation to the battle array is N, can be detected
Number of targets is up to N-1.Corresponding relatively prime battle array arrangement is as shown in Figure 1, the position of submatrix A is [0, Mc,2Mc,…,(Nc-1)Mc]
The position of d, submatrix B are [0, Nc,2Nc,…,(2Mc-1)Nc] d, d is the half-wavelength of Antenna Operation centre frequency.Utilize this side
The minimum 2M of freedom degree for the relatively prime battle array of difference that formula generatescNc+ 1, such as work as Mc=3, NcWhen=8, the corresponding relatively prime battle array of difference
The position arrangement of unit is as shown in Figure 2.It can be seen that the number of sequential cells is 26 × 2+1=53, i.e. freedom degree is greater than for 53
2McNc+ 1=49.The relatively prime battle array of the difference shares 4McNc-2NcIt is i.e. dark-grey as shown in Figure 2 that+1=81 apertures contain 14 holes
Shown in the circle of color.Therefore, effective aperture number achieved by the relatively prime battle array of difference generated using traditional approach is 81-14=67
It is a.
In order under the premise of not increasing array element number, obtain sequential cells number as much as possible i.e. freedom degree and
Effective aperture number proposes the method for the present invention.Obtaining its new position after a submatrix such as A in relatively prime battle array is inverted is
[-(Nc-1)Mc,…,-2Mc,-Mc,0]d；By unit maximum position after transformation from traditional max ((Nc-1)Mcd,(2Mc-1)
NcD) i.e. from (Nc-1)McD, (2Mc-1)NcBeing maximized in d becomes [(Nc-1)Mc+(2Mc-1)Nc]d.It can be seen that [(Nc-1)Mc
+(2Mc-1)Nc] d is consistently greater than max ((Nc-1)Mcd,(2Mc-1)Ncd).Therefore, difference caused by the relatively prime battle array after reversion is mutual
Matter array element maximum position perseverance is greater than the unit maximum position of the relatively prime battle array of difference caused by conventional method.
But simple reversion not can solve the optimal problem of sequential cells number, therefore propose optimization shift factor FA
Universal method.The specific method is as follows:
In formula: AdisFor the location matrix of subarray A, BdisFor the location matrix of subarray B, function g (a, b) is matrix a
With maximum value continuous in the difference of b；A=-A in the present inventiondis+ FAd replaces original mutually after subarray A is carried out reversion translation transformation
Proton battle array A, and submatrix B does not do any operation.Therefore, the method only needs to convert a submatrix, relatively simple.Separately
Outside, the result that optimization method obtains is to obtain the optimal translation of relatively prime battle array corresponding to the relatively prime battle array of the most difference of sequential cells number
Factor FA.
As shown in figure 3, for the relatively prime battle array for having set number of unit, setting shift factor FA's takes the optimization process of FA
It is worth range, which can be very big, and step-length can be with very little；FA value range and step-length are determined according to practical application scene；?
The initial value of FA is used to acquire the relatively prime battle array of difference, and obtained sequential cells number is Q1, second value of FA is with similarly
Method obtains the relatively prime battle array of difference, and obtained sequential cells number is Q2；Compare Q1And Q2Size, if Q1< Q2It sets up, explanation
Initial value is not that optimal solution is not optimal shift factor number FA, then excludes the initial value of FA value range, and second is worth
It is set as initial value, and obtains the relatively prime battle array of corresponding difference, calculates and obtains sequential cells number, be defined as Q1, former value model
The third value enclosed becomes second value of new value range, and the relatively prime battle array of corresponding difference and sequential cells number is calculated,
It is defined as Q2.Later with it is discussed above identical, until finding maximum Q2The value of corresponding shift factor FA, then this value is FA
Optimal value.If Q1< Q2It is invalid, illustrate using the sequential cells number of the relatively prime battle array of the obtained difference of initial value than second
The sequential cells number of a obtained relatively prime battle array of difference of value is more, then third is worth to second substituted in former value range
Value calculates and obtains the relatively prime battle array of corresponding difference and sequential cells number is defined as Q2, then compare Q1And Q2Size, in repetition
Step is stated, if still obtaining Q after having traversed all values1< Q2It is invalid.Illustrate that initial value is the optimal solution of shift factor FA.
Q under normal circumstances1< Q2It permanent can not set up or permanent invalid, need to select above-mentioned corresponding processing side after determining
Formula, to obtain the optimal solution of shift factor FA.Battle array A of relatively prime a period of time after reversion translate as shown in Figure 4.By new son
Battle array A and atom battle array B carries out alternate position spike calculating, obtains the relatively prime battle array of new difference, this gust is the relatively prime battle array of difference of optimal arrangement, connects
Continuous unit number is most.Equally with Mc=3, NcFor=8, by optimal shift factor FA=8 obtained by the above method, accordingly
Difference relatively prime battle array distribution it is as shown in Figure 5.As it can be seen that sequential cells number is 39 × 2+1=79, i.e., freedom degree is 79, than passing
System method obtains 53 big.In addition, the relatively prime battle array of the difference, which shares 107 apertures, contains 14 holes, depth as illustrated in fig. 5
Gray circles.Effective aperture number is 93.In the case where no change relatively prime array element number, the optimization that is generated by it
The relatively prime battle array of difference increases the number, freedom degree and effective aperture number of sequential cells compared with the relatively prime battle array of conventional differential.
In order to further be compared, it regard the relatively prime battle array of difference after the relatively prime battle array of conventional differential and optimization as dual-mode antenna
Battle array, launching centre frequency is in the Gaussian modulation pulse that the 3.5GHz duration is 0.5ns as direct impulse.Assuming that environment noise
Than being respectively 37 °, 39 ° and 41 ° for -10dB, target DOA, first antenna horizontal distance of extraterrestrial target distance is 1km, is utilized
Logarithm normal distribution shadow model simulates communication environments, which is according in environment between signal sending end and receiving end
The difference of obstacle density realizes different channels.The model is widely used in the signal processing, utilizes traditional Capon
Algorithm carries out DOA estimation, as a result as shown in Figure 6 and Figure 7.Wherein Capon is common beam forming algorithm.
Since target DOA is at a distance of relatively closely, the relatively prime battle array of traditional difference can only tell the position of two DOA, respectively at 38 °
With 41.5 °, wherein 38 ° are center that target DOA is 37 ° and 39 °, error is 1 °, and the error of 41.5 ° of DOA of estimation is
0.5°.But using the relatively prime battle array of optimization difference caused by method proposed by the present invention, it can clearly obtain the position of 3 DOA
It sets respectively at 36.9 °, 39.2 ° and 41 °, error is respectively 0.1 °, 0.2 ° and 0 °.It significantly improves the positioning accuracy of DOA and divides
Resolution.Herein in the power spectrum comprising DOA information, using the optimization obtained waveform of the relatively prime battle array of difference than utilizing conventional differential
The waveform secondary lobe that relatively prime battle array obtains is lower, convenient for the further identification of DOA, therefore, relatively prime battle array arrangement side proposed by the invention
Method under the premise of not changing DOA evaluation method, can realize the DOA positioning of high-precision low sidelobe.
Arrangement method of the present invention can effectively increase number, freedom degree and the effective aperture of the relatively prime battle array sequential cells of difference
Number.Two submatrixs for any relatively prime battle array, that is, relatively prime battle array are that any relatively prime unit number is applicable in, and do not need to calculate repeatedly, only
It needs the unit number input of two submatrixs optimal arrangement can be obtained.From conventional method for different relatively prime battle arrays need it is different
Optimization method is different, and invention increases the universalities of method.
Claims (3)
1. a kind of arrangement method of the relatively prime battle array of adaptive optimization, which comprises the following steps:
Step 1: any one submatrix position in relatively prime battle array being inverted, i.e., coordinate position is inverted to negative axis from positive axis；
Step 2: the submatrix after inverting in step 1 is translated into the relatively prime battle array after being optimized, shift factor FA.
2. a kind of arrangement method of the relatively prime battle array of adaptive optimization according to claim 1, which is characterized in that the translation because
The optimization method of sub- FA is as follows:
In formula: AdisFor the location matrix of subarray A, BdisFor the location matrix of subarray B, function g (a, b) is matrix a and b
Continuous maximum value in difference.
3. a kind of obtain the target DOA evaluation method of the relatively prime battle array of difference using such as claim 1 method, which is characterized in that including
Following steps:
S1: the relatively prime battle array of difference that the relatively prime battle array after the optimization obtained according to step 2 is optimized, using the relatively prime battle array of the difference as receipts
Send out antenna array；
S2: penetrating direct impulse using the dual-mode antenna paroxysm in step S1, and the echo letter of target is received after extraterrestrial target reflects
Number；
S3: the echo-signal that step S2 is obtained carries out matrixing processing；
S4: the DOA from different path targets is calculated using Capon method, completes the estimation of target DOA.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910122620.5A CN109613475B (en) | 2019-02-18 | 2019-02-18 | Self-adaptive optimization co-prime matrix arrangement method and target DOA estimation method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910122620.5A CN109613475B (en) | 2019-02-18 | 2019-02-18 | Self-adaptive optimization co-prime matrix arrangement method and target DOA estimation method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109613475A true CN109613475A (en) | 2019-04-12 |
| CN109613475B CN109613475B (en) | 2022-11-15 |
Family
ID=66022121
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910122620.5A Active CN109613475B (en) | 2019-02-18 | 2019-02-18 | Self-adaptive optimization co-prime matrix arrangement method and target DOA estimation method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109613475B (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110133576A (en) * | 2019-05-23 | 2019-08-16 | 成都理工大学 | Biradical relatively prime MIMO array orientation algorithm for estimating based on cascade residual error network |
| CN112014792A (en) * | 2020-08-28 | 2020-12-01 | 陕西理工大学 | Improved two-dimensional co-prime array DOA estimation method |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105182293A (en) * | 2015-08-25 | 2015-12-23 | 西安电子科技大学 | Method for estimating DOA and DOD of MIMO radar based on co-prime array |
| CN106324558A (en) * | 2016-08-30 | 2017-01-11 | 东北大学秦皇岛分校 | Broadband signal DOA estimation method based on co-prime array |
-
2019
- 2019-02-18 CN CN201910122620.5A patent/CN109613475B/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105182293A (en) * | 2015-08-25 | 2015-12-23 | 西安电子科技大学 | Method for estimating DOA and DOD of MIMO radar based on co-prime array |
| CN106324558A (en) * | 2016-08-30 | 2017-01-11 | 东北大学秦皇岛分校 | Broadband signal DOA estimation method based on co-prime array |
Non-Patent Citations (2)
| Title |
|---|
| SI QIN ET AL.: ""Generalized Coprime Array Configurations for Direction-of-Arrival Estimation"", 《IEEE TRANSACTIONS ON SIGNAL PROCESSING》 * |
| TAREK HASAN AL MAHMUD ET AL.: ""Interpolating Coprime Arrays With Translocated and Axis Rotated Compressed Subarrays by Iterative Power Factorization for DOA Estimation"", 《IEEE ACCESS》 * |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110133576A (en) * | 2019-05-23 | 2019-08-16 | 成都理工大学 | Biradical relatively prime MIMO array orientation algorithm for estimating based on cascade residual error network |
| CN110133576B (en) * | 2019-05-23 | 2020-11-27 | 成都理工大学 | Double-base co-prime MIMO array azimuth estimation algorithm based on cascade residual error network |
| CN112014792A (en) * | 2020-08-28 | 2020-12-01 | 陕西理工大学 | Improved two-dimensional co-prime array DOA estimation method |
| CN112014792B (en) * | 2020-08-28 | 2023-03-31 | 陕西理工大学 | Improved two-dimensional co-prime array DOA estimation method |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109613475B (en) | 2022-11-15 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105259552B (en) | A kind of synthetic aperture radar image-forming method and apparatus based on NLFM signal | |
| US5359329A (en) | Jammer reference target measurement system | |
| CN109188387B (en) | Target parameter estimation method for distributed coherent radar based on interpolation compensation | |
| CN105974389B (en) | Machine based on iterative processing sweeps metre wave radar Monopulse estimation method | |
| CN110412514B (en) | Linear frequency modulation continuous wave waveform optimization method under MIMO system | |
| CN106842210A (en) | A kind of new multiple submatrixes synthetic aperture sonar fast imaging algorithm | |
| CN105974390B (en) | Ji Sao metre wave radars mass centre angle-measuring method based on doppler information | |
| CN112834980B (en) | Vortex electromagnetic wave super-resolution direction finding method based on propagation operator | |
| CN110967692A (en) | Imaging method | |
| CN106950565A (en) | Space-borne SAR Imaging jitter compensation method, imaging method | |
| CN109613475A (en) | A kind of arrangement method and target DOA evaluation method of the relatively prime battle array of adaptive optimization | |
| CN111649803B (en) | Three-dimensional radar level meter based on vertical linear array and design method thereof | |
| CN113314832B (en) | Millimeter wave vehicle-mounted MIMO radar antenna array device and design method | |
| CN109490881A (en) | Interference SAR measurement of higher degree system and measurement method based on vortex electromagnetic wave | |
| CN104237886B (en) | High-precision synthetic aperture radar imaging method | |
| CN109738856A (en) | A kind of noise suppressed Capon active target DOA evaluation method based on time reversal | |
| CN110988864A (en) | MTI radar speed measuring method with frequency agility | |
| Feng et al. | MIMO–monopulse target localisation for automotive radar | |
| CN109061599B (en) | STAP method based on cyclostationarity and symmetric prior knowledge | |
| CN111736129B (en) | Space-time modulation bistatic radar system, implementation method and medium thereof | |
| CN108490425A (en) | A kind of angle-measuring method of bistatic MIMO radar | |
| Chen et al. | Multiple-frequency CW radar and the array structure for uncoupled angle-range indication | |
| CN116500620A (en) | Data processing method and device of millimeter wave radar, storage medium and unmanned vehicle | |
| CN110456342A (en) | Single-shot penetrates the far field multiple mobile object detection method of antenna radar | |
| CN112198510B (en) | Method for inhibiting azimuth ambiguity of dual-receiving-channel hybrid full-polarization radar system |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |