CN105974389B - Machine based on iterative processing sweeps metre wave radar Monopulse estimation method - Google Patents
Machine based on iterative processing sweeps metre wave radar Monopulse estimation method Download PDFInfo
- Publication number
- CN105974389B CN105974389B CN201610504677.8A CN201610504677A CN105974389B CN 105974389 B CN105974389 B CN 105974389B CN 201610504677 A CN201610504677 A CN 201610504677A CN 105974389 B CN105974389 B CN 105974389B
- Authority
- CN
- China
- Prior art keywords
- angle
- pulse
- target
- signal
- wave
- 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.)
- Active
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/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
-
- 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
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/02—Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
- G01S13/50—Systems of measurement based on relative movement of target
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 invention belongs to Radar Technology field, discloses a kind of machine based on iterative processing and sweep metre wave radar Monopulse estimation method, influence of the noise to signal can be reduced, reduce angle error;Including:Antenna is divided into two submatrixs, emits pulse signal;Coherent accumulation, the data after being accumulated are carried out to the echo-signal of reception;To the echo-signal of reception into horizontal phasing control, maximum amplitude is taken, and find the angle value θ of respective pulsesmax;In the angle, θmaxLeft and right take (L 1)/2 angle value θl, and corresponding amplitude is found out, to obtain two groups of new data;It is obtained and wave beam and difference beam by this two groups of data;Pair and wave beam and difference beam carry out and difference beam Monopulse estimation, obtain the off-axis angle of target;Off-axis angle value is assigned to θmax;It steps be repeated alternatively until that the estimation angle of target meets preset error condition.
Description
Technical field
The present invention relates to Radar Technology fields more particularly to a kind of machine based on iterative processing to sweep the survey of metre wave radar pulse
Angle method can be used for carrying out pulse method angle measurement to target when center of antenna normal in mechanical scanning radar is not aligned with target.
Background technology
For metre wave radar since signal decaying is small, detection range is remote, has in over-the-horizon detection, Anti-amyloid-β antibody etc. only
Special advantage, but simultaneously because the wavelength of metre wave radar is longer, wave beam is wider so that and its angular resolution is poor, and angle measurement accuracy is low.
Mechanical scanning radar is due at low cost, and the factors such as realization is simple are always by people's extensive use.Mechanic scan radar angle-measuring method is general
Using maximum-signal method, but the precision of this method is poor, and metre wave radar is swept especially for machine.Its electric size is smaller, wave beam
Wider, maximum-signal method angle measurement accuracy is not high.In order to improve angle measurement accuracy, Monopulse estimation technology may be used.
Zhang Guangyi academician is in " phased-array radar principle [M], Beijing:National Defense Industry Press gives in 1994 " one books
Monopulse estimation method.Monopulse estimation method refers to that can measure target actual position using the echo-signal of single pulse
A kind of angle-measuring method.In machine sweeps metre wave radar using width phase and difference beam single pulse method can be quickly obtain the standard of target
True position.Width phase and difference beam single pulse method refer to each pulse echo that radar is received after treatment, obtain each
Pulse and wave beam and difference beam, by calculating difference and ratio, table look-up known to each pulse target off-axis angle.Along with antenna
Reference angle just obtains the target angle that single pulse measures, and it is exactly the true of target to be averaged to the angle that all pulses measure
Position.
Since width phase and difference beam single pulse method realize that conveniently calculation amount is small, angle measurement accuracy is high, has in practical applications
Have great advantage.But the shortcomings of being swept wider metre wave radar wave beam, central beam misalignment target by machine is influenced, width phase and poor wave
Although beam Monopulse estimation method precision is higher than maximum-signal method, its angle measurement accuracy is affected by noise very big, and precision needs
It further increases, especially when antenna centre normal is not aligned with target, angle measurement accuracy is still very low, to influence radar
Resolving power and detection, track target accuracy.
Invention content
For the deficiency of above-mentioned prior art, the purpose of the present invention is to provide a kind of machines based on iterative processing to sweep metric wave
Radar Monopulse estimation method reduces angle error to reduce influence of the noise to signal, improves angle measurement accuracy.
The technical scheme is that:It is equivalent at the antenna surface being made of two submatrixs that machine is swept antenna, in antenna scanning
Subarray processing is carried out to the echo-signal of reception in the process, to the data formation obtained after processing and difference beam, is swept with opportunity
Array radar takes maximum amplitude, and find correspondence into horizontal phasing control to the echo-signal of reception to the data obtained after processing
The angle value of pulse, taking (L-1)/2 angle value in the left and right of the angle, (L is the pulse that beam pattern three dB bandwidth includes
Number), and corresponding and difference beam amplitude is found out, angle estimation is carried out using width phase and difference beam single burst algorithm, obtains one
A angle estimation value θe.Again in θeLeft and right take (L-1)/2 angle value, and find out corresponding and difference beam amplitude, reuse
Width phase and difference beam single burst algorithm carry out angle estimation, and an angle estimation value of getting back then proceedes to repeat above behaviour
Make, carries out multiple width phase and difference beam pulse iterative processing.
In order to achieve the above objectives, the embodiment of the present invention adopts the following technical scheme that:
A kind of machine based on iterative processing sweeps metre wave radar Monopulse estimation method, and described method includes following steps:
Step 1, it sets machine and sweeps the antenna of metre wave radar as the uniform line-array that is made of N number of array element, by the uniform line-array
It is divided into the left submatrix that array number is N/2 and the right submatrix that array number is N/2, wherein N is the even number more than or equal to 4;
Step 2, -90 ° of measurable angle range~90 ° that machine is swept to metre wave radar are divided into multiple wave positions, between adjacent wave position between
Every the first preset angle;Determine the mirror angular curve of each wave position, composition mirror angular curve table;
Step 3, machine sweeps metre wave radar K pulse of transmitting, is spaced the second preset angle between adjacent pulse, determines i-th of arteries and veins
The reference angle of punchingI=1,2 ..., K;Wherein, reference angle is the angle of center of antenna normal and horizontal reference plane;
Step 4, machine sweeps the echo-signal X that metre wave radar receives K pulse by antenna array, determines the mesh of i-th of pulse
Off-axis angle, i=1,2 ..., K are marked, and to echo-signal X into horizontal phasing control, obtains the echo-signal Y after phase adjustment, obtains
The corresponding pulse of maximum amplitude and the corresponding target off-axis angle of the pulse in the echo-signal Y after phase adjustment are taken, is denoted as
One target off-axis angle;Angle of the off-axis angle between target and antenna normal;
Step 5, the echo-signal T that left submatrix receives is obtained respectively1The echo-signal T received with right submatrix2, in the first mesh
(L-1)/2 angle value is taken on the left of mark off-axis angle, (L-1)/2 angle value, and adjacent angular are taken on the right side of first object off-axis angle
The second preset angle is divided between degree;And the echo-signal T received from left submatrix1In obtain and first object off-axis angle on the left of
(L-1)/2 angle value and the corresponding range value of (L-1)/2 angle value on right side form the first new echo-signal Z1, from
The echo-signal T that right submatrix receives2In obtain on the left of first object off-axis angle (L-1)/2 angle value and right side
(L-1)/2 the corresponding range value of angle value forms the second new echo-signal Z2;L is the pulse for including in 3dB beam angles
Number;
Step 6, according to the described first new echo-signal Z1With the second new echo-signal Z2Obtain with wave beam and difference beam, from
And difference and ratio are obtained, according to the difference and than obtaining the target off-axis angle θ of i-th of pulse with the mirror angular curve tablei, i=1,
2,...,K;By the target off-axis angle θ of i-th of pulseiWith the reference angle of i-th of pulseIt is added, obtains i-th of pulse to target
Measurement angle ψi, i=1,2 ... L;
Step 7, step 6 is repeated, until obtaining measurement angle of the L pulse to target, L measurement angle is made even
, the estimation angle of target is obtained;
Step 8, the first object off-axis angle estimation angle of the target being assigned in step 4;
Step 9, it is repeated in and executes step 5 to step 8, until the estimation angle of target meets preset error condition,
And by the estimation angle-determining obtained for the last time be target angle estimation value.
The present invention has the advantage that compared with prior art:The present invention carries out multiple pulse method to receiving data, will be single
The advantages of impulse method successive ignition, to which by pulse method, to center of antenna, whether the very sensitive disadvantage of alignment target weakens
, to improve the angle measurement accuracy of pulse method.
Description of the drawings
In order to more clearly explain the embodiment of the invention or the technical proposal in the existing technology, to embodiment or will show below
There is attached drawing needed in technology description to be briefly described, it should be apparent that, the accompanying drawings in the following description is only this
Some embodiments of invention for those of ordinary skill in the art without creative efforts, can be with
Obtain other attached drawings according to these attached drawings.
Fig. 1 is that a kind of machine based on iterative processing provided in an embodiment of the present invention sweeps metre wave radar Monopulse estimation method
Implementation process schematic diagram;
Fig. 2 is traditional width phase Monopulse estimation method and difference beam curve synoptic diagram;
Fig. 3 is the mirror angular curve schematic diagram of traditional width phase Monopulse estimation method;
Fig. 4 is the root-mean-square error that angle measurement is carried out using target of the method for the present invention and traditional width phase single pulse method pair
Correlation curve schematic diagram.
Specific implementation mode
Following will be combined with the drawings in the embodiments of the present invention, and technical solution in the embodiment of the present invention carries out clear, complete
Site preparation describes, it is clear that described embodiments are only a part of the embodiments of the present invention, instead of all the embodiments.It is based on
Embodiment in the present invention, it is obtained by those of ordinary skill in the art without making creative efforts every other
Embodiment shall fall within the protection scope of the present invention.
The embodiment of the present invention provides a kind of machine based on iterative processing and sweeps metre wave radar Monopulse estimation method, referring to Fig.1,
Described method includes following steps:
Step 1, it sets machine and sweeps the antenna of metre wave radar as the uniform line-array that is made of N number of array element, by the uniform line-array
It is divided into the left submatrix that array number is N/2 and the right submatrix that array number is N/2, wherein N is the even number more than or equal to 4.
Step 2, -90 ° of measurable angle range~90 ° that machine is swept to metre wave radar are divided into multiple wave positions, between adjacent wave position between
Every the first preset angle;Determine the mirror angular curve of each wave position, composition mirror angular curve table.
Step 2 is specially:
Measurable angle range is divided into a wave position for -90 °~90 ° by (2a) at interval of the first preset angle △ θ, wherein first is pre-
If angle θ=50.8/ △ (N-d/ λ), d is array element interval, and λ is signal wavelength, λ/2 d=;
(2b) for any one wave position, obtain the wave position left submatrix synthesis wave beam BlWith the wave position in right submatrix
Synthesize wave beam Br, according to the wave position left submatrix synthesis wave beam BlWith the wave position right submatrix synthesis wave beam BrObtain the wave
Position and wave beam and difference beam;
(2c) obtains the difference and ratio of the wave position according to the wave position and wave beam and difference beamAnd it is somebody's turn to do
The mirror angular curve of wave position;Wherein imag () expressions take imaginary-part operation;Using the angular range of the wave position as abscissa, with difference and ratio
For ordinate, the mirror angular curve of the wave position is drawn by the corresponding difference of the wave position measurable angle range and ratio;
(2d) obtains the mirror angular curve of all wave positions, and the mirror angular curve of all wave positions is formed mirror angular curve table.
Step 3, machine sweeps metre wave radar K pulse of transmitting, is spaced the second preset angle between adjacent pulse, determines i-th of arteries and veins
The reference angle of punchingI=1,2 ..., K;Wherein, reference angle is the angle of center of antenna normal and horizontal reference plane.
In step 3, machine sweeps metre wave radar and emits K pulse, and being spaced the second preset angle between adjacent pulse is specially:
Machine sweeps metre wave radar when antenna scanning works, and emits a pulse at interval of the second preset angle △ Φ, and in 3dB
Include L pulse in beam angle;Wherein, △ Φ=360*t/T, t are the pulse repetition period, and T is that machine sweeps turning for metre wave radar
Speed.
The reference angle of center of antenna normal and the angle of horizontal reference plane as i-th of pulseI=1,2 ..., K, by
The metre wave radar moment is swept in machine rotating, so center of antenna normal is variation.
Step 4, machine sweeps the echo-signal X that metre wave radar receives K pulse by antenna array, determines the mesh of i-th of pulse
Off-axis angle, i=1,2 ..., K are marked, and to echo-signal X into horizontal phasing control, obtains the echo-signal Y after phase adjustment, obtains
The corresponding pulse of maximum amplitude and the corresponding target off-axis angle of the pulse in the echo-signal Y after phase adjustment are taken, is denoted as
One target off-axis angle;Angle of the target off-axis angle between target and antenna normal.
Step 4 is specially:
(4a) machine sweeps the echo-signal X=AS+n that metre wave radar receives K pulse by antenna array, wherein A=[a
(θ1),...a(θi)...,a(θK)] indicate the target phase information that echo-signal includes, a (θi)=[1, exp (j2 π d/ λ sin
θi),...,exp(j2π(N-1)d/λsinθi)]TFor the target off-axis angle vector of i-th of pulse, i=1,2 ..., K;D is battle array
Member interval, λ is signal wavelength, indicates that the transposition of vector, exp indicate that j represents imaginary unit using e as the exponential depth at bottom;
S=[S1,...Si...,SK]TIndicate the complex envelope of echo-signal, SiIndicate the complex envelope information of i-th of pulse,Wherein fdIndicate the Doppler frequency of target, fd=2Vf0/ c, V indicate that target sweeps metre wave radar with respect to machine
Radial velocity, f0Indicate that the centre frequency of radar emission signal, c represent the light velocity, t is the pulse repetition period;N represent mean value as
0, N × K rank white Gaussian noise matrixes that variance is 1;
(4b), into horizontal phasing control, obtains the echo-signal Y=a (θ after phase adjustment to echo-signal X0) X, it is assumed that
The actual angle of target is θ0a(θ0)=[1, exp (j2 π d/ λ sin (θ0)),…,exp(j2π(N-1)d/λsin(θ0))]T;
It should be noted that the actual angle of the target described in the present embodiment is the mesh estimated in advance according to prior information
Target angle, belongs to the scope of rough estimate, and the present invention is accurately estimated the angle of target.
(4c) obtains the corresponding pulse of maximum amplitude and pulse corresponding first in the echo-signal Y after the adjustment of position
Target off-axis angle θmax。
Step 5, the echo-signal T that left submatrix receives is obtained respectively1The echo-signal T received with right submatrix2, in the first mesh
(L-1)/2 angle value is taken on the left of mark off-axis angle, (L-1)/2 angle value, and adjacent angular are taken on the right side of first object off-axis angle
Second preset angle is divided between degree;And the echo-signal T received from left submatrix1In obtain and first object off-axis angle
The corresponding range value of (L-1)/2 angle value in left side forms the first new echo-signal Z1, from the echo-signal T of right submatrix reception2
In obtain corresponding with (L-1)/2 angle value on the right side of first object off-axis angle range value and form the second new echo-signal Z2;L
For the pulse number for including in 3dB beam angles.
Wherein, T1For 1 × kth moment battle array, T2For 1 × kth moment battle array.
Step 6, according to the described first new echo-signal Z1With the second new echo-signal Z2Obtain with wave beam and difference beam, from
And difference and ratio are obtained, according to the difference and than obtaining the target off-axis angle θ of i-th of pulse with the mirror angular curve tablei, i=1,
2,...,K;By the target off-axis angle θ of i-th of pulseiWith the reference angle of i-th of pulseIt is added, obtains i-th of pulse to target
Measurement angle ψi, i=1,2 ... L.
Step 6 is specially:
(6a) is according to the described first new echo-signal Z1With the second new echo-signal Z2It obtains and wave beam PΣWith difference beam PΔ,
PΣ=Z1+Z2, PΔ=Z1-Z2;
(6b) is utilized and wave beam PΣWith difference beam PΔ, calculate difference and than P=imag (PΛ/PΣ), according to described poor and ratio
Value obtains the target off-axis angle θ of i-th of pulse by inquiry mirror angular curve tablei, i=1,2 ..., K;
(6c) is by the target off-axis angle θ of i-th of pulseiWith the reference angle of i-th of pulseIt is added, obtains i-th of pulse pair
The measurement angle of targetI=1,2 ... L.
Step 7, step 6 is repeated, until obtaining measurement angle of the L pulse to target, L measurement angle is made even
, the estimation angle of target is obtained.
Step 8, the first object off-axis angle estimation angle of the target being assigned in step 4;
Step 9, it is repeated in and executes step 5 to step 7, until the estimation angle of target meets preset error condition,
And by the estimation angle-determining obtained for the last time be target angle estimation value.
It is specially until the estimation angle of target meets preset error condition in step 9:The estimation angle and reality of target
The root-mean-square error of border angle is less than preset value, and wherein preset value is 10。
The effect of the present invention can be verified by following Computer Simulation:
One, simulated conditions
Simulated conditions 1:Assuming that the antenna surface that antenna is made of 10 array elements, is divided into two submatrixs, 5 battle arrays of each submatrix
The beam angle of member, antenna is about 10 °, and antenna carries out rotation sweep with the speed of 10s/r, and radar emits a frequency every 10ms
Rate is the pulse of 300MHz, and antenna receives 21 pulse signals every time, and antenna impulse center is not aligned with target, antenna impulse
4 pulses of center deviation target;Assuming that there are a target, the real angle of target is 10 °;Signal-to-noise ratio takes -10dB to arrive 10dB,
Here signal-to-noise ratio refers to single array element, single pulse;Iterations M=4.
Two, emulation content
Emulation 1, is generated respectively using simulated conditions 1 traditional single pulse angle-measuring method and difference beam curve and mirror angle song
Line, as shown in Figures 2 and 3 respectively, wherein Fig. 2 abscissas are angle, and ordinate is amplitude;Fig. 3 abscissas are angle, ordinate
For difference and ratio.
From Figure 2 it can be seen that when antenna centre normal alignment target, echo-signal and beam amplitude in maximum position, it is poor
Beam amplitude is not in the above situation if antenna does not have alignment target in minimum position.
As seen from Figure 3, when antenna centre normal alignment target and the difference of difference beam and than be 0, this with shown in Fig. 2
Situation match, compares figure 2 it can be found that in antennas the heart normal alignment target when, the difference beam amplitude of echo is 0, then
Its difference and ratio are necessarily 0;Similarly, if center of antenna discovery is not aligned with target, poor and ratio has deviation, in addition,
Mirror angular curve shown in Fig. 3 can be as the foundation of tabling look-up of the method for the present invention.
Traditional and difference beam single pulse method is respectively adopted to target using simulated conditions 1 and the method for the present invention is surveyed
Angle obtains the root-mean-square error curve that two methods change with signal-to-noise ratio, as shown in figure 4, abscissa is signal-to-noise ratio in Fig. 4, indulges
Coordinate is root-mean-square error.
From fig. 4, it can be seen that the method for the present invention is higher with difference beam single pulse method precision than traditional, especially in signal-to-noise ratio
When low, clearly, Fig. 4 fully demonstrates the validity of the method for the present invention to effect, in practical applications center of antenna normal pair
The case where inaccurate target, is relatively common, and the method for the present invention can reduce the influence of noise, improves angle measurement accuracy.
The above description is merely a specific embodiment, but scope of protection of the present invention is not limited thereto, any
Those familiar with the art in the technical scope disclosed by the present invention, can easily think of the change or the replacement, and should all contain
Lid is within protection scope of the present invention.Therefore, protection scope of the present invention should be based on the protection scope of the described claims.
Claims (6)
1. a kind of machine based on iterative processing sweeps metre wave radar Monopulse estimation method, which is characterized in that the method includes such as
Lower step:
Step 1, antenna that machine sweeps metre wave radar is set as the uniform line-array being made of N number of array element, and the uniform line-array is divided into
The right submatrix that the left submatrix and array number that array number is N/2 are N/2, wherein N is the even number more than or equal to 4;
Step 2, -90 ° of measurable angle range~90 ° that machine is swept to metre wave radar are divided into multiple wave positions, and is spaced between adjacent wave position
One preset angle;Determine the mirror angular curve of each wave position, composition mirror angular curve table;
Step 3, machine sweeps metre wave radar K pulse of transmitting, is spaced the second preset angle between adjacent pulse, determines i-th of pulse
Reference angleI=1,2 ..., K;Wherein, reference angle is the angle of center of antenna normal and horizontal reference plane;
Step 4, machine sweeps the echo-signal X that metre wave radar receives K pulse by antenna array, determines that the target of i-th of pulse is inclined
Shaft angle, i=1,2 ..., K, and to echo-signal X into horizontal phasing control, the echo-signal Y after phase adjustment is obtained, obtain phase
The corresponding pulse of maximum amplitude and the corresponding target off-axis angle of the pulse in echo-signal Y after the adjustment of position, are denoted as the first mesh
Mark off-axis angle;Angle of the target off-axis angle between target and center of antenna normal;
Step 5, the echo-signal T that left submatrix receives is obtained respectively1The echo-signal T received with right submatrix2, inclined in first object
(L-1)/2 angle value is taken on the left of shaft angle, takes (L-1)/2 angle value on the right side of first object off-axis angle, and adjacent angular it
Between between be divided into second preset angle;And the echo-signal T received from left submatrix1In obtain and first object off-axis angle on the left of
(L-1)/2 angle value corresponding range value composition first echo signal Z1, from the echo-signal T of right submatrix reception2In obtain
Range value composition second echo signal Z corresponding with (L-1)/2 angle value on the right side of first object off-axis angle2;L is wave beam side
The pulse number for including to Fig. 3 dB width;
Step 6, according to the first echo signal Z1With second echo signal Z2Obtain with wave beam and difference beam, to obtain
Difference and ratio, according to the difference and than obtaining the target off-axis angle θ of i-th of pulse with the mirror angular curve tablei, i=1,2 ...,
K;By the target off-axis angle θ of i-th of pulseiWith the reference angle of i-th of pulseIt is added, obtains measurement of i-th of pulse to target
Angle ψi, i=1,2 ... L;
Step 7, step 6 is repeated, until obtaining measurement angle of the L pulse to target, L measurement angle is averaged,
Obtain the estimation angle of target;
Step 8, the first object off-axis angle estimation angle of the target being assigned in step 4;
Step 9, it is repeated in and executes step 5 to step 8, until the estimation angle of target meets preset error condition, and general
The estimation angle-determining that last time obtains is the angle estimation value of target.
2. a kind of machine based on iterative processing according to claim 1 sweeps metre wave radar Monopulse estimation method, feature
It is, step 2 is specially:
Measurable angle range is divided into a wave position for -90 °~90 ° by (2a) at interval of the first preset angle △ θ, wherein the first preset angle
θ=50.8/ △ (N-d/ λ), d are array element interval, and λ is signal wavelength, λ/2 d=;
(2b) for any one wave position, obtain the wave position left submatrix synthesis wave beam BlWith the wave position right submatrix synthesis
Wave beam Br, according to the wave position left submatrix synthesis wave beam BlWith the wave position right submatrix synthesis wave beam BrObtain the wave position
With wave beam and difference beam;
(2c) obtains the difference and ratio of the wave position according to the wave position and wave beam and difference beam, and obtains the mirror angular curve of the wave position;
(2d) obtains the mirror angular curve of all wave positions, and the mirror angular curve of all wave positions is formed mirror angular curve table.
3. a kind of machine based on iterative processing according to claim 1 sweeps metre wave radar Monopulse estimation method, feature
It is, in step 3, machine sweeps metre wave radar and emits K pulse, and being spaced the second preset angle between adjacent pulse is specially:
Machine sweeps metre wave radar when antenna scanning works, and emits a pulse at interval of the second preset angle △ Φ, and in 3dB wave beams
Include L pulse in width;Wherein, △ Φ=360*t/T, t are the pulse repetition period, and T is the rotating speed that machine sweeps metre wave radar.
4. a kind of machine based on iterative processing according to claim 1 sweeps metre wave radar Monopulse estimation method, feature
It is, step 4 is specially:
(4a) machine sweeps the echo-signal X=AS+n that metre wave radar receives K pulse by antenna array, wherein A=[a (θ1),...a
(θi)...,a(θK)] indicate the target phase information that echo-signal includes, a (θi)=[1, exp (j2 π d/ λ sin θsi),...,
exp(j2π(N-1)d/λsinθi)]TFor the target off-axis angle vector of i-th of pulse, i=1,2 ..., K;D is array element interval, λ
For signal wavelength, []TIndicate the transposition of vector;
S=[S1,...Si...,SK]TIndicate the complex envelope of echo-signal, SiIndicate the complex envelope information of i-th of pulse,Wherein fdIndicate the Doppler frequency of target, fd=2Vf0/ c, V indicate that target sweeps metre wave radar with respect to machine
Radial velocity, f0Indicate that the centre frequency of radar emission signal, c represent the light velocity, t is the pulse repetition period;N represent mean value as
0, N × K rank white Gaussian noise matrixes that variance is 1;
(4b), into horizontal phasing control, obtains the echo-signal Y=a (θ after phase adjustment to echo-signal X0) X, wherein θ0For
The actual angle of target, a (θ0)=[1, exp (j2 π d/ λ sin (θ0)),…,exp(j2π(N-1)d/λsin(θ0))]T;
(4c) obtains the corresponding pulse of maximum amplitude and the corresponding first object of the pulse in the echo-signal Y after the adjustment of position
Off-axis angle θmax。
5. a kind of machine based on iterative processing according to claim 1 sweeps metre wave radar Monopulse estimation method, feature
It is, step 6 is specially:
(6a) is according to the first echo signal Z1With second echo signal Z2It obtains and wave beam PΣWith difference beam PΔ, PΣ=Z1+Z2,
PΔ=Z1-Z2;
(6b) is utilized and wave beam PΣWith difference beam PΔ, calculate difference and than P=imag (PΔ/PΣ), according to the value of the difference and ratio, lead to
It crosses inquiry mirror angular curve table and obtains the target off-axis angle θ of i-th of pulsei, i=1,2 ..., K;
(6c) is by the target off-axis angle θ of i-th of pulseiWith the reference angle of i-th of pulseIt is added, obtains i-th of pulse to target
Measurement angleI=1,2 ... L.
6. a kind of machine based on iterative processing according to claim 1 sweeps metre wave radar Monopulse estimation method, feature
It is, is specially until the estimation angle of target meets preset error condition in step 9:The estimation angle and actual corners of target
The root-mean-square error of degree is less than preset value, and wherein preset value is 1 °.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610504677.8A CN105974389B (en) | 2016-06-30 | 2016-06-30 | Machine based on iterative processing sweeps metre wave radar Monopulse estimation method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610504677.8A CN105974389B (en) | 2016-06-30 | 2016-06-30 | Machine based on iterative processing sweeps metre wave radar Monopulse estimation method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105974389A CN105974389A (en) | 2016-09-28 |
CN105974389B true CN105974389B (en) | 2018-09-25 |
Family
ID=56953605
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610504677.8A Active CN105974389B (en) | 2016-06-30 | 2016-06-30 | Machine based on iterative processing sweeps metre wave radar Monopulse estimation method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105974389B (en) |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108449098B (en) | 2018-03-22 | 2020-08-21 | 京东方科技集团股份有限公司 | Signal receiving method and device |
CN108594195B (en) * | 2018-04-20 | 2021-07-27 | 西安电子科技大学 | Low-repetition frequency modulation continuous wave radar aircraft target classification method based on division mixing |
CN109061594B (en) * | 2018-08-03 | 2022-10-28 | 中国航空工业集团公司雷华电子技术研究所 | Angle measuring method based on subarray space-time adaptive processing in strong clutter environment |
CN109521418B (en) * | 2018-12-28 | 2022-12-02 | 西安电子科技大学 | Foundation radar angle measurement method based on interference field |
CN110082750B (en) * | 2019-03-25 | 2021-05-07 | 西安电子科技大学 | Amplitude comparison angle measurement method capable of eliminating phase errors among channels |
CN110018439B (en) * | 2019-03-27 | 2020-09-25 | 中国电子科技集团公司第三十八研究所 | Direction finding method combining sum-difference beam and beam MUSIC |
CN110596692B (en) * | 2019-08-19 | 2022-10-14 | 电子科技大学 | Self-adaptive monopulse direction finding method based on joint constraint |
CN110988835B (en) * | 2019-11-27 | 2023-03-31 | 中国船舶重工集团公司第七二四研究所 | Distributed coherent radar angle measurement method |
CN112083382B (en) * | 2020-09-15 | 2022-08-02 | 四川九洲空管科技有限责任公司 | High-precision small phased array secondary radar azimuth compensation calculation method |
CN113189579B (en) * | 2021-03-24 | 2022-08-16 | 四川九洲空管科技有限责任公司 | Hybrid monopulse angle measurement system and method |
CN114675257B (en) * | 2022-05-30 | 2022-08-02 | 中国人民解放军国防科技大学 | Non-rectangular array polarized monopulse radar double-target resolution method |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7619556B1 (en) * | 2008-02-28 | 2009-11-17 | Rockwell Collins, Inc. | System and method for synthesizing localizer and glide slope deviations from weather radar |
CN101825707A (en) * | 2010-03-31 | 2010-09-08 | 北京航空航天大学 | Monopulse angular measurement method based on Keystone transformation and coherent integration |
CN101887120A (en) * | 2010-06-10 | 2010-11-17 | 哈尔滨工业大学 | Method for measuring target azimuth by single-beam mechanical scanning radar |
CN103389493A (en) * | 2013-06-25 | 2013-11-13 | 西安电子科技大学 | Multi-beam single-pulse angle measuring method based on beam selection method |
CN103728614A (en) * | 2014-01-15 | 2014-04-16 | 西安电子科技大学 | Mechanical scanning meter wave radar based method for improving single pulse angle measurement |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP5093298B2 (en) * | 2010-06-04 | 2012-12-12 | 株式会社デンソー | Direction detection device |
-
2016
- 2016-06-30 CN CN201610504677.8A patent/CN105974389B/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7619556B1 (en) * | 2008-02-28 | 2009-11-17 | Rockwell Collins, Inc. | System and method for synthesizing localizer and glide slope deviations from weather radar |
CN101825707A (en) * | 2010-03-31 | 2010-09-08 | 北京航空航天大学 | Monopulse angular measurement method based on Keystone transformation and coherent integration |
CN101887120A (en) * | 2010-06-10 | 2010-11-17 | 哈尔滨工业大学 | Method for measuring target azimuth by single-beam mechanical scanning radar |
CN103389493A (en) * | 2013-06-25 | 2013-11-13 | 西安电子科技大学 | Multi-beam single-pulse angle measuring method based on beam selection method |
CN103728614A (en) * | 2014-01-15 | 2014-04-16 | 西安电子科技大学 | Mechanical scanning meter wave radar based method for improving single pulse angle measurement |
Non-Patent Citations (4)
Title |
---|
Adaptive Monopulse Angle Measurement for Meter-Wave Radar Based on Differential Constraints Technique;ZY Jun et.al;《Advanced Materials Research》;20141231;第812-817页 * |
一种多目标情况下的单脉冲测角方法;赵永波等;《西安电子科技大学学报》;20050630;第32卷(第3期);第383-386页 * |
一种改进的单脉冲测角方法;赵永波等;《第九届全国雷达学术年会论文集》;20040820;第513-516页 * |
基于鉴角曲线拟合的自适应单脉冲测角方法;陈成增等;《系统工程与电子技术》;20130731;第35卷(第7期);第1403-1408页 * |
Also Published As
Publication number | Publication date |
---|---|
CN105974389A (en) | 2016-09-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105974389B (en) | Machine based on iterative processing sweeps metre wave radar Monopulse estimation method | |
CN103728614B (en) | The improvement Monopulse estimation method of metre wave radar is swept based on machine | |
CN105974390B (en) | Ji Sao metre wave radars mass centre angle-measuring method based on doppler information | |
CN107976660B (en) | Missile-borne multi-channel radar ultra-low-altitude target analysis and multi-path echo modeling method | |
CN108549059B (en) | Low-altitude target elevation angle estimation method under complex terrain condition | |
CN106066473B (en) | The multi-beam joint angle measurement of MIMO radar target and Contact fusion method under orthogonal waveforms | |
CN103744077B (en) | The angle-measuring method of metre wave radar under multi-target condition swept by machine | |
CN105445701B (en) | The pulse angle estimating method of DDMA MIMO radar targets | |
CN104155647B (en) | The method sweeping the estimation azimuth of target of meter wave array radar based on machine | |
CN101887120B (en) | Method for measuring target azimuth by single-beam mechanical scanning radar | |
CN104931948B (en) | A kind of improved method of reception scheme of FDA radars based on conventional beams scanning | |
CN104237844B (en) | Distributed meter wave planar array radar azimuth angle measuring method based on phase compensation | |
CN103197294B (en) | Elevation angle estimating method of multi-frequency fusion maximum likelihood low-altitude target | |
CN110488255A (en) | A kind of phased-array radar pulse high-resolution angle measuring system and method | |
CN108845325A (en) | Towed linear-array sonar submatrix error misfits estimation method | |
CN112147593B (en) | Four-dimensional parameter estimation method for high-speed dense explosive fragment target | |
CN110879391B (en) | Radar image data set manufacturing method based on electromagnetic simulation and missile-borne echo simulation | |
CN105824016B (en) | The steady space-time adaptive processing method of motion platform detections of radar treetop level target | |
CN109521418B (en) | Foundation radar angle measurement method based on interference field | |
CN104122548B (en) | Meter wave array radar angle-measuring method swept by machine | |
CN106125058B (en) | The improvement mass centre angle-measuring method of metre wave radar is swept based on machine | |
CN111708023B (en) | Phase inconsistency correction method and imaging system for millimeter wave multi-transmitting multi-receiving linear array | |
CN104199005B (en) | The distributed meter wave array radar optimum length of base determines method | |
CN112068087B (en) | Shore-based multi-channel radar aircraft-simulated sea clutter measurement method | |
CN111175705B (en) | Radar communication integrated design method based on spatial waveform two-dimensional weighting |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |