CN107294534A - 用于窄带信号采样的双通道tiadc频响失配实时校正方法 - Google Patents

用于窄带信号采样的双通道tiadc频响失配实时校正方法 Download PDF

Info

Publication number
CN107294534A
CN107294534A CN201710339024.3A CN201710339024A CN107294534A CN 107294534 A CN107294534 A CN 107294534A CN 201710339024 A CN201710339024 A CN 201710339024A CN 107294534 A CN107294534 A CN 107294534A
Authority
CN
China
Prior art keywords
frequency response
mrow
error
narrow band
mismatch
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
CN201710339024.3A
Other languages
English (en)
Other versions
CN107294534B (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.)
Sun Yat Sen University
SYSU CMU Shunde International Joint Research Institute
Original Assignee
Sun Yat Sen University
SYSU CMU Shunde International Joint Research Institute
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 Sun Yat Sen University, SYSU CMU Shunde International Joint Research Institute filed Critical Sun Yat Sen University
Priority to CN201710339024.3A priority Critical patent/CN107294534B/zh
Publication of CN107294534A publication Critical patent/CN107294534A/zh
Application granted granted Critical
Publication of CN107294534B publication Critical patent/CN107294534B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • HELECTRICITY
    • H03ELECTRONIC CIRCUITRY
    • H03MCODING; DECODING; CODE CONVERSION IN GENERAL
    • H03M1/00Analogue/digital conversion; Digital/analogue conversion
    • H03M1/10Calibration or testing
    • H03M1/1009Calibration

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Compression, Expansion, Code Conversion, And Decoders (AREA)
  • Analogue/Digital Conversion (AREA)

Abstract

本发明提供了一种用于窄带信号采样的双通道时间交织模数转换器(TIADC)频响失配实时校正方法,该方法采用多项式形式的频率响应函数表征系统的线性失配特性,对窄带输入信号人为叠加高斯白噪声作为TIADC的输入信号并利用轻微过采样获取失配信息,基于变步长最小均方误差(VSSLMS)算法对失配误差进行实时地边估计边校正。VSSLMS算法相对于基本的LMS算法,前者在大的误差范围内有快速收敛性,在小的误差范围内有较小的失调量,提高了跟踪性能。得到校正后的输出再通过设计的一个滤波器将白噪声滤除从而得到期望输出。该方法解决了窄带输入信号无法保证在过采样带上出现误差能量的问题,保证系统参数能够收敛到精确值从而实现有效校正,并且简单易行,校正效果好。

Description

用于窄带信号采样的双通道TIADC频响失配实时校正方法
技术领域
本发明涉及信号采样与处理技术领域,更具体地,涉及一种用于窄带信号采样的双通道TIADC频响失配实时校正方法。
背景技术
随着集成电路技术的不断发展,数字化技术的推广,对模数转换器件ADC的采样速率以及采样精度的要求越来越高,不仅要求数据采集系统有高的采样率,还要有高的采样精度。在实际的运用中,对实时采样速率以及采样精度有极高的依赖性。然而ADC的最大采样速率受限于它的分辨率,分辨率与采样速率之间是一对矛盾体,高采样速率要求较短的转换时间,而高分辨率则要求较长的转换时间。根据目前的IC设计工艺,要实现更高速的采样速率,需要开发一种基于新结构和新方法的ADC模块。现有技术所提供的能够实现超高速采样的系统就是利用时间交织(Time-interleaved)结构的ADC系统。
这种结构的ADC系统利用M片有着相同采样率fs的单个ADC模块,采用并行的结构,每片ADC模块以相隔1/(M*fs)的时间间隔进行采样,以达到采样率为M*fs(总采样率f=M*fs)的效果。理论上,这种M通道并行交替采样的ADC系统能够使得整个系统的采样率达到单个ADC模块的M倍。但是由于制造工艺本身固有的缺点,不可能使得每一片ADC模块完全一模一样,所以必然会使得各个通道的ADC模块之间存在失配误差,且每片ADC自身带有微分和积分非线性特性,从而严重降低了整个ADC系统的信噪比。
目前,大多数方法主要针对线性失配,例如增益误差,时间误差等进行估计和校正,部分方法针对模数转换器(ADC)自身的积分和微分非线性造成的失配进行估计和校正。由于基于通道传递函数的校正方法能够把任何线性误差的效果转移为频域响应失配误差(frequency-response mismatch errors),对线性失配的校正效果更优。
为了在频谱上获取只含有失配误差信号的频带,基于多项式频响失配的校正结构需要对输入信号进行过采样,然而对于窄带信号不能保证失配误差出现在过采样带上,从而使校正算法失效。
发明内容
本发明为解决现有技术在对窄带信号的采样结果进行校正时由于窄带信号不能保证失配误差出现在过采样带上而导致的校正算法失效的缺陷,提供了一种用于窄带信号采样的双通道TIADC频响失配实时校正方法,该校正方法通过在窄带信号上人为叠加高斯白噪声作为2-TIADCs系统的输入信号,从而利用轻微过采样获取失配误差信号,同时基于变步长最小均方误差(VSSLMS)算法对失配误差信号进行实时地边估计边校正,得到校正后的输出再设计一个滤波器将高斯白噪声滤除从而得到期望输出。
为实现以上发明目的,采用的技术方案是:
用于窄带信号采样的双通道TIADC频响失配实时校正方法,包括以下步骤:
S1.设置输入的窄带信号s(t)满足奈奎斯特采样定理,并且使系统可对其进行轻微的过采样;
S2.将高斯白噪声通过一个低通滤波器对其采用过采样,并叠加到窄带信号s(t)上作为2-TIADCs系统的输入;
S3.设计一级微分器d1(n),通过时域卷积获取所需的更高级的微分器dp(n),2≤p≤P;
S4.采用多项式形式的频率响应函数对线性频响失配进行建模:
确定通道频率响应函数的阶数P,令cp为p阶多项式系数,系统的归一化频率响应函数为线性频响失配误差
S5.将2-TIADCs系统的输出y(n)通过步骤S3中设计的各级微分滤波器,再经过一个(-1)n的乘法器,所得信号为Yd(n)=[y0(n),...,yp(n),...,yP(n)],其中:
yp(n)=dp(n)*y(n)(-1)n
S6.令某时刻的误差参数为设计一个VSSLMS算法对误差参数进行迭代更新;
S7.利用步骤S6中迭代更新所得的参数对线性频响失配误差进行重构得到
其中Yd(n)T为步骤S5中所得的信号Yd(n)的转置;
S8.将2-TIADCs系统的输出y(n)减去步骤S7中重构的估计误差得到校正输出
S9.设计一个LMS滤波器将步骤S8所得的校正输出中的高斯白噪声滤除,得到期望输出
优选地,所述步骤S3中使用到的微分器为线性相位数字微分器。
优选地,所述步骤S6利用VSSLMS算法对误差参数进行迭代更新的具体过程如下:
设计相应的高通滤波器f(n),使高通滤波器f(n)的截止频率高于采样信号的截止频率,令迭代更新的公式如下:
μ(n+1)=α·μ(n)+γ·ε2(n)
其中Yd f(n)=Yd(n)*f(n)=[y0(n)*f(n)...,yp(n)*f(n)...,yP(n)*f(n)]T,α是一个接近1的系数,γ是一个接近零的系数。
与现有技术相比,本发明的有益效果是:
本发明提供了一种用于窄带信号采样的双通道时间交织模数转换器(TIADC)频响失配实时校正方法,该方法采用多项式形式的频率响应函数表征系统的线性失配特性,对窄带输入信号人为叠加高斯白噪声作为TIADC的输入信号并利用轻微过采样获取失配信息,基于变步长最小均方误差(VSSLMS)算法对失配误差进行实时地边估计边校正。VSSLMS算法相对于基本的LMS算法,前者在大的误差范围内有快速收敛性,在小的误差范围内有较小的失调量,提高了跟踪性能。得到校正后的输出再通过设计的一个滤波器将白噪声滤除从而得到期望输出。该方法解决了窄带输入信号无法保证在过采样带上出现误差能量的问题,保证系统参数能够收敛到精确值从而实现有效校正,并且简单易行,校正效果好。
附图说明
图1为时间交织模数转换器的结构示意图。
图2为基于通道频率响应的双通道TIADC模型示意图。
图3为误差校正的原理框图。
图4为基于VSSLMS算法的自适应边估计边校正的示意图。
图5为校正方法的流程图。
具体实施方式
附图仅用于示例性说明,不能理解为对本专利的限制;
以下结合附图和实施例对本发明做进一步的阐述。
实施例1
如图5所示,用于窄带信号采样的双通道TIADC频响失配实时校正方法,包括以下步骤:
S1.设置输入的窄带信号s(t)满足奈奎斯特采样定理,并且使系统可对其进行轻微的过采样;
S2.将高斯白噪声通过一个低通滤波器对其采用过采样,并叠加到窄带信号s(t)上作为2-TIADCs系统的输入;
S3.设计一级微分器d1(n),通过时域卷积获取所需的更高级的微分器dp(n),2≤p≤P;
S4.采用多项式形式的频率响应函数对线性频响失配进行建模:
确定通道频率响应函数的阶数P,令cp为p阶多项式系数,系统的归一化频率响应函数为线性频响失配误差
S5.将2-TIADCs系统的输出y(n)通过步骤S3中设计的各级微分滤波器,再经过一个(-1)n的乘法器,所得信号为Yd(n)=[y0(n),...,yp(n),...,yP(n)],其中:
yp(n)=dp(n)*y(n)(-1)n
S6.令某时刻的误差参数为设计一个VSSLMS算法对误差参数进行迭代更新;
S7.利用步骤S6中迭代更新所得的参数对线性频响失配误差进行重构得到
其中Yd(n)T为步骤S5中所得的信号Yd(n)的转置;
S8.将2-TIADCs系统的输出y(n)减去步骤S7中重构的估计误差得到校正输出
S9.设计一个LMS滤波器将步骤S8所得的校正输出中的高斯白噪声滤除,得到期望输出
在具体的实施过程中,所述步骤S3中使用到的微分器为线性相位数字微分器。
在具体的实施过程中,所述步骤S6利用VSSLMS算法对误差参数进行迭代更新的具体过程如下:
设计相应的高通滤波器f(n),使高通滤波器f(n)的截止频率高于采样信号的截止频率,令迭代更新的公式如下:
μ(n+1)=α·μ(n)+γ·ε2(n)
其中Yd f(n)=Yd(n)*f(n)=[y0(n)*f(n)...,yp(n)*f(n)...,yP(n)*f(n)]T,α是一个接近1的系数,γ是一个接近零的系数。
实施例1
本实施例在实施例1的基础上,进行了具体的实验:
本实施例的实验采用的理想输入信号是多正弦信号,频率为f1=0.1fs,f2=0.45fs,f3=0.7fs,其中fs为采样频率,叠加一个均值为0,方差为1的高斯白噪声。如图1所示为时间交织模数转换器的结构示意图,输入信号以M通道输入,每条通道以相同的采样率但不同的采样时刻(相邻通道相差Ts时刻)对高速输入信号采样,最终合并出输出信号,以此实现高速采样的模数转化。图2为基于通道频率响应的双通道TIADC模型示意图,该种模型能够把任何线性误差转移到通道传递函数的参数上,可用统一的方法进行补偿。设置阶数P=3,系统参数为c=[-0.025,0.005,-0.0015,-0.0001]。
如实施例1所述,本发明提供的校正方法需要用到微分器。本实施例中,所使用的微分器为线性相位数字微分器。设计一级的微分滤波器,长度为N=41,在通过时域卷积运算获得长度相等的二级和三级微分滤波器,对应的滤波器系数为:
利用TIADC的输出和上述的微分滤波器构造
Yd(n)=[y0(n),y1(n),y2(n),y3(n)],
其中:
如图4所示的自适应边估计边校正示意图,结合实施例1中的步骤S6和S8,对系统的参数进行迭代更新,系统收敛后获得估计参数为:
可见频响失配误差参数在误差允许范围内可以稳定收敛到真值。信号未经过校正前,存在大量的噪声毛刺,其幅度最高可以达到-40dB。而通过校正且滤除高斯白噪声后,噪声频谱受到抑制,其最高的毛刺的幅度可降低到-80dB左右。
从以上的实验结果可以得出,本发明解决了窄带输入信号无法保证在过采样带上出现误差能量的问题,保证系统参数能够收敛到精确值从而实现有效补偿,并且简单易行,补偿效果好。而且选用VSSLMS算法相对于基本的LMS算法,前者在大的误差范围内有快速收敛性,在小的误差范围内有较小的失调量,提高了跟踪性能。
显然,本发明的上述实施例仅仅是为清楚地说明本发明所作的举例,而并非是对本发明的实施方式的限定。对于所属领域的普通技术人员来说,在上述说明的基础上还可以做出其它不同形式的变化或变动。这里无需也无法对所有的实施方式予以穷举。凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明权利要求的保护范围之内。

Claims (3)

1.用于窄带信号采样的双通道TIADC频响失配实时校正方法,其特征在于:包括以下步骤:
S1.设置输入的窄带信号s(t)满足奈奎斯特采样定理,并且使系统可对其进行轻微的过采样;
S2.将高斯白噪声通过一个低通滤波器对其采用过采样,并叠加到窄带信号s(t)上作为2-TIADCs系统的输入;
S3.设计一级微分器d1(n),通过时域卷积获取所需的更高级的微分器dp(n),2≤p≤P;
S4.采用多项式形式的频率响应函数对线性频响失配进行建模:
确定通道频率响应函数的阶数P,令cp为p阶多项式系数,系统的归一化频率响应函数为线性频响失配误差
S5.将2-TIADCs系统的输出y(n)通过步骤S3中设计的各级微分滤波器,再经过一个(-1)n的乘法器,所得信号为Yd(n)=[y0(n),...,yp(n),...,yP(n)],其中:
yp(n)=dp(n)*y(n)(-1)n
S6.令某时刻的误差参数为设计一个VSSLMS算法对误差参数进行迭代更新;
S7.利用步骤S6中迭代更新所得的参数对线性频响失配误差进行重构得到
其中Yd(n)T为步骤S5中所得的信号Yd(n)的转置;
S8.将2-TIADCs系统的输出y(n)减去步骤S7中重构的估计误差得到校正输出
S9.设计一个LMS滤波器将步骤S8所得的校正输出中的高斯白噪声滤除,得到期望输出
2.根据权利要求1所述的用于窄带信号采样的双通道TIADC频响失配实时校正方法,其特征在于:所述步骤S3中使用到的微分器为线性相位数字微分器。
3.根据权利要求1所述的用于窄带信号采样的双通道TIADC频响失配实时校正方法,其特征在于:所述步骤S6利用VSSLMS算法对误差参数进行迭代更新的具体过程如下:
设计相应的高通滤波器f(n),使高通滤波器f(n)的截止频率高于采样信号的截止频率,令迭代更新的公式如下:
<mrow> <mover> <mi>c</mi> <mo>^</mo> </mover> <mrow> <mo>(</mo> <mi>n</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> <mo>=</mo> <mover> <mi>c</mi> <mo>^</mo> </mover> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>+</mo> <mi>&amp;mu;</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> <mo>&amp;CenterDot;</mo> <mi>&amp;epsiv;</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> <mo>&amp;CenterDot;</mo> <msup> <msub> <mi>Y</mi> <mi>d</mi> </msub> <mi>f</mi> </msup> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </mrow>
μ(n+1)=α·μ(n)+γ·ε2(n)
其中Yd f(n)=Yd(n)*f(n)=[y0(n)*f(n)...,yp(n)*f(n)...,yP(n)*f(n)]T,α是一个接近1的系数,γ是一个接近零的系数。
CN201710339024.3A 2017-05-15 2017-05-15 用于窄带信号采样的双通道tiadc频响失配实时校正方法 Active CN107294534B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710339024.3A CN107294534B (zh) 2017-05-15 2017-05-15 用于窄带信号采样的双通道tiadc频响失配实时校正方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710339024.3A CN107294534B (zh) 2017-05-15 2017-05-15 用于窄带信号采样的双通道tiadc频响失配实时校正方法

Publications (2)

Publication Number Publication Date
CN107294534A true CN107294534A (zh) 2017-10-24
CN107294534B CN107294534B (zh) 2020-10-23

Family

ID=60094625

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710339024.3A Active CN107294534B (zh) 2017-05-15 2017-05-15 用于窄带信号采样的双通道tiadc频响失配实时校正方法

Country Status (1)

Country Link
CN (1) CN107294534B (zh)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108366037A (zh) * 2018-01-24 2018-08-03 佛山市顺德区中山大学研究院 基于iq通信收发系统的多通道tiadc失真校正方法
CN109557512A (zh) * 2018-12-06 2019-04-02 航天南湖电子信息技术股份有限公司 一种高灵敏度和高动态范围的雷达接收机
CN110557122A (zh) * 2019-09-25 2019-12-10 电子科技大学 一种tiadc系统频响非一致性误差的校正方法
CN111384950A (zh) * 2018-12-31 2020-07-07 特克特朗尼克公司 用于时间交织数模转换器的线性和非线性校准
CN115459770A (zh) * 2022-09-19 2022-12-09 中国电子科技集团公司第十四研究所 一种多通道时间交织adc采样时间失配的校正方法
CN118399959A (zh) * 2024-06-28 2024-07-26 西安欣创电子技术有限公司 系统函数的数据处理方法、电子设备及存储介质

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101783683A (zh) * 2008-12-29 2010-07-21 英特赛尔美国股份有限公司 双通道时间交错型模数转换器中的误差估计与校正
US20130015990A1 (en) * 2010-08-17 2013-01-17 Texas Instruments Incorporated Track and hold architecture with tunable bandwidth
CN105024696A (zh) * 2015-07-02 2015-11-04 大唐微电子技术有限公司 多通道并行模数转换系统采样时间误差的校准装置及方法
CN106374920A (zh) * 2016-09-05 2017-02-01 中山大学 一种基于多项式模型的tiadc系统的估计与补偿实现方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101783683A (zh) * 2008-12-29 2010-07-21 英特赛尔美国股份有限公司 双通道时间交错型模数转换器中的误差估计与校正
US20130015990A1 (en) * 2010-08-17 2013-01-17 Texas Instruments Incorporated Track and hold architecture with tunable bandwidth
CN105024696A (zh) * 2015-07-02 2015-11-04 大唐微电子技术有限公司 多通道并行模数转换系统采样时间误差的校准装置及方法
CN106374920A (zh) * 2016-09-05 2017-02-01 中山大学 一种基于多项式模型的tiadc系统的估计与补偿实现方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
SAIHUA XU: "Recursive Filters of Time-Interleaved ADC Mismatch Compensation", 《IEEE TRANSACTIONS ON CIRCUITS AND SYSTEMS》 *
聂伟 等: "一种TIADC系统时钟偏斜误差的全并行校正方法", 《实验技术与管理》 *

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108366037A (zh) * 2018-01-24 2018-08-03 佛山市顺德区中山大学研究院 基于iq通信收发系统的多通道tiadc失真校正方法
CN108366037B (zh) * 2018-01-24 2020-09-15 佛山市顺德区中山大学研究院 基于iq通信收发系统的多通道tiadc失真校正方法
CN109557512A (zh) * 2018-12-06 2019-04-02 航天南湖电子信息技术股份有限公司 一种高灵敏度和高动态范围的雷达接收机
CN111384950A (zh) * 2018-12-31 2020-07-07 特克特朗尼克公司 用于时间交织数模转换器的线性和非线性校准
CN110557122A (zh) * 2019-09-25 2019-12-10 电子科技大学 一种tiadc系统频响非一致性误差的校正方法
CN110557122B (zh) * 2019-09-25 2022-04-19 电子科技大学 一种tiadc系统频响非一致性误差的校正方法
CN115459770A (zh) * 2022-09-19 2022-12-09 中国电子科技集团公司第十四研究所 一种多通道时间交织adc采样时间失配的校正方法
CN115459770B (zh) * 2022-09-19 2024-04-16 中国电子科技集团公司第十四研究所 一种多通道时间交织adc采样时间失配的校正方法
CN118399959A (zh) * 2024-06-28 2024-07-26 西安欣创电子技术有限公司 系统函数的数据处理方法、电子设备及存储介质
CN118399959B (zh) * 2024-06-28 2024-08-23 西安欣创电子技术有限公司 系统函数的数据处理方法、电子设备及存储介质

Also Published As

Publication number Publication date
CN107294534B (zh) 2020-10-23

Similar Documents

Publication Publication Date Title
CN107294534A (zh) 用于窄带信号采样的双通道tiadc频响失配实时校正方法
CN107302357A (zh) 一种双通道tiadc线性频响失配和非线性失配的联合校正方法
CN108111169A (zh) 一种四通道tiadc线性失配和非线性失配的联合校正方法
CN105024696B (zh) 多通道并行模数转换系统采样时间误差的校准装置及方法
CN107124141B (zh) 基于复数域矩阵数值求解自适应误差验证的数字预失真的方法
CN104393872A (zh) 一种多通道并行adc系统的采样时间误差校正方法
JP3752237B2 (ja) A/d変換装置
CN101656538A (zh) 一种基于拉格朗日插值方法的时间交替模拟数字转换(tiadc)系统时间失配实时补偿算法
CN107342750B (zh) 适用于多奈奎斯特区的分数延迟优化方法及其实现结构
CN102801425A (zh) 一种基于双压控振荡器环路的Sigma-Delta模数转换器
Han et al. An All-Digital Background Calibration Technique for M-Channel Downsampling Time-Interleaved ADCs Based on Interpolation
CN110708039B (zh) 一种farrow滤波器的系数拟合方法
CN105207646A (zh) 交替dft调制滤波器组的极小化极大设计方法
Ta et al. Fully digital background calibration technique for channel mismatches in TIADCs
CN103326724A (zh) 一种新的针对tiadc系统时间误差的估计方法
Li et al. Two-stage difference-based estimation method for timing skew in TI-ADCs
CN107171665A (zh) 带通信号的双通道tiadc非线性系统参数估计方法
CN114189229B (zh) 一种基于自适应分段算法的小数延时滤波方法
CN112422103A (zh) 减少升余弦滤波器乘法器数量的方法及fir升余弦滤波器
CN107154804B (zh) 低通信号的双通道tiadc非线性系统参数估计方法
Liu et al. Adaptive nonlinear mismatch calibration technique for TIADC based on memory polynomial model
CN108736891A (zh) 一种iq和tiadc失配校正系统及方法
Liu et al. A novel oversampling scheme for design of hybrid filter bank based ADCs
Zhang et al. Analysis and simulation of distributed T/H with averaging technique in folding and interpolating A/D converters
Guiqing et al. Real-time blind calibration for frequency response mismatches in two-channel TI-ADCs

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