CN113553997A - 一种联合改进小波阈值的信号处理方法、存储介质和*** - Google Patents
一种联合改进小波阈值的信号处理方法、存储介质和*** Download PDFInfo
- Publication number
- CN113553997A CN113553997A CN202111092643.XA CN202111092643A CN113553997A CN 113553997 A CN113553997 A CN 113553997A CN 202111092643 A CN202111092643 A CN 202111092643A CN 113553997 A CN113553997 A CN 113553997A
- Authority
- CN
- China
- Prior art keywords
- wavelet
- signal
- threshold
- noise
- processing
- 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
- 238000003672 processing method Methods 0.000 title claims abstract description 9
- 238000000034 method Methods 0.000 claims abstract description 68
- 238000012545 processing Methods 0.000 claims abstract description 22
- 238000009825 accumulation Methods 0.000 claims abstract description 18
- 230000006872 improvement Effects 0.000 claims abstract description 8
- 238000000354 decomposition reaction Methods 0.000 claims description 21
- 230000000694 effects Effects 0.000 claims description 15
- 238000013139 quantization Methods 0.000 claims description 11
- 230000009466 transformation Effects 0.000 claims description 5
- 238000005516 engineering process Methods 0.000 abstract description 4
- 239000013307 optical fiber Substances 0.000 abstract description 4
- 238000012935 Averaging Methods 0.000 description 6
- 238000004364 calculation method Methods 0.000 description 6
- 238000010586 diagram Methods 0.000 description 6
- 238000005259 measurement Methods 0.000 description 5
- 230000001186 cumulative effect Effects 0.000 description 4
- 238000004043 dyeing Methods 0.000 description 4
- 238000010835 comparative analysis Methods 0.000 description 3
- 230000007547 defect Effects 0.000 description 3
- 230000009467 reduction Effects 0.000 description 3
- 238000004088 simulation Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 238000012067 mathematical method Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000010355 oscillation Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000002269 spontaneous effect Effects 0.000 description 1
- 238000010186 staining Methods 0.000 description 1
Images
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/02—Preprocessing
- G06F2218/04—Denoising
- G06F2218/06—Denoising by applying a scale-space analysis, e.g. using wavelet analysis
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01D—MEASURING NOT SPECIALLY ADAPTED FOR A SPECIFIC VARIABLE; ARRANGEMENTS FOR MEASURING TWO OR MORE VARIABLES NOT COVERED IN A SINGLE OTHER SUBCLASS; TARIFF METERING APPARATUS; MEASURING OR TESTING NOT OTHERWISE PROVIDED FOR
- G01D5/00—Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable
- G01D5/26—Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable characterised by optical transfer means, i.e. using infrared, visible, or ultraviolet light
- G01D5/32—Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable characterised by optical transfer means, i.e. using infrared, visible, or ultraviolet light with attenuation or whole or partial obturation of beams of light
- G01D5/34—Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable characterised by optical transfer means, i.e. using infrared, visible, or ultraviolet light with attenuation or whole or partial obturation of beams of light the beams of light being detected by photocells
- G01D5/353—Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable characterised by optical transfer means, i.e. using infrared, visible, or ultraviolet light with attenuation or whole or partial obturation of beams of light the beams of light being detected by photocells influencing the transmission properties of an optical fibre
- G01D5/35338—Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable characterised by optical transfer means, i.e. using infrared, visible, or ultraviolet light with attenuation or whole or partial obturation of beams of light the beams of light being detected by photocells influencing the transmission properties of an optical fibre using other arrangements than interferometer arrangements
- G01D5/35354—Sensor working in reflection
- G01D5/35358—Sensor working in reflection using backscattering to detect the measured quantity
- G01D5/35364—Sensor working in reflection using backscattering to detect the measured quantity using inelastic backscattering to detect the measured quantity, e.g. using Brillouin or Raman backscattering
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Signal Processing (AREA)
- General Engineering & Computer Science (AREA)
- Image Processing (AREA)
Abstract
本发明公开了一种联合改进小波阈值的信号处理方法、存储介质和***,涉及光电信息技术中光纤传感的信号处理领域,具体涉及一种提高信号信噪比方法,在D.L.Donoho提出的小波阈值去噪方法的基础上,对其阈值函数进行改进,提出了一种新的阈值函数,该阈值函数不仅能够有效的克服硬阈值函数不连续的缺点,还可以解决软阈值函数存在的恒定偏差问题,能有效提高去噪效果。并在此方法上,将累加平均法联合改进小波阈值的方法来进行信号去噪处理,能够进一步提升输出信号的信噪比,进一步减小了均方差。
Description
技术领域
本发明属于光电信息技术中光纤传感的信号处理领域,具体涉及一种联合改进小波阈值的信号处理方法、存储介质和***。
背景技术
布里渊光时域反射仪(BOTDR)是利用布里渊频移与应力/应变以及温度的线性关系来实现光纤传感的技术,可以应用于许多大型的建设,如桥梁大坝、隧道交通、传输电网等的安全监测,具有广大的应用前景。但是由于自发的布里渊散射信号十分微弱,并且有着各种噪声的干扰,在长距离测量时,传感光纤中的有用信号就很容易淹没在这些噪声之中,使得仪器难以准确探测,所以提高BOTDR***采集信号的信噪比显得十分重要。
传统去噪是建立在傅里叶变换基础上的线性滤波方法,这种方法虽然简单并且易于实现,但是在提高信噪比和分辨率上存在着矛盾。而小波变换能将时频特征同时局部化,是一种多分辨率的分析方法,被广泛应用于信号去噪。在1995年,D.L.Donoho等人在小波变换的基础上提出了一种小波阈值去噪算法,引起了广泛的关注。小波阈值去噪是在小波变换的基础上对含噪信号进行分解和重构,经小波分解后,有用信号的系数要大于噪声的系数,于是通过设定合适的阈值限定分解后的小波系数,认为绝对值比该阈值小的小波系数主要由噪声引起,直接置为零,把绝对值比该阈值大的系数认为是由信号引起的,通过某个阈值函数处理将其收缩(软阈值方法)或者保留(硬阈值方法)得到预估的小波系数,最后通过对预估的小波系数重构得到去噪信号。
小波阈值去噪的去噪效果取决于阈值的选取以及阈值函数的设计,常用的阈值去噪方法有硬阈值和软阈值两种去噪方法。采用硬阈值方法能更多的保留真实信号的尖峰等特点,但因其在阈值处不连续,去噪效果会出现震荡;软阈值方法克服了硬阈值法在阈值处不连续的缺点,能在去噪后产生光滑的结果,但原始小波系数与阈值量化处理后的小波系数之间总存在恒定偏差,容易丢失信号的高频信息,一定程度上影响了去噪效果。
累加平均法是一种通过多次累加测量数据来提升信噪比的方法,曾被广泛应用于信号降噪处理中;但当累加次数达到一定数值后,信噪比提升速度明显变慢,而所需存储空间和测量时间依然线性增加,这不利于硬件实现,也无法满足测量实时性要求。
发明内容
本发明要解决的问题是,克服常用的阈值去噪方法中存在不连续和恒定偏差的缺陷,提出一种联合改进小波阈值函数的信号处理方法,该方法能够抑制震荡,消除恒定偏差,减少信号丢失,有效提高信号的信噪比。
本发明解决其问题所采用的技术方案是:一种联合改进小波阈值的信号处理方法,包含以下步骤:
S1 模拟BOTDR***数据采集模块采集到的布里渊散射信号,模型如下:
其中x(t)为含噪信号,s(t)为原始信号,n(t)为噪声信号。则其信噪比即可表示为:
进一步地,模拟BOTDR***数据采集模块采集到的布里渊散射信号,包括:
构造正弦函数,并采用randn函数对其添加白噪声,得到一个一维染噪信号;
选取sym6、sym7、sym8和db1、db3、db4小波基对小波进行小波分解,分解层数设定为3~5层,通过matlab仿真对比分析,采用db4作为小波基,对所述一维染噪信号进行5层小波分解,得到高频小波系数;
采用改进的阈值计算方法计算出阈值,其中改进的阈值计算方法表达式为:
其中μ为改进阈值函数的调整因子,且μ>0;
进一步地,所述方法还包括引入输出信噪比(SNR)、均方误差(RMSE)来验证函数的去噪效果,其计算公式分别为:
当信号去噪后,信号信噪比越高、均方误差越小,表明去噪后的信号越接近原始信号,去噪效果越好。
一种存储介质,所述存储介质记录有上述一种联合改进小波阈值的信号处理方法。
一种计算执行***所述计算执行***用于计算有上述一种联合改进小波阈值的信号处理方法。
采用本申请的技术方法的有益效果如下:
本发明的一种联合改进小波阈值的信号处理方法,新提出的小波阈值函数不仅克服了传统小波阈值函数不连续的缺点,还有效解决了传统小波阈值函数存在的恒定偏差问题,同时采用累加平均法联合本文改进小波阈值方法对染噪信号进行信号处理,能够有效提高信噪比,减小均方误差,使得处理后的信号尽可能接近原始信号。
附图说明
图1为本发明信号处理流程图;
图2为染噪信号信噪比-5.88db的原始信号以及染噪信号;
图3为染噪信号信噪比-11.98db的原始信号以及染噪信号;
图4为改进阈值函数的示意图;
图5为染噪信号信噪比-5.88db的去噪效果图;
图6为染噪信号信噪比-5.88db的累加联合去噪效果图;
图7为染噪信号信噪比-11.98db去噪效果图;
图8为染噪信号信噪比-11.98db累加联合去噪效果图。
具体实施方式
下面将结合本公开实施例中的附图,对本公开实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本公开一部分实施例,而不是全部的实施例。基于本公开中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其它实施例,都属于本公开保护的范围。
如图1至图8所示,一种联合改进小波阈值的信号处理方法,包括以下几个步骤:
步骤S1:模拟原始信号;
在本实施例中,步骤S1中所述的模拟的原始信号表达式为:
步骤S2:在原始信号中添加白噪声,得到一维染噪信号;
步骤S3:将得到的一维染噪信号用累加平均法进行去噪;
累加次数一般设定为2的指数倍,但当累加次数达到一定数值后,信噪比提升速度明显变慢,所需存储空间和测量时间却依然线性增加,在本实施例中,对一维染噪信号进行512次的累加处理;
在本实施例中,选取sym6、sym7、sym8和db1、db3、db4小波基对小波进行小波分解,通过matlab仿真对比分析,选取db4作为小波基,并且在选取db4为小波基后,还对比分析了3~5层的分解层数,最终选取了5层的分解层数,即,选取db4小波基对所述一维染噪信号进行5层小波分解,得到高频小波系数;
在本实施例中,步骤S5具体包括以下步骤:
步骤S51:采用改进阈值规则计算出阈值,其中改进的阈值表达式为:
其中μ为改进阈值函数的调整因子,且μ>0;这个调整因子的实际值可以根据使用情况自适应的调节,拥有更强的自适应性;
下面利用数学方法证明阈值函数的连续性和改进现有固定偏差问题;
综上:本方案构造的改进阈值函数在±th处均连续,改进阈值函数的连续性表明改进的阈值函数能够避免吉布斯震荡的产生;
偏差性:
进一步地,本发明上述方法中还包括引入输出信噪比(SNR)、均方误差(RMSE)来验证去噪效果,其计算公式分别为:
利用模拟的染噪信号在不同信噪比的情况下,采用本发明的去噪方法与传统阈值去噪方法进行信号去噪,得到本发明与传统阈值去噪方法的输出信噪比和均方误差;
表1:染噪信号信噪比-5.88db下不同方法去噪效果对比结果
去噪方法 | 信噪比(SNR) | 均方误差(RMSE) |
一维染噪信号 | -5.8754 | 5.0120 |
软阈值方法 | 5.4018 | 1.3682 |
硬阈值方法 | 4.6773 | 1.4872 |
改进阈值方法 | 5.7091 | 1.3206 |
累加平均法 | 18.0596 | 0.3186 |
累加联合软阈值方法 | 19.3931 | 0.2733 |
累加联合硬阈值方法 | 23.3763 | 0.1728 |
累加联合改进阈值方法 | 24.7935 | 0.1467 |
表2:染噪信号信噪比-11.98db下不同方法去噪效果对比结果
去噪方法 | 信噪比(SNR) | 均方误差(RMSE) |
一维染噪信号 | -11.9763 | 10.1170 |
软阈值方法 | 1.9412 | 2.0379 |
硬阈值方法 | 1.9412 | 2.0379 |
改进阈值方法 | 2.0334 | 2.0164 |
累加平均法 | 12.4703 | 0.6064 |
累加联合软阈值方法 | 15.7114 | 0.4175 |
累加联合硬阈值方法 | 19.0580 | 0.2840 |
累加联合改进阈值方法 | 20.5497 | 0.2392 |
在本实施例中,如图5至图8所示为不同阈值函数去噪效果图,结合表1、表2不同方法去噪效果对比分析结果可以看出,虽然传统的小波阈值方法也能提高信号的信噪比,减小均方误差,但是本方案提出的联合改进小波阈值的方法能够进一步提高信噪比,减小均方误差,最大程度的有效还原信号。
工作原理
通过采用累加平均法联合本文改进小波阈值的方法对染噪信号进行信号处理,能够有效提高信噪比,减小均方误差,使得处理后的信号尽可能接近原始信号。
原文以上的实施例仅为说明发明的技术思想,不能以此限定本发明的保护范围,凡是按照本发明提出的技术思想,在技术方案基础上所做的任何改动,均落入本发明保护范围之内。本发明未涉及的技术均可通过现有的技术加以实现。
Claims (10)
1.一种联合改进小波阈值的信号处理方法,其特征在于,所述方法包括以下步骤:对累加处理的染噪信号进行多层小波分解,得到高频小波系数;对高频小波系数进行阈值量化处理,得到低频小波系数;将阈值量化处理之后的每层的高频小波系数和低频小波系数实行小波逆变换,进行信号重构,得到重构信号。
2.根据权利要求1所述的一种联合改进小波阈值的信号处理方法,其特征在于,所述染噪信号通过构造正弦函数,并对其添加白噪声得到。
3.根据权利要求1所述的一种联合改进小波阈值的信号处理方法,其特征在于,所述染噪信号进行多层小波分解采用sym6、sym7、sym8、db1、db3、db4小波基中的一种进行小波分解。
4.根据权利要求1所述的一种联合改进小波阈值的信号处理方法,其特征在于,所述小波分解的分解层数设定为3~5层。
5.根据权利要求1所述的一种联合改进小波阈值的信号处理方法,其特征在于,所述低频小波系的获取方法,包括以下步骤:
第一步:采用改进阈值规则计算出阈值;
第二步:利用该阈值与改进后的小波阈值函数对高频小波系数进行阈值量化处理,得到低频小波系数。
6.根据权利要求1所述的一种联合改进小波阈值的信号处理方法,其特征在于,所述染噪信号的累加处理次数为2的指数倍。
7.根据权利要求1所述的一种联合改进小波阈值的信号处理方法,其特征在于,所述一种联合改进小波阈值的信号处理方法包括引入输出信噪比、均方误差来验证去噪效果。
8.根据权利要求1所述的一种联合改进小波阈值的信号处理的方法,其特征在于,所述高频小波系数通过db4小波基对所述染噪信号进行5层小波分解得到。
9.一种存储介质,其特征在于,所述存储介质记录有权利要求1所述的一种联合改进小波阈值的信号处理方法。
10.一种计算执行***,其特征在于,所述计算执行***用于执行权利要求1所述的一种联合改进小波阈值的信号处理方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111092643.XA CN113553997A (zh) | 2021-09-17 | 2021-09-17 | 一种联合改进小波阈值的信号处理方法、存储介质和*** |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111092643.XA CN113553997A (zh) | 2021-09-17 | 2021-09-17 | 一种联合改进小波阈值的信号处理方法、存储介质和*** |
Publications (1)
Publication Number | Publication Date |
---|---|
CN113553997A true CN113553997A (zh) | 2021-10-26 |
Family
ID=78134645
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111092643.XA Pending CN113553997A (zh) | 2021-09-17 | 2021-09-17 | 一种联合改进小波阈值的信号处理方法、存储介质和*** |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113553997A (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114577419A (zh) * | 2022-04-24 | 2022-06-03 | 南京信息工程大学 | 一种提高地下连续墙渗漏安全监测有效性方法 |
CN116108336A (zh) * | 2023-04-13 | 2023-05-12 | 吉林省百皓科技有限公司 | 一种基于小波变换的二氧化氯传感器信号去噪方法 |
CN116399379A (zh) * | 2023-06-07 | 2023-07-07 | 山东省科学院激光研究所 | 分布式光纤声波传感***及其测量方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108985179A (zh) * | 2018-06-22 | 2018-12-11 | 福建和盛高科技产业有限公司 | 一种基于改进小波阈值函数的电能质量信号去噪方法 |
CN112395992A (zh) * | 2020-11-18 | 2021-02-23 | 云南电网有限责任公司电力科学研究院 | 一种基于改进小波阈值的电力谐波信号去噪方法 |
AU2021101814A4 (en) * | 2021-04-08 | 2021-06-17 | Achanta, Sampath Dakshina Murthy MR | A novel image denoising method with hybrid dual tree complex wavelet transform |
CN113221746A (zh) * | 2021-05-13 | 2021-08-06 | 西南科技大学 | 基于改进小波阈值函数的微震信号去噪方法 |
-
2021
- 2021-09-17 CN CN202111092643.XA patent/CN113553997A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108985179A (zh) * | 2018-06-22 | 2018-12-11 | 福建和盛高科技产业有限公司 | 一种基于改进小波阈值函数的电能质量信号去噪方法 |
CN112395992A (zh) * | 2020-11-18 | 2021-02-23 | 云南电网有限责任公司电力科学研究院 | 一种基于改进小波阈值的电力谐波信号去噪方法 |
AU2021101814A4 (en) * | 2021-04-08 | 2021-06-17 | Achanta, Sampath Dakshina Murthy MR | A novel image denoising method with hybrid dual tree complex wavelet transform |
CN113221746A (zh) * | 2021-05-13 | 2021-08-06 | 西南科技大学 | 基于改进小波阈值函数的微震信号去噪方法 |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114577419A (zh) * | 2022-04-24 | 2022-06-03 | 南京信息工程大学 | 一种提高地下连续墙渗漏安全监测有效性方法 |
CN116108336A (zh) * | 2023-04-13 | 2023-05-12 | 吉林省百皓科技有限公司 | 一种基于小波变换的二氧化氯传感器信号去噪方法 |
CN116399379A (zh) * | 2023-06-07 | 2023-07-07 | 山东省科学院激光研究所 | 分布式光纤声波传感***及其测量方法 |
CN116399379B (zh) * | 2023-06-07 | 2023-11-03 | 山东省科学院激光研究所 | 分布式光纤声波传感***及其测量方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113553997A (zh) | 一种联合改进小波阈值的信号处理方法、存储介质和*** | |
CN112395992B (zh) | 一种基于改进小波阈值的电力谐波信号去噪方法 | |
CN109557429B (zh) | 基于改进小波阈值去噪的gis局部放电故障检测方法 | |
CN103630808B (zh) | 一种基于提升小波变换的局部放电信号去噪方法 | |
CN113378661A (zh) | 一种基于改进小波阈值和相关检测的直流电能信号去噪方法 | |
CN112084845B (zh) | 基于多尺度小波系数自相关的低频1/f噪声消除方法 | |
CN113723171B (zh) | 基于残差生成对抗网络的脑电信号去噪方法 | |
CN109709585B (zh) | 去除gps坐标时间序列中有色噪声的方法 | |
CN113297987B (zh) | 一种基于双目标函数优化的变分模态分解信号降噪方法 | |
CN107886078A (zh) | 一种基于分层自适应阈值函数的小波阈值降噪方法 | |
CN106569034A (zh) | 一种基于小波和高阶pde的局部放电信号去噪方法 | |
CN114577419B (zh) | 一种提高地下连续墙渗漏安全监测有效性方法 | |
CN115165274A (zh) | 工程机械设备振动状态自适应智能监测装置及方法 | |
CN105550998A (zh) | 基于二代小波整数变换的图像增强方法及图像增强*** | |
CN109724693B (zh) | 一种基于平稳小波的融合光谱去噪方法 | |
CN113553960B (zh) | 一种基于小波包方差熵的风电爬坡不确定性评估方法 | |
CN114330445A (zh) | 一种基于变压器振动信号敏感imf的小波阈值去噪方法 | |
CN114690003A (zh) | 一种基于eemd的局放信号降噪方法 | |
CN110333054A (zh) | 一种针对白车身焊接设备的缓变微小故障检测方法 | |
CN113702666A (zh) | 一种光纤陀螺惯性测量单元信号联合降噪方法 | |
CN109558857B (zh) | 一种混沌信号降噪方法 | |
CN110287853B (zh) | 一种基于小波分解的暂态信号去噪方法 | |
CN116383605A (zh) | 基于小波分解的车辆震动信号去噪方法 | |
CN113837042B (zh) | 一种基于局部均值分解的小波阈值自动化变形监测信号去噪方法 | |
CN110703089B (zh) | 一种用于低频振荡Prony分析的小波阈值去噪方法 |
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: 20211026 |
|
RJ01 | Rejection of invention patent application after publication |