CN107167658A - 一种高精度抗干扰的电力系统基波频率及相位测量方法 - Google Patents

一种高精度抗干扰的电力系统基波频率及相位测量方法 Download PDF

Info

Publication number
CN107167658A
CN107167658A CN201710331779.9A CN201710331779A CN107167658A CN 107167658 A CN107167658 A CN 107167658A CN 201710331779 A CN201710331779 A CN 201710331779A CN 107167658 A CN107167658 A CN 107167658A
Authority
CN
China
Prior art keywords
mrow
msub
mfrac
high accuracy
frequency
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
Application number
CN201710331779.9A
Other languages
English (en)
Other versions
CN107167658B (zh
Inventor
李帅
陈福锋
赵谦
龚世敏
孙伟
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nanjing SAC Automation Co Ltd
Original Assignee
Nanjing SAC Automation Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Nanjing SAC Automation Co Ltd filed Critical Nanjing SAC Automation Co Ltd
Priority to CN201710331779.9A priority Critical patent/CN107167658B/zh
Publication of CN107167658A publication Critical patent/CN107167658A/zh
Application granted granted Critical
Publication of CN107167658B publication Critical patent/CN107167658B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R23/00Arrangements for measuring frequencies; Arrangements for analysing frequency spectra
    • G01R23/02Arrangements for measuring frequency, e.g. pulse repetition rate; Arrangements for measuring period of current or voltage
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R25/00Arrangements for measuring phase angle between a voltage and a current or between voltages or currents

Abstract

本发明公开了一种高精度抗干扰的电力系统基波频率及相位测量方法,通过复向量提取算法获取基波的实部ure,基波的虚部uim;构造复向量并计算相角θ1;结合上次相角计算结果θ2,计算基波频率f;对θ1进行延时相位补偿,计算基波相角θ;本发明采用复向量提取算法直接获取基波的实部、虚部,避免了数据干扰、谐波、直流分量等因素的影响;同时,参数设计时,实部、虚部具有相同的增益与延时,保证了计算结果与实际波形的一致性,因而具有高精度、抗干扰的特点。

Description

一种高精度抗干扰的电力系统基波频率及相位测量方法
技术领域
本发明涉及一种高精度抗干扰的电力系统基波频率及相位测量方法,属于电力系统测量技术领域。
背景技术
基波的频率、相位等信息作为电力系统交流采样的重要特征量,与保护、测量、控制等环节密切相关。目前,现有的测量方法主要基于傅里叶分析、过零点检测等,存在一定的局限。傅里叶分析算法较为成熟,但由于频谱泄露等问题,存在一定的测量偏差;过零点算法原理简单,但当叠加谐波或直流分量时,该算法误差较大。此外,电力设备的运行工况较为复杂,采样数据极易受到干扰,直接影响频率、相位的测量。因此,构造一种高精度抗干扰的基波频率及相位测量方法,是当前急需解决的问题。
发明内容
目的:为了克服现有技术中存在的不足,本发明提供一种高精度抗干扰的电力系统基波频率及相位测量方法。
技术方案:为解决上述技术问题,本发明采用的技术方案为:
一种高精度抗干扰的电力系统基波频率及相位测量方法,包括如下步骤:
步骤1,确定复向量提取算法的初始长度N1及截断长度N;
步骤2,设定频率参数fc1、fc2
步骤3,根据公式1,计算实部提取算法的初始迭代系数an,其中fs为采样频率,n为迭代系数的序号;
步骤4,根据公式2,计算虚部提取算法的初始迭代系数bn
步骤5,通过窗函数加权初始迭代系数an、bn
步骤6,根据公式3、4,截取初始迭代公式,得出迭代系数cn、dn
步骤7,判断实部、虚部的提取公式的幅值增益是否满足如下指标:
(a)40~60Hz的幅值增益的相对偏差≤0.5%;
(b)直流分量及100Hz以上的幅值增益≤0.05;
若不满足,重复步骤2,调整频率参数,重新计算初始迭代系数;
步骤8,根据公式5、6,计算基波的实部ure、虚部uim,其中xn为输入波形的采样序列;
步骤9,根据公式7、8,构造复向量并计算相角θ1
步骤10,设上次相角计算的结果为θ2,每个工频周波即20ms计算一次相角,则基波频率f,如公式9所示
步骤11,对θ1进行延时相位补偿,则基波相位θ如下
步骤12,数据更新,令θ2=θ1,重复步骤8,进行下一次运算。
有益效果:本发明提供的一种高精度抗干扰的电力系统基波频率及相位测量方法,通过复向量提取算法直接获取基波的实部、虚部,避免了数据干扰、谐波、直流分量等因素的影响;并且参数设计时,实部、虚部具有相同的增益与延时,确保了计算结果与原始数据一致,因此本方法计算精度高、抗干扰能力强,具有良好的应用前景。
附图说明
图1为本发明基波频率及相位测量方法的原理框图;
图2为实部、虚部提取算法的参数设计流程;
图3为满足指标的幅频特性图;
图4为仿真模型结构图;
图5为仿真结果示意图。
具体实施方式
下面结合附图对本发明作更进一步的说明。
如图1、图2所示,本发明的电力系统基波频率及相位测量方法,包括以下步骤,其中步骤1~7为实部、虚部提取算法的参数计算过程,步骤8~12为频率相位计算过程。
步骤1,确定复向量提取算法的初始长度N1及截断长度N;
步骤2,设定频率参数fc1、fc2
步骤3,根据公式1,计算实部提取算法的初始迭代系数an,其中fs为采样频率,n为迭代系数的序号;
步骤4,根据公式2,计算虚部提取算法的初始迭代系数bn
步骤5,通过窗函数加权初始迭代系数an、bn
步骤6,根据公式3、4,截取初始迭代公式,得出迭代系数cn、dn
步骤7,判断实部、虚部的提取公式的幅值增益是否满足如下指标:
(a)40~60Hz的幅值增益的相对偏差≤0.5%;
(b)直流分量及100Hz以上的幅值增益≤0.05;
若不满足,重复步骤2,调整频率参数,重新计算初始迭代系数;
步骤8,根据公式5、6,计算基波的实部ure、虚部uim,其中xn为输入波形的采样序列;
步骤9,根据公式7、8,构造复向量并计算相角θ1
步骤10,设上次相角计算的结果为θ2,每个工频周波即20ms计算一次相角,则基波频率f,如公式9所示
步骤11,对θ1进行延时相位补偿,则基波相位θ如下
步骤12,数据更新,令θ2=θ1,重复步骤8,进行下一次运算。
下面给出算法实例。针对采样频率为4000Hz的系统,按步骤1-7设计了一组实部、虚部提取参数,如下表所示:
如图3所示,其幅频特性满足步骤7的指标要求。
如图4所示,根据步骤8-11,用MATLAB/Simulink搭建了频率及相角计算的仿真模型。
为了验证算法的抗干扰能力及精确性,输入波形采用基波叠加谐波的方式,如下式所示。其中,基波频率设为55Hz,幅值为100;谐波频率设为100Hz及150Hz,谐波幅值为20,如下式所示。
y=100sin(2πt·55)+20sin(2πt·100+π/3)+20sin(2πt·150+π/2)(11)
如图5所示,给出了0.1~0.3s的仿真结果。计算出的频率为55Hz,与实际一致;相角的计算结果亦与基波波形吻合。
以上所述仅是本发明的优选实施方式,应当指出:对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。

Claims (7)

1.一种高精度抗干扰的电力系统基波频率及相位测量方法,其特征在于:包括如下步骤:
步骤1,确定复向量提取算法的初始长度N1及截断长度N;
步骤2,设定频率参数fc1、fc2
步骤3,计算实部提取算法的初始迭代系数an
步骤4,计算虚部提取算法的初始迭代系数bn
步骤5,通过窗函数加权初始迭代系数an、bn
步骤6,根据截取初始迭代公式,得出迭代系数cn、dn
步骤7,计算基波的实部ure、虚部uim
步骤8,构造复向量并计算相角θ1
<mrow> <mover> <mi>u</mi> <mo>&amp;RightArrow;</mo> </mover> <mo>=</mo> <msub> <mi>u</mi> <mrow> <mi>r</mi> <mi>e</mi> </mrow> </msub> <mo>+</mo> <mi>j</mi> <mo>&amp;CenterDot;</mo> <msub> <mi>u</mi> <mrow> <mi>i</mi> <mi>m</mi> </mrow> </msub> </mrow>
<mrow> <msub> <mi>&amp;theta;</mi> <mn>1</mn> </msub> <mo>=</mo> <mi>a</mi> <mi>n</mi> <mi>g</mi> <mi>l</mi> <mi>e</mi> <mrow> <mo>(</mo> <mover> <mi>u</mi> <mo>&amp;RightArrow;</mo> </mover> <mo>)</mo> </mrow> </mrow>
步骤9,设上次相角计算的结果为θ2,每个工频周波计算一次相角,则基波频率f计算如下:
<mrow> <mi>f</mi> <mo>=</mo> <mn>50</mn> <mo>&amp;CenterDot;</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>+</mo> <mfrac> <mrow> <msub> <mi>&amp;theta;</mi> <mn>1</mn> </msub> <mo>-</mo> <msub> <mi>&amp;theta;</mi> <mn>2</mn> </msub> </mrow> <mrow> <mn>2</mn> <mi>&amp;pi;</mi> </mrow> </mfrac> <mo>)</mo> </mrow> </mrow>
步骤10,对θ1进行延时相位补偿,则基波相位θ计算如下,其中fs为采样频率;
<mrow> <mi>&amp;theta;</mi> <mo>=</mo> <msub> <mi>&amp;theta;</mi> <mn>1</mn> </msub> <mo>+</mo> <mfrac> <mrow> <mi>f</mi> <mo>&amp;CenterDot;</mo> <mrow> <mo>(</mo> <mi>N</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <mo>&amp;CenterDot;</mo> <mi>&amp;pi;</mi> </mrow> <msub> <mi>f</mi> <mi>S</mi> </msub> </mfrac> </mrow>
步骤11,数据更新,令θ2=θ1,重复步骤7,进行下一次运算。
2.根据权利要求1所述的一种高精度抗干扰的电力系统基波频率及相位测量方法,其特征在于:所述计算实部提取算法的初始迭代系数an如下:
<mrow> <msub> <mi>a</mi> <mi>n</mi> </msub> <mo>=</mo> <mfrac> <mn>1</mn> <mrow> <mi>&amp;pi;</mi> <mo>&amp;lsqb;</mo> <mi>n</mi> <mo>-</mo> <mrow> <mo>(</mo> <msub> <mi>N</mi> <mn>1</mn> </msub> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <mo>/</mo> <mn>2</mn> <mo>&amp;rsqb;</mo> </mrow> </mfrac> <mrow> <mo>(</mo> <mi>cos</mi> <mfrac> <mrow> <mn>2</mn> <msub> <mi>&amp;pi;f</mi> <mrow> <mi>c</mi> <mn>1</mn> </mrow> </msub> <mo>&amp;lsqb;</mo> <mi>n</mi> <mo>-</mo> <mrow> <mo>(</mo> <msub> <mi>N</mi> <mn>1</mn> </msub> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <mo>/</mo> <mn>2</mn> <mo>&amp;rsqb;</mo> </mrow> <msub> <mi>f</mi> <mi>s</mi> </msub> </mfrac> <mo>-</mo> <mi>cos</mi> <mfrac> <mrow> <mn>2</mn> <msub> <mi>&amp;pi;f</mi> <mrow> <mi>c</mi> <mn>2</mn> </mrow> </msub> <mo>&amp;lsqb;</mo> <mi>n</mi> <mo>-</mo> <mrow> <mo>(</mo> <msub> <mi>N</mi> <mn>1</mn> </msub> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <mo>/</mo> <mn>2</mn> <mo>&amp;rsqb;</mo> </mrow> <msub> <mi>f</mi> <mi>s</mi> </msub> </mfrac> <mo>)</mo> </mrow> <mo>.</mo> </mrow>
3.根据权利要求1所述的一种高精度抗干扰的电力系统基波频率及相位测量方法,其特征在于:所述计算虚部提取算法的初始迭代系数bn如下:
<mrow> <msub> <mi>b</mi> <mi>n</mi> </msub> <mo>=</mo> <mfrac> <mn>1</mn> <mrow> <mi>&amp;pi;</mi> <mo>&amp;lsqb;</mo> <mi>n</mi> <mo>-</mo> <mrow> <mo>(</mo> <msub> <mi>N</mi> <mn>1</mn> </msub> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <mo>/</mo> <mn>2</mn> <mo>&amp;rsqb;</mo> </mrow> </mfrac> <mrow> <mo>(</mo> <mi>sin</mi> <mfrac> <mrow> <mn>2</mn> <msub> <mi>&amp;pi;f</mi> <mrow> <mi>c</mi> <mn>1</mn> </mrow> </msub> <mo>&amp;lsqb;</mo> <mi>n</mi> <mo>-</mo> <mrow> <mo>(</mo> <msub> <mi>N</mi> <mn>1</mn> </msub> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <mo>/</mo> <mn>2</mn> <mo>&amp;rsqb;</mo> </mrow> <msub> <mi>f</mi> <mi>s</mi> </msub> </mfrac> <mo>-</mo> <mi>sin</mi> <mfrac> <mrow> <mn>2</mn> <msub> <mi>&amp;pi;f</mi> <mrow> <mi>c</mi> <mn>2</mn> </mrow> </msub> <mo>&amp;lsqb;</mo> <mi>n</mi> <mo>-</mo> <mrow> <mo>(</mo> <msub> <mi>N</mi> <mn>1</mn> </msub> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <mo>/</mo> <mn>2</mn> <mo>&amp;rsqb;</mo> </mrow> <msub> <mi>f</mi> <mi>s</mi> </msub> </mfrac> <mo>)</mo> </mrow> <mo>.</mo> </mrow>
4.根据权利要求1所述的一种高精度抗干扰的电力系统基波频率及相位测量方法,其特征在于:所述截取初始迭代公式如下:
<mrow> <msub> <mi>c</mi> <mi>n</mi> </msub> <mo>=</mo> <msub> <mi>a</mi> <mrow> <mi>n</mi> <mo>+</mo> <mrow> <mo>(</mo> <msub> <mi>N</mi> <mn>1</mn> </msub> <mo>-</mo> <mi>N</mi> <mo>)</mo> </mrow> <mo>/</mo> <mn>2</mn> </mrow> </msub> </mrow>
<mrow> <msub> <mi>c</mi> <mi>n</mi> </msub> <mo>=</mo> <msub> <mi>a</mi> <mrow> <mi>n</mi> <mo>+</mo> <mrow> <mo>(</mo> <msub> <mi>N</mi> <mn>1</mn> </msub> <mo>-</mo> <mi>N</mi> <mo>)</mo> </mrow> <mo>/</mo> <mn>2</mn> </mrow> </msub> <mo>.</mo> </mrow>
5.根据权利要求1所述的一种高精度抗干扰的电力系统基波频率及相位测量方法,其特征在于:所述步骤6后还包括判断实部、虚部的提取公式的幅值增益是否满足如下指标:
(a)40~60Hz的幅值增益的相对偏差≤0.5%;
(b)直流分量及100Hz以上的幅值增益≤0.05;
若不满足,重复步骤2,调整频率参数,重新计算初始迭代系数。
6.根据权利要求1所述的一种高精度抗干扰的电力系统基波频率及相位测量方法,其特征在于:所述计算基波的实部ure、虚部uim,其中xn为输入波形的采样序列;
<mrow> <msub> <mi>u</mi> <mrow> <mi>r</mi> <mi>e</mi> </mrow> </msub> <mo>=</mo> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>n</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>N</mi> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <msub> <mi>c</mi> <mi>n</mi> </msub> <msub> <mi>x</mi> <mi>n</mi> </msub> </mrow>
<mrow> <msub> <mi>u</mi> <mrow> <mi>i</mi> <mi>m</mi> </mrow> </msub> <mo>=</mo> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>n</mi> <mo>=</mo> <mn>0</mn> </mrow> <mrow> <mi>N</mi> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <msub> <mi>d</mi> <mi>n</mi> </msub> <msub> <mi>x</mi> <mi>n</mi> </msub> <mo>.</mo> </mrow>
7.根据权利要求1所述的一种高精度抗干扰的电力系统基波频率及相位测量方法,其特征在于:所述步骤9中频率及相角计算的周期为一个工频周波即20ms。
CN201710331779.9A 2017-05-11 2017-05-11 一种高精度抗干扰的电力系统基波频率及相位测量方法 Active CN107167658B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710331779.9A CN107167658B (zh) 2017-05-11 2017-05-11 一种高精度抗干扰的电力系统基波频率及相位测量方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710331779.9A CN107167658B (zh) 2017-05-11 2017-05-11 一种高精度抗干扰的电力系统基波频率及相位测量方法

Publications (2)

Publication Number Publication Date
CN107167658A true CN107167658A (zh) 2017-09-15
CN107167658B CN107167658B (zh) 2019-06-28

Family

ID=59815623

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710331779.9A Active CN107167658B (zh) 2017-05-11 2017-05-11 一种高精度抗干扰的电力系统基波频率及相位测量方法

Country Status (1)

Country Link
CN (1) CN107167658B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107797452A (zh) * 2017-10-25 2018-03-13 北京富力通达科技有限公司 多通道协调加载系统相位和幅值的控制方法
CN108333426A (zh) * 2017-12-25 2018-07-27 南京丰道电力科技有限公司 基于傅氏算法的电力系统频率测量方法
CN109884388A (zh) * 2019-04-02 2019-06-14 华中科技大学 一种基于半周期移相法的电网频率计量、测量装置及方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101216512A (zh) * 2007-12-29 2008-07-09 湖南大学 一种非正弦周期信号实时高精度检测方法
CN101726690A (zh) * 2009-11-30 2010-06-09 上海电力学院 电力系统短路全电流各分量分解方法
CN102236048A (zh) * 2011-06-16 2011-11-09 国网电力科学研究院 一种电力系统相量频率测量方法
CN102798748A (zh) * 2012-07-13 2012-11-28 中冶南方工程技术有限公司 一种基于迭代傅里叶变换计算交流电压信号幅值和相位的方法
CN103543331A (zh) * 2013-10-24 2014-01-29 佟晓白 一种计算电信号谐波和间谐波的方法
CN104833851A (zh) * 2015-04-24 2015-08-12 合肥工业大学 基于分布式相关卡尔曼滤波的电力系统谐波估计方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101216512A (zh) * 2007-12-29 2008-07-09 湖南大学 一种非正弦周期信号实时高精度检测方法
CN101726690A (zh) * 2009-11-30 2010-06-09 上海电力学院 电力系统短路全电流各分量分解方法
CN102236048A (zh) * 2011-06-16 2011-11-09 国网电力科学研究院 一种电力系统相量频率测量方法
CN102798748A (zh) * 2012-07-13 2012-11-28 中冶南方工程技术有限公司 一种基于迭代傅里叶变换计算交流电压信号幅值和相位的方法
CN103543331A (zh) * 2013-10-24 2014-01-29 佟晓白 一种计算电信号谐波和间谐波的方法
CN104833851A (zh) * 2015-04-24 2015-08-12 合肥工业大学 基于分布式相关卡尔曼滤波的电力系统谐波估计方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
朱骏 等: "基波幅值、相位和频率的软件检测技术", 《电工技术》 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107797452A (zh) * 2017-10-25 2018-03-13 北京富力通达科技有限公司 多通道协调加载系统相位和幅值的控制方法
CN108333426A (zh) * 2017-12-25 2018-07-27 南京丰道电力科技有限公司 基于傅氏算法的电力系统频率测量方法
CN109884388A (zh) * 2019-04-02 2019-06-14 华中科技大学 一种基于半周期移相法的电网频率计量、测量装置及方法
CN109884388B (zh) * 2019-04-02 2021-05-18 华中科技大学 一种基于半周期移相法的电网频率计量、测量装置及方法

Also Published As

Publication number Publication date
CN107167658B (zh) 2019-06-28

Similar Documents

Publication Publication Date Title
CN103018534B (zh) 确定谐波电压的方法及系统
CN103401238B (zh) 一种基于总体测辨法的电力负荷建模方法
CN102868183B (zh) 单相并网逆变器的基于多谐振滑模面的滑模变结构控制方法
CN107590317A (zh) 一种计及模型参数不确定性的发电机动态估计方法
CN107016489A (zh) 一种电力系统抗差状态估计方法和装置
CN103838959A (zh) 偏最小二乘回归应用于配电网谐波源定位与检测中的方法
CN108155643B (zh) 一种基于滑模观测器的单相电网电压参数的鲁棒估计方法
CN107167658A (zh) 一种高精度抗干扰的电力系统基波频率及相位测量方法
CN104597320B (zh) 一种适用于多个频率交流信号计算的方法
CN103616613B (zh) 一种利用输电线路两端行波自然频率的故障定位方法
CN109828163A (zh) 一种用于电网的三相不平衡检测方法
CN106841778A (zh) 基于pmu实现的次同步和超同步谐波参数的处理方法
CN113675850B (zh) 一种基于非线性鲁棒估计的电网信息快速精准感知方法
CN105911341A (zh) 一种谐波无功功率的测量方法
Wang et al. Modeling and simulation of full-bridge series resonant converter based on generalized state space averaging
CN104502707A (zh) 一种基于三次样条插值的电力系统同步相量测量方法
CN104833852A (zh) 一种基于人工神经网络的电力系统谐波信号估计测量方法
CN106154037A (zh) 一种基于校验的同步相量自适应计算方法
CN104849552A (zh) 一种基于自适应陷波滤波器的谐波提取方法
CN111884218A (zh) 一种双馈入vsc输电系统稳定性评估方法及系统
CN102854422B (zh) 一种变压器支路三相不对称故障分析方法
CN103995180A (zh) 一种考虑不等式约束的电力系统频率估计方法
CN105203877B (zh) 能够消除量测误差影响的单回输电线路零序参数辨识方法
CN110320400B (zh) 准同步采样和改进能量算子的电压闪变包络参数提取方法
CN104931777A (zh) 一种基于两条dft复数谱线的信号频率测量方法

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