CN113050043A - Ground penetrating radar ultra wide band Gaussian pulse FRI sampling method based on non-ideal LPF - Google Patents
Ground penetrating radar ultra wide band Gaussian pulse FRI sampling method based on non-ideal LPF Download PDFInfo
- Publication number
- CN113050043A CN113050043A CN202110323328.7A CN202110323328A CN113050043A CN 113050043 A CN113050043 A CN 113050043A CN 202110323328 A CN202110323328 A CN 202110323328A CN 113050043 A CN113050043 A CN 113050043A
- Authority
- CN
- China
- Prior art keywords
- sampling
- signal
- ideal
- expressed
- sample
- 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.)
- Pending
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/28—Details of pulse systems
-
- 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
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/02—Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
- G01S13/0209—Systems with very large relative bandwidth, i.e. larger than 10 %, e.g. baseband, pulse, carrier-free, ultrawideband
-
- 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
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/885—Radar or analogous systems specially adapted for specific applications for ground probing
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
A ground penetrating radar ultra wide band Gaussian pulse FRI sampling method based on a non-ideal LPF comprises the following steps: generating an original signal; step two, constructing a sampling kernel function; step three, low-speed sampling modeling; step four, sampling the sample y [ n ]]And obtaining and processing; step five, sampling a sample h [ n ]]And obtaining and processing; step six, representing the signal relation; step seven, in the data processing stage, the formula (9) only contains unknown parametersIs a typical parameter estimation problem, solving equation (9) yields the result. The invention provides a ground penetrating radar ultra wide band Gaussian pulse FRI sampling method based on a non-ideal LPF, which reconstructs the processed sampling informationThe non-ideal effects brought by the filter can be eliminated, thereby improving the reconstruction accuracy of the system.
Description
Technical Field
The invention relates to the technical field of signal processing, in particular to a ground penetrating radar ultra wide band Gaussian pulse FRI sampling method based on a non-ideal LPF.
Background
With the development of communication technology in recent years, the bandwidth of signals used in the field of radar is continuously increasing. According to the traditional Nyqiust sampling theorem, the sampling rate of a signal is required to be more than twice of the highest frequency of the signal, so that the original signal can be reconstructed without distortion. In order to effectively reduce the sampling frequency of the signal and effectively reconstruct the original signal, many experts and scholars have conducted a series of researches on the undersampling method. The conventional under-sampling methods include a Compressed Sensing (CS) method and a Finite Innovation Rate (FRI) method. FRI sampling theory was first proposed in 2002 by Vetterli et al. According to the theory, various parameters of the parameterized signal can be effectively reconstructed by acquiring partial frequency domain information of the parameterized signal and processing the frequency domain information by using a specific algorithm, so that the signal is reconstructed. The classic FRI sampling structure is shown in fig. 1.
In fig. 1, the most important part is the sampling kernel g (t), which is the unit impulse response of the filter in practical application. After the signal is processed by the sampling kernel function, partial frequency domain information of the signal can be obtained, and then the sampling sample y [ n ] can be obtained by low-speed sampling, wherein the sampling sample comprises partial frequency domain information of the original signal x (t). The frequency domain information may then be processed according to existing nulling filter methods or subspace estimation methods, which may complete the reconstruction of the original signal. The existing common sampling kernel functions include a sinc sampling kernel, an SOS sampling kernel, a B spline function and the like.
Taking a basic sinc sampling kernel as an example, modeling is carried out on the process of carrying out FRI sampling reconstruction on the ultra-wideband Gaussian pulse of the ground penetrating radar. Using a ground penetrating radar ultra wide band gaussian pulse as an original signal x (t), in a radar system, a complete echo signal can be represented by superposition of received echoes, and the radar echo signal is represented without considering noise and other interference as follows:
where h (t) is a known Gaussian pulse shape,respectively corresponding to the amplitude parameter and the time delay parameter of the detected pulse, wherein T is the pulse repetition interval of the radar, and after the original signal is processed by a sampling kernel function, namely y (T) is expressed as the convolution process of the original signal and the sampling kernel function:
the sampled values y [ n ] can be obtained through low-speed sampling, and y [ n ] is expressed as:
wherein T isSFor the sampling period of the signal, it can be seen from the formula that the sampling samples contain partial fourier coefficients of the original signal, and then a typical parameter estimation problem can be obtained by simplification, and by solving the problem, the amplitude and delay parameters of the signal can be obtained, i.e. the signal is successfully reconstructed.
A key component in FRI sampling systems is the selection of the sampling kernel, i.e., the selection of the filter. However, in the hardware implementation process, because the filter cannot achieve the effect of simulation, a non-ideal effect is brought, and reconstruction accuracy is influenced. How to eliminate the non-ideal effect of the filter is an important problem influencing the reconstruction precision of the ground penetrating radar ultra wide band Gaussian pulse FRI sampling system.
Disclosure of Invention
In order to overcome the defects of the prior art and solve the problem of non-ideal effect of a filter in a ground penetrating radar ultra wide band Gaussian pulse FRI sampling system, the invention provides a non-ideal LPF-based ground penetrating radar ultra wide band Gaussian pulse FRI sampling method, an original signal can be expressed as a combination of base signals with different time delays and amplitudes, the original signal and the base signals are processed by the method of the invention, so that two groups of sampling information containing the non-ideal effect of the filter can be obtained, the sampling information without the non-ideal effect of the filter can be obtained after processing, the non-ideal effect brought by the filter can be eliminated by reconstructing the processed sampling information, and the reconstruction precision of the system is improved.
The technical scheme adopted by the invention for solving the technical problems is as follows:
a ground penetrating radar ultra wide band Gaussian pulse FRI sampling method based on a non-ideal LPF comprises the following steps:
step one, generating an original signal: parameter setting of an original signal, the original signal being represented asWherein T ∈ [0, T) is the observation time, for unknown delay and amplitude parameters, L is the number of signals x (t) including base signals h (t), which are expressed as:
step two, constructing a sampling kernel function: the basic sinc kernel is taken as a sampling kernel function, and non-ideal effects existing in the hardware implementation process are considered, so that the sampling kernel function is expressed asIts continuous time Fourier transform can be expressed as
Step three, low-speed sampling modeling: using a delay interval of t-nTsThe continuous impulse string as the sampling function needs to satisfy the requirement that the low-speed sampling frequency domain needs more than twice the cut-off frequency of the filter, namely 1/TsF is more than 2f, and f is the cut-off frequency of the non-ideal filter;
step four, sampling the sample y [ n ]]And acquiring and processing: the original signal x (T) is sampled at a low speed after passing through a non-ideal filter, the sampling rate being 1/Ts> 2f, obtaining sample y [ n ]]Then, fourier transform is performed to obtain:
from equation (2)For a signal extended by one period, one period is taken, that is, the signal has complete frequency domain information, and if n is 0, the signal is expressed as:
step five, sampling a sample h [ n ]]And acquiring and processing: the base signal h (T) is passed through a non-ideal filter and then sampled at a low rate of 1/Ts> 2f, obtaining a sample h [ n ]]Then, fourier transform is performed to obtain:
in the same step four, let n be 0, formula (4) is expressed as:
step six, representing the signal relation: between the original signal x (t) and the base signal h (t)In at restThe relationship, in the frequency domain, is expressed as:
Step seven, a data processing stage: and D, processing the frequency domain information obtained in the fourth step and the fifth step to obtain:
after simplification, it is written as:
the formula (9) only contains unknown parametersIs a typical parameter estimation problem, solving equation (9) yields the result.
Further, in the sixth step, equation (9) is solved by using a nulling filter method or a subspace estimation method.
The method of the invention is different from the traditional FRI undersampling system, and considers and solves the influence of non-ideal effect on reconstruction precision. The method of the invention processes the base signal h (t) as the original signal, and the original signal and the base signal both contain non-ideal effects. Then the non-ideal effect is eliminated through processing, and finally the signal is reconstructed by using the frequency domain information with the non-ideal effect eliminated. The method eliminates the non-ideal effect caused by sampling kernel function, and can improve the reconstruction precision to a great extent.
The invention has the following beneficial effects: the reconstruction precision of the system is improved, and the anti-noise capability is strong.
Drawings
Fig. 1 is a diagram of a classic FRI sampling architecture.
Fig. 2 is a block diagram of the FRI sampling system structure of the method of the present invention.
Fig. 3 is a graph showing the reconstruction effect of experiment two, in which (a) shows the reconstruction effect of the FRI system and (b) shows the reconstruction effect of the system according to the method of the present invention.
Fig. 4 is a comparison graph of system noise resistance.
Detailed Description
The invention is further described below with reference to the accompanying drawings.
Referring to fig. 2 to 4, a structural block diagram of a ground penetrating radar ultra wide band gaussian pulse FRI sampling method based on a non-ideal LPF is shown in fig. 2, and the method comprises the following steps:
step one, generating an original signal: parameter settings of the original signal, which can be represented asWherein T ∈ [0, T) is the observation time, for unknown delay and amplitude parameters, L is the number of signals x (t) including base signals h (t), which are expressed as:
step two, constructing a sampling kernel function: by the basic sincThe kernel is a sampling kernel function, and takes into account non-ideal effects existing during hardware implementation, so the sampling kernel function can be expressed asIts continuous time Fourier transform is expressed as
Step three, low-speed sampling modeling: using a delay interval of t-nTsAs a sampling function, so that the requirement of a low-speed sampling frequency domain for more than twice the cut-off frequency of the filter, i.e. 1/T, is satisfiedsF is more than 2f, and f is the cut-off frequency of the non-ideal filter;
step four, sampling the sample y [ n ]]And acquiring and processing: the original signal x (T) is sampled at a low speed after passing through a non-ideal filter, the sampling rate being 1/Ts> 2f, obtaining sample y [ n ]]Then, fourier transform is performed to obtain:
the formula (1.14) showsFor a signal extended by one period, one period is taken, that is, the signal has complete frequency domain information, and if n is 0, the signal is expressed as:
step five, sampling a sample h [ n ]]And acquiring and processing: the base signal h (T) is passed through a non-ideal filter and then sampled at a low rate of 1/Ts> 2f, obtaining a sample h [ n ]]Then, fourier transform is performed to obtain:
in the same step four, let n be 0, and formula (4) is expressed as:
step six, representing the signal relation: between the original signal x (t) and the base signal h (t)The relationship, in the frequency domain, is expressed as:
Step seven, a data processing stage: and D, processing the frequency domain information obtained in the fourth step and the fifth step to obtain:
after simplification, it is written as:
the formula (8) contains only unknown parametersIs a typical parameter estimation problem, using a nulling filter approach or subspace estimationAnd (6) counting and solving.
In order to verify the reconstruction effect and the anti-noise performance of the method, MATLAB software is used for simulating the method. Use signalAs the original signal, the amplitude parameter is set to al=[0.9,0.6,0.8,0.5]With time delay parameter set to tl=[0.2,0.4,0.6,0.8]Using the formulaGenerating a base signal, t0For the initial delay of the base signal, set to 0.5, σ is set to 0.001.
Experiment one: the use of an idealized sampling kernel, i.e. the impulse effect of the filter, is represented as a typical cut-off characteristic. By using the above signals to perform experiments on the reconstruction effect of the FRI sampling system and the method of the present invention, it can be observed that the reconstructed amplitude parameters and time delay parameters are shown in table 1.
TABLE 1
It can be seen from table 1 that both the FRI sampling system and the method of the present invention can reconstruct the original signal with high accuracy when using an idealized filter.
And secondly, performing simulation experiments on the FRI sampling system and the method by using the non-ideal sampling kernel function. The non-ideal filter used in this experiment was a Chebyshev type I filter, and the results of the experiment are shown in fig. 3(a) and 3 (b).
It can be seen from fig. 3(a) that the reconstruction effect is adversely affected when the non-ideal sampling kernel is used, and the original signal can be reconstructed with high precision by using the method of the present invention.
And thirdly, performing an anti-noise experiment on the FRI sampling system and the method. Gaussian white noise was added to both systems, gradually increasing from 0dB to 100dB in steps of 5 dB. For the convenience of comparison, two different parameter reconstruction methods are introduced here to reconstruct the signal, including a nulling filter method and a subspace estimation method. Reconstruction accuracy the reconstruction delay accuracy is measured using Normalized Mean-Square Error (NMSE), which is expressed as follows:
the experimental result is shown in fig. 4, and it can be seen from the graph that the anti-noise performance of the method of the present invention is significantly better than that of the FRI sampling system.
The embodiments described in this specification are merely illustrative of implementations of the inventive concepts, which are intended for purposes of illustration only. The scope of the present invention should not be construed as being limited to the particular forms set forth in the examples, but rather as being defined by the claims and the equivalents thereof which can occur to those skilled in the art upon consideration of the present inventive concept.
Claims (2)
1. A ground penetrating radar ultra wide band Gaussian pulse (FRI) sampling method based on a non-ideal LPF is characterized by comprising the following steps:
step one, generating an original signal: parameter setting of an original signal, the original signal being represented asWherein T ∈ [0, T) is the observation time,for unknown delay and amplitude parameters, L is the number of signals x (t) including base signals h (t), which are expressed as:
step two, constructing a sampling kernel function: taking a basic sinc kernel as a sampling kernel function,and takes into account non-ideal effects that exist during hardware implementation, the sampling kernel is expressed asIts continuous time Fourier transform can be expressed as
Step three, low-speed sampling modeling: using a delay interval of t-nTsThe continuous impulse string as the sampling function needs to satisfy the requirement that the low-speed sampling frequency domain needs more than twice the cut-off frequency of the filter, namely 1/TsF is more than 2f, and f is the cut-off frequency of the non-ideal filter;
step four, sampling the sample y [ n ]]And acquiring and processing: the original signal x (T) is sampled at a low speed after passing through a non-ideal filter, the sampling rate being 1/Ts> 2f, obtaining sample y [ n ]]Then, fourier transform is performed to obtain:
from equation (2)For a signal extended by one period, one period is taken, that is, the signal has complete frequency domain information, and if n is 0, the signal is expressed as:
step five, sampling a sample h [ n ]]And acquiring and processing: the base signal h (T) is passed through a non-ideal filter and then sampled at a low rate of 1/Ts> 2f, obtaining a sample h [ n ]]Then, fourier transform is performed to obtain:
in the same step four, let n be 0, formula (4) is expressed as:
step six, representing the signal relation: between the original signal x (t) and the base signal h (t)The relationship, in the frequency domain, is expressed as:
Step seven, a data processing stage: and D, processing the frequency domain information obtained in the fourth step and the fifth step to obtain:
after simplification, it is written as:
2. The FRI sampling method for ultra-wideband Gaussian pulses of ground penetrating radar based on non-ideal LPF as claimed in claim 1, wherein in the sixth step, equation (9) is solved by using the nulling filter method or subspace estimation method.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110323328.7A CN113050043A (en) | 2021-03-26 | 2021-03-26 | Ground penetrating radar ultra wide band Gaussian pulse FRI sampling method based on non-ideal LPF |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110323328.7A CN113050043A (en) | 2021-03-26 | 2021-03-26 | Ground penetrating radar ultra wide band Gaussian pulse FRI sampling method based on non-ideal LPF |
Publications (1)
Publication Number | Publication Date |
---|---|
CN113050043A true CN113050043A (en) | 2021-06-29 |
Family
ID=76515297
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110323328.7A Pending CN113050043A (en) | 2021-03-26 | 2021-03-26 | Ground penetrating radar ultra wide band Gaussian pulse FRI sampling method based on non-ideal LPF |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113050043A (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117572416A (en) * | 2024-01-15 | 2024-02-20 | 深圳安德空间技术有限公司 | Centimeter-level ultra-wideband ground penetrating radar system based on time expansion framework |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20070143078A1 (en) * | 2001-03-26 | 2007-06-21 | Martin Vetterli | Sampling method, reconstruction method, and device for sampling and/or reconstructing signals |
US20100259428A1 (en) * | 2007-10-30 | 2010-10-14 | Technion City | Method and apparatus for reconstructing digitized distorted signals |
CN102143108A (en) * | 2011-03-17 | 2011-08-03 | 电子科技大学 | Improved self-adaption predistortion technology |
US20130245471A1 (en) * | 2012-03-16 | 2013-09-19 | Qualcomm Incorporated | System and method for analysis and reconstruction of variable pulse-width signals having low sampling rates |
CN105814603A (en) * | 2013-05-30 | 2016-07-27 | 帝国创新有限公司 | Method and apparatus |
CN106788334A (en) * | 2016-12-07 | 2017-05-31 | 哈尔滨工业大学 | A kind of compensation method for modulating the non-ideal filter of wide-band transducer |
CN108984474A (en) * | 2018-06-04 | 2018-12-11 | 哈尔滨工业大学 | A kind of lack sampling method of non-ideal piecewise polynomial signal |
CN109889231A (en) * | 2019-02-01 | 2019-06-14 | 哈尔滨工业大学 | Burst signal lack sampling method based on random demodulation and the limited new fixed rate of interest |
CN112468114A (en) * | 2020-10-14 | 2021-03-09 | 浙江工业大学 | FRI sampling system and method based on non-ideal sinc kernel |
-
2021
- 2021-03-26 CN CN202110323328.7A patent/CN113050043A/en active Pending
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20070143078A1 (en) * | 2001-03-26 | 2007-06-21 | Martin Vetterli | Sampling method, reconstruction method, and device for sampling and/or reconstructing signals |
US20100259428A1 (en) * | 2007-10-30 | 2010-10-14 | Technion City | Method and apparatus for reconstructing digitized distorted signals |
CN102143108A (en) * | 2011-03-17 | 2011-08-03 | 电子科技大学 | Improved self-adaption predistortion technology |
US20130245471A1 (en) * | 2012-03-16 | 2013-09-19 | Qualcomm Incorporated | System and method for analysis and reconstruction of variable pulse-width signals having low sampling rates |
CN105814603A (en) * | 2013-05-30 | 2016-07-27 | 帝国创新有限公司 | Method and apparatus |
CN106788334A (en) * | 2016-12-07 | 2017-05-31 | 哈尔滨工业大学 | A kind of compensation method for modulating the non-ideal filter of wide-band transducer |
CN108984474A (en) * | 2018-06-04 | 2018-12-11 | 哈尔滨工业大学 | A kind of lack sampling method of non-ideal piecewise polynomial signal |
CN109889231A (en) * | 2019-02-01 | 2019-06-14 | 哈尔滨工业大学 | Burst signal lack sampling method based on random demodulation and the limited new fixed rate of interest |
CN112468114A (en) * | 2020-10-14 | 2021-03-09 | 浙江工业大学 | FRI sampling system and method based on non-ideal sinc kernel |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117572416A (en) * | 2024-01-15 | 2024-02-20 | 深圳安德空间技术有限公司 | Centimeter-level ultra-wideband ground penetrating radar system based on time expansion framework |
CN117572416B (en) * | 2024-01-15 | 2024-03-29 | 深圳安德空间技术有限公司 | Centimeter-level ultra-wideband ground penetrating radar system based on time expansion framework |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Tur et al. | Innovation rate sampling of pulse streams with application to ultrasound imaging | |
CN204495996U (en) | broadband sonar receiver | |
CN102508206B (en) | Linear frequency modulation (FM) signal parameter estimation method based on small-wave-packet denoising and power spectral entropy | |
CN110456351B (en) | Maneuvering target ISAR imaging method based on time-varying amplitude LFM signal parameter estimation | |
EP2367293A1 (en) | Low-rate sampling of pulse streams | |
CN108984474B (en) | Under-sampling method of non-ideal piecewise polynomial signal | |
CN106199532B (en) | Based on mixing Fourier-wavelet analysis Gpr Signal noise-reduction method | |
CN108844617B (en) | The adaptive sparse tree construction noise-reduction method of main reducing gear very noisy vibration signal | |
CN109738878A (en) | Radar range profile's recognition methods based on compressed sensing and frequency stepping waveform | |
CN111562597A (en) | Beidou satellite navigation interference source identification method based on BP neural network | |
CN109361376A (en) | A kind of high-accuracy time delay estimation method of Higher Order Cumulants | |
CN111224672A (en) | Multi-harmonic signal undersampling method based on multi-channel time delay | |
CN111159888A (en) | Covariance matrix sparse iteration time delay estimation method based on cross-correlation function | |
CN113608259A (en) | Seismic thin layer detection method based on ICEEMDAN constraint generalized S transformation | |
CN113050043A (en) | Ground penetrating radar ultra wide band Gaussian pulse FRI sampling method based on non-ideal LPF | |
CN110032968B (en) | Denoising method based on dual-tree complex wavelet and self-adaptive semi-soft threshold method | |
CN114690003A (en) | EEMD-based partial discharge signal noise reduction method | |
CN112652290B (en) | Method for generating reverberation audio signal and training method of audio processing model | |
CN112468114B (en) | FRI sampling system and method based on non-ideal sinc core | |
CN113589253A (en) | Method for detecting weak echo signal based on wavelet transform algorithm of pseudo time domain | |
Zhu et al. | Research on NMR noise reduction method based on improved CEEMD | |
CN110808929A (en) | Real-complex conversion type signal-to-noise ratio estimation algorithm of subtraction strategy | |
CN114545353B (en) | Pulse Doppler signal undersampling and parameter estimation method based on FRI sampling and Chinese remainder theorem | |
CN103971012B (en) | Lamb wave signal denoising method based on John Saris model and fractional differentiation | |
CN116184333A (en) | Linear frequency modulation signal parameter estimation method based on local iterative filtering |
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 |