CN108957503A - Positioning and false target recognition methods based on ADS-B message step-by-step counting - Google Patents
Positioning and false target recognition methods based on ADS-B message step-by-step counting Download PDFInfo
- Publication number
- CN108957503A CN108957503A CN201810586091.XA CN201810586091A CN108957503A CN 108957503 A CN108957503 A CN 108957503A CN 201810586091 A CN201810586091 A CN 201810586091A CN 108957503 A CN108957503 A CN 108957503A
- Authority
- CN
- China
- Prior art keywords
- base station
- ads
- pulse
- message
- positioning
- 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
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/38—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
- G01S19/39—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/42—Determining position
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)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
The invention discloses a kind of positioning and false target recognition methods based on ADS-B message step-by-step counting, the counting for receiving ADS-B message pulse signal and its preamble pulse of the base station received from the aircraft transmitting with S mode answering machine is detectd in the invention by difference, estimates each step-by-step counting deviation for detecing and receiving and receiving between base station and refer to aircraft ADS-B;Then the ADS-B signal or secondary radar answer signal of unknown other aircrafts transmitting of the location information that receives of base station are received for multiple detect, false target can be effectively verified using the range estimation method for containing space.The present invention can be realized the MDS system based on ADS-B, without receiving synchronization time between base station in each detect, the time synchronization cost of multipoint location system is reduced, the deficiency of secondary radar and ADS-B positioning is effectively made up, false target information is excluded, ensures civil aviaton's flight safety.
Description
Technical field
The present invention relates to a kind of positioning and false target recognition methods based on ADS-B message step-by-step counting.
Background technique
With the fast development of civil aviaton's industry, the quantity of aircraft significantly increases, in order to preferably ensure the safety of flight
Property, ADS-B technology are come into being, and compensate for the deficiency of secondary radar to a certain extent.But due to technical restriction, ADS-B
It is not loaded in above all aircraft, while ICAO is easy to be forged, such ADS-B receiver is difficult to discriminate whether
The authenticity of target.At the same time, multipoint location system now is spread between receiving terminal to solve receiving station's clocking error
If wire link, its cost is caused to be up to thousands of ten thousand yuan, local airport can not all purchase the set facility.
Summary of the invention
To solve the deficiencies in the prior art, the purpose of the present invention is to provide a kind of based on ADS-B message step-by-step counting
Positioning and false target recognition methods solve this problem using model and programming is established, greatly reduce cost, fill
Divide and existing technological means is utilized to make up the defect of ADS-B technology and secondary radar, enhances airport scene monitoring
Problem has greatly ensured the safety of aircraft, has identified false target in time.
In order to achieve the above objectives, the present invention adopts the following technical scheme that:
A kind of positioning and false target recognition methods based on ADS-B message step-by-step counting, characterized in that including walking as follows
It is rapid:
The aircraft that step 1) carries S mode answering machine sends data pulse, respectively detects and receives base station calculating preamble pulse starting
Number;
Step 2) is each to detect and receives base station and match identical pulse signal, is counted by pulse deviation, calculates and each detects receipts
The measurement error of base station;
Step 3) introduces measurement pulse deviation correction mechanism accurately to obtain the measurement error of each base station;
Step 4) given threshold is compared the base station number of the ADS-B signal received with threshold value, if being not less than threshold
Value, then go to step 5), otherwise, go to step 6);
Signal station according to star-like arrangement, while to objective emission signal, is considered the measurement error of each base station by step 5),
Intersection point, which is generated, by three unilateral curved surfaces determines position;
Step 6) receives base station and obtains target real time position or divide flight range by detecing, and is consistent verification method using purpose
Exclude spurious signal interference.
A kind of positioning and false target recognition methods based on ADS-B message step-by-step counting above-mentioned, characterized in that described
The particular content of step 1) are as follows:
The answer signal of S mode answering machine is carried as signal radiation source using aircraft, ground, which is respectively detectd, receives base station from receipts
Pulse number before top guide starts to calculate before to message;
The preamble pulse Base Serial Number received is calculated, is matched to identical pulse signal, and send it to main website, main website
It is uniformly processed, ADS-B message should carry out pulse transmission since first sends pulse with the pulse distance of 0.5 microsecond, when detecing
After receipts base station receives same message preamble pulse, the step-by-step counting received before is calculated.
A kind of positioning and false target recognition methods based on ADS-B message step-by-step counting above-mentioned, characterized in that described
The particular content of step 2) are as follows:
Step-by-step counting before when receiving station i receives preface information is denoted as Pi, the time of echo signal transmitting is consistent
, start to calculate to receive the same inceptive impulse signal, calculates the intermediate propagation by road and receiving station connects
Pulse signal error is received, when receiving station begins search for ADS-B preamble pulse, pulse signal received before is collected
Get up;
Wherein, the error for receiving base station and center base station is each detectdΔPi1It is each detect
The pulse number for receiving pulse and center base station that base station receives subtracts each other,It is Average pulse width, Δ l is S mode answering machine hair
Signal out carries latitude and longitude information, the aircraft obtained by calculation with S mode answering machine to each receiving station away from
From;C is the speed of pulse transmitting.
A kind of positioning and false target recognition methods based on ADS-B message step-by-step counting above-mentioned, characterized in that described
The particular content of step 3) are as follows:
Using the method for weighted moving average come update the system error, specific formula is as follows:Wherein, Yi
It is step-by-step counting difference between base station that i-th observes, XiThe weight of i-th of the observation artificially assigned, according to it is pre-
The distance of the observed value of survey phase gives flexible strategy value, and closer observed value gives smaller flexible strategy value, the sum of all flexible strategy values
Sum to 1;The sum of observation is n;
Definition set { Y1, Y2, Y3... ..., Y100, wherein by newest observation Δ Pi1It is defined as Y1, by Y100
Set is eliminated, the lower target number of other all values all increases by one;
It gives different observations different flexible strategy, gives larger flexible strategy value for the observed value close to time span of forecast, for
Range prediction phase farther away observed value then accordingly gives lesser flexible strategy value.
A kind of positioning and false target recognition methods based on ADS-B message step-by-step counting above-mentioned, characterized in that definition
Flexible strategy are e-i, Y is all used by calculating the value that reality obtains every time101It indicates, It new every time will obtain
The Y arrived101It feeds back in main website, as Δ Pi1。
A kind of positioning and false target recognition methods based on ADS-B message step-by-step counting above-mentioned, characterized in that described
4 are set the threshold in step 4).
A kind of positioning and false target recognition methods based on ADS-B message step-by-step counting above-mentioned, characterized in that described
The particular content of step 5) are as follows:
The real time position of other airbound targets is determined using the method for multipoint positioning, is at least needed four to detect and is received base station,
Three unilateral hyperboloids are formed to generate intersection point;Model is established, the measurement error for each receiving station that step 3) obtains is included in
Model calculates, and calculates target position by the method for multipoint positioning.
A kind of positioning and false target recognition methods based on ADS-B message step-by-step counting above-mentioned, characterized in that including
Following steps:
Step 51) is since the multipoint positioning in three-dimensional space at least needs four base stations, by receiving station according to star cloth station:
Main website occupies center, chooses other three lines between secondary station and main website and forms a star, main website position coordinates are set as
(x1, y1, z1), it is assumed that spatial position coordinate is (x, y, z);RiIt is to detect the distance for receiving base station to echo signal for i-th, i=1,2,
3,4, R1iIt indicates that echo signal arrival the distance between main website and i-th of secondary station are poor, is indicated with equation: t1iEqual to each time for detecing and receiving base station and receiving the unknown aircraft of geography information
Subtract each measurement error detectd and receive base station and main website, i.e. t1i=T1i-Δε1i, wherein T1iIt is same that be the base station i receive with main website
The time difference of the display of signal, Δ ε1iIt is the measurement error of the base station i and main website that are calculated in step 2);
Step 52) is by the initial position (x of target0, y0, z0) it is scheduled on the position for being no more than 200km apart from main website distance at random
It sets;
Formula inside step 51) is carried out Taylor expansion by step 53), and the removal above component of second order is available: Ψ=
ht-GtD, wherein each matrix vector is expressed as follows: Utilize weighted least-squares method solving target position
Set deviation are as follows:Wherein, Ri' indicate
Target initial position (the x of hypothesis0, y0, z0The distance between) and each base station;Q is to indicate that practicing measuring unit in each base station draws
The constant error of the relative error entered, each base station is identical;
Initial position is substituted into (x by step 54)0, y0, z0) distance (x is calculated in formulan, yn, zn), position deviation is calculated, it will
Thresholding is arranged to ε=150m, when | xn-x0|+|yn-y0|+|zn-z0| < ε terminates to calculate and result is fed back (xn, yn,
zn);It is such as unsatisfactory for, then in calculate next time, enables (xn, yn, zn)=(xn+Δx0, yn+Δy0, zn+Δz0) compute repeatedly;
Step 55) is by (xn, yn, zn) coordinate is converted into longitude and latitude height coordinate, by with inside included ADS-B message
Location information is compared, with to determine whether being an effective target.
A kind of positioning and false target recognition methods based on ADS-B message step-by-step counting above-mentioned, characterized in that described
The particular content of step 6) are as follows:
A) when it is multiple detect receive base stations receive echo signal, utilize multipoint positioning method calculate target real time position (xn, yn,
zn), then compared with the location information in the ADS-B message that target is sent, judge whether the target is spurious signal;
B) when detect receive base station receive it is bad or when there is base station crash, find out airplane motion using correlation function
The degree of association being intended to ADS-B message, formula are as follows:
Wherein, Δ ε is two relative pulse countings for detecing receipts base station for the same airbound target,It is Average pulse width, (xa,
ya,za) and (xb,yb,zb) it is two corresponding longitude and latitude height of position coordinates for detecing receipts base station, (xk,yk,zk) be airbound target from
The longitude and latitude elevation information of band, Δ P are the pulse number deviations being previously set, and contain space for estimating, judge aircraft position
The benchmark of range;
It is calculated by the way that aircraft longitude and latitude height is substituted into formula, judges the position and ADS- that multipoint location system detects
The distance between the included location information of B message error, if the location information that the aircraft carries is shown in the distance of setting
In band, then it is assumed that be an effective target;If location information is shown not in range, then it is assumed that be a false target.
Advantageous effects of the invention: using the aircraft of included S mode answering machine as existence conditions, it can be effective
In the case where not being equipped with wired link, calculate the clocking error between each receiving station;Then do not consider to eliminate error
But error is included in mathematical model by the means by programming, and is flown finally by continuous conclusion
The real time position of device.This method is greatly saved the cost of construction multipoint location system, and it is right around airport effectively to have ensured
The monitoring of airbound target compensates for the deficiency of secondary radar and ADS-B monitoring.
Detailed description of the invention
Fig. 1 is the total flow diagram of the present invention;
Fig. 2 is that different detect receives the step-by-step counting distribution that base station receives same ADS-B signal;
Fig. 3 is that base station receives the range difference schematic diagram for referring to aircraft;
Fig. 4 is arrangement of base stations schematic diagram.
Specific embodiment
The invention will be further described below in conjunction with the accompanying drawings.Following embodiment is only used for clearly illustrating the present invention
Technical solution, and not intended to limit the protection scope of the present invention.
For this method using the automatic broadcast data pulse signal of aircraft for carrying S mode answering machine, the ground base station Zhen Shou is logical
The Base Serial Number for collecting with reference to aircraft ADS-B message pulse signal and calculating preamble pulse is crossed, identical pulse is matched to
Collected information unification is sent to center base station and handled by signal, and main website is corrected by the method for weighted moving average and counted
Deviation establishes model and each measurement error detectd and receive base station and center base station is obtained by calculation.When there are numbers to reach multiple spot
When base station necessary to positioning, each detect receives base station by sending request signal to other unknown aircrafts of location information, leads to
It crosses the signal that feedback obtains and false target identification and multipoint positioning is carried out to aircraft, form the multipoint positioning phase based on ADS-B
Close monitoring system;When number of base stations necessary to multipoint positioning less than requiring, based on the location information being calculated, using packet
The range estimation method for holding space can effectively verify false target.Letter definition in attached drawing is as follows: Δ P21No. 2 base stations with
The received step-by-step counting difference in No. 1 base station, Δ P31It is No. 3 base stations and the received step-by-step counting difference in No. 1 base station.L1It is included S
The distance of No. 1 base station of aircraft distance of mode answering machine;LiIt is the base station aircraft distance i of included S mode answering machine
Distance.BSiIt is the code name of i-th of base station.
As shown in Figure 1, comprising the following steps:
Step 1, in daily life, each to detect the aircraft progress multipoint positioning for receiving base station to unknown geographical location information
When, the last positioning result deviation caused by clocking error is larger.The present invention proposes a kind of same without being laid with fixed optical cable
The time difference is walked, the method for cost is greatly reduced.The answer signal that the S mode answering machine carried using aircraft is sent is as signal amplitude
Source is penetrated, the pulse number received before since base station calculate top guide before receiving message is respectively detectd on ground.Calculate the preamble pulse received
Base Serial Number is matched to identical pulse signal, and sends it to main website, and main website is uniformly processed.Its same ADS-B message
The reception condition that pulse signal and its preamble pulse detect receipts base station in difference is as shown in Figure 2.
Step 2, when main website receive it is each detect receive base station be matched to the pulse with the same preamble pulses when, use
Following formula is calculated:ΔPi1It is each pulse and center base detectd and receive that base station receives
The pulse number stood subtracts each other,It is Average pulse width, Δ l is that the signal that S mode answering machine issues carries latitude and longitude information, is led to
Cross the S mode answering machine that is calculated to each receiving station distance;C is the speed of pulse transmitting, the i.e. light velocity.
So each at this time detect receives the measurement error of base station and main website and is exactlyThis
When it is each detect receive base station at a distance from main website erect-position as shown in figure 3, distance, delta l at this timei1Exactly use l1-li, the difference of umber of pulse
It is worth Δ Pi1Equal to the Y calculated in step 3 by formula101。
Step 3, over time, the error of the increase of target measurement number, whole system will over time
It is linear to increase.In order to solve the continuous growing concern of deviation, we are using the method for weighted moving average come update the system error.Tool
Body formula is as follows:
We define a set { Y1, Y2, Y3... ..., Y100, wherein by newest observation Δ Pi1It is defined as
Y1, by Y100Set is eliminated, the lower target number of other all values all increases by one.It is different that we give different observations
Flexible strategy, give larger flexible strategy value for the observed value close to time span of forecast, observed value farther away for the range prediction phase is then corresponding
Lesser flexible strategy value is given, to be equal to 1 since the sum of all flexible strategy adds up, so it is e that we, which might as well define flexible strategy,-x。
So the value that we actually obtain every time all uses Y by calculating101It indicates,It will
The Y newly obtained every time101It feeds back in main website, becomes and calculate each parameter detectd and receive base station and main website measurement error in step 2.
With reference to step 2, the pulse deviation being calculated is put into calculation formula and calculates by us to be detectd the measurement for receiving base station and misses
Difference.In measurement process, larger situation is differed with a upper measured value if there is measured value, generally it is considered that being some
It detects receipts base station not receiving signal or the case where system in case of system halt occur, we should cast out { Y1, Y2, Y3... ..., Y100, weight
Newly observation is included in start to be calculated again.
Step 4, according to the testing result of step 3, we determine other airbound targets with the method for multipoint positioning, judgement
Whether it is an effective target and its real time position.We set needs at least four and detect receipts base station thus, are formed
Three unilateral hyperboloids generate intersection point.
Model is established, the clocking error for each receiving station that step 3 obtains is included in model and is calculated, multipoint positioning is passed through
Method calculates target position.According to the base station number for receiving signal, the content of next step is selected.
Step 5, if the base station number received is not less than 4, specific implementation process is as follows:
Step 51) is since the multipoint positioning in three-dimensional space at least needs four base stations, by receiving station according to star cloth station:
Main website occupies center, chooses other three lines between secondary station and main website and forms a star, main website position coordinates are set as
(x1, y1, z1), it is assumed that spatial position coordinate is (x, y, z);RiIt is to detect the distance for receiving base station to echo signal for i-th, i=1,2,
3,4, R1iIt indicates that echo signal arrival the distance between main website and i-th of secondary station are poor, is indicated with equation: t1iEqual to each time for detecing and receiving base station and receiving the unknown aircraft of geography information
Subtract each measurement error detectd and receive base station and main website, i.e. t1i=T1i-Δε1i, wherein T1iIt is same that be the base station i receive with main website
The time difference of the display of signal, Δ ε1iIt is the measurement error of the base station i and main website that are calculated in step 2);
Step 52) is by the initial position (x of target0, y0, z0) it is scheduled on the position for being no more than 200km apart from main website distance at random
It sets;
Above formula is carried out Taylor expansion by step 53), and the removal above component of second order is available: Ψ=ht-Gtd.It is wherein each
Matrix vector is expressed as follows:
Utilize weighted least-squares method (WLS) solving target position deviation are as follows:
Wherein, Ri' indicate the target initial position (x assumed0, y0, z0The distance between) and each base station;C indicates pulse
The speed transmitted in space, the i.e. light velocity;t1iEqual to it is each detect to receive base station and receive time of the unknown aircraft of geography information subtract
Go each measurement error detectd and receive base station and main website;Q is the relative error for indicating to practice measuring unit introducing in each base station, generally
It is considered that the constant error of each base station is identical.
Step 54) we by initial position substitute into (x0, y0, z0) distance is calculated in formula, position deviation is then calculated, I
Thresholding is arranged to ε=150m, as | Δ x |+| Δ y |+| Δ z | < ε terminates to calculate result is simultaneously fed back (xn, yn,
zn);It is such as unsatisfactory for, then computes repeatedly.
Step 6, when the base station number received is less than 4, in order to identify the authenticity of ADS-B message, false flight is excluded
Interference of the target for monitoring station takes the be consistent method of verifying of purpose to judge that message is true and false.
A) when it is multiple detect receive base stations receive echo signal, utilize multipoint positioning method calculate target real time position (xn, yn,
zn), then compared with the location information in the ADS-B message that target is sent, judge whether the target is spurious signal;
B) when detect receive base station receive it is bad or when there is base station crash, find out airplane motion using correlation function
The degree of association being intended to ADS-B message, formula are as follows:
Wherein, Δ ε is two relative pulse countings for detecing receipts base station for the same airbound target,It is Average pulse width, (xa,
ya,za) and (xb,yb,zb) it is two corresponding longitude and latitude height of position coordinates for detecing receipts base station, (xk,yk,zk) be airbound target from
The longitude and latitude elevation information of band, Δ P are the pulse number deviations being previously set, and contain space for estimating, judge aircraft position
The benchmark of range;
It is calculated by the way that aircraft longitude and latitude height is substituted into formula, judges the position and ADS- that multipoint location system detects
The distance between the included location information of B message error.If the included location information of the aircraft is shown in what we were arranged
In distance band, then it is assumed that be an effective target;If location information is shown not in range, then it is assumed that be a false mesh
Mark.
Final result: the case where if it is a, according to calculating (xn, yn, zn) coordinate is converted into longitude and latitude height coordinate, then
By being compared with the location information inside included ADS-B message, with to determine whether being an effective target.
The case where if it is b, we obtain detecing the pulse deviation Y for receiving base station and main website by step 3101With it is each in step 2
It is a to detect the measurement error for receiving base station and main website, it has obtained arbitrarily detecing and has received base station and main website to the relative pulse of the same airbound target
Δ ε is counted, the longitude and latitude height coordinate that aircraft is carried substitutes into formula and calculates, and sees whether in the containing space as defined in us
Face.If judgement is an effective target, if not existing, judgement is a false target.
The above is only a preferred embodiment of the present invention, it is noted that for the ordinary skill people of the art
For member, without departing from the technical principles of the invention, several improvement and deformations can also be made, these improvement and deformations
Also it should be regarded as protection scope of the present invention.
Claims (9)
1. a kind of positioning and false target recognition methods based on ADS-B message step-by-step counting, characterized in that including walking as follows
It is rapid:
The aircraft that step 1) carries S mode answering machine sends data pulse, respectively detects and receives base station calculating preamble pulse Base Serial Number;
Each detect of step 2) receives the identical pulse signal of base station matching, is counted by pulse deviation, calculates each detect and receive base station
Measurement error;
Step 3) introduces measurement pulse deviation correction mechanism accurately to obtain the measurement error of each base station;
Step 4) given threshold is compared the base station number of the ADS-B signal received with threshold value, if being not less than threshold value,
Step 5) is then gone to, otherwise, goes to step 6);
Signal station according to star-like arrangement, while to objective emission signal, is considered the measurement error of each base station, passed through by step 5)
Three unilateral curved surfaces generate intersection point and determine position;
Step 6) by detects receive base station obtain target real time position or divide flight range, using purpose be consistent verification method exclude
Spurious signal interference.
2. a kind of positioning and false target recognition methods based on ADS-B message step-by-step counting according to claim 1,
It is characterized in, the particular content of the step 1) are as follows:
The answer signal of S mode answering machine is carried as signal radiation source using aircraft, ground, which is respectively detectd, receives base station and disappear from receiving
Pulse number before top guide starts to calculate before ceasing;
The preamble pulse Base Serial Number received is calculated, is matched to identical pulse signal, and send it to main website, main website is unified
Processing, ADS-B message should carry out pulse transmission since first sends pulse with the pulse distance of 0.5 microsecond, when detecing receipts base
After station receives same message preamble pulse, the step-by-step counting received before is calculated.
3. a kind of positioning and false target recognition methods based on ADS-B message step-by-step counting according to claim 2,
It is characterized in, the particular content of the step 2) are as follows:
Step-by-step counting before when receiving station i receives preface information is denoted as Pi, echo signal transmitting time be consistent, with
It receives the same inceptive impulse signal to start to calculate, calculates the intermediate propagation by road and receiving station receives pulse
Signal errors collects pulse signal received before when receiving station begins search for ADS-B preamble pulse;
Wherein, the error for receiving base station and center base station is each detectdΔPi1It is that each detect receives base station
The pulse number of the pulse and center base station that receive subtracts each other,It is Average pulse width, Δ l is the letter that S mode answering machine issues
Number included latitude and longitude information, the distance of the aircraft obtained by calculation with S mode answering machine to each receiving station;C is
The speed of pulse transmitting.
4. a kind of positioning and false target recognition methods based on ADS-B message step-by-step counting according to claim 2,
It is characterized in, the particular content of the step 3) are as follows:
Using the method for weighted moving average come update the system error, specific formula is as follows:Wherein, YiIt is
Step-by-step counting difference between the i base station observed, XiThe weight of i-th of the observation artificially assigned, according to time span of forecast
The distance of observed value give flexible strategy value, closer observed value gives smaller flexible strategy value, all flexible strategy values and add
To be equal to 1;The sum of observation is n;
Definition set { Y1, Y2, Y3... ..., Y100, wherein by newest observation Δ Pi1It is defined as Y1, by Y100It rejects
Gather out, the lower target number of other all values all increases by one;
It gives different observations different flexible strategy, larger flexible strategy value is given for the observed value close to time span of forecast, for distance
The farther away observed value of time span of forecast then accordingly gives lesser flexible strategy value.
5. a kind of positioning and false target recognition methods based on ADS-B message step-by-step counting according to claim 4,
It is characterized in, definition flexible strategy are e-i, Y is all used by calculating the value that reality obtains every time101It indicates,
The Y that will newly obtain every time101It feeds back in main website, as Δ Pi1。
6. a kind of positioning and false target recognition methods based on ADS-B message step-by-step counting according to claim 1,
It is characterized in, sets the threshold to 4 in the step 4).
7. a kind of positioning and false target recognition methods based on ADS-B message step-by-step counting according to claim 6,
It is characterized in, the particular content of the step 5) are as follows:
The real time position of other airbound targets is determined using the method for multipoint positioning, is at least needed four to detect and is received base station, is formed
Three unilateral hyperboloids generate intersection point;Model is established, the measurement error for each receiving station that step 3) obtains is included in model
It calculates, target position is calculated by the method for multipoint positioning.
8. a kind of positioning and false target recognition methods based on ADS-B message step-by-step counting according to claim 7,
It is characterized in, includes the following steps:
Step 51) is since the multipoint positioning in three-dimensional space at least needs four base stations, by receiving station according to star cloth station: main website
Center is occupied, other three lines between secondary station and main website is chosen and forms a star, main website position coordinates are set as (x1, y1,
z1), it is assumed that spatial position coordinate is (x, y, z);RiIt is to detect the distance for receiving base station to echo signal for i-th, i=1,2,3,4, R1i
It indicates that echo signal arrival the distance between main website and i-th of secondary station are poor, is indicated with equation: t1iEqual to each time for detecing and receiving base station and receiving the unknown aircraft of geography information
Subtract each measurement error detectd and receive base station and main website, i.e. t1i=T1i-Δε1i, wherein T1iIt is same that be the base station i receive with main website
The time difference of the display of signal, Δ ε1iIt is the measurement error of the base station i and main website that are calculated in step 2);
Step 52) is by the initial position (x of target0, y0, z0) it is scheduled on the position for being no more than 200km apart from main website distance at random;
Formula inside step 51) is carried out Taylor expansion by step 53), and the removal above component of second order is available: Ψ=ht-
GtD, wherein each matrix vector is expressed as follows: Utilize weighted least-squares method solving target
Position deviation are as follows:Wherein, R 'iTable
Show the target initial position (x of hypothesis0, y0, z0The distance between) and each base station;Q is to indicate to practice measuring unit in each base station
The constant error of the relative error of introducing, each base station is identical;
Initial position is substituted into (x by step 54)0, y0, z0) distance (x is calculated in formulan, yn, zn), position deviation is calculated, by thresholding
It is arranged to ε=150m, when | xn-x0|+|yn-y0|+|zn-z0| < ε terminates to calculate and by result (xn, yn, zn) feed back;Such as
It is unsatisfactory for, then in calculate next time, enables (xn, yn, zn)=(xn+Δx0, yn+Δy0, zn+Δz0) compute repeatedly;
Step 55) is by (xn, yn, zn) coordinate is converted into longitude and latitude height coordinate, by with the position inside included ADS-B message
Information is compared, with to determine whether being an effective target.
9. a kind of positioning and false target recognition methods based on ADS-B message step-by-step counting according to claim 1,
It is characterized in, the particular content of the step 6) are as follows:
A) when it is multiple detect receive base stations receive echo signal, utilize multipoint positioning method calculate target real time position (xn, yn, zn),
Then compared with the location information in the ADS-B message that target is sent, judge whether the target is spurious signal;
B) when detect receive base station receive it is bad or when there is base station crash, using correlation function find out airplane motion with
The degree of association that ADS-B message is intended to, formula are as follows:
Wherein, Δ ε is two relative pulse countings for detecing receipts base station for the same airbound target,It is Average pulse width, (xa,
ya,za) and (xb,yb,zb) it is two corresponding longitude and latitude height of position coordinates for detecing receipts base station, (xk,yk,zk) be airbound target from
The longitude and latitude elevation information of band, Δ P are the pulse number deviations being previously set, and contain space for estimating, judge aircraft position
The benchmark of range;
It is calculated by the way that aircraft longitude and latitude height is substituted into formula, judges the position and ADS-B report that multipoint location system detects
The distance between the included location information of text error, if the location information that the aircraft carries is shown in the distance band of setting
It is interior, then it is assumed that be an effective target;If location information is shown not in range, then it is assumed that be a false target.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810586091.XA CN108957503B (en) | 2018-06-08 | 2018-06-08 | Positioning and false target identification method based on ADS-B message pulse counting |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810586091.XA CN108957503B (en) | 2018-06-08 | 2018-06-08 | Positioning and false target identification method based on ADS-B message pulse counting |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108957503A true CN108957503A (en) | 2018-12-07 |
CN108957503B CN108957503B (en) | 2022-07-26 |
Family
ID=64493814
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810586091.XA Active CN108957503B (en) | 2018-06-08 | 2018-06-08 | Positioning and false target identification method based on ADS-B message pulse counting |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108957503B (en) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109886116A (en) * | 2019-01-21 | 2019-06-14 | 中国人民解放军战略支援部队信息工程大学 | Specific emitter identification method and device based on ADS-B |
CN110061801A (en) * | 2019-04-25 | 2019-07-26 | 电子科技大学 | Abnormal object discrimination system and method in a kind of aerial surveillance systems |
CN110940972A (en) * | 2019-12-09 | 2020-03-31 | 中国民航大学 | Method for extracting S-mode signal arrival time of multi-preamble pulse combined filtering detection |
WO2021097680A1 (en) * | 2019-11-19 | 2021-05-27 | 华为技术有限公司 | Positioning method and apparatus |
CN113050137A (en) * | 2021-03-09 | 2021-06-29 | 江西师范大学 | Multi-point cooperative measurement spatial information acquisition method |
CN113050052A (en) * | 2021-03-08 | 2021-06-29 | 四川九洲空管科技有限责任公司 | Multi-target response simulation method for friend or foe identification system |
CN113589225A (en) * | 2021-07-29 | 2021-11-02 | 四川九洲电器集团有限责任公司 | Aerial target positioning method and device |
CN114279466A (en) * | 2021-12-23 | 2022-04-05 | 中国电子科技集团公司第十四研究所 | Sensor error correction method |
CN114626409A (en) * | 2022-02-21 | 2022-06-14 | 中铁第四勘察设计院集团有限公司 | Near-fault acceleration pulse identification method, storage medium and computer equipment |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102323567A (en) * | 2011-06-13 | 2012-01-18 | 四川大学 | A kind of method that detects the ADS-B false target |
EP2738972A2 (en) * | 2012-11-30 | 2014-06-04 | Honeywell International Inc. | Parallel-frequency partially-coherent reception of pulse-position modulated ADS-B messages |
CN104378253A (en) * | 2014-09-23 | 2015-02-25 | 北京航空航天大学 | Header locating method of S-mode acknowledgement signal in ADS-B system |
CN106301413A (en) * | 2016-08-16 | 2017-01-04 | 中电科现代导航(西安)科技有限公司 | The airborne anti-interference reception of ADS B and processing method |
CN106846919A (en) * | 2017-01-16 | 2017-06-13 | 南京航空航天大学 | A kind of four-dimensional dynamic track prediction method based on ADS B information updatings |
-
2018
- 2018-06-08 CN CN201810586091.XA patent/CN108957503B/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102323567A (en) * | 2011-06-13 | 2012-01-18 | 四川大学 | A kind of method that detects the ADS-B false target |
EP2738972A2 (en) * | 2012-11-30 | 2014-06-04 | Honeywell International Inc. | Parallel-frequency partially-coherent reception of pulse-position modulated ADS-B messages |
CN104378253A (en) * | 2014-09-23 | 2015-02-25 | 北京航空航天大学 | Header locating method of S-mode acknowledgement signal in ADS-B system |
CN106301413A (en) * | 2016-08-16 | 2017-01-04 | 中电科现代导航(西安)科技有限公司 | The airborne anti-interference reception of ADS B and processing method |
CN106846919A (en) * | 2017-01-16 | 2017-06-13 | 南京航空航天大学 | A kind of four-dimensional dynamic track prediction method based on ADS B information updatings |
Non-Patent Citations (2)
Title |
---|
MING-SHIH HUANG,ET AL: "Tracking of Noncooperative Airborne Targets Using ADS-B Signal and Radar Sensing", 《INTERNATIONAL JOURNAL OF AEROSPACE ENGINEERING》 * |
周建红: "基于MDS系统的定位算法与同步技术研究", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109886116A (en) * | 2019-01-21 | 2019-06-14 | 中国人民解放军战略支援部队信息工程大学 | Specific emitter identification method and device based on ADS-B |
CN110061801A (en) * | 2019-04-25 | 2019-07-26 | 电子科技大学 | Abnormal object discrimination system and method in a kind of aerial surveillance systems |
WO2021097680A1 (en) * | 2019-11-19 | 2021-05-27 | 华为技术有限公司 | Positioning method and apparatus |
CN110940972A (en) * | 2019-12-09 | 2020-03-31 | 中国民航大学 | Method for extracting S-mode signal arrival time of multi-preamble pulse combined filtering detection |
CN110940972B (en) * | 2019-12-09 | 2023-04-07 | 中国民航大学 | Method for extracting S-mode signal arrival time of multi-preamble pulse combined filtering detection |
CN113050052B (en) * | 2021-03-08 | 2022-08-09 | 四川九洲空管科技有限责任公司 | Multi-target response simulation method for friend or foe identification system |
CN113050052A (en) * | 2021-03-08 | 2021-06-29 | 四川九洲空管科技有限责任公司 | Multi-target response simulation method for friend or foe identification system |
CN113050137A (en) * | 2021-03-09 | 2021-06-29 | 江西师范大学 | Multi-point cooperative measurement spatial information acquisition method |
CN113589225A (en) * | 2021-07-29 | 2021-11-02 | 四川九洲电器集团有限责任公司 | Aerial target positioning method and device |
CN114279466A (en) * | 2021-12-23 | 2022-04-05 | 中国电子科技集团公司第十四研究所 | Sensor error correction method |
CN114279466B (en) * | 2021-12-23 | 2024-02-27 | 中国电子科技集团公司第十四研究所 | Sensor error correction method |
CN114626409A (en) * | 2022-02-21 | 2022-06-14 | 中铁第四勘察设计院集团有限公司 | Near-fault acceleration pulse identification method, storage medium and computer equipment |
CN114626409B (en) * | 2022-02-21 | 2023-09-26 | 中铁第四勘察设计院集团有限公司 | Near fault acceleration pulse identification method, storage medium and computer equipment |
Also Published As
Publication number | Publication date |
---|---|
CN108957503B (en) | 2022-07-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108957503A (en) | Positioning and false target recognition methods based on ADS-B message step-by-step counting | |
CN109541538B (en) | Method for realizing distributed ADS-B passive positioning anti-cheating function based on opportunity reference source | |
CA2144291C (en) | Control target surveillance system | |
ES2361776T3 (en) | USE OF PRECISION GEOMETRIC DILUTION (GDOP) TO SELECT THE BEST RECEIVER GROUP IN A MULTILATERATION SYSTEM. | |
CN107484136A (en) | A kind of UWB multi-tag localization methods based on polling dispatching | |
CN103562742B (en) | Synchronous and the forwarding of determination signal of geographic registration application and the method and system of the time variations in propagating for range observation, driver | |
CN105510876A (en) | Electromagnetic wave propagation characteristic-based indoor distance measurement positioning method | |
CN106872942A (en) | For the positioning precision calculation method of Distributed Multi positioning monitoring system | |
CN104684014B (en) | Mobile communications network transmission testing method | |
CN108882149A (en) | NLOS apart from dependent probability compensates localization method | |
CN111781570B (en) | Radar online precision analysis method based on real-time ADS-B data | |
CN110764050A (en) | Simulation test method for multipoint positioning system | |
RU2333538C2 (en) | Method of indication of observed object position | |
CN104501807A (en) | Indoor location method based on geomagnetic field and historical localization track | |
CN104318094B (en) | A kind of radar network composite Data Association of formation flight | |
CN117075151A (en) | Navigation decoy interference system and method based on artificial intelligence | |
CN107561530A (en) | A kind of target identification method based on collaboration passive detection multidimensional information | |
CN110389353A (en) | A kind of forest fire point precise positioning method and system | |
CN113391267A (en) | Frequency spectrum detection system positioning method based on ATDOA algorithm | |
RU2453995C1 (en) | Method to receive radio signals from sources of radio radiations | |
TWI745153B (en) | Method, device and system of ultra-wideband localization | |
CN108535690A (en) | A kind of signal matching method of multipoint positioning scene monitoring system | |
CN112394320B (en) | Indoor high-precision centroid positioning method based on support vector machine | |
CN111866725B (en) | People stream detection method based on WIFI probe technology | |
CN104181549A (en) | Method for positioning aircraft based on 1090ES data chain ground stations |
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 |