CN104714226B - Tracking before a kind of Dynamic Programming detection based on phase - Google Patents

Tracking before a kind of Dynamic Programming detection based on phase Download PDF

Info

Publication number
CN104714226B
CN104714226B CN201510133826.XA CN201510133826A CN104714226B CN 104714226 B CN104714226 B CN 104714226B CN 201510133826 A CN201510133826 A CN 201510133826A CN 104714226 B CN104714226 B CN 104714226B
Authority
CN
China
Prior art keywords
frame
target
resolution cell
targetpath
flight path
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
Application number
CN201510133826.XA
Other languages
Chinese (zh)
Other versions
CN104714226A (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 CN201510133826.XA priority Critical patent/CN104714226B/en
Publication of CN104714226A publication Critical patent/CN104714226A/en
Application granted granted Critical
Publication of CN104714226B publication Critical patent/CN104714226B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Systems 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/66Radar-tracking systems; Analogous systems

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 provides tracking before a kind of Dynamic Programming detection based on phase, belong to Radar Targets'Detection tracking technique field, more particularly to fluctuating target detecting and tracking technical field under multiple Gauss noise background.After radar return data are received, multiple likelihood ratio is calculated using phase information, all possible flight path is scanned for by the multiple likelihood ratio of accumulation, dbjective state is estimated.This method by using phase information, can preferably embody target and the difference of noise, improve the detecting and tracking performance to fluctuating target than track algorithm before existing Dynamic Programming detection.

Description

Tracking before a kind of Dynamic Programming detection based on phase
Technical field
The present invention relates to Radar Targets'Detection tracking technique, fluctuating weak target radar is examined more particularly under multiple Gauss background Survey tracking technique.
Background technology
With reaching its maturity for stealth technology, the radar cross section of the target such as aircraft, guided missile reduces one to two quantity Level, while including helicopter, SUAV, a large amount of Military applications for the Small objects at a slow speed of low-latitude flying such as device that fly at low speed, So that all suffering to dim targets detection is difficult asking in each field such as airborne fire control, airborne early warning, ground information, shipborne radar Topic.Therefore, the detection and tracking to weak target turn into the technical problem that radar is badly in need of solving.
Tracking is a kind of effective dim target detection tracking before Dynamic Programming detection, and it is returned to multiframe radar Wave number is searched for, so as to estimate target according to Combined Treatment is carried out by accumulating targetpath value function to all possible flight path Time of day.Tracking targetpath value function has two kinds before existing Dynamic Programming detection:(1) amplitude;(2) envelope logarithm is seemingly So compare.Both targetpath value functions do not utilize phase information, and not utilizing for phase information causes target detection tracing property It can decline, be unfavorable for the detecting and tracking of weak target.To improve the detecting and tracking performance to constant amplitude weak target, Davey Et al. propose tracking before the detection based on multiple likelihood ratio, referring to document " S.~J.Davey, M.~G.Rutten, and B.~Cheung, ``Using phase to improve track-before-detect, " IEEE Trans.Aerosp.Electron.Syst., vol.48, no.1, pp.832-849, Jan.2012 ", as a result show phase information The detecting and tracking performance to constant amplitude target can effectively be improved.
But the radar cross section of real goal can not possibly use a simple constant.Generally even for simple mesh Mark, radar cross section is also the angle of sight, frequency, the complicated function of plan.The fluctuating target models of Swerling 3 are conventional fluctuating Object module, it is 4 ° of χ that it, which characterizes radar cross section probability density function,2The target of distribution, is that many small scattering objects and one are strong The radar cross section of the approximate solution of scattering object, wherein strong scatterer is equal to small scattering object radar cross section sumTimes.
The content of the invention
The purpose of the present invention is in the tracking fluctuating targets of Swerling 3 for tracking before existing Dynamic Programming detection When the detecting and tracking probability that exists it is low the problem of, tracking before a kind of Dynamic Programming detection based on phase of Curve guide impeller, The detecting and tracking performance and the purpose of the stability of a system to the fluctuating weak targets of Swerling 3 are improved using phase information.
The invention provides tracking before a kind of Dynamic Programming detection based on phase, this method includes:
Step 1:Initialization systematic parameter includes:Datum plane size Nx×Ny, frame number K, state transfer number q are handled, is utilized The thresholding V that Monte Carlo simulation experimental calculation goes outT, target point spread function h (xk), noise average power σn, target mean power σt
Step 2:Resolution cell each to each frame carries out targetpath value function assignment:Each each resolution cell pair of frame The targetpath value function answered is the corresponding multiple likelihood ratio of the resolution cell;
Step 3, by dynamic programming algorithm all possible targetpath in preceding K frames is scanned for, and record each mesh Positional information of the flight path in the corresponding resolution cell of every frame is marked, the resolution cell correspondence of same targetpath will be belonged in preceding K frames Targetpath value function be overlapped, obtain target discrimination value of the target in k-th frame;
Step 4, by each target discrimination value of k-th frame with setting thresholding VTIt is compared, if higher than thresholding, assert should With the presence of target in decision content correspondence resolution cell;If all target discrimination values of k-th frame are below thresholding, announce that target is not deposited ;
If step 5, step 4 assertive goal are present, K frame location informations before the target recorded using step 3 recover mesh Target flight path;
Step 6:False track is deleted:All targetpaths for recovering to obtain in step 5 are contrasted, if a plurality of Flight path has same position in a certain frame, then the flight path with maximum target decision content in these flight paths navigates as real goal Mark, deletes remaining flight path;
Step 7, output targetpath.
The step 2 is concretely comprised the following steps:
The corresponding multiple likelihood ratio of step 2.1, each each resolution cell of frame of calculating:To resolution cell (i, j), it is assumed that wherein Comprising target, then dbjective state is xk=[i, vx,j,vy], wherein the position of target is (i, j), and the speed of target is (vx,vy), Then the corresponding multiple likelihood ratio of the resolution cell is
Wherein
zk(i, j) represents that resolution cell (i, j) is corresponding
Measuring value;Zk=[zk(1,1),…,zk(1,Ny),zk(2,1),…,zk(1,Ny)]TRepresent the collection of kth frame measuring value Close;(·)TRepresent transposition;Represent target point Spread function, wherein h(m,n)(xk);1≤m≤Nx,1≤n≤NyRepresent contribution of the target energy to resolution cell (m, n);NxTable Show the line number of datum plane, NyRepresent the columns of datum plane, hH(xk) represent h (xk) conjugate transposition, M-1.5,0Represent parameter For -1.5,0 Whittaker functions.
Step 2.2, the corresponding targetpath value function of the individual resolution cell of kth frame (i, j) is entered as step 2.1 fallen into a trap The obtained corresponding multiple likelihood ratio of the resolution cell.
The step 3 is concretely comprised the following steps:
Step 3.1, to each Range resolution unit object decision content assignment of the 1st frame:1st frame each Range resolution unit correspondence Target discrimination value obtained targetpath value function is calculated in step 2.3 for the Range resolution unit;
Step 3.2, to each Range resolution unit object decision content assignment of the 2nd frame:It is assumed that being deposited in the 2nd each resolution cell of frame In a target, judge the 2nd each target of frame region that may be present in the 1st frame, search out target discrimination value in the region Maximum resolution cell and the position for recording the unit, the target discrimination value of the corresponding resolution cell of the 2nd frame are to find out in the 1st frame Maximum target decision content targetpath value function corresponding with the 2nd frame resolution cell and;
Step 3.3, using with step 3.2 identical method the 3rd frame is calculated to the target discrimination value of k-th frame.
K value is 3~20 in the step 3.
The step 6 is concretely comprised the following steps:
In step 6.1, all targetpaths for judging no process, there is the mesh of maximum target decision content in k-th frame Mark flight path is determined as real goal flight path;
Step 6.2, by each without through judgement targetpath and real goal flight path be compared, if the boat Mark has same position with real goal flight path in certain frame, then the flight path is determined as into false track;
If step 6.3, the targetpath still judged either with or without process, repeat step 6.1, step 6.2 are until all mesh Flight path is marked all by judging.
Tracking before a kind of Dynamic Programming detection based on phase of the present invention, calculates each frame using phase information and respectively differentiates The corresponding multiple likelihood ratio of unit, chooses multiple likelihood ratio as targetpath value function and carries out Dynamic Programming value accumulation.Utilize phase Information can effectively improve the detecting and tracking performance to the fluctuating targets of Swerling 3, reduce missing inspection.
Brief description of the drawings
Fig. 1 is FB(flow block) of the invention.
Fig. 2 is the present invention (Nx=30, Nx=30, q=9, K=5, σn=1, if m=i-1, i, i+1 then h(m,n)(xk)= exp(-(m-i)2/ 8), otherwise h(m,n)(xk)=0) and traditional single frame detection (not using phase, Nx=30, Nx=30, q=9, K =1, σn=1, if m=i-1, i, i+1 then h(m,n)(xk)=exp (- (m-i)2/ 8), otherwise h(m,n)(xkDetect performance in)=0) Comparison diagram.The solid line of wherein band " o " represents the performance of the present invention, and the solid line of band " " represents traditional single frame detection performance.
Fig. 3 is the present invention (Nx=30, Nx=30, q=9, K=5, σn=1, if m=i-1, i, i+1 then h(m,n)(xk)= exp(-(m-i)2/ 8), otherwise h(m,n)(xkTracking before the planning detection of)=0) and conventional dynamic (not using phase, Nx=30, Nx=30, q=9, K=5, σn=1, if m=i-1, i, i+1 then h(m,n)(xk)=exp (- (m-i)2/ 8), otherwise h(m,n)(xk) =0) detect performance comparison figure.The solid line of wherein band " o " represents the performance of the present invention, and the solid line of band " " represents that conventional dynamic is advised Draw tracking detection performance before detection.
Embodiment
The present invention chooses multiple likelihood ratio as targetpath value function, and the selection of the targetpath value function takes full advantage of Target amplitude statistical property, noise statisticses and target phase information.Tracking is used before conventional dynamic planning detection Three kinds of information of amplitude target flight path value function are not all utilized, and envelope log-likelihood ratio make use of target amplitude statistical property with Noise statisticses, phase information is not utilized but.The present invention can further improve Dynamic Programming detection using phase information The detecting and tracking performance of preceding tracking.
The main method for using l-G simulation test of the invention is verified that all steps, conclusion are all on MATLAB R2012b Checking is correct.Specific implementation step is as follows:
Step 1:Initialization systematic parameter includes:(the i.e. N of datum plane size 30 × 30x=Ny=30), processing frame number K= 5, state transfer number q=9, the thresholding V gone out using Monte Carlo simulation experimental calculationT=[- 0.0124,0.6802,1.2511, 1.7538,1.9351,2.3633,1.9839,1.6518,1.4338,0.2654, -0.5999, -1.8316], target point diffusion Function h (xk):If m=i-1, i, i+1 then h(m,n)(xk)=exp (- (m-i)2/ 8), otherwise h(m,n)(xk)=0, noise is average Power σn=1, target mean power σt=[1.0000,1.2589,1.5849,1.9953,2.5119,3.1623,3.9811, 5.0119,6.3096,7.9433,10.0000,12.5893] (correspondence signal to noise ratio is 0-11dB).
Step 2:Resolution cell each to each frame carries out targetpath value function assignment:Each each resolution cell pair of frame The targetpath value function answered is the corresponding multiple likelihood ratio of the resolution cell;
The corresponding multiple likelihood ratio of step 2.1, each each resolution cell of frame of calculating:To resolution cell (i, j), it is assumed that wherein Comprising target, (then dbjective state is xk=[i, vx,j,vy], wherein the position of target is (i, j), and the speed of target is (vx, vy)), then the corresponding multiple likelihood ratio of the resolution cell is
Wherein
zk(i, j) represents the corresponding measuring value of resolution cell (i, j);Zk=[zk(1,1),…,zk(1,Ny), zk(2,1),…,zk(1,Ny)]TRepresent the set of kth frame measuring value;(·)TRepresent transposition;Represent target point spread function, wherein h(m,n) (xk);1≤m≤Nx,1≤n≤NyRepresent contribution of the target energy to resolution cell (m, n);NxRepresent the line number of datum plane, Ny Represent the columns of datum plane, hH(xk) represent h (xk) conjugate transposition, M-1.5,0Represent the Whittaker that parameter is -1.5,0 Function.
Step 2.2, the corresponding targetpath value function of the individual resolution cell of kth frame (i, j) is entered as step 2.1 fallen into a trap The obtained corresponding multiple likelihood ratio of the resolution cell.
Step 3, by dynamic programming algorithm all possible targetpath in preceding K frames is scanned for, and record each mesh Positional information of the flight path in the corresponding resolution cell of every frame is marked, the resolution cell correspondence of same targetpath will be belonged in preceding K frames Targetpath value function be overlapped, obtain target discrimination value of the target in k-th frame;
Step 3.1, to each Range resolution unit object decision content assignment of the 1st frame:1st frame each Range resolution unit correspondence Target discrimination value obtained targetpath value function is calculated in step 2.3 for the Range resolution unit;
Step 3.2, to each Range resolution unit object decision content assignment of the 2nd frame:It is assumed that being deposited in the 2nd each resolution cell of frame In a target, judge the 2nd each target of frame region that may be present in the 1st frame, search out target discrimination value in the region Maximum resolution cell and the position for recording the unit, the target discrimination value of the corresponding resolution cell of the 2nd frame are to find out in the 1st frame Maximum target decision content targetpath value function corresponding with the 2nd frame resolution cell and;Specific processing mode is as follows:2nd The individual Range resolution unit of frame (i, j), it is assumed that wherein in the presence of a target, then the target may be deposited in former frame (the 1st frame) Region be (i-1-vx,j-1-vy),(i-1-vx,j-vy),(i-1-vx,j+1-vy),(i-vx,j-1-vy),(i-vx,j- vy), (i-vx,j+1-vy),(i+1-vx,j-1-vy),(i+1-vx,j-vy),(i+1-vx,j+1-vy);
Step 3.3, using with step 3.2 identical method the 3rd frame is calculated to the target discrimination value of k-th frame.
Step 4, by each target discrimination value of k-th frame with setting thresholding VTIt is compared, if higher than thresholding, assert should With the presence of target in decision content correspondence resolution cell;If all target discrimination values of k-th frame are below thresholding, announce that target is not deposited ;
If step 5, step 4 assertive goal are present, K frame location informations before the target recorded using step 3 recover mesh Target flight path;
Step 6:False track is deleted:All targetpaths for recovering to obtain in step 5 are contrasted, if a plurality of Flight path has same position in a certain frame, then the flight path with maximum target decision content in these flight paths navigates as real goal Mark, deletes remaining flight path;
In step 6.1, all targetpaths for judging no process, there is the mesh of maximum target decision content in k-th frame Mark flight path is determined as real goal flight path;
Step 6.2, by each without through judgement targetpath and real goal flight path be compared, if the boat Mark has same position with real goal flight path in certain frame, then the flight path is determined as into false track;
If step 6.3, the targetpath still judged either with or without process, repeat step 6.1, step 6.2 are until all mesh Flight path is marked all by judging.
Step 7, output targetpath.
It can be seen that method proposed by the present invention is than tracking and traditional single frame detection before conventional dynamic planning detection Method detection performance is respectively provided with very big lifting, effectively can carry out detecting and tracking to the faint fluctuating targets of Swerling 3.

Claims (4)

1. a kind of tracking before Dynamic Programming detection based on phase, this method includes:
Step 1:Initialization systematic parameter includes:Datum plane size Nx×Ny, frame number K, state transfer number q are handled, illiteracy spy is utilized The thresholding V that Carlow emulation experiment is calculatedT, target point spread function h (xk), noise average power σn, target mean power σt
Step 2:Resolution cell each to each frame carries out targetpath value function assignment:Each each resolution cell of frame is corresponding Targetpath value function is the corresponding multiple likelihood ratio of the resolution cell;
Step 3, by dynamic programming algorithm all possible targetpath in preceding K frames is scanned for, and record each target boat Mark will belong to the corresponding mesh of resolution cell of same targetpath in the positional information of the corresponding resolution cell of every frame in preceding K frames Mark flight path value function is overlapped, and obtains target discrimination value of the target in k-th frame;
Step 4, by each target discrimination value of k-th frame with setting thresholding VTIt is compared, if higher than thresholding, assert the judgement With the presence of target in value correspondence resolution cell;If all target discrimination values of k-th frame are below thresholding, announce that target is not present;
If step 5, step 4 assertive goal are present, K frame location informations before the target recorded using step 3 recover target Flight path;
Step 6:False track is deleted:All targetpaths for recovering to obtain in step 5 are contrasted, if a plurality of flight path There is same position in a certain frame, then using the flight path with maximum target decision content in these flight paths as real goal flight path, Delete remaining flight path;
Step 7, output targetpath;
It is characterized in that the step 2 is concretely comprised the following steps:
The corresponding multiple likelihood ratio of step 2.1, each each resolution cell of frame of calculating:To resolution cell (i, j), it is assumed that wherein include Target, then dbjective state is xk=[i, vx,j,vy], wherein the speed of target is (vx,vy), then the resolution cell is corresponding answers seemingly So ratio is
Wherein
zk(i, j) represents the corresponding measuring value of resolution cell (i, j);Zk=[zk(1,1),...,zk(1,Ny),zk(2,1),...,zk(1, Ny)]TRepresent the set of kth frame measuring value;(·)TRepresent transposition; Represent target point spread function, wherein h(m,n)(xk);1≤m≤Nx,1≤n≤NyRepresent target energy to resolution cell (m, n) Contribution;NxRepresent the line number of datum plane, NyRepresent the columns of datum plane, hH(xk) represent h (xk) conjugate transposition, M-1.5,0 Represent the Whittaker functions that parameter is -1.5,0;
Step 2.2, the corresponding targetpath value function of the individual resolution cell of kth frame (i, j) is entered as calculating in step 2.1 The corresponding multiple likelihood ratio of the resolution cell arrived.
2. tracking before a kind of Dynamic Programming detection based on phase as claimed in claim 1, it is characterised in that the step Rapid 3 concretely comprise the following steps:
Step 3.1, to each Range resolution unit object decision content assignment of the 1st frame:Each corresponding mesh of Range resolution unit of 1st frame Mark decision content is that the Range resolution unit calculates obtained targetpath value function in step 2.3;
Step 3.2, to each Range resolution unit object decision content assignment of the 2nd frame:It is assumed that having one in the 2nd each resolution cell of frame Individual target, judges the 2nd each target of frame region that may be present in the 1st frame, searches out target discrimination value in the region maximum Resolution cell and record the position of the unit, the target discrimination value of the corresponding resolution cell of the 2nd frame is the maximum found out in the 1st frame Target discrimination value targetpath value function corresponding with the 2nd frame resolution cell and;
Step 3.3, using with step 3.2 identical method the 3rd frame is calculated to the target discrimination value of k-th frame.
3. tracking before a kind of Dynamic Programming detection based on phase as claimed in claim 2, it is characterised in that the step K value is 3~20 in rapid 3.
4. tracking before a kind of Dynamic Programming detection based on phase as claimed in claim 1, it is characterised in that the step Rapid 6 concretely comprise the following steps:
In step 6.1, all targetpaths for judging no process, the target for having maximum target decision content in k-th frame is navigated Mark is determined as real goal flight path;
Step 6.2, each is compared without targetpath and the real goal flight path by judging, if the flight path and Real goal flight path has same position in certain frame, then the flight path is determined as into false track;
If step 6.3, the targetpath still judged either with or without process, repeat step 6.1, step 6.2 are navigated until all targets Mark is all by judging.
CN201510133826.XA 2015-03-25 2015-03-25 Tracking before a kind of Dynamic Programming detection based on phase Active CN104714226B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510133826.XA CN104714226B (en) 2015-03-25 2015-03-25 Tracking before a kind of Dynamic Programming detection based on phase

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510133826.XA CN104714226B (en) 2015-03-25 2015-03-25 Tracking before a kind of Dynamic Programming detection based on phase

Publications (2)

Publication Number Publication Date
CN104714226A CN104714226A (en) 2015-06-17
CN104714226B true CN104714226B (en) 2017-08-11

Family

ID=53413712

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510133826.XA Active CN104714226B (en) 2015-03-25 2015-03-25 Tracking before a kind of Dynamic Programming detection based on phase

Country Status (1)

Country Link
CN (1) CN104714226B (en)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109523482B (en) * 2018-11-14 2021-04-30 太原理工大学 Deep neural network-based restoration method for degraded image containing texture
CN109669180B (en) * 2018-12-14 2020-09-15 河海大学 Continuous wave radar unmanned aerial vehicle detection method

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4214239A (en) * 1963-04-02 1980-07-22 The United States Of America As Represented By The Secretary Of The Navy Magnetic drum storage for the track before detect system
CN102901956A (en) * 2012-09-27 2013-01-30 电子科技大学 Method for detecting weak target by radar
CN104007422A (en) * 2014-05-21 2014-08-27 西安电子科技大学 Complex likelihood ratio track-before-detect method based on dynamic planning
CN104062651A (en) * 2014-06-30 2014-09-24 电子科技大学 Track-before-detect method based on G0 clutter background and constant target amplitude
CN104076355A (en) * 2014-07-04 2014-10-01 西安电子科技大学 Method for conducting before-detection tracking on weak and small target in strong-clutter environment based on dynamic planning

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4214239A (en) * 1963-04-02 1980-07-22 The United States Of America As Represented By The Secretary Of The Navy Magnetic drum storage for the track before detect system
CN102901956A (en) * 2012-09-27 2013-01-30 电子科技大学 Method for detecting weak target by radar
CN104007422A (en) * 2014-05-21 2014-08-27 西安电子科技大学 Complex likelihood ratio track-before-detect method based on dynamic planning
CN104062651A (en) * 2014-06-30 2014-09-24 电子科技大学 Track-before-detect method based on G0 clutter background and constant target amplitude
CN104076355A (en) * 2014-07-04 2014-10-01 西安电子科技大学 Method for conducting before-detection tracking on weak and small target in strong-clutter environment based on dynamic planning

Also Published As

Publication number Publication date
CN104714226A (en) 2015-06-17

Similar Documents

Publication Publication Date Title
Uhlmann Algorithms for multiple-target tracking
CN104714225B (en) Dynamic programming tracking-before-detection method based on generalized likelihood ratios
CN105277923B (en) A kind of single channel Radar Signal Sorting Method
CN108280395B (en) Efficient identification method for flight control signals of low-small-slow unmanned aerial vehicle
CN103728599B (en) The method of false targets interference is suppressed with the active radar and passive radar net of other place configure
CN111079859B (en) Passive multi-station multi-target direction finding cross positioning and false point removing method
CN104076355A (en) Method for conducting before-detection tracking on weak and small target in strong-clutter environment based on dynamic planning
CN106383339B (en) A kind of mirror target suppressing method of multi-site radar signal grade joint-detection
US10317574B2 (en) System and method for identification and/or nowcasting of hail events
CN109975780A (en) Helicopter model recognizer based on pulse Doppler radar time domain echo
CN106842194A (en) A kind of adaptive targets detection method and device
CN103728598A (en) Method for restraining flight path deception jamming through remotely configured initiative and passive radar net
CN106934324A (en) Based on the radar data correlating methods for simplifying many hypothesis algorithms
CN103901408A (en) CFAR detection method according to uniformly divided reference signals and radar system
Schultz et al. Automated storm tracking and the lightning jump algorithm using GOES-R Geostationary Lightning Mapper (GLM) proxy data
CN104714226B (en) Tracking before a kind of Dynamic Programming detection based on phase
CN105259540B (en) A kind of optimization method of multistation radar anti-active cheating formula interference
CN105445699A (en) Range measurement method and range measurement system based on non line-of-sight error elimination
CN104407345A (en) Improved dynamic planning method for weak target detection tracking
CN113484838B (en) Multi-base radar active false target identification method and system
CN104035084A (en) Dynamic planning pre-detection tracking method for heterogeneous clutter background
CN104880696B (en) One kind is based on the preferable multi-channel target detection method of passage
CN105259554B (en) A kind of method and device of multiple target classification tracking
CN106093891A (en) Radar fence anti-crowding measure false target jamming profile method based on doppler velocity inspection
CN106291472B (en) A kind of disturbance restraining method of external illuminators-based radar

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