CN102419430B - Parallel-baseline-based two-dimensional direction finding method of round array phase interferometer - Google Patents

Parallel-baseline-based two-dimensional direction finding method of round array phase interferometer Download PDF

Info

Publication number
CN102419430B
CN102419430B CN 201110235023 CN201110235023A CN102419430B CN 102419430 B CN102419430 B CN 102419430B CN 201110235023 CN201110235023 CN 201110235023 CN 201110235023 A CN201110235023 A CN 201110235023A CN 102419430 B CN102419430 B CN 102419430B
Authority
CN
China
Prior art keywords
baseline
parallel
phase differential
phase
vector
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Fee Related
Application number
CN 201110235023
Other languages
Chinese (zh)
Other versions
CN102419430A (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.)
University of Electronic Science and Technology of China
Original Assignee
University of Electronic Science and Technology of China
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 University of Electronic Science and Technology of China filed Critical University of Electronic Science and Technology of China
Priority to CN 201110235023 priority Critical patent/CN102419430B/en
Publication of CN102419430A publication Critical patent/CN102419430A/en
Application granted granted Critical
Publication of CN102419430B publication Critical patent/CN102419430B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Preparing Plates And Mask In Photomechanical Process (AREA)

Abstract

The invention belongs to the technical field of radio monitoring technologies, and provides a parallel baseline based method for realizing two-dimensional wideband direction finding for a phase interferometer. The method comprises the steps of: figuring out a possible fuzzy number combination of a phase difference between two groups of parallel baselines by utilizing a linear relationship among fuzzy numbers of the phase difference between the two groups of parallel baselines in a uniform round array, thereby estimating the direction cosine of a possible incident signal; subsequently, calculating the corresponding phase differences among all the longest baselines; correlating actually measured phase difference vectors of all the longest baselines of the uniform round array; finding out a phase difference vector corresponding to the maximum correlation coefficient for estimation of a theoretical phase difference vector; and acquiring non-fuzzy phase difference vectors of all the longest baselines by resolving the phase ambiguity of all the longest baselines. According to the method provided by the invention, the deficiencies of other ambiguity resolving methods can be overcome; the estimation of the direction cosine approaches a CRLB (Cremer-Rao Lower Bound), so that the direction estimation of the incident signal achieves a very high direction finding precision; in addition, the calculation quantity in the invention is less.

Description

A kind of circle battle array phase-interferometer two dimension direction-finding method based on parallel baseline
Technical field
The invention belongs to the radio monitoring technical field, particularly the two dimension of the wide-band phase-interferometer in radio monitoring direction-finding method.
Background technology
Interferometer direction finding has the advantages such as algorithm is simple, highly sensitive, real-time is good, applicable antenna array is various informative, is widely used in the direction-finding system in electronic reconnaissance field.In order to improve the ability of direction finding precision and anti-multipath effect, require antenna aperture enough large, yet this must cause the fuzzy of direction finding.So in the phase-interferometer direction finding, the ambiguity of phase place is that ambiguity is to affect whether successful key issue of direction finding.
For addressing the above problem, various ambiguity solution methods are arisen at the historic moment.Present existing interferometer ambiguity solution method mainly contains: utilize long, the method ambiguity solution that short baseline combines (is seen document: Chen Qi, yellow brilliant, Song Shiqiong. the design studies [C] of circle battle array in nine yuan of uniform circular array interferometer direction finding systems. the 14 Annual Conference collection of thesis of electronic countermeasure branch of Chinese Institute of Electronics, 2005, (1): 717-721.), (see document: Gong Xiang iridium based on the phase differential ambiguity solution of diversity distance, Yuan Junquan, Sun Xiaochang. based on the ambiguity solution method research [J] of the phase differential variation value of diversity distance. signal is processed, 2003,19(4): 308-311) and the method (Chen Li of many baselines group cluster, Chen Hao, Xiao grants first. five yuan of uniform circular array interferometer weighting Direction Finding Algorithms and separate the condition [J] of phase ambiguity. and electronic countermeasure: 2004, (1): 8-12.).Although long and short baseline combined method simple and flexible requires the shortest baseline to be less than the half-wavelength of incoming signal, this has limited the maximum operation frequency of antenna.Phase differential ambiguity solution method based on diversity distance requires array element distance to satisfy certain irregular relation, and signal to noise ratio (S/N ratio) also there is certain restriction (Zhou Yaqiang, old flying, Huangfu may, Sun Zhongkang. Algorithm of Solving Multi-baseline Interferometer Phase DifferenceAmbiguity in Noisy Circumstance [J]. electronics and information journal: 2005,27(2) 259-261.).The method calculated amount of many baselines group cluster is larger, and can't provide a clear and definite cluster thresholding, has a strong impact on the understanding fuzzy performance.
In the various direction-finding methods of prior art, all have the problem that usable range is limited or calculated amount is large, the precision of simultaneous direction finding also is difficult to guarantee.
Summary of the invention
In various direction-finding methods of the prior art, exist usable range limited or calculated amount is large, the precision of simultaneous direction finding also is difficult to the technical matters that guarantees, and therefore an a kind of circle battle array phase-interferometer two dimension direction-finding method based on parallel baseline is provided.
The invention discloses a kind of circle battle array phase-interferometer two dimension direction-finding method based on parallel baseline, described method specifically comprises following steps:
The first step: choose two groups of parallel baselines in the planar array
Figure 2011102350237100002DEST_PATH_IMAGE002
With
Figure 2011102350237100002DEST_PATH_IMAGE004
, base length is respectively
Figure 2011102350237100002DEST_PATH_IMAGE006
, and
Figure 2011102350237100002DEST_PATH_IMAGE008
, the angle between two groups of parallel baselines is
Figure 2011102350237100002DEST_PATH_IMAGE010
Second step: the phase differential that calculates first group of parallel baseline
Figure 2011102350237100002DEST_PATH_IMAGE012
Phase differential with second group of parallel baseline
Figure 2011102350237100002DEST_PATH_IMAGE014
, and calculate
Figure 2011102350237100002DEST_PATH_IMAGE016
Root is the phase differential vector of long baseline
Figure 2011102350237100002DEST_PATH_IMAGE018
, wherein
Figure 38906DEST_PATH_IMAGE016
Be element number of array;
The 3rd step: the phase differential that utilizes second step to obtain
Figure 956047DEST_PATH_IMAGE012
With
Figure 676397DEST_PATH_IMAGE014
Calculate the possible fuzzy number of two groups of parallel baselines
Figure 2011102350237100002DEST_PATH_IMAGE020
,
Figure 2011102350237100002DEST_PATH_IMAGE022
, ,
Figure 2011102350237100002DEST_PATH_IMAGE026
:
Figure 2011102350237100002DEST_PATH_IMAGE028
；
Figure 2011102350237100002DEST_PATH_IMAGE030
；
Wherein
Figure 2011102350237100002DEST_PATH_IMAGE032
,
Figure 2011102350237100002DEST_PATH_IMAGE034
,
Figure 2011102350237100002DEST_PATH_IMAGE036
Be the incoming signal wavelength,
Figure 2011102350237100002DEST_PATH_IMAGE038
Expression rounds up,
Figure 2011102350237100002DEST_PATH_IMAGE040
Expression rounds downwards;
The 4th step: go on foot the fuzzy number that obtains by the 3rd
Figure 445507DEST_PATH_IMAGE020
,
Figure 259880DEST_PATH_IMAGE022
,
Figure 20025DEST_PATH_IMAGE024
,
Figure 224742DEST_PATH_IMAGE026
, have
Figure 2011102350237100002DEST_PATH_IMAGE042
Right
Figure 2011102350237100002DEST_PATH_IMAGE044
,
Figure 2011102350237100002DEST_PATH_IMAGE046
Right
Figure 2011102350237100002DEST_PATH_IMAGE048
, that is: Group , calculate respectively the direction cosine that N organizes incoming signal;
The 5th step: obtain according to the 4th step
Figure 2011102350237100002DEST_PATH_IMAGE054
The direction cosine of group incoming signal calculate the N group
Figure 2072DEST_PATH_IMAGE016
The longest baseline phase differential vector of root
Figure 2011102350237100002DEST_PATH_IMAGE056
The 6th step: the 5th step was obtained With the M root of the actual measurement phase differential vector of long baseline
Figure 477364DEST_PATH_IMAGE018
Carry out related calculation and select related coefficient when maximum corresponding phase differential vector be designated as
Figure DEST_PATH_IMAGE058
The 7th step: vectorial according to the phase differential that the 6th step obtained
Figure 169376DEST_PATH_IMAGE058
Obtain
Figure 841141DEST_PATH_IMAGE016
Root is the fuzzy number vector of long baseline
Figure DEST_PATH_IMAGE060
The 8th step: the fuzzy number vector that obtains with the 7th step obtains without fuzzy phase differential vector
Figure DEST_PATH_IMAGE062
Preferably, said method also further comprises:
The 9th step: find the solution the least square solution of direction cosine by the fuzzy phase differential vector of the nothing calculated, calculate the estimation of incoming signal direction cosine.
Preferably, said method also further comprises:
The tenth step: utilize the 9th to go on foot the estimation computer azimuth angle of the direction cosine that obtain and the estimation of the angle of pitch.
Preferably, above-mentioned planar array is uniform circular array.
Preferably, in above-mentioned the 6th step, the function of related operation is:
Figure DEST_PATH_IMAGE064
。
Preferably, obtain in above-mentioned 6 the 7th steps Root is the fuzzy number vector of long baseline Function is:
Figure DEST_PATH_IMAGE066
。
Beneficial effect of the present invention is: proposed a kind of circle battle array phase-interferometer two dimension direction-finding method based on parallel baseline, the method is by choosing two groups of parallel baselines in the uniform circular array, utilize the linear relationship between the fuzzy number of parallel baseline phase differential, calculate the combination of possible fuzzy number, corresponding with it phase differential vector is obtained in combination according to fuzzy number, then utilize related operation to find out the poor estimation of notional phase, by solution the longest baseline the longest baseline phase differential of phase ambiguity vector without fuzzy value, calculate at last the estimation of incoming signal position angle and the angle of pitch with least square method.The present invention not only can overcome the deficiency of other ambiguity solution algorithm, and high-precision incoming signal direction estimated value can be provided, and the direction cosine that obtain are estimated to approach preferably a carat Metro lower limit.In addition, calculated amount of the present invention is less, has guaranteed the real-time of broadband direction finding.
Description of drawings
Fig. 1 is the process flow diagram based on the circle battle array phase-interferometer two dimension direction-finding method of parallel baseline.
Fig. 2 is
Figure 724412DEST_PATH_IMAGE016
Unit's uniform circular array model.
Fig. 3 is two groups of parallel baseline Selection Models.
Fig. 4 is direction cosine
Figure DEST_PATH_IMAGE068
The comparison of estimation and carat Metro lower limit.
Fig. 5 is direction cosine The comparison of estimation and carat Metro lower limit.
Fig. 6 is direction finding standard deviation in position angle of the present invention.
Fig. 7 is angle of pitch direction finding standard deviation of the present invention.
Embodiment
Describe the specific embodiment of the present invention in detail below in conjunction with Figure of description.
The process flow diagram of justifying battle array phase-interferometer two dimension direction-finding method based on parallel baseline as shown in Figure 1, described method specifically comprises following steps:
The first step: choose two groups of parallel baselines in the uniform circular array
Figure DEST_PATH_IMAGE072
With
Figure DEST_PATH_IMAGE074
, base length is respectively
Figure 753198DEST_PATH_IMAGE006
, and
Figure 334352DEST_PATH_IMAGE008
, the angle between two groups of parallel baselines is
Figure 607201DEST_PATH_IMAGE010
Second step: the phase differential that calculates first group of parallel baseline
Figure 273806DEST_PATH_IMAGE012
Phase differential with second group of parallel baseline
Figure 352620DEST_PATH_IMAGE014
, and calculate Root is the phase differential vector of long baseline , wherein
Figure 385933DEST_PATH_IMAGE016
Be element number of array;
The 3rd step: the phase differential that utilizes second step to obtain With
Figure 493358DEST_PATH_IMAGE014
Calculate:
Figure 435907DEST_PATH_IMAGE028
（1）
Figure 811524DEST_PATH_IMAGE030
（2）
Wherein
Figure 169824DEST_PATH_IMAGE032
, ,
Figure 631210DEST_PATH_IMAGE036
Be the incoming signal wavelength,
Figure 556440DEST_PATH_IMAGE038
Expression rounds up,
Figure 718431DEST_PATH_IMAGE040
Expression rounds downwards;
The 4th step: go on foot and to obtain by the 3rd
Figure 717611DEST_PATH_IMAGE042
Right
Figure 930856DEST_PATH_IMAGE044
, in like manner can get
Figure 281066DEST_PATH_IMAGE046
Right
Figure 309065DEST_PATH_IMAGE048
, therefore can obtain Group Or Combination, will
Figure 330876DEST_PATH_IMAGE054
Group
Figure 834670DEST_PATH_IMAGE052
Respectively substitution:
Figure DEST_PATH_IMAGE078
（3）
Thereby can obtain
Figure 74021DEST_PATH_IMAGE054
The direction cosine of group incident direction , wherein , and
Figure DEST_PATH_IMAGE084
Be integer;
The 5th step: obtain according to the 4th step
Figure 244715DEST_PATH_IMAGE054
The direction cosine of group incoming signal are updated to the theoretical calculation formula that calculates phase differential, and are available
Figure 569517DEST_PATH_IMAGE054
Group
Figure 877002DEST_PATH_IMAGE016
The longest baseline phase differential vector of root
Figure 705280DEST_PATH_IMAGE056
The 6th step: the 5th step was obtained
Figure 440018DEST_PATH_IMAGE056
Substitution:
（4）
Corresponding phase differential vector was the estimation of notional phase difference vector when modus ponens (4) was maximum
Figure 160029DEST_PATH_IMAGE058
The 7th step: go on foot the poor estimation of notional phase that obtains with the 6th
Figure 111324DEST_PATH_IMAGE058
Substitution:
Figure 16963DEST_PATH_IMAGE066
（5）
Separate
Figure 50778DEST_PATH_IMAGE016
Root is phase differential fuzzy of long baseline, obtains fuzzy number vector
Figure 762382DEST_PATH_IMAGE060
The 8th step: with the fuzzy number substitution:
Figure DEST_PATH_IMAGE086
（6）
Thereby obtain without fuzzy phase differential vector
The 9th step: find the solution the least square solution of direction cosine by the fuzzy phase differential vector of the nothing calculated, it is as follows to find the solution formula:
Figure DEST_PATH_IMAGE088
（7）
Wherein
Figure DEST_PATH_IMAGE090
,
Figure DEST_PATH_IMAGE092
It is one
Figure DEST_PATH_IMAGE094
Matrix, its every delegation chooses combination corresponding to the corresponding array element of the longest baseline, supposes that certain row is corresponding to array element
Figure DEST_PATH_IMAGE096
And array element
Figure DEST_PATH_IMAGE098
Combination, then the 1st element of this row is
Figure DEST_PATH_IMAGE100
, the 2nd element is
Figure DEST_PATH_IMAGE102
,
Figure DEST_PATH_IMAGE104
Combination one total
Figure 140328DEST_PATH_IMAGE016
Kind,
Figure 723756DEST_PATH_IMAGE016
Be element number of array;
The tenth step: utilize the 9th to go on foot the estimation computer azimuth angle of the direction cosine that obtain and the estimation of the angle of pitch:
Figure DEST_PATH_IMAGE106
（8）
Principle of work of the present invention is as follows:
Consider M unit uniform circular array as shown in Figure 2, the array element radius is , as a reference point with the center of circle.If the incoming signal direction is , its frequency is
Figure DEST_PATH_IMAGE112
, wavelength is
Figure DEST_PATH_IMAGE114
, wherein,
Figure DEST_PATH_IMAGE116
Be the light velocity, then Individual array element with respect to the phase place of reference point is:
（9）
So the
Figure DEST_PATH_IMAGE122
Individual array element and
Figure DEST_PATH_IMAGE124
Phase differential between the individual array element can be expressed as:
Figure DEST_PATH_IMAGE126
（10）
Order ,
Figure DEST_PATH_IMAGE130
Be
Figure 115135DEST_PATH_IMAGE124
Individual array element and Base length between the individual array element, as shown in Figure 3, wherein the center of circle is true origin, direct north is
Figure DEST_PATH_IMAGE132
Axle is parallel to
Figure DEST_PATH_IMAGE134
Base direction is
Figure DEST_PATH_IMAGE136
Axle.Therefore, formula (10) can be written as:
Figure DEST_PATH_IMAGE138
（11）
Wherein,
Figure DEST_PATH_IMAGE140
Thereby can be poor in the hope of two groups of baseline notional phases , wherein Be two groups of parallel baseline sequence numbers.When the ratio (baseline wavelength ratio) of the longest base length and signal wavelength is larger, phase ambiguity can appear, so:
The phase differential of first group of parallel baseline can be expressed as:
Figure DEST_PATH_IMAGE146
（12）
The phase differential of second group of parallel baseline can be expressed as:
Figure DEST_PATH_IMAGE148
（13）
Wherein
Figure DEST_PATH_IMAGE150
Be the phase difference measurement value,
Figure DEST_PATH_IMAGE152
Be fuzzy number,
Figure DEST_PATH_IMAGE154
It is the angle of two groups of parallel baselines.
For first group of baseline, ideally have:
Figure DEST_PATH_IMAGE156
（14）
Put in order:
（15）
As seen
Figure 600047DEST_PATH_IMAGE044
Linear, wherein It is all known,
Figure 608455DEST_PATH_IMAGE032
,
Figure 599544DEST_PATH_IMAGE036
Be the incoming signal wavelength.By formula (15) can with
Figure 377007DEST_PATH_IMAGE022
Corresponding
Figure 792421DEST_PATH_IMAGE020
, because the impact of noise,
Figure 22545DEST_PATH_IMAGE020
May be not integer, do following processing:
Figure DEST_PATH_IMAGE162
（16）
In like manner can with
Figure 817326DEST_PATH_IMAGE026
Corresponding
Figure 449295DEST_PATH_IMAGE024
Got by formula (16)
Figure 38540DEST_PATH_IMAGE042
Right , in like manner can get
Figure 354431DEST_PATH_IMAGE046
Right
Figure 840908DEST_PATH_IMAGE048
, therefore can obtain
Figure 928949DEST_PATH_IMAGE050
Group
Figure 124877DEST_PATH_IMAGE076
Or
Figure DEST_PATH_IMAGE164
Combination.Utilize formula (12) second formulas and formula (13) second formulas to estimate
Figure 199143DEST_PATH_IMAGE054
The direction cosine of group incoming signal will
Figure 540126DEST_PATH_IMAGE054
Direction cosine substitution formula (11) is obtained
Figure 799069DEST_PATH_IMAGE054
Group The longest baseline phase differential of root
Figure DEST_PATH_IMAGE166
, wherein
Figure 103459DEST_PATH_IMAGE082
, and Be integer.Wherein have and only have one group of phase differential vector with
Figure 994372DEST_PATH_IMAGE016
The longest baseline phase differential of root
Figure 905171DEST_PATH_IMAGE018
Differ
Figure DEST_PATH_IMAGE168
The relation of integral multiple, for finding out this group phase differential vector, we will
Figure 321240DEST_PATH_IMAGE166
With
Figure 902394DEST_PATH_IMAGE018
Do suc as formula the related operation shown in (4), select so that formula (4) when maximum corresponding phase differential vector be the estimation of notional phase difference vector
Figure 237561DEST_PATH_IMAGE058
Will
Figure 904166DEST_PATH_IMAGE058
Substitution formula (5) is separated
Figure 920663DEST_PATH_IMAGE016
Root is phase differential fuzzy of long baseline, obtains fuzzy number vector , with the fuzzy number substitution:
Figure 800074DEST_PATH_IMAGE086
（17）
Thereby obtain
Figure 956905DEST_PATH_IMAGE016
Root the longest baseline without the Fuzzy Phase difference vector
Figure 777093DEST_PATH_IMAGE062
Utilize formula (7) to obtain the least-squares estimation of direction cosine
Figure DEST_PATH_IMAGE170
, utilize formula (8) thus calculate the position angle and the estimation of the angle of pitch
Figure DEST_PATH_IMAGE172
Thereby, finished the estimation of incoming signal position angle and the angle of pitch.
The below illustrates concrete effect of the present invention: adopt a kind of circle battle array phase-interferometer two dimension direction-finding method based on parallel baseline provided by the invention, at first be to choose that two groups of parallel baselines carry out rough measure in the uniform circular array, find out and separate all phase ambiguities of the longest baseline after the fuzzy number, then utilize the phase differential vector of the longest baseline behind the ambiguity solution by the estimation of least square method travel direction cosine, and then obtain the estimation of position angle and the angle of pitch.
Consider 9 yuan of uniform circular arrays, choose 81,72 with 67,40 liang the group parallel baselines, the angle of two groups of parallel baselines is Be 50 meters at the array radius, signal source is simple signal, and signal incident direction is
Figure DEST_PATH_IMAGE176
, the ratio (baseline wavelength ratio) of the longest base length and wavelength changes to 13 from 0.5, and signal to noise ratio (S/N ratio) is respectively in the situation of 10dB, 20dB, 30dB has carried out emulation experiment, wherein carries out Monte Carlo Experiment 1000 times under each baseline wavelength ratio.Fig. 4 with Figure 5 shows that different signal to noise ratio (S/N ratio)s under the curve that changes with the baseline wavelength ratio of direction cosine estimated value and carat Metro lower limit.Fig. 6 is under the different signal to noise ratio (S/N ratio)s with Fig. 7, the curve that the direction finding standard deviation changes with the baseline wavelength ratio.Can be seen by Fig. 4 and Fig. 5, under the different signal to noise ratio (S/N ratio)s of choosing and baseline wavelength ratio, the estimation of the direction cosine of method provided by the present invention approaches a carat Metro lower limit, thereby guaranteed the precision that incoming signal position angle and the angle of pitch are estimated, as shown in Figure 6 and Figure 7, when the baseline wavelength ratio greater than 1 the time, the angle measurement error of position angle and the angle of pitch can guarantee in 1 °.
The algorithm that the present invention proposes is not only applicable to the uniform circular array of interferometer direction finding system practical application, is applicable to too other planar arraies.Only need to find out two groups of parallel baselines in other planar arraies and get final product, antenna is structured the formation does not have special requirement.
The present invention expands to any new feature or any combination that discloses in this manual, and the either method that discloses or step or any combination of process.

Claims (4)

1. circle battle array phase-interferometer two dimension direction-finding method based on parallel baseline, described method specifically comprises following steps:
The first step: choose two groups of parallel baselines in the planar array
Figure DEST_PATH_IMAGE002
With
Figure DEST_PATH_IMAGE004
, base length is respectively
Figure DEST_PATH_IMAGE006
, and
Figure DEST_PATH_IMAGE008
, the angle between two groups of parallel baselines is
Figure DEST_PATH_IMAGE010
Second step: the phase differential that calculates first group of parallel baseline
Figure DEST_PATH_IMAGE012
Phase differential with second group of parallel baseline , and calculate Root is the phase differential of long baseline Vector, wherein
Figure 525784DEST_PATH_IMAGE016
Be element number of array;
The 3rd step: the phase differential that utilizes second step to obtain
Figure 380607DEST_PATH_IMAGE012
With
Figure 599492DEST_PATH_IMAGE014
, calculate two groups of fuzzy numbers that parallel baseline is possible
Figure DEST_PATH_IMAGE020
,
Figure DEST_PATH_IMAGE022
,
Figure DEST_PATH_IMAGE024
, :
Figure DEST_PATH_IMAGE028
；
Figure DEST_PATH_IMAGE030
；
Wherein
Figure DEST_PATH_IMAGE032
,
Figure DEST_PATH_IMAGE034
,
Figure DEST_PATH_IMAGE036
Be the incoming signal wavelength,
Figure DEST_PATH_IMAGE038
Expression rounds up,
Figure DEST_PATH_IMAGE040
Expression rounds downwards;
The 4th step: go on foot the fuzzy number that obtains by the 3rd
Figure 433850DEST_PATH_IMAGE020
,
Figure 248222DEST_PATH_IMAGE022
,
Figure 8367DEST_PATH_IMAGE024
,
Figure 213084DEST_PATH_IMAGE026
, calculate respectively the direction cosine that N organizes possible incoming signal, wherein
Figure DEST_PATH_IMAGE042
The 5th step: obtain according to the 4th step
Figure DEST_PATH_IMAGE044
The direction cosine of group incoming signal calculate
Figure 723568DEST_PATH_IMAGE044
Group
Figure 330130DEST_PATH_IMAGE016
The longest baseline phase differential vector of root
Figure DEST_PATH_IMAGE046
The 6th step: the 5th step was obtained Group
Figure 442936DEST_PATH_IMAGE016
The longest baseline phase differential vector of root
Figure 117631DEST_PATH_IMAGE046
With the actual measurement
Figure 641016DEST_PATH_IMAGE016
Root is the phase differential vector of long baseline
Figure 8543DEST_PATH_IMAGE018
Carry out related calculation and select related coefficient when maximum corresponding phase differential vector be designated as
Figure DEST_PATH_IMAGE048
The 7th step: vectorial according to the phase differential that the 6th step obtained Obtain Root is the fuzzy number vector of long baseline
Figure DEST_PATH_IMAGE050
The 8th step: the fuzzy number vector that obtains with the 7th step obtains without fuzzy phase differential vector
The 9th step: find the solution the least square solution of direction cosine by the fuzzy phase differential vector of the nothing calculated, calculate the estimation of incoming signal direction cosine;
The tenth step: utilize the 9th to go on foot the estimation computer azimuth angle of the direction cosine that obtain and the estimation of the angle of pitch.
2. the circle battle array phase-interferometer two dimension direction-finding method based on parallel baseline as claimed in claim 1 is characterized in that described planar array is uniform circular array.
3. the circle battle array phase-interferometer two dimension direction-finding method based on parallel baseline as claimed in claim 1 is characterized in that the function of related operation is in described the 6th step:
Figure DEST_PATH_IMAGE054
。
4. the circle battle array phase-interferometer two dimension direction-finding method based on parallel baseline as claimed in claim 1 is characterized in that obtaining in described the 7th step
Figure 794862DEST_PATH_IMAGE016
Root is the fuzzy number vector of long baseline
Figure 67711DEST_PATH_IMAGE050
Function is:
Figure DEST_PATH_IMAGE056
。
CN 201110235023 2011-08-17 2011-08-17 Parallel-baseline-based two-dimensional direction finding method of round array phase interferometer Expired - Fee Related CN102419430B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN 201110235023 CN102419430B (en) 2011-08-17 2011-08-17 Parallel-baseline-based two-dimensional direction finding method of round array phase interferometer

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN 201110235023 CN102419430B (en) 2011-08-17 2011-08-17 Parallel-baseline-based two-dimensional direction finding method of round array phase interferometer

Publications (2)

Publication Number Publication Date
CN102419430A CN102419430A (en) 2012-04-18
CN102419430B true CN102419430B (en) 2013-04-17

Family

ID=45943931

Family Applications (1)

Application Number Title Priority Date Filing Date
CN 201110235023 Expired - Fee Related CN102419430B (en) 2011-08-17 2011-08-17 Parallel-baseline-based two-dimensional direction finding method of round array phase interferometer

Country Status (1)

Country Link
CN (1) CN102419430B (en)

Families Citing this family (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102798839B (en) * 2012-09-13 2013-12-11 广州新软计算机技术有限公司 Device for positioning active to-be-measured-object in real time by using synchronous antenna arrays
CN103235281B (en) * 2013-04-03 2015-01-21 电子科技大学 Correlation interferometer direction-finding method based on phase difference increment
CN104122527B (en) * 2014-07-14 2016-08-17 中国人民解放军国防科学技术大学 A kind of round battle array phase-interferometer broadband based on look-up table instantaneous direction finding method
CN104749555B (en) * 2014-12-19 2017-05-03 中国航天科技集团公司第五研究院第五一三研究所 Phase difference direction finding and spatial spectrum direction finding combined direction-finding positioning system
CN107861095A (en) * 2017-10-10 2018-03-30 上海交通大学 A kind of single radio-frequency channel two dimensional wireless electricity direction-finding system
CN108872933B (en) * 2018-07-16 2022-02-08 电子科技大学 Single-station blind interferometer positioning method
CN109376422B (en) * 2018-10-18 2023-06-06 中国电子科技集团公司第三十六研究所 Uniform circular array optimal design evaluation method and device
CN110007267B (en) * 2019-01-29 2020-08-18 杭州电子科技大学 Uniform circular array interferometer direction finding ambiguity resolving method based on mixed base line
CN109633526B (en) * 2019-01-29 2020-09-01 杭州电子科技大学 Direction finding ambiguity resolving method of non-uniform circular array phase interferometer based on direction function
CN109814064B (en) * 2019-02-28 2023-04-14 中国电子科技集团公司第三十六研究所 Direction finding method and device based on three-array-element L-shaped rectangular array interferometer
CN112433192A (en) * 2020-11-05 2021-03-02 中国电子科技集团公司第二十九研究所 Low-cost high-precision direction-finding method for non-fixed-frequency pulse signal
CN114355280B (en) * 2022-03-18 2022-05-17 中国电子科技集团公司第二十九研究所 Multi-sensor composite array antenna arraying and multi-information fusion sorting angle measuring method
CN117075035B (en) * 2023-08-15 2024-04-30 湖南红船科技有限公司 Spin short baseline high-precision direction finding method, system, equipment, medium and terminal
CN117110980B (en) * 2023-10-23 2024-01-12 中国航天科工集团八五一一研究所 FPGA-based self-adaptive monopulse direction finding method
CN117347945B (en) * 2023-12-04 2024-03-22 中国航天科工集团八五一一研究所 Interferometer system direction finding method based on antenna array three-dimensional layout

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0438488A (en) * 1990-06-04 1992-02-07 Mitsubishi Heavy Ind Ltd Measuring method for azimuth of radio wave source utilizing flight data
US6784840B2 (en) * 2002-12-23 2004-08-31 Itt Manufacturing Enterprises, Inc. Method for determining azimuth and elevation angles using a single axis direction finding system
CN101109798B (en) * 2007-07-06 2011-04-20 哈尔滨工程大学 Direction finding method of accurate direction finding device for P/L waveband radiation source
CN101105525A (en) * 2007-07-06 2008-01-16 哈尔滨工程大学 Pure phase type broad frequency band microwave radiation source direction finding system and method
US7961147B1 (en) * 2008-07-25 2011-06-14 Rockwell Collins, Inc. Long baseline phase interferometer ambiguity resolution using frequency differences
CN101980043B (en) * 2010-09-15 2013-04-10 电子科技大学 Anti-receiver phase jump method for measuring directions of interference sources

Also Published As

Publication number Publication date
CN102419430A (en) 2012-04-18

Similar Documents

Publication Publication Date Title
CN102419430B (en) Parallel-baseline-based two-dimensional direction finding method of round array phase interferometer
CN102411136B (en) Phase interferometer direction finding method for ambiguity resolution by extension baselines
CN102419432B (en) Round array phase interferometer two-dimensional (2D) direction-finding method based on virtual base line
CN106487358B (en) A kind of maneuvering target turning tracking
CN107290709B (en) The relatively prime array Wave arrival direction estimating method decomposed based on vandermonde
CN103501538B (en) Based on the indoor orientation method of multipath energy fingerprint
CN107390199B (en) A kind of radar maneuvering target tracking waveform design method
CN106569171B (en) Direction of arrival angle estimation method based on two-layer hybrid array
CN103018730A (en) Distributed sub-array wave arrival direction estimation method
CN106569172B (en) Arrival direction estimation method
CN102209382A (en) Wireless sensor network node positioning method based on received signal strength indicator (RSSI)
CN104330767A (en) Interferometer direction-finding method based on phase fuzzy number search and least square fit
CN105911521B (en) Associated wireless electric signal complex envelope and the direct localization method of the over-the-horizon target of carrier phase information
CN109738861A (en) A kind of three-dimensional combined estimation method based on Wi-Fi channel state information
CN111601253B (en) Passive passive intelligent tracking and positioning method and system, storage medium and tracking and positioning terminal
CN109195110B (en) Indoor positioning method based on hierarchical clustering technology and online extreme learning machine
CN103353588B (en) Two-dimensional DOA (direction of arrival) angle estimation method based on antenna uniform planar array
CN106102163A (en) WLAN fingerprint positioning method based on RSS linear correlation Yu secondary weighted centroid algorithm
CN102540138A (en) Multi-base-line phase searching type two-dimensional spatial spectrum direction-measuring method
CN107749627B (en) Improved matching pursuit based intelligent power distribution network power flow Jacobian matrix estimation method
CN104793177B (en) Microphone array direction-finding method based on least square method
CN104931923A (en) Grid iterative estimation of signal parameters via rotational invariance techniques (ESPRIT), namely, extensible rapid estimation algorithm capable of being used for uniform circular array 2-dimensional direction of arrival (2D DOA)
CN104820204A (en) Weighted least square positioning method with reduced deviation
CN107132514A (en) Improve ESPRIT bistatic MIMO radar angle evaluation method
CN104535987A (en) Amplitude phase error self-correcting method applicable to uniform circular array acoustic susceptance 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

Granted publication date: 20130417

Termination date: 20140817

EXPY Termination of patent right or utility model