CN110032707B - 一种基于fpga的四进小波多通道信号处理方法 - Google Patents
一种基于fpga的四进小波多通道信号处理方法 Download PDFInfo
- Publication number
- CN110032707B CN110032707B CN201910285083.6A CN201910285083A CN110032707B CN 110032707 B CN110032707 B CN 110032707B CN 201910285083 A CN201910285083 A CN 201910285083A CN 110032707 B CN110032707 B CN 110032707B
- Authority
- CN
- China
- Prior art keywords
- wavelet
- way
- coefficient
- way wavelet
- coefficients
- 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.)
- Active
Links
Images
Classifications
-
- 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
- G06F17/148—Wavelet transforms
-
- 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/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Computational Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Theoretical Computer Science (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Computing Systems (AREA)
- Complex Calculations (AREA)
Abstract
本发明公开了一种基于FPGA的四进小波多通道信号处理方法,包括以下步骤:求取四进小波系数,四进小波系数包括四进小波分解系数和四进小波重构系数;对四进小波系数进行优化以降低四进小波重构的误差;将优化后的四进小波系数转化为整数乘法结合移位运算的结构;结合四进小波系数,根据整数乘法结合移位运算的结构,利用FPGA平台对待处理信号进行并行运算,包括四进小波分解、抽取插值和四进小波重构处理,在四进小波重构过程中根据实际情况进行置零、压缩、设置阈值、级联滤波。本发明能对数字信号进行多个频率域的特征性处理,可以通过实际情况提升四进小波滤波器分解层次以达到不同分辨率要求,且运算复杂度低、运算速度快、节省硬件运算资源。
Description
技术领域
本发明属于信号处理领域,特别是一种基于FPGA的四进小波多通道信号处理方法。
背景技术
信号处理的核心是数字滤波,数字滤波器可以分为经典滤波器和现代滤波器两大部分。经典滤波器假设噪声和信号的频谱属于不同的频率相应范围,设置频率域通带和阻带就可以实现噪声去除,常用的高通、低通、带通、带阻滤波器都属于经典滤波器。因为经典滤波系统主要滤除反应电子系统特征的固定频段信号且真实信号往往在该频段也有分布,所以对于高分辨率要求的微弱信号处理,经典滤波无法达到要求,基于统计特性估计出有用信号和噪声信号的现代滤波方法被应用到数字滤波领域。
现代滤波方法的质量与所处理信号本身频率特性密不可分,例如LMS最小均方自适应滤波算法,其本身收敛速度较慢,不适用于快速变化的信号;RLS递推最小二乘算法则具有快速收敛和稳定的特性,对快变信号处理效果优越。既然算法对不同频率域信号的处理效果具有明显差异性,那么将信号分解到不同频率域处理就能够显著改善信号滤波效果,对于微弱信号处理提供有效的算法支持。
FPGA采用并行处理工作模式,多个通道并行计算,解决了将信号分解到不同频率域处理的延时问题,小波变换理论的正变换分解和逆变换重构为多频率通道处理提供了理论支持,完全重构条件解决了实际应用中将信号分解成不同频率通道的通带和过渡带彼此间有重叠的问题,保证信号在不同频率通道处理后可以还原真实信号。
从信号采样角度划分,小波滤波器主要划分为双通道和M(M>2)通道滤波器。二进小波在对称性、线性相位和紧支撑上不能同时兼顾(哈尔小波除外,但是哈尔小波函数失去了光滑性)。在数字信号处理时,滤波器不满足线性相位,则会造成相位失真。为了消除相位失真,一直以来都采用牺牲正交性、用双正交小波的方法来完成滤波。这些局限性使小波逐渐向多小波和四进小波两个方向发展,多小波在发展中遇到和二进小波一样的矛盾。多进小波扩充了高通滤波器子带数目,提升了小波构造的自由度,可以做到同时兼顾对称性、线性相位和双正交,为小波滤波器发展提供了新的空间。但多进小波应用到多通道滤波器存在的问题有:多进小波的系数求取困难大,且除了哈尔阵中的系数其余系数均为无理数;但是哈尔阵缺乏光滑性,难以完成实时信号处理的任务;且在基于定点数运算的FPGA上实现无理数乘法存在字长效应,无法保证完全重构。
发明内容
本发明的目的在于提供一种运算速度快、运算复杂度低、节省硬件运算资源的基于FPGA的四进小波多通道信号处理方法。
实现本发明目的的技术解决方案为:一种基于FPGA的四进小波多通道信号处理方法,包括以下步骤:
步骤1、求取四进小波系数,四进小波系数包括四进小波分解系数和四进小波重构系数;
步骤2、对步骤1获得的四进小波系数进行优化以降低四进小波重构的误差;
步骤3、将优化后的四进小波系数转化为整数乘法结合移位运算的结构;
步骤4、结合步骤1的四进小波系数,根据步骤3的整数乘法结合移位运算的结构,利用FPGA平台对待处理信号进行并行运算,包括四进小波分解、抽取插值和四进小波重构处理,在四进小波重构过程中根据实际情况进行置零、压缩、设置阈值、级联滤波。
本发明与现有技术相比,其显著优点为:1)本发明通过增加约束条件和一个可变因子λ,将繁杂的四进小波求取条件简化至3组方程,降低运算复杂度;2)本发明通过改进二进小波的提升方法,能够获得避免字长效应且保证完全重构的可硬件实现的系数组合,该系数组合为简单的求和与移位运算,不必调用乘法器核,提高了运算速度,节省了硬件运算资源。
下面结合附图对本发明作进一步详细的描述。
附图说明
图1为本发明基于FPGA的四进小波多通道信号处理系统的结构示意图。
图2为本发明实施例的系统结构示意图。
图3为本发明实施例模拟信号处理集成箱直接输出的模拟信号波形图。
图4为本发明实施例中应用本发明信号处理装置和方法获得的数据拟合结果图。
图5为本发明实施例FPGA芯片采用各滤波器滤波效果分析图。
具体实施方式
结合图1,本发明一种基于FPGA的四进小波多通道信号处理方法,包括以下步骤:
步骤1、求取四进小波系数,四进小波系数包括四进小波分解系数和四进小波重构系数;
步骤2、对步骤1获得的四进小波系数进行优化以降低四进小波重构的误差;
步骤3、将优化后的四进小波系数转化为整数乘法结合移位运算的结构;
步骤4、结合步骤1的四进小波系数,根据步骤3的整数乘法结合移位运算的结构,利用FPGA平台对待处理信号进行并行运算,包括四进小波分解、抽取插值和四进小波重构处理,在四进小波重构过程中根据实际情况进行置零、压缩、设置阈值、级联滤波。
进一步地,步骤1中所述求取四进小波系数,具体为:
在现有四进小波系数求取条件的基础上增加约束条件:
1)四进小波分解和重构系数完全相同;
2)弱化双正交条件,四进小波分解或重构的通道长度为4L,L∈Z+;
3)令尺度通道系数呈中心对称,小波通道系数通过改变尺度通道系数的符号和顺序获取;
4)将四进小波分解和重构系数放大λ倍,且使得λ2可以写为2n',n'为整数;
步骤1-1、在现有的四进小波系数的错位型结构的基础上,增加约束条件1),简化双正交条件方程组为:
∑ihihi+4j=λ2δj
其中,j为自然数,j=0时δj=1,j≠0时δj=0;
最终获得简化后的四进小波系数方程组:
∑ihi=2×λ
∑ihihi+4j=λ2δj
由简化后的四进小波系数方程组即可求取出尺度通道系数;
步骤1-2、在步骤1-1四进小波系数方程组的基础上增加消失矩条件,获得新的四进小波系数方程组,其中消失矩条件为:
式中,m=0,1,…,r-1,则{gi (k)}拥有r阶消失矩,其中r的取值等于小波通道系数所需的消失矩数;
步骤1-3、在步骤1-2的基础上,求出的尺度通道系数并非一组解,而是带有变量的一个集合,因此增加约束条件2)约束自由度变量,具体为:去除双正交三个条件中的以弱化双正交条件,由此去除尺度通道系数的非0元素,将4L+2长度的尺度通道系数变为4L长度的尺度通道系数;
步骤1-4、在步骤1-3的基础上,增加约束条件3),具体为:
第一小波通道系数gi (1):
g2p (1)=(-1)L-ph2p+1
g2p+1 (1)=(-1)L-p(-h2p)
式中,p=0、1、2、…、2L-1;
第二小波通道系数gi (2)通过对尺度通道系数的奇数位符号进行取反获得;
第三小波通道系数gi (3)通过对第一小波通道系数的奇数位符号进行取反获得。
进一步地,步骤2所述对步骤1获得的四进小波系数进行优化以降低四进小波重构的误差,具体为:
步骤2-1、由四进小波系数的频率域形式获得四进小波提升算法的频率域形式,具体为:
四进小波系数的频率域形式为:
式中,H(ω)为尺度通道系数频率域,H(ω)=H0(ω)+H1(ω)+H2(ω)+H3(ω),H0(ω)、H1(ω)、H2(ω)、H3(ω)分别为i=4q、4q+1、4q+2、4q+3的尺度通道系数hi的频率域形式,q=0、1、…、L-1;G(k)(ω)为小波通道系数频率域, 分别为i=4q、4q+1、4q+2、4q+3的小波通道系数gi (k)的频率域形式,ω为频率;
用Aan表示四进小波系数的频率域形式的代数余子式,对于尺度通道a=1,由代数余子式代替二进小波提升算法中的G(ω),获得四进小波提升算法的频域形式Hnew(ω)为:
式中,Tn-1(ω)为eiω的多项式;
步骤2-2、利用Laurant多项式,将步骤2-1获得的四进小波提升算法频域形式转换为时域形式Hnew:
式中,H1、H2、……、Hb为不同的尺度通道系数,cs为调节因子且cs 2为整数,s=1,2,…,b;此形式成立则在FPGA硬件平台可以建立H1、H2、……、Hb不同尺度通道的并行模块,b个模块同时进行,最后移位获得输出,避免浮点数运算带来数据丢失。
步骤2-3、对步骤2-1的四进小波提升算法的频域形式Hnew(ω)进行存在性证明,具体为:
步骤2-3-1、令Tn-1(ω)=Hn-1(ω),则
Hnew(ω)=2H(ω)
由此可知Hnew(ω)符合步骤1中的约束条件4);
步骤2-3-2、证明Tn-1(ω)的时域形式Tn-1(z)是否必然存在,具体为:
将步骤2-1中的G(ω)写为z频率域形式:
根据步骤1约束条件1)可得:
GG*=I
则矩阵G的秩为4,因此矩阵G中的任意一个元素的余子式的模不为0,则必然存在:
Tn-1(z)=λn-1H(z)Hn-1(z)/|H(z)|
由此证明Tn-1(z)必然存在,则Hnew(ω)存在。
进一步地,步骤3所述将优化后的四进小波系数转化为整数乘法结合移位运算的结构,具体为:
将步骤2-2公式中的H1、H2、……、Hb转换为与整数相乘的形式,并重新整合该公式,获得新的系数序列H1、H2、……、Hb',由此获得整数乘法结合移位运算的结构。以二进小波系数为例,如果调节因子选择c为c2为2满足要求,新的H2=[11],此时调节因子所以只需经过和的并行运算分解和重构后,还原系数放大倍数λ2=2,即重构分解至滤波结束时,进行一次右移位就可以避免的浮点数系数1.414213562373乘法,保证完全重构且避免字长效应。
下面结合实施例对本发明作进一步详细的描述。
实施例
本实施例对位移传感器的信号处理为例进行分析,系统结构如图2所示,位移传感器输出模拟信号至模拟信号处理集成箱,获得的模拟信号如图3所示,由图可知电压值被完全淹没于噪声中无法读取出位移对应的电压值;AD7606将模拟信号处理集成箱输出位移信号的模拟量转换为16位数字信号;该数字信号经FPGA芯片AC616中完成电压值还原和四进小波滤波后经过串口输出在电脑端显示。其中AD7606电压转换范围±5V,经过四进小波滤波后串口发送回电脑的结果如表1所示:
表1位移信号采样数据
位移/10μm | 1 | 2 | 3 | 4 | 5 |
电压/mV | 1791.06 | 1803.52 | 1817.66 | 1820.02 | 1831.36 |
位移编号 | 6 | 7 | 8 | 9 | 10 |
电压/10μm | 1843.92 | 1850.16 | 1857.75 | 1873.04 | 1881.59 |
图4为本实施例获得的数据拟合结果图,相对于图3未经本发明装置和方法处理直接输出的模拟信号,本发明能够读取出位移传感器直接输出的被噪声淹没的电压值,且由图可知,电压值精度9.726mV,能保证最大量程±5V千分之一的精度。
本实施例中FPGA芯片采用不同滤波器处理固定信噪比数字信号的滤波结果如下表2所示,四进小波采用了2层分解,第一层小波通道置零,第二层小波通道1压缩比例1/8,小波通道2置零,小波通道3设置阈值为8。
表2各滤波器滤波效果
滤波器长度=8 | 含噪信号 | 中值滤波 | 均值滤波 | LMS | 本发明四进小波滤波 |
SNR | 1.0043 | 4.1579 | 5.8617 | 7.1549 | 10.2795 |
SNR | 3.0043 | 5.1979 | 5.7793 | 7.3144 | 11.7722 |
SNR | 5.0043 | 6.6511 | 6.9713 | 8.8359 | 12.1933 |
SNR | 7.0043 | 8.3721 | 7.0476 | 9.3503 | 13.8261 |
SNR | 9.0043 | 10.0189 | 7.4366 | 9.9656 | 14.6712 |
各滤波器滤波效果分析图如图5所示,可以得出本发明多通道滤波在第二层分解时滤波效果明显提升,对于信噪比小于4的高噪信号,处理效果是常规滤波器的2倍。
本发明能对数字信号进行多个频率域的特征性处理,可以通过实际情况提升四进小波滤波器分解层次以达到不同分辨率要求,且运算复杂度低、运算速度快、节省硬件运算资源。
Claims (3)
1.一种基于FPGA的四进小波多通道信号处理方法,其特征在于,包括以下步骤:
步骤1、求取四进小波系数,四进小波系数包括四进小波分解系数和四进小波重构系数;
步骤2、对步骤1获得的四进小波系数进行优化以降低四进小波重构的误差;
步骤3、将优化后的四进小波系数转化为整数乘法结合移位运算的结构;
步骤4、结合步骤1的四进小波系数,根据步骤3的整数乘法结合移位运算的结构,利用FPGA平台对待处理信号进行并行运算,包括四进小波分解、抽取插值和四进小波重构处理,在四进小波重构过程中根据实际情况进行置零、压缩、设置阈值、级联滤波;步骤1中所述求取四进小波系数,具体为:
在现有四进小波系数求取条件的基础上增加约束条件:
1)四进小波分解和重构系数完全相同;
2)弱化双正交条件,四进小波分解或重构的通道长度为4L,L∈Z+;
3)令尺度通道系数呈中心对称,小波通道系数通过改变尺度通道系数的符号和顺序获取;
4)将四进小波分解和重构系数放大λ倍,且使得λ2可以写为2n',n'为整数;
步骤1-1、在现有四进小波系数的错位型结构的基础上,增加约束条件1),简化双正交条件方程组为:
∑ihihi+4j=λ2δj
其中,j为自然数,j=0时δj=1,j≠0时δj=0;
最终获得简化后的四进小波系数方程组:
∑ihi=2×λ
∑ihihi+4j=λ2δj
由简化后的四进小波系数方程组即可求取出尺度通道系数;
步骤1-2、在步骤1-1四进小波系数方程组的基础上增加消失矩条件,获得新的四进小波系数方程组,其中消失矩条件为:
式中,m=0,1,…,r-1,则{gi (k)}拥有r阶消失矩,其中r的取值等于小波通道系数所需的消失矩数;
步骤1-4、在步骤1-3的基础上,增加约束条件3),具体为:
第一小波通道系数gi (1):
g2p (1)=(-1)L-ph2p+1
g2p+1 (1)=(-1)L-p(-h2p)
式中,p=0、1、2、…、2L-1;
第二小波通道系数gi (2)通过对尺度通道系数的奇数位符号进行取反获得;
第三小波通道系数gi (3)通过对第一小波通道系数的奇数位符号进行取反获得。
2.根据权利要求1所述基于FPGA的四进小波多通道信号处理方法,其特征在于,步骤2中所述对步骤1获得的四进小波系数进行优化以降低四进小波重构的误差,具体为:
步骤2-1、由四进小波系数的频率域形式获得四进小波提升算法的频率域形式,具体为:
四进小波系数的频率域形式为:
式中,H(ω)为尺度通道系数频率域,H(ω)=H0(ω)+H1(ω)+H2(ω)+H3(ω),H0(ω)、H1(ω)、H2(ω)、H3(ω)分别为i=4q、4q+1、4q+2、4q+3的尺度通道系数hi的频率域形式,q=0、1、…、L-1;G(k)(ω)为小波通道系数频率域,分别为i=4q、4q+1、4q+2、4q+3的小波通道系数gi (k)的频率域形式,ω为频率;
用Aan表示四进小波系数的频率域形式的代数余子式,对于尺度通道a=1,由代数余子式代替二进小波提升算法中的G(ω),获得四进小波提升算法的频域形式Hnew(ω)为:
式中,Tn-1(ω)为eiω的多项式;
步骤2-2、利用Laurant多项式,将步骤2-1获得的四进小波提升算法频域形式转换为时域形式Hnew:
式中,H1、H2、…、Hb为不同的尺度通道系数,cs为调节因子且cs 2为整数,s=1,2,…,b;
步骤2-3、对步骤2-1的四进小波提升算法的频域形式Hnew(ω)进行存在性证明,具体为:
步骤2-3-1、令Tn-1(ω)=Hn-1(ω),则
Hnew(ω)=2H(ω)
由此可知Hnew(ω)符合步骤1中的约束条件4);
步骤2-3-2、证明Tn-1(ω)的时域形式Tn-1(z)是否必然存在,具体为:
将步骤2-1中的G(ω)写为z频率域形式:
根据步骤1约束条件1)可得:
GG*=I
则矩阵G的秩为4,因此矩阵G中的任意一个元素的余子式的模不为0,则必然存在:
Tn-1(z)=λn-1H(z)Hn-1(z)/|H(z)|
由此证明Tn-1(z)必然存在,则Hnew(ω)存在。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910285083.6A CN110032707B (zh) | 2019-04-10 | 2019-04-10 | 一种基于fpga的四进小波多通道信号处理方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910285083.6A CN110032707B (zh) | 2019-04-10 | 2019-04-10 | 一种基于fpga的四进小波多通道信号处理方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110032707A CN110032707A (zh) | 2019-07-19 |
CN110032707B true CN110032707B (zh) | 2023-05-05 |
Family
ID=67237986
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910285083.6A Active CN110032707B (zh) | 2019-04-10 | 2019-04-10 | 一种基于fpga的四进小波多通道信号处理方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110032707B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113238132B (zh) * | 2021-04-27 | 2024-05-03 | 平顶山学院 | 一种频域介电谱测试仪的检测装置及检测方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102594472A (zh) * | 2012-03-22 | 2012-07-18 | 北京邮电大学 | 一种基于小波分解阈值去噪的无线信道测量的方法及其装置 |
CN108195932A (zh) * | 2017-11-30 | 2018-06-22 | 北京宇航系统工程研究所 | 一种飞行器管路损伤超声导波定量评估方法 |
-
2019
- 2019-04-10 CN CN201910285083.6A patent/CN110032707B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102594472A (zh) * | 2012-03-22 | 2012-07-18 | 北京邮电大学 | 一种基于小波分解阈值去噪的无线信道测量的方法及其装置 |
CN108195932A (zh) * | 2017-11-30 | 2018-06-22 | 北京宇航系统工程研究所 | 一种飞行器管路损伤超声导波定量评估方法 |
Non-Patent Citations (1)
Title |
---|
A Digital Filtering Algorithm Based on Four-channel Wavelet and its Application in Active Optical System;Qiu Guangbo等;《SPIE》;20190118;第1-11页 * |
Also Published As
Publication number | Publication date |
---|---|
CN110032707A (zh) | 2019-07-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Kormylo et al. | Two-pass recursive digital filter with zero phase shift | |
WO1984002019A1 (en) | Image processing method using a block overlap transformation procedure | |
CN110784229B (zh) | 一种基于快速傅里叶变换的带有模拟滤波器补偿的mwc后端信号重构方法 | |
US4791597A (en) | Multiplierless FIR digital filter with two to the Nth power coefficients | |
US3872290A (en) | Finite impulse response digital filter with reduced storage | |
EP0612148B1 (en) | Digital filtering circuit operable as a three-stage moving average filter | |
CN110032707B (zh) | 一种基于fpga的四进小波多通道信号处理方法 | |
CN107294512B (zh) | 一种基于树型结构的非均匀滤波器组滤波方法 | |
CN106921367A (zh) | 一种sigma delta ADC的数字抽取滤波器 | |
CN101242168B (zh) | 一种fir数字滤波器直接型实现方法及实现装置 | |
CN113037248A (zh) | 一种基于分段频域优化的小数延时滤波器设计方法 | |
US20060224649A1 (en) | Interpolation and decimation using newton polyphase filters | |
CN105099398B (zh) | 基于相位调制的非均匀dft调制滤波器组的构建方法 | |
CN115296647A (zh) | 基于frm技术的wola滤波器组及子带分割设计方法 | |
CN101072018B (zh) | 数字信号的分频滤波方法及其系统 | |
CN107302358B (zh) | 一种四通道tiadc的非线性失配补偿方法 | |
CN102571031B (zh) | 一种用于实现小波变换的宽带低插损声表面波滤波器组 | |
King et al. | A unified approach to scrambling filter design | |
WO2014064408A2 (en) | Method and apparatus for processing a signal | |
CN114584109B (zh) | 一种对高速串行接口传输的并行信号进行实时滤波的方法 | |
CN111010146B (zh) | 一种基于快速滤波器组的信号重构结构及其设计方法 | |
CN113708741B (zh) | 一种系数可配置数字滤波器的设计方法及装置 | |
CN107454030B (zh) | 一种电力线宽带载波半并行发射机及其实现方法 | |
Asemani et al. | Time-division multiplexing architecture for hybrid filter bank A/D converters | |
Wang et al. | Adaptive analog signal sampling through lagrange interpolation |
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 |