CN1625088A - Blind estimating method for error of array antenna channel - Google Patents
Blind estimating method for error of array antenna channel Download PDFInfo
- Publication number
- CN1625088A CN1625088A CN 200310111076 CN200310111076A CN1625088A CN 1625088 A CN1625088 A CN 1625088A CN 200310111076 CN200310111076 CN 200310111076 CN 200310111076 A CN200310111076 A CN 200310111076A CN 1625088 A CN1625088 A CN 1625088A
- Authority
- CN
- China
- Prior art keywords
- array antenna
- theta
- signal source
- error
- data
- 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
Images
Landscapes
- Variable-Direction Aerials And Aerial Arrays (AREA)
Abstract
This invention discloses a blind estimation method of array antenna channel error. It adopts a marrow-band signal resource used for correction, changes the azimuth of the signal resource through moving the movable signal resource, and gets the data of the signal resource pre-motion and after-motion. Then adopt MUSIC computation method to roughly estimate the azimuth of the signal resource pre-motion and after-motion, define a target function, get the accurate azimuth of the signal resource pre-motion and after-motion through getting the extremum of the target function, and to obtain the amplitude and phase errors of every channel of the array antenna.
Description
Technical field:
The invention belongs to communication technical field, it is particularly related to the estimation technique of array antenna channel error.
Background technology:
The growth of multimedia mobile communication business has improved the requirement to the communication system data transmission rate, for this reason, needs to solve some mis-behave problems that cause owing to reasons such as multipath fading, common-channel interference.(SDMA) can effectively address this problem to division multiple access technology, and the core of SDMA technology is an array antenna, and the error of array antenna can have a strong impact on the performance of SDMA technology, also has identical problem in field of radar.
Based on the high-resolution DOA algorithm for estimating of characteristic value decomposition,, have a wide range of applications in fields such as radar, sonar, mobile communication as the MUSIC algorithm.Theory analysis and evidence, this class algorithm is under desirable array bar spare, and performance is fine, but this class algorithm is all very sensitive to noise disturbance and systematic error, and the performance of noise disturbance and this class algorithm of systematic error meeting severe exacerbation descends its resolution.When noise disturbance mainly showed little number of samples, low signal-to-noise ratio, there was a thresholding in this class algorithm, and when signal to noise ratio was lower than this thresholding, the resolution capability of algorithm sharply descended.Noise disturbance can reduce its influence with the increasing number of samples usually.Systematic error can make a big impact to the performance of algorithm equally.Strengthen the systematic error that number of samples can not reduce array.Propose a lot of algorithms at present and addressed this problem, can see in the literature, mainly contained three class bearing calibrations at present.The first kind is to place a signal source (to see document: Matrix-constructioncalibration method for antenna arrays, Hung, E.K.L. in the far field of array antenna to be corrected; Aerospace and ElectronicSystems, IEEE Transactions on, Volume:36 Issue:3, July 2000, Page (s): 819-828, " the matrix correction method of array antenna "), its azimuth is known, uses the signal of this signal source emission that array antenna is proofreaied and correct.This class methods superior performance, but it is very big to implement difficulty in actual environment.Second class methods need two RF signal sources, its azimuth can be unknown, (see document: Array antenna errorscalibration using signal sources in unknown directions, Jia Yongkang but the differential seat angle between the two must be known; Bao Zheng; Wu Huan; Radar, 1996.Proceedings., CIE International Conference of, 8-10 Oct.1996, Page (s): 503-506 " the error in array antenna bearing calibration of the unknown signaling angle of arrival ").The shortcoming of this method also is that the enforcement difficulty is big.The 3rd class be correction signal with power splitter be separated into M road homophase, etc. the signal of amplitude, directly inject each array element passage, analyze the output signal of array, then the channel error of correction array antenna.These class methods can effectively solve the channel transfer error problem, but have improved the hardware complexity of system, and will do change on the hardware to system.
The MUSIC algorithm: the core concept of MUSIC algorithm is exactly that the covariance matrix of array output signal is made feature decomposition, feature space be divided into signal subspace and noise subspace.Signal subspace and noise subspace quadrature define MUSIC spectrum, search speed spectrum peak, the angle of spectrum peak correspondence be the direction of arrival of incoming signal [see document: open the prominent personage, protect polished, " Array Signal Processing ", 2000.12].
Summary of the invention:
The invention provides a kind of blind estimating method of array antenna channel error, adopt method of the present invention accurately to estimate the channel amplitude sum of errors phase error of array antenna not needing measurement update with under azimuthal situation of radio-frequency signal source.
For convenience of description, at first define:
Term definition:
Far field condition: generally be meant 10 times of wavelength that the distance between signal source and the array antenna transmits greater than signal source.
Direction of arrival (DOA): the angle between the line between the reference array element of signal source and array antenna and the normal of array antenna.
The blind estimating method of a kind of array antenna channel error of the present invention is characterized in that adopting following step:
Step 1 (signalization source):
And array antenna between distance satisfy under the situation of far field condition, place a narrow band signal source s (t), the direction of arrival of incoming signal (signal source emission signal) is θ
1, as shown in Figure 1.The start signal source is in the data of the individual snap of the output collection N (N>400) of array antenna, shutdown signal source then.The k that collects fast beat of data is:
x
1k=[x
1k(1)…x
1k(M)]
H (1)
X in the formula
1k(i) in N the snap that obtains before the expression displaced signal sources in k snap, i the data that array element is exported.()
HExpression Hermitian transposition.
Step 2 (displaced signal sources):
Displaced signal sources s (t) is placed on the position that another one satisfies far field condition, and this moment, the azimuth of signal source was θ
2(θ
1≠ θ
2).Start signal source once more is in the data that the output of array antenna is gathered the individual snap of N (N>400), shutdown signal source.The k that collects fast beat of data is:
x
2k=[x
2k(1)…x
2k(M)]
H (2)
X in the formula
2k(i) in N the snap that obtains after the expression displaced signal sources in k snap, i the data that array element is exported;
Step 3 (data are synthetic):
Data according to twice collection before and after moving obtain two covariance matrix R
1And R
2
Step 4 (asking intermediate variable):
To R
1And R
2Make feature decomposition respectively, can obtain the characteristic vector u of two big characteristic values of correspondence
M1And u
M2, these two characteristic vectors are obtained two intermediate variables with its first element do normalization,
Step 5 (estimation azimuth):
Adopt MUSIC algorithm estimated signal source to move the azimuth of front and back, just can obtain
With
Step 6 (DOA estimation domain):
Define two independent variable territories, with
With
Being the center, is that radius is expanded to both sides with δ (desirable 0.3 degree of reference value), obtains two DOA estimated value intervals
With
Step 7 (the estimation flow process of array antenna channel error)
At first, the phase error of k the passage of data rough estimate array antenna that obtains before moving according to signal source:
In the formula (7), the angle of ∠ () operator representation calculated complex, d is the spacing of adjacent array element, λ is the wavelength of incoming signal.
The phase error of k the passage of data rough estimate array antenna that obtains after moving according to signal source:
Then, introduce a parameter:
At last, determine the true incidence angle of signal source:
And then according to formula (9) calculating ξ (0); Get suitable step-length ε (reference value gets 0.01).Make l=0.
Ψ
1(l)=Ψ
1(0)+and l ε bl=l+1, Ψ
2(l)=Ψ
2(0)+and l ε, calculate ξ (l);
If c
Then change b; D seeks the minimum of ξ.If minimum does not exist, then suitably increase δ (reference value is got δ=δ+0.5), return a; The angle value Ψ of the minimum correspondence of e ξ
1(l) and Ψ
2(l) be exactly real azimuth angle theta before and after signal source moves
1, θ
2
Step 8 (range error of computing array antenna channels and phase error)
The u that obtains according to step 4
S1, the range error of k passage of computing array antenna is estimated
For:
In the formula (10), u
S1(k) expression vector u
S1In k element
The u that obtains according to step 4 again
S1The θ that obtains with step 7
1, k channel phases estimation error of computing array antenna
For:
Through after the above-mentioned steps, just can accurately estimate the channel amplitude error of array antenna
And phase error
Operation principle of the present invention is:
Covariance matrix to the equidistant linear array antenna output signal of M unit carries out feature decomposition, with dominant eigenvalue characteristic of correspondence vector u
m, serve as that reference obtains do normalization with first element,
The direction vector of array antenna is a (θ),
In the formula (13), θ is the angle of arrival of incoming signal, and d is the spacing of adjacent array element, and λ is the wavelength of incoming signal.
When there was channel amplitude sum of errors phase error in array antenna, the direction vector of incoming signal was:
a′(θ)=∧a(θ)
ρ in the formula
KdThe amplitude gain of representing k passage, φ
kThe phase error of representing k passage.Can prove:
u
s=a′(θ) (15)
Then have
Following formula is got phase information, then being estimated as of the phase error of k passage of array antenna:
To the following formula delivery, being estimated as of the range error of k passage of array antenna then:
In sum, the present invention use the narrow band signal source by adopting a correction, and mobile displaced signal sources changes the azimuth of signal source, before move in the picked up signal source and respectively one group of data after moving.And then adopt MUSIC algorithm rough estimate signal source to move the azimuth of front and back, define a target function, by asking the extreme value of target function, the accurate azimuth of front and back is moved in the picked up signal source, just can obtain the range error and the phase error of each passage of array antenna.
In the channel error alignment technique of array antenna, a correction radio-frequency signal source is set in the place of the condition that satisfies the far field, and correction accuracy is relevant with the certainty of measurement of the direction of arrival of radio-frequency signal source with correction.This algorithm can accurately be estimated the channel amplitude sum of errors phase error of array antenna not needing measurement update with under azimuthal situation of radio-frequency signal source.
Description of drawings:
Fig. 1 be among the present invention signal source schematic diagram is set
Wherein, s (t) represents testing source, θ
1The expression signal source just incides the direction of arrival of the signal of array antenna, in the array antenna shown in the figure with respect to the azimuth of array antenna, 1 represents first array element, and m represents m array element, and m (1≤m≤M), M represents M array element, and M is a natural number.
Fig. 2 is the schematic diagram that is provided with after signal source moves among the present invention
Wherein, s (t) represents testing source, θ
2The expression signal source just incides the direction of arrival of the signal of array antenna with respect to the azimuth of array antenna.In the array antenna shown in the figure, 1 represents first array element, m (1≤m≤M) represent the array element, M to represent M array element.M is a natural number.
Fig. 3 is the actual error table of array antenna passage
Wherein, K represents the channel position of array antenna, ρ
kThe physical channel range error of k passage of expression array antenna, φ
kThe physical channel phase error of k passage of expression array antenna.
Fig. 4 is the array antenna channel error result who adopts this algorithm to estimate
Wherein, K represents the channel position of array antenna, ρ
kThe range error of k the passage of array antenna that expression employing method of the present invention obtains is estimated φ
kThe phase error estimation and phase error of k the passage of array antenna that expression employing method of the present invention obtains.
Embodiment:
The blind estimation of 8 element array antenna channels errors:
We adopt equidistantly line array of 8 (M=8) unit, and array element distance is a half-wavelength.In the placed around of aerial array a signal source, and the distance between the array satisfies far field condition, the signal to noise ratio of signal is 20dB.Azimuth before and after signal source moves is respectively :-10 °, 20 °, the fast umber of beats of image data is 400, Fig. 3 is the actual margin sum of errors phase error of each passage of array antenna, utilizes method of the present invention to estimate, just can obtain amplitude gain sum of errors phase error shown in Figure 4.
Data shown in comparison diagram 3 and Fig. 4 adopt method of the present invention can effectively estimate the channel error of array antenna as can be seen.
Claims (1)
1, a kind of blind estimating method of array antenna channel error is characterized in that adopting following step:
Step 1 (signalization source):
And array antenna between distance satisfy under the situation of far field condition, place a narrow band signal source s (t), the direction of arrival of incoming signal (signal source emission signal) is θ
1The start signal source, the data of gathering the individual snap of N (N=>400) at the output of array antenna, shutdown signal source then, the k that collects fast beat of data is:
x
1k=[x
1k(1)…x
1k(M)]
H (1)
X in the formula
1k(i) in 400 snaps that obtain before the expression displaced signal sources in k snap, i the data that array element is exported, ()
HExpression Hermitian transposition;
Step 2 (displaced signal sources):
Displaced signal sources s (t) is placed on the position that another one satisfies far field condition, and this moment, the azimuth of signal source was θ
2(θ
1≠ θ
2), start signal source once more is in the data that the output of array antenna is gathered the individual snap of N (N>400), shutdown signal source; The k that collects fast beat of data is:
x
2k=[x
2k(1)…x
2k(M)]
H (2)
X in the formula
2k(i) in N the snap that obtains after the expression displaced signal sources in k snap, i the data that array element is exported;
Step 3 (data are synthetic):
Data according to twice collection before and after moving obtain two covariance matrix R
1And R
2,
Step 4 (asking intermediate variable):
To R
1And R
2Make feature decomposition respectively, can obtain the characteristic vector u of two big characteristic values of correspondence
M1And u
M2, these two characteristic vectors are obtained two intermediate variables with its first element do normalization,
Step 5 (estimation azimuth):
Adopt MUSIC algorithm estimated signal source to move the azimuth of front and back, just can obtain
With
Step 6 (DOA estimation domain):
Define two independent variable territories, with
With
Being the center, is that radius is expanded to both sides with δ (desirable 0.3 degree of reference value), obtains two DOA estimated value intervals
With
Step 7 (the estimation flow process of array antenna channel error)
At first, the phase error of k the passage of data rough estimate array antenna that obtains before moving according to signal source:
In the formula (7), the angle of ∠ () operator representation calculated complex, d is the spacing of adjacent array element, λ is the wavelength of incoming signal; The phase error of k the passage of data rough estimate array antenna that obtains after moving according to signal source:
Then, introduce a parameter:
At last, determine the true incidence angle of signal source:
A with
With
Be starting point,
Calculate according to formula (7) and (8)
With
And then according to formula (9) calculating ζ (0); Step-length ε (its reference value can get 0.01) makes l=0;
ψ
1(l)=ψ
1(0)+lε
B l=l+1, ψ
2(l)=ψ
2(0)+and l ε, calculate ξ (l);
If c
Then change b;
D seeks the minimum of ξ, if minimum does not exist, then suitably increases δ (reference value is got δ=δ+0.5), returns a;
The angle value ψ of the minimum correspondence of e ξ
1(l) and ψ
2(l) be exactly real azimuth angle theta before and after signal source moves
1, θ
2
Step 8 (range error of computing array antenna channels and phase error)
The u that obtains according to step 4
S1, the range error of k passage of computing array antenna is estimated
For:
In the formula (10), u
S1(k) expression vector u
S1In k element
The u that obtains according to step 4 again
S1The θ that obtains with step 7
1, k channel phases estimation error of computing array antenna
For:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2003101110763A CN100399730C (en) | 2003-12-03 | 2003-12-03 | Blind estimating method for error of array antenna channel |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2003101110763A CN100399730C (en) | 2003-12-03 | 2003-12-03 | Blind estimating method for error of array antenna channel |
Publications (2)
Publication Number | Publication Date |
---|---|
CN1625088A true CN1625088A (en) | 2005-06-08 |
CN100399730C CN100399730C (en) | 2008-07-02 |
Family
ID=34759317
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNB2003101110763A Expired - Fee Related CN100399730C (en) | 2003-12-03 | 2003-12-03 | Blind estimating method for error of array antenna channel |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN100399730C (en) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101938305A (en) * | 2010-08-13 | 2011-01-05 | 四川九洲电器集团有限责任公司 | Amplitude and phase calibration method of phased array system receiving channel |
CN102544755A (en) * | 2011-12-31 | 2012-07-04 | 哈尔滨工业大学 | Uniform linear array calibration method based on strong scattering points |
CN109828249A (en) * | 2019-02-22 | 2019-05-31 | 北京遥感设备研究所 | A kind of calibrating method for synthetic aperture radiometer array |
CN110095658A (en) * | 2018-01-31 | 2019-08-06 | 罗克韦尔柯林斯公司 | Method and system for ESA measurement |
CN111060871A (en) * | 2019-12-26 | 2020-04-24 | 南京长峰航天电子科技有限公司 | Five-element array positioning method and device based on improved genetic algorithm |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1120858B1 (en) * | 1999-12-15 | 2007-04-04 | Nippon Telegraph and Telephone Corporation | Adaptive array transceiver apparatus |
US6765529B2 (en) * | 2000-07-14 | 2004-07-20 | Sanyo Electric Co., Ltd. | Calibration device, adaptive array device, calibration method, program recording medium and program |
JP3651430B2 (en) * | 2001-09-17 | 2005-05-25 | 日本電気株式会社 | Array antenna calibration apparatus and calibration method |
-
2003
- 2003-12-03 CN CNB2003101110763A patent/CN100399730C/en not_active Expired - Fee Related
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101938305A (en) * | 2010-08-13 | 2011-01-05 | 四川九洲电器集团有限责任公司 | Amplitude and phase calibration method of phased array system receiving channel |
CN101938305B (en) * | 2010-08-13 | 2012-12-26 | 四川九洲电器集团有限责任公司 | Amplitude and phase calibration method of phased array system receiving channel |
CN102544755A (en) * | 2011-12-31 | 2012-07-04 | 哈尔滨工业大学 | Uniform linear array calibration method based on strong scattering points |
CN102544755B (en) * | 2011-12-31 | 2013-12-11 | 哈尔滨工业大学 | Uniform linear array calibration method based on strong scattering points |
CN110095658A (en) * | 2018-01-31 | 2019-08-06 | 罗克韦尔柯林斯公司 | Method and system for ESA measurement |
CN109828249A (en) * | 2019-02-22 | 2019-05-31 | 北京遥感设备研究所 | A kind of calibrating method for synthetic aperture radiometer array |
CN111060871A (en) * | 2019-12-26 | 2020-04-24 | 南京长峰航天电子科技有限公司 | Five-element array positioning method and device based on improved genetic algorithm |
Also Published As
Publication number | Publication date |
---|---|
CN100399730C (en) | 2008-07-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107329108B (en) | The relatively prime array Wave arrival direction estimating method rebuild based on interpolation virtual array covariance matrix Toeplitzization | |
CN1262130C (en) | Method of estimating space reaching direction of fixed wave beam | |
CN107505602A (en) | DOA estimation method based on DFT under nested battle array | |
CN1131645C (en) | Method and system for creating RF signature for signature region in wireless communication system | |
CN107436421A (en) | Mixed signal DOA estimation method under a kind of management loading framework | |
CN1131644C (en) | Method and system for comparing measured radio frequency signal propagation characteristics in wireless communication system | |
CN103546221B (en) | A kind of wide-band coherent signal direction of arrival estimation method | |
CN106324559A (en) | Large-baseline four-element array broadband signal direction finding system and method | |
CN112051540B (en) | Quick high-precision direction finding method | |
CN109471063B (en) | Uniform linear array high-resolution direction-of-arrival estimation method based on delayed snapshot | |
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) | |
CN106443569A (en) | Robust adaptive beamforming method based on steering vector correction | |
CN113835063B (en) | Unmanned aerial vehicle array amplitude and phase error and signal DOA joint estimation method | |
CN105572629B (en) | A kind of two-dimentional direction-finding method of low computational complexity suitable for General Cell structure | |
CN1625088A (en) | Blind estimating method for error of array antenna channel | |
CN107144811A (en) | A kind of cyclic subspace direction-finding method of single channel receiving array signal | |
CN108614234B (en) | Direction-of-arrival estimation method based on multi-sampling snapshot co-prime array received signal fast Fourier inverse transformation | |
CN112731280B (en) | ESPRIT-DOA estimation method in inter-mass array mixed noise environment | |
CN113391261A (en) | Extended co-prime array rapid DOA estimation method based on Taylor expansion | |
CN104656055A (en) | Single signal arrival angle estimating method based on large-scale multi-antenna system | |
CN106507952B (en) | A kind of quick spatial spectrum computational methods based on circle battle array | |
CN1580813A (en) | Method for conducting array corrction and information source direction finding by using high precision auxiliary array element | |
CN1199052C (en) | Method of raising accuracy estimation of received signal wave arrival direction | |
CN1297075C (en) | Interference blocking spatial spectrum estimation method | |
CN109239648A (en) | Spectrum correlation subspaces direction-finding method based on symmetrical cycle frequency |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20080702 Termination date: 20101203 |