CN107153174B - GJB661-89 standard beacon signal identification method - Google Patents
GJB661-89 standard beacon signal identification method Download PDFInfo
- Publication number
- CN107153174B CN107153174B CN201710353779.9A CN201710353779A CN107153174B CN 107153174 B CN107153174 B CN 107153174B CN 201710353779 A CN201710353779 A CN 201710353779A CN 107153174 B CN107153174 B CN 107153174B
- Authority
- CN
- China
- Prior art keywords
- pulse
- beacon
- information
- bcn
- queue
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/28—Details of pulse systems
Abstract
The invention discloses a method for identifying GJB661-89 standard beacon signals. The implementation scheme is as follows: 1. firstly, obtaining beacon information fed back by a group of beacon machines, and then inquiring beacon pulse positions for decoding; 2. comparing the pulse coding information obtained by decoding with the existing pulse coding information in the pulse coding queue, and updating the pulse coding queue; 3. judging whether the antenna scans the boundary or not through the radar beam pointing angle, acquiring a group of new beacon position information if the antenna does not scan the boundary, calculating the pulse coding information in the pulse coding queue and storing the pulse coding information in the beacon machine information queue if the antenna scans the boundary, and clearing the information in the pulse coding queue; 4. and displaying the corresponding beacon machine position and the corresponding beacon machine code according to the decoding instruction and clearing the beacon machine information queue. The invention further improves the identification capability and the positioning precision of the beacon machine, can effectively identify a plurality of standard beacon machine signals, meets the higher requirement of positioning and navigation, and can be used for various radars.
Description
Technical Field
The invention belongs to the technical field of radars, and particularly relates to a method for identifying a double-pulse coded signal and a multi-pulse coded signal in a beacon machine which meet the requirements of GJB661-89 standard, which can be used for airborne radars, ground radars, carrier-borne radars and shore-based radars.
Background
The beacon machine is a positioning navigation equipment, it can receive radar inquiry signal of a certain frequency band, then use another frequency of the same frequency band to make response, according to the difference of carrying mode, the beacon machine can be divided into portable type, machine-mounted type and ship-mounted type.
The beacon response signal is divided into a single pulse signal, a double pulse coding signal and a multi-pulse coding signal, wherein the double pulse coding signal is coded by the time interval of two pulses; the multi-pulse coded signal contains at most six pulses per group, the distance between each pulse is 3us, wherein the first pulse and the sixth pulse must exist in different codes, the distance between the first pulse and the sixth pulse is 15us, the existence or non-existence of the fifth pulse, the fourth pulse, the third pulse and the second pulse represents the value of the corresponding bit in 4-bit binary data, the existence of the pulse represents 1, the non-existence represents 0, and the multi-pulse coded signal is used for representing codes 1-15, for example: the presence of the fourth and second pulses represents binary data 0101, encoded 5.
In practical use, a plurality of beacons are usually arranged in a target area, each beacon has a unique code, a radar transmits a query signal meeting the GJB661-89 standard to the target area, the beacon responds to the coded signal after receiving the standard query signal, and the beacon code of interest is marked on a screen after the radar receives the standard query signal and decodes and identifies the standard query signal.
The prior art identifies beacon signals by: the operator firstly inputs a concerned beacon machine code, and then matches the beacon echo signal aiming at the code, the processing flow can only identify one beacon machine, and the positioning precision of the beacon machine is lower due to the lower data sampling rate, and the higher requirement of positioning navigation can not be met. Therefore, how to improve the identification capability and the positioning accuracy of the beacon is an important problem which needs to be solved urgently at present.
Disclosure of Invention
The invention aims to provide a method for identifying GJB661-89 standard beacon signals aiming at the defects of the prior art, so as to improve the identification capability and the positioning accuracy of a beacon and meet higher requirements of positioning and navigation.
In order to achieve the purpose, the technical scheme of the invention comprises the following steps:
1) carrying out high-frequency sampling on a beacon echo fed back by a GJB661-89 standard beacon machine to obtain digital beacon information, wherein the digital beacon information comprises a radar beam pointing angle theta and beacon pulse position information, and storing the beacon pulse position information into a beacon pulse position information array bcn _ position [ ];
2) inquiring the value of each element in the beacon pulse position information array bcn _ position [ ] one by one from the element 0 of the beacon pulse position information array bcn _ position [ ] and judging whether the value is 0, if so, jumping to 4), and if not, decoding the beacon pulse position bcn _ position [ i ] to obtain new pulse coding information;
3) updating the information in the pulse coding queue according to the new pulse coding information;
4) judging whether the radar beam scans the boundary according to the radar beam pointing angle theta: when the absolute value of theta is less than 55 degrees, the radar beam is not scanned to the boundary, and the step 1) is returned; otherwise, the radar beam is scanned to the boundary, and step 5) is executed;
5) calculating each parameter in each pulse code information of the updated pulse code queue to obtain the beacon machine distance RBCNBeacon orientation thetaBCNAnd beacon encoding BBCNThree items of beacon machine information are stored in a beacon machine information queue, and all information in the pulse coding queue is eliminated;
6) and extracting the distance and the direction of the beacon corresponding to the beacon code concerned by the operator from the beacon information queue, displaying the distance and the direction, and removing all information in the beacon information queue.
Compared with the prior art, the invention has the following advantages:
first, the invention saves all the pulse coding information when updating the information in the pulse coding queue in the processing process, and calculates each parameter in each pulse coding information of the updated pulse coding queue, thereby not only identifying a plurality of beacons simultaneously, but also improving the identification capability of the beacons.
Secondly, the invention carries out high-frequency sampling on the beacon echo fed back by the GJB661-89 standard beacon machine to obtain the digital beacon information, thereby improving the processing precision of the beacon pulse position and the positioning precision of the beacon machine.
Drawings
FIG. 1 is a flow chart of the present invention for decoding and identifying GJB661-89 standard beacon signals.
Detailed Description
Referring to fig. 1, the implementation steps of the invention are as follows:
step 1, acquiring data to be processed.
The specific implementation of this step is as follows:
1a) setting a beacon pulse position information array bcn _ position [ ], and initializing each element of the beacon pulse position information array to 0;
1b) carrying out high-frequency sampling and quantization on the beacon echo fed back by the GJB661-89 standard beacon machine to obtain digital beacon information, wherein the digital beacon information comprises a radar beam pointing angle theta and beacon pulse position information,
1c) the beacon pulse position information is saved into the beacon pulse position information array bcn _ position [ ].
And 2, resolving pulse codes.
The specific implementation of this step is as follows:
2a) inquiring the value of each element in the beacon pulse position information array bcn _ position [ ] one by one from the element 0 of the beacon pulse position information array bcn _ position [ ], and judging whether the value is 0 or not; if the value is 0, skipping to the step 4, and if the value is not 0, executing the step 2b) to decode;
2b) the current position element bcn _ position [ i ] of the beacon pulse position array bcn _ position [ ] is processed according to the decoding mode selected by the operator: when the decoding mode selected by the operator is a double-pulse decoding mode, executing the step 2c), and when the decoding mode selected by the operator is a multi-pulse decoding mode, executing the step 2 d);
2c) querying the beacon pulse position array for the current position element bcn _ position [ i ] and the next beacon pulse position element bcn _ position [ i +1 ];
if bcn _ position [ i +1]]If not 0, the time interval T corresponding to the two pulses is calculatedbcn：
Wherein Sample is the beacon location data sampling rate, different TbcnRepresenting double-pulse codes 1-9, the obtained double-pulse codes are shown in table 1:
TABLE 1 double-pulse coding table
If bcn _ position [ i +1] is 0, no dipulse encoded signal is present, dipulse encoding is 0; then step 2f) is executed;
2d) query beacon pulse position array Current position element bcn _ position [ i]And no more than 5 beacon pulse position elements bcn _ position [ i + n ] behind]N is less than or equal to 5, and the time interval T between the n + i pulse and the i pulse is calculatedbcn：
2e) According to TbcnJudging whether the five predicted pulse positions of the multi-pulse code have pulses or not;
if n is less than or equal to 5, T existsbcn15us, then there is a multi-pulse encoded signal, recording TbcnThe value of n is k when the power is 15 us; setting 4-bit multi-pulse code as binary b 4b 2b 1b 0, initializing the binary b 4b 2b 0 to 0, and calculating the obtained T according to the value of n-1, … and k-1bcnThe binary b 4b 2b 1b 0 is assigned a value of 1, namely:
if Tbcn3us, the binary b0 position is 1;
if Tbcn6us, the binary b1 position is 1;
if Tbcn9us, the binary b2 position is 1;
if Tbcn12us, the binary b4 position is 1, thus obtaining multi-pulse code;
if n is less than or equal to 5, T does not existbcnIn case of 15us, there is no multi-pulse encoded signal, the multi-pulse encoding is 0, and then step 2f) is performed;
2f) the decoded double pulse code or multi-pulse code is used as the parameter of the pulse code B in the new pulse code information, and the current position element bcn _ position [ i ] of the beacon pulse position array bcn _ position [ ] is used as the parameter of the pulse position R in the new pulse code information, thereby obtaining the new pulse code information.
And 3, updating the pulse coding queue.
The specific implementation of this step is as follows:
3a) setting a position error D, and inquiring the existing pulse coding information in a pulse coding queue, wherein the pulse coding information comprises: pulse code B, pulse angle cumulative sum thetasumPulse position R and the number SUM of successful decoding of the pulse code, the position error D is determined according to the speed characteristics of the platform loaded by the equipment, for example: if the conventional speed of a certain platform is Vm/s, the position error D is equal to Sample multiplied by V/150, and the Sample is the beacon position data sampling rate;
3b) comparing the existing pulse code information in the pulse code queue with the new pulse code information obtained from step 2f), if the pulse codes B in the two parameters are the same and the difference R between the pulse positions R in the two parametersdiffIf not, executing step 3c), otherwise, executing step 3D);
3c) changing the existing pulse coding information in the pulse coding queue, namely performing the following operations:
sum of theta at pulse anglesumAdding a current radar beam pointing angle theta on the basis of the current radar beam pointing angle theta;
replacing the existing pulse position R with the current pulse position R;
adding 1 on the basis of the SUM number of successful pulse code decoding;
3d) setting new pulse coding information and adding the new pulse coding information into a pulse coding queue, wherein the parameters of the new pulse coding information are as follows:
the pulse code B is the current code and,
pulse angle cumulative sum thetasumAs the current angle, the angle of the mirror is,
the pulse position R is the current position,
the number SUM of successful decoding of the pulse code is 1.
And 4, judging whether to start to solve the beacon machine information.
Judging whether the radar beam scans the boundary according to the radar beam pointing angle theta: when the absolute value of theta is less than 55 degrees, the radar beam is not scanned to the boundary, and the step 1b) is returned; otherwise, the radar beam is scanned to the boundary, and step 5 is executed.
And 5, resolving beacon machine information.
Calculating each parameter in each pulse code information of the updated pulse code queue to obtain the beacon machine distance RBCNBeacon orientation thetaBCNAnd beacon encoding BBCNAnd storing the three items of beacon machine information into a beacon machine information queue, and then clearing all information in the pulse coding queue.
The steps are specifically realized as follows:
5a) acquiring each pulse code information of the updated pulse code queue, and executing the step 5b) to the step 5d) on each pulse code information;
5b) calculating beacon distance RBCN：
Wherein R is the beacon machine position, c is the light speed, and Sample is the beacon position data sampling rate;
5c) calculating beacon orientation thetaBCN：
Wherein theta issumThe SUM is the pulse angle accumulated SUM, and the SUM is the number of successful decoding of the pulse code;
5d) calculating beacon code BBCN：
BBCN＝B
Wherein B is a pulse code;
5e) the calculated beacon machine distance RBCNBeacon orientation thetaBCNAnd beacon encoding BBCNStoring the three items of beacon machine information into a beacon machine information queue;
5f) all information in the pulse code queue is cleared.
And 6, displaying beacon machine information.
And firstly extracting the distance and the direction of the beacon corresponding to the beacon code concerned by the operator from the beacon information queue, displaying the distance and the direction, then clearing all information in the beacon information queue, and waiting for the generation of the next group of beacon information.
The effects of the present invention can be further illustrated by the following test results:
laboratory tests: the beacon echo fed back by the GJB661-89 standard beacon machine is simulated and generated by using a beacon signal generator, the receiving processing is carried out by using the method of the invention, the positioning error of the beacon machine is not more than 100m, and the identification number of the beacon machine is not less than 5;
external field test: the method is applied to an X-type radar, a GJB661-89 standard beacon machine is arranged at a position 115km away from the X-type radar, and the X-type radar receives and processes the beacon echo fed back by the GJB661-89 standard beacon machine by using the method, so that the beacon machine code can be rapidly identified, and the distance and the direction of the GJB661-89 standard beacon machine can be accurately calculated.
The foregoing description is only an example of the present invention and should not be construed as limiting the invention, as it will be apparent to those skilled in the art that various modifications and variations in form and detail can be made without departing from the principle and structure of the invention after understanding the present disclosure and the principles, but such modifications and variations are considered to be within the scope of the appended claims.
Claims (6)
- The method for identifying the GJB661-89 standard beacon signals comprises the following steps:1) carrying out high-frequency sampling on a beacon echo fed back by a GJB661-89 standard beacon machine to obtain digital beacon information, wherein the digital beacon information comprises a radar beam pointing angle theta and beacon pulse position information, and storing the beacon pulse position information into a beacon pulse position information array bcn _ position [ ];2) inquiring the value of each element in the beacon pulse position information array bcn _ position [ ] one by one from the element 0 of the beacon pulse position information array bcn _ position [ ] and judging whether the value is 0, if the value is 0, jumping to 4), and if the value is not 0, decoding the beacon pulse position element bcn _ position [ i ] to obtain new pulse coding information;3) updating the information in the pulse coding queue according to the new pulse coding information;4) judging whether the radar beam scans the boundary according to the radar beam pointing angle theta: when the absolute value of theta is less than 55 degrees, the radar beam is not scanned to the boundary, and the step 1) is returned; otherwise, the radar beam is scanned to the boundary, and step 5) is executed;5) calculating each parameter in each pulse code information of the updated pulse code queue to obtain the beacon machine distance RBCNBeacon orientation thetaBCNAnd beacon encoding BBCNThree items of beacon machine information are stored in a beacon machine information queue, and all information in the pulse coding queue is eliminated;6) and extracting the distance and the direction of the beacon corresponding to the beacon code concerned by the operator from the beacon information queue, displaying the distance and the direction, and removing all information in the beacon information queue.
- 2. The method according to claim 1, wherein the beacon pulse position bcn _ position [ i ] is decoded in step 2) to obtain new pulse code information, and the decoding is performed according to a decoding mode selected by an operator;when the selected decoding mode is a double-pulse decoding mode, inquiring the current position element bcn _ position [ i ] of the beacon pulse position array]And the subsequent beacon pulse position element bcn _ position [ i +1 [ ]]If bcn _ position [ i +1]]If not 0, the time interval T corresponding to two pulses is calculatedbcn：Wherein Sample is the beacon location data sampling rate, different TbcnRepresents the double-pulse code 1 ~ 9, if bcn _ position [ i + 1%]If the value is 0, no double-pulse coding signal exists, and the double-pulse coding is 0;when the selected decoding mode is a multi-pulse decoding mode, inquiring the beacon pulse position array current position element bcn _ position [ i]And no more than 5 beacon pulse position elements bcn _ position [ i + n ] behind]N is less than or equal to 5, and the time interval T between the n + i pulse and the i pulse is calculatedbcn：If T isbcn15us, then there is a multi-pulse encoded signal, recording TbcnWhen the value of n is k, 15us, 4-bit multi-pulse coding is set as binary b 4b 2b 1b 0; calculating T when n is 1, …, k-1bcnIf T is a value ofbcn3us, bit b0 is binary 1; if T isbcn6us, bit b1 is binary 1; if T isbcn9us, bit b2 is binary 1; if T isbcn12us, the binary b 4bit is 1, thus obtaining multi-pulse code;if n is less than or equal to 5, T does not existbcnIn the case of 15us, the multi-pulse code signal is not present, and the multi-pulse code is 0.
- 3. The method of claim 1, wherein the step 3) of updating the information in the pulse code queue according to the new pulse code information is performed by the steps of:(3a) setting a position error D, and inquiring the existing pulse coding information in a pulse coding queue, wherein the pulse coding information comprises: pulse code B, pulse angle cumulative sum thetasumThe pulse position R and the pulse code decoding success number SUM, wherein the position error is determined according to the platform speed loaded by the equipment;(3b) comparing the existing pulse code information in the pulse code queue with the new pulse code information, if the pulse codes B in the two parameters are the same and the difference R between the pulse positions R in the two parametersdiffNot greater than positionError D, perform (3c), otherwise perform (3D);(3c) changing the existing pulse coding information in the pulse coding queue: i.e. the sum theta is accumulated over the pulse anglesumAdding a current radar beam pointing angle theta on the basis of the current radar beam pointing angle theta; replacing the existing pulse position R with the current pulse position R; adding 1 on the basis of the SUM number of successful pulse code decoding;(3d) adding a new pulse code information in the pulse code queue, wherein the parameters of the new pulse code information are as follows: the pulse code B is the current code, the sum of pulse angles is addedsumAnd the current angle is the pulse position R, and the number SUM of successful pulse code decoding is 1.
- 4. Method according to claim 1, wherein in step 5) a beacon distance R is calculatedBCNCalculated by the following formula:where R is the pulse position, c is the speed of light, and Sample is the beacon position data sampling rate.
- 5. The method of claim 1 wherein the beacon orientation θ is calculated in step 5)BCNCalculated by the following formula:wherein theta issumAnd SUM is the pulse angle accumulated SUM and the number of successful decoding of the pulse code.
- 6. Method according to claim 1, wherein in step 5) a beacon code B is calculatedBCNCalculated by the following formula:BBCN＝Bwherein B is a pulse code.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710353779.9A CN107153174B (en) | 2017-05-18 | 2017-05-18 | GJB661-89 standard beacon signal identification method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710353779.9A CN107153174B (en) | 2017-05-18 | 2017-05-18 | GJB661-89 standard beacon signal identification method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107153174A CN107153174A (en) | 2017-09-12 |
| CN107153174B true CN107153174B (en) | 2019-12-31 |
Family
ID=59793203
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710353779.9A Active CN107153174B (en) | 2017-05-18 | 2017-05-18 | GJB661-89 standard beacon signal identification method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107153174B (en) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109085593B (en) * | 2018-08-08 | 2022-10-25 | 中国航空工业集团公司雷华电子技术研究所 | Meteorological radar beacon detection method based on Kmeans clustering algorithm |
| CN113176451B (en) * | 2021-04-22 | 2023-08-15 | 浙江威星智能仪表股份有限公司 | Pulse signal detection method based on table look-up method |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5367303A (en) * | 1993-09-17 | 1994-11-22 | Alliedsignal Inc. | Parallel observer spatial evaluator |
| EP1640743A2 (en) * | 2004-09-15 | 2006-03-29 | Kabushiki Kaisha Toshiba | Secondary surveillance radar system, and ground system for use therein |
| CN101702024A (en) * | 2009-10-16 | 2010-05-05 | 四川九洲电器集团有限责任公司 | Method and device for processing signals based on broadcast mode automatic dependent surveillance |
| CN102298146A (en) * | 2011-05-18 | 2011-12-28 | 中国航空无线电电子研究所 | Method for decoding secondary radar signal |
| CN106199587A (en) * | 2016-06-27 | 2016-12-07 | 芜湖航飞科技股份有限公司 | A kind of secondary radar response decoding system and device |
-
2017
- 2017-05-18 CN CN201710353779.9A patent/CN107153174B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5367303A (en) * | 1993-09-17 | 1994-11-22 | Alliedsignal Inc. | Parallel observer spatial evaluator |
| EP1640743A2 (en) * | 2004-09-15 | 2006-03-29 | Kabushiki Kaisha Toshiba | Secondary surveillance radar system, and ground system for use therein |
| CN101702024A (en) * | 2009-10-16 | 2010-05-05 | 四川九洲电器集团有限责任公司 | Method and device for processing signals based on broadcast mode automatic dependent surveillance |
| CN102298146A (en) * | 2011-05-18 | 2011-12-28 | 中国航空无线电电子研究所 | Method for decoding secondary radar signal |
| CN106199587A (en) * | 2016-06-27 | 2016-12-07 | 芜湖航飞科技股份有限公司 | A kind of secondary radar response decoding system and device |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107153174A (en) | 2017-09-12 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN1157091C (en) | Selection of location measurement units for determining position of mobile communication station | |
| CN107153174B (en) | GJB661-89 standard beacon signal identification method | |
| EP3896873A1 (en) | Field intensity prediction method and apparatus, and device and storage medium | |
| CN109188443B (en) | Passive target tracking method based on interactive multiple models | |
| CA3057729A1 (en) | A communication method and apparatus | |
| CN106990395B (en) | Radar target calibration method and device | |
| CN113516037A (en) | Marine vessel track segment association method, system, storage medium and equipment | |
| CN103634881B (en) | Frequency searching method and device, terminal | |
| CN111200456B (en) | Fast and low-consumption 3D beam forming method based on joint autonomous positioning | |
| CN111770436B (en) | Indoor Wi-Fi positioning method and device based on CSI and storage medium | |
| EP2627086A1 (en) | Method and device for encoding a block of an image and corresponding reconstructing method and device | |
| CN111160204A (en) | Geological radar image identification method and system based on principal component analysis BP neural network | |
| JP6545369B2 (en) | Lookup table generation method | |
| CN107479044B (en) | Adaptive track starting method based on point track density real-time statistics | |
| CN109270520B (en) | Processing method for acquiring secondary radar response target identity code based on amplitude information | |
| CN114216434A (en) | Target confirmation method, system, equipment and storage medium for mobile measurement and control station | |
| CN109669199B (en) | Positioning method, terminal, computer and storage medium | |
| CN108549061B (en) | Signal clustering method | |
| CN110426696B (en) | Pulse defect radar signal characteristic sequence searching method | |
| CN109490885B (en) | SAR deblurring imaging method based on multi-pulse combination | |
| CN113625249A (en) | Reflector positioning method, robot and computer readable storage medium | |
| CN109143219B (en) | Method for improving TWS radar pitch angle measurement precision | |
| CN111917440A (en) | Wave beam selection and limited feedback method for offshore over-the-horizon large-scale antenna array | |
| CN111901747A (en) | Indoor accurate positioning method and system based on LANDMAC | |
| JP2017067504A (en) | Meteorological radar signal processing equipment, processing program and processing method |
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 |