CN104007413B - Consider the array position error calibration method of information source azimuthal error - Google Patents
Consider the array position error calibration method of information source azimuthal error Download PDFInfo
- Publication number
- CN104007413B CN104007413B CN201410174084.0A CN201410174084A CN104007413B CN 104007413 B CN104007413 B CN 104007413B CN 201410174084 A CN201410174084 A CN 201410174084A CN 104007413 B CN104007413 B CN 104007413B
- Authority
- CN
- China
- Prior art keywords
- array
- phase
- mrow
- error
- information source
- 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.)
- Expired - Fee Related
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/04—Details
- G01S3/10—Means for reducing or compensating for quadrantal, site, or like errors
-
- 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
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
The invention belongs to array signal processing field, there is the method for the array position error correction of deviation, array position error accurately corrected in more particularly to information source azimuth information.The present invention using single information source correction, but receiving array can precision rotation, the effect of multiple information source independence timesharing corrections is reached, so as to obtain great amount of samples；First use least square fitting estimates the phase of information source steering vector, then rejects the larger array element data of some phase-fitting errors, and then use least square fitting estimation phase again by the way of open country is picked, and obtains information source azimuth information；Sensor position uncertainties are finally corrected using least square method.Using the method for the correction sensor position uncertainties of the present invention, sensor position uncertainties can be accurately corrected very much, method is simple, and be well suited for Practical Project and used.
Description
Technical field
The present invention relates to array signal processing neck, and in particular to a kind of array position error school for considering information source azimuthal error
Correction method.
Background technology
Signal direction of arrival (DOA) estimation has extremely wide application in various fields such as radar, communication, sonars.With
Multiple signal classification (MUSIC) algorithm is the Power estimation algorithm of representative all with very high resolving power and estimated accuracy, its premise
It is accurately known array manifold.By various non-ideal factors (machining error, receiving channel inconsistency, array element mutual coupling etc.)
Influence, array manifold often occurs deviation and disturbance to a certain extent, will make the performance of high resolution remote sensing data algorithm
Severe exacerbation, or even failure.It is limited to current technological level, goes to overcome model error to exist merely in terms of hardware design making
Very big difficulty, the estimation of array error turns into the problem that angular surveying is badly in need of solving.
Existing array correcting method can be generally divided into self-correcting class (see document：Direction finding in
the presence of mutual coupling,Friedlander B,Weiss A J；IEEE Trans.on AP,
1991,39(3):273-284；A Bayesian approach to autocalibration for parametric
Array signal processing, Viberg M；Swindlehurst A L.；IEEE Trans.on SP,1994,42
(12):3495-3507) with active correction class (see document：A kind of array antenna element position, amplitude and phase error it is active
Bearing calibration, Jia Yongkang protects polished, Wu the Huanshui River；Electronic letters, vol, 1996,24 (3):47-52；Sensor array calibration
In the presence of mutual coupling and unknown sensor gains and phases, See C
M S；Electronics Letters,1994,30(5):373-374).Self-correcting class method is generally by the orientation of space information source
Disturbance parameter with array is according to certain majorized function progress Combined estimator.Array calibration can not need the known auxiliary in orientation
Information source, and can online be completed when true bearing is estimated.But, it is often mutual between information source orientation and array error parameter
Coupling so that the identifiability of parameter Estimation hardly results in guarantee (Spatial signature estimation in self-correcting
for uniform linear arrays with unknown receiver gains and phases[J],ASTELY D,
SWINDLEHURST A,OTTERSTEN B.IEEE Transactions on Signal Processing,1999,47(8):
2128-2138).Active correction is by setting the accurately known auxiliary information source in orientation to estimate array disturbance parameter offline
Meter, is method practical at present.
The optimal function that least square method finds data by minimizing the quadratic sum of error is matched, and is a kind of practical number
Learn optimisation technique.When posterior probability density function is the unimodal density function for being symmetrical with Posterior Mean, least mean-square error is estimated
Meter is the best estimate of numerous type cost functions (see document：The statistic mixed-state of signal and estimation theory (second edition), Li Dao
This).The utilization of least square method widely, especially in terms of straight line or curve matching.Least square method also applies to
(see document in the estimation of array element amplitude phase error or site error：Sensor position uncertainties correction and performance point based on least square
Analysis, Yang Zhiwei, Liao Guisheng；System engineering and electronic technology 2007,29 (2):167-169.).
Sensor position uncertainties make array steering vector orientation dependence characteristics occur, therefore the method based on active correction is often
Need multiple independent correction information sources；The accuracy of auxiliary information source position also can impact position error correction result；In addition, working as
Relatively large deviation occurs in information source DOA estimate when some element positions have larger error, and the presence of these problems makes existing
Sensor position uncertainties bearing calibration is difficult to meet the actual demand of engineering.
Notice many phased array arrays be mounted on can with the pedestal of precision rotation, fixed single calibration source and itself
Rotated, it is possible to change the relative angle of front normal and calibration source, be equivalent to a large amount of virtual time-sharing works of acquisition many
Calibration source receives data, and chance is brought to improve sensor position calibration precision.Certainly, the angular accuracy of determining of rotating basis also can shadow
Ring to final correction accuracy, it is necessary to eliminate its influence brought using special measure.
The content of the invention
The technical problems to be solved by the invention are：Propose a kind of array position error correction for considering information source azimuthal error
Method, is accurately corrected to array position error.
The present invention solves the technical scheme that is used of above-mentioned technical problem：
Consider the array position error calibration method of information source azimuthal error, comprise the following steps：
A. by revolving-turret, the array antenna being located on turntable is continuously rotated J angle, obtain under each angle
The sample data that array antenna received is arrived；
B. the phase sequence for obtaining signal source steering vector is calculated according to sample data；
C. the phase sequence of the signal source steering vector obtained to calculating is pre-processed, and to the orientation angles of signal source
It is corrected；
D. sensor position uncertainties are corrected using the orientation angles of the signal source after correction, calculate sensor position uncertainties
Caused phase error；
E. each array element of phase error compensation calculating obtained to array antenna.
Further, step a is specifically included：
The even linear array of M array element is placed in can be on the turntable of precision rotation, revolving-turret, antenna is continuously rotated J
Angle, θj(j=1 ... J), obtains the sample data X of each anglej(t)。
Further, step b is specifically included：
The covariance matrix R for receiving signal is calculated using sample dataj：
Rj=E [Xj(t)Xj H(t)]
The covariance matrix obtained to above-mentioned calculating carries out feature decomposition, obtains estimating for normalized signal source steering vector
EvaluationWherein e is the R corresponding characteristic vector of eigenvalue of maximum, e1For e first element, further calculate
To the phase sequence φ of signal source steering vectorj=angle (a (θj))。
Further, the side that the phase sequence of the signal source steering vector obtained described in step c to calculating is pre-processed
Method is：
φj=[φ1j,φ2j,…,φij,…,φMj] (i=1 ..., M), using first member as reference array element, φ1j=0；
Work as i>When 1, the average of phase intervals is between array element：
Further, it is to the method that the orientation angles of signal source are corrected described in step c：
Phase data to each pretreated passage carries out least square method solution, obtains phase change rate kj：
Further according to phase change rate kjCorrection signal source azimuth angle degree information：
Then the error of fitting of each array element phase is calculated：
φij_error=(2 π (i-1) x0·sinθj/λ-φij)2, then reject by the way of open country is picked some phases and intend
The larger array element data of error are closed, and then phase is estimated using least square fitting again, the side of signal source is calculated
Position
Angle information.
Further, step d is specifically included：
For i-th of array element, with by pretreated phase data and it is corrected after information source orientation angles θj
(j=1 ... J) obtains the coefficient matrix of overdetermined equationWith
Solved using least square method：
Then sensor position uncertainties are calculated：
△xi=xi-(i-1)x0,△yi=yi；
Caused phase error is calculated finally according to sensor position uncertainties：
The beneficial effects of the invention are as follows：Only need to single calibration source, and need not accurately know the azimuth information of the information source,
By the precision rotation of receiving array to obtain mass data sample, the situation of multiple independent source timesharing corrections is equivalent to；This
The least square fitting algorithm that invention first picks open country using band obtains the information source azimuth information of degree of precision, rejects site error larger
Array number avoid its influence to angle estimation according to this, sensor position uncertainties are finally corrected again.This method is simple, error correction
Precision is also higher, and is well suited for Practical Project and uses.
Brief description of the drawings
Fig. 1 is aerial array model schematic in the present invention.
Fig. 2 is that 8 array-element antennas have sensor position uncertainties and the phase curve comparison diagram without sensor position uncertainties.
Fig. 3 is the schematic flow sheet of the present invention.
During Fig. 4 is the present invention, the angle after information source angle correct (adds with real angle i.e. on the basis of desired angle
Angle after 2 degree of random perturbation) difference and real angle and the relativity figure of desired angle difference.
Fig. 5 (a) is that the sensor position uncertainties of information source angle correct and information source angle not timing in the present invention are corrected in X side
To comparison diagram, Fig. 5 (b) is comparison diagram in the Y direction.
Fig. 6 is 7 passage information source angle correct of the invention and information source angle do not correct the phase errors of two kinds of situations with it is true
The contrast curve of phase error.
Fig. 7 be information source angle ideal value of the present invention be 60 ° when, information source angle correct and sensor position uncertainties correction, information source
Angle do not correct and sensor position uncertainties correction and sensor position uncertainties not timing Music spectrum comparison diagram.
Fig. 8 is that information source angle ideal value of the present invention spends the Music angle measurement accuracys that each different angle is obtained 20 to 60
Comparison diagram, compared for information source angle correct and sensor position uncertainties correction, information source angle are not corrected and sensor position uncertainties are corrected
With sensor position uncertainties not three kinds of situations of timing.
Embodiment
For ease of description, the technical term in the present invention is annotated first as follows：
Sensor position uncertainties：Alignment error, measurement error, array working environment and airborne or hip-based platform shake in practice
The bay position caused by factor such as dynamic changes, and causes the mistake existed between preferable element position and actual element position
Difference.
Active correction：Array error parameter is estimated offline by the auxiliary information source for setting orientation accurately known in space
Meter.
DOA estimates：Estimation of Spatial Spectrum main purpose is the variously-shaped array that the multiple sensors of utilization space are constituted,
Parameter and information source position to spatial distribution radiation source are estimated.To signal direction of arrival (DireetionofArrival,
DOA) estimation is the core content of Estimation of Spatial Spectrum.
Steering vector：It is all array elements of array antenna to the response with unit energy arrowband information source.Because array rings
Should be different in different directions, the direction of steering vector and information source is to be mutually related, the unique dependence of this association
In the geometry of array.
The present invention operation principle be：
The even linear array of M array element is in can be on the turntable of precision rotation, array element spacing x0For half-wavelength, in array element far field
In, signal source is at using linear array axis normal as the θ of reference.Using first array element as reference array element, it is assumed that preferable array element position
Put and be arranged in x-axis, then the physical location and ideal position of array element are respectively (xi,yi) and ((i-1) x0, 0) (i=1,2 ...
M), sensor position uncertainties are (△ xi,△yi).The snapshot data received can be expressed as：
X (t)=A (θ) S (t)+N (t), t=1,2 ... K. (1)
S (t) is incoming signal complex envelope, and N (t) is the array noise vector of M × 1, and A (θ) is array received steering vector, K
For fast umber of beats.The covariance matrix R of array is defined as
R=E [X (t) XH(t)]=ARsAH+σ2I (2)
Wherein Rs=E [S (t) SH(t) it is] covariance matrix of signal source.I is unit matrix, σ2For noise power.By son
Space Principles are understood, are carried out feature decomposition to R, can be obtained the estimate of signal source steering vectorWherein e is R's
The corresponding characteristic vector of eigenvalue of maximum, e1For e first element.So can obtain the phase sequence of signal source steering vector
For φ=angle (a (θ)).
As seen from Figure 1, using the first array element as reference position, positioned at (xi,yi) the array element phase of position is：
φi=2 π [xi·sinθ+yi·cosθ]/λ (3)
Wherein φ=[φ1…φM].Similarly, if element position has deviation (△ xi,△yi), then led by position deviation
The phase error of cause can be expressed as:
△φi=2 π [△ xi·sinθ+△yi·cosθ]/λ (4)
Because preferable element position is arranged in x-axis, i.e., the position of array element is in ideal array：
xi0=(i-1) x0, (i=1,2 ... M) (5)
Then have：
φi0=2 π [(i-1) x0·sinθ]/λ (6)
From (6) as can be seen that any array element is the linear function of sine value sin θ relative to the phase of reference array element.If each
The location interval of individual array element is uniform, then phasei0Between interval also will be uniform.If being spaced uneven or array element
There is error in position, then phase becomes uneven, so as to influence DOA angle measurements.As shown in Fig. 2 when sensor position uncertainties are not present,
The phase of 8 antennas is linear relationship, when there are sensor position uncertainties, and the phase of 8 antennas does not just meet linear relationship strictly
, consider and be generally much smaller than array element interval x to sensor position uncertainties △ x, △ y0, therefore it is approximately line that each array element phase sequence, which remains unchanged,
Property form, it is possible to use least square fitting phase straight line, but be due to some sensor position uncertainties significantly greater than its
His array element (4 array elements in such as figure), then in the element position phase deviation straight line farther out, if using least square fitting when
Time can reject the larger array element of these element position deviations, then when the straight line being fitted will be more nearly no sensor position uncertainties
Situation.
Because sensor position uncertainties △ x, △ y are generally much smaller than array element interval x0, each array element phase sequence, which remains unchanged, is approximately
Linear form, you can use formula (7) approximate expression into following form：
φi=k (i-1) (7)
WhereinIt will of course be noted that the span of phase should expand to [- ∞ ∞], i.e.,
It should make the phase of M passage that the trend for being incremented by or successively decreasing substantially is presented by adding and subtracting 360 degree, meet the approximately linear of phase
Change, needs to carry out necessary pretreatment to the original phase data dealt for this.
Because phase approximately meets linear change, then [φ1,φ2,…,φi,…,φM] processing method be：First
Array element is reference array element, φ1=0；Work as i>When 1, the average of phase intervals is between array element：
The data of each passage are substituted into (7), it is possible to obtain over-determined systems.
Then the coefficient matrix of this overdetermined equation is represented by A=[0 1 ... M-1]T, φ=[φ1 φ2…φM]T, so
Least square solution to k is：
And then can obtain：
So the error of fitting of each array element phase can be obtained：
φi_error=(2 π (i-1) x0·sinθ/λ-φi)2(i=1,2 ... M) (11)
If fruit part sensor position uncertainties are than larger, or there is the influence of other interference, part array element there may be than larger
Error of fitting, i.e. part φi_errorThe significantly greater than error of fitting of other array element phases, therefore, when the fitting of some array elements
When application condition is big, a least square can be first carried out with all data and solves operation, by error φi_errorSome very big battle arrays
Metadata reject, by remaining data again in the same way carry out once fitting operation, the k higher with regard to ratio of precision can be obtained,
It further can be obtained by information source angle.
After signal source angular deviation correction, then do sensor position uncertainties correction, the method for being also based on least square.
Because aerial array is placed on graduated rotating disk, information source can be rotated continuously with aerial array normal angle θ, have rotated altogether
J (J >=3) individual angle, only considers the position deviation of single array element, and obtained equation has following shape (12) formula：
φij=2 π [xi·sinθj+yi·cosθj]/λ (12)
Consider that various errors may cause phase change, have：
Wherein i represents array element sequence, and j then represents angle sequence number,It is the phase place change as caused by various errors.
The phase data of different angles is substituted into formula (13), J equation can be obtained altogether：
But unknown number only has three, it is respectively：xi、yiWith △ φi, these three can be just solved using least square method
Parameter.Wherein coefficient matrix is Then least square
Xie Wei：
Then sensor position uncertainties are：
△xi=xi-(i-1)x0,△yi=yi (16)
Then obtaining the phase error as caused by sensor position uncertainties is：
Again △ φijEach array element is compensated, can thus eliminate that sensor position uncertainties bring have impact on.The present invention
In bearing calibration idiographic flow referring to Fig. 3.
Embodiment：
Assuming that the array number M=10 of even linear array, information source goes to 60 ° of degree from 20 ° with aerial array normal angle θ ideals,
2 ° of random perturbations are added in each angle at intervals of 1 °, but during emulation, the signal of signal source transmitting is linear FM signal, is carried
Frequency f0=5300MHz, array element spacing x0For half-wavelength, bandwidth B=500kHz, pulse repetition period T=400us, sample rate fs=
50MHz, signal to noise ratio is 20dB, the sensor position uncertainties △ x=[0.0020-0.0015- of 0 0.0022-0.0021 0.008
0.0023 0.0020], △ y=[0.0021-0.0018 of 0-0.0017 0.0031 0.0015-0.0026 0.0023].
First, information source is with aerial array normal angle θ from 20 ° of scannings to 60 °, and each sweep spacing is 1 °, so as to obtain
Substantial amounts of sample number.But plus 2 ° of random perturbations in each angle during emulation, so the azimuth information of information source is to exist
Deviation.
Secondly, the phase of information source steering vector is obtained, and φ is pre-processed, it is met approximate linear change
Trend.Table 1 describes the result of phase before and after the processing when information source angle is 20 °.
The processing of the information source steering vector phase of table 1
Again, linear fit is carried out to phase using least square method, obtains the corrected value of information source angle, θ, and obtain
The error of fitting of phase, then weeds out error of fitting than larger array element data, then is carried out with remaining data once minimum
Two multiply fitting, obtain the more accurately corrected value of information source angle, θ.
The angle that Fig. 4 compared for after correction (has added 2 ° of random perturbations with real angle on the basis of desired angle
Angle afterwards) difference and real angle and desired angle difference (i.e. 2 ° of random perturbations), from fig. 4, it can be seen that emulating
There is positive and negative 2 degree of deviation in Shi Caiyong real angle and desired angle, and the information source angle, θ after correcting is closer to true angle
Degree, the difference with real angle is basic in the range of 0.1 degree.Therefore, information source azimuth information is more accurate, to follow-up array element position
Put error correction and Music angle measurements advantageously.
Finally, sensor position uncertainties are corrected on the basis of correction information source angle, and with not correcting information source angle
And directly the result that sensor position uncertainties are corrected is compared using step 4.
Fig. 5 compared for information source angle correct and information source angle does not correct the sensor position uncertainties correction of two kinds of situations.From Fig. 5
As can be seen that sensor position uncertainties correction is very accurate after information source angle correct, it is coincide substantially with real sensor position uncertainties,
And the deviation of information source angle not timing, obtained sensor position uncertainties and true sensor position uncertainties is very big.
Fig. 6 compared for 7 passage information source angle corrects and information source angle does not correct the phase error and true phase of two kinds of situations
Position error.
From fig. 6, it can be seen that the phase error after information source angle correct can be very good to fit real phase error.
And believe
The non-timing of source angle, obtained phase error can not reflect the changing rule of true phase error.
When Fig. 7 compared for information source angle ideal value for 60 ° (due to adding random perturbation, real angle is 59.755 °)
Music angle measurements.From figure 7 it can be seen that before not by sensor position uncertainties correction, the angle that Music Power estimations are obtained is
59.3 °, use information source desired angle directly to correct the angle obtained after sensor position uncertainties for 60.6 °, carried using this patent
The angle that method is obtained after being compensated to sensor position uncertainties is 59.8 °, angle measurement accuracy highest, and the peak value of Music spectrums
More sharp, peak value is bigger, to the detection of weak signal target advantageously.Fig. 8 is information source to be obtained in 20 to 60 degree each different angle
Music angle measurement accuracys comparison diagram, compared for not carrying out sensor position uncertainties correction respectively, carry out array element with desired angle
Site error is corrected and this patent institute extracting method carries out the contrast of Music angle measurement accuracys after sensor position uncertainties correction.It is overall next
Say, the ratio of precision obtained using the bearing calibration of the present invention is higher.
Claims (5)
1. consider the array position error calibration method of information source azimuthal error, it is characterised in that comprise the following steps：
A. by revolving-turret, the array antenna being located on turntable is continuously rotated J angle, obtain the array under each angle
The sample data that antenna is received；
B. the phase sequence for obtaining signal source steering vector is calculated according to sample data；
C. the phase sequence of the signal source steering vector obtained to calculating is pre-processed, and the orientation angles of signal source are carried out
Correction；
D. sensor position uncertainties are corrected using the orientation angles of the signal source after correction, calculate sensor position uncertainties and cause
Phase error；
E. each array element of phase error compensation calculating obtained to array antenna；
It is to the method that the orientation angles of signal source are corrected described in step c：
Phase data to each pretreated passage carries out least square method solution, obtains phase change rate kj：
Wherein, A is coefficient matrix；
Further according to phase change rate kjCorrection signal source azimuth angle degree information：
<mrow>
<msub>
<mi>&theta;</mi>
<mi>j</mi>
</msub>
<mo>=</mo>
<mi>argsin</mi>
<mrow>
<mo>(</mo>
<mfrac>
<mi>&lambda;</mi>
<mrow>
<mn>2</mn>
<msub>
<mi>&pi;x</mi>
<mn>0</mn>
</msub>
</mrow>
</mfrac>
<mfrac>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>M</mi>
</munderover>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<msub>
<mi>&phi;</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
</mrow>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>M</mi>
</munderover>
<msup>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
φijIt is i to represent array element sequence, and angle sequence is j array element phase；
Then the error of fitting of each array element phase is calculated：
φij_error=(2 π (i-1) x0.sinθj/λ-φij)2, x0Represent array element interval；Reject some by the way of open country is picked again
The larger array element data of phase-fitting error, and then phase is estimated using least square fitting again, calculate signal source azimuth angle
Spend information.
2. the array position error calibration method of information source azimuthal error is considered as claimed in claim 1, it is characterised in that step
A is specifically included：
The even linear array of M array element is placed in can be on the turntable of precision rotation, and revolving-turret makes antenna continuously rotate J angle, θj
(j=1 ... J), obtains the sample data X of each anglej(t)。
3. the array position error calibration method of information source azimuthal error is considered as claimed in claim 2, it is characterised in that
Step b is specifically included：
The covariance matrix R for receiving signal is calculated using sample dataj：
Rj=E [Xj(t)Xj H(t)]
The covariance matrix obtained to above-mentioned calculating carries out feature decomposition, obtains the estimate of normalized signal source steering vectorWherein e is the R corresponding characteristic vector of eigenvalue of maximum, e1For e first element, further calculate and obtain
The phase sequence φ of signal source steering vectorj=angle (a (θj))。
4. the array position error calibration method of information source azimuthal error is considered as claimed in claim 3, it is characterised in that step
It is to the method that the phase sequence for calculating obtained signal source steering vector is pre-processed described in c：
φj=[φ1j,φ2j,…,φij,…,φMj] (i=1 ..., M), using first member as reference array element, φ1j=0；Work as i>When 1, battle array
The average of phase intervals is between member：
5. the array position error calibration method of information source azimuthal error is considered as claimed in claim 4, it is characterised in that step
D is specifically included：
For i-th of array element, with by pretreated phase data and it is corrected after information source orientation angles θj(j=
1 ... J) obtain the coefficient matrix of overdetermined equationWithProfit
Solved with least square method：
Then sensor position uncertainties are calculated：
Δxi=xi-(i-1)x0,Δyi=yi；xiRepresent the abscissa of i-th of element position, yiRepresent i-th element position
Ordinate；
Caused phase error is calculated finally according to sensor position uncertainties：
Represent the array element phase place change as caused by various errors.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410174084.0A CN104007413B (en) | 2014-04-28 | 2014-04-28 | Consider the array position error calibration method of information source azimuthal error |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410174084.0A CN104007413B (en) | 2014-04-28 | 2014-04-28 | Consider the array position error calibration method of information source azimuthal error |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104007413A CN104007413A (en) | 2014-08-27 |
| CN104007413B true CN104007413B (en) | 2017-08-25 |
Family
ID=51368157
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410174084.0A Expired - Fee Related CN104007413B (en) | 2014-04-28 | 2014-04-28 | Consider the array position error calibration method of information source azimuthal error |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104007413B (en) |
Families Citing this family (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106646345A (en) * | 2016-09-09 | 2017-05-10 | 西北工业大学 | Direction finding error compensation method based on data domain |
| CN108375752B (en) * | 2018-02-05 | 2022-06-24 | 中国人民解放军战略支援部队信息工程大学 | Amplitude-phase error single-radiation-source direction finding method based on all-angle search |
| CN110187304B (en) * | 2019-05-21 | 2021-05-04 | 泰凌微电子(上海)股份有限公司 | Signal arrival angle estimation method and device |
| CN112068072B (en) * | 2020-08-07 | 2023-02-03 | 北京万集科技股份有限公司 | Multi-antenna combined positioning method and device, storage medium and electronic device |
| CN112305526B (en) * | 2020-10-22 | 2023-05-26 | 电子科技大学 | Distributed array system synchronization method based on external calibration source |
| CN113126047B (en) * | 2021-04-15 | 2023-06-16 | 西安电子科技大学 | Self-adaptive channel correction method based on far-field point source |
| CN113518358B (en) * | 2021-04-25 | 2022-05-31 | 长春理工大学 | Method and device for correcting location errors in wireless sensor network deployment |
| CN113438738B (en) * | 2021-08-27 | 2021-10-29 | 网络通信与安全紫金山实验室 | Positioning parameter estimation method, positioning parameter estimation device, computer equipment and storage medium |
| CN113768541B (en) * | 2021-10-27 | 2024-02-13 | 之江实验室 | Method for correcting position error of array element of complex curved surface ultrasonic array transducer |
| CN114814721B (en) * | 2022-03-16 | 2024-10-15 | 西北工业大学 | Acoustic vector array azimuth estimation method based on constraint tensor decomposition |
| CN114814719A (en) * | 2022-05-06 | 2022-07-29 | 成都天奥测控技术有限公司 | Orientation estimation method, module, equipment and system based on adaptive feedback |
| CN115508775B (en) * | 2022-10-20 | 2024-10-15 | 电子科技大学 | Node positioning method using incoming wave azimuth difference measurement |
| CN116760437B (en) * | 2023-08-17 | 2023-10-20 | 四川省华盾防务科技股份有限公司 | Broadband microwave receiving and transmitting control method and system applied to phased array system |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101644765B (en) * | 2009-04-23 | 2012-03-14 | 中国科学院声学研究所 | Amplitude and phase error correction method used for linear array of underwater acoustic transducer |
| CN102621532B (en) * | 2012-03-31 | 2013-05-22 | 华中科技大学 | Synthetic aperture radiometer visibility phase error correction method based on array rotation |
-
2014
- 2014-04-28 CN CN201410174084.0A patent/CN104007413B/en not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| CN104007413A (en) | 2014-08-27 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104007413B (en) | Consider the array position error calibration method of information source azimuthal error | |
| CN104297718B (en) | Interferometer array integrated correction method | |
| Davis et al. | Angle estimation with adaptive arrays in external noise fields | |
| CN104076348B (en) | Radar beyond visual range base line passive cooperative localization method | |
| CN109270486B (en) | Two-dimensional direction finding method for frequency-varying signal by adopting rotation of two sensors | |
| CN109507635A (en) | Utilize the array amplitude phase error evaluation method of two unknown orientation auxiliary sources | |
| CN109856605A (en) | A kind of while formation of the digital multiple beam quadratic fit curve is directed toward modification method | |
| CN110196414A (en) | A kind of antenna beam pointing method based on compensation Antenna error | |
| CN115061156A (en) | Array antenna satellite navigation deception resisting method and system based on integrated navigation | |
| CN109975840B (en) | Positioning correction method for zero-setting antenna satellite navigation receiver | |
| CN103616685B (en) | Based on the ISAR image geometry calibrating method of characteristics of image | |
| CN104181513A (en) | Array element position correcting method of radar antenna | |
| CN108983169B (en) | Meter wave radar terrain correction method based on digital elevation model | |
| CN117724035A (en) | Interferometer direction finding positioning method based on two-stage correction | |
| CN109324309A (en) | A kind of rotation single antenna angle-measuring equipment and its measurement method | |
| RU2610150C1 (en) | Method of determining ground radio-frequency source coordinates when performing on-board radio-direction finding | |
| Yuan et al. | Correction method for magnitude and phase variations in acoustic arrays based on focused beamforming | |
| CN112240957B (en) | Method for correcting amplitude-phase characteristics of antenna in satellite navigation interference direction finding | |
| RU2711341C1 (en) | Two-dimensional direction finding method | |
| CN102222230B (en) | Multi-input multi-output imaging radar target identification method based on parameterization dictionary | |
| CN108387868B (en) | Scanning type radiation source positioning method based on signal arrival angle difference measurement | |
| KR101610051B1 (en) | A radio direction finder | |
| CN109856619B (en) | Radar direction finding relative system error correction method | |
| CN102305927A (en) | 2-dimensional (2-D) radar error registration two-step algorithm under earth centered earth fixed (ECEF) coordinate system | |
| CN112697129B (en) | Quick method for determining two-dimensional attitude of antenna linear array by utilizing one-dimensional angle measurement |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170825 Termination date: 20180428 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |