CN103903609B - A kind of circular array design of beamformer with constant beamwidth method - Google Patents
A kind of circular array design of beamformer with constant beamwidth method Download PDFInfo
- Publication number
- CN103903609B CN103903609B CN201410146364.0A CN201410146364A CN103903609B CN 103903609 B CN103903609 B CN 103903609B CN 201410146364 A CN201410146364 A CN 201410146364A CN 103903609 B CN103903609 B CN 103903609B
- Authority
- CN
- China
- Prior art keywords
- array
- wave beam
- beamformer
- weight vector
- frequency
- 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
Links
- 239000011159 matrix material Substances 0.000 claims abstract description 14
- 230000017105 transposition Effects 0.000 claims description 4
- 230000021615 conjugation Effects 0.000 claims description 3
- 125000004122 cyclic group Chemical group 0.000 claims description 3
- 230000015572 biosynthetic process Effects 0.000 description 15
- 230000002194 synthesizing Effects 0.000 description 13
- 238000003786 synthesis reaction Methods 0.000 description 12
- 238000005755 formation reaction Methods 0.000 description 3
- 238000000034 method Methods 0.000 description 3
- 238000001514 detection method Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000004891 communication Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 230000000875 corresponding Effects 0.000 description 1
- 230000003247 decreasing Effects 0.000 description 1
- 230000001419 dependent Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 230000004304 visual acuity Effects 0.000 description 1
Abstract
The present invention relates to a kind of circular array design of beamformer with constant beamwidth method, propose a kind of simple and accurate design of beamformer with constant beamwidth method for uniform rings shape array.The method utilizes the characteristic of circular matrix, gives the accurate solution of ring array least squares error beam pattern synthtic price index, and the wave beam finally synthesized and least squares error are expressed simply as subcomponent superposition sum.When expecting that wave beam is converted to appropriate format, the closed expression of available beamformer with constant beamwidth weight vector.Instant invention overcomes prior art operation complicated and accurate not enough not: optimal solution be exactly represented as subcomponent superposition and, and have the expression formula of enclosed.
Description
Technical field
The invention belongs to a kind of Beam-former method for designing, relate to a kind of circular array beamformer with constant beamwidth
Method for designing, it is adaptable to the broadband target detection of circular array and broadband target Bearing Estimation, belongs to marine acoustics, array letter
Number process and the field such as sonar technique.
Background technology
Circular array is simple due to formation, do not have port and starboard ambiguity and can all-round in the range of form substantially constant
Wave beam, be therefore widely used in fields such as sonar, radar, communication and voice engineerings.Signal processing side about this array
Method also emerges in an endless stream, and wherein Broadband Beamforming Method has been a great concern especially.
Owing to the wave beam of conventional arrays output changes with frequency, especially when signal is from other direction incidence of main lobe,
Different frequency signals frequency spectrum will be distorted, and as having carried out low-pass filtering, is unfavorable for subsequent treatment.For solving this
Problem, needs design beamformer with constant beamwidth, i.e. wave beam not with frequency shift, particularly ensures that main lobe is in given bandwidth
Interior invariable.This type of method has been proposed that a lot, is based primarily upon two class thoughts: (1) utilizes Wave beam forming medium frequency and aperture
Relation, invariable by changing the effective aperture wave beam that makes different frequency corresponding of array, typically have document 1
“Constant beamwidth receiving arrays for broad band sonar systems,Acustica,
1970, vol.23 (1), p.21-26 " disclosed combination battle array method, but the method requires that formation strictly meets expansion structure, basic matrix
Size is relatively big, needs array element number many.(2) wave beam needing synthesis is made to approach, such as with expectation wave beam under certain criterion
Document 2 " arbitrary geometry array wideband constant beamwidth Wave beam forming new method, acoustic journal, 2001, vol.26 (1), p.55-58 "
Disclosed Bessel series method and document 3 " any sensor array time domain constant beam-width wave beam based on Second-order cone programming
Formed, acoustic journal, 2005, vol.30 (4), p.309-316 " disclosed in convex optimization method.Although this kind of method is the suitableeest
For any formation, but they fail to make full use of the characteristic of circular array, do not provide Exact.Document 4 " On
the design of digital broadband beamformer for uniform circular array with
frequency invariant characteristics,IEEE International Symposium on Circuits
And Systems (ISCAS), 2002, vol.1, p.I-693-I-696 " it is specifically designed for ring array and discloses a kind of simple permanent
Determining beam-width beam forming method, ultimate principle is to utilize phase place modal theory, not amount with frequency shift is separated, then
Matching expectation wave beam, obtains required constant beam-width wave beam.Although the method is simple, but phase place modal theory relates to spatial domain adopts
Sample and truncated series, the final result obtained is not accurate.
Summary of the invention
Solve the technical problem that
In place of the deficiencies in the prior art, the present invention proposes a kind of circular array beamformer with constant beamwidth
Method for designing, solves prior art operation complicated and accurate not enough.
Technical scheme
A kind of circular array design of beamformer with constant beamwidth method, it is characterised in that step is as follows:
Step 1: determine that expectation wave beam is the two dimensional beam in horizontal extent, by circular array at reference frequency theory meter
Obtaining, its expression formula is:Wherein M is element number of array, ()*Represent and ask altogether
Yoke, ps(ka, φ) is the unit amplitude plane wave signal incident from direction φ that the s array element receives, ωm(krar) for giving
Fixed expectation weight vector element, and need to meetM=1,2 ..., M/2-1;
DescribedIn formula:arFor the radius of circle of reference annular array, kr=2 π/
λr, λrRepresent the wavelength with reference to incident plane wave, φs=s β, β=2 π/M;
Described
Step 2: calculate expectation wave beam and synthesize the least mean-square error of wave beam:
Wherein:And λmIt is respectively circular matrix ρ in reference frequency krarWith eigenvalue during frequency ka to be calculated,For Cyclic Moment
Battle arrayEigenvalue, meet relation With
Described
Described
Described
In formula: J0() is the 0th rank cylindricality Bessel function, Δ rsr=2arSin (s β/2) and Δ rs=2asin (s β/
2) it is with reference to the distance between ring array and ring array to be studied each the m and m' sensor, s=| m-m ' | respectively;
Step 3: calculate beamformer with constant beamwidth weight vector in the frequency band range determined, the of weight vector
M element expression is:
Composition weight vector is ω=[ω0,ω1,…,ωM-1]T, subscript T represents transposition;
Step 4: according to the multi-form of weight vector, uses following two method to synthesize final wave beam:
1, directly calculate:
2, Element space calculates:
Wherein w=V ω, V=[v0,v1,…,vM-1], P (ka, φ)=[p0(ka,φ),p1(ka,φ),…,pM-1(ka,
φ)]TFor array manifold vector, vm=M-1/2[1 eimβ … ei(M-1)mβ]TFor circular matrix characteristic vector and meet relationKa is value in the frequency range of wave beam, and this frequency range is according to design accuracy, step 2 minimum obtained is all
Side curve of error δ (ka) determines.
Beneficial effect
A kind of circular array design of beamformer with constant beamwidth method that the present invention proposes, for uniform rings shape battle array
Row propose a kind of simple and accurate design of beamformer with constant beamwidth method.The method utilizes the characteristic of circular matrix, gives
Go out the accurate solution of ring array least squares error beam pattern synthtic price index, and the wave beam finally synthesized and least square have been missed
Difference is expressed simply as subcomponent superposition sum.When expecting that wave beam is converted to appropriate format, available constant beam-width wave beam shape
Grow up to be a useful person the closed expression of weight vector.Instant invention overcomes prior art operation complicated and accurate not enough: optimal solution
Be exactly represented as subcomponent superposition and, and have the expression formula of enclosed.
Have the beneficial effect that:
1. the present invention is without the biggest array aperture, and array number does not the most require a lot, simpler than the method for document 1.
2. The present invention gives the accurate solution making synthesis wave beam minimum with the square error of expectation wave beam, principle is not made
Any approximation, is not artificially induced other error, more more reliable than the mode territory method of document 4.
3. when expectation wave beam be transformed to appropriate format, the present invention obtained beamformer with constant beamwidth weight vector and
The closed expression of final error, and other document is all difficult to obtain such result.
Accompanying drawing explanation
Fig. 1: circular array schematic diagram.
Fig. 2: expectation beam pattern.
Fig. 3: expectation wave beam and the least squares error variation relation with frequency synthesizing wave beam.
The beam pattern of the constant beam-width broadband beams figure of Fig. 4: synthesis is with the change of frequency.
The Overlapping display of the constant beam-width broadband beams figure of Fig. 5: synthesis.
Fig. 6: the directional gain of the conventional constant beam-width wave beam with synthesis.
Detailed description of the invention
In conjunction with embodiment, accompanying drawing, the invention will be further described:
The present invention is applicable to the design of beamformer with constant beamwidth method of uniform rings shape array, gives uniform rings
The accurate solution of the least squares error beam pattern synthtic price index of battle array, and expectation wave beam is converted to appropriate format, obtain constant
The closed expression of beam-width beam shaper weight vector.Its process is:
(1) given expectation wave beam.Owing to circular array more attention is the performance of Wave beam forming in horizontal plane, so
Expect the two dimensional beam that wave beam is only taken as in horizontal extent.Expect that wave beam is obtained in reference frequency Theoretical Calculation by circular array
Arriving, its expression formula is:
Wherein M is element number of array, ()*Represent and ask conjugation, ps(krar, φ) be m-th array element receive from direction φ
Incident unit amplitude plane wave signal, expression formula is:
In formulaarFor the radius of circle of reference annular array, kr=2 π/λr, λrRepresent the ripple with reference to incident plane wave
Long.Additionally have:
β=2 π/M, φs=s β.ωm(krar) it is given expectation weight vector element, and need to meet
(2) relation that least squares error changes is calculated with frequency, to determine applicable according to required synthesis precision
Frequency band range.Expect that wave beam with the least squares error of synthesis wave beam is:
WhereinAnd λmIt is respectively circular matrix ρ in reference frequency krarWith eigenvalue during frequency ka to be calculated,For
Circular matrixEigenvalue, be satisfied by relationλm=λM-mWithCalculating formula is followed successively by:
J in formula0() is the 0th rank cylindricality Bessel function, Δ rsr=2arSin (s β/2) and Δ rs=2asin (s β/2)
It is with reference to the distance between ring array and ring array to be studied each the m and m' sensor, s=| m-m ' | respectively.
(3) beamformer with constant beamwidth weight vector in the frequency band range determined is calculated.Weight vector m-th unit
The expression formula of element is:
Weight vector is ω=[ω0,ω1,…,ωM-1]T, subscript T represents transposition.
(4) synthesize final wave beam, according to the multi-form of weight vector, be divided into following two situation:
A) direct formula for calculating is:
Wherein E=[E0,E1,…,EM-1]T, subscript H represents conjugate transpose;
B) Element space calculating formula is:
Wherein w=V ω, V=[v0,v1,…,vM-1], P (ka, φ)=[p0(ka,φ),p1(ka,φ),…,pM-1(ka,
φ)]TFor array manifold vector, vm=M-1/2[1 eimβ … ei(M-1)mβ]TFor circular matrix characteristic vector and meet relationKa is value in the frequency range of wave beam, and this frequency range is according to design accuracy, step 2 minimum obtained is all
Side curve of error δ (ka) determines.
Noting, a) result with b) two kinds of method synthesis is consistent, and first method is relatively simple, and second method is more
It is beneficial to practical operation.
Specific embodiment is as follows:
(1) with reference to Fig. 1 and 2.This circular array comprises M equally distributed array element.Due to the more pass of circular array
Note is the performance of Wave beam forming in horizontal plane, it is desirable to wave beam is only taken as the two dimensional beam in horizontal extent.Expect wave beam
Being obtained in reference frequency Theoretical Calculation by circular array, its expression formula is:
Wherein M is element number of array, ()*Represent and ask conjugation, ps(krar, φ) be m-th array element receive from direction φ
Incident unit amplitude plane wave signal, expression formula is:
In formulaarFor the radius of circle of reference annular array, kr=2 π/λr, λrRepresent the ripple with reference to incident plane wave
Long.Additionally have:
β=2 π/M, φs=s β.ωm(krar) it is given expectation weight vector element, and need to meet
Assuming to expect that wave beam is calculated by 16 yuan of uniform rings battle arrays, reference frequency is krar=3, its weight vector element
ω0~ω8It is listed in table 1, notices that other element is by relationDirectly obtain.The expectation ripple obtained by formula (11)
The side lobe levels that is mainly characterized by of bundle figure is below-20dB, and the linear trend successively decreased.
Table 1
ω0 | ω1 | ω2 | ω3 | ω4 | ω5 | ω6 | ω7 | ω8 |
-0.0721 | -O.O548i | -0.0365 | 0.0538i | 0.1129 | -O.3079i | -0.9794 | 4.2747i | 25.9878 |
(2) with reference to Fig. 3.Calculate the relation that changes with frequency of least mean-square error, with true according to required synthesis precision
The fixed frequency band range being suitable for.
Expect that wave beam is expressed as δ=< B with the square error of design wave beamd(φ)-B(φ),Bd(φ)-B (φ) >, symbol
<>representsIt is derived by further:
WhereinAnd λmIt is respectively circular matrix ρ in reference frequency krarWith eigenvalue during frequency ka to be calculated,
For circular matrixEigenvalue, be satisfied by relationλm=λM-mWithCalculating formula is followed successively by:
J in formula0() is the 0th rank cylindricality Bessel function, Δ rsr=2arSin (s β/2) and Δ rs=2asin (s β/2)
It is with reference to the distance between ring array and ring array to be studied each the m and m' sensor, s=| m-m ' | respectively.By
Formula (14) is known, square error is solely dependent upon expectation weight vector and the eigenvalue of associated cyclic matrix, and equal to M sub-error
Sum.
A minimum is had with frequency variation curve at ka=3, meaning by formula (14) calculated least squares error
Now synthesis wave beam to be completely superposed with expectation wave beam.On minimizing right side, error raises with frequency and becomes big, left side then phase
Instead.If ensureing that the error of synthesis is less thanThen suitably frequency band range is ka ∈ [1,6].
(3) beamformer with constant beamwidth weight vector in the frequency band range determined is calculated.Weight vector m-th unit
The expression formula of element is:
Weight vector is ω=[ω0,ω1,…,ωM-1]T, subscript T represents transposition.
(4) synthesize final wave beam, according to the multi-form of weight vector, be divided into following two situation:
A) direct formula for calculating is:
Wherein E=[E0,E1,…,EM-1]T, subscript H represents conjugate transpose;
B) Element space calculating formula is:
Wherein w=V ω, V=[v0,v1,…,vM-1], P (ka, φ)=[p0(ka,φ),p1(ka,φ),…,pM-1(ka,
φ)]TFor array manifold vector, vm=M-1/2[1 eimβ … ei(M-1)mβ]TFor circular matrix characteristic vector and meet relationKa is value in the frequency range of wave beam, and this frequency range is according to design accuracy, step 2 minimum obtained is all
Side curve of error δ (ka) determines.
Noting, a) result with b) two kinds of method synthesis is consistent, and first method is relatively simple, and second method is more
It is beneficial to practical operation.
With reference to Fig. 4 and Fig. 5 and Fig. 6.Consider 16 yuan of uniform rings battle arrays, formula (19) wave beam synthesized is with expectation wave beam all
Coincideing very well, their beam angle is constant, especially main lobe region.It is emphasized that obtain here is constant
Beam-width beam is minimum with the square error of expectation wave beam in all getable wave beams.It addition, the finger of constant beam-width wave beam
Directional index is all at about 12dB, much larger than conventional method in given frequency range.This means that the present invention is had
There is the constant beam-width wave beam of super directivity, significant for the target detection probability and resolving power improving system.
Claims (1)
1. a circular array design of beamformer with constant beamwidth method, it is characterised in that step is as follows:
Step 1: determine that expectation wave beam is the two dimensional beam in horizontal extent, circular array obtain in reference frequency Theoretical Calculation
Arriving, its expression formula is:Wherein M is element number of array, ()*Represent and seek conjugation,
pn(krar, φ) and it is the unit amplitude plane wave signal incident from direction φ that receive of the n-th array element, ωm(krar) it is given
Expectation weight vector element, and need to meetM=1,2 ..., M/2-1;
DescribedIn formula:arFor the radius of circle of reference annular array, kr=2 π/λr, λr
Represent the wavelength with reference to incident plane wave, φn=n β, β=2 π/M;
Described
Step 2: calculate expectation wave beam and synthesize the least mean-square error of wave beam:
Wherein:And λmIt is respectively circular matrix ρ in reference frequency krarWith eigenvalue during frequency ka to be calculated,For Cyclic Moment
Battle arrayEigenvalue, meet relationλm=λM-mWith
Describedρs(krar)=J0(kr·Δrsr);
Describedρs(ka)=J0(k·Δrs);
Described
In formula: J0() is the 0th rank cylindricality Bessel function, Δ rsr=2arSin (s β/2) and Δ rs=2a sin (s β/2) point
It not with reference to the distance between ring array and ring array to be studied each the m and m' sensor, s=| m-m ' |;
Step 3: calculate beamformer with constant beamwidth weight vector in the frequency band range determined, the m-th of weight vector
Element expression is:
Composition weight vector is ω=[ω0,ω1,…,ωM-1]T, subscript T represents transposition;
Step 4: according to the multi-form of weight vector, uses following two method to synthesize final wave beam:
1, directly calculate:
2, Element space calculates:
Wherein w=V ω, V=[v0,v1,…,vM-1], P (ka, φ)=[p0(ka,φ),p1(ka,φ),…,pM-1(ka,φ)]T
For array manifold vector, vm=M-1/2[1 eimβ … ei(M-1)mβ]TFor circular matrix characteristic vector and meet relationKa is value in the frequency range of wave beam, and this frequency range is according to design accuracy, step 2 minimum obtained is all
Side curve of error δ (ka) determines.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410146364.0A CN103903609B (en) | 2014-04-11 | A kind of circular array design of beamformer with constant beamwidth method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410146364.0A CN103903609B (en) | 2014-04-11 | A kind of circular array design of beamformer with constant beamwidth method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103903609A CN103903609A (en) | 2014-07-02 |
CN103903609B true CN103903609B (en) | 2016-11-30 |
Family
ID=
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104730491B (en) | A kind of virtual array DOA estimation method based on L-type battle array | |
Nannuru et al. | Sparse Bayesian learning for beamforming using sparse linear arrays | |
CN103018730B (en) | Distributed sub-array wave arrival direction estimation method | |
Ma et al. | Theoretical and practical solutions for high-order superdirectivity of circular sensor arrays | |
CN101644773B (en) | Real-time frequency domain super-resolution direction estimation method and device | |
CN101995574B (en) | Near field focusing beam forming positioning method | |
CN101813772B (en) | Array beamforming method by quickly expanding and dragging broadband frequency domain | |
CN103630887B (en) | The multi-beam sonar of a kind of Multibeam synthesis method and use the method | |
CN103902830B (en) | A kind of sane Sidelobe control of circular array surpasses directional wave beam forming method | |
CN103353596A (en) | Wave beam space domain meter wave radar height measurement method based on compressed sensing | |
CN104166120B (en) | A kind of acoustic vector justifies battle array robust wideband MVDR direction estimation methods | |
JP2013174498A (en) | Monopulse angle measurement apparatus and monopulse angle measurement method | |
CN105182285A (en) | Target direction-finding method based on acoustic vector two-dimensional nested array | |
CN105005038A (en) | Improved acoustic vector array coherent source DOA estimation algorithm | |
CN104931923A (en) | Grid iterative estimation of signal parameters via rotational invariance techniques (ESPRIT), namely, extensible rapid estimation algorithm capable of being used for uniform circular array 2-dimensional direction of arrival (2D DOA) | |
CN107170441B (en) | Circular ring array optimal frequency invariant response super-directivity beam forming method | |
CN109085556A (en) | A kind of high-frequency ground wave radar wave field forming method based on a second order p-ratio | |
CN105277932B (en) | A kind of Doppler correction method in the Wave beam forming based on down coversion | |
CN107238829B (en) | Double-layer circular ring array stable super-directivity beam forming method | |
CN108595758B (en) | Method for synthesizing optimal broadband beam pattern of sensor array in any form | |
CN103837871A (en) | Inverse beamforming method and system | |
CN103903609B (en) | A kind of circular array design of beamformer with constant beamwidth method | |
CN101825715A (en) | Method for forming wave beam of satellite navigation receiving equipment | |
CN109375197B (en) | Small-size vector array low-frequency scattering correction method | |
CN109061597B (en) | Ionospheric clutter suppression method based on blind source separation and time-frequency ridge domain filtering |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant |