CN107462873B - Radar interference rapid identification method - Google Patents
Radar interference rapid identification method Download PDFInfo
- Publication number
- CN107462873B CN107462873B CN201710533754.7A CN201710533754A CN107462873B CN 107462873 B CN107462873 B CN 107462873B CN 201710533754 A CN201710533754 A CN 201710533754A CN 107462873 B CN107462873 B CN 107462873B
- Authority
- CN
- China
- Prior art keywords
- interference
- target
- points
- radar
- determining
- 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
Images
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/36—Means for anti-jamming, e.g. ECCM, i.e. electronic counter-counter measures
-
- 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/021—Auxiliary means for detecting or identifying radar signals or the like, e.g. radar jamming signals
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
A radar interference rapid identification method is based on AGC mutation, multi-target detection is carried out, then the type of a signal entering a range gate at present is comprehensively judged according to the relation among the multi-targets or the change of waveform entropy, and conventional interference types are rapidly identified. The method realizes the rapid identification of various interferences by using a method with small calculation amount and without adding a circuit module, provides a basis for subsequent interference countermeasures, reserves more sufficient processing time, is suitable for complex scenes with various interferences such as foil strip interference, noise interference, distance dragging interference, outboard dragging and the like, has low calculation complexity, high real-time performance and quick and effective interference identification, does not need adding a circuit module on the hardware of a radar signal processing system, has simple algorithm and is convenient for application in engineering practice.
Description
Technical Field
The invention relates to a method for quickly identifying radar interference.
Background
The anti-interference measure that the radar can call is determined by the recognition ability of the signal processing system to the interference, and the conventional radar has the ability of resisting various active or passive interferences such as dragging, foil strips, noise and the like. An interference suppression algorithm is usually effective only for one interference type, and due to the fact that interference means are complex and various in actual communication countermeasure, one or more unknown interference types may exist in a received signal. If all the interference countermeasures are added to each target detection blindly, it will undoubtedly cause a great burden to the signal processing circuit, and the real-time performance of the signal processing is seriously affected. Therefore, the interference type can be accurately and quickly identified, a basis is provided for the radar to call a targeted anti-interference measure, and the method has great significance.
At present, in the domestic patent CN105866749A ("a radar interference identification method with synchronous range and speed pull") a broadband radar interference identification method is introduced, which performs interference or target identification by setting a characteristic interference signal error angle and the type of baseband echo signal received by a radar. The method only aims at synchronous dragging of distance and speed, the types of interference identification are less, the input signals of the interference machine need to be set, the prior information of the interference is depended, and the scene adaptability and the engineering realizability are poor.
Patent CN104267379A ("an active and passive radar cooperative anti-interference method based on waveform design") introduces an active and passive radar cooperative anti-interference method, which can estimate interference related parameters by using a passive part through detecting echo signals, so as to design transmitted waveform parameters of an active radar, thereby achieving the purposes of interference identification and countermeasure. Although the method can identify more types of interference, the framework of the method belongs to a multi-mode composite radar, the engineering is complex to realize, and compared with an active radar, an additional set of hardware circuit needs to be added.
Patent CN104678368B ("one-dimensional phase-scanning three-coordinate radar spatial adaptive interference suppression method") introduces a method of extracting interference angle features through a frequency domain, and then forming nulls in the interference direction by an adaptive beam forming technique to suppress interference. The method has high interference identification and countermeasure effect. But the method only has obvious effect on interference mainly based on angle spoofing, and the self-adaptive beam forming algorithm is complex, has higher requirement on the real-time performance of a signal processing system, and can not identify other types of interference.
An interference identification method based on frequency domain sparsity is introduced in an intermittent sampling forwarding interference identification based on frequency domain sparsity, which is disclosed in journal of air force early warning college academic proceedings of 5 th year 2015. According to the method, the intermittent sampling forwarding interference and the echo signals are subjected to Fourier transform, the sparsity difference on the signal spectrum is mined, and the sparsity is used as a characteristic parameter of interference identification, so that the intermittent sampling forwarding interference and the real echo can be effectively distinguished. However, the method is limited in that only the interference-generated mechanism can distinguish the intermittent sampling forwarding interference from the real echo, and cannot identify the non-forwarding interference, and cannot finally obtain the specific type of the interference, so that a basis cannot be provided for a signal processing system to invoke which anti-interference measure.
A master thesis "radar interference signal analysis technology research" of the science and technology university of western's electronics in 2014 introduces a time-frequency analysis method and a high-order statistic method, extracts time-frequency domain characteristic parameters in interference signals, compares the time-frequency domain characteristic parameters with theoretical interference signal model parameters, and provides theoretical basic parameter information for identification of interference patterns. The method can accurately identify various interference types, but the wavelet analysis method and the bispectrum analysis method of the third-order statistics have huge calculation amount, greatly increase the burden of a signal processing circuit, occupy a large amount of processing time, and seriously influence the implementation of subsequent anti-interference measures, so the engineering realization significance is not great.
Disclosure of Invention
The invention provides a method for quickly identifying radar interference, which does not need to add an additional processing circuit, has small calculation amount, low calculation complexity and high real-time performance, can simultaneously identify more interference types, and is quick and effective in interference identification.
In order to achieve the above object, the present invention provides a method for quickly identifying radar interference, comprising the following steps:
step S1, in the routine detection of the tracked target of each frame of radar echo signal, carrying out mutation judgment on the AGC value of the current frame, and if the AGC value has mutation, carrying out step S2;
step S2, performing multi-target detection on the current two-dimensional detection plane after the tracking points and the protection area thereof are removed, if multi-target is detected, recording distance and frequency dimension coordinates of the multi-target points, and performing step S3, if no target is detected, performing step S4;
step S3, performing Doppler dimension diffusion detection on the detected multiple target points, if Doppler expansion phenomenon occurs, determining that a target with higher polymerization degree and wider frequency dimension distribution appears on the plane at the moment, conforming to the basic characteristics of foil strip interference, confirming that the interference type is foil strip interference, and if expansion phenomenon does not occur, performing step S5;
step S4, performing waveform entropy mutation judgment on the current two-dimensional detection plane, if the waveform entropy has mutation, determining that the plane noise is not pure Gaussian white noise at the moment, the radar is subjected to aiming type or frequency sweeping type noise suppression interference, determining that the interference type is noise interference, and if the waveform entropy has no mutation, determining that other situations beyond the range of the identification method occur, and determining that the situation is non-interference;
and step S5, comparing the detected relative distances between the multiple target points and the tracking points one by one, if the relative distances are stable, considering that even if other frequency dimension non-extended targets except the target point exist, the relative distances between the target and the original tracking points are kept unchanged, conforming to a plurality of angle inverses of target tail towing or ship-borne active false target jammers, confirming that the jamming type is outboard towing jamming, and if the relative distances are changed, belonging to the basic characteristic of typical distance towing jamming, and confirming that the jamming type is distance towing jamming.
In step S1, if the AGC value is increased or decreased from the original value by an amount greater than the threshold value, it is determined that the AGC value is suddenly changed.
In step S2, the protection area of the tracking point includes: the distance around the tracking point is n points and the frequency is m points, where the value of n is determined according to the volume and resolution of the object and the value of m is determined according to the motion velocity and resolution of the object.
In step S3, the distance dimension and the frequency dimension of the multiple target points are compared with each other, and if there are 3 points in a circular area with a radius of 1 or a square area with a side length of l, it is determined that the doppler spread phenomenon occurs.
In step S4, the waveform entropy of the current frame is compared with the waveform entropy of the previous frame, and if the change is greater than the threshold value, the waveform entropy jump is considered to occur.
The method is based on AGC mutation, carries out multi-target detection, comprehensively judges the type of the current signal entering the range gate according to the relation between the multiple targets or the change of the waveform entropy, realizes the rapid identification of the conventional interference type, realizes the rapid identification of multiple interferences by using a method with smaller calculation amount and without adding a circuit module, provides a basis for subsequent interference countermeasures, reserves more sufficient processing time, is suitable for complex scenes with multiple interferences such as foil strip interference, noise interference, distance dragging interference, outboard dragging and the like, is low in calculation complexity, high in instantaneity, rapid and effective in interference identification, does not need to add a circuit module on the hardware of a radar signal processing system, is simple in algorithm, and is convenient to apply in engineering practice.
Drawings
Fig. 1 is a flowchart of a method for quickly identifying radar interference according to the present invention.
Detailed Description
The preferred embodiment of the present invention is described in detail below with reference to fig. 1.
As shown in fig. 1, the present invention provides a method for quickly identifying radar interference, which includes the following steps:
step S1, in the routine detection of the tracked target of each frame of radar echo signal, carrying out mutation judgment on the AGC value of the current frame, and if the AGC value has mutation, carrying out step S2;
step S2, performing multi-target detection on the current two-dimensional detection plane after the tracking points and the protection area thereof are removed, if multi-target is detected, recording distance and frequency dimension coordinates of the multi-target points, and performing step S3, if no target is detected, performing step S4;
step S3, performing Doppler dimension diffusion detection on the detected multiple target points, if Doppler expansion phenomenon occurs, determining that a target with higher polymerization degree and wider frequency dimension distribution appears on the plane at the moment, conforming to the basic characteristics of foil strip interference, confirming that the interference type is foil strip interference, and if expansion phenomenon does not occur, performing step S5;
step S4, performing waveform entropy mutation judgment on the current two-dimensional detection plane, if the waveform entropy has mutation, determining that the plane noise is not pure Gaussian white noise at the moment, the radar is subjected to aiming type or frequency sweeping type noise suppression interference, determining that the interference type is noise interference, and if the waveform entropy has no mutation, determining that other situations beyond the range of the identification method occur, and determining that the situation is non-interference;
and step S5, comparing the detected relative distances between the multiple target points and the tracking points one by one, if the relative distances are stable, considering that even if other frequency dimension non-extended targets except the target point exist, the relative distances between the target and the original tracking points are kept unchanged, conforming to a plurality of angle inverses of target tail towing or ship-borne active false target jammers, confirming that the jamming type is outboard towing jamming, and if the relative distances are changed, belonging to the basic characteristic of typical distance towing jamming, and confirming that the jamming type is distance towing jamming.
The AGC value mainly reflects the intensity of the signal energy entering the receiver at present, and if the AGC changes suddenly, the signal energy entering the receiver at the moment is considered to generate obvious changes. It is generally believed that the interfering signal should interfere with the real echo quickly and efficiently, and the signal energy should be significantly greater than the real echo. In step S1, the abrupt change of the AGC value includes a rising edge of the AGC value and a falling edge of the AGC value; increasing the AGC from an original value, wherein the increasing amount is larger than a threshold value, and considering that the AGC has a rising edge; reducing the value from the original value, wherein the reduction is larger than a threshold value, and considering that the AGC has a falling edge; and the jumping is considered to occur when the AGC is increased or reduced from the original value by more than a threshold value. Depending on the AGC voltage definition of different receivers, AGC up may reflect an increase in signal energy, and a decrease may also be defined as an increase in signal energy. Here, the AGC mutation is only used to indicate that the energy entering the receiver becomes stronger, and it is considered that there is a possibility of an interference signal entering.
In step S2, the protection area of the tracking point includes: the distance around the tracking point is n (determined by the target volume and resolution, assuming that the target length is 300 m and the resolution is 15 m, n is 5) points and the frequency is m (determined by the target moving speed and resolution, assuming that the target moving speed is 122 m/s and the resolution is 15 m, m is 5) points. The tracking points and the protection area thereof are due to the characteristics of the target body, so that a plurality of scattering points appear, and the tracking points and the protection area thereof are not considered as multiple targets.
The multi-target detection adopts a conventional target detection method, such as constant false alarm and the like.
In step S3, the distance dimension and the frequency dimension of the multiple target points are compared with each other, and if there are 3 points in a circular area with a radius of 1 or a square area with a side length of l, it is determined that the doppler spread phenomenon occurs.
In step S4, the waveform entropy reflects the noise level of the entire detection plane. And comparing the waveform entropy of the current frame with the waveform entropy of the previous frame, and if the change is larger than a threshold value, determining that the waveform entropy mutation occurs.
The method is based on AGC mutation, carries out multi-target detection, comprehensively judges the type of the current signal entering the range gate according to the relation between the multiple targets or the change of the waveform entropy, realizes the rapid identification of the conventional interference type, realizes the rapid identification of multiple interferences by using a method with smaller calculation amount and without adding a circuit module, provides a basis for subsequent interference countermeasures, reserves more sufficient processing time, is suitable for complex scenes with multiple interferences such as foil strip interference, noise interference, distance dragging interference, outboard dragging and the like, is low in calculation complexity, high in instantaneity, rapid and effective in interference identification, does not need to add a circuit module on the hardware of a radar signal processing system, is simple in algorithm, and is convenient to apply in engineering practice.
While the present invention has been described in detail with reference to the preferred embodiments, it should be understood that the above description should not be taken as limiting the invention. Various modifications and alterations to this invention will become apparent to those skilled in the art upon reading the foregoing description. Accordingly, the scope of the invention should be determined from the following claims.
Claims (4)
1. A radar interference rapid identification method is characterized by comprising the following steps:
step S1, in the routine detection of the tracked target of each frame of radar echo signal, carrying out mutation judgment on the AGC value of the current frame, and if the AGC value has mutation, carrying out step S2;
step S2, performing multi-target detection on the current two-dimensional detection plane after the tracking points and the protection area thereof are removed, if multi-target is detected, recording distance and frequency dimension coordinates of the multi-target points, and performing step S3, if no target is detected, performing step S4;
step S3, performing Doppler dimension diffusion detection on the detected multiple target points, if a Doppler expansion phenomenon occurs, determining that a target with higher polymerization degree and wider frequency dimension distribution appears on the plane at the moment, conforming to the basic characteristics of foil strip interference, confirming that the interference type is foil strip interference, and if the Doppler expansion phenomenon does not occur, performing step S5; comparing the distance dimension and the frequency dimension of the multiple target points, and if 3 points are simultaneously positioned in a circular area with the radius of 1 or a square area with the side length of l, determining that the Doppler expansion phenomenon occurs;
step S4, performing waveform entropy mutation judgment on the current two-dimensional detection plane, if the waveform entropy has mutation, determining that the plane noise is not pure Gaussian white noise at the moment, the radar is subjected to aiming type or frequency sweeping type noise suppression interference, determining that the interference type is noise interference, and if the waveform entropy has no mutation, determining that other situations beyond the range of the identification method occur, and determining that the situation is non-interference;
and step S5, comparing the detected relative distances between the multiple target points and the tracking points one by one, if the relative distances are stable, considering that even if other frequency dimension non-extended targets except the target point exist, the relative distances between the target and the original tracking points are kept unchanged, conforming to a plurality of angle inverses of target tail towing or ship-borne active false target jammers, confirming that the jamming type is outboard towing jamming, and if the relative distances are changed, belonging to the basic characteristic of typical distance towing jamming, and confirming that the jamming type is distance towing jamming.
2. The method for rapidly identifying radar interference according to claim 1, wherein in step S1, if the AGC value is increased or decreased from the original value by an amount greater than a threshold value, the AGC value is considered to have a sudden change.
3. The method for fast detecting radar interference according to claim 1, wherein in step S2, the protection area of the tracking point includes: the distance around the tracking point is n points and the frequency is m points, where the value of n is determined according to the volume and resolution of the object and the value of m is determined according to the motion velocity and resolution of the object.
4. The method for rapidly identifying radar interference according to claim 1, wherein in step S4, the waveform entropy of the current frame is compared with the waveform entropy of the previous frame, and if the change is greater than a threshold value, it is determined that the sudden change of the waveform entropy has occurred.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710533754.7A CN107462873B (en) | 2017-07-03 | 2017-07-03 | Radar interference rapid identification method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710533754.7A CN107462873B (en) | 2017-07-03 | 2017-07-03 | Radar interference rapid identification method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107462873A CN107462873A (en) | 2017-12-12 |
| CN107462873B true CN107462873B (en) | 2020-06-23 |
Family
ID=60546548
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710533754.7A Active CN107462873B (en) | 2017-07-03 | 2017-07-03 | Radar interference rapid identification method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107462873B (en) |
Families Citing this family (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108919203B (en) * | 2018-05-21 | 2022-04-22 | 西安电子科技大学 | Radar active interference identification method and system |
| US11585919B2 (en) * | 2018-07-19 | 2023-02-21 | Qualcomm Incorporated | Multi-radar coexistence using slow rate interference identification and suppression |
| CN109633558A (en) * | 2018-10-25 | 2019-04-16 | 上海无线电设备研究所 | A kind of DOA estimation algorithm based on polarization time-frequency distributions |
| CN109521433A (en) * | 2018-10-31 | 2019-03-26 | 歌尔股份有限公司 | Determination method, the processing method and processing device of the abnormal frame point cloud data of laser radar |
| CN109283507B (en) * | 2018-12-03 | 2023-03-21 | 北京遥感设备研究所 | Radar target identification method and system based on time-frequency domain characteristics |
| CN110531327B (en) * | 2019-10-11 | 2021-08-31 | 上海无线电设备研究所 | Radar anti-interference detection method based on automatic gain control |
| CN112859116B (en) * | 2020-12-28 | 2022-11-15 | 杭州电子科技大学 | Single-antenna GNSS deception jamming detection method based on main feature blind source separation |
| CN112731309B (en) * | 2021-01-06 | 2022-09-02 | 哈尔滨工程大学 | Active interference identification method based on bilinear efficient neural network |
| CN113759324B (en) * | 2021-07-21 | 2024-01-30 | 西安电子科技大学 | Countermeasure method based on foil interference actual measurement data |
| CN114660562A (en) * | 2022-03-17 | 2022-06-24 | 成都美数科技有限公司 | Adaptive filtering method and filter |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5483240A (en) * | 1994-09-28 | 1996-01-09 | Rockwell International Corporation | Radar terrain bounce jamming detection using ground clutter tracking |
| JP2005331259A (en) * | 2004-05-18 | 2005-12-02 | Mitsubishi Electric Corp | Radar jamming method |
| CN102707266A (en) * | 2012-05-24 | 2012-10-03 | 北京理工大学 | Radar with anti-interference and multi-target identification functions and detection method thereof |
| CN104215938A (en) * | 2014-09-30 | 2014-12-17 | 中国人民解放军63892部队 | Polarization detection and recognition method and device for terrain rebound interference |
| CN106019248A (en) * | 2016-06-30 | 2016-10-12 | 北京润科通用技术有限公司 | Chaff-jammed echo simulation method and system |
| CN106371082A (en) * | 2016-08-26 | 2017-02-01 | 上海无线电设备研究所 | Linear velocity pull-off jamming identification method |
-
2017
- 2017-07-03 CN CN201710533754.7A patent/CN107462873B/en active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5483240A (en) * | 1994-09-28 | 1996-01-09 | Rockwell International Corporation | Radar terrain bounce jamming detection using ground clutter tracking |
| JP2005331259A (en) * | 2004-05-18 | 2005-12-02 | Mitsubishi Electric Corp | Radar jamming method |
| CN102707266A (en) * | 2012-05-24 | 2012-10-03 | 北京理工大学 | Radar with anti-interference and multi-target identification functions and detection method thereof |
| CN104215938A (en) * | 2014-09-30 | 2014-12-17 | 中国人民解放军63892部队 | Polarization detection and recognition method and device for terrain rebound interference |
| CN106019248A (en) * | 2016-06-30 | 2016-10-12 | 北京润科通用技术有限公司 | Chaff-jammed echo simulation method and system |
| CN106371082A (en) * | 2016-08-26 | 2017-02-01 | 上海无线电设备研究所 | Linear velocity pull-off jamming identification method |
Non-Patent Citations (2)
| Title |
|---|
| 脉冲多普勒雷达特征分析及干扰识别方法研究;户锋刚等;《舰船电子对抗》;20070630;第30卷;第30-32、36页 * |
| 识别箔条干扰的一种实用方法;李伟等;《现代雷达》;20000630(第3期);第35-38、43页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107462873A (en) | 2017-12-12 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107462873B (en) | Radar interference rapid identification method | |
| US9157992B2 (en) | Knowledge aided detector | |
| US10690748B2 (en) | System and method for interference detection in a RF receiver | |
| CN108627807B (en) | Anti-interference method for airborne radar | |
| CN108120976B (en) | Ground clutter spectrum leakage suppression method based on Doppler channel characteristics | |
| CN107607916B (en) | Self-defense type speed and distance joint deception jamming resisting method | |
| CN112014806B (en) | Unintentional interference suppression method for airborne radar under complex interference scene | |
| CN107064903B (en) | GLRT-based target detection method under multiple heterogeneous satellites | |
| CN111948618B (en) | Forward scattering target detection method and system based on satellite external radiation source | |
| CN111896928B (en) | Multi-station radar target detection method based on active deception interference elimination | |
| CN105891799B (en) | Suitable for the active interference detection method of mechanical scanning radar | |
| CN113567934B (en) | Radar detection method based on cooperation | |
| CN109031221B (en) | Data layer multi-dimensional feature extraction-based radar dense false target interference resisting method | |
| CN109143184B (en) | Double-threshold detection method for scanning radar | |
| CN109633624B (en) | RGPO interference identification method based on filtering data processing | |
| CN109655794B (en) | Detection and identification method for suppressing interference by narrow-band self-defense noise | |
| Raja Abdullah et al. | Improvement in detection with forward scattering radar | |
| CN106019250B (en) | Based on angle scintillations relay type decoy discrimination method | |
| CN112285656B (en) | Time-sensitive target deception jamming prevention method based on track characteristics | |
| CN108508413B (en) | Target detection method based on probability statistics under low signal-to-noise ratio condition | |
| KR102421093B1 (en) | Apparatus for Cancelling Interference Plot based on Signal Processing Gain, and Method thereof | |
| CN104898097A (en) | FPGA-based phase demodulation constant false alarm rate (CFAR) radar signal detection method | |
| Novillo et al. | Performance evaluation of radar systems in noise jamming environments | |
| CN116953655B (en) | Forward scattering radar moving target detection method and device based on Durbin test | |
| Abdullah et al. | Target prediction in forward scattering radar |
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 |