CN111624561A - Filtering method for discrete multi-target signal equalization - Google Patents
Filtering method for discrete multi-target signal equalization Download PDFInfo
- Publication number
- CN111624561A CN111624561A CN202010337759.4A CN202010337759A CN111624561A CN 111624561 A CN111624561 A CN 111624561A CN 202010337759 A CN202010337759 A CN 202010337759A CN 111624561 A CN111624561 A CN 111624561A
- Authority
- CN
- China
- Prior art keywords
- frequency
- target
- filter
- echo
- signal
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
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/40—Means for monitoring or calibrating
- G01S7/4052—Means for monitoring or calibrating by simulation of echoes
- G01S7/4056—Means for monitoring or calibrating by simulation of echoes specially adapted to FMCW
-
- 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/50—Systems of measurement based on relative movement of target
- G01S13/52—Discriminating between fixed and moving objects or between objects moving at different speeds
- G01S13/536—Discriminating between fixed and moving objects or between objects moving at different speeds using transmission of continuous unmodulated waves, amplitude-, frequency-, or phase-modulated waves
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
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
The embodiment of the invention provides a filtering method for discrete multi-target signal equalization, which is characterized in that through designing parameters of a filter, beat signal frequencies corresponding to all targets are in a transition section of the filter, and the passband frequency of the filter is selected to be close to the maximum beat signal frequency, so that the attenuation of a near target is large, the attenuation of a far target is reduced, and the effect that echo amplitudes of targets with different distances are nearly consistent is realized.
Description
Technical Field
The invention relates to the field of radar signal processing, in particular to a filtering method for discrete multi-target signal equalization.
Background
Linear Frequency Modulated Continuous Wave (LFMCW) systems have the characteristics of small volume, low cost, extremely high resolution, no distance blind zone and being particularly suitable for near field applications, and are receiving continuous attention.
Conventional methods can distinguish well for single objects or objects where speed distinction is significant. However, as the number of targets increases, the targets no longer appear as separate peaks on the spectrum but are instead numerous; due to the influences of signal attenuation, noise interference and the like, the difference of the target echo amplitudes at different distances is large, namely the target echo amplitude at a close distance is large, the target echo amplitude at a far distance is small, and how to avoid annihilation of a long-distance target by a short-distance target is a key problem in multi-target signal processing. However, the literature has been very little concerned with this problem, and it is therefore necessary to develop a method for equalizing discrete multi-target echo signals.
Disclosure of Invention
The invention aims to provide a filtering method for discrete multi-target signal equalization, which is characterized by comprising the following steps of:
step 1: according to the characteristics of the linear frequency modulation continuous wave radar, the target distance R can be determinediDetermining the frequency of a beat signal of a target echo
Wherein f isiIs the beat signal frequency, B is the frequency modulation bandwidth, T is the frequency modulation period, and c is the speed of light.
According to the maximum target distance RmaxAnd a minimum target distance RminCalculating maximum and minimum beat frequencies
Step 2: according to different target distances RiDetermining the amplitude A of the target echo signalr,i
Assuming radar scattering cross-sections σ of different targetsiEquality, determining that the echo signal amplitude is inversely proportional to the square of the frequency, i.e. A, according to the relation between the target distance and the echo signal amplitude and the beat signal frequencyr,i∝1/fi 2(ii) a Assuming N targets, let minimum distance (R)min) The echo amplitude of the target is 1, and the echo amplitudes of other targets are
Then, the echo amplitudes of different frequencies are expressed as
And step 3: determining the stop band frequency f of a high pass filterstopAnd pass band frequency fpass。
Designing the stop band frequency of the high-pass filter to be less than the minimum beat signal frequency, i.e. fstop<fmin(ii) a Designing the passband frequency f of a high pass filterpassApproaching the maximum beat frequency as much as possible, e.g. satisfying | fpass-fmax|/fmax≤;
And 4, step 4: determining that the attenuation characteristic of the filter satisfies
p is the attenuation coefficient of the direct current signal (f is 0), and p < 1. The echo signal amplitude equalization of a plurality of targets with dispersed distance can be realized by passing the target echo through the designed filter.
Preferably, in step 3, the setting of the passband frequency includes: pass band frequency greater than maximum beat signal frequency fmaxI.e. all target beat signal frequenciesAre all located in the transition section (f) of the filterstop<fmin<…<fmax<fpass)。
Preferably, in step 3, the setting of the passband frequency includes: to ensure that the attenuation of the echo signal of the farthest target in the filter is as small as possible, the passband frequency is set at the maximum beat signal frequency fmaxFrequency f of sub-maximum beat signalN-1In between, i.e. fstop<fmin<…<fN-1<fpass≤fmax。
Preferably, in order to ensure the best filter effect, a cascaded filter design is adopted in consideration of the fact that the frequency range of beat signals of all targets exceeds one order of magnitude.
Preferably, step 3 is followed by:
step 3 a: assuming that N filters are designed, the linear segments of the N filter characteristics are represented as:
after a plurality of filters are cascaded, the total filter characteristics are respectively
Substituting (22) into (23) yields the overall filter characteristic
Preferably, step 3a is followed by step 3 b:
and (3) cascading the filters determined in the step (3 a), and processing the signal after the target echo passes through the filters, so that the amplitude equalization of the echo signals of a plurality of targets with dispersed distances can be realized.
In the process of processing the echo signals of the triangular frequency-modulated continuous waves, the invention provides a filtering method for discrete multi-target signal equalization, which can better ensure the equalization of the echo signals of various discrete targets and is convenient for realizing the measurement of a plurality of target parameters.
Drawings
FIG. 1 is a schematic diagram of echo signals of targets with different distances according to the present invention.
Fig. 2 is a schematic diagram of echo signal spectra and filter attenuation characteristics corresponding to different range targets in the present invention.
Fig. 3 is a schematic diagram of the amplitude-frequency characteristics of N high-pass filters according to the present invention.
Fig. 4 is a schematic diagram of the amplitude-frequency characteristic of the cascaded high-pass filter of the present invention.
FIG. 5 is a spectrum diagram of a echo signal before filtering in the present invention.
Fig. 6 is an amplitude-frequency characteristic curve of the cascaded high-pass filter of the present invention.
FIG. 7 is a spectrum diagram of echo signals before and after filtering in the present invention.
FIG. 8 is a flow chart of a method of the present invention.
Detailed Description
In order to make the aforementioned objects, features and advantages of the present invention comprehensible, embodiments accompanied with figures are described in further detail below.
The final implementation of the scheme is shown in figure 4. According to the characteristics of the linear frequency modulation continuous wave radar, the frequency of beat signals of target echoes and the amplitude of echo signals corresponding to targets with different distances can be determined according to the target distance. And designing parameters of a high-pass filter to enable the beat signal frequencies corresponding to all targets to be in a transition section of the filter, so that the echo amplitudes of different distances to the targets are approximately consistent. If the frequency ranges of beat signals of all targets are wide, a cascade filter design is adopted to ensure the optimal filter effect. The target echo is processed by the cascade filter, and the echo signal amplitude equalization of a plurality of targets with dispersed distance can be realized.
With reference to fig. 1 to 8, the specific implementation of the present invention is as follows:
step 1: according to the characteristics of chirped continuous wave radar, the firstThe frequency of beat signal of i target echoes can be determined by the target distance RiDetermining
Wherein f isiIs the beat signal frequency, B is the frequency modulation bandwidth, T is the frequency modulation period, and c is the speed of light.
According to the maximum target distance RmaxAnd a minimum target distance RminCalculating maximum and minimum beat frequencies
Step 2: according to different target distances RiDetermining the amplitude A of the target echo signalr,i
Assuming that the radar scattering cross sections sigma of different targets are equal, according to the relation between the target distance and the amplitude of the echo signal and the frequency of the beat signal, the inverse proportion of the amplitude of the echo signal to the square of the frequency can be determined, namely Ar,i∝1/fi 2See figure 1. Assuming N targets, let minimum distance target RminThe echo amplitude is 1, the echo amplitudes of other targets are
Then, the echo amplitudes of different frequencies can be expressed as
Therefore, a high-pass filter can be designed, so that the beat signal frequencies corresponding to all targets are in the transition section of the filter, and the echo amplitudes of different distances to the targets are basically approximately consistent.
And 3, step 3: determining the stop band frequency f of a high pass filterstopAnd pass band frequency fpass。
Setting the stop band frequency of the high-pass filter to be less than the minimum beat signal frequency, i.e. fstop<fmin. There are 2 schemes for setting the passband frequency: 1) pass band frequency greater than maximum beat signal frequency fmaxI.e. all target beat signal frequencies are located in the transition section (f) of the filterstop<fmin<…<fmax<fpass) (ii) a 2) To ensure that the attenuation of the echo signal of the farthest target in the filter is as small as possible, the passband frequency is set at the maximum beat signal frequency fmaxFrequency f of sub-maximum beat signalN-1In between, i.e. fstop<fmin<…<fN-1<fpass≤fmaxSee fig. 2.
And 4, step 4: considering that the frequency range of beat signals of all targets may exceed an order of magnitude, in order to ensure the best effect of the filter, a cascade filter design is adopted. Assume that N filters are designed (fig. 3). Note that f in the figurekIs the pass band frequency of each filter, is not the beat frequency corresponding to target k, and fN=fpass。
The linear segments of the N filter characteristic curves can be expressed as:
after a plurality of filters are cascaded, the total filter characteristics are respectively
Substituting (23) into (24) yields the total filter characteristic (FIG. 4)
And 5, step 5: and (4) cascading the filters determined in the step (4), and processing the target echo after passing through the filters to realize the amplitude equalization of the echo signals of a plurality of targets with dispersed distances.
The technical solution of the present invention will be described below with specific examples.
A filtering method for discrete multi-target signal equalization realizes discrete distribution target echo signal equalization through filter design. The following takes a triangular sweep signal with a center frequency of 24.125GHz and a bandwidth of 250MHz as an example to further describe the technical scheme of the present invention in detail. However, it should be noted that the present invention is not limited to 24GHz signals, but is a general method for designing a discrete multi-target signal equalization filter. The simulation parameters are shown in table 1.
TABLE 1
The realization process is as follows:
step 1: according to linear frequency modulation continuous wave radar parameters and target distance RiCalculating the frequency f of the target echo beat signali. According to the parameters in Table 1, the frequency of the echo beat signal of each target is [ f [ ]1,f2,f3, f4]=[16.67kHz,33.33kHz,83.33kHz,166.67kHz]. Thus, the minimum and maximum beat signal frequencies fmin=16.67kHz,fmax=166.67kHz。
Step 2: according to the distance R of different targetsiDetermining the amplitude A of the echo signalr,i. Let the minimum distance target RminThe echo amplitude is 1, and the echo amplitudes of other 3 targets are Ar,2=0.25,Ar,3=0.04, Ar,40.01. The superimposed signal of the echoes of different targets is fourier transformed, the spectrum of which is shown in fig. 5. It can be seen that the corresponding spectral amplitude of the target decreases with increasing frequency, and the near (low frequency) target has a much larger amplitude than the far (high frequency) target (by a factor of 80), which results in the far target being annihilated in the noise, resulting in a false alarm.
And 3, step 3: here we set fstop10kHz, for convenience, scheme 1) is employed to set the passband frequency, f pass180 kHz. According to the calculation, the maximum target echo amplitude is different from the minimum value by 40dB, so other parameters of the filter are as follows: the stopband amplitude is-50 dB, and the passband amplitude is 0 dB.
And 4, step 4: because the frequency range of the beat signal exceeds one order of magnitude, a segmented filter design is adopted to ensure the optimal effect of the filter. Here, we design 5 filters with filter parameters: (1) p is a radical of1=0.1,f1=20kHz;(2)p2=0.1,f2=40kHz;(3)p2=0.12,f2=80kHz;(4)p2=0.15, f2=120kHz;(5)p2=0.15,f 2180 kHz. The amplitude-frequency characteristic of the high-pass filter after the cascade connection is shown in figure 6.
And 5, step 5: and the echo signals of the targets pass through the filter and then are processed, so that the amplitude equalization of the echo signals of a plurality of targets with dispersed distances can be realized. The filtered echo signals are fourier transformed, see fig. 7. Compared with the frequency spectrum before filtering, the frequency spectrum amplitudes of the 4 filtered targets are relatively close, and the method is favorable for distinguishing different targets and calculating target information.
The beneficial effects of the invention include:
(1) can filter out the low frequency signal who reveals
Aiming at the low-frequency triangular wave frequency component leaked in the echo signal, the high-pass filter can effectively filter the low-frequency signal, and is convenient for amplification and target observation of the subsequent useful signal.
(2) Facilitating simultaneous measurement of multiple range discrete targets
Through designing a high-pass filter, the beat signal frequency corresponding to all targets is in the transition section of the filter, namely, the attenuation of a near target is large, and the attenuation of a far target is small, so that the echo amplitudes of different distance targets are approximately consistent, and the simultaneous measurement of a plurality of distance discrete targets is facilitated.
(3) The hardware circuit can be repeatedly configured, and various target test conditions can be met.
It will be apparent to those skilled in the art that various changes and modifications may be made in the invention without departing from the spirit and scope of the invention. Thus, if such modifications and variations of the present invention fall within the scope of the claims of the present invention and their equivalents, the present invention is also intended to include such modifications and variations.
Claims (6)
1. A filtering method for discrete multi-target signal equalization is characterized by comprising the following steps:
step 1: according to the characteristics of the linear frequency modulation continuous wave radar, the target distance R can be determinediDetermining the frequency of a beat signal of a target echo
Wherein f isiIs the beat signal frequency, B is the frequency modulation bandwidth, T is the frequency modulation period, and c is the speed of light.
According to the maximum target distance RmaxAnd a minimum target distance RminCalculating maximum and minimum beat frequencies
Step 2: according to different target distances RiDetermining the amplitude A of the target echo signalr,i
Assuming radar scattering cross-sections σ of different targetsiEquality, determining the inverse ratio of the echo signal amplitude to the square of the frequency, i.e. A, according to the relation between the target distance and the echo signal amplitude and the beat signal frequencyr,i∝1/fi 2(ii) a Assuming N targets, let minimum distance (R)min) The echo amplitude of the target is 1, and the echo amplitudes of other targets are
The amplitude of the echo at different frequencies is expressed as
And step 3: determining the stop band frequency f of a high pass filterstopAnd pass band frequency fpass。
Designing the stop band frequency of the high-pass filter to be less than the minimum beat signal frequency, i.e. fstop<fmin(ii) a Designing the passband frequency f of a high pass filterpassApproaching maximum beat frequency f as much as possiblemaxIf satisfy | fpass-fmax|/fmax≤;
And 4, step 4: determining that the attenuation characteristic of the filter satisfies
p is the attenuation coefficient of the direct current signal (f is 0), and p < 1. The echo signal amplitude equalization of a plurality of targets with dispersed distance can be realized by the target echo through the designed filter.
2. The filtering method for discrete multi-target signal equalization as claimed in claim 1, wherein the step 3, the setting of the passband frequency comprises: high pass band frequencyAt the maximum beat signal frequency fmaxI.e. all target beat signal frequencies are located in the transition section (f) of the filterstop<fmin<…<fmax<fpass)。
3. The filtering method for discrete multi-target signal equalization as claimed in claim 1, wherein the step 3, the setting of the passband frequency comprises: to ensure that the attenuation of the echo signal of the farthest target in the filter is as small as possible, the passband frequency is set at the maximum beat signal frequency fmaxFrequency f of sub-maximum beat signalN-1In between, i.e. fstop<fmin<…<fN-1<fpass≤fmax。
4. A method of discrete multi-target signal equalization filtering as claimed in claim 1, characterized in that, in order to ensure optimum filter effect, a cascaded filter design is used, taking into account the fact that the range of the beat signal frequency of all targets exceeds an order of magnitude.
5. The filtering method for discrete multi-target signal equalization as claimed in claim 4, further comprising after said step 3:
step 3 a: assuming N filters, the linear segments of the N filter characteristics are represented as:
after a plurality of filters are cascaded, the total filter characteristics are respectively
Substituting (7) into (8) to obtain the total filter characteristic;
6. the filtering method for discrete multi-target signal equalization as claimed in claim 5, wherein said step 3a is followed by the step 3b of:
and (3) cascading the filters determined in the step (3 a), and processing the target echo after passing through the filter bank to realize the amplitude equalization of the echo signals of a plurality of targets with dispersed distances.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010337759.4A CN111624561B (en) | 2020-04-26 | 2020-04-26 | Filtering method for discrete multi-target signal equalization |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010337759.4A CN111624561B (en) | 2020-04-26 | 2020-04-26 | Filtering method for discrete multi-target signal equalization |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111624561A true CN111624561A (en) | 2020-09-04 |
CN111624561B CN111624561B (en) | 2023-08-11 |
Family
ID=72257924
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010337759.4A Active CN111624561B (en) | 2020-04-26 | 2020-04-26 | Filtering method for discrete multi-target signal equalization |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111624561B (en) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108693524A (en) * | 2017-12-06 | 2018-10-23 | 深圳大学 | Multiple mobile object matching process based on Modulation Continuous Wave Radar and its system |
CN109270511A (en) * | 2018-11-28 | 2019-01-25 | 深圳大学 | Radar Signal Processing System and method based on FPGA and DSP |
US20190113600A1 (en) * | 2017-10-13 | 2019-04-18 | Infineon Technologies Ag | Radar sensing with interference suppression |
CN110118957A (en) * | 2019-05-06 | 2019-08-13 | 上海航天电子有限公司 | A kind of target matching method of the upper and lower frequency sweep of triangle CW with frequency modulation |
-
2020
- 2020-04-26 CN CN202010337759.4A patent/CN111624561B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20190113600A1 (en) * | 2017-10-13 | 2019-04-18 | Infineon Technologies Ag | Radar sensing with interference suppression |
CN108693524A (en) * | 2017-12-06 | 2018-10-23 | 深圳大学 | Multiple mobile object matching process based on Modulation Continuous Wave Radar and its system |
CN109270511A (en) * | 2018-11-28 | 2019-01-25 | 深圳大学 | Radar Signal Processing System and method based on FPGA and DSP |
CN110118957A (en) * | 2019-05-06 | 2019-08-13 | 上海航天电子有限公司 | A kind of target matching method of the upper and lower frequency sweep of triangle CW with frequency modulation |
Non-Patent Citations (4)
Title |
---|
MINGFEI LIU ET AL.: ""Multi-target Detection of FMCW Radar Based on Width Filtering"", 《EIDWT 2017: ADVANCES IN INTERNTWORKING, DATA & WEB TECHNOLOGY》 * |
张亚楠等: "一种脉压体制雷达检测多目标技术的研究", 《火控雷达技术》 * |
陆威: ""毫米波调频连续波雷达机动目标检测与跟踪技术研究"", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
颜卫忠 等: ""多目标识别LFMCW雷达的算法及DSP设计"", 《2018年全国微波毫米波会议论文集(上册)》 * |
Also Published As
Publication number | Publication date |
---|---|
CN111624561B (en) | 2023-08-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108415010B (en) | Radar multi-target detection method based on trapezoidal LFMCW modulation | |
Schroeder | Frequency‐correlation functions of frequency responses in rooms | |
Varshney et al. | Sidelobe reduction for matched filter range processing | |
Duk et al. | Target detection in sea-clutter using stationary wavelet transforms | |
JP2005530164A (en) | Method for suppressing interference in an object detection system | |
CN112882006B (en) | Millimeter wave short-range target detection method based on composite frequency modulation continuous wave | |
CN103344944A (en) | Radar pulse compression filter optimization design method applied to random signal waveforms | |
CN112859014A (en) | Radar interference suppression method, device and medium based on radar signal sorting | |
CN111896926A (en) | Low-altitude target detection method and system based on strong clutter suppression | |
Vizitiu | Some aspects of sidelobe reduction in pulse compression radars using NLFM signal processing | |
CN111812640A (en) | Through-wall imaging method for inhibiting clutter interference of wall body | |
CN111624561A (en) | Filtering method for discrete multi-target signal equalization | |
Galushko | On application of taper windows for sidelobe suppression in LFM pulse compression | |
CN109061626B (en) | Method for detecting low signal-to-noise ratio moving target by step frequency coherent processing | |
CN110632573A (en) | Airborne broadband radar space-time two-dimensional keystone transformation method | |
CN108680782B (en) | Voltage flicker parameter detection method based on extreme point symmetric mode decomposition | |
CN107769815A (en) | Linear frequency modulation short-range detecting system noise AM interference suppressing method | |
CN107689933B (en) | MFSK signal modulation identification method under Alpha stable distributed noise | |
CN112462343A (en) | Method for extracting radar signal intra-pulse characteristic parameters through frequency domain transformation processing | |
CN109085568B (en) | Frequency modulation continuous wave multi-target detection method based on secondary frequency mixing | |
Limin et al. | Low probability of intercept radar signal recognition based on the improved AlexNet model | |
CN109541556B (en) | Method for identifying frequency shift interference of linear frequency modulation signal | |
CN109460614B (en) | Signal time-frequency decomposition method based on instantaneous bandwidth | |
CN113176541B (en) | Self-adaptive anti-spectrum dispersion interference method and system | |
CN113820703B (en) | Unmanned aerial vehicle target rotor wing parameter estimation method based on scattering transformation |
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 |