CN109271855A - A kind of feature extracting method of industry control signal - Google Patents
A kind of feature extracting method of industry control signal Download PDFInfo
- Publication number
- CN109271855A CN109271855A CN201810902497.4A CN201810902497A CN109271855A CN 109271855 A CN109271855 A CN 109271855A CN 201810902497 A CN201810902497 A CN 201810902497A CN 109271855 A CN109271855 A CN 109271855A
- Authority
- CN
- China
- Prior art keywords
- signal
- sample
- physical signal
- industry control
- control 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
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/08—Feature extraction
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/12—Classification; Matching
Abstract
A kind of feature extracting method of industry control signal, comprising: S1, acquisition simultaneously pre-process industry control signal to obtain physical signal sample；S2, at least two groups physical signal sample is chosen from the physical signal sample to do the variogram of the attribute byte of different classes of physical signal sample；S3, related coefficient curve is obtained based on the variogram；S4, the characteristic threshold value that different classes of physical signal sample is obtained based on the related coefficient curve.Implement the feature extracting method of industry control signal of the invention, the feature of industry control signal can be extracted by succinct operation, and high for the feature recognition efficiency of industry control signal, had a wide range of application.
Description
Technical field
The present invention relates to signal detection fields, more specifically to a kind of feature extracting method of industry control signal.
Background technique
As the requirement to industry control safety is higher and higher, existing industry control safety detection is all realized in network protocol detection layers
No matter face in network protocol level constructs the protection of which kind of degree, all not can avoid the various network attack means to emerge one after another
Intrusion.Industrial control equipment physical signal layer is undoubtedly the most fundamental detection basis, and foundation is more in the detection means of signals layer
A kind of reliable detection method.However the different classes of corresponding function of industry control signal is also different, their feature also has certainly
It is distinguished.To rapidly identify the type of the industry control signal received, judge that it is to belong to normal signal, fault-signal also
It is other intrusive viruses signals, it is necessary to which enough signal characteristic knowledge is arrived in study, it is desirable that acquires enough samples to mention
Take the feature of inhomogeneity signal.Therefore, it is necessary to a kind of operations succinctly, and recognition efficiency is high, the spy for the industry control signal having a wide range of application
Levy extracting method.
Summary of the invention
The technical problem to be solved in the present invention is that in view of the above drawbacks of the prior art, a kind of industry control signal is provided
Feature extracting method.
The technical solution adopted by the present invention to solve the technical problems is: constructing a kind of feature extraction side of industry control signal
Method, comprising:
S1, acquisition simultaneously pre-process industry control signal to obtain physical signal sample；
S2, at least two groups physical signal sample is chosen from the physical signal sample to do different classes of physical signal
The variogram of the attribute byte of sample；
S3, related coefficient curve is obtained based on the variogram；
S4, the characteristic threshold value that different classes of physical signal sample is obtained based on the related coefficient curve.
In the feature extracting method of industry control signal of the present invention, the step S1 further comprises:
S11, the sort instructions signal for obtaining industrial control equipment；
The physics letter of the corresponding mark label of sort instructions signal described in S12, the excitation sort instructions signal and recording
Number；
S13, the physical signal is normalized to obtain the physical signal sample.
In the feature extracting method of industry control signal of the present invention, the step S11 further comprises:
The instruction classification of S111, the selected industrial control equipment that will be detected and the selected industrial control equipment that will be detected
Scheme；
S112, based on described instruction classification schemes by way of labelling to the industry control signal of the industrial control equipment
Classify to obtain sort instructions signal.
In the feature extracting method of industry control signal of the present invention, the step S2 further comprises:
S21, selection at least two groups quantity is identical from the physical signal sample and covers the physics of codomain variation range
Sample of signal；
S22, variance point is obtained according to following formula based on the physical signal sample, and variance is done based on the variance point
Figure:
Wherein σ2Indicate the variance of the corresponding attribute byte of different types of physical signal sample, N indicates the physics of acquisition
The sample size of sample of signal, xiIndicate the amplitude size of physics sample of signal, μ indicates the equal of the amplitude of physics sample of signal
Value.
In the feature extracting method of industry control signal of the present invention, the step S3 further comprises:
S31, the waveform that the physical signal sample where the attribute byte is obtained based on the variogram；
S32, the waveform mean value of each classification is calculated using the reference signal of the attribute byte as the category；
S33, related coefficient curve is obtained based on the reference signal and the range of waveforms of the physical signal sample.
In the feature extracting method of industry control signal of the present invention, the step S33 further comprises, based on as follows
Formula calculates the correlation coefficient value of the related coefficient curve:
Wherein ρXYIndicate correlation coefficient value, Cov (X, Y) indicates the covariance of stochastic variable X and Y, and D (X) and D (Y) are respectively
Indicate the variance of stochastic variable X and Y, μxIndicate the mean value of stochastic variable X, μyThe mean value of stochastic variable Y is indicated, wherein random become
Measuring X indicates the reference signal, and stochastic variable Y indicates the physical signal sample in the range of waveforms of the physical signal sample
Point.
Another technical solution that the present invention solves the use of its technical problem is to construct a kind of computer readable storage medium,
It is stored thereon with computer program, the feature extracting method of the industry control signal is realized when described program is executed by processor.
Implement the feature extracting method of industry control signal of the invention, the spy of industry control signal can be extracted by succinct operation
Sign, and it is high for the feature recognition efficiency of industry control signal, have a wide range of application.
Detailed description of the invention
Present invention will be further explained below with reference to the attached drawings and examples, in attached drawing:
Fig. 1 is the flow chart of the feature extracting method of the industry control signal of first embodiment according to the present invention；
Fig. 2 is the flow chart of the feature extracting method of the industry control signal of first embodiment according to the present invention；
Fig. 3 is the schematic diagram of sort instructions signal of the invention；
Fig. 4 is the schematic diagram that attribute byte of the invention extracts.
Specific embodiment
In order to make the objectives, technical solutions, and advantages of the present invention clearer, with reference to the accompanying drawings and embodiments, right
The present invention is further elaborated.It should be appreciated that the specific embodiments described herein are merely illustrative of the present invention, and
It is not used in the restriction present invention.
The present invention relates to a kind of feature extracting methods of industry control signal, comprising: S1, acquisition simultaneously pre-process industry control signal to obtain
Obtain physical signal sample；S2, selection at least two groups physical signal sample is different classes of to do from the physical signal sample
The variogram of the attribute byte of physical signal sample；S3, related coefficient curve is obtained based on the variogram；S4, based on described
Related coefficient curve obtains the characteristic threshold value of different classes of physical signal sample.The feature for implementing industry control signal of the invention mentions
Method is taken, the feature of industry control signal can be extracted by succinct operation, and high for the feature recognition efficiency of industry control signal,
Have a wide range of application.
Fig. 1 is the flow chart of the feature extracting method of the industry control signal of first embodiment according to the present invention.In step S1
In, it obtains and pre-processes industry control signal to obtain physical signal sample.In a preferred embodiment of the invention, it can first obtain
The command signal for taking the industry control signal of the industrial control equipment, then divides command signal by way of labelling
Class, to obtain sort instructions signal.Then can be by executing sort instructions signal, then recording obtains the physical signal
Sample.In further preferred embodiment of the invention, physical signal sample can be normalized.
In step s 2, it is different classes of to do that at least two groups physical signal sample is chosen from the physical signal sample
The variogram of the attribute byte of physical signal sample.In a preferred embodiment of the invention, it is selected from the physical signal sample
It takes at least two groups quantity identical and covers the physical signal sample of codomain variation range, be then based on the physical signal sample and obtain
Take variance point.Certainly, in other preferred embodiments of the invention, those skilled in the art can carry out respectively according to actual needs
Kind signal behavior, and variogram can be made using any method as known in the art.
In step s3, related coefficient curve is obtained based on the variogram.In a preferred embodiment of the invention,
The waveform of the physical signal sample where the attribute byte is obtained based on the variogram.Further, it calculates each
The waveform mean value of classification is using the reference signal of the attribute byte as the category.Subsequently, based on the reference signal and the object
The range of waveforms for managing sample of signal obtains related coefficient curve.In the present invention, curved section where attribute byte known to analysis
Variance yields is larger, and size is 0.5 or so, by the range of waveforms of the classification data of the available this feature of this characteristic.It sees
Examine this feature byte includes how many a classifications, then calculates each classification and correspond to mean value in dimension, the reference as the category
Signal.After obtaining reference signal, the point of itself and normalized signal character pair range is sought into related coefficient, finally obtains the ginseng
Examine the related coefficient curve under signal.
In step s 4, the characteristic threshold value of different classes of physical signal sample is obtained based on the related coefficient curve.
In a preferred embodiment of the invention it is found that it is bent that a plurality of related coefficient will be obtained if this feature byte includes multiple classifications
Line, the distribution tendency by observing curve, which can determine, identifies the threshold value for belonging to different types of classification data.
Implement the feature extracting method of industry control signal of the invention, the spy of industry control signal can be extracted by succinct operation
Sign, and it is high for the feature recognition efficiency of industry control signal, have a wide range of application.
Fig. 2 is the flow chart of the feature extracting method of the industry control signal of first embodiment according to the present invention.In step S1
In, obtain the sort instructions signal of industrial control equipment.It is selected first to be detected in of the invention one preferred city embodiment
Industrial control equipment and the industrial control equipment that will be detected instruction classification scheme, be then based on described instruction classification side
Case classifies to obtain sort instructions signal to the industry control signal of the industrial control equipment by way of labelling.Fig. 3
Show sort instructions signal of the invention.
In step s 2, the corresponding mark label of sort instructions signal described in the sort instructions signal and recording is motivated
Physical signal.In a preferred embodiment of the invention, the sort instructions signal can be motivated and based on triggering-suspension record
The physical signal of the corresponding mark label of sort instructions signal described in wave.In further preferred embodiment of the invention
In, the command function of sort instructions signal described in the playback sort instructions signal authentication can be passed through.By excitation, institute is active
Its available described physical signal for having tag along sort of sort instructions signal that can classify.Those skilled in the art know,
Known any triggering-suspension condition setting may be used to the present invention in signal transmission field.Those skilled in the art can be with
Correlation selection is carried out according to actual needs.The present invention is not limited by specific triggering-suspension condition.
In step s3, the physical signal is normalized to obtain the physical signal sample.Usually may be used
To be handled using the positive and negative upper and lower limits normalization of voltage value the waveform of the physical signal.
In step s 4, selection at least two groups quantity is identical from the physical signal sample and covers codomain variation range
Physical signal sample.
In step s 5, variance point is obtained according to setting formula based on the physical signal sample, and is based on the variance
Point does variogram.In other preferred embodiments of the invention, those skilled in the art can use ability according to actual needs
Known any method makes variogram in domain.
In a preferred embodiment of the invention, each variance point σ on variogram can be calculated by following formula2
Wherein σ2Indicate the variance of different types of sort instructions signal character pair byte, N indicates that the classification of acquisition refers to
Enable the sample size of signal, xiIndicate the amplitude size of sample point, μ indicates the mean value of sample point amplitude.
In step s 6, the wave of the physical signal sample where the attribute byte is obtained based on the variogram
Shape.Those skilled in the art know that the waveform of the physical signal sample can also be using described where the attribute byte
The serial number of physical signal sample indicates.In a preferred embodiment of the invention, curved section where attribute byte known to analysis
Variance yields is larger, and size is 0.5 or so, by the range of waveforms of the classification data of the available this feature of this characteristic.It sees
Examine this feature byte includes how many a classifications.
In the step s 7, the waveform mean value of each classification is calculated using the reference signal of the attribute byte as the category.Fig. 4
Show attribute byte of the invention.In the present invention, can observe this feature byte includes how many a classifications, then calculate every
A classification corresponds to the mean value in dimension, the reference signal as the category.In other embodiments of the invention, it can also use
Other reference signals.
In step s 8, related coefficient song is obtained based on the reference signal and the range of waveforms of the physical signal sample
Line.
In further preferred embodiment of the invention, the visualization formula of related coefficient algorithm are as follows:
Wherein ρXYIndicate correlation coefficient value, Cov (X, Y) indicates the covariance of stochastic variable X and Y, and D (X) and D (Y) are respectively
Indicate the variance of stochastic variable X and Y, abscissa X is the line number of the characteristic point, and ordinate Y is the related coefficient.
In the present invention, after obtaining reference signal, the point of itself and normalized signal character pair range is asked into phase relation
Number, finally obtains the related coefficient curve under the reference signal.If this feature byte includes multiple classifications, will obtain a plurality of
Related coefficient curve, the distribution tendency by observing curve, which can determine, to be identified and belongs to different types of physical signal sample
Threshold value.
In further preferred embodiment of the invention, (it can be tested test sample is received using the threshold value
Industry control signal) it is detected.For example, can by calculate the correlation coefficient value of test sample and various contrast signals come with threshold value
Compare, threshold value can be used as us and distinguish class signal another characteristic, see whether it belongs to corresponding classification.
The present invention is based on statistics to ask variance and relevance algorithms, first determines the interval range of attribute byte, then calculates phase
Closing property obtains the threshold value of identification various types of signal, and operation is succinct, and recognition efficiency is high, has a wide range of application, makes signal characteristic abstraction technology
Obtain huge development.
Description of the invention process also described by the mode of method and step specific function implementation procedure and its mutually
Relationship.For ease of description, having carried out special definition to the boundary of these method and steps and sequence in text.Make these functions and
Under the premise of its relationship can work normally, also their boundary of redefinable and sequence.But these weights to boundary and sequence
New definition is fallen among purport of the invention and the protection scope stated.
The present invention can also be implemented by computer program product, and program includes that can be realized the complete of the method for the present invention
Method of the invention may be implemented when it is installed in computer system in portion's feature.Computer program in this document is signified
: system can be made using any expression formula for one group of instruction that any program language, code or symbol are write, the instruction group
With information processing capability, to be directly realized by specific function, or after carrying out one or two following step specific function is realized
Can: a) it is converted into other Languages, coding or symbol；B) it reproduces in a different format.
Although the present invention be illustrated by specific embodiment, it will be appreciated by those skilled in the art that, do not departing from
In the case where the scope of the invention, various transformation and equivalent substitute can also be carried out to the present invention.In addition, being directed to particular condition or material
Material, can do various modifications to the present invention, without departing from the scope of the present invention.Therefore, the present invention is not limited to disclosed tool
Body embodiment, and should include the whole embodiments fallen within the scope of the appended claims.
Claims (7)
1. a kind of feature extracting method of industry control signal characterized by comprising
S1, acquisition simultaneously pre-process industry control signal to obtain physical signal sample；
S2, at least two groups physical signal sample is chosen from the physical signal sample to make different classes of physical signal sample
Attribute byte variogram；
S3, related coefficient curve is obtained based on the variogram；
S4, the characteristic threshold value that different classes of physical signal sample is obtained based on the related coefficient curve.
2. the feature extracting method of industry control signal according to claim 1, which is characterized in that the step S1 is further wrapped
It includes:
S11, the sort instructions signal for obtaining industrial control equipment；
The physical signal of the corresponding mark label of sort instructions signal described in S12, the excitation sort instructions signal and recording；
S13, the physical signal is normalized to obtain the physical signal sample.
3. the feature extracting method of industry control signal according to claim 2, which is characterized in that the step S11 is further
Include:
The instruction classification side of S111, the selected industrial control equipment that will be detected and the selected industrial control equipment that will be detected
Case；
S112, the industry control signal of the industrial control equipment is carried out by way of labelling based on described instruction classification schemes
Classification is to obtain sort instructions signal.
4. the feature extracting method of industry control signal according to claim 1, which is characterized in that the step S2 is further wrapped
It includes:
S21, selection at least two groups quantity is identical from the physical signal sample and covers the physical signal of codomain variation range
Sample；
S22, variance point is obtained according to following formula based on the physical signal sample, and variogram is done based on the variance point；
Wherein σ2Indicate the variance of the corresponding attribute byte of different types of physical signal sample, N indicates the physical signal sample of acquisition
This sample size, xiIndicate the amplitude size of physics sample of signal, μ indicates the mean value of the amplitude of physics sample of signal.
5. the feature extracting method of industry control signal according to claim 1, which is characterized in that the step S3 is further wrapped
It includes:
S31, the waveform that the physical signal sample where the attribute byte is obtained based on the variogram；
S32, the waveform mean value of each classification is calculated using the reference signal of the attribute byte as the category；
S33, related coefficient curve is obtained based on the reference signal and the range of waveforms of the physical signal sample.
6. the feature extracting method of industry control signal according to claim 5, which is characterized in that the step S33 is further
Including calculating the correlation coefficient value of the related coefficient curve based on following formula:
Wherein ρXYIndicate correlation coefficient value, Cov (X, Y) indicates the covariance of stochastic variable X and Y, and D (X) and D (Y) are respectively indicated
The variance of stochastic variable X and Y, μxIndicate the mean value of stochastic variable X, μyIndicate the mean value of stochastic variable Y, wherein stochastic variable X table
Show the reference signal, stochastic variable Y indicates the physical signal sample point in the range of waveforms of the physical signal sample.
7. a kind of computer readable storage medium, is stored thereon with computer program, which is characterized in that described program is by processor
The feature extracting method of industry control signal described in any one of -6 claims according to claim 1 is realized when execution.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810902497.4A CN109271855B (en) | 2018-08-09 | 2018-08-09 | Method for extracting characteristics of industrial control signals |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810902497.4A CN109271855B (en) | 2018-08-09 | 2018-08-09 | Method for extracting characteristics of industrial control signals |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109271855A true CN109271855A (en) | 2019-01-25 |
| CN109271855B CN109271855B (en) | 2021-07-09 |
Family
ID=65153484
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810902497.4A Active CN109271855B (en) | 2018-08-09 | 2018-08-09 | Method for extracting characteristics of industrial control signals |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109271855B (en) |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8019702B1 (en) * | 2007-12-07 | 2011-09-13 | Google Inc. | Supervised learning with multi-scale time intervals using a statistical classification model to classify unlabeled events |
| CN102780700A (en) * | 2012-07-27 | 2012-11-14 | 天津市亚安科技股份有限公司 | Method and device for realizing baud rate and protocol self-adaption in video monitoring field |
| US20130304739A1 (en) * | 2012-05-10 | 2013-11-14 | Samsung Electronics Co., Ltd. | Computing system with domain independence orientation mechanism and method of operation thereof |
| CN106209843A (en) * | 2016-07-12 | 2016-12-07 | 工业和信息化部电子工业标准化研究院 | A kind of data flow anomaly towards Modbus agreement analyzes method |
| CN107181724A (en) * | 2016-03-11 | 2017-09-19 | 华为技术有限公司 | A kind of recognition methods for cooperateing with stream, system and the server using this method |
| CN107689847A (en) * | 2017-09-30 | 2018-02-13 | 电子科技大学 | Mark auxiliary radio communication system signal method of sending and receiving |
| CN107809310A (en) * | 2017-11-15 | 2018-03-16 | 深圳大学 | Information detecting method, device and receiving device based on packet variance |
| CN108108712A (en) * | 2017-12-29 | 2018-06-01 | 哈尔滨工业大学 | A kind of Emitter Fingerprint feature extracting method based on variance dimension |
-
2018
- 2018-08-09 CN CN201810902497.4A patent/CN109271855B/en active Active
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8019702B1 (en) * | 2007-12-07 | 2011-09-13 | Google Inc. | Supervised learning with multi-scale time intervals using a statistical classification model to classify unlabeled events |
| US20130304739A1 (en) * | 2012-05-10 | 2013-11-14 | Samsung Electronics Co., Ltd. | Computing system with domain independence orientation mechanism and method of operation thereof |
| CN102780700A (en) * | 2012-07-27 | 2012-11-14 | 天津市亚安科技股份有限公司 | Method and device for realizing baud rate and protocol self-adaption in video monitoring field |
| CN107181724A (en) * | 2016-03-11 | 2017-09-19 | 华为技术有限公司 | A kind of recognition methods for cooperateing with stream, system and the server using this method |
| CN106209843A (en) * | 2016-07-12 | 2016-12-07 | 工业和信息化部电子工业标准化研究院 | A kind of data flow anomaly towards Modbus agreement analyzes method |
| CN107689847A (en) * | 2017-09-30 | 2018-02-13 | 电子科技大学 | Mark auxiliary radio communication system signal method of sending and receiving |
| CN107809310A (en) * | 2017-11-15 | 2018-03-16 | 深圳大学 | Information detecting method, device and receiving device based on packet variance |
| CN108108712A (en) * | 2017-12-29 | 2018-06-01 | 哈尔滨工业大学 | A kind of Emitter Fingerprint feature extracting method based on variance dimension |
Non-Patent Citations (2)
| Title |
|---|
| KHALED ASSALEH ET AL.: "Features extraction and analysis for classifying causable patterns in control charts", 《COMPUTERS & INDUSTRIAL ENGINEERING》 * |
| 杨凡: "雷达辐射源特征提取和个体识别", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109271855B (en) | 2021-07-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Petrovic | A comparison between the silhouette index and the davies-bouldin index in labelling ids clusters | |
| CN110826059B (en) | Method and device for defending black box attack facing malicious software image format detection model | |
| CN108030494B (en) | Electrocardiosignal error flag training sample recognition methods based on cross validation | |
| CN103488941B (en) | Hardware Trojan horse detection method and system | |
| CN106096411B (en) | A kind of Android malicious code family classification methods based on bytecode image clustering | |
| CN112491796B (en) | Intrusion detection and semantic decision tree quantitative interpretation method based on convolutional neural network | |
| CN108682007B (en) | JPEG image resampling automatic detection method based on depth random forest | |
| CN112729688B (en) | Oil-gas pipeline leakage detection method based on vibration and temperature double parameters | |
| CN105516127A (en) | Internal threat detection-oriented user cross-domain behavior pattern mining method | |
| CN111398679B (en) | Sub-synchronous oscillation identification and alarm method based on PMU (phasor measurement Unit) | |
| CN111915437A (en) | RNN-based anti-money laundering model training method, device, equipment and medium | |
| CN106792883A (en) | Sensor network abnormal deviation data examination method and system | |
| CN109976308A (en) | A kind of extracting method of the fault signature based on Laplce's score value and AP cluster | |
| CN109002810A (en) | Model evaluation method, Radar Signal Recognition method and corresponding intrument | |
| CN106874760A (en) | A kind of Android malicious code sorting techniques based on hierarchy type SimHash | |
| Gamboa-Medina et al. | Feature extraction in pressure signals for leak detection in water networks | |
| CN109088744A (en) | Powerline network abnormal intrusion detection method, device, equipment and storage medium | |
| CN108280357A (en) | Data leakage prevention method, system based on semantic feature extraction | |
| CN108900538A (en) | A kind of industry control signal detecting method and device | |
| CN106973051A (en) | Set up method, device, storage medium and the processor of detection Cyberthreat model | |
| Petrovic et al. | Labelling clusters in an intrusion detection system using a combination of clustering evaluation techniques | |
| CN109271855A (en) | A kind of feature extracting method of industry control signal | |
| Sipola et al. | Dimensionality reduction framework for detecting anomalies from network logs | |
| US9323987B2 (en) | Apparatus and method for detecting forgery/falsification of homepage | |
| CN113190851B (en) | Active learning method of malicious document detection model, electronic equipment and storage medium |
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 |