CN109657340A - A kind of 1090ES signal overlaps the estimation method of probability - Google Patents
A kind of 1090ES signal overlaps the estimation method of probability Download PDFInfo
- Publication number
- CN109657340A CN109657340A CN201811541773.5A CN201811541773A CN109657340A CN 109657340 A CN109657340 A CN 109657340A CN 201811541773 A CN201811541773 A CN 201811541773A CN 109657340 A CN109657340 A CN 109657340A
- Authority
- CN
- China
- Prior art keywords
- signal
- overlapping
- probability
- averagely
- factor
- 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
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G08—SIGNALLING
- G08G—TRAFFIC CONTROL SYSTEMS
- G08G5/00—Traffic control systems for aircraft, e.g. air-traffic control [ATC]
- G08G5/0073—Surveillance aids
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/08—Probabilistic or stochastic CAD
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Aviation & Aerospace Engineering (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
The invention discloses the estimation methods that a kind of 1090ES signal overlaps probability, emulation and modeled analysis are carried out by the arrival time of secondary radar Mode A/C answer signal and 1090ES signal to radio frequency operation on 1090MHz frequency point, proposes the algorithm for estimating that the 1090ES signal based on Poisson distribution model overlaps probability.The algorithm can be according to secondary radar portion number, frame-scan period and inquiry period in scene, and aerial target number parameter, when it is assumed that Mode A/C answer signal, 1090ES signal reach certain ADS-B earth station when arrival time random distribution, can directly it be estimated by overlapping probability of the parsing mathematic(al) representation to 1090ES signal.
Description
Technical field
The invention belongs to blank pipe surveillance technology fields, and in particular to a kind of 1090ES signal overlaps the estimation method of probability.
Background technique
Blank pipe surveillance technology is one of most important technology in air traffic control system (ATCS), and the task of monitoring is to aircraft
And other targets are reliably detected, and position, state and the warning information of accurate aircraft and other detection targets are provided,
It provides the real-time dynamic information of target for air traffic control system information management.Currently, commonly monitoring that means are main in air traffic control system
Once surveillance radar (PSR), scene surveillance radar, secondary surveillance radar (SSR), Automatic dependent surveillance broadcast (ADS-B)
With multipoint positioning (MLAT).In global sail plan and Chinese chunking upgrading plan, within current and future suitable long-term,
It still can be main monitoring means with above-mentioned several surveillance technologies.
ADS-B is the aircraft operation monitoring system based on GPS GPS and vacant lot, the communication of absolutely empty data-link
System.ADS-B system is realized ground by way of itself relevant information of aircraft active transmission or in the air to the monitoring of target
Function.International Civil Aviation Organization (ICAO) determines it as the Main way of the following surveillance technology development, and international airline circle is being accumulated
Pole promotes the application of this technology.Compared with traditional secondary radar system, ADS-B system is capable of providing real-time and accurate aviation
The information such as device position, speed, and the system has construction investment small, the advantages that maintenance cost is low, long service life.1090ES
(1090Extended Squitter, 1090 extension texts) are one of the three kinds of data-links for realizing ADS-B function, radio frequency work
Make in 1090MHz.Identical as the format of secondary radar S mode answer signal, one frame signal is made of a series of pulse, this
The important informations such as position, speed, course, the state of aircraft are contained in a little pulses.
With the fast development of air-transport industry, number of vehicles continues to increase, secondary radar, multipoint positioning and 1090ES
This frequency point of 1090MHz is used etc. main monitoring means.For busy terminal airspace or air route, with the increasing of destination number
Add, the problem of co-channel interference of 1090MHz frequency point is got worse, and co-channel interference can make signal missing inspection and lead to track jump point, sternly
Targetpath interrupts when weight.Since surveillance technology is in the critical role of air traffic control system, 1090MHz signal is studied, especially
The probability and interference mitigation technology that 1090ES signal is disturbed have ADS-B earth station tentative plan of construction program extremely important
Directive function.When the overlapping probability of 1090ES signal is excessive, it is considered as connecing using multiple fan antenna point airspaces progress signal
It receives, overlaps probability to reduce signal.
Summary of the invention
It is an object of the invention to overcome the deficiencies of the prior art and provide the estimations that a kind of 1090ES signal overlaps probability
Method is overlapped the probability of crosstalk in the signal of 1090MHz frequency point by other radio frequency operations applied to estimation 1090ES signal, is
The design of ADS-B earth station construction scheme provides foundation.
In order to solve the above technical problems, the present invention provides the estimation method that a kind of 1090ES signal overlaps probability, it is special
Sign is, including following procedure:
S1 calculates the signal dutyfactor factor of the work on 1090MHz frequency point；
S2 is calculated according to the relationship between the emulation averagely overlapping number factor fitted and the signal dutyfactor factor
The averagely overlapping number factor；
S3 calculates averagely overlapping frequency according to the averagely overlapping number factor；
S4 calculates 1090ES signal and overlaps probability according to averagely overlapping frequency.
Further, in S1, signal dutyfactor of the work on 1090MHz frequency point is calculated because of subformula are as follows:
Wherein, fS-112And fS-56Respectively S-112 and S-56 average every batch of objective emission number per second；tS-112With
tS-56Respectively S-112 and S-56 signal durations；It is averaged every batch of target response per second for corresponding i-th secondary radar
Number；tA/CFor A/C signal duration；NSSRFor secondary radar portion number；kDDimension be the second.
Further, in S2, averagely overlapping number factor kλWith signal dutyfactor factor kDBetween relationship are as follows:
kλ=6.78kD-0.0032
Wherein, two coefficients in above formula are obtained by Monte Carlo simulation.
Further, simulation process are as follows:
It is assumed that the aerial target transmitting various types of signal that ADS-B ground station reception arrives is in the unit time in 1090MHz channel
It is reached at random in section；
In conjunction with the duration of each signal, to judge whether each signal has with other signals overlapping, and statistics has overlapping signal
Number calculates signal and overlaps probability；
The averagely overlapping number factor of probability calculation is overlapped by signal, draws out the averagely overlapping number factor and signal dutyfactor
Relationship between the factor.
Further, averagely overlapping frequency process is calculated are as follows:
Scheduled target lot number N, which is multiplied by, according to the resulting averagely overlapping number factor of calculating calculates average overlapping generation time
Number.
Further, it using classical Poisson distribution model, calculates signal and overlaps probability.
Further, the calculation formula of the overlapping probability of signal is as follows:
P=1-e-λ
Wherein, λ is averagely overlapping frequency.
Compared with prior art, the beneficial effects obtained by the present invention are as follows being: evaluation method proposed by the present invention is based on Poisson
Distributed model according to aerial target sum, secondary radar portion number and its frame-scan period in scene and inquires the parameters such as period, by
Digital analytical expression directly calculates the probability that 1090ES signal is overlapped by other signals.
Detailed description of the invention
Fig. 1 is the linear relationship between the averagely overlapping number factor and the signal dutyfactor factor；
Fig. 2 is the flow diagram of the method for the present invention；
Fig. 3 is the corresponding relationship between target lot number and secondary radar portion number and the overlapping probability of signal.
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.
By the aeronautical telecommunication specification attachment 10 of International Civil Aviation Organization (ICAO) publication it is found that 1090ES signal and secondary radar
Mode S answer signal format is identical, and 1090ES signal is divided into two kinds of long format (S-112) and short format (S-56).Each
S-112 signal duration 120us, transmitting 5.2 times per second of each target；Each S-56 signal duration 64us, each target are every
Second transmitting 1 time.Each target will not only reply the inquiry of the various modes from secondary radar, and active broadcast is also wanted to emit
1090ES signal.Mode A is identical with mode C answer signal format, when each Mode A/C answer signal a length of 20.75us, often
Second response times inquire that the parameters such as period are determined by secondary radar portion number, frame-scan period, Mode A/C in scene.Due to
The important parameters such as position, speed, the S mould code of target are all contained in S-112 contents, and therefore, the present invention is directed to 1090ES
Its overlapping probability of S-112 signal estimations.
A kind of 1090ES signal of the invention overlaps the estimation method of probability, shown in Figure 2, including following procedure:
Step S1 calculates the signal dutyfactor factor of the work on 1090MHz frequency point；
Target emergence pattern A/C is answered when according to secondary radar portion number in scene, every secondary radar is inswept certain batch of target
The parameters such as number, the target total number answered calculate signal dutyfactor factor kD, i.e., each target emission mode A/C response per second
The signals of all radio frequency operations on 1090MHz frequency point such as S-56 signal, S-112 signals of 1090ES, 1090ES signals
Total duration kD；
Wherein, fS-112And fS-56Respectively S-112 and S-56 average every batch of objective emission number per second；tS-112With
tS-56Respectively S-112 and S-56 signal durations；It is averaged every batch of target response per second for corresponding i-th secondary radar
Number；tA/CFor A/C signal duration；NSSRFor secondary radar portion number；kDDimension be the second.
Step S2, according to the linear pass between the emulation averagely overlapping number factor fitted and the signal dutyfactor factor
System calculates the averagely overlapping number factor.
Because each signal reach ADS-B earth station time it is irregular follow, most common hypothesis is exactly: it is assumed that
In 1090MHz channel, ADS-B ground station reception to aerial target transmitting various types of signal reached at random in unit interval.
Under the assumption, it is overlapping to judge whether each signal has with other signals for the duration in conjunction with each signal, and statistics has
Overlapping signal number calculates signal and overlaps probability.The averagely overlapping number factor of probability calculation is overlapped by signal again, can be drawn
Out as shown in Figure 1 as a result, and fitting fixed linear equation.
Resulting kDAveragely overlapping number factor k is calculated with following formulaλ；
kλ=6.78kD-0.0032
Wherein, two coefficients in above formula are obtained by Monte Carlo simulation.Coefficient 6.78 is from (tS-112+tA/C)/
tA/CRatio, this is because 1090MHz channel internal schema A/C answer signal number is more, 1090ES occur it is overlapping be nearly all with
Mode A/C answer signal generates, and there is a situation where overlapping considerably less between 1090ES signal；Coefficient -0.0032 comes from Meng Teka
Lip river emulation.
Step S3, according to the resulting k of calculatingλThe target lot number N assumed in advance calculates averagely overlapping event frequency
λ。
Overlapping event be averaged frequency calculation formula it is as follows:
λ=kλN
Step S4 calculates 1090ES S-112 signal and overlaps probability according to averagely overlapping number.
Because each time of arrival (toa) is random within the unit time, it is applicable in Poisson distribution model.Utilize classical Poisson distribution
Model, the calculation formula that signal overlaps probability are as follows:
P=1-e-λ
Embodiment
In certain scene, there are 4 secondary radars, every secondary radar run-down time is 4 seconds, per an inswept target
20 responses can be caused in the process, then estimate that the specific implementation step of overlapping probability is as follows in different target number:
(1) signal dutyfactor factor k is calculated according to above-mentioned condition and following formulaD, i.e., each target emission mode A/C per second answers
Answer the letters of all radio frequency operations on 1090MHz frequency point such as S-56 signal, S-112 signals of 1090ES, 1090ES signals
Number total duration；
(2) resulting k is calculated according to step (1)DAveragely overlapping number factor k is calculated with following formulaλ；
kλ=6.78kD-0.0032
=0.00428
(3) initial value of target number is set as N=50, calculates resulting k according to step (2)λIt is average overlapping with calculating
Number λ；
λ=kλ× N=0.00428 × 50=0.214
(4) 1090ES S-112 signal is calculated using Poisson distribution model overlap probability
P=1-e-λ=1-e-0.214=0.19
It in general, is all that overlapping probability accordingly is calculated according to target lot number and secondary radar portion number, for blank pipe
It is relatively more intuitive for personnel.
(5) different target number: 100-1000 (be divided into 50) is set, when repeating step (3) and calculating different target number
Overlapping probability and record data；
(6) again assume that the secondary radar portion number in scene is respectively 3,2,1,0, repeat the above steps (1)-(5), calculates
Overlapping probability in different secondary radar portion number situation simultaneously records data；Overlapping probability is closely related with secondary radar portion number, because
Aerial target is constantly inquired for secondary radar, so that the A/C response of target can be caused, and the biggish A/C answer signal of density
It can be generated with 1090ES signal overlapping.
(7) curve is drawn according to the data that step (5) and (6) are recorded, (SSR is secondary radar in figure as shown in Figure 3
English abbreviation), with the increase of target lot number and the increase of secondary radar portion number, 1090ES signal overlaps probability and significantly rises.
It should be understood by those skilled in the art that, embodiments herein can provide as method, system or computer program
Product.Therefore, complete hardware embodiment, complete software embodiment or reality combining software and hardware aspects can be used in the application
Apply the form of example.Moreover, it wherein includes the computer of computer usable program code that the application, which can be used in one or more,
The computer program implemented in usable storage medium (including but not limited to magnetic disk storage, CD-ROM, optical memory etc.) produces
The form of product.
The application is referring to method, the process of equipment (system) and computer program product according to the embodiment of the present application
Figure and/or block diagram describe.It should be understood that every one stream in flowchart and/or the block diagram can be realized by computer program instructions
The combination of process and/or box in journey and/or box and flowchart and/or the block diagram.It can provide these computer programs
Instruct the processor of general purpose computer, special purpose computer, Embedded Processor or other programmable data processing devices to produce
A raw machine, so that being generated by the instruction that computer or the processor of other programmable data processing devices execute for real
The device for the function of being specified in present one or more flows of the flowchart and/or one or more blocks of the block diagram.
These computer program instructions, which may also be stored in, is able to guide computer or other programmable data processing devices with spy
Determine in the computer-readable memory that mode works, so that it includes referring to that instruction stored in the computer readable memory, which generates,
Enable the manufacture of device, the command device realize in one box of one or more flows of the flowchart and/or block diagram or
The function of being specified in multiple boxes.
These computer program instructions also can be loaded onto a computer or other programmable data processing device, so that counting
Series of operation steps are executed on calculation machine or other programmable devices to generate computer implemented processing, thus in computer or
The instruction executed on other programmable devices is provided for realizing in one or more flows of the flowchart and/or block diagram one
The step of function of being specified in a box or multiple boxes.
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 improvements and modifications, these improvements and modifications can also be made
Also it should be regarded as protection scope of the present invention.
Claims (7)
1. the estimation method that a kind of 1090ES signal overlaps probability, characterized in that including following procedure:
S1 calculates the signal dutyfactor factor of the work on 1090MHz frequency point；
S2 is calculated average according to the relationship between the emulation averagely overlapping number factor fitted and the signal dutyfactor factor
The overlapping number factor；
S3 calculates averagely overlapping frequency according to the averagely overlapping number factor；
S4 calculates signal and overlaps probability according to averagely overlapping frequency.
2. the estimation method that a kind of 1090ES signal according to claim 1 overlaps probability, characterized in that in S1, calculate
The signal dutyfactor to work on 1090MHz frequency point is because of subformula are as follows:
Wherein, fS-112And fS-56Respectively S-112 and S-56 average every batch of objective emission number per second；tS-112And tS-56Point
It Wei not S-112 and S-56 signal durations；It is averaged every batch of target response number per second for corresponding i-th secondary radar；
tA/CFor A/C signal duration；NSSRFor secondary radar portion number；kDDimension be the second.
3. the estimation method that a kind of 1090ES signal according to claim 1 overlaps probability, characterized in that average in S2
Overlapping number factor kλWith signal dutyfactor factor kDBetween relationship are as follows:
kλ=6.78kD-0.0032
Wherein, two coefficients in above formula are obtained by Monte Carlo simulation.
4. the estimation method that a kind of 1090ES signal according to claim 3 overlaps probability, characterized in that simulation process
Are as follows:
It is assumed that the aerial target transmitting various types of signal that ADS-B ground station reception arrives is in unit interval in 1090MHz channel
It is random to reach；
In conjunction with the duration of each signal, to judge whether each signal has with other signals overlapping, and statistics has overlapping signal number, counts
It calculates signal and overlaps probability；
The averagely overlapping number factor of probability calculation is overlapped by signal, draws out the averagely overlapping number factor and the signal dutyfactor factor
Between relationship.
5. the estimation method that a kind of 1090ES signal according to claim 1 overlaps probability, characterized in that calculate average hand over
Folded frequency process are as follows:
The averagely overlapping frequency of scheduled target lot number calculating is multiplied by according to the resulting averagely overlapping number factor is calculated.
6. the estimation method that a kind of 1090ES signal according to claim 1 overlaps probability, characterized in that utilize classical pool
Loose distributed model calculates signal and overlaps probability.
7. the estimation method that a kind of 1090ES signal according to claim 6 overlaps probability, characterized in that signal is overlapping general
The calculation formula of rate is as follows:
P=1-e-λ
Wherein, λ is averagely overlapping frequency.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811541773.5A CN109657340B (en) | 2018-12-17 | 2018-12-17 | 1090ES signal overlapping probability estimation method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811541773.5A CN109657340B (en) | 2018-12-17 | 2018-12-17 | 1090ES signal overlapping probability estimation method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109657340A true CN109657340A (en) | 2019-04-19 |
| CN109657340B CN109657340B (en) | 2023-03-14 |
Family
ID=66113237
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811541773.5A Active CN109657340B (en) | 2018-12-17 | 2018-12-17 | 1090ES signal overlapping probability estimation method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109657340B (en) |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104394107A (en) * | 2014-10-24 | 2015-03-04 | 中国电子科技集团公司第二十八研究所 | A method of decoding common-frequency overlapping of 1090 ES signals based on secondary inflexion-point method |
-
2018
- 2018-12-17 CN CN201811541773.5A patent/CN109657340B/en active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104394107A (en) * | 2014-10-24 | 2015-03-04 | 中国电子科技集团公司第二十八研究所 | A method of decoding common-frequency overlapping of 1090 ES signals based on secondary inflexion-point method |
Non-Patent Citations (2)
| Title |
|---|
| 孙清清等: "1090MHz信号同频干扰与窜扰概率探究", 《通信技术》 * |
| 胡来招等: "多信号环境下脉冲重叠概率的研究", 《电子对抗》 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109657340B (en) | 2023-03-14 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US9558670B1 (en) | Method and system for air traffic rerouting for airspace constraint resolution | |
| US9171473B1 (en) | Method and system for dynamic automated corrections to weather avoidance routes for aircraft in en route airspace | |
| US9842506B2 (en) | Systems and methods for conflict detection using dynamic thresholds | |
| US8437949B2 (en) | Method and systems for determining required interval management performance (RIMP) | |
| US20180240348A1 (en) | Methods and systems for probabilistic spacing advisory tool (psat) | |
| US8775062B2 (en) | Terminal aircraft sequencing and conflict resolution | |
| CN107656266B (en) | S-mode secondary radar information fusion system based on ADS-B and fusion method thereof | |
| CN104077928A (en) | Position information updating method suitable for irregular-flying navigation through points out of airways | |
| US11442163B2 (en) | Device and method for an aircraft bird congestion indicator system | |
| CN104966408A (en) | GPS positioning data compensation method | |
| CN101377887A (en) | Statistical method and apparatus for airliner delay | |
| CN101630446A (en) | Method for evaluating aircraft state based on broadcast type automatic correlative monitoring data and system thereof | |
| CN113240212B (en) | Data processing method, electronic device and medium for generating flight trajectory | |
| CN102831487A (en) | Flow predicted result verification method based on historical scheduled flight running data analysis | |
| CN101393691B (en) | Airport arriving flow alarm implementing method and system | |
| CN101465066A (en) | Method for obtaining airplane vertical occupation ratio | |
| Gariel et al. | An algorithm for conflict detection in dense traffic using ADS-B | |
| JP6234141B2 (en) | Movement plan management system, movement plan management method, and movement plan management program | |
| Mott | Estimation of aircraft distances using transponder signal strength information | |
| CN109657340A (en) | A kind of 1090ES signal overlaps the estimation method of probability | |
| Andrés et al. | Probabilistic analysis of air traffic in adverse weather scenarios | |
| US20230221173A1 (en) | Airport sound and noise management methods and systems | |
| US20190122566A1 (en) | Method for securing a provisional itinerary for an aircraft, corresponding system and computer program | |
| Zhu et al. | Simulation and Analysis for Overlapping Probability of ADS-B 1090ES Signal | |
| CN111816005A (en) | Remote piloted aircraft environment monitoring optimization method based on ADS-B |
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 |