CN112526430A - Aircraft power supply characteristic frequency transient parameter calibration method - Google Patents
Aircraft power supply characteristic frequency transient parameter calibration method Download PDFInfo
- Publication number
- CN112526430A CN112526430A CN202011447473.8A CN202011447473A CN112526430A CN 112526430 A CN112526430 A CN 112526430A CN 202011447473 A CN202011447473 A CN 202011447473A CN 112526430 A CN112526430 A CN 112526430A
- Authority
- CN
- China
- Prior art keywords
- frequency
- hilbert transform
- power supply
- uncertainty
- original 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
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R35/00—Testing or calibrating of apparatus covered by the other groups of this subclass
- G01R35/005—Calibrating; Standards or reference devices, e.g. voltage or resistance standards, "golden" references
-
- 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/14—Fourier, Walsh or analogous domain transformations, e.g. Laplace, Hilbert, Karhunen-Loeve, transforms
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/30—Circuit design
- G06F30/36—Circuit design at the analogue level
- G06F30/367—Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods
-
- 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
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Mathematical Physics (AREA)
- Data Mining & Analysis (AREA)
- General Engineering & Computer Science (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Evolutionary Computation (AREA)
- Microelectronics & Electronic Packaging (AREA)
- Algebra (AREA)
- Geometry (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Measuring Frequencies, Analyzing Spectra (AREA)
Abstract
The invention discloses a method for calibrating frequency transient parameters of power supply characteristics of an aircraft power supply, and belongs to the technical field of electrical measurement. The invention evaluates the uncertainty evaluation introduced by a discrete Hilbert transform algorithm in the process of calibrating the frequency transient parameters, and the realization method comprises the following steps: the data acquisition system acquires the output frequency transient signal and then obtains the instantaneous frequency based on a discrete Hilbert transform algorithm. The uncertainty introduced by the error of the discrete hilbert transform due to the end-point effect caused by insufficient sampling rate when the frequency changes instantaneously is analyzed. The method comprises the steps of establishing a relation between instantaneous frequency and Hilbert transform of an original signal according to the principle that the instantaneous frequency is obtained through discrete Hilbert transform, calculating errors introduced by an end point effect under different sampling rates based on a numerical simulation technology, evaluating uncertainty, and perfecting uncertainty analysis of a frequency transient parameter calibration result of an aircraft power supply characteristic test system, so that the calibration result is more accurate and reliable.
Description
Technical Field
The invention relates to a method for calibrating frequency transient parameters of power supply characteristics of an aircraft power supply, and belongs to the technical field of electrical measurement.
Background
At present, power supply systems of large-scale military and civil aircraft at home and abroad generally adopt a constant-frequency alternating-current or variable-frequency alternating-current power supply system, and meanwhile, the research focus of electrical parameter testing of an aircraft power grid is gradually transferred to the research of transient parameters by steady-state parameters such as voltage, current, distortion and modulation, wherein the steady-state parameters comprise frequency transient parameters. According to the regulation of relevant standards, the calibration of frequency transient parameters requires the acquisition of the whole variation process by a data acquisition system, and the transient time and the transient frequency value are obtained by data analysis.
The mainstream instantaneous frequency extraction method comprises the modes of short-time Fourier transform, wavelet transform, Hilbert transform and the like, and the time resolution and the frequency resolution of the wavelet transform and the short-time Fourier transform algorithm cannot be considered at the same time, so that the Hilbert transform is a better choice for the situation of the text. Based on a discrete Hilbert transform algorithm in Matlab, a time domain signal is transformed into a frequency domain by means of an FFT algorithm, a Hilbert transform frequency spectrum is obtained through recombination of the frequency domain, and finally, the signal is transformed into a time domain by means of an IFFT algorithm, so that Hilbert transform is completed. According to the method, an end point effect is generated at a position where an original signal waveform is suddenly changed due to insufficient sampling rate, so that a discrete Hilbert transform result is deviated, and uncertainty is introduced.
Disclosure of Invention
In order to solve the problems of incomplete uncertainty evaluation process and unreliable calibration result in the conventional aircraft power supply characteristic frequency transient parameter calibration method, the invention aims to provide the aircraft power supply characteristic frequency transient parameter calibration method.
The purpose of the invention is realized by the following technical scheme:
the invention discloses a method for calibrating frequency transient parameters of power supply characteristics of an aircraft power supply, which is used for evaluating uncertainty introduced by a discrete Hilbert transform algorithm in the process of calibrating the frequency transient parameters of a power supply characteristic test system of the aircraft power supply and comprises the following steps: and acquiring the output frequency transient signal by a data acquisition system, and then obtaining the instantaneous frequency based on a discrete Hilbert transform algorithm. The uncertainty introduced by the error of the discrete hilbert transform due to the end-point effect caused by insufficient sampling rate when the frequency changes instantaneously is analyzed. The method comprises the steps of obtaining an instantaneous frequency according to discrete Hilbert transform, establishing a relation between the instantaneous frequency and Hilbert transform of an original signal and the original signal, calculating errors introduced by an end point effect under different sampling rates based on a numerical simulation technology, evaluating uncertainty, and perfecting uncertainty analysis of a frequency transient parameter calibration result of an aircraft power supply characteristic test system so that the calibration result is more accurate and reliable.
Preferably, the numerical simulation technology is preferably Matlab simulation technology.
The invention discloses a method for calibrating frequency transient parameters of power supply characteristics of an aircraft power supply.
Step one, determining the relation between the transient frequency and the original signal and the Hilbert transform thereof.
Firstly, obtaining transient frequency f (n), original signal s (n) obtained by collection and discrete Hilbert transform of the original signalThe relationship between:
in an analog system, the instantaneous frequency of the signal is obtained by equations (1), (2) and (3).
Wherein the content of the first and second substances,for the hilbert transform of the original signal s (t),to analyze the phase of the signal, f (t) is the time-frequency characteristic curve of the original signal.
The formula (2) is brought into the formula (3) and converted to obtain the formula (4).
Because the actual data is a sequence obtained by a data acquisition system, discrete Hilbert transform is also adopted in the processing process, and therefore, the formula (4) is discretized, and difference is used for replacing differentiation to obtain a formula (5), namely, the relation between the transient frequency and the original signal and the Hilbert transform is determined.
Wherein s (n) is the acquisition sequence of the original signal,is discrete Hilbert transform of the original signal, and f (n) is a time-frequency characteristic sequence of the original signal.
And step two, calculating uncertainty u caused by errors caused by insufficient sampling rate, namely realizing evaluation of uncertainty introduced by a discrete Hilbert transform algorithm, and perfecting uncertainty analysis of a frequency transient parameter calibration result of the aircraft power supply characteristic test system, so that the calibration result is more accurate and reliable.
As is clear from equation (5), the instantaneous frequency can be calculated from the original signal and the hilbert transform sequence of the original signal, and the uncertainty u can be obtained by the following procedure.
The sampling rate of the signal is fs, the instantaneous frequency f (n) of the signal s (n) is calculated, and the instantaneous frequency value f of the transient moment is recorded0;
Let the new sampling rate be 2kFs (k ═ 1,2,3, …), and new instantaneous frequency values are calculated until f is foundkIf the difference between the current frequency value and the instantaneous frequency value obtained by the previous updating is smaller than a preset threshold value, the error introduced by the end point effect cannot be obviously reduced by increasing the sampling rate;
let Δ f be f0-fkCalculating the uncertainty contribution of the part according to uniform distribution
And the uncertainty u is obtained, namely, the uncertainty introduced by the discrete Hilbert transform algorithm is evaluated, and the uncertainty analysis of the frequency transient parameter calibration result of the aircraft power supply characteristic test system is perfected, so that the calibration result is more accurate and reliable.
Preferably, the difference in instantaneous frequency values being less than a predetermined threshold value means
Has the advantages that:
the invention discloses a method for calibrating the frequency transient parameters of the power supply characteristics of an aircraft power supply, which is characterized in that errors introduced by an end point effect are calculated under different sampling rates by a numerical simulation technology, instantaneous frequency can be calculated by a Hilbert transform sequence of an original signal and an original signal, uncertainty introduced by a discrete Hilbert transform algorithm is evaluated, uncertainty analysis of the frequency transient parameter calibration result of a power supply characteristic test system of the aircraft power supply is perfected, uncertainty introduced by each link of the calibration method can be more comprehensively evaluated in the calibration process of the frequency transient parameters of the power supply characteristics of the aircraft power supply, and the calibration result is more accurate and reliable.
Drawings
FIG. 1 is a flow chart of a frequency transient parameter calibration process of an aircraft power supply characteristic test system.
Detailed Description
The following examples illustrate the invention in detail.
As shown in fig. 1, in the method for calibrating the frequency transient parameter of the power supply characteristic of the aircraft disclosed in the present embodiment, the primary source of uncertainty u of the frequency transient parameter calibration result is an error introduced by an "end effect" due to insufficient sampling rate at a sudden change of the waveform.
Step one, determining the relation between the transient frequency and the original signal and the Hilbert transform thereof.
And simulating an original signal s (n) according to the frequency transient signal characteristic of the power supply characteristic of the airplane power supply, wherein the process is as follows. Setting the sampling rate fs of the signal to be 10MS/s and the total sampling point number to be 10M according to the highest acquisition rate of the used high-precision data acquisition system, creating two sine wave sequences s1 and s2, wherein the frequencies are respectively the lower limit and the upper limit f of the frequency conversion alternating current power supply system in the aircraft power supply system1360Hz and f2=800Hz。
Signal s is defined as equation (8).
And step two, calculating uncertainty u caused by errors caused by insufficient sampling rate, namely realizing evaluation of uncertainty introduced by a discrete Hilbert transform algorithm, and perfecting uncertainty analysis of a frequency transient parameter calibration result of the aircraft power supply characteristic test system, so that the calibration result is more accurate and reliable.
Calculating the instantaneous frequency f (n) of the signal s (n) according to the formula (5), and recordingInstantaneous frequency value f of time0=800.835Hz;
Let the new sampling rate be 2kFs (k ═ 1,2,3, …), and new instantaneous frequency values are calculated until f is foundkThe difference of the instantaneous frequency value obtained in the previous step is less than 0.0001 time;
let Δ f be f0-fk0.825Hz, the partial uncertainty contribution is calculated in a uniform distribution, then
And the uncertainty u is obtained, namely, the uncertainty introduced by the discrete Hilbert transform algorithm is evaluated, and the uncertainty analysis of the frequency transient parameter calibration result of the aircraft power supply characteristic test system is perfected, so that the calibration result is more accurate and reliable.
The above detailed description is intended to illustrate the objects, aspects and advantages of the present invention, and it should be understood that the above detailed description is only exemplary of the present invention and is not intended to limit the scope of the present invention, and any modifications, equivalents, improvements and the like made within the spirit and principle of the present invention should be included in the scope of the present invention.
Claims (4)
1. A method for calibrating frequency transient parameters of power supply characteristics of an aircraft power supply is characterized by comprising the following steps: the method for evaluating the uncertainty introduced by a discrete Hilbert transform algorithm in the frequency transient parameter calibration process of the aircraft power supply characteristic test system comprises the following steps: acquiring the output frequency transient signal by a data acquisition system and then obtaining instantaneous frequency based on a discrete Hilbert transform algorithm; analyzing the uncertainty introduced by errors caused by an end point effect due to insufficient sampling rate when the frequency of the discrete Hilbert transform is instantaneously changed; the method comprises the steps of obtaining an instantaneous frequency according to discrete Hilbert transform, establishing a relation between the instantaneous frequency and Hilbert transform of an original signal and the original signal, calculating errors introduced by an end point effect under different sampling rates based on a numerical simulation technology, evaluating uncertainty, and perfecting uncertainty analysis of a frequency transient parameter calibration result of an aircraft power supply characteristic test system so that the calibration result is more accurate and reliable.
2. An aircraft power supply characteristic frequency transient parameter calibration method as defined in claim 1, wherein: the numerical simulation technology selects Matlab simulation technology.
3. An aircraft power supply characteristic frequency transient parameter calibration method as defined in claim 1 or 2, characterized in that: the realization method is that,
determining the relation between transient frequency and original signal and Hilbert transform thereof;
firstly, obtaining transient frequency f (n), original signal s (n) obtained by collection and discrete Hilbert transform of the original signalThe relationship between:
in an analog system, the instantaneous frequency of the signal is given by equations (1), (2) and (3);
wherein the content of the first and second substances,for the hilbert transform of the original signal s (t),to analyze the phase of the signal, f (t) is the time-frequency characteristic curve of the original signal;
carrying the formula (2) into the formula (3), and converting to obtain a formula (4);
because the actual data is a sequence obtained by a data acquisition system, discrete Hilbert transform is also adopted in the processing process, and therefore, the formula (4) is discretized, and differential is replaced by difference to obtain a formula (5), namely, the relation between the transient frequency and the original signal and the Hilbert transform is determined;
wherein s (n) is the acquisition sequence of the original signal,discrete Hilbert transform of the original signal, f (n) time-frequency characteristic sequence of the original signal;
step two, calculating uncertainty u caused by errors caused by insufficient sampling rate, namely realizing evaluation of uncertainty introduced by a discrete Hilbert transform algorithm, and perfecting uncertainty analysis of a frequency transient parameter calibration result of an aircraft power supply characteristic test system, so that the calibration result is more accurate and reliable;
as shown in the formula (5), the instantaneous frequency can be calculated from the original signal and the hilbert transform sequence of the original signal, and the uncertainty u can be obtained by the following steps;
the sampling rate of the signal is fs, the instantaneous frequency f (n) of the signal s (n) is calculated, and the instantaneous frequency value f of the transient moment is recorded0;
Let the new sampling rate be 2kFs (k ═ 1,2,3, …), and new instantaneous frequency values are calculated until f is foundkIf the difference between the current frequency value and the instantaneous frequency value obtained by the previous updating is smaller than a preset threshold value, the error introduced by the end point effect cannot be obviously reduced by increasing the sampling rate;
let Δ f be f0-fkCalculating the uncertainty contribution of the part according to uniform distribution
And the uncertainty u is obtained, namely, the uncertainty introduced by the discrete Hilbert transform algorithm is evaluated, and the uncertainty analysis of the frequency transient parameter calibration result of the aircraft power supply characteristic test system is perfected, so that the calibration result is more accurate and reliable.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011447473.8A CN112526430B (en) | 2020-12-09 | 2020-12-09 | Aircraft power supply characteristic frequency transient parameter calibration method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011447473.8A CN112526430B (en) | 2020-12-09 | 2020-12-09 | Aircraft power supply characteristic frequency transient parameter calibration method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112526430A true CN112526430A (en) | 2021-03-19 |
CN112526430B CN112526430B (en) | 2022-08-09 |
Family
ID=75000479
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011447473.8A Active CN112526430B (en) | 2020-12-09 | 2020-12-09 | Aircraft power supply characteristic frequency transient parameter calibration method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112526430B (en) |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6901353B1 (en) * | 2003-07-08 | 2005-05-31 | The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration | Computing Instantaneous Frequency by normalizing Hilbert Transform |
EP1804334A1 (en) * | 2005-12-27 | 2007-07-04 | Nederlandse Organisatie voor Toegepast-Natuuurwetenschappelijk Onderzoek TNO | Phased array antenna apparatus |
CN101887407A (en) * | 2010-07-16 | 2010-11-17 | 哈尔滨工业大学 | Hilbert-Huang transform-based equipment or system built-in test signal characteristic extraction method |
CN101988935A (en) * | 2009-12-18 | 2011-03-23 | 西安电子科技大学 | Digital down converter-Hilbert-Huang transform based instantaneous frequency measuring method |
CN106224224A (en) * | 2016-07-13 | 2016-12-14 | 北京航空航天大学 | A kind of based on Hilbert-Huang transform and quality the Hydraulic pump fault feature extracting method away from entropy |
CN106650218A (en) * | 2016-10-19 | 2017-05-10 | 上海电机学院 | Harmonic analysis method based on complementary ensemble empirical mode decomposition (CEEMD) algorithm and Hilbert transform |
CN109815581A (en) * | 2019-01-21 | 2019-05-28 | 西北工业大学 | The appraisal procedure of aircraft electrical supply parameters test analysis algorithm software |
CN110187310A (en) * | 2019-05-13 | 2019-08-30 | 北京遥感设备研究所 | A kind of radar equipment LFM signal transient frequency curve approximating method and system |
CN110865238A (en) * | 2019-11-22 | 2020-03-06 | 北京东方计量测试研究所 | Alternating current resistance measurement method and device based on quasi-harmonic model sampling algorithm |
CN110927822A (en) * | 2019-12-03 | 2020-03-27 | 吉林大学 | Method for evaluating accuracy of magnetic gradient tensor obtained by Hilbert transform algorithm |
-
2020
- 2020-12-09 CN CN202011447473.8A patent/CN112526430B/en active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6901353B1 (en) * | 2003-07-08 | 2005-05-31 | The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration | Computing Instantaneous Frequency by normalizing Hilbert Transform |
EP1804334A1 (en) * | 2005-12-27 | 2007-07-04 | Nederlandse Organisatie voor Toegepast-Natuuurwetenschappelijk Onderzoek TNO | Phased array antenna apparatus |
CN101988935A (en) * | 2009-12-18 | 2011-03-23 | 西安电子科技大学 | Digital down converter-Hilbert-Huang transform based instantaneous frequency measuring method |
CN101887407A (en) * | 2010-07-16 | 2010-11-17 | 哈尔滨工业大学 | Hilbert-Huang transform-based equipment or system built-in test signal characteristic extraction method |
CN106224224A (en) * | 2016-07-13 | 2016-12-14 | 北京航空航天大学 | A kind of based on Hilbert-Huang transform and quality the Hydraulic pump fault feature extracting method away from entropy |
CN106650218A (en) * | 2016-10-19 | 2017-05-10 | 上海电机学院 | Harmonic analysis method based on complementary ensemble empirical mode decomposition (CEEMD) algorithm and Hilbert transform |
CN109815581A (en) * | 2019-01-21 | 2019-05-28 | 西北工业大学 | The appraisal procedure of aircraft electrical supply parameters test analysis algorithm software |
CN110187310A (en) * | 2019-05-13 | 2019-08-30 | 北京遥感设备研究所 | A kind of radar equipment LFM signal transient frequency curve approximating method and system |
CN110865238A (en) * | 2019-11-22 | 2020-03-06 | 北京东方计量测试研究所 | Alternating current resistance measurement method and device based on quasi-harmonic model sampling algorithm |
CN110927822A (en) * | 2019-12-03 | 2020-03-27 | 吉林大学 | Method for evaluating accuracy of magnetic gradient tensor obtained by Hilbert transform algorithm |
Non-Patent Citations (5)
Title |
---|
李小舟等: "基于希尔伯特变换的信号解调算法及其在飞机供电特性参数测试***中的应用", 《计量学报》 * |
王建强 等: "飞机变频交流电源供电特性频率瞬变参数测量方法研究", 《中北大学学报(自然科学版)》 * |
王慧等: "测试高分辨率ADC有效位数的HHT方法", 《电子学报》 * |
王文健等: "飞机直流供电特性测试***畸变频谱参数校准方法", 《计测技术》 * |
王静等: "25 Hz相敏轨道电路相位差测量非整周期采样误差校正算法", 《计算机应用》 * |
Also Published As
Publication number | Publication date |
---|---|
CN112526430B (en) | 2022-08-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Gutiérrez et al. | Flicker measurement system for wind turbine certification | |
CN110161421B (en) | Method for reconstructing battery impedance in set frequency range on line | |
CN107402334B (en) | Electric system time/supersynchronous coupled impedance model frequency domain identification method and system | |
CN109782206B (en) | Calibration compensation method for broadband transient voltage measuring device | |
CN108535613B (en) | Voltage flicker parameter detection method based on combined window function | |
CN105954695B (en) | Synchronization-based homogeneous sensor mutation parameter identification method and device | |
CN111289796B (en) | Method for detecting subsynchronous oscillation of high-proportion renewable energy power system | |
CN103257273A (en) | Method for measuring phase difference of signals with same frequency period | |
CN104410360A (en) | Safe operation method of photovoltaic power generation system, training method for artificial neural network and real-time detection method in safe operation method, and real-time detection device | |
CN110967658B (en) | Analog input merging unit calibrator tracing method based on digital differential method | |
CN110389312B (en) | Calibrator phasor measurement method suitable for field PMU test | |
CN108896944A (en) | A kind of synchronous measuring apparatus laboratory investment instrument and its synchronous phasor measuring method | |
CN112986775A (en) | Active calibration method for partial discharge pulse phase | |
US7206340B2 (en) | Characterizing jitter of repetitive patterns | |
CN112526430B (en) | Aircraft power supply characteristic frequency transient parameter calibration method | |
CN111812573B (en) | Method and system for analyzing harmonic test error characteristics of extra-high voltage CVT | |
CN112505391B (en) | Frequency-adaptive alternating current signal effective value acquisition method | |
CN109283391A (en) | A kind of synchronized phasor method for measurement based on nonlinear fitting | |
KR101807280B1 (en) | Apparatus and Method for detecting accidents in Power Systems using PMU signal | |
CN106053936B (en) | A kind of method and system obtaining electrical signal instantaneous frequency | |
Cho et al. | Oscillation recognition using a geometric feature extraction process based on periodic time-series approximation | |
Chen | A calibration test platform of power quality instruments for grid integration of wind energy system | |
CN108982954B (en) | Method and system for calculating phase voltage amplitude and phase suitable for feeder line terminal | |
CN113884966B (en) | AC/DC sensor broadband calibrator and performance detection calibration method thereof | |
CN104360297B (en) | A kind of PMU dynamic property multi objective method of testings compared based on instantaneous value |
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 |