CN104297718A - Interferometer array integrated correction method - Google Patents
Interferometer array integrated correction method Download PDFInfo
- Publication number
- CN104297718A CN104297718A CN201410513334.9A CN201410513334A CN104297718A CN 104297718 A CN104297718 A CN 104297718A CN 201410513334 A CN201410513334 A CN 201410513334A CN 104297718 A CN104297718 A CN 104297718A
- Authority
- CN
- China
- Prior art keywords
- interferometer
- pitching
- baseline
- angle
- error
- 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/023—Monitoring or calibrating
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Instruments For Measurement Of Length By Optical Means (AREA)
Abstract
The invention provides an interferometer array integrated correction method. The method includes the steps of firstly, parameterizing a target signal initial incidence direction, a processing and installing comprehensive error and a feeder line passageway comprehensive error; secondly, obtaining a measuring equation set according to a known track rotation antenna array, and obtaining error parameters and incidence signal directions through parameter estimation; thirdly, modifying the error parameters into an interference angle measurement formula to obtain a precisely-corrected angle through the iteration technology. By means of the method, the incidence signal directions do not need to be accurately calibrated, the basic known quantities of the method are the rotation angle and the phase difference measurement value, the interferometer array processing and installing comprehensive error and the feeder line passageway comprehensive error can be obtained at the same time through parameter estimation, the problem that the direction error and the pitching error are coupled can be solved through the iteration technology when the angle measurement error correction is executed, and finally the precisely-corrected angle is obtained. The interferometer array integrated correction method is accurate in error model, simple in execution process, high in correction precision and capable of improving the correction efficiency and the angle measurement accuracy of an interferometer.
Description
Technical field
The present invention relates to a kind of interferometer array integrated correction method, be applied in the direction finding of radiation source with in locating, belong to radio direction finding (RDF) field.This technology is widely applied in fields such as radar, sonar, mobile communication, geophysical surveys, and particularly in electronic countermeasure, interferometer can timely and accurately to target localization.
Background technology
The ultimate principle of interferometer direction finding is the incident angle being calculated information source by the difference of the different array element of measuring radiation source arrival interferometer.The current raising along with radiation source working frequency range, in order to ensure the precision of interferometry, to the processing of the phase equalization of each array element of interferometer, array and installation accuracy, strict requirement must be proposed, and the passage consistency of array element rear end is strictly controlled, to reduce each interchannel Sensor gain and phase perturbations to greatest extent, but be limited to the impact of processing installation accuracy, calibration measure must be taked to calibrate processing composition error, feeder line passage composition error are installed.
By Wanfang Database, domestic open academic journal, academic dissertation, meeting paper, patent, standard, regulation etc. are retrieved, search condition is: cql: // interfere and (angle measurement or direction finding) and (calibration or corrects), retrieve the document that 45 sections of conditions are relevant altogether; By to defense S & T information Resource service system (DTIRSS, Defence Technology Information Resources Service System) retrieve, this service system comprises IEL, Elsevier, EI, Springer, AD database, NASA database, DE database etc., search condition is interferometer+ quadratic search (Direction Finding or angle measuring)+quadratic search calibration, retrieves 7 sections of condition pertinent literatures altogether.
Recognize from domestic and international open source literature retrieval, interferometer array calibration steps is broadly divided into following a few class:
1) only consider passage sensor gain and phase uncertainties, or regard center of antenna position deviation and the error synthesis of passage amplitude phase unbalance as passage sensor gain and phase uncertainties and carry out measuring or estimating, form calibration parameter or correction card.Concrete, the method comprises and arranges reference channel and compare method, Measurement and calibration table method, the width phase parameter estimation technique, and wherein the width phase parameter estimation technique comprises again correlation method, FFT method, least square method etc.Although above-mentioned concrete grammar is different, but common ground only considers channel amplitude and phase information, when antenna misalignment ideal position and after can not ignore, inter-channel phase error can change along with the change of signal incident angle, now, this method only considering calibration of amplitude and phase is just no longer accurate; These class methods can not assess the processing installation composition error of bay in addition.Typical case's document is as follows:
[1] Cui Yuxin. multi-channel digital interferometer direction finding receiver phase Design of Calibration System [J]. marine electronic resists, and 2010,33 (1)
[2] Ding Yikai. interferometer channel error corrects and signal acquisition circuit design [D]. Xian Electronics Science and Technology University, 2010.
[3] Zhang Haiyan. the research of Five-channel phase-interfer-ometer direction-finding and realization [D] (correlation method calibration, the calibration of FFT method, the calibration of LS method). Chengdu University of Technology, 2004.
[4] Chen Junhua. circle battle array interferometer direction finding technology [D]. Shanghai Communications University, 2008.
[5] season weighs. the realization [J] of reconnaissance system for radar near filed phase bearing calibration. and aerospace electronic warfare, 2013,29 (5)
[6] Wen Fuzhong, Li Bo. the phase equalization bearing calibration [J] of Multi-channel interferometer direction finding. electronic information countermeasure techniques, 2012,27 (2)
[7] hair tiger, Yang Jianbo, Liu Peng. the interferometer direction finding state of the art and developmental research [J]. electronic information countermeasure techniques, 2010,25 (6)
[8] Guo Fang. many channel receivers broadband alignment technique research and implementation [D]. University of Electronic Science and Technology, 1999.
[9]Schneider,Stephen?W.The?Calibration?of?Four-Arm?Spiral?Modal?Measurements?for?Angle-of-Arrival?Determination?(Preprint).AD?Report,2006?JUL.
[10]HAN?Yue-tao.Multi-baseline?interferometer?broadband?real-time?phase?calibration?method.2011?4th?International?Congress?on?Image?and?Signal?Processing.
2) center of antenna and passage amplitude phase unbalance error are demarcated respectively.First to antenna phase center independent measurement, obtain the position of its phase center, then postpone to calibrate to channel phases.These class methods depend on the Accurate Calibration to bay phase center, and after antenna array completes, the relative position of each antenna also depends on its installation site, if group is after battle array completes, then demarcate, are often limited by antenna feeder and passage, are unfavorable for implementing.Typical case's document is as follows:
[1] Feng Xiaoyu. the Research on Calibration Technology [D] of interferometer direction finding. Xian Electronics Science and Technology University, 2012.
[2] Wang Yuwen. wide-band ratio width adds interferometer direction finding antenna array [D]. University of Electronic Science and Technology, 1999.
3) signal source that incident direction is known is set, calculates lubber line error, realize calibrating interferometer.The background of pertinent literature revises the baseline vector error in the application of Long baselines (km magnitude) interferometer, and method is by arranging fixing derived reference signal, and this method needs the incident direction by miscellaneous equipment or means precision calibration signal.Typical case's document is as follows:
[1] Zhong Danxing, Yang Zhengbin, Monday the .LBI direction detecting positioning system such as space many calibrations source correcting algorithm [J]. systems engineering and electronic technology, 2008,30 (5)
[2] Li Yuxing, Li Yanjie, Wang Pei. interferometer system error analysis and correction [A]. the tenth academic nd Annual Meeting collection [C] .2008. of national radar
4) a kind of special linear array interferometer is described in the document of the people such as Kawase, this interferometer is made up of three pairs of antennas, first pair of astronomical cycle is on rotating mechanism, for echo signal being reflexed to middle second pair of antenna, signal reflex is entered last on receiving antenna by second pair of antenna again, rotate first pair of antenna with certain stepping, rotate a circle the phase difference measurements system of equations that just can obtain containing error, can obtain interferometer system error by solving it.When document method implements calibration, first need the position angle and the angle of pitch that are obtained incident target by optical telescope; Secondly the difference error that document is considered comprises the droop error of the error in length of first pair of antenna pivot arm, anglec of rotation error and the surfaces of revolution, due to three pairs of antennas successively will be experienced, incoming signal transmission path model is different from general interferometer, and therefore its error model is not suitable for general interferometer; Do not provide further angle measurement correction algorithm after Access to publication error term, do not point out how to eliminate orientation and pitch error coupled problem in interference angle measuring algorithm.Pertinent literature is as follows:
[1]Kawase,Seiichiro.Radio-interferometer?for?monitoring?orbital-resource?congestion.ICSSC.Japan,April?15-19,2003.
5) based on the array calibration technology of Estimation of Spatial Spectrum, obtain in " Estimation of Spatial Spectrum theoretical and algorithm " book that this technology is shown at Wang Yongliang etc. and comparatively comprehensively sum up, these calibration algorithms include source calibration, self calibration and the collimation technique based on auxiliary array element, these algorithms can by Ro-vibrational population, phase position error and array element mutual coupling error etc. are brought in Estimation of Spatial Spectrum model, subspace and signal space characteristic etc. is utilized to carry out estimation calibration, these collimation technique basic theories are Estimation of Spatial Spectrum technology, the data used are sampling snaps of fixing information source incident direction, in addition for non-homogeneous array element of structuring the formation, Power estimation technology could be used after generally needing pre-service.Typical case's document is as follows:
[1] work such as Wang Yongliang, Chen Hui. Estimation of Spatial Spectrum theory and algorithm. publishing house of Tsing-Hua University, 2004.
[2] Chen Juan. based on modeling and interferometer direction finding algorithm research (passive calibration) [D] of the satellite reconaissance link of ADS. University of Electronic Science and Technology, 2009.
[3] Wang Lei. based on the radio communication reconnaissance system key technology research [D] of Array Signal Processing. China Science & Technology University, 2005.
[4] first duckweed. based on the passive radar direction finding technology [D] of broadband system. Harbin Engineering University, 2011.
[5] burnt strong. satellite communication system interference source DOA estimation method research [D]. Harbin Institute of Technology, 2008.
[6] Zhang Weihua. communication direction finding modeling and simulation research [D]. the National University of Defense Technology, 2005.
[7] Zhang Yaguang. array received signal imitation and DOA algorithm for estimating are studied [D]. University of Electronic Science and Technology, 2009.
Summary of the invention
The technical matters that the present invention solves is: install composition error and feeder line passage composition error to obtain interferometer array processing simultaneously, seek a kind of without Accurate Calibration incoming signal direction, without the interferometer array calibration steps of Estimation of Spatial Spectrum technology; Meanwhile, also need to seek a kind of angle calibration system method that effectively can utilize interferometer array calibration result, the method needs to overcome azimuthal error and pitch error when angle measurement is resolved and to intercouple problem, eliminates interferometer array processing and installs composition error and feeder line passage composition error to the impact of angle measurement accuracy.The invention provides a kind of interferometer array integrated correction method, can solve the technical problem, the present invention does not need the incident direction of accurate known signal, do not need to adopt Estimation of Spatial Spectrum technology yet, the present invention is basic known quantity with the anglec of rotation, first by signal original incident direction, composition error and feeder line passage composition error parametrization are simultaneously installed in processing, secondly rotary antenna battle array obtains and measures system of equations, then error parameter and incoming signal direction is obtained by parameter estimation, finally error parameter correction is entered to interfere angle measurement formula, pass through iterative technique, overcome azimuthal error and pitch error to intercouple problem, thus obtain the angle after close adjustment.
The technical scheme that the present invention solves is: a kind of interferometer array integrated correction method, comprises the error model stage of setting up, obtains the measurement equation stage, solve the error parameter stage and obtain the angle measurement calibration result stage;
It is described that to set up error model stage etch as follows:
(1) on interferometer front, set up surving coordinate system, to choose on interferometer front an array element as interferometer measurement coordinate origin O, using desirable azimuth reference datum as Y-axis, using desirable pitching baseline as Z axis, by right-hand rule determination X-axis; The angle in objective definition incident direction and XOY face is pitching angle theta, and the angle in the projection of objective definition incident direction in XOY plane and XOZ face is position angle φ;
(2) in the interferometer measurement coordinate system set up in described step (1), if target original incident position angle is φ
0, the original incident angle of pitch is θ
0; If total M array element in interferometer array, 0th array element is positioned at initial point O place, i-th array element is arranged in desirable orientation baseline, and a jth array element is arranged in desirable pitching baseline, i, j=1,2 ..., M-1 and i ≠ j, there is the processing installation composition error in X-axis, Y-axis, Z axis three directions relative to ideal position in the i-th array element, is respectively Δ x
i, Δ y
i, Δ z
i, jth array element relative to ideal position exist X-axis, Y-axis, Z axis three directions processing install composition error be respectively Δ x
j, Δ y
j, Δ z
j; In addition, the Received signal strength of each array element enters interferometer signal processor by respective feeder line passage, strictly not isometric between feeder line passage, if with the feeder line passage of the 0th array element for reference length, then by i-th, the feeder line passage of j array element is set to Δ d respectively relative to the feeder line passage composition error of the 0th array element
i, Δ d
j; Respectively by Δ x
i, Δ y
i, Δ z
i, Δ d
i, θ
0, φ
0, Δ x
j, Δ y
j, Δ z
j, Δ d
j, θ
0, φ
0as six unknown parameters of azimuth reference datum, pitching baseline, set up the phase eikonal equation of interferometer azimuth reference datum, pitching baseline:
In formula, i, j=1,2 ..., M-1 and i ≠ j, L
ibe the i-th array element azimuth reference datum design length, L
jfor jth array element pitching baseline design length,
for the interferometer phase aberration measurements of azimuth reference datum,
for the interferometer phase aberration measurements of pitching baseline, λ is the wavelength of incoming signal;
(3) interferometer be arranged on elaborate servo mechanism or be erected on precise rotating platform, by adjusting the position of interferometer, the Z axis of the interferometer measurement coordinate system that described step (1) is set up overlaps with the direction rotating shaft of elaborate servo mechanism or precise rotating platform, the Y-axis of interferometer measurement coordinate system that described step (1) is set up and the pitching shaft parallel of elaborate servo mechanism or precise rotating platform; Keep target incident direction constant, make interferometer respectively around surving coordinate system Z axis and Y-axis anglec of rotation δ φ, δ θ, after coordinate conversion, by the interferometer azimuth reference datum phase eikonal equation in described step (2) and pitching baseline phase eikonal equation, by rotation of coordinate, obtaining take rotation amount δ φ, δ θ as azimuth reference datum, the pitching baseline six parameter phase difference model of basic known quantity:
In formula, i, j=1,2 ..., M-1 and i ≠ j;
It is as follows that equation stage etch is measured in described acquisition:
(4) control servo control mechanism or turntable rotate δ φ successively
k, δ θ
kaccurate step angle, k represents that kth time rotates, and 0≤k≤N-1, N represents number of revolutions, N>=0, resident certain hour on each step angle, gathers the interferometer phase aberration measurements that described current accurate step angle is corresponding;
(5) the interferometer phase aberration measurements gathered in every section of residence time in described step (4) is averaged, obtain corresponding step angle upper position baseline phase aberration measurements
pitching baseline phase aberration measurements
(6) described step (5) is processed the phase aberration measurements obtained
substitute into azimuth reference datum, pitching baseline six parameter phase difference model in described step (3), obtain azimuth reference datum, the non-linear six parameter phase difference measurements system of equations of pitching baseline:
In formula, i, j=1,2 ..., M-1 and i ≠ j, x
i=[Δ x
i, Δ y
i, Δ z
i, Δ d
i, θ
0, φ
0]
trepresent the parameter that the i-th array element is to be asked, x
j=[Δ x
j, Δ y
j, Δ z
j, Δ d
j, θ
0, φ
0]
trepresent the parameter that jth array element is to be asked, w
ik, w
jkrepresent the measurement noises that difference collection is introduced;
represent respectively and walk postrotational azimuth reference datum, pitching baseline six parameter phase eikonal equation through k, form is as follows:
In formula, 0≤k≤N-1, l represents the l time rotation, 0≤l≤k, δ φ
lbe the l time orientation to accurate step angle, δ θ
lbe that the l time pitching is to accurate step angle;
It is described that to solve error parameter stage etch as follows:
(7) non-linear six parameter phase difference measurements system of equations step (6) obtained, at initial value x
i0, x
j0place's linearization, to x in the i-th array element
iinitialize x
i0, x
i0=[Δ x
i0, Δ y
i0, Δ z
i0, Δ d
i0, θ
00, φ
00]
t, wherein Δ x
i0, Δ y
i0, Δ z
i0the initial value of composition error is installed in processing for X-axis, Y-axis, Z axis three directions, Δ d
i0for the initial value of feeder line passage composition error, θ
00for the initial value of pitching initial incidence angle, φ
00for the initial value of orientation initial incidence angle, at initial value x
i0place carries out first order Taylor expansion to the non-linear six parameter phase difference measurements system of equations of azimuth reference datum in described step (6), obtains azimuth reference datum lienarized equation group, is write as vector form:
Α in formula
ifor first-order partial derivative matrix, form is as follows:
for difference residual error, form is as follows:
Dx
ifor initial value x
i0correction; w
i=[w
i0, w
i1... w
iN-1]
tfor measurement noises;
In like manner, to x in jth array element
jinitialize x
j0, x
j0=[Δ x
j0, Δ y
j0, Δ z
j0, Δ d
j0, θ
00, φ
00]
t, wherein Δ x
j0, Δ y
j0, Δ z
j0the initial value of composition error is installed in processing for X-axis, Y-axis, Z axis three directions, Δ d
j0for the initial value of feeder line passage composition error, θ
00for the initial value of pitching initial incidence angle, φ
00for the initial value of orientation initial incidence angle, at x
j0place carries out first order Taylor expansion to the non-linear six parameter phase difference measurements system of equations of pitching baseline in described step (6), obtains pitching baseline lienarized equation group, is write as vector form:
In formula, Α
jfor first-order partial derivative matrix, form is as follows:
for difference residual error, form is as follows:
Dx
jfor initial value x
j0correction, i, j=1,2 ..., M-1 and i ≠ j; w
j=[w
j0, w
j1... w
jN-1]
tfor measurement noises;
(8) obtain the least square solution of azimuth reference datum, pitching baseline lienarized equation group in step (7) according to weighted least-squares method, this solution is initial value x
i0, x
j0correction dx
i, dx
j:
P in formula
i, P
jfor least square weighting matrix, i, j=1,2...M-1 and i ≠ j;
(9) according to the linearization six parameter phase difference measurements solution of equations that described step (8) obtains, judge to resolve end condition, if the absolute value of least square solution | dx
i|, | dx
j| be less than the interferometer array calibration accuracy of setting, then initial value x
i0, x
j0be non-linear six parameter phase difference measurements solution of equations; If the absolute value of least square solution | dx
i|, | dx
j| be more than or equal to the interferometer array calibration accuracy of setting, utilize least square solution to initial value x
i0, x
j0revise, obtain correction x '
i0, x '
j0, utilize revised x '
i0, x '
j0x in replacement step (7)
i0, x
j0, re-execute described step (7), (8), (9), until dx
i, dx
jmeet the requirement of interferometer array calibration accuracy, initial value x
i0, x
j0correction formula as follows:
x′
i0=x
i0+dx
i
x′
j0=x
j0+dx
j
I, j=1 in formula, 2 ..., M-1 and i ≠ j;
This step finally obtains composition error item corresponding to i-th array element
X
i0=[Δ x
i0, Δ y
i0, Δ z
i0, Δ d
i0, θ
00, φ
00]
t, the composition error item that a jth array element is corresponding
X
j0=[Δ x
j0, Δ y
j0, Δ z
j0, Δ d
j0, θ
00, φ
00]
t, namely the processing of the i-th array element in X-axis, Y-axis, Z axis three directions is installed composition error and is respectively Δ x
i0, Δ y
i0, Δ z
i0, feeder line passage composition error is Δ d
i, the processing of jth array element in X-axis, Y-axis, Z axis three directions is installed composition error and is respectively Δ x
j0, Δ y
j0, Δ z
j0, feeder line passage composition error is Δ d
j;
Described acquisition angle measurement calibration result stage etch is as follows:
(10) the whole composition error item Δ x will obtained in step (9)
i0, Δ y
i0, Δ z
i0, Δ d
i0, Δ x
j0, Δ y
j0, Δ z
j0, Δ d
j0, substitute into the phase difference measurements equation of the orientation of described step (3), pitching baseline, introduce noise item w during azimuth reference datum interferometry simultaneously
i, noise item w when pitching base line interference is measured
j, obtaining being with noisy take rotation amount δ φ, δ θ as azimuth reference datum, the pitching baseline six parameter phase eikonal equation of basic known quantity:
In formula,
for band noisy orientation difference, pitching difference that interferometry obtains;
(11) by carrying out computing to base line interference phase level difference measurements multiple in step (10), synthesize a virtual baseline and interfere phase level difference measurements, virtual baseline after synthesis is generally not more than half wavelength lambda/2 of incoming signal, utilize the phase level difference measurements after synthesis, realize resolving without fuzzy angle measurement, for azimuth reference datum phase difference measurements equation in step (10), if i=m and i=n, m<n, 1≤m≤M-1, during 1≤n≤M-1, synthesize a short virtual baseline by computing, computing operator is H
i(), then corresponding orientation virtual baseline phase difference measurements equation is written as:
Δ x
n0, Δ y
n0, Δ z
n0, Δ d
n0, L
n, w
nrespectively during corresponding i=n
Δ x
i0, Δ y
i0, Δ z
i0, Δ d
i0, L
i, w
i;
Δ x
m0, Δ y
m0, Δ z
m0, Δ d
m0, L
m, w
mrespectively during corresponding i=m
Δ x
i0, Δ y
i0, Δ z
i0, Δ d
i0, L
i, w
i;
In like manner, for pitching baseline phase difference measurements equation in described step (10), if j=p and j=q, p<q, synthesize a short virtual baseline when 1≤p≤M-1,1≤q≤M-1, computing operator is H
j(), then write out pitching virtual baseline phase difference measurements equation:
Δ x
q0, Δ y
q0, Δ z
q0, Δ d
q0, L
q, w
qrespectively during corresponding j=q
Δ x
j0, Δ y
j0, Δ z
j0, Δ d
j0, L
j, w
j;
Δ x
p0, Δ y
p0, Δ z
p0, Δ d
p0, L
p, w
prespectively during corresponding j=p
Δ x
j0, Δ y
j0, Δ z
j0, Δ d
j0, L
j, w
j;
(12) in the orientation that step (11) obtains, pitching virtual baseline phase difference measurements equation, differ from orientation
pitching differs
middlely deduct 2 π H respectively
i(Δ d
n0, Δ d
m0)/λ, 2 π H
j(Δ d
q0, Δ d
p0)/λ, and will subtract each other result reduction to (-π, π] principal value interval, obtain following without integral circumference ambiguity orientation, pitching baseline phase difference measurements system of equations:
In formula
<> symbol represent that reduction is arrived (-π, π] principal value interval; w
1=H
i(w
n, w
m), w
2=H
j(w
q, w
p); F in formula
1(θ, φ), f
2(θ, φ) is defined as follows:
(13) vector y=[θ, φ] is defined
t, to y initialize y
0,
herein
for angle of pitch initial value,
for position angle initial value, at y
0near to carrying out first order Taylor expansion without integral circumference ambiguity orientation, pitching baseline phase difference measurements system of equations in step (12), obtain linearization without fuzzy phase difference measurements system of equations, write as vector form:
In formula, u is that form is as follows without fuzzy phase difference measurements residual error:
B is that form is as follows without fuzzy phase difference measurements system of equations first-order partial derivative matrix:
Dy=[d θ, d φ]
trepresent y
0correction variable; W=[w
1, w
2]
tfor measurement noises;
(14) what obtain described step (13) directly inverts without fuzzy phase difference measurements system of equations, can obtain initial value y in step (13)
0correction variable dy:
dy=B
-1u
(15) linearization obtained according to described step (14) is without fuzzy phase difference measurements solution of equations, judge to resolve end condition, if linearization is without the absolute value of fuzzy phase difference measurements solution of equations | dy| is less than the Interferometer angle measurement precision of permission, then initial value y
0be non-linear without fuzzy phase difference measurements solution of equations; If linearization is without the absolute value of fuzzy phase difference measurements solution of equations | dy| is more than or equal to the Interferometer angle measurement precision of setting, utilizes dy to initial value y
0revise, obtain modified value y '
0, with y '
0replace the initial value y in described step (13)
0, re-execute described step (13), (14), (15), until dy meets the Interferometer angle measurement precision of setting, to initial value y
0correction formula as follows:
y′
0=y
0+dy
This step finally obtains the Interferometer angle measurement value after calibrating
namely the angle of pitch after Interferometer angle measurement calibration is
position angle after calibration is
The present invention's advantage is compared with prior art:
(1) the present invention is in array error item, using array element phase heart deviation and installation deviation as a composition error item, using array element feeder line error and channel error as a composition error item, carry out parametrization simultaneously, by measuring and parameter estimation in interferometer phase, calculate composition error item, for the interferometer that frequency range is higher, as millimeter wave interferometer, its wavelength is shorter, composition error item on angle measurement accuracy impact significantly, is ignored approximate calibration method that is some or certain several error component if adopt again, can not be met the requirement of angle measurement accuracy;
(2) the present invention only needs to be arranged on turntable or servo control mechanism when starting by interferometer front, then by by certain track angle measuring interferometer front, obtains phase difference measurements equation, finally carries out mathematics manipulation and obtain composition error; This process does not need to carry out independent measurement to target incident direction, phase heart position, installation site, feeder line time delay and channel time delay etc., and manual intervention rate is low, is easy to be designed to interferometer array auto-calibration device, improves the efficiency of calibration;
(3) the present invention avoids when array calibration using nonuniform noise Estimation of Spatial Spectrum technology, and theory and the computation complexity of the array calibration technology of Estimation of Spatial Spectrum are all higher;
(4) the present invention is in order to improve angle measurement accuracy, do not need to form look-up table, but processing is installed in composition error, feeder line passage composition error substitution interference angle measurement equation, utilize alternative manner, overcome azimuthal error and pitch error to intercouple problem, and then eliminate composition error item to the impact of angle measurement accuracy.
Accompanying drawing explanation
Fig. 1 is interferometer array integrated correction method flow diagram of the present invention;
Fig. 2 is that composition error schematic diagram is installed in a kind of L-type interferometer array processing of the present invention;
Fig. 3 is a kind of L-type interferometer array of the present invention feeder line passage composition error schematic diagram;
Fig. 4 (a) is front angle Error Graph for the present invention calibrates;
Fig. 4 (b) is rear angle Error Graph for the present invention calibrates.
Embodiment
Be example below in conjunction with 5 array elements " L " the type interferometer shown in accompanying drawing 2, Fig. 3, the present invention is described in further detail, a kind of interferometer array integrated correction method as shown in Figure 1 comprises the error model stage of setting up, obtains the measurement equation stage, solves the error parameter stage and obtain the angle measurement calibration result stage, and concrete embodiment is described below:
Set up error model stage etch as follows:
(1) on interferometer front, surving coordinate system is set up.In " L " type interferometer shown in Fig. 2, to choose on interferometer front array element 0 as interferometer measurement coordinate origin O, using desirable azimuth reference datum O2 as Y-axis, using desirable pitching baseline O4 as Z axis, by right-hand rule determination X-axis; The angle in objective definition incident direction and XOY face is pitching angle theta, and the angle in the projection of objective definition incident direction in XOY plane and XOZ face is position angle φ;
(2) set up containing all orientation of composition error item, the phase eikonal equation of pitching baseline.In the interferometer measurement coordinate system that step (1) is set up, if target original incident position angle is φ
0, the original incident angle of pitch is θ
0; In " L " type interferometer array shown in Fig. 2, total M=5 array element, the 0th array element is positioned at initial point O place, and i-th array element is arranged in the azimuth reference datum of Theoretical Design, and a jth array element is arranged in the pitching baseline of Theoretical Design, i=1,2, j=3,4.The i-th, there is the processing installation composition error in X-axis, Y-axis, Z axis three directions in j array element relative to Theoretical Design position, be respectively Δ x
i, Δ y
i, Δ z
i, Δ x
j, Δ y
j, Δ z
j; In addition, as shown in Figure 3, the Received signal strength of each array element enters interferometer signal processor by respective feeder line passage, strictly not isometric between feeder line passage, if with the feeder line passage of the 0th array element for reference length, then can by i-th, the feeder line passage of j array element is designated as Δ d respectively relative to the feeder line passage composition error of the 0th array element
i, Δ d
j, then can set up the phase eikonal equation of interferometer azimuth reference datum, pitching baseline:
In formula, i=1,2, j=3,4, L
ibe the i-th array element azimuth reference datum design length, L
jfor jth array element pitching baseline design length,
for the interferometer phase aberration measurements of azimuth reference datum,
for the interferometer phase aberration measurements of pitching baseline, λ is the wavelength of incoming signal;
(3) foundation take rotation amount as six parameter phase difference models of basic known quantity.Certain interferometer shown in Fig. 2 is arranged in elaborate servo mechanism, by adjusting the position of interferometer, the Z axis of the interferometer measurement coordinate system that step (1) is set up overlaps with elaborate servo mechanism direction rotating shaft, the Y-axis of the interferometer measurement coordinate system that step (1) is set up and elaborate servo mechanism pitching shaft parallel; Keep target incident direction constant, make interferometer respectively around surving coordinate system Z axis and Y-axis anglec of rotation δ φ, δ θ, after coordinate conversion, according to the interferometer azimuth reference datum phase eikonal equation described in step (2) and pitching baseline phase eikonal equation, obtaining take rotation amount δ φ, δ θ as azimuth reference datum, the pitching baseline six parameter phase difference model of basic known quantity:
In formula, i=1,2 j=3,4.
Obtain measurement equation stage etch as follows:
(4) control servo control mechanism or turntable rotate δ φ successively
k, δ θ
kaccurate step angle, k represents that kth time rotates, and 0≤k≤N-1, N represents number of revolutions, N>=0, resident certain hour on each step angle, gathers the interferometer phase aberration measurements that described current accurate step angle is corresponding;
(5) to phase aberration measurements average value processing.The interferometer phase aberration measurements gathered in every section of residence time in step (4) is averaged, obtains corresponding step angle upper position baseline phase aberration measurements
pitching baseline phase aberration measurements
(6) non-linear six parameter phase difference measurements system of equations are obtained.Step (5) is processed the phase aberration measurements obtained
substitute into azimuth reference datum, pitching baseline six parameter phase difference model in step (3), obtain azimuth reference datum, the non-linear six parameter phase difference measurements system of equations of pitching baseline:
In formula, i=1,2 j=3,4, x
i=[Δ x
i, Δ y
i, Δ z
i, Δ d
i, θ
0, φ
0]
trepresent the parameter that the i-th array element is to be asked, x
j=[Δ x
j, Δ y
j, Δ z
j, Δ d
j, θ
0, φ
0]
trepresent the parameter that jth array element is to be asked, w
ik, w
jkrepresent the measurement noises that difference collection is introduced;
represent respectively and walk postrotational azimuth reference datum, pitching baseline six parameter phase eikonal equation through k, form is as follows:
In formula, 0≤k≤N-1, l represents the l time rotation, 0≤l≤k, δ φ
lbe the l time orientation to accurate step angle, δ θ
lbe that the l time pitching is to accurate step angle;
It is described that to solve error parameter stage etch as follows:
(7) at initial value x
i0, x
j0place is by non-linear six parameter phase difference measurements system of equations linearizations.Non-linear six parameter phase difference measurements system of equations step (6) obtained, at initial value x
i0, x
j0place's linearization, to x in the i-th array element
iinitialize x
i0, x
i0=[Δ x
i0, Δ y
i0, Δ z
i0, Δ d
i0, θ
00, φ
00]
t, wherein Δ x
i0, Δ y
i0, Δ z
i0the initial value of composition error is installed in processing for X-axis, Y-axis, Z axis three directions, Δ d
i0for the initial value of feeder line passage composition error, θ
00for the initial value of pitching initial incidence angle, φ
00for the initial value of orientation initial incidence angle, at initial value x
i0place carries out first order Taylor expansion to the non-linear six parameter phase difference measurements system of equations of azimuth reference datum in described step (6), obtains azimuth reference datum lienarized equation group, is write as vector form:
Α in formula
ifor first-order partial derivative matrix, form is as follows:
for difference residual error, form is as follows:
Dx
ifor initial value x
i0correction; w
i=[w
i0, w
i1... w
iN-1]
tfor measurement noises.
In like manner, to x in jth array element
jinitialize x
j0, x
j0=[Δ x
j0, Δ y
j0, Δ z
j0, Δ d
j0, θ
00, φ
00]
t, wherein Δ x
j0, Δ y
j0, Δ z
j0the initial value of composition error is installed in processing for X-axis, Y-axis, Z axis three directions, Δ d
j0for the initial value of feeder line passage composition error, θ
00for the initial value of pitching initial incidence angle, φ
00for the initial value of orientation initial incidence angle, at x
j0place carries out first order Taylor expansion to the non-linear six parameter phase difference measurements system of equations of pitching baseline in described step (6), obtains pitching baseline lienarized equation group, is write as vector form:
In formula, Α
jfor first-order partial derivative matrix, form is as follows:
for difference residual error, form is as follows:
Dx
jfor initial value x
j0correction, w
j=[w
j0, w
j1... w
jN-1]
tfor measurement noises, i=1,2j=3 in formula, 4;
(8) least square solution of linearization six parameter phase difference measurements system of equations is asked.Initial value x in (7) is obtained according to weighted least-squares method
i0, x
j0correction dx
i, dx
j:
P in formula
i, P
jfor least square weighting matrix, i=1,2 j=3,4.
(9) if the absolute value of least square solution | dx
i|, | dx
j| be less than the interferometer array calibration accuracy of setting, then initial value x
i0, x
j0be non-linear six parameter phase difference measurements solution of equations; If the absolute value of least square solution | dx
i|, | dx
j| be more than or equal to the interferometer array calibration accuracy of setting, utilize least square solution to initial value x
i0, x
j0revise, obtain correction x '
i0, x '
j0, utilize revised x '
i0, x '
j0x in replacement step (7)
i0, x
j0, re-execute described step (7), (8), (9), until dx
i, dx
jmeet the requirement of interferometer array calibration accuracy.Initial value x
i0, x
j0correction formula as follows:
x′
i0=x
i0+dx
i
x′
j0=x
j0+dx
j
In formula, i=1,2 j=3,4;
This step finally obtains composition error item x corresponding to i-th array element
i0=[Δ x
i0, Δ y
i0, Δ z
i0, Δ d
i0, θ
00, φ
00]
t, the composition error item x that a jth array element is corresponding
j0=[Δ x
j0, Δ y
j0, Δ z
j0, Δ d
j0, θ
00, φ
00]
t, namely the processing of the i-th array element in X-axis, Y-axis, Z axis three directions is installed composition error and is respectively Δ x
i0, Δ y
i0, Δ z
i0, feeder line passage composition error is Δ d
i, the processing of jth array element in X-axis, Y-axis, Z axis three directions is installed composition error and is respectively Δ x
j0, Δ y
j0, Δ z
j0, feeder line passage composition error is Δ d
j;
Described calibration angle measurement result stage etch is as follows:
(10) the whole composition error item Δ x will obtained in step (9)
i0, Δ y
i0, Δ z
i0, Δ d
i0, Δ x
j0, Δ y
j0, Δ z
j0, Δ d
j0, substitute into the phase difference measurements equation of the orientation of described step (3), pitching baseline, introduce noise item w during azimuth reference datum interferometry simultaneously
i, noise item w when pitching base line interference is measured
j, obtaining being with noisy take rotation amount δ φ, δ θ as azimuth reference datum, the pitching baseline six parameter phase eikonal equation of basic known quantity:
In formula,
for orientation difference, pitching difference that interferometry obtains;
(11) orientation, pitching virtual baseline phase difference measurements system of equations is obtained.By interfering phase level difference measurements to carry out arithmetical operation to multiple Long baselines, synthesizing a Short baseline and interfering phase level difference measurements, utilizing this result, can realize resolving without fuzzy angle measurement.
As shown in Figure 2, for azimuth reference datum, if by arithmetical operation H when i=m and i=n, m=1, n=2
i(L
n, L
m) a synthesis short virtual baseline L
n-2L
m, now can according to operator H
i() writes out orientation virtual baseline phase difference measurements equation:
In like manner, for pitching baseline, if by arithmetical operation H when j=p and j=q, p=3, q=4
j(L
q, L
p) a synthesis short virtual baseline L
q-2L
p, now can according to operator H
j() writes out pitching virtual baseline phase difference measurements equation:
(12) deduct the impact of feeder line channel error, and will differ reduction to (-π, π] principal value interval, obtain non-linear without fuzzy phase difference measurements system of equations.From the phase difference measurements equation of the orientation that step (11) obtains, pitching baseline, deduction feeder line passage composition error Δ d
n0, Δ d
m0, Δ d
q0, Δ d
p0impact, namely differ from orientation
pitching differs
middlely deduct 2 π (Δ d respectively
n0-2 Δ d
m0)/λ, 2 π (Δ d
q0-2 Δ d
p0)/λ, and by result of calculation reduction to (-π, π] principal value interval, obtain following without integral circumference ambiguity orientation, pitching baseline phase difference measurements system of equations:
In formula
<> symbol represent that reduction is arrived (-π, π] principal value interval; w
1=w
n-2w
m, w
2=w
q-2w
p; F in formula
1(θ, φ), f
2(θ, φ) is defined as follows:
(13) vector y=[θ, φ] is defined
t, to y initialize y
0,
herein
for angle of pitch initial value,
for position angle initial value, at y
0near to carrying out first order Taylor expansion without integral circumference ambiguity orientation, pitching baseline phase difference measurements system of equations in step (12), obtain linearization without fuzzy phase difference measurements system of equations, write as vector form:
In formula, u is that form is as follows without fuzzy phase difference measurements residual error:
B is that form is as follows without fuzzy phase difference measurements system of equations first-order partial derivative matrix:
Dy=[d θ, d φ]
trepresent y
0correction variable; W=[w
1, w
2]
tfor measurement noises;
(14) what obtain described step (13) directly inverts without fuzzy phase difference measurements system of equations, can obtain initial value y in step (13)
0correction variable dy:
dy=B
-1u
(15) linearization obtained according to described step (14) is without fuzzy phase difference measurements solution of equations, judge to resolve end condition, if linearization is without the absolute value of fuzzy phase difference measurements solution of equations | dy| is less than the Interferometer angle measurement precision of permission, then initial value y
0be non-linear without fuzzy phase difference measurements solution of equations; If linearization is without the absolute value of fuzzy phase difference measurements solution of equations | dy| is more than or equal to the Interferometer angle measurement precision of setting, utilizes dy to initial value y
0revise, obtain modified value y '
0, with y '
0replace the initial value y in described step (13)
0, re-execute described step (13), (14), (15), until dy meets the Interferometer angle measurement precision of setting.To initial value y
0correction formula as follows:
y′
0=y
0+dy
This step finally obtains the Interferometer angle measurement value after calibrating
namely the angle of pitch after Interferometer angle measurement calibration is
position angle after calibration is
The random survey phase error of certain millimeter wave interferometer is about 20 °, the data phase aberration measurements of each stepping point being performed to about 5 seconds are average, the difference is smoothly utilized to carry out array calibration, after obtaining interferometer composition error item, utilize described angle measurement calibration algorithm, eliminate array element processing and composition error, systematic error caused by passage feeder line disparity error are installed, interfered angle measurement result accurately.When not calibrating interferometer, as shown in accompanying drawing 4a, there is obvious system deviation (1 ° of magnitude) in angle measurement result, this angle measurement result cannot be applied in real work; After calibrating interferometer as depicted in fig. 4b, angle measurement result is without obvious system deviation, and show as very little stochastic error (0.01 ° of magnitude), this angle measurement accuracy can meet interferometer real work demand.
Non-elaborated part of the present invention belongs to techniques well known.
Claims (1)
1. an interferometer array integrated correction method, is characterized in that: comprise the error model stage of setting up, obtain the measurement equation stage, solve the error parameter stage and obtain the angle measurement calibration result stage;
It is described that to set up error model stage etch as follows:
(1) on interferometer front, set up surving coordinate system, to choose on interferometer front an array element as interferometer measurement coordinate origin O, using desirable azimuth reference datum as Y-axis, using desirable pitching baseline as Z axis, by right-hand rule determination X-axis; The angle in objective definition incident direction and XOY face is pitching angle theta, and the angle in the projection of objective definition incident direction in XOY plane and XOZ face is position angle φ;
(2) in the interferometer measurement coordinate system set up in described step (1), if target original incident position angle is φ
0, the original incident angle of pitch is θ
0; If total M array element in interferometer array, 0th array element is positioned at initial point O place, i-th array element is arranged in desirable orientation baseline, and a jth array element is arranged in desirable pitching baseline, i, j=1,2 ..., M-1 and i ≠ j, there is the processing installation composition error in X-axis, Y-axis, Z axis three directions relative to ideal position in the i-th array element, is respectively Δ x
i, Δ y
i, Δ z
i, jth array element relative to ideal position exist X-axis, Y-axis, Z axis three directions processing install composition error be respectively Δ x
j, Δ y
j, Δ z
j; In addition, the Received signal strength of each array element enters interferometer signal processor by respective feeder line passage, strictly not isometric between feeder line passage, if with the feeder line passage of the 0th array element for reference length, then by i-th, the feeder line passage of j array element is set to Δ d respectively relative to the feeder line passage composition error of the 0th array element
i, Δ d
j; Respectively by Δ x
i, Δ y
i, Δ z
i, Δ d
i, θ
0, φ
0, Δ x
j, Δ y
j, Δ z
j, Δ d
j, θ
0, φ
0as six unknown parameters of azimuth reference datum, pitching baseline, set up the phase eikonal equation of interferometer azimuth reference datum, pitching baseline:
In formula, i, j=1,2 ..., M-1 and i ≠ j, L
ibe the i-th array element azimuth reference datum design length, L
jfor jth array element pitching baseline design length,
for the interferometer phase aberration measurements of azimuth reference datum,
for the interferometer phase aberration measurements of pitching baseline, λ is the wavelength of incoming signal;
(3) interferometer be arranged on elaborate servo mechanism or be erected on precise rotating platform, by adjusting the position of interferometer, the Z axis of the interferometer measurement coordinate system that described step (1) is set up overlaps with the direction rotating shaft of elaborate servo mechanism or precise rotating platform, the Y-axis of interferometer measurement coordinate system that described step (1) is set up and the pitching shaft parallel of elaborate servo mechanism or precise rotating platform; Keep target incident direction constant, make interferometer respectively around surving coordinate system Z axis and Y-axis anglec of rotation δ φ, δ θ, after coordinate conversion, by the interferometer azimuth reference datum phase eikonal equation in described step (2) and pitching baseline phase eikonal equation, by rotation of coordinate, obtaining take rotation amount δ φ, δ θ as azimuth reference datum, the pitching baseline six parameter phase difference model of basic known quantity:
In formula, i, j=1,2 ..., M-1 and i ≠ j;
It is as follows that equation stage etch is measured in described acquisition:
(4) control servo control mechanism or turntable rotate δ φ successively
k, δ θ
kaccurate step angle, k represents that kth time rotates, and 0≤k≤N-1, N represents number of revolutions, N>=0, resident certain hour on each step angle, gathers the interferometer phase aberration measurements that described current accurate step angle is corresponding;
(5) the interferometer phase aberration measurements gathered in every section of residence time in described step (4) is averaged, obtain corresponding step angle upper position baseline phase aberration measurements
pitching baseline phase aberration measurements
(6) described step (5) is processed the phase aberration measurements obtained
substitute into azimuth reference datum, pitching baseline six parameter phase difference model in described step (3), obtain azimuth reference datum, the non-linear six parameter phase difference measurements system of equations of pitching baseline:
In formula, i, j=1,2 ..., M-1 and i ≠ j, x
i=[Δ x
i, Δ y
i, Δ z
i, Δ d
i, θ
0, φ
0]
trepresent the parameter that the i-th array element is to be asked, x
j=[Δ x
j, Δ y
j, Δ z
j, Δ d
j, θ
0, φ
0]
trepresent the parameter that jth array element is to be asked, w
ik, w
jkrepresent the measurement noises that difference collection is introduced;
represent respectively and walk postrotational azimuth reference datum, pitching baseline six parameter phase eikonal equation through k, form is as follows:
In formula, 0≤k≤N-1, l represents the l time rotation, 0≤l≤k, δ φ
lbe the l time orientation to accurate step angle, δ θ
lbe that the l time pitching is to accurate step angle;
It is described that to solve error parameter stage etch as follows:
(7) non-linear six parameter phase difference measurements system of equations step (6) obtained, at initial value x
i0, x
j0place's linearization, to x in the i-th array element
iinitialize x
i0, x
i0=[Δ x
i0, Δ y
i0, Δ z
i0, Δ d
i0, θ
00, φ
00]
t, wherein Δ x
i0, Δ y
i0, Δ z
i0the initial value of composition error is installed in processing for X-axis, Y-axis, Z axis three directions, Δ d
i0for the initial value of feeder line passage composition error, θ
00for the initial value of pitching initial incidence angle, φ
00for the initial value of orientation initial incidence angle, at initial value x
i0place carries out first order Taylor expansion to the non-linear six parameter phase difference measurements system of equations of azimuth reference datum in described step (6), obtains azimuth reference datum lienarized equation group, is write as vector form:
Α in formula
ifor first-order partial derivative matrix, form is as follows:
for difference residual error, form is as follows:
dx
ifor initial value x
i0correction;
W
i=[w
i0, w
i1... w
iN-1]
tfor measurement noises;
In like manner, to x in jth array element
jinitialize x
j0, x
j0=[Δ x
j0, Δ y
j0, Δ z
j0, Δ d
j0, θ
00, φ
00]
t, wherein Δ x
j0, Δ y
j0, Δ z
j0the initial value of composition error is installed in processing for X-axis, Y-axis, Z axis three directions, Δ d
j0for the initial value of feeder line passage composition error, θ
00for the initial value of pitching initial incidence angle, φ
00for the initial value of orientation initial incidence angle, at x
j0place carries out first order Taylor expansion to the non-linear six parameter phase difference measurements system of equations of pitching baseline in described step (6), obtains pitching baseline lienarized equation group, is write as vector form:
In formula, Α
jfor first-order partial derivative matrix, form is as follows:
for difference residual error, form is as follows:
Dx
jfor initial value x
j0correction, i, j=1,2 ..., M-1 and i ≠ j; w
j=[w
j0, w
j1... w
jN-1]
tfor measurement noises;
(8) obtain the least square solution of azimuth reference datum, pitching baseline lienarized equation group in step (7) according to weighted least-squares method, this solution is initial value x
i0, x
j0correction dx
i, dx
j:
P in formula
i, P
jfor least square weighting matrix, i, j=1,2...M-1 and i ≠ j;
(9) according to the linearization six parameter phase difference measurements solution of equations that described step (8) obtains, judge to resolve end condition, if the absolute value of least square solution | dx
i|, | dx
j| be less than the interferometer array calibration accuracy of setting, then initial value x
i0, x
j0be non-linear six parameter phase difference measurements solution of equations; If the absolute value of least square solution | dx
i|, | dx
j| be more than or equal to the interferometer array calibration accuracy of setting, utilize least square solution to initial value x
i0, x
j0revise, obtain correction x '
i0, x '
j0, utilize revised x '
i0, x '
j0x in replacement step (7)
i0, x
j0, re-execute described step (7), (8), (9), until dx
i, dx
jmeet the requirement of interferometer array calibration accuracy, initial value x
i0, x
j0correction formula as follows:
x′
i0=x
i0+dx
i
x′
j0=x
j0+dx
j
I, j=1 in formula, 2 ..., M-1 and i ≠ j;
This step finally obtains composition error item corresponding to i-th array element
X
i0=[Δ x
i0, Δ y
i0, Δ z
i0, Δ d
i0, θ
00, φ
00]
t, the composition error item that a jth array element is corresponding
X
j0=[Δ x
j0, Δ y
j0, Δ z
j0, Δ d
j0, θ
00, φ
00]
t, namely the processing of the i-th array element in X-axis, Y-axis, Z axis three directions is installed composition error and is respectively Δ x
i0, Δ y
i0, Δ z
i0, feeder line passage composition error is Δ d
i, the processing of jth array element in X-axis, Y-axis, Z axis three directions is installed composition error and is respectively Δ x
j0, Δ y
j0, Δ z
j0, feeder line passage composition error is Δ d
j;
Described acquisition angle measurement calibration result stage etch is as follows:
(10) the whole composition error item Δ x will obtained in step (9)
i0, Δ y
i0, Δ z
i0, Δ d
i0, Δ x
j0, Δ y
j0, Δ z
j0, Δ d
j0, substitute into the phase difference measurements equation of the orientation of described step (3), pitching baseline, introduce noise item w during azimuth reference datum interferometry simultaneously
i, noise item w when pitching base line interference is measured
j, obtaining being with noisy take rotation amount δ φ, δ θ as azimuth reference datum, the pitching baseline six parameter phase eikonal equation of basic known quantity:
In formula,
for band noisy orientation difference, pitching difference that interferometry obtains;
(11) by carrying out computing to base line interference phase level difference measurements multiple in step (10), synthesize a virtual baseline and interfere phase level difference measurements, virtual baseline after synthesis is generally not more than half wavelength lambda/2 of incoming signal, utilize the phase level difference measurements after synthesis, realize resolving without fuzzy angle measurement, for azimuth reference datum phase difference measurements equation in step (10), if i=m and i=n, m<n, 1≤m≤M-1, during 1≤n≤M-1, synthesize a short virtual baseline by computing, computing operator is H
i(), then corresponding orientation virtual baseline phase difference measurements equation is written as:
Δ x
n0, Δ y
n0, Δ z
n0, Δ d
n0, L
n, w
nrespectively during corresponding i=n
Δ x
i0, Δ y
i0, Δ z
i0, Δ d
i0, L
i, w
i;
Δ x
m0, Δ y
m0, Δ z
m0, Δ d
m0, L
m, w
mrespectively during corresponding i=m
Δ x
i0, Δ y
i0, Δ z
i0, Δ d
i0, L
i, w
i;
In like manner, for pitching baseline phase difference measurements equation in described step (10), if j=p and j=q, p<q, synthesize a short virtual baseline when 1≤p≤M-1,1≤q≤M-1, computing operator is H
j(), then write out pitching virtual baseline phase difference measurements equation:
Δ x
q0, Δ y
q0, Δ z
q0, Δ d
q0, L
q, w
qrespectively during corresponding j=q
Δ x
j0, Δ y
j0, Δ z
j0, Δ d
j0, L
j, w
j;
Δ x
p0, Δ y
p0, Δ z
p0, Δ d
p0, L
p, w
prespectively during corresponding j=p
Δ x
j0, Δ y
j0, Δ z
j0, Δ d
j0, L
j, w
j;
(12) in the orientation that step (11) obtains, pitching virtual baseline phase difference measurements equation, differ from orientation
pitching differs
middlely deduct 2 π H respectively
i(Δ d
n0, Δ d
m0)/λ, 2 π H
j(Δ d
q0, Δ d
p0)/λ, and will subtract each other result reduction to (-π, π] principal value interval, obtain following without integral circumference ambiguity orientation, pitching baseline phase difference measurements system of equations:
In formula
<> symbol represent that reduction is arrived (-π, π] principal value interval; w
1=H
i(w
n, w
m), w
2=H
j(w
q, w
p); F in formula
1(θ, φ), f
2(θ, φ) is defined as follows:
(13) vector y=[θ, φ] is defined
t, to y initialize y
0,
herein
for angle of pitch initial value,
for position angle initial value, at y
0near to carrying out first order Taylor expansion without integral circumference ambiguity orientation, pitching baseline phase difference measurements system of equations in step (12), obtain linearization without fuzzy phase difference measurements system of equations, write as vector form:
u≈B|
y=y0dy+w
In formula, u is that form is as follows without fuzzy phase difference measurements residual error:
B is that form is as follows without fuzzy phase difference measurements system of equations first-order partial derivative matrix:
Dy=[d θ, d φ]
trepresent y
0correction variable; W=[w
1, w
2]
tfor measurement noises;
(14) what obtain described step (13) directly inverts without fuzzy phase difference measurements system of equations, can obtain initial value y in step (13)
0correction variable dy:
dy=B
-1u
(15) linearization obtained according to described step (14) is without fuzzy phase difference measurements solution of equations, judge to resolve end condition, if linearization is without the absolute value of fuzzy phase difference measurements solution of equations | dy| is less than the Interferometer angle measurement precision of permission, then initial value y
0be non-linear without fuzzy phase difference measurements solution of equations; If linearization is without the absolute value of fuzzy phase difference measurements solution of equations | dy| is more than or equal to the Interferometer angle measurement precision of setting, utilizes dy to initial value y
0revise, obtain modified value y '
0, with y '
0replace the initial value y in described step (13)
0, re-execute described step (13), (14), (15), until dy meets the Interferometer angle measurement precision of setting, to initial value y
0correction formula as follows:
y′
0=y
0+dy
This step finally obtains the Interferometer angle measurement value after calibrating
namely the angle of pitch after Interferometer angle measurement calibration is
position angle after calibration is
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410513334.9A CN104297718B (en) | 2014-09-29 | 2014-09-29 | Interferometer array integrated correction method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410513334.9A CN104297718B (en) | 2014-09-29 | 2014-09-29 | Interferometer array integrated correction method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104297718A true CN104297718A (en) | 2015-01-21 |
CN104297718B CN104297718B (en) | 2017-04-12 |
Family
ID=52317520
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410513334.9A Active CN104297718B (en) | 2014-09-29 | 2014-09-29 | Interferometer array integrated correction method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104297718B (en) |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106680769A (en) * | 2017-02-20 | 2017-05-17 | 上海微小卫星工程中心 | Property inspection method of positioning system for rotary baseline interferometer |
CN106769132A (en) * | 2016-11-30 | 2017-05-31 | 上海卫星工程研究所 | For the passive direction finding method of testing based on two-dimensional quadrature interferometer under the conditions of whole star |
CN107560553A (en) * | 2017-10-26 | 2018-01-09 | 清华大学 | The installation deviation scaling method of interferometer in multiaxis laser displacement measurement system |
CN108089441A (en) * | 2017-12-08 | 2018-05-29 | 中国科学院长春光学精密机械与物理研究所 | Clap machine secondary mirror six degree of freedom accurate adjusting mechanism calibration algorithm and storage medium in space |
CN108107424A (en) * | 2017-11-24 | 2018-06-01 | 北京遥感设备研究所 | A kind of missile-borne volume-search coverage method being resident based on servo |
CN108254718A (en) * | 2017-12-04 | 2018-07-06 | 西安电子科技大学 | Multiple target angle ambiguity solution method based on angle measuring interferometer |
CN109814063A (en) * | 2019-02-27 | 2019-05-28 | 中国电子科技集团公司第三十六研究所 | A kind of interferometer direction finding method and apparatus |
CN111025333A (en) * | 2019-12-05 | 2020-04-17 | 中国人民解放军63620部队 | Installation deviation calibration and correction method based on short baseline attitude determination of navigation satellite signal |
CN113552534A (en) * | 2021-08-05 | 2021-10-26 | 中国人民解放军战略支援部队航天工程大学 | Phase calibration method of rotating baseline interferometer based on pulse signal |
CN116148784A (en) * | 2023-02-24 | 2023-05-23 | 武汉大学 | Automatic phase calibration system and method for single-station lightning positioning system |
CN116359835A (en) * | 2023-05-15 | 2023-06-30 | 中国人民解放***箭军工程大学 | Y-type baseline interferometer direction finding device and direction finding method |
CN117347945A (en) * | 2023-12-04 | 2024-01-05 | 中国航天科工集团八五一一研究所 | Interferometer system direction finding method based on antenna array three-dimensional layout |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20070273576A1 (en) * | 2006-05-27 | 2007-11-29 | Bae Systems Information And Electronic Systems Integration, Inc. | Multipath resolving correlation interferometer direction finding |
CN102419432A (en) * | 2011-08-25 | 2012-04-18 | 电子科技大学 | Round array phase interferometer two-dimensional (2D) direction-finding method based on virtual base line |
CN102788975A (en) * | 2012-06-20 | 2012-11-21 | 西安空间无线电技术研究所 | Multi-array-element single-channel interference angle measurement device |
-
2014
- 2014-09-29 CN CN201410513334.9A patent/CN104297718B/en active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20070273576A1 (en) * | 2006-05-27 | 2007-11-29 | Bae Systems Information And Electronic Systems Integration, Inc. | Multipath resolving correlation interferometer direction finding |
CN102419432A (en) * | 2011-08-25 | 2012-04-18 | 电子科技大学 | Round array phase interferometer two-dimensional (2D) direction-finding method based on virtual base line |
CN102788975A (en) * | 2012-06-20 | 2012-11-21 | 西安空间无线电技术研究所 | Multi-array-element single-channel interference angle measurement device |
Non-Patent Citations (2)
Title |
---|
HAN YUE-TAO ET AL.: "Multi-baseline interferometer broadband real-time phase calibration method", 《IMAGE AND SINAL PROCESSING(CISP),2011 4TH INTERNATIONALCONGRESS》 * |
崔旭: "基于虚拟基线的干涉仪测向改进方法", 《通信技术》 * |
Cited By (23)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106769132B (en) * | 2016-11-30 | 2019-09-17 | 上海卫星工程研究所 | For under the conditions of whole star based on the passive direction finding test method of two-dimensional quadrature interferometer |
CN106769132A (en) * | 2016-11-30 | 2017-05-31 | 上海卫星工程研究所 | For the passive direction finding method of testing based on two-dimensional quadrature interferometer under the conditions of whole star |
CN106680769A (en) * | 2017-02-20 | 2017-05-17 | 上海微小卫星工程中心 | Property inspection method of positioning system for rotary baseline interferometer |
CN107560553B (en) * | 2017-10-26 | 2019-11-22 | 清华大学 | The installation deviation scaling method of interferometer in multiaxis laser displacement measurement system |
WO2019080888A1 (en) * | 2017-10-26 | 2019-05-02 | 清华大学 | Installation deviation calibration method for interferometer in multi-axis laser displacement measurement system |
CN107560553A (en) * | 2017-10-26 | 2018-01-09 | 清华大学 | The installation deviation scaling method of interferometer in multiaxis laser displacement measurement system |
US11022423B2 (en) | 2017-10-26 | 2021-06-01 | Tsinghua University | Method for calibrating an error of installation of an interferometer in a multi-axis laser displacement measurement system |
CN108107424A (en) * | 2017-11-24 | 2018-06-01 | 北京遥感设备研究所 | A kind of missile-borne volume-search coverage method being resident based on servo |
CN108107424B (en) * | 2017-11-24 | 2021-05-07 | 北京遥感设备研究所 | Missile-borne airspace search method based on servo dwell |
CN108254718A (en) * | 2017-12-04 | 2018-07-06 | 西安电子科技大学 | Multiple target angle ambiguity solution method based on angle measuring interferometer |
CN108254718B (en) * | 2017-12-04 | 2021-07-20 | 西安电子科技大学 | Multi-target angle ambiguity resolution method based on rotary interferometer |
CN108089441A (en) * | 2017-12-08 | 2018-05-29 | 中国科学院长春光学精密机械与物理研究所 | Clap machine secondary mirror six degree of freedom accurate adjusting mechanism calibration algorithm and storage medium in space |
CN108089441B (en) * | 2017-12-08 | 2021-01-15 | 中国科学院长春光学精密机械与物理研究所 | Calibration algorithm and storage medium for six-degree-of-freedom precision adjustment mechanism of secondary mirror of space shooting machine |
CN109814063A (en) * | 2019-02-27 | 2019-05-28 | 中国电子科技集团公司第三十六研究所 | A kind of interferometer direction finding method and apparatus |
CN111025333A (en) * | 2019-12-05 | 2020-04-17 | 中国人民解放军63620部队 | Installation deviation calibration and correction method based on short baseline attitude determination of navigation satellite signal |
CN113552534A (en) * | 2021-08-05 | 2021-10-26 | 中国人民解放军战略支援部队航天工程大学 | Phase calibration method of rotating baseline interferometer based on pulse signal |
CN113552534B (en) * | 2021-08-05 | 2022-02-01 | 中国人民解放军战略支援部队航天工程大学 | Phase calibration method of rotating baseline interferometer based on pulse signal |
CN116148784A (en) * | 2023-02-24 | 2023-05-23 | 武汉大学 | Automatic phase calibration system and method for single-station lightning positioning system |
CN116148784B (en) * | 2023-02-24 | 2024-04-12 | 武汉大学 | Automatic phase calibration system and method for single-station lightning positioning system |
CN116359835A (en) * | 2023-05-15 | 2023-06-30 | 中国人民解放***箭军工程大学 | Y-type baseline interferometer direction finding device and direction finding method |
CN116359835B (en) * | 2023-05-15 | 2023-08-15 | 中国人民解放***箭军工程大学 | Y-type baseline interferometer direction finding device and direction finding method |
CN117347945A (en) * | 2023-12-04 | 2024-01-05 | 中国航天科工集团八五一一研究所 | Interferometer system direction finding method based on antenna array three-dimensional layout |
CN117347945B (en) * | 2023-12-04 | 2024-03-22 | 中国航天科工集团八五一一研究所 | Interferometer system direction finding method based on antenna array three-dimensional layout |
Also Published As
Publication number | Publication date |
---|---|
CN104297718B (en) | 2017-04-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104297718A (en) | Interferometer array integrated correction method | |
EP2824479B1 (en) | Method and system for locating a moving vehicle | |
CN108375752B (en) | Amplitude-phase error single-radiation-source direction finding method based on all-angle search | |
CN108037374B (en) | Array antenna near field calibration method | |
CN104007413B (en) | Consider the array position error calibration method of information source azimuthal error | |
CN103064089B (en) | Method for calibrating satellite navigation digital multi-beam launching array antenna phase center | |
CN105445709B (en) | A kind of thinning array near field passive location amplitude and phase error correction method | |
CN110673196B (en) | Time difference positioning method based on multidimensional calibration and polynomial root finding | |
CN109975840B (en) | Positioning correction method for zero-setting antenna satellite navigation receiver | |
CN106980104B (en) | Signal direction of arrival self-correction method for sensor array | |
CN109507635A (en) | Utilize the array amplitude phase error evaluation method of two unknown orientation auxiliary sources | |
Lösler et al. | Terrestrial monitoring of a radio telescope reference point using comprehensive uncertainty budgeting: Investigations during CONT14 at the Onsala Space Observatory | |
CN107783156B (en) | Satellite tracking method | |
CN103926548A (en) | Method for quickly measuring precision of reflection face of radiotelescope | |
CN104793177B (en) | Microphone array direction-finding method based on least square method | |
Li et al. | Agricultural machinery GNSS/IMU-integrated navigation based on fuzzy adaptive finite impulse response Kalman filtering algorithm | |
Xue et al. | A closed-form of Newton method for solving over-determined pseudo-distance equations | |
Lindgren et al. | Distributed localization using acoustic Doppler | |
CN109541324A (en) | Element antenna phase center measurement method in a kind of battle array based on Compact Range | |
CN112698310A (en) | Multi-emission-channel calibration method and system for vehicle-mounted millimeter wave radar | |
Petrov et al. | Auto-calibration of automotive radars in operational mode using simultaneous localisation and mapping | |
CN109085556A (en) | A kind of high-frequency ground wave radar wave field forming method based on a second order p-ratio | |
CN112240957B (en) | Method for correcting amplitude-phase characteristics of antenna in satellite navigation interference direction finding | |
CN110445557B (en) | Deep space measurement and control interferometry large-aperture antenna pointing calibration method and device | |
CN109001670B (en) | Distributed passive positioning method and device combining time difference and angle |
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 |