CN102540153A - Array amplitude and phase error correcting method based on interstation direct wave interference - Google Patents

Array amplitude and phase error correcting method based on interstation direct wave interference Download PDF

Info

Publication number
CN102540153A
CN102540153A CN2011104390900A CN201110439090A CN102540153A CN 102540153 A CN102540153 A CN 102540153A CN 2011104390900 A CN2011104390900 A CN 2011104390900A CN 201110439090 A CN201110439090 A CN 201110439090A CN 102540153 A CN102540153 A CN 102540153A
Authority
CN
China
Prior art keywords
direct
station
wave jamming
wave
frame
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
Application number
CN2011104390900A
Other languages
Chinese (zh)
Other versions
CN102540153B (en
Inventor
吴雄斌
刘斌
李伦
徐新安
沈志奔
陈骁峰
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Wuhan University WHU
Original Assignee
Wuhan University WHU
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Wuhan University WHU filed Critical Wuhan University WHU
Priority to CN 201110439090 priority Critical patent/CN102540153B/en
Publication of CN102540153A publication Critical patent/CN102540153A/en
Application granted granted Critical
Publication of CN102540153B publication Critical patent/CN102540153B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Abstract

Disclosed is an array amplitude and phase error correcting method based on interstation direct wave interference. Certain conditions are required if a high-frequency ground wave radar needs to receive the interstation direct wave interference with high signal to interference and noise ratio at a base band, and the array amplitude and phase error correcting method includes steps of setting working parameters of the high-frequency ground wave radar according to the conditions; and determining features of interstation direct wave interference; and searching multiframe distance echo spectra according to the features, determining information of distance elements with the interstation direct wave interference, constructing interstation direct wave interference receiving snapshot data and correcting array amplitude and phase errors. The array amplitude and phase error correcting method based on the interstation direct wave interference has the advantages that the array amplitude and phase error correcting method can realize active correction without extra auxiliary information sources, has excellent accuracy and robustness, can realize correction by the aid of asynchronous interstation direct wave interference, is simple to be realized, is low in calculation quantity, can realize long-term stable operation, and greatly reduces development cost and maintenance cost for the radar while improving exploration performances.

Description

A kind of array amplitude phase error calibration method based on direct-wave jamming between station
Technical field
The method that direct-wave jamming carries out the calibration of high-frequency ground wave radar array amplitude-phase consistency between station is utilized the present invention relates to a kind of.
Background technology
High-frequency ground wave radar uses MUSIC（Multiple Signal Classification）Ocean current DOA is carried out Deng super resolution algorithm（Direction of Arrival）Estimation, its general principle is that the single order marine echo Electromagnetic Scattering Theory proposed according to Barrick carries out Estimation of Spatial Spectrum extraction ocean current DOA information to the backscattering echo received.But this kind of super resolution algorithm is all based on certain signal model, only when model error is smaller, estimated result could accurately and stably, otherwise, and its performance is or even entirely ineffective by degradation.And in actual applications, it is difficult to realize preferably reception system, therefore model error is inevitable completely from hardware.For high-frequency ground wave radar, inconsistent each receiving channel magnitude-phase characteristics is to influence the key factor of MUSIC algorithm performances.Accordingly, it would be desirable to before ocean surface kinetic parameter inverting is carried out, be calibrated to the amplitude phase error of array.
Array calibration method has been broadly divided into source calibration and passive calibration.
In active calibration method, a known signal source is put to from array remote open area enough, transmission signal measures the amplitude and phase of each receiving channel output signal, deducts phase difference caused by array manifold position, you can obtain channel error information.This calibration method principle is simple, works well, is widely applied in practice.But the antenna array of high-frequency ground wave radar is nearer away from seawater, calibrated with auxiliary source, somewhat expensive and extremely inconvenience, it is difficult to long-term work.
In passive calibration methods, the signal source accurately known without direction directly utilizes the measured data and some prioris received（Such as array format）The disturbance parameter of the orientation of space information source and array is subjected to Combined estimator, orientation estimation can be completed simultaneously and is calibrated for error；In addition, passage element position and array element coupling error are equivalent to calibrate array disturbance parameter dependent on the amplitude phase error model in orientation by the amplitude phase error calibration method that orientation is relied on.Write in Wang Yongliang, Chen Hui etc.《Estimation of Spatial Spectrum theory and algorithm》（Publishing house of Tsing-Hua University 2004）The method is elaborated in one book.These methods need repeatedly complicated interative computation, and amount of calculation is very big, not necessarily meets requirement of real-time, and is possible to converge on local minimum, rather than global minimum, so as to obtain the result of mistake.
Wuhan University's wave propagation laboratory proposes the technology that list DOA spectrums in a kind of utilization marine echo click through row of channels calibration, this technology uses statistical method to a large amount of echo-signals, without complicated interative computation, amount of calculation is small, real-time, the Stability and veracity of passage calibration are improved, specific implementation details refer to No. 200610071360.6 Chinese invention patent applications " a kind of passive channel calibration method based on non-linear antenna array ".But the invention, for active calibration method, the Stability and veracity of its calibration result is still a bit weaker.
High-frequency ground wave radar typically uses FMICW（Frequency modulated interrupted continuous wave, abridge FMICW）System.Entitled " high-frequency ground wave radar target detection and tracking " that Rafaat Khan et al. are delivered（Target Detection and Tracking With a High Frequency Ground Wave Radar,IEEE Journal of Oceanic Engineering,1994,19(4):540~548）Paper in this is had a detailed description.Under FMICW waveform systems, marine echo（Including direct-wave jamming between sea clutter and station）Into after receiver, through mixing, LPF, A/D conversion and FFT can obtain as shown in Figure 2 apart from echo spectrum.As shown in figure 1, the local oscillation signal that A stations are produced is stood by A and launched, the local oscillation signal that B stations are produced is stood by B to be launched, A station transmission signals and B stations transmission signal（Reflected thing reflection）Stand and receive in B, B stations receive signal and obtained by mixing, LPF, A/D conversions and FFT apart from echo spectrum.Shown in Fig. 2 in echo spectrum, abscissa is distance element, and ordinate is power（Unit is dB）, which show direct-wave jamming between station.
For the situation of prior art, a kind of effective scheme of high-frequency ground wave radar array amplitude-phase consistency calibration is proposed, is the great difficult problem that the art is badly in need of solving.
The content of the invention
For the limitation of existing method, the purpose of the present invention is simple there is provided a kind of implementation method using direct-wave jamming between station, and with low cost, operand is small, the accurate and stable array amplitude-phase consistency calibration method of calibration result.
The technical scheme is that a kind of array amplitude phase error calibration method based on direct-wave jamming between station, comprises the following steps：
Step 1, the condition setting running parameter produced according to direct-wave jamming between station, produces direct-wave jamming between station；
Step 2, the feature of direct-wave jamming between station is determined according to step 1 gained running parameter；
Step 3, using step 2 determine station between the feature of direct-wave jamming the multiframe of reception is searched for frame by frame apart from echo spectrum, it is determined that between standing where direct-wave jamming apart from metamessage；
Step 4, snapshot data is received according to step 3 gained distance element information structuring, carries out array amplitude phase error calibration.
Moreover, the condition that direct-wave jamming is produced between being stood described in step 1 is as follows,
1）                                                ；
2）
Figure 355850DEST_PATH_IMAGE002
；
3）
Figure 2011104390900100002DEST_PATH_IMAGE003
,
Figure 822735DEST_PATH_IMAGE004
；
Wherein
Figure 2011104390900100002DEST_PATH_IMAGE005
,
Figure 533333DEST_PATH_IMAGE006
,
Figure 2011104390900100002DEST_PATH_IMAGE007
,,
Figure DEST_PATH_IMAGE009
,
Figure 715364DEST_PATH_IMAGE010
For the frequency sweep slope at the station of A, B two,
Figure DEST_PATH_IMAGE011
For the initial frequency at the station of A, B two,
Figure 985939DEST_PATH_IMAGE012
For the pulse period,For pulse frequency,
Figure 675678DEST_PATH_IMAGE014
For frequency sweep cycle,
Figure DEST_PATH_IMAGE015
Radar start time difference and satisfaction for the station of A, B two
Figure 373506DEST_PATH_IMAGE016
,
Figure DEST_PATH_IMAGE017
For the propagation time of direct-wave jamming between station,
Figure 515906DEST_PATH_IMAGE018
For direct wave propagation delay frequency difference between station,
Figure DEST_PATH_IMAGE019
For constant；CoefficientAnd coefficient
Figure DEST_PATH_IMAGE021
It is the condition that meets 1）And two minimum natural numbers of value；Coefficient
Figure 134417DEST_PATH_IMAGE022
,
Figure DEST_PATH_IMAGE023
Expression is rounded.
Moreover, the feature of direct-wave jamming is as follows between being stood described in step 2,
1）There is a direct-wave jamming between two stations at each station, two stand between the frequency of direct-wave jamming be respectivelyWith, wherein
Figure 508209DEST_PATH_IMAGE026
,
Figure DEST_PATH_IMAGE027
, and
Figure 409604DEST_PATH_IMAGE028
,
Figure DEST_PATH_IMAGE029
For pulse frequency；If
Figure 808355DEST_PATH_IMAGE030
, then direct-wave jamming is all present between two stations；If
Figure DEST_PATH_IMAGE031
, when
Figure 847987DEST_PATH_IMAGE032
When only one of which station between direct-wave jamming exist, direct-wave jamming is filtered out between another station, when
Figure DEST_PATH_IMAGE033
When two station between direct-wave jamming all exist；If
Figure 840344DEST_PATH_IMAGE034
, then direct-wave jamming between a station is only existed；Wherein
Figure DEST_PATH_IMAGE035
For the bandwidth of baseband filter；
2）Direct-wave jamming is on direct wave propagation delay frequency difference between station between the station at two stations
Figure 256413DEST_PATH_IMAGE018
Symmetrically；
3）Under the conditions of asynchronous, direct-wave jamming change over time between the station at two stations and gradually move, and between the station at two stations direct-wave jamming move it is in opposite direction, speed is identical；
4）Under the conditions of asynchronous, direct-wave jamming is at a time met the direct wave propagation delay frequency difference between station between the station at two stations
Figure 634305DEST_PATH_IMAGE018
Place.
Moreover, set the station of A, B two, determined using A stations as main website, described in step 3 between station where direct-wave jamming apart from metamessage when, A stations perform following steps,
Step 3.1, our station is read
Figure 110417DEST_PATH_IMAGE036
Frame pitch from echo spectrum,
Figure 573759DEST_PATH_IMAGE036
For the frame number of sampling；
Step 3.2, the bandwidth of baseband filter is judged
Figure 793519DEST_PATH_IMAGE035
And pulse frequency
Figure 25917DEST_PATH_IMAGE029
Relation；
WhenWhen, search for frame by frame
Figure DEST_PATH_IMAGE037
In the range of apart
Figure 498935DEST_PATH_IMAGE029
And it is dry to make an uproar than more than the default dry spectral peak made an uproar than threshold value PG, and the position of spectral peak is recorded, obtain data
Figure 381441DEST_PATH_IMAGE038
, wherein
Figure DEST_PATH_IMAGE039
；From data
Figure 484657DEST_PATH_IMAGE038
Middle one group data of the selection with gradual change feature are designated as apart from metamessage
Figure 427205DEST_PATH_IMAGE040
, terminate flow；
When
Figure DEST_PATH_IMAGE041
When, search for frame by frame
Figure 740506DEST_PATH_IMAGE042
With
Figure DEST_PATH_IMAGE043
In the range of apart
Figure 177435DEST_PATH_IMAGE029
And it is dry to make an uproar than more than the default dry spectral peak made an uproar than threshold value PG, searching for frame by frame
Figure 384425DEST_PATH_IMAGE044
In the range of it is dry make an uproar than more than the default dry spectral peak made an uproar than threshold value PG, and record the position of spectral peak, obtain data
Figure 373241DEST_PATH_IMAGE038
；Into step 3.3；
When
Figure DEST_PATH_IMAGE045
When, search for frame by frame
Figure 173838DEST_PATH_IMAGE044
In the range of it is dry make an uproar than more than the default dry spectral peak made an uproar than threshold value PG, and record the position of spectral peak, obtain data
Figure 398146DEST_PATH_IMAGE038
；Into step 3.3；
Step 3.3, the data that spectrum peak position is recorded in same time period are read at B stations, data
Figure 619360DEST_PATH_IMAGE046
Step 3.1 is performed by B stations and step 3.2 is obtained；
Step 3.4, data are compared
Figure 907253DEST_PATH_IMAGE038
And data
Figure 810618DEST_PATH_IMAGE046
, from dataIn select on direct wave propagation delay frequency difference between station
Figure 322819DEST_PATH_IMAGE018
Symmetrical and with gradual change feature data are designated as apart from metamessage
Figure 957062DEST_PATH_IMAGE040
, from dataIn select on direct wave propagation delay frequency difference between station
Figure 700207DEST_PATH_IMAGE018
Symmetrical and with gradual change feature data are designated as apart from metamessage
Figure DEST_PATH_IMAGE047
, apart from metamessage
Figure 201727DEST_PATH_IMAGE047
It is sent to B stations.
Moreover, it is as follows according to distance element information structuring reception snapshot data described in step 4,
Receive snapshot data,
WhereinFor frame number,
Figure 443801DEST_PATH_IMAGE050
ForiIndividual passagekDirect-wave jamming between the station that frame is received,NFor channel number.
Advantage of the invention is that：Extra auxiliary information source is not needed to can be achieved with source calibration, with good precision and robustness；It is adjustable using direct-wave jamming between asynchronous station, implementation method is simple；Operand is small；Can long-term stable operation；While detection performance is improved, the development cost and maintenance cost of radar are significantly reduced.
Brief description of the drawings
Fig. 1 is high-frequency ground wave radar fundamental diagram；
Fig. 2 is high-frequency ground wave radar apart from echo spectrogram； 
Fig. 3 is the FMICW waveform diagrams of the embodiment of the present invention；
Fig. 4 is where direct-wave jamming between the station of the embodiment of the present invention apart from metamessage schematic diagram；
Fig. 5 is that the embodiment of the present invention determines the flow chart where direct-wave jamming apart from metamessage between standing.
Embodiment
Technical solution of the present invention is described in detail below in conjunction with drawings and Examples.The array amplitude phase error calibration method that the embodiment of the present invention provides high-frequency ground wave radar is as follows：
The key of the present invention is the condition setting running parameter produced according to direct-wave jamming between station, then multiframe is searched for frame by frame apart from echo spectrum according to the feature of direct-wave jamming between station, it is determined that the direct-wave jamming between metamessage, build station between station where direct-wave jamming receives snapshot data.
Assuming that there is two high-frequency ground wave radar stations, the distance between A stations and B stations are, respectively have a high-frequency ground wave radar, using FMICW waveform systems, the initial frequency of A stations radar is
Figure 475342DEST_PATH_IMAGE052
, swept bandwidth is
Figure DEST_PATH_IMAGE053
, B station radar initial frequency be
Figure 147763DEST_PATH_IMAGE054
, swept bandwidth is
Figure DEST_PATH_IMAGE055
, A stand radar than B station radar shift to an earlier date
Figure 631965DEST_PATH_IMAGE015
（
Figure 415244DEST_PATH_IMAGE016
）Time starts, and remaining running parameter is all identical, and the wherein frequency sweep time is
Figure 425925DEST_PATH_IMAGE056
, frequency sweep cycle is
Figure 269248DEST_PATH_IMAGE014
, the pulse period is
Figure 240746DEST_PATH_IMAGE012
, pulse width is
Figure DEST_PATH_IMAGE057
, baseband filter band is a width of
Figure 827716DEST_PATH_IMAGE035
, as shown in figure 3, Fig. 3 abscissa is the time
Figure DEST_PATH_IMAGE059
。
The output frequency of direct-wave jamming after mixed is between the station that B stations are received
     
Figure 544833DEST_PATH_IMAGE060
                              （1）        
Wherein coefficient
Figure DEST_PATH_IMAGE061
,
Figure 559056DEST_PATH_IMAGE062
For positive integer,For pulse period corresponding pulse frequency,For mixer combination frequency interferences exponent number,
Figure 533145DEST_PATH_IMAGE064
Value is bigger, and the loss of direct-wave jamming is bigger between standing, in order that direct-wave jamming has stronger dry make an uproar than and for constant at base band between station, it is necessary to meet
Figure 128206DEST_PATH_IMAGE001
,
Figure DEST_PATH_IMAGE065
（Constant S size is determined by radar system and distance between sites）Requirement.
Order
Figure 313331DEST_PATH_IMAGE005
,
Figure 384055DEST_PATH_IMAGE006
, （Direct wave propagation delay frequency difference between standing）,
Figure 152608DEST_PATH_IMAGE008
,（1）Formula can be written as
Figure 508634DEST_PATH_IMAGE066
          （2）
Because, so（1）Only one of which expression formula may be met in two expression formulas of formula
Figure 942020DEST_PATH_IMAGE068
,
Figure 408905DEST_PATH_IMAGE065
Requirement.
The output frequency of direct-wave jamming after mixed is between the station that A stations are received can similarly be obtained
     （3）
By（2）Formula and（3）Knowable to formula, when
Figure 978558DEST_PATH_IMAGE070
When,
Figure DEST_PATH_IMAGE071
, due to
Figure 771064DEST_PATH_IMAGE017
Very little, under normal circumstances
Figure 550802DEST_PATH_IMAGE072
, so direct-wave jamming is also met between the station that A stations are received
Figure DEST_PATH_IMAGE073
,
Figure 821377DEST_PATH_IMAGE065
Requirement.
In summary, be in base band
Figure 511115DEST_PATH_IMAGE037
Place produce be more than it is certain it is dry make an uproar than station between the condition of direct-wave jamming be：1）
Figure 208944DEST_PATH_IMAGE001
；2）；3）
Figure 550244DEST_PATH_IMAGE003
,
Figure 219123DEST_PATH_IMAGE004
.Wherein。For the frequency sweep slope at the station of A, B two,
Figure 454560DEST_PATH_IMAGE011
For the initial frequency at the station of A, B two,
Figure 977945DEST_PATH_IMAGE012
For the pulse period,For pulse frequency,
Figure 400147DEST_PATH_IMAGE014
For frequency sweep cycle,
Figure 675271DEST_PATH_IMAGE015
Radar start time difference and satisfaction for the station of A, B two
Figure 194108DEST_PATH_IMAGE016
,For the propagation time of direct-wave jamming between station,For direct wave propagation delay frequency difference between station,For constant.The probability of direct-wave jamming appearance is between now standing
Figure 585720DEST_PATH_IMAGE074
.Coefficient
Figure 357367DEST_PATH_IMAGE020
And coefficient
Figure 183372DEST_PATH_IMAGE021
It is to meetAnd two minimum natural numbers of value；Coefficient
Figure DEST_PATH_IMAGE075
And
Figure 169094DEST_PATH_IMAGE016
, i.e.,
Figure 111642DEST_PATH_IMAGE076
,
Figure 424943DEST_PATH_IMAGE023
Expression is rounded.
Step 1, the condition setting running parameter produced according to direct-wave jamming between above-mentioned station, produces direct-wave jamming between station.The running parameter of required setting includes the frequency sweep slope at the station of A, B two
Figure 845560DEST_PATH_IMAGE010
, the initial frequency at the station of A, B two
Figure 193496DEST_PATH_IMAGE011
.Because the probability that direct-wave jamming occurs between standing is
Figure 306945DEST_PATH_IMAGE074
Running parameter can be arranged as required to during specific implementation, then start radar system, direct-wave jamming between standing whether is produced by design software procedure judges, radar system direct-wave jamming between producing station is restarted if direct-wave jamming between not producing station.
Step 2, the feature of direct-wave jamming between station is determined according to step 1 gained running parameter.
In embodiment, direct-wave jamming has following feature between the station of reception：
1）There is a direct-wave jamming between two stations at each station, two stand between the frequency of direct-wave jamming be respectively
Figure 107542DEST_PATH_IMAGE024
With
Figure 331850DEST_PATH_IMAGE025
, wherein
Figure 268713DEST_PATH_IMAGE026
,
Figure 553064DEST_PATH_IMAGE027
, and
Figure 840957DEST_PATH_IMAGE028
,
Figure 744322DEST_PATH_IMAGE029
For pulse frequency；If, then direct-wave jamming is all present between two stations；If
Figure 990944DEST_PATH_IMAGE031
, when
Figure 890767DEST_PATH_IMAGE032
When only one of which station between direct-wave jamming exist, direct-wave jamming is filtered out between another station, when
Figure 355681DEST_PATH_IMAGE033
When two station between direct-wave jamming all exist；If
Figure 657349DEST_PATH_IMAGE034
, then direct-wave jamming between a station is only existed；Wherein
Figure 893290DEST_PATH_IMAGE035
For the bandwidth of baseband filter；
2）Direct-wave jamming is on direct wave propagation delay frequency difference between station between the station at two stations
Figure 280409DEST_PATH_IMAGE018
Symmetrically；
3）Under the conditions of asynchronous, direct-wave jamming change over time between the station at two stations and gradually move, and between the station at two stations direct-wave jamming move it is in opposite direction, speed is identical；
4）Under the conditions of asynchronous, direct-wave jamming is at a time met the direct wave propagation delay frequency difference between station between the station at two stations
Figure 525576DEST_PATH_IMAGE018
Place.
Wherein
Figure 416172DEST_PATH_IMAGE035
For the bandwidth of baseband filter.
Step 3, using step 2 determine station between the feature of direct-wave jamming the multiframe of reception is searched for frame by frame apart from echo spectrum, it is determined that between standing where direct-wave jamming apart from metamessage.
Multiframe echo distance spectrum is scanned for according to the feature of direct-wave jamming between station, apart from metamessage where direct-wave jamming between the two station stations presented as shown in Figure 4, wherein abscissa is the time（Unit is hour, h）, ordinate is distance element, where the direct-wave jamming between metamessage and B stations station between the station of search gained A stations where direct-wave jamming apart from metamessage be in symmetrical structure.According to determine between station where direct-wave jamming apart from metamessage
Figure 88593DEST_PATH_IMAGE040
, construction reception snapshot data
Figure DEST_PATH_IMAGE077
.If the station of A, B two, using A stations as main website.Obtained in such as Fig. 5, embodiment apart from metamessage
Figure 572795DEST_PATH_IMAGE040
Detailed process be, by A station perform following steps：
Step 3.1, our station is read
Figure 621654DEST_PATH_IMAGE036
Frame pitch from echo spectrum,
Figure 632335DEST_PATH_IMAGE036
For the frame number of sampling；
Step 3.2, the bandwidth of baseband filter is judged
Figure 475657DEST_PATH_IMAGE035
And pulse frequency
Figure 571789DEST_PATH_IMAGE029
Relation；
When
Figure 158759DEST_PATH_IMAGE030
When, search for frame by frame
Figure 23947DEST_PATH_IMAGE037
In the range of apart
Figure 38171DEST_PATH_IMAGE029
And it is dry to make an uproar than more than the default dry spectral peak made an uproar than threshold value PG, and the position of spectral peak is recorded, obtain data
Figure 621599DEST_PATH_IMAGE038
, wherein
Figure 746681DEST_PATH_IMAGE039
；From data
Figure 731954DEST_PATH_IMAGE038
Middle selection has gradual change feature（It is in slow consecutive variations）One group of data be designated as apart from metamessage
Figure 917079DEST_PATH_IMAGE040
, terminate flow；
When
Figure 863169DEST_PATH_IMAGE041
When, search for frame by frame
Figure 916576DEST_PATH_IMAGE042
With
Figure 631722DEST_PATH_IMAGE043
In the range of apart
Figure 112382DEST_PATH_IMAGE029
And it is dry to make an uproar than more than the default dry spectral peak made an uproar than threshold value PG, searching for frame by frameIn the range of it is dry make an uproar than more than the default dry spectral peak made an uproar than threshold value PG, and record the position of spectral peak, obtain data
Figure 137287DEST_PATH_IMAGE038
；Into step 3.3；
When
Figure 706940DEST_PATH_IMAGE045
When, search for frame by frame
Figure 358501DEST_PATH_IMAGE044
In the range of it is dry make an uproar than more than the default dry spectral peak made an uproar than threshold value PG, and record the position of spectral peak, obtain data
Figure 154550DEST_PATH_IMAGE038
；Into step 3.3；
Step 3.3, the data that spectrum peak position is recorded in same time period are read at B stations
Figure 549759DEST_PATH_IMAGE046
, data
Figure 973918DEST_PATH_IMAGE046
Obtained by B stations according to the same manner, that is, perform step 3.1 and step 3.2, data are designated as obtained by simply recording the position of spectral peak
Figure 61960DEST_PATH_IMAGE046
；
Step 3.4, data are compared
Figure 204360DEST_PATH_IMAGE038
And data, from data
Figure 681925DEST_PATH_IMAGE038
In select on direct wave propagation delay frequency difference between station
Figure 940868DEST_PATH_IMAGE018
Symmetrical and with gradual change feature data are designated as apart from metamessage, from data
Figure 307576DEST_PATH_IMAGE046
In select on direct wave propagation delay frequency difference between station
Figure 706327DEST_PATH_IMAGE018
Symmetrical and with gradual change feature data are designated as apart from metamessage
Figure 11538DEST_PATH_IMAGE047
, apart from metamessage
Figure 987584DEST_PATH_IMAGE047
It is sent to B stations.B stations directly can be used apart from metamessagez ’(k) voluntarily build reception snapshot data, array amplitude phase error calibration is carried out, concrete mode is consistent with A stations.
During specific implementation, preset dry make an uproar and empirical value can be preset according to experiment than threshold value PG.
Step 4, snapshot data is received according to step 3 gained distance element information structuring, carries out array amplitude phase error calibration.
General each radar station has a receiving array, and each receiving array has multiple receiving channels.Apart from metamessage
Figure 47124DEST_PATH_IMAGE040
Represent thekPosition between frame station where direct-wave jamming, according to apart from metamessage
Figure 257657DEST_PATH_IMAGE040
The reception snapshot data of construction
Figure 986578DEST_PATH_IMAGE077
It is as follows：
Figure 65393DEST_PATH_IMAGE048
Wherein
Figure 438736DEST_PATH_IMAGE049
For frame number,
Figure 820170DEST_PATH_IMAGE040
ForkPosition between frame station where direct-wave jamming,
Figure 36388DEST_PATH_IMAGE050
ForiIndividual passagekDirect-wave jamming between the station that frame is received,NFor channel number.B stations are using apart from metamessage
Figure 794260DEST_PATH_IMAGE047
Bring above formula into.
According to existing estimate covariance Matrix Formula
Figure 146743DEST_PATH_IMAGE078
Obtain estimate covariance matrix
Figure DEST_PATH_IMAGE079
（Wherein
Figure 964658DEST_PATH_IMAGE080
Represent conjugate transposition）, to estimate covariance matrix
Figure 277959DEST_PATH_IMAGE079
Carry out feature decomposition and obtain characteristic vector
Figure DEST_PATH_IMAGE081
；Finally according to existing array amplitude phase error value formula
Figure 573942DEST_PATH_IMAGE082
Obtain array amplitude phase error value
Figure DEST_PATH_IMAGE083
, realize that array amplitude phase error is calibrated according to array amplitude phase error value.
Wherein：
Figure 921878DEST_PATH_IMAGE084
Figure DEST_PATH_IMAGE085
Figure 621677DEST_PATH_IMAGE086
Wherein,
Figure 422274DEST_PATH_IMAGE036
For the fast umber of beats for estimate covariance matrix（The number of arrayed data sampling in a period of time）,
Figure DEST_PATH_IMAGE087
For the orientation of direct-wave jamming between station,
Figure 521948DEST_PATH_IMAGE088
For complex constant,
Figure DEST_PATH_IMAGE089
For steering vector,
Figure 724390DEST_PATH_IMAGE090
For array amplitude phase error matrix.
Figure DEST_PATH_IMAGE091
Represent that 1 column vector is converted to diagonal matrix.
Figure 884107DEST_PATH_IMAGE092
It is characteristic vector
Figure 172000DEST_PATH_IMAGE081
'sNIndividual element,It is steering vector
Figure 809786DEST_PATH_IMAGE089
'sNIndividual element,
Figure 991369DEST_PATH_IMAGE083
It is array amplitude phase error matrixDiagonal onNIndividual element, these elements respectively withNIndividual passage correspondence.
Although the array amplitude-phase consistency calibration method that the present invention is described obtains success on high-frequency ground wave radar first, in essence, this method is also possible to the detection system applied to other use FMICW waveform systems.Specific embodiment described herein is only to spirit explanation for example of the invention.Those skilled in the art can be made various modifications or supplement to described specific embodiment or be substituted using similar mode, but without departing from spirit of the invention or surmount scope defined in appended claims.Direct-wave jamming is that calibration can be achieved between technical scheme utilizes asynchronous station, can similarly be calibrated simultaneously for direct-wave jamming between synchronous station, and implement simpler.

Claims (5)

1. a kind of array amplitude phase error calibration method based on direct-wave jamming between station, it is characterised in that comprise the following steps：
Step 1, the condition setting running parameter produced according to direct-wave jamming between station, produces direct-wave jamming between station；
Step 2, the feature of direct-wave jamming between station is determined according to step 1 gained running parameter；
Step 3, using step 2 determine station between the feature of direct-wave jamming the multiframe of reception is searched for frame by frame apart from echo spectrum, it is determined that between standing where direct-wave jamming apart from metamessage；
Step 4, snapshot data is received according to step 3 gained distance element information structuring, carries out array amplitude phase error calibration.
2. the array amplitude phase error calibration method according to claim 1 based on direct-wave jamming between station, it is characterised in that：The condition that direct-wave jamming is produced between being stood described in step 1 is as follows,
1）                                                ；
2）
Figure 848177DEST_PATH_IMAGE002
；
3）
Figure 884267DEST_PATH_IMAGE003
,
Figure 385786DEST_PATH_IMAGE004
；
Wherein
Figure 772905DEST_PATH_IMAGE005
,,
Figure 908668DEST_PATH_IMAGE007
,
Figure 456456DEST_PATH_IMAGE008
,
Figure 65291DEST_PATH_IMAGE009
,
Figure 114150DEST_PATH_IMAGE010
For the frequency sweep slope at the station of A, B two,
Figure 124831DEST_PATH_IMAGE011
For the initial frequency at the station of A, B two,
Figure 968154DEST_PATH_IMAGE012
For the pulse period,
Figure 64286DEST_PATH_IMAGE013
For pulse frequency,
Figure 385677DEST_PATH_IMAGE014
For frequency sweep cycle,
Figure 516444DEST_PATH_IMAGE015
Radar start time difference and satisfaction for the station of A, B two,For the propagation time of direct-wave jamming between station,
Figure 239177DEST_PATH_IMAGE018
For direct wave propagation delay frequency difference between station,
Figure 224451DEST_PATH_IMAGE019
For constant；Coefficient
Figure 409576DEST_PATH_IMAGE020
And coefficient
Figure 480300DEST_PATH_IMAGE021
It is the condition that meets 1）And two minimum natural numbers of value；Coefficient
Figure 409073DEST_PATH_IMAGE022
,
Figure 248853DEST_PATH_IMAGE023
Expression is rounded.
3. the array amplitude phase error calibration method according to claim 2 based on direct-wave jamming between station, it is characterised in that：The feature of direct-wave jamming is as follows between being stood described in step 2,
1）There is a direct-wave jamming between two stations at each station, two stand between the frequency of direct-wave jamming be respectively
Figure 729513DEST_PATH_IMAGE024
With, wherein
Figure 754417DEST_PATH_IMAGE026
,
Figure 324070DEST_PATH_IMAGE027
, and
Figure 850998DEST_PATH_IMAGE028
,
Figure 896314DEST_PATH_IMAGE029
For pulse frequency；If
Figure 166889DEST_PATH_IMAGE030
, then direct-wave jamming is all present between two stations；If
Figure 715682DEST_PATH_IMAGE031
, when
Figure 679090DEST_PATH_IMAGE032
When only one of which station between direct-wave jamming exist, direct-wave jamming is filtered out between another station, when
Figure 946124DEST_PATH_IMAGE033
When two station between direct-wave jamming all exist；If, then direct-wave jamming between a station is only existed；Wherein
Figure 423690DEST_PATH_IMAGE035
For the bandwidth of baseband filter；
2）Direct-wave jamming is on direct wave propagation delay frequency difference between station between the station at two stations
Figure 682633DEST_PATH_IMAGE018
Symmetrically；
3）Under the conditions of asynchronous, direct-wave jamming change over time between the station at two stations and gradually move, and between the station at two stations direct-wave jamming move it is in opposite direction, speed is identical；
4）Under the conditions of asynchronous, direct-wave jamming is at a time met the direct wave propagation delay frequency difference between station between the station at two stations
Figure 46749DEST_PATH_IMAGE018
Place.
4. the array amplitude phase error calibration method according to claim 1 based on direct-wave jamming between station, it is characterised in that：If the station of A, B two, determined using A stations as main website, described in step 3 between station where direct-wave jamming apart from metamessage when, A stations perform following steps,
Step 3.1, our station is read
Figure 49340DEST_PATH_IMAGE036
Frame pitch from echo spectrum,For the frame number of sampling；
Step 3.2, the bandwidth of baseband filter is judged
Figure 612357DEST_PATH_IMAGE035
And pulse frequency
Figure 729348DEST_PATH_IMAGE029
Relation；
WhenWhen, search for frame by frame
Figure 523309DEST_PATH_IMAGE037
In the range of apart
Figure 124055DEST_PATH_IMAGE029
And it is dry to make an uproar than more than the default dry spectral peak made an uproar than threshold value PG, and the position of spectral peak is recorded, obtain data, wherein
Figure 807157DEST_PATH_IMAGE039
；From data
Figure 180500DEST_PATH_IMAGE038
Middle one group data of the selection with gradual change feature are designated as apart from metamessage
Figure 686568DEST_PATH_IMAGE040
, terminate flow；
When
Figure 778152DEST_PATH_IMAGE041
When, search for frame by frame
Figure 536024DEST_PATH_IMAGE042
WithIn the range of apart
Figure 440843DEST_PATH_IMAGE029
And it is dry to make an uproar than more than the default dry spectral peak made an uproar than threshold value PG, searching for frame by frameIn the range of it is dry make an uproar than more than the default dry spectral peak made an uproar than threshold value PG, and record the position of spectral peak, obtain data；Into step 3.3；
WhenWhen, search for frame by frame
Figure 636146DEST_PATH_IMAGE044
In the range of it is dry make an uproar than more than the default dry spectral peak made an uproar than threshold value PG, and record the position of spectral peak, obtain data
Figure 561377DEST_PATH_IMAGE038
；Into step 3.3；
Step 3.3, the data that spectrum peak position is recorded in same time period are read at B stations
Figure 661051DEST_PATH_IMAGE046
, data
Figure 988127DEST_PATH_IMAGE046
Step 3.1 is performed by B stations and step 3.2 is obtained；
Step 3.4, data are compared
Figure 147844DEST_PATH_IMAGE038
And data
Figure 294792DEST_PATH_IMAGE046
, from data
Figure 221594DEST_PATH_IMAGE038
In select on direct wave propagation delay frequency difference between station
Figure 403177DEST_PATH_IMAGE018
Symmetrical and with gradual change feature data are designated as apart from metamessage
Figure 733795DEST_PATH_IMAGE040
, from dataIn select on direct wave propagation delay frequency difference between stationSymmetrical and with gradual change feature data are designated as apart from metamessage
Figure 111184DEST_PATH_IMAGE047
, apart from metamessage
Figure 612703DEST_PATH_IMAGE047
It is sent to B stations.
5. the array amplitude phase error calibration method according to claim 4 based on direct-wave jamming between station, it is characterised in that：It is as follows according to distance element information structuring reception snapshot data described in step 4,
Receive snapshot data
Figure 734243DEST_PATH_IMAGE048
,
WhereinFor frame number,
Figure 10952DEST_PATH_IMAGE050
ForiIndividual passagekDirect-wave jamming between the station that frame is received,NFor channel number.
CN 201110439090 2011-12-26 2011-12-26 Array amplitude and phase error correcting method based on interstation direct wave interference Expired - Fee Related CN102540153B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN 201110439090 CN102540153B (en) 2011-12-26 2011-12-26 Array amplitude and phase error correcting method based on interstation direct wave interference

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN 201110439090 CN102540153B (en) 2011-12-26 2011-12-26 Array amplitude and phase error correcting method based on interstation direct wave interference

Publications (2)

Publication Number Publication Date
CN102540153A true CN102540153A (en) 2012-07-04
CN102540153B CN102540153B (en) 2013-03-06

Family

ID=46347488

Family Applications (1)

Application Number Title Priority Date Filing Date
CN 201110439090 Expired - Fee Related CN102540153B (en) 2011-12-26 2011-12-26 Array amplitude and phase error correcting method based on interstation direct wave interference

Country Status (1)

Country Link
CN (1) CN102540153B (en)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103245933A (en) * 2013-04-17 2013-08-14 武汉大学 Motion compensation method for float-type high-frequency over-the-horizon radar
CN104020455A (en) * 2014-05-30 2014-09-03 中国电子科技集团公司第二十二研究所 Ground and sky wave array calibration method based on direct waves
CN105891792A (en) * 2015-02-16 2016-08-24 松下知识产权经营株式会社 Radar device
CN106291470A (en) * 2016-07-28 2017-01-04 中国船舶重工集团公司第七〇九研究所 A kind of based on the disturbance restraining method of feature during high-frequency ground wave radar ocean current result sky
CN107643514A (en) * 2017-09-12 2018-01-30 武汉大学 A kind of array calibration method of buoy base/shipborne radar based on direct wave
CN108051772A (en) * 2017-10-25 2018-05-18 西南电子技术研究所（中国电子科技集团公司第十研究所） Width be combined measurement incoming wave azimuth information method

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20040242176A1 (en) * 2001-06-25 2004-12-02 Vromans Johannes Sophia Device for down-transforming the frequency of signals
CN101349741A (en) * 2008-08-29 2009-01-21 西安电子科技大学 Phased array digital multi-beam forming machine for electron reconnaissance
CN101483273A (en) * 2009-02-24 2009-07-15 中国航天科技集团公司第五研究院第五○四研究所 Calibration method for amplitude and phase variable array antenna
CN101644765A (en) * 2009-04-23 2010-02-10 中国科学院声学研究所 Amplitude and phase error correction method used for linear array of underwater acoustic transducer

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20040242176A1 (en) * 2001-06-25 2004-12-02 Vromans Johannes Sophia Device for down-transforming the frequency of signals
CN101349741A (en) * 2008-08-29 2009-01-21 西安电子科技大学 Phased array digital multi-beam forming machine for electron reconnaissance
CN101483273A (en) * 2009-02-24 2009-07-15 中国航天科技集团公司第五研究院第五○四研究所 Calibration method for amplitude and phase variable array antenna
CN101644765A (en) * 2009-04-23 2010-02-10 中国科学院声学研究所 Amplitude and phase error correction method used for linear array of underwater acoustic transducer

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
严盟等: "阵列天线在近场条件下的幅相校正和阵元位置估计", 《雷达科学与技术》, vol. 03, no. 06, 31 December 2005 (2005-12-31) *
王利众等: "直接数字式频率合成幅相校正技术", 《舰船科学技术》, vol. 26, no. 03, 30 June 2004 (2004-06-30) *
贺志毅等: "一种宽带雷达幅相误差分析与校正方法", 《系统工程与电子技术》, vol. 25, no. 11, 31 December 2003 (2003-12-31) *

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103245933A (en) * 2013-04-17 2013-08-14 武汉大学 Motion compensation method for float-type high-frequency over-the-horizon radar
CN104020455A (en) * 2014-05-30 2014-09-03 中国电子科技集团公司第二十二研究所 Ground and sky wave array calibration method based on direct waves
CN104020455B (en) * 2014-05-30 2016-09-07 中国电子科技集团公司第二十二研究所 A kind of sky based on direct wave earthwave array calibration method
CN105891792A (en) * 2015-02-16 2016-08-24 松下知识产权经营株式会社 Radar device
CN106291470A (en) * 2016-07-28 2017-01-04 中国船舶重工集团公司第七〇九研究所 A kind of based on the disturbance restraining method of feature during high-frequency ground wave radar ocean current result sky
CN106291470B (en) * 2016-07-28 2018-08-24 中国船舶重工集团公司第七一九研究所 The disturbance restraining method of feature when a kind of ocean current result sky based on high-frequency ground wave radar
CN107643514A (en) * 2017-09-12 2018-01-30 武汉大学 A kind of array calibration method of buoy base/shipborne radar based on direct wave
CN107643514B (en) * 2017-09-12 2019-12-10 武汉大学 array calibration method of buoy-based/shipborne radar based on direct wave
CN108051772A (en) * 2017-10-25 2018-05-18 西南电子技术研究所（中国电子科技集团公司第十研究所） Width be combined measurement incoming wave azimuth information method
CN108051772B (en) * 2017-10-25 2021-12-07 西南电子技术研究所（中国电子科技集团公司第十研究所） Amplitude-phase combined method for measuring incoming wave azimuth information

Also Published As

Publication number Publication date
CN102540153B (en) 2013-03-06

Similar Documents

Publication Publication Date Title
CN103616679B (en) Based on difference beam modulation and the PD radar range finding angle-measuring method of wave form analysis
CN110412558A (en) The vehicle-mounted fmcw radar velocity ambiguity method of solution based on TDM MIMO
CN106405541B (en) Full-coherent continuous wave Doppler radar and distance and speed measuring method thereof
CN101825707B (en) Monopulse angular measurement method based on Keystone transformation and coherent integration
CN102540153A (en) Array amplitude and phase error correcting method based on interstation direct wave interference
CN104155648B (en) The high-frequency ground wave radar single snap MUSIC direction-finding method reset based on array data
CN101644773B (en) Real-time frequency domain super-resolution direction estimation method and device
CN103412301B (en) Metrewave radar angle measurement method
CN109581352B (en) Super-resolution angle measurement system based on millimeter wave radar
CN103713285B (en) Based on the distributed meter wave array radar angle-measuring method of information fusion
CN102914772B (en) Precession target two-dimensional imaging method based on equivalent scattering points
CN109510787A (en) Underwater acoustic channel lower linear FM signal method for parameter estimation and device
CN103901429A (en) Inverse synthetic aperture radar imaging method for maneuvering targets on basis of sparse aperture
CN104020438A (en) Direction-of-arrival estimation method based on sparse representation
CN107271955B (en) Time difference and scale difference estimation method for broadband linear frequency modulation signal
CN108490443A (en) Multiple submatrixes synthetic aperture sonar ω k imaging algorithms based on analytic solutions and NUFFT
CN104168232A (en) Method for determining multipath time delay and Doppler frequency shift in underwater acoustic channel
CN103558584A (en) Method and device for detecting direction of arrival
CN102394686B (en) Device and method for estimating angle of high-precision array antenna receiving system
JP2010175457A (en) Radar apparatus
CN108490425A (en) A kind of angle-measuring method of bistatic MIMO radar
CN104122548A (en) Angle measuring method for mechanical-scanning meter-wave array radar
CN105527624B (en) A kind of method and weather radar system of radar return dynamic estimation noise
CN102957472A (en) Method and device for receiving communication signal of thinned array antenna system
CN114035149B (en) Method for improving sensitivity of interferometer direction-finding system

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
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20130306

Termination date: 20201226