CN110954746A - 一种基于四项Nuttall余弦窗的六插值FFT算法 - Google Patents
一种基于四项Nuttall余弦窗的六插值FFT算法 Download PDFInfo
- Publication number
- CN110954746A CN110954746A CN201911185279.4A CN201911185279A CN110954746A CN 110954746 A CN110954746 A CN 110954746A CN 201911185279 A CN201911185279 A CN 201911185279A CN 110954746 A CN110954746 A CN 110954746A
- Authority
- CN
- China
- Prior art keywords
- frequency
- nuttall
- harmonic
- signal
- cosine
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 230000003595 spectral effect Effects 0.000 claims abstract description 64
- 238000005070 sampling Methods 0.000 claims abstract description 31
- 238000001228 spectrum Methods 0.000 claims abstract description 23
- 238000004364 calculation method Methods 0.000 claims abstract description 11
- 238000000034 method Methods 0.000 claims abstract description 9
- 238000012937 correction Methods 0.000 claims description 12
- 238000001514 detection method Methods 0.000 abstract description 16
- 230000008901 benefit Effects 0.000 abstract description 4
- 230000000694 effects Effects 0.000 abstract description 4
- 230000004888 barrier function Effects 0.000 abstract description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 230000001360 synchronised effect Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000010183 spectrum analysis Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R23/00—Arrangements for measuring frequencies; Arrangements for analysing frequency spectra
- G01R23/16—Spectrum analysis; Fourier analysis
Landscapes
- Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- General Physics & Mathematics (AREA)
- Complex Calculations (AREA)
Abstract
本申请提供了一种基于四项Nuttall余弦窗的六插值FFT算法,包括以下步骤:将谐波信号以采样频率fc进行均匀采样,得到的离散时间信号x(n),对离散时间信号x(n)加四项Nuttall余弦窗截断数据,然后进行FFT计算得到谐波信号频谱;分析谐波信号频谱,选择峰值频点k0附近左右的各3条谱线即6条谱线,建立谱线关联式,通过谱线关联式获得各次谐波信号的幅值、频率和相角。本申请中通过四项Nuttall余弦窗的优势及峰值频点附近6条谱线包含不同程度的幅值信息,降低常规FFT算法在非同步采样时存在的频谱泄露和栅栏效应产生的误差,从而实现更高精度的谐波幅值、频率及初相角检测;有效解决了目前常见的插值未完全利用峰值频点k0附近各谱线所包含的信息,使得检测精度不够高的问题。
Description
技术领域
本申请涉及电力谐波计量检测技术领域,尤其涉及一种基于四项Nuttall余弦窗的六插值FFT算法。
背景技术
电力***中,由于大量电力电子元器件的使用,造成电网中的谐波含量越来越多,严重影响电能质量,这样的局势威胁到电力***的安全运行。如果能够精确的检测到谐波信号,并对其采取抑制措施,便能降低谐波对电力***的影响。自从离散傅里叶变换(discrete Fourier transform,DFT)提出以来,离散频谱分析实现了信号分析处理从时域到频率的转变,且其快速算法(fast Fourier transform,FFT)以它计算快速、方便的优点已成为了数字信号分析的基础,被广泛应用于电网谐波检测的领域中。
在现阶段应用于谐波检测领域中,FFT算法仍然作为主要算法,且FFT算法能在同步采样和整周期采样情况下,进行谐波的幅值、频率及初相角的精确测量。但电网是时变的***,频率并非恒值,常规FFT算法很难再做到对谐波信号的同步和整周期采样。为降低检测误差,目前主要采取加窗插值的FFT算法。
常用的窗函数为Hanning窗、Blackman窗、Blackman-Harris窗,但均难以兼顾到主瓣宽度大、第一旁瓣衰减大且旁瓣峰值渐近衰减快的综合要求。实验证明,在插值谱线数未超过6前,计算准确度呈现持续增高,常见的插值为单谱线插值、双谱线插值、三谱线插值及四谱线插值,均未完全利用峰值频点k0附近各谱线所包含的信息,使得检测精度不够高。
发明内容
本申请提供了一种基于四项Nuttall余弦窗的六插值FFT算法,与现有FFT算法相比,平衡了准确性和计算量,降低了因频率偏移造成谐波检测出错的风险,提高了谐波检测的精度,有效解决了目前常见的插值为单谱线插值、双谱线插值、三谱线插值及四谱线插值,均未完全利用峰值频点k0附近各谱线所包含的信息,使得检测精度不够高的问题。
本申请解决上述技术问题所采取的技术方案如下:
一种基于四项Nuttall余弦窗的六插值FFT算法,包括以下步骤:
将谐波信号以采样频率fc进行均匀采样,得到的离散时间信号x(n),对离散时间信号x(n)加四项Nuttall余弦窗截断数据,然后进行FFT计算得到谐波信号频谱;
其中,根据各谱线的贡献度不同,取系数Bm,m=1,2,3,4,5,6,y1,y2,y3,y4,y5,y6分别为6条谱线的幅值绝对值;
通过所述谱线关联式获得各次谐波信号的幅值、频率和相角。
可选的,所述将谐波信号以采样频率fc进行均匀采样,得到的离散时间信号x(n),包括:
四项Nuttall余弦窗满足组合窗项数不大于4,其在时域表示形式为:
将谐波信号以采样频率fc进行均匀采样,得到的离散时间信号为:
可选的,所述对离散时间信号x(n)加四项Nuttall余弦窗截断数据,然后进行FFT计算得到谐波信号频谱,包括:
对离散时间信号x(n)加四项Nuttall余弦窗,令W(λ)为w(n)的离散傅里叶变换、X(λ)为x(n)加窗函数的离散傅里叶变换,如下:
其中:Δf=fc/N,λ=1,2,3,...,N-1,ft为t次谐波;
由于负频点谱峰的影响较小,由上式可得t次谐波分析如下式:
引入变量α,取α=k1-k0-0.5,由于0<k1-k0<1,则α∈(-0.5,+0.5);
根据各谱线的贡献度不同,取系数Bm,m=1,2,3,4,5,6;
引入变量β,使得β与各谱线产生关联,得到:
求出频率修正式为:
ft=k0Δf=(α+k1+0.5)Δf
其中,k1可以利用峰值索引的方式得出;
由于采用点数N足够大,幅值修正式为:
At=N-1(B1y1+B2y2+B3y3+B4y4+B5y5+B6y6)n(α)
频率修正式为:
可选的,所述求出α的拟合曲线为n(a)包括采用反函数及拟合曲线的算法得出的n(a)。
本申请提供的技术方案包括以下有益技术效果:
本申请提供了一种基于四项Nuttall余弦窗的六插值FFT算法,包括以下步骤:将谐波信号以采样频率fc进行均匀采样,得到的离散时间信号x(n),对离散时间信号x(n)加四项Nuttall余弦窗截断数据,然后进行FFT计算得到谐波信号频谱;分析谐波信号频谱,选择峰值频点k0附近左右的各3条谱线即6条谱线,建立谱线关联式,通过谱线关联式获得各次谐波信号的幅值、频率和相角。本申请中通过四项Nuttall余弦窗的优势及峰值频点附近6条谱线包含不同程度的幅值信息,降低常规FFT算法在非同步采样时存在的频谱泄露和栅栏效应产生的误差,从而实现更高精度的谐波幅值、频率及初相角检测;有效解决了目前常见的插值为单谱线插值、双谱线插值、三谱线插值及四谱线插值,均未完全利用峰值频点k0附近各谱线所包含的信息,使得检测精度不够高的问题。
附图说明
为了更清楚地说明本申请的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,显而易见地,对于本领域普通技术人员而言,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为本申请实施例提供的非同步采样6谱线视图。
具体实施方式
为了使本领域技术人员更好地理解本申请中的技术方案,下面将结合本申请实施例中的附图,对申请实施例中的技术方案进行清楚、完整地描述;显然,所描述的实施例仅是本申请的一部分实施例,而不是全部的实施例。基于本申请中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都应当属于本申请保护的范围。
本申请实施例提供的一种基于四项Nuttall余弦窗的六插值FFT算法,包括以下步骤:
将谐波信号以采样频率fc进行均匀采样,得到的离散时间信号x(n),对离散时间信号x(n)加四项Nuttall余弦窗截断数据,然后进行FFT计算得到谐波信号频谱;
分析所述谐波信号频谱,在非同步采样时,峰值频点k0为非整数,使得其不会正好处于离散谱线频点上,由于以k0为中心,其附近左右的各3条谱线均含有不同程度的幅值信息,故插值谱线应考虑这6条谱线,为保证完全利用目标谱线附件的谱线信息,因此选择峰值频点k0附近左右的各3条谱线即6条谱线,即如图1中所示k1,k2,k3,k4,k5,k6,建立谱线关联式,
其中,根据6条谱线对求取修正值的贡献度不同,取系数Bm,m=1,2,3,4,5,6,y1,y2,y3,y4,y5,y6分别为6条谱线的幅值绝对值;目标谱线附近的6条谱线因距离目标谱线距离不同,其在算法中对目标谱线的影响不一,因此对取修正值的贡献度不同。
通过所述谱线关联式获得各次谐波信号的幅值、频率和相角。
本申请实施例提供的基于四项Nuttall余弦窗的六插值FFT算法,针对电网频率偏移时,在FFT算法的基础上利用四项Nuttall余弦窗,降低常规FFT算法因非同步采样时存在的频谱泄露而产生的误差;接着利用离散谱线频点上目标谱线附件的6条谱线,降低常规FFT算法因非同步采样时存在的栅栏效应而产生的误差。本发明提供的基于四项余弦窗的六插值FFT算法,与现有FFT算法相比,平衡了准确性和计算量,降低了因频率偏移造成谐波检测出错的风险,提高了谐波检测的精度。
Hanning窗、Blackman窗、Blackman-Harris窗的主瓣宽度分别为8π/N、12π/N、16π/N,Hanning窗、Blackman窗、Blackman-Harris窗的旁瓣衰减幅值分别为-32dB、-52dB、-92dB,Hanning窗和Blackman窗的衰减速率为18dB/oct,Blackman-Harris窗为6dB/oct。Blackman-Harris窗虽然有很好的主瓣宽度及旁瓣衰减幅值,但其衰减速率仍有不足。
单谱线插值算法简单、运算快,但校正精度低;双谱线、三谱线、四谱线虽利用了峰谱频点附近的多个谱线信息,但仍未完全考虑更多的谱线信息。
为满足一般电网中整数次谐波的检测,结合窗函数的特点,窗函数应选择主瓣宽度大、第一旁瓣衰减大且旁瓣峰值渐近衰减快的四项Nuttall余弦窗。结合余弦信号在经过加窗后出现泄漏现象,即在非同步采样时,峰值频点k0为非整数,使得其不会正好处于离散谱线频点上,由于以k0为中心,其附近左右的各3条谱线均含有不同程度的幅值信息,故插值谱线应考虑这6条谱线。
进一步的,所述将谐波信号以采样频率fc进行均匀采样,得到的离散时间信号x(n),包括:
四项Nuttall余弦窗满足组合窗项数不大于4,其在时域表示形式为:
将谐波信号以采样频率fc进行均匀采样,得到的离散时间信号为:
可选的,所述对离散时间信号x(n)加四项Nuttall余弦窗截断数据,然后进行FFT计算得到谐波信号频谱,包括:
对离散时间信号x(n)加四项Nuttall余弦窗,令W(λ)为w(n)的离散傅里叶变换、X(λ)为x(n)加窗函数的离散傅里叶变换,如下:
其中:Δf=fc/N,λ=1,2,3,...,N-1,ft为t次谐波;
由于负频点谱峰的影响较小,由上式可得t次谐波分析如下式:
引入变量α,取α=k1-k0-0.5,由于0<k1-k0<1,则α∈(-0.5,+0.5);
根据各谱线的贡献度不同,取系数Bm,m=1,2,3,4,5,6;
引入变量β,使得β与各谱线产生关联,得到:
求出频率修正式为:
ft=k0Δf=(α+k1+0.5)Δf
其中,k1可以利用峰值索引的方式得出;
由于采用点数N足够大,幅值修正式为:
At=N-1(B1y1+B2y2+B3y3+B4y4+B5y5+B6y6)n(α)
频率修正式为:
可选的,所述求出α的拟合曲线为n(a)包括采用反函数及拟合曲线的算法得出的n(a)。
本申请实施例旨在通过四项Nuttall余弦窗的优势及峰值频点附近6条谱线包含不同程度的幅值信息,实现更高精度的谐波幅值、频率及初相角检测。
需要说明的是,术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的物品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下,由语句“包括一个……”限定的要素,并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。
以上所述仅是本申请的具体实施方式,使本领域技术人员能够理解或实现本申请。对这些实施例的多种修改对本领域的技术人员来说将是显而易见的,本文中所定义的一般原理可以在不脱离本申请的精神或范围的情况下,在其它实施例中实现。因此,本申请将不会被限制于本文所示的这些实施例,而是要符合与本文所公开的原理和新颖特点相一致的最宽的范围。
应当理解的是,本申请并不局限于上面已经描述并在附图中示出的内容,并且可以在不脱离其范围进行各种修改和改变。本申请的范围仅由所附的权利要求来限制。
Claims (5)
4.根据权利要求3所述的基于四项Nuttall余弦窗的六插值FFT算法,其特征在于,所述分析所述谐波信号频谱,选择峰值频点k0附近左右的各3条谱线即6条谱线,即k1,k2,k3,k4,k5,k6,建立谱线关联式,包括:
引入变量α,取α=k1-k0-0.5,由于0<k1-k0<1,则α∈(-0.5,+0.5);
根据各谱线的贡献度不同,取系数Bm,m=1,2,3,4,5,6;
引入变量β,使得β与各谱线产生关联,得到:
求出频率修正式为:
ft=k0Δf=(α+k1+0.5)Δf
其中,k1可以利用峰值索引的方式得出;
由于采用点数N足够大,幅值修正式为:
At=N-1(B1y1+B2y2+B3y3+B4y4+B5y5+B6y6)n(α)
频率修正式为:
5.根据权利要求4所述的基于四项Nuttall余弦窗的六插值FFT算法,其特征在于,所述求出α的拟合曲线为n(a)包括采用反函数及拟合曲线的算法得出的n(a)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911185279.4A CN110954746A (zh) | 2019-11-27 | 2019-11-27 | 一种基于四项Nuttall余弦窗的六插值FFT算法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911185279.4A CN110954746A (zh) | 2019-11-27 | 2019-11-27 | 一种基于四项Nuttall余弦窗的六插值FFT算法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN110954746A true CN110954746A (zh) | 2020-04-03 |
Family
ID=69978685
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911185279.4A Pending CN110954746A (zh) | 2019-11-27 | 2019-11-27 | 一种基于四项Nuttall余弦窗的六插值FFT算法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110954746A (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111929532A (zh) * | 2020-08-06 | 2020-11-13 | 国网北京市电力公司 | 一种10kV配电电缆中间接头定位方法 |
CN114185308A (zh) * | 2021-11-29 | 2022-03-15 | 华中科技大学 | 基于频谱稀疏和校正的主轴误差测量方法及*** |
CN115356534A (zh) * | 2022-10-21 | 2022-11-18 | 成都特隆美储能技术有限公司 | 一种减少电网电流fft谐波提取时频谱泄露的方法 |
CN117949964A (zh) * | 2024-03-27 | 2024-04-30 | 青岛科技大学 | 一种激光多普勒测速***信号处理方法 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103197141A (zh) * | 2013-03-21 | 2013-07-10 | 广东电网公司电力科学研究院 | 一种测量电力***信号频率及谐波参数的方法 |
CN104897960A (zh) * | 2015-06-15 | 2015-09-09 | 中南民族大学 | 基于加窗四谱线插值fft的谐波快速分析方法及*** |
CN105137180A (zh) * | 2015-08-24 | 2015-12-09 | 三峡大学 | 基于六项余弦窗四谱线插值的高精度谐波分析方法 |
CN105353215A (zh) * | 2015-12-23 | 2016-02-24 | 合肥工业大学 | 基于Nuttall窗四谱线插值FFT的谐波检测方法 |
CN106814230A (zh) * | 2016-12-15 | 2017-06-09 | 贵州电网有限责任公司贵阳供电局 | 基于Kaiser窗FFT四峰插值的电能计量方法 |
CN107179439A (zh) * | 2017-06-26 | 2017-09-19 | 贵州电网有限责任公司贵阳供电局 | 基于五项余弦窗fft三峰插值的电能计量方法 |
CN107643446A (zh) * | 2017-08-11 | 2018-01-30 | 中南民族大学 | 一种基于主瓣宽度的多谱线插值谐波分析方法及*** |
CN109541304A (zh) * | 2018-11-21 | 2019-03-29 | 三峡大学 | 基于六项最小旁瓣窗插值的电网高次弱幅值谐波检测方法 |
CN109782063A (zh) * | 2018-10-23 | 2019-05-21 | 国网安徽省电力有限公司芜湖供电公司 | 一种基于纳托尔自卷积窗三谱线插值fft的动态间谐波分析方法 |
-
2019
- 2019-11-27 CN CN201911185279.4A patent/CN110954746A/zh active Pending
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103197141A (zh) * | 2013-03-21 | 2013-07-10 | 广东电网公司电力科学研究院 | 一种测量电力***信号频率及谐波参数的方法 |
CN104897960A (zh) * | 2015-06-15 | 2015-09-09 | 中南民族大学 | 基于加窗四谱线插值fft的谐波快速分析方法及*** |
CN105137180A (zh) * | 2015-08-24 | 2015-12-09 | 三峡大学 | 基于六项余弦窗四谱线插值的高精度谐波分析方法 |
CN105353215A (zh) * | 2015-12-23 | 2016-02-24 | 合肥工业大学 | 基于Nuttall窗四谱线插值FFT的谐波检测方法 |
CN106814230A (zh) * | 2016-12-15 | 2017-06-09 | 贵州电网有限责任公司贵阳供电局 | 基于Kaiser窗FFT四峰插值的电能计量方法 |
CN107179439A (zh) * | 2017-06-26 | 2017-09-19 | 贵州电网有限责任公司贵阳供电局 | 基于五项余弦窗fft三峰插值的电能计量方法 |
CN107643446A (zh) * | 2017-08-11 | 2018-01-30 | 中南民族大学 | 一种基于主瓣宽度的多谱线插值谐波分析方法及*** |
CN109782063A (zh) * | 2018-10-23 | 2019-05-21 | 国网安徽省电力有限公司芜湖供电公司 | 一种基于纳托尔自卷积窗三谱线插值fft的动态间谐波分析方法 |
CN109541304A (zh) * | 2018-11-21 | 2019-03-29 | 三峡大学 | 基于六项最小旁瓣窗插值的电网高次弱幅值谐波检测方法 |
Non-Patent Citations (3)
Title |
---|
华敏 等: "一种高精度六谱线插值FFT谐波与间谐波分析方法", 《电力***保护与控制》 * |
李晓辉 等: "基于四项余弦窗四谱线插值FFT的电力谐波算法分析", 《沈阳农业大学学报》 * |
王保帅 等: "一种用于谐波分析的高精度多谱线插值算法", 《电工技术学报》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111929532A (zh) * | 2020-08-06 | 2020-11-13 | 国网北京市电力公司 | 一种10kV配电电缆中间接头定位方法 |
CN114185308A (zh) * | 2021-11-29 | 2022-03-15 | 华中科技大学 | 基于频谱稀疏和校正的主轴误差测量方法及*** |
CN115356534A (zh) * | 2022-10-21 | 2022-11-18 | 成都特隆美储能技术有限公司 | 一种减少电网电流fft谐波提取时频谱泄露的方法 |
CN117949964A (zh) * | 2024-03-27 | 2024-04-30 | 青岛科技大学 | 一种激光多普勒测速***信号处理方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110954746A (zh) | 一种基于四项Nuttall余弦窗的六插值FFT算法 | |
CN104897960B (zh) | 基于加窗四谱线插值fft的谐波快速分析方法及*** | |
CN102818930B (zh) | 一种高精度快速计算电力谐波参数的方法 | |
CN110837001B (zh) | 一种电力***中谐波和间谐波的分析方法与装置 | |
CN102539915B (zh) | 时延傅立叶变换测频法精确计算电力谐波参数方法 | |
CN105137180B (zh) | 基于六项余弦窗四谱线插值的高精度谐波分析方法 | |
CN101113995A (zh) | 基于Nuttall窗双峰插值FFT的基波与谐波检测方法 | |
CN107643446B (zh) | 一种基于主瓣宽度的多谱线插值谐波分析方法及*** | |
CN103308766A (zh) | 一种基于凯撒自卷积窗双谱线插值fft谐波分析方法及其装置 | |
CN101261292A (zh) | 基于5项Rife-Vincent(I)窗双谱线插值FFT的基波与谐波检测方法 | |
CN110837003B (zh) | 一种基于三角窗的双窗全相位dft同步相量测量方法及*** | |
CN105307095B (zh) | 一种基于fft的高分辨率音频频率测量方法 | |
Luo et al. | Frequency estimation of the weighted real tones or resolved multiple tones by iterative interpolation DFT algorithm | |
Huibin et al. | Energy based signal parameter estimation method and a comparative study of different frequency estimators | |
Borkowski et al. | Comparison of sine-wave frequency estimation methods in respect of speed and accuracy for a few observed cycles distorted by noise and harmonics | |
CN109655665A (zh) | 基于布莱克曼窗的全相位傅里叶谐波分析方法 | |
CN102495285B (zh) | 对称窗函数功率重心估计电力谐波参数的方法 | |
Li et al. | Frequency estimation based on modulation FFT and MUSIC algorithm | |
CN104931777B (zh) | 一种基于两条dft复数谱线的信号频率测量方法 | |
CN105372492B (zh) | 基于三条dft复数谱线的信号频率测量方法 | |
CN113129912B (zh) | 一种单音信号的检测方法 | |
CN111579868B (zh) | 一种高次谐波的测量方法及装置 | |
CN114487589A (zh) | 电网宽频信号自适应测量方法、装置及*** | |
Lewandowski et al. | Optimal base frequency estimation of an electrical signal based on Prony’s estimator and a FIR filter | |
Yue et al. | Modified algorithm of sinusoid signal frequency estimation based on Quinn and Aboutanios iterative algorithms |
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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20200403 |
|
RJ01 | Rejection of invention patent application after publication |