CN108398670A - A kind of pulse signal direction-finding method and device based on angle measuring interferometer - Google Patents
A kind of pulse signal direction-finding method and device based on angle measuring interferometer Download PDFInfo
- Publication number
- CN108398670A CN108398670A CN201810118224.0A CN201810118224A CN108398670A CN 108398670 A CN108398670 A CN 108398670A CN 201810118224 A CN201810118224 A CN 201810118224A CN 108398670 A CN108398670 A CN 108398670A
- Authority
- CN
- China
- Prior art keywords
- pulse
- angle measuring
- measuring interferometer
- fuzzy
- finding
- 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
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/40—Means for monitoring or calibrating
-
- 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
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/021—Auxiliary means for detecting or identifying radar signals or the like, e.g. radar jamming signals
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
The present invention relates to a kind of pulse signal direction-finding method and device based on angle measuring interferometer, the direction-finding method include:The mathematical model of angle measuring interferometer is established, and according to the mathematical model, determines the boundary of pulse repetition frequency;Judge whether pulse signal meets the boundary inequality of pulse recurrence frequency, if being unsatisfactory for boundary inequality, is determined as low repetition pulse, and carry out low repetition pulse direction finding;If meeting boundary inequality, it is determined as high repetition pulse, and carry out high repetition pulse direction finding.The boundary that the present invention passes through calculating discrete pulses signal repetition rate, it may determine that whether pulse signal is low repetition pulse, if low repetition pulse, then no fuzzy interval is calculated according to two first pulse signals to gather, then the element in gathering without fuzzy interval is screened using posterior pulse signal again, and then high-precision orientation, pitching angular estimation are obtained, effectively increase the accuracy of direction finding result.
Description
Technical field
The present invention relates to a kind of pulse signal direction-finding method and device based on angle measuring interferometer, belongs to emitter Signals
Passive direction finding technical field.
Background technology
Important answer is suffered from electronic reconnaissance, radio control and communication to the passive direction finding of emitter Signals
With.In recent years, high-resolution direction finding technology has obtained significant progress, receives the extensive concern of many scholars both at home and abroad, and
It has carried out specializing in and proposing some effective methods for direction finding technology.
Typical direction-finding system is usually made of multiple channels, higher to obtain using MUSIC, ESPRIT, ML scheduling algorithm
Precision and resolving power.However, so many port number also gives volume, weight, power consumption and the back end signal processing of system
Harsh requirement is proposed, is not easy to realize in portable system.Another typical direction-finding system is interferometer, only due to it
A small amount of channel is needed, realizes simple and is widely used.However, there are lances for the direction finding precision reconciliation hazy condition of interferometer
Shield, higher direction finding precision needs longer baseline length, but also will appear fuzzy problem simultaneously.
In recent years, scholar proposes carries out high-precision direction finding in the way of angle measuring interferometer to target emanation source.Rotation
Interferometer can be equivalent to multiple baseline length and carry out direction finding, thus ambiguity solution easy to implement to target during rotation,
To under the conditions of Long baselines, ensure that the measurement accuracy to target bearing.However, existing method is believed both for continuous time
It number provides, preferable effect can be reached for the positioning of signal of communication.When incoming signal is discrete pulse signal,
Such as pulsed radar signal, existing method does not consider that the repetition situation of pulse signal, direction finding effect are poor usually at present.Therefore
For pulsed radar signal, the problem of how realizing accurate direction finding, be our urgent need to resolve.
Invention content
The object of the present invention is to provide a kind of pulse signal direction-finding method and device based on angle measuring interferometer, for solving
The repetition situation of pulse signal, the poor problem of direction finding effect are not considered.
In order to solve the above technical problems, the present invention provides a kind of pulse signal direction-finding method based on angle measuring interferometer,
Including following scheme:
Method scheme one:Steps are as follows:
The mathematical model of angle measuring interferometer is established, and according to the mathematical model, determines the boundary of pulse repetition frequency;
Judge whether pulse signal meets the boundary inequality of pulse recurrence frequency and sentence if being unsatisfactory for boundary inequality
It is set to low repetition pulse, and carries out low repetition pulse direction finding;If meeting boundary inequality, it is determined as high repetition pulse, goes forward side by side
Row high repetition pulse direction finding;
Wherein, the direction finding process for carrying out low repetition pulse is as follows:
Azimuth and Fuzzy Phase of two pulses of acquisition when reaching where corresponding angle measuring interferometer are poor, and paint
Make corresponding all equiphase difference curves;
The intersection point of the equiphase difference curve corresponding to two pulses is solved, and calculates the corresponding no fuzzy interval of each intersection point
Set;
When several subsequent pulses reach, the azimuth where corresponding to angle measuring interferometer and Fuzzy Phase are obtained
Difference, and draw out all equiphase difference curves;
Each section in gathering for no fuzzy interval judges the equiphase difference curve corresponding to succeeding impulse successively
Whether the section is passed through;If passing through the section, retains the section, otherwise reject the section in no fuzzy interval is gathered, directly
Next section is only remained in gathering to no fuzzy interval;
It handles to obtain no fuzzifying equation without fuzzy interval set progress deblurring to finally obtaining, according to the fuzzy side of the nothing
Journey obtains orientation and pitch angle estimated value, realizes low repetition pulse direction finding.
Method scheme two:On the basis of method scheme one, the step of mathematical model for establishing angle measuring interferometer, includes:
Using the center of angle measuring interferometer as coordinate origin, using angle measuring interferometer baseline initial time direction as x-axis
Direction obtains xoy Plane of rotations using the normal direction of initial baseline as y-axis direction, with perpendicular to the direction of xoy Plane of rotations
For z-axis direction;
According to the turning rate of angle measuring interferometer, baseline length, azimuth and the elevation angle, two days of angle measuring interferometer are calculated
Position vector and target vector of the line in space;
It is formed by time delay on two antennas of angle measuring interferometer in conjunction with pulse signal, calculates two channels of angle measuring interferometer
Between phase difference.
Method scheme three:On the basis of method scheme two, when determining the boundary of pulse repetition frequency, according to rotation
Phase difference between two channels of interferometer, is removed fuzzy by the way of difference ambiguity solution, obtains pulse recurrence frequency
Boundary inequality, corresponding mathematic(al) representation are:
Wherein, F is discrete pulses repetition rate, and d is the baseline length of angle measuring interferometer, wrFor the rotation of angle measuring interferometer
Corner rate, λ are the wavelength of discrete pulses.
Method scheme four:On the basis of method scheme three, the mathematic(al) representation of no fuzzifying equation is:
Wherein, φ=[φ1,φ2,...,φN]TIndicate that N number of pulse angle measuring interferometer within each impulse sampling time exists
, the ε poor without Fuzzy Phase of current locationiIndicate the azimuth where angle measuring interferometer, ε=[ε when i-th of pulse reaches1,
ε2,...,εN]TIndicate the vector of this N number of azimuth composition,Indicate that the pitch angle of angle measuring interferometer, θ indicate where target
Azimuth, β indicate the inconsistent phase between channel, ν=[ν1,ν2,...,νN]TIndicate the measurement error of phase difference.
Method scheme five:Further include being carried out to no fuzzifying equation using least square method on the basis of method scheme four
Processing is to reduce influence of the measurement error of inconsistent phase and phase difference between channel to orientation and pitch angle estimated value.
Method scheme six:On the basis of method scheme one, two, three, four or five, the direction finding process of high repetition pulse is carried out
It is as follows:
According to the N number of pulse received in pulse observation time, corresponding phase difference sequence is sought;
Ambiguity solution processing is carried out to phase difference sequence, obtains no Fuzzy Phase difference sequence;
Phase difference in no Fuzzy Phase difference sequence is changed according to trigonometric function form, obtain orientation and pitches angular estimation
Value realizes high repetition pulse direction finding.
Method scheme seven:On the basis of method scheme six, phase difference sequence is carried out at ambiguity solution using difference method
Reason, obtains no Fuzzy Phase difference sequence.
The present invention also provides a kind of pulse signal direction-finding device based on angle measuring interferometer, including following scheme:
Device scheme one:Including processor and memory, the processor is stored in for handling in the memory
Instruction is to realize following method:
The mathematical model of angle measuring interferometer is established, and according to the mathematical model, determines the boundary of pulse repetition frequency;
Judge whether pulse signal meets the boundary inequality of pulse recurrence frequency and sentence if being unsatisfactory for boundary inequality
It is set to low repetition pulse, and carries out low repetition pulse direction finding;If meeting boundary inequality, it is determined as high repetition pulse, goes forward side by side
Row high repetition pulse direction finding;
Wherein, the direction finding process for carrying out low repetition pulse is as follows:
Azimuth and Fuzzy Phase of two pulses of acquisition when reaching where corresponding angle measuring interferometer are poor, and paint
Make corresponding all equiphase difference curves;
The intersection point of the equiphase difference curve corresponding to two pulses is solved, and calculates the corresponding no fuzzy interval of each intersection point
Set;
When several subsequent pulses reach, the azimuth where corresponding to angle measuring interferometer and Fuzzy Phase are obtained
Difference, and draw out all equiphase difference curves;
Each section in gathering for no fuzzy interval judges the equiphase difference curve corresponding to succeeding impulse successively
Whether the section is passed through;If passing through the section, retains the section, otherwise reject the section in no fuzzy interval is gathered, directly
Next section is only remained in gathering to no fuzzy interval;
It handles to obtain no fuzzifying equation without fuzzy interval set progress deblurring to finally obtaining, according to the fuzzy side of the nothing
Journey obtains orientation and pitch angle estimated value, realizes low repetition pulse direction finding.
Device scheme two:On the basis of device scheme one, the step of mathematical model for establishing angle measuring interferometer, includes:
Using the center of angle measuring interferometer as coordinate origin, using angle measuring interferometer baseline initial time direction as x-axis
Direction obtains xoy Plane of rotations using the normal direction of initial baseline as y-axis direction, with perpendicular to the direction of xoy Plane of rotations
For z-axis direction;
According to the turning rate of angle measuring interferometer, baseline length, azimuth and the elevation angle, two days of angle measuring interferometer are calculated
Position vector and target vector of the line in space;
It is formed by time delay on two antennas of angle measuring interferometer in conjunction with pulse signal, calculates two channels of angle measuring interferometer
Between phase difference.
Device scheme three:On the basis of device scheme two, when determining the boundary of pulse repetition frequency, according to rotation
Phase difference between two channels of interferometer, is removed fuzzy by the way of difference ambiguity solution, obtains pulse recurrence frequency
Boundary inequality, corresponding mathematic(al) representation are:
Wherein, F is discrete pulses repetition rate, and d is the baseline length of angle measuring interferometer, wrFor the rotation of angle measuring interferometer
Corner rate, λ are the wavelength of discrete pulses.
Device scheme four:On the basis of device scheme three, the mathematic(al) representation of no fuzzifying equation is:
Wherein, φ=[φ1,φ2,...,φN]TIndicate that N number of pulse angle measuring interferometer within each impulse sampling time exists
, the ε poor without Fuzzy Phase of current locationiIndicate the azimuth where angle measuring interferometer, ε=[ε when i-th of pulse reaches1,
ε2,...,εN]TIndicate the vector of this N number of azimuth composition,Indicate that the pitch angle of angle measuring interferometer, θ indicate where target
Azimuth, β indicate the inconsistent phase between channel, ν=[ν1,ν2,...,νN]TIndicate the measurement error of phase difference.
Device scheme five:Further include being carried out to no fuzzifying equation using least square method on the basis of device scheme four
Processing is to reduce influence of the measurement error of inconsistent phase and phase difference between channel to orientation and pitch angle estimated value.
Device scheme six:On the basis of device scheme one, two, three, four or five, the direction finding process of high repetition pulse is carried out
It is as follows:
According to the N number of pulse received in pulse observation time, corresponding phase difference sequence is sought;
Ambiguity solution processing is carried out to phase difference sequence, obtains no Fuzzy Phase difference sequence;
Phase difference in no Fuzzy Phase difference sequence is changed according to trigonometric function form, obtain orientation and pitches angular estimation
Value realizes high repetition pulse direction finding.
Device scheme seven:On the basis of device scheme six, phase difference sequence is carried out at ambiguity solution using difference method
Reason, obtains no Fuzzy Phase difference sequence.
The beneficial effects of the invention are as follows:By the boundary for calculating discrete pulses signal repetition rate, it can be determined that go out arteries and veins
Rush whether signal is low repetition pulse, if low repetition pulse, then calculates no confusion region according to two first pulse signals
Between gather, then the element in gathering without fuzzy interval is screened using posterior pulse signal again, and then obtains height
The orientation of precision, pitching angular estimation effectively increase the accuracy of direction finding result by judging pulse repetition situation.
Description of the drawings
Fig. 1 is the flow diagram of the pulse signal direction-finding method based on angle measuring interferometer;
Fig. 2 is the mathematical model schematic diagram of angle measuring interferometer;
Fig. 3 is phase difference and angle measuring interferometer azimuth correspondence schematic diagram.
Specific implementation mode
In order to make the purpose , technical scheme and advantage of the present invention be clearer, below in conjunction with attached drawing and specific implementation
The present invention will be described in further detail for example.
The present invention provides a kind of pulse signal direction-finding device based on angle measuring interferometer, which includes processor
And memory, the processor are used to handle instruction stored in memory to realize that a kind of pulse based on angle measuring interferometer is believed
Number direction-finding method, the flow diagram of the direction-finding method is as shown in Figure 1, specifically include following steps:
(1) mathematical model of angle measuring interferometer is established.
Wherein, angle measuring interferometer usually utilizes two antennas along center rotation to complete direction finding.Establish space coordinates
As shown in Figure 2:Coordinate origin is located at the center of angle measuring interferometer, using angle measuring interferometer baseline initial time direction as x-axis
Direction, using the normal direction of initial baseline as y-axis direction, Plane of rotation is xoy planes, and z-axis is perpendicular to xoy planes.If rotation
The turning rate of interferometer is wr, baseline length d;The frequency of incoming signal is f, and azimuth where target is θ, pitch angle
ForThen position vector r of two antennas of interferometer in space1And r2It is expressed as:
r1=[dcos (wrt)/2,dsin(wrt)/2,0]T
r2=[- dcos (wrt)/2,-dsin(wrt)/2,0]T
The direction vector e in target emanation source is expressed as:
Then signal is formed by delay, τ on two antennas of interferometer1And τ2Respectively:
τ1=r1 Te/c
τ2=r2 Te/c
Wherein, c is the light velocity.
It follows that the phase difference (t) between two channels of interferometer is expressed as:
Wherein, λ is the wavelength of discrete pulses, and β indicates the inconsistent phase between channel.
The mathematic(al) representation of phase difference (t) describes the relationship of phase difference and time, correspondence schematic diagram such as Fig. 3
It is shown.The formula can also be expressed as phase difference and interferometer space different angle phase difference, i.e.,:
Wherein, εiIndicate the azimuth where angle measuring interferometer when i-th of pulse reaches.In actual measurement, pass through
Observe the phase difference of different locationi, to complete the estimation at azimuth and pitch angle, i=0,1 ..., N.
(2) according to the mathematical model of angle measuring interferometer, the boundary of radar pulse repetition frequency is determined.
Specifically, determining that the detailed process on the boundary of radar pulse repetition frequency is:
It is removed by way of difference ambiguity solution fuzzy, obtains correct phase, i.e.,:
Namely when the angle interval delta ε of the spatial sampling of angle measuring interferometer meets:
When the adjacent angle interval of angle measuring interferometer meets above formula, system can be removed by way of difference ambiguity solution
It is fuzzy, namely:
The repetition rate F of radar pulse needs to meet:
(3) boundary inequality is utilized, judges whether pulse signal meets the boundary inequality of pulse recurrence frequency, this time-division
For two kinds of situations:
Situation 1:When pulse recurrence frequency is unsatisfactory for repetition inequality, then it is determined as low repetition pulse.At this point, according to rotation
Become a cadre interferometer model the characteristics of, the direction finding process for carrying out low repetition pulse is as follows:
1) when the 1st discrete pulses reach, the azimuth where angle measuring interferometer is ε1, the fuzzy phase that measures
Potential difference is φ* 1, draw out all possible equiphase difference curve.
2) when the 2nd discrete pulses reach, the azimuth where angle measuring interferometer is ε2, the fuzzy phase that measures
Potential difference is φ* 2, draw out all possible equiphase difference curve.Solve all φ* 1And φ* 2Corresponding equiphase difference curve
Intersection point, and calculate the corresponding no fuzzy interval set of each intersection point.
If intersection point integrates as Φ=[χ1,χ2,...,χM], no fuzzy interval integrates as Ψ=[ρ1,ρ2,...,ρM], for these
Intersection point, φ1Corresponding fuzzy number is respectively n1=[n11,n12,...,n1M], φ2Corresponding fuzzy number is respectively n2=[n21,
n22,...,n2M], wherein M indicates that intersection point concentrates the number of intersection point, φ* i(i=1,2 ..., N) representative has fuzzy phase difference,
φi(i=1,2 ..., N) it represents without fuzzy phase difference.
Certainly, as other embodiments, corresponding no mould can also be sought out based on any other two pulses
Paste section set Ψ.
3) when i-th (i >=3), a discrete pulses reach, the azimuth where angle measuring interferometer is εi, measure
Fuzzy Phase difference is φ* i, draw out equiphase difference curve.For each section ρ in no fuzzy interval collection Ψi, sentence successively
Disconnected φ* iEquiphase difference curve whether pass through ρi, if so, retaining ρi, if it is not, then by ρiIt is rejected from Ψ.Update is without fuzzy
Interval Set Ψ updates fuzzy number.
4) and so on, after element number M=1 in Ψ, ambiguity removal, obtain no fuzzifying equation:
Wherein, φ=[φ1,φ2,...,φN]TIndicate that N number of pulse angle measuring interferometer within each impulse sampling time exists
, the ε poor without Fuzzy Phase of current locationiIndicate the azimuth where angle measuring interferometer, ε=[ε when i-th of pulse reaches1,
ε2,...,εN]TIndicate the vector of this N number of azimuth composition,Indicate that the pitch angle of angle measuring interferometer, θ indicate the side where target
Parallactic angle, β indicate the inconsistent phase between channel, ν=[ν1,ν2,...,νN]TIndicate the measurement error of phase difference.
Using the methods of least square, inconsistent phase (systematic error) β and phase measurement that reduce between channel are missed
The influence of poor ν azimuthals and pitch angle measurement, to finally obtain high-precision orientation, pitching angular estimation.Due to using most
The process that small square law handles fuzzifying equation belongs to the prior art, and details are not described herein again.
Situation 2:When pulse recurrence frequency meets repetition inequality, then it is determined as high repetition pulse.For Gao Zhongying arteries and veins
Punching, within the sampling time of each pulse, it is believed that angle measuring interferometer unchanged direction, therefore within each pulse period, it can be with
Phase difference in the position is obtained with traditional digital interference instrument method measurement, the specific direction finding process for carrying out high repetition pulse is such as
Under:
A) assume to have received N number of pulse in observation time, seek phase difference sequence φ*=[φ* 1,φ* 2,...,φ* N]T。
Specifically, for i-th of the pulse signal received in pulse observation time, two channel signals carry out FFT, obtain respectively
Phase difference onto incoming signal frequency point* i, phase difference sequence φ*=[φ* 1,φ* 2,...,φ* N]TMathematic(al) representation be:
Wherein, φ* 1,φ* 2,...,φ* NIndicate fuzzy phase difference, εiIndicate rotational interference when i-th of pulse reaches
Azimuth where instrument, ε=[ε1,ε2,...,εN]TIndicate the vector of this N number of azimuth composition,Indicate bowing for angle measuring interferometer
The elevation angle, θ indicate that the azimuth where target, β indicate that the inconsistent phase between channel, 2 π of β mod indicate β pairs of 2 π remainders.
B) under the conditions of Gao Zhongying, since differential phase is less than π, fuzzy problem is not present, uses difference ambiguity solution at this time
Mode to phase difference sequence φ*It is handled, restores original phase, obtained without fuzzy phase sequence φ=[φ1,
φ2,...,φN]T, have at this time:
Wherein, φ1,φ2,…,φNIndicate the phase difference after deblurring, after ambiguity solution, φ1,φ2,…,φNWith ε1,
ε2,...,εNCorrespondence can be indicated with Fig. 3.
C) phase difference after will be fuzzy changes according to trigonometric function form, obtains the estimated value of azimuth and pitch angle.
Specifically, sequence φ1,φ2,…,φNChange according to trigonometric function form, initial phase and amplitude separately include
The azimuth of incoming signal and the information of pitch angle.According toAzimuth and pitch angle
It is estimated as:
Wherein, N is pulse number, W (wr) it is φnSequence and exp (- j εn) sequence related coefficient, expression formula is:
Wherein, due to exp (- j εn) sequence periodicity, have:
Through the above steps, the direction finding process of high repetition pulse may be implemented, which can remove systematic error
Interferometer direction finding result influence, improve the accuracy of direction finding result.
The present invention is directed to measurement problem of the angle measuring interferometer to the arrival bearing of pulse signal, by marginal analysis, by it
It is divided into Gao Zhongying and low repetition both of which, and describes high repetition pulse and low repetition pulse in detail and realize high-precision direction finding
Detailed process.It the method achieve the high-precision direction finding of radar pulse signal, can be used for the angle estimation to radar signal, improve
Applicability of the angle measuring interferometer to radar signal.
Claims (10)
1. a kind of pulse signal direction-finding method based on angle measuring interferometer, which is characterized in that steps are as follows:
The mathematical model of angle measuring interferometer is established, and according to the mathematical model, determines the boundary of pulse repetition frequency;
Judge whether pulse signal meets the boundary inequality of pulse recurrence frequency and be determined as if being unsatisfactory for boundary inequality
Low repetition pulse, and carry out low repetition pulse direction finding;If meeting boundary inequality, it is determined as high repetition pulse, and carry out height
Repetition pulse direction finding;
Wherein, the direction finding process for carrying out low repetition pulse is as follows:
Azimuth and Fuzzy Phase of two pulses of acquisition when reaching where corresponding angle measuring interferometer are poor, and draw out
Corresponding all equiphase difference curves;
The intersection point of the equiphase difference curve corresponding to two pulses is solved, and calculates the corresponding no fuzzy interval collection of each intersection point
It closes;
When several subsequent pulses reach, the azimuth and Fuzzy Phase where the corresponding angle measuring interferometer of acquisition are poor, and
Draw out all equiphase difference curves;
Each section in gathering for no fuzzy interval judges whether is equiphase difference curve corresponding to succeeding impulse successively
Across the section;If passing through the section, retain the section, otherwise reject the section in no fuzzy interval is gathered, until nothing
Next section is only remained in the set of fuzzy interval;
It handles to obtain no fuzzifying equation without fuzzy interval set progress deblurring to finally obtaining, be obtained without fuzzifying equation according to this
To orientation and pitch angle estimated value, low repetition pulse direction finding is realized.
2. the pulse signal direction-finding method according to claim 1 based on angle measuring interferometer, which is characterized in that establish rotation
The step of mathematical model of interferometer includes:
Using the center of angle measuring interferometer as coordinate origin, using angle measuring interferometer baseline initial time direction as x-axis direction,
Using the normal direction of initial baseline as y-axis direction, obtain xoy Plane of rotations, using perpendicular to the direction of xoy Plane of rotations as z-axis
Direction;
According to the turning rate of angle measuring interferometer, baseline length, azimuth and the elevation angle, calculates two antennas of angle measuring interferometer and exist
The position vector and target vector in space;
It is formed by time delay on two antennas of angle measuring interferometer in conjunction with pulse signal, calculates between two channels of angle measuring interferometer
Phase difference.
3. the pulse signal direction-finding method according to claim 2 based on angle measuring interferometer, which is characterized in that determine pulse
When the boundary of signal repetition rate, according to the phase difference between the two of angle measuring interferometer channels, using the side of difference ambiguity solution
Formula removal is fuzzy, obtains the boundary inequality of pulse recurrence frequency, corresponding mathematic(al) representation is:
Wherein, F is discrete pulses repetition rate, and d is the baseline length of angle measuring interferometer, wrFor the rotation angle of angle measuring interferometer
Rate, λ are the wavelength of discrete pulses.
4. the pulse signal direction-finding method according to claim 3 based on angle measuring interferometer, which is characterized in that without fuzzy side
The mathematic(al) representation of journey is:
Wherein, φ=[φ1,φ2,...,φN]TIndicate N number of pulse within each impulse sampling time angle measuring interferometer current
, the ε poor without Fuzzy Phase of positioniIndicate the azimuth where angle measuring interferometer, ε=[ε when i-th of pulse reaches1,ε2,...,
εN]TIndicate the vector of this N number of azimuth composition,Indicate that the pitch angle of angle measuring interferometer, θ indicate the azimuth where target, β
Indicate the inconsistent phase between channel, ν=[ν1,ν2,...,νN]TIndicate the measurement error of phase difference.
5. the pulse signal direction-finding method according to claim 4 based on angle measuring interferometer, which is characterized in that further include adopting
No fuzzifying equation is handled with least square method to reduce the measurement error of inconsistent phase and phase difference between channel
Influence to orientation and pitch angle estimated value.
6. the pulse signal direction-finding method according to any one of claims 1-5 based on angle measuring interferometer, feature exist
In the direction finding process for carrying out high repetition pulse is as follows:
According to the N number of pulse received in pulse observation time, corresponding phase difference sequence is sought;
Ambiguity solution processing is carried out to phase difference sequence, obtains no Fuzzy Phase difference sequence;
Phase difference in no Fuzzy Phase difference sequence is changed according to trigonometric function form, obtains orientation and pitch angle estimated value,
Realize high repetition pulse direction finding.
7. the pulse signal direction-finding method according to claim 6 based on angle measuring interferometer, which is characterized in that use difference
Method carries out ambiguity solution processing to phase difference sequence, obtains no Fuzzy Phase difference sequence.
8. a kind of pulse signal direction-finding device based on angle measuring interferometer, which is characterized in that described including processor and memory
Processor is for handling the instruction being stored in the memory to realize following method:
The mathematical model of angle measuring interferometer is established, and according to the mathematical model, determines the boundary of pulse repetition frequency;
Judge whether pulse signal meets the boundary inequality of pulse recurrence frequency and be determined as if being unsatisfactory for boundary inequality
Low repetition pulse, and carry out low repetition pulse direction finding;If meeting boundary inequality, it is determined as high repetition pulse, and carry out height
Repetition pulse direction finding;
Wherein, the direction finding process for carrying out low repetition pulse is as follows:
Azimuth and Fuzzy Phase of two pulses of acquisition when reaching where corresponding angle measuring interferometer are poor, and draw out
Corresponding all equiphase difference curves;
The intersection point of the equiphase difference curve corresponding to two pulses is solved, and calculates the corresponding no fuzzy interval collection of each intersection point
It closes;
When several subsequent pulses reach, the azimuth and Fuzzy Phase where the corresponding angle measuring interferometer of acquisition are poor, and
Draw out all equiphase difference curves;
Each section in gathering for no fuzzy interval judges whether is equiphase difference curve corresponding to succeeding impulse successively
Across the section;If passing through the section, retain the section, otherwise reject the section in no fuzzy interval is gathered, until nothing
Next section is only remained in the set of fuzzy interval;
It handles to obtain no fuzzifying equation without fuzzy interval set progress deblurring to finally obtaining, be obtained without fuzzifying equation according to this
To orientation and pitch angle estimated value, low repetition pulse direction finding is realized.
9. the pulse signal direction-finding device according to claim 8 based on angle measuring interferometer, which is characterized in that establish rotation
The step of mathematical model of interferometer includes:
Using the center of angle measuring interferometer as coordinate origin, using angle measuring interferometer baseline initial time direction as x-axis direction,
Using the normal direction of initial baseline as y-axis direction, obtain xoy Plane of rotations, using perpendicular to the direction of xoy Plane of rotations as z-axis
Direction;
According to the turning rate of angle measuring interferometer, baseline length, azimuth and the elevation angle, calculates two antennas of angle measuring interferometer and exist
The position vector and target vector in space;
It is formed by time delay on two antennas of angle measuring interferometer in conjunction with pulse signal, calculates between two channels of angle measuring interferometer
Phase difference.
10. the pulse signal direction-finding device according to claim 9 based on angle measuring interferometer, which is characterized in that determine arteries and veins
When rushing the boundary of signal repetition rate, according to the phase difference between the two of angle measuring interferometer channels, using difference ambiguity solution
Mode removes fuzzy, obtains the boundary inequality of pulse recurrence frequency, corresponding mathematic(al) representation is:
Wherein, F is discrete pulses repetition rate, and d is the baseline length of angle measuring interferometer, wrFor the rotation angle of angle measuring interferometer
Rate, λ are the wavelength of discrete pulses.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810118224.0A CN108398670B (en) | 2018-02-06 | 2018-02-06 | Pulse signal direction finding method and device based on rotary interferometer |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810118224.0A CN108398670B (en) | 2018-02-06 | 2018-02-06 | Pulse signal direction finding method and device based on rotary interferometer |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108398670A true CN108398670A (en) | 2018-08-14 |
CN108398670B CN108398670B (en) | 2020-06-05 |
Family
ID=63095224
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810118224.0A Active CN108398670B (en) | 2018-02-06 | 2018-02-06 | Pulse signal direction finding method and device based on rotary interferometer |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108398670B (en) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109444838A (en) * | 2018-09-12 | 2019-03-08 | 上海无线电设备研究所 | One kind being based on the dual frequency solution velocity ambiguity method and system of pulse accumulation frame |
CN109596900A (en) * | 2018-12-06 | 2019-04-09 | 中国科学院遥感与数字地球研究所 | The test method and system of multifrequency antenna electric axis deviation |
CN109814064A (en) * | 2019-02-28 | 2019-05-28 | 中国电子科技集团公司第三十六研究所 | One kind being based on three array element L-type right angle battle array interferometer direction finding method and apparatus |
CN113552534A (en) * | 2021-08-05 | 2021-10-26 | 中国人民解放军战略支援部队航天工程大学 | Phase calibration method of rotating baseline interferometer based on pulse signal |
CN114487987A (en) * | 2022-02-23 | 2022-05-13 | 西安迅尔电子有限责任公司 | Efficient miniature ultra-wideband interferometer ambiguity resolution method |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100033377A1 (en) * | 2008-08-11 | 2010-02-11 | The Government of the United States as represented by the U. S. Navy | Sparse Interferometric Array for Radio Frequency Azimuth 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 |
CN103323811A (en) * | 2013-05-21 | 2013-09-25 | 西安电子科技大学 | Parameter estimation method based on virtual concentric annulus array |
CN105334490A (en) * | 2015-11-19 | 2016-02-17 | 西安空间无线电技术研究所 | Intervening and orientating method for broadband signal |
CN108196221A (en) * | 2017-12-20 | 2018-06-22 | 北京遥感设备研究所 | A kind of elimination of burst noise method based on multi-baseline interference instrument direction ambiguity section |
-
2018
- 2018-02-06 CN CN201810118224.0A patent/CN108398670B/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100033377A1 (en) * | 2008-08-11 | 2010-02-11 | The Government of the United States as represented by the U. S. Navy | Sparse Interferometric Array for Radio Frequency Azimuth 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 |
CN103323811A (en) * | 2013-05-21 | 2013-09-25 | 西安电子科技大学 | Parameter estimation method based on virtual concentric annulus array |
CN105334490A (en) * | 2015-11-19 | 2016-02-17 | 西安空间无线电技术研究所 | Intervening and orientating method for broadband signal |
CN108196221A (en) * | 2017-12-20 | 2018-06-22 | 北京遥感设备研究所 | A kind of elimination of burst noise method based on multi-baseline interference instrument direction ambiguity section |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109444838A (en) * | 2018-09-12 | 2019-03-08 | 上海无线电设备研究所 | One kind being based on the dual frequency solution velocity ambiguity method and system of pulse accumulation frame |
CN109596900A (en) * | 2018-12-06 | 2019-04-09 | 中国科学院遥感与数字地球研究所 | The test method and system of multifrequency antenna electric axis deviation |
CN109596900B (en) * | 2018-12-06 | 2020-06-19 | 中国科学院遥感与数字地球研究所 | Method and system for testing deviation of electric axis of multi-frequency antenna |
CN109814064A (en) * | 2019-02-28 | 2019-05-28 | 中国电子科技集团公司第三十六研究所 | One kind being based on three array element L-type right angle battle array interferometer direction finding method and apparatus |
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 |
CN114487987A (en) * | 2022-02-23 | 2022-05-13 | 西安迅尔电子有限责任公司 | Efficient miniature ultra-wideband interferometer ambiguity resolution method |
Also Published As
Publication number | Publication date |
---|---|
CN108398670B (en) | 2020-06-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108398670A (en) | A kind of pulse signal direction-finding method and device based on angle measuring interferometer | |
CN110375730B (en) | Indoor positioning navigation system based on IMU and UWB fusion | |
CN104076353B (en) | A kind of Area Objects echo beam center speed measurement method | |
CN108614258B (en) | Underwater positioning method based on single underwater sound beacon distance measurement | |
CN105842676B (en) | A kind of radar shadown analysis method and device | |
CN109934855A (en) | A kind of livewire work scene power components three-dimensional rebuilding method based on cloud | |
CN106597363A (en) | Pedestrian location method in indoor WLAN environment | |
CN108445444A (en) | Correct the direct localization method of volume Kalman filtering | |
CN108303043B (en) | Multi-sensor information fusion plant leaf area index detection method and system | |
CN109471101A (en) | A kind of Missile-borne SAR method for correcting image | |
CN104396321B (en) | Method of locating terminal and positioning device | |
CN108663654B (en) | 360-degree all-dimensional dynamic direction finding method based on continuous quantum pigeon group | |
CN105467370A (en) | Cross-range scaling method for precession object ISAR image of composite bistatic radar | |
CN110231592A (en) | Indoor orientation method, device, computer readable storage medium and terminal device | |
CN107861096A (en) | Least square direction-finding method based on voice signal reaching time-difference | |
CN108983169B (en) | Meter wave radar terrain correction method based on digital elevation model | |
CN114964259A (en) | Attitude calculation method and system based on infrared focal plane array time domain signal and infrared image | |
CN108562872A (en) | A kind of method of exceptional value detection when ultra-short baseline under-water acoustic locating calibration | |
CN110672115A (en) | Moving target track acquisition method based on multi-observation whistle digital telescope | |
CN110703205A (en) | Ultrashort baseline positioning method based on adaptive unscented Kalman filtering | |
CN116774264B (en) | Moving target positioning method based on low orbit satellite opportunistic signal Doppler | |
CN115826004B (en) | Three-star cooperative direct positioning method based on two-dimensional angle and time difference combination | |
CN111289944A (en) | Unmanned ship position and course measuring method based on UWB positioning | |
CN111127334A (en) | Real-time geometric correction method and system for SAR image based on RD plane pixel mapping | |
Chen et al. | IMM tracking of a 3D maneuvering target with passive TDOA system |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |