CN117892072B - Signal detection method based on Gaussian rank correlation coefficient - Google Patents
Signal detection method based on Gaussian rank correlation coefficient Download PDFInfo
- Publication number
- CN117892072B CN117892072B CN202410051774.0A CN202410051774A CN117892072B CN 117892072 B CN117892072 B CN 117892072B CN 202410051774 A CN202410051774 A CN 202410051774A CN 117892072 B CN117892072 B CN 117892072B
- Authority
- CN
- China
- Prior art keywords
- rank
- gaussian
- signal
- statistic
- echo signals
- 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
- 238000001514 detection method Methods 0.000 title claims abstract description 33
- 238000006243 chemical reaction Methods 0.000 claims abstract description 12
- 230000001186 cumulative effect Effects 0.000 claims description 7
- 238000005315 distribution function Methods 0.000 claims description 7
- 230000009466 transformation Effects 0.000 claims description 4
- 238000005070 sampling Methods 0.000 claims description 3
- 238000010276 construction Methods 0.000 claims description 2
- 238000000034 method Methods 0.000 abstract description 6
- 238000012360 testing method Methods 0.000 description 3
- 238000000342 Monte Carlo simulation Methods 0.000 description 2
- 230000002159 abnormal effect Effects 0.000 description 2
- 239000000654 additive Substances 0.000 description 2
- 230000000996 additive effect Effects 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 238000005314 correlation function Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000004891 communication Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000006866 deterioration Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 238000010295 mobile communication Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000001629 suppression Effects 0.000 description 1
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/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
-
- 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/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
- G01S7/539—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/10—Pre-processing; Data cleansing
- G06F18/15—Statistical pre-processing, e.g. techniques for normalisation or restoring missing data
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L25/00—Baseband systems
- H04L25/02—Details ; arrangements for supplying electrical power along data transmission lines
- H04L25/03—Shaping networks in transmitter or receiver, e.g. adaptive shaping networks
- H04L25/03006—Arrangements for removing intersymbol interference
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L25/00—Baseband systems
- H04L25/02—Details ; arrangements for supplying electrical power along data transmission lines
- H04L25/03—Shaping networks in transmitter or receiver, e.g. adaptive shaping networks
- H04L25/03828—Arrangements for spectral shaping; Arrangements for providing signals with specified spectral properties
- H04L25/03834—Arrangements for spectral shaping; Arrangements for providing signals with specified spectral properties using pulse shaping
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Computer Networks & Wireless Communication (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Power Engineering (AREA)
- Remote Sensing (AREA)
- Probability & Statistics with Applications (AREA)
- Evolutionary Biology (AREA)
- Bioinformatics & Computational Biology (AREA)
- Mathematical Physics (AREA)
- Radar, Positioning & Navigation (AREA)
- Signal Processing (AREA)
- General Engineering & Computer Science (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Algebra (AREA)
- Evolutionary Computation (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Artificial Intelligence (AREA)
- Software Systems (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Databases & Information Systems (AREA)
- Operations Research (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
The invention provides a signal detection method based on Gaussian rank correlation coefficients, which comprises the following steps: acquiring echo signals received by two paths of receivers; calculating rank statistics of the echo signals, and performing Gaussian rank conversion on the echo signals; constructing statistics based on the Gaussian rank transformed echo signals; calculating a corresponding threshold value according to the given false alarm probability; and comparing the statistic with the threshold value to finish the detection process. The method has the advantages of simple structure, easy engineering realization and strong robustness, and can be used as a powerful tool for signal detection under the background noise containing pulse components.
Description
Technical Field
The invention belongs to the technical field of communication, and particularly relates to a signal detection method based on a Gaussian rank correlation coefficient.
Background
The signal detection is an important signal analysis technology and is widely applied to the fields of phased array radars, underwater sonars, mobile communication and the like. The purpose is to determine whether and when a useful signal is present, and facilitate the subsequent extraction of more useful high quality signals by the system through techniques such as parameter estimation, code recognition, modulation analysis, and the like.
Currently, matched filters are classic signal detection algorithms in the literature. It is popular because of its simple form, complete theory, easy implementation. The matched filter is a linear correlator that maximizes the signal-to-noise ratio of the signal output with optimal performance in an additive white gaussian noise background. However, due to factors such as atmospheric radiation and channel interference, noise sources in actual scenes may also generate some noise following non-gaussian distribution in addition to additive gaussian noise. These non-gaussian distributed noise features a short duration but a large amplitude, commonly referred to as impulse noise. The signal-to-noise ratio of the received echo signal is rapidly reduced over a period of time due to the presence of impulse noise, resulting in rapid deterioration or even failure of the matched filter.
Disclosure of Invention
Aiming at the problems, the invention uses Gaussian rank transformation of signals to realize suppression of large abnormal values in impulse noise, and then calculates correlation functions of the transformed signals to design statistics, thereby realizing high-precision signal detection in impulse noise environment.
The technical scheme of the invention is as follows:
A signal detection method based on Gaussian rank correlation coefficient comprises the following steps:
acquiring echo signals received by two paths of receivers;
Respectively calculating rank statistics of the two paths of echo signals, and carrying out Gaussian rank conversion on the two paths of echo signals;
constructing statistics based on the two paths of echo signals after Gaussian rank conversion;
Calculating a corresponding threshold value according to the given false alarm probability;
Comparing the statistic with the threshold value to determine the detection result of the source signal.
Further, the two echo signals specifically are:
x(i)=θ1s(i)+z1(i)
y(i)=θ2s(i)+z2(i)
i=1,2…,n
Wherein { x (i), y (i) } is an echo signal received by two paths of receivers at a sampling time i, s (i) is a source signal, θ l is less than or equal to 1, z l (i) is background noise, l=1, 2 is a first path of receiver, and n is a signal length.
Further, the formula for calculating the rank statistics of the two echo signals respectively is as follows:
i=1,2…,n
Wherein R i represents the rank statistic of the first echo signal x (i), and Q i represents the rank statistic of the second echo signal y (i);
wherein H (·) represents heaviside step functions.
Further, the gaussian rank conversion for the two echo signals is specifically:
i=1,2…,n
Wherein Φ -1 (·) represents the inverse of the standard normal distribution cumulative distribution function Φ.
Further, the construction statistics based on the two paths of echo signals after Gaussian rank transformation is specifically as follows:
Wherein the coefficient is
Wherein the coefficient isRepresenting the variance value of the signal after gaussian rank transformation, Φ -1 (·) represents the inverse of the standard normal distribution cumulative distribution function Φ.
Further, the calculating the corresponding threshold according to the given false alarm probability specifically includes:
λ=σΦ-1(1-Pf)
Where Φ -1 (·) represents the inverse of the standard normal distribution cumulative distribution function Φ, σ represents the variance of the statistics, and P f is the given false alarm probability.
Further, the comparing statistics and the threshold value are specifically:
Comparing the statistic with the threshold value, and if the statistic is larger than the threshold value, the source signal exists; and if the statistic is smaller than or equal to the threshold value, the source signal is considered to be absent.
The invention has the technical effects that:
The invention fully suppresses the negative influence of a large abnormal value in impulse noise by carrying out Gaussian rank conversion on the signals, and shows excellent robustness in the impulse noise environment. Therefore, under the environment noise containing pulse components, the Gaussian rank correlation coefficient is an effective tool for signal detection, and the performance is quite good.
Drawings
The accompanying drawings illustrate various embodiments by way of example in general and not by way of limitation, and together with the description and claims serve to explain the inventive embodiments. Wherever possible, the same reference numbers will be used throughout the drawings to refer to the same or like parts. Such embodiments are illustrative and not intended to be exhaustive or exclusive of the present apparatus or method.
FIG. 1 is a schematic flow chart of a signal detection method of the present invention;
FIG. 2 shows a schematic diagram of the detection structure of the detection method of the present invention;
fig. 3 shows a comparison graph of detection probabilities of signal detection in a pulse noise environment for a matched filter and a gaussian rank correlation detector.
Detailed Description
It should be noted that, without conflict, the embodiments of the present application and features of the embodiments may be combined with each other. The application will be described in detail below with reference to the drawings in connection with embodiments.
The embodiment provides a signal detection method based on a gaussian rank correlation coefficient, as shown in fig. 1, comprising the following steps:
step S10: collecting echo signals received by two paths of receivers;
Step S20: calculating rank statistics of two paths of received signals, and acquiring corresponding nonlinear Gaussian rank conversion;
step S30: calculating correlation functions of the two paths of transformed received signals, and constructing test statistics T;
step S40: calculating a corresponding threshold lambda according to the given false alarm probability P f;
Step S50: comparing the statistic with the threshold value, and if the statistic T is larger than the threshold value lambda, the source signal exists. Otherwise, the source signal is considered to be absent, thereby completing the detection process.
As a preferred implementation manner of this embodiment, the echo signals received by the two-way receiver at the ith sampling time have the following form:
x(i)=θ1s(i)+z1(i)
y(i)=θ2s(i)+z2(i),i=1,2…,n
where { x (i), y (i) } is an echo signal received by two paths of receivers, s (i) is a source signal, θ l +.1 (l=1, 2) represents a signal attenuation factor, z l (i) (l=1, 2) represents background noise, and n is a signal length. Our goal is to determine whether a signal s is present by measuring the degree of correlation of two received signals.
As a preferred implementation manner of this embodiment, calculating rank statistics of two paths of received signals may be implemented by the following formula:
wherein H (·) represents heaviside step functions.
As a preferred implementation manner of this embodiment, the gaussian rank transform of the two-path received signal is:
As a preferred implementation of this example, the test statistic is constructed by:
Wherein the coefficient is
As a preferred implementation manner of this embodiment, the formula for calculating the threshold value by the false alarm probability is:
λ=σΦ-1(1-Pf)
Where Φ -1 (·) represents the inverse of the standard normal distribution cumulative distribution function Φ, σ represents the variance of the statistic, which is defined as:
as a preferred implementation of this embodiment, the signal detection process is implemented by the following decision:
Theoretical basis for calculating threshold value.
According to the statistic expression and in combination with the definition of the rank statistic, the average value of the statistic can be obtained as follows:
when n is an even number, it is possible to obtain
And when n is an odd number
To sum up, the mean value of the statistics isAccording to the relation between the variance and the mean, the statistic variance has the following calculation formula:
It has been demonstrated in the literature that for independent co-distributed binary samples { { x (1), y (1) },.+ -. X (n), y (n) }, the corresponding rank statistic is { (R 1,Q1),...,(Rn,Qn) }, assuming J, K, L, M are real valued functions
The expression of the statistic variance σ 2 can be found by combining the above equations.
When the sample length n is sufficiently large, the false alarm probability P f is expressed by:
at this time, a calculation formula of the threshold value can be obtained by inversion.
Example 1
In order to analyze the performance of the Gaussian rank correlation detector and the matched filter in signal detection under pulse noise, the invention is verified by a Monte Carlo experiment. Fig. 2 is a schematic diagram of a detection structure of a gaussian rank correlation detector, where x (1), x (2) …, x (n) and y (1), y (2) …, y (n) are echo signals received by two receivers, R 1,R2…,Rn and Q 1,Q2…,Qn are rank statistics of two received signals, x '(1), x' (2) …, x '(n) and y' (1), y '(2) …, y' (n) are signals obtained by performing gaussian rank conversion on two received signals, t is a statistic obtained by performing correlation operation on the two converted signals, and finally comparing the statistic t with a threshold lambda to obtain a binary hypothesis test result (i.e. 0: the source signal is absent or i 1: the source signal is present).
The experimental parameters were set as follows:
The source signal is randomly generated from a length n=200 of signal subject to a standard normal distribution.
Impulse noise is modeled by a mixed gaussian distribution:
Where ε=0.01 denotes the probability that the pulse component occurs throughout the impulse noise environment and δ 2=100>>δ1 denotes the standard deviation of the pulse component. The signal-to-noise ratio of the received signal at this time can be defined as:
Through Monte Carlo experiments, performance performances of the Gaussian rank correlation detector and the matched filter under different signal to noise ratios are compared and analyzed, and the Gaussian rank correlation detector can be verified to have robustness under the impulse noise environment. The number of experiments is 10 4, the false alarm probability is P f =0.1, and the signal attenuation factor is θ 1=θ2 =0.5. The experimental results are shown in FIG. 3.
As shown in the experimental result of fig. 3, the detection probability curve of the matched filter is close to a horizontal line with P f =0.1 due to the existence of the pulse component, so that the detection effect is completely lost, and the gaussian rank correlation detector has higher detection probability and shows robustness to the impulse noise, so that the gaussian rank correlation coefficient can be used as a powerful tool for signal detection in the impulse noise environment.
The foregoing is only a preferred embodiment of the present invention, but the scope of the present invention is not limited thereto, and any person skilled in the art, who is within the scope of the present invention, should make equivalent substitutions or modifications according to the technical solution of the present invention and the inventive concept thereof, and should be covered by the scope of the present invention.
Claims (6)
1. A signal detection method based on Gaussian rank correlation coefficients is characterized by comprising the following steps:
acquiring echo signals received by two paths of receivers;
Respectively calculating rank statistics of the two paths of echo signals, and carrying out Gaussian rank conversion on the two paths of echo signals;
constructing statistics based on the two paths of echo signals after Gaussian rank conversion;
Calculating a corresponding threshold value according to the given false alarm probability;
Comparing the statistic with the threshold value to determine the detection result of the source signal;
the construction statistics of the two paths of echo signals based on Gaussian rank transformation is specifically as follows:
Wherein the coefficient is
Wherein the coefficient isThe variance value of the signal after the gaussian rank conversion is represented, Φ -1 (·) represents an inverse function of a standard normal distribution cumulative distribution function Φ, n is a signal length, R i represents rank statistics of the first echo signal x (i), Q i represents rank statistics of the second echo signal y (i), x '(i) is a signal after the first echo signal x (i) is subjected to the gaussian rank conversion, and y' (i) is a signal after the second echo signal y (i) is subjected to the gaussian rank conversion.
2. The signal detection method based on a gaussian rank correlation coefficient according to claim 1, wherein the two echo signals are specifically:
x(i)=θ1s(i)+z1(i)
y(i)=θ2s(i)+z2(i)
i=1,2…,n
Wherein { x (i), y (i) } is an echo signal received by two paths of receivers at a sampling time i, s (i) is a source signal, θ l is less than or equal to 1, z l (i) is background noise, l=1, 2 is a first path of receiver, and n is a signal length.
3. The signal detection method based on a gaussian rank correlation coefficient according to claim 1, wherein the formula for calculating rank statistics of two echo signals respectively is:
i=1,2…,n
Wherein R i represents the rank statistic of the first echo signal x (i), and Q i represents the rank statistic of the second echo signal y (i);
wherein H (·) represents heaviside step functions.
4. The signal detection method based on a gaussian rank correlation coefficient according to claim 1, wherein performing a gaussian rank transform on two echo signals is specifically:
i=1,2…,n
Wherein Φ -1 (·) represents the inverse of the standard normal distribution cumulative distribution function Φ.
5. The signal detection method based on a gaussian rank correlation coefficient according to claim 1, wherein said calculating a corresponding threshold value according to a given false alarm probability is specifically:
λ=σΦ-1(1-Pf)
Where Φ -1 (·) represents the inverse of the standard normal distribution cumulative distribution function Φ, σ represents the variance of the statistics, and P f is the given false alarm probability.
6. The signal detection method based on a gaussian rank correlation coefficient according to claim 1, wherein the magnitude of the comparison statistic and the threshold value is specifically:
Comparing the statistic with the threshold value, and if the statistic is larger than the threshold value, the source signal exists; and if the statistic is smaller than or equal to the threshold value, the source signal is considered to be absent.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202410051774.0A CN117892072B (en) | 2024-01-12 | 2024-01-12 | Signal detection method based on Gaussian rank correlation coefficient |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202410051774.0A CN117892072B (en) | 2024-01-12 | 2024-01-12 | Signal detection method based on Gaussian rank correlation coefficient |
Publications (2)
Publication Number | Publication Date |
---|---|
CN117892072A CN117892072A (en) | 2024-04-16 |
CN117892072B true CN117892072B (en) | 2024-06-07 |
Family
ID=90648601
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202410051774.0A Active CN117892072B (en) | 2024-01-12 | 2024-01-12 | Signal detection method based on Gaussian rank correlation coefficient |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN117892072B (en) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20120032232A (en) * | 2010-09-28 | 2012-04-05 | 동아대학교 산학협력단 | Signal detection device on generalized normal-laplace distributed noise environments |
CN105761216A (en) * | 2016-01-25 | 2016-07-13 | 西北大学 | Image de-noising processing method and device |
CN110568415A (en) * | 2019-07-22 | 2019-12-13 | 广东工业大学 | signal detection method based on Arctan function under Gaussian mixture model |
-
2024
- 2024-01-12 CN CN202410051774.0A patent/CN117892072B/en active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20120032232A (en) * | 2010-09-28 | 2012-04-05 | 동아대학교 산학협력단 | Signal detection device on generalized normal-laplace distributed noise environments |
CN105761216A (en) * | 2016-01-25 | 2016-07-13 | 西北大学 | Image de-noising processing method and device |
CN110568415A (en) * | 2019-07-22 | 2019-12-13 | 广东工业大学 | signal detection method based on Arctan function under Gaussian mixture model |
Also Published As
Publication number | Publication date |
---|---|
CN117892072A (en) | 2024-04-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Conte et al. | Recursive estimation of the covariance matrix of a compound-Gaussian process and its application to adaptive CFAR detection | |
Flam et al. | On MMSE estimation: A linear model under Gaussian mixture statistics | |
Cui et al. | Multiple-input multiple-output radar detectors design in non-Gaussian clutter | |
Tian et al. | DOA estimation for CD sources by complex cyclic correntropy in an impulsive noise environment | |
CN110568415B (en) | Signal detection method based on Arctan function under Gaussian mixture model | |
CN117892072B (en) | Signal detection method based on Gaussian rank correlation coefficient | |
Roenko et al. | An overview of the adaptive robust DFT | |
Wang et al. | Hopping time estimation of frequency-hopping signals based on HMM-enhanced Bayesian compressive sensing with missing observations | |
CN111090089A (en) | Space-time adaptive detection method based on two types of auxiliary data | |
Chen et al. | New approximate distributions for the generalized likelihood ratio test detection in passive radar | |
CN116736342A (en) | GNSS spoofing detection method and system based on integrated CNN | |
Guven et al. | Classifying LPI radar waveforms with time-frequency transformations using multi-stage CNN system | |
CN108718223B (en) | Blind spectrum sensing method for non-cooperative signals | |
CN107271978B (en) | Target detection method based on Rao detection under multiple heterogeneous satellites | |
Oliinyk et al. | Center weighted median filter application to time delay estimation in non-Gaussian noise environment | |
CN114740501A (en) | Satellite navigation deception signal detection method based on frequency domain parameters | |
CN114089307A (en) | Radar detection and classification method and system under target and interference conditions | |
Stringer et al. | Application of Generative Machine Learning for Adaptive Detection with Limited Sample Support | |
Zhuo et al. | A subbandadaptive filter for DTV based passive radar | |
JP5328210B2 (en) | Incident wave number estimation device and incident wave number estimation method | |
JP2009250833A (en) | Device and method for estimating incident wave number | |
CN113759316B (en) | FFT-based linear frequency modulation continuous wave radar pre-detection parameter estimation method | |
Zhang et al. | A Low Complexity Algorithm for Time-Frequency Joint Estimation of CAF Based on Numerical Fitting | |
GOVALKAR | SIAMESE NETWORK BASED ONE-SHOT LEARNING FOR PULSE AND SIGNAL IDENTIFICATION APPLICATIONS A Thesis By | |
He et al. | CFAR detection of weak target in clutter using chaos synchronization |
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 |