CN108809882A - 一种基于相位幅度特性的软判决方法 - Google Patents

一种基于相位幅度特性的软判决方法 Download PDF

Info

Publication number
CN108809882A
CN108809882A CN201810612762.5A CN201810612762A CN108809882A CN 108809882 A CN108809882 A CN 108809882A CN 201810612762 A CN201810612762 A CN 201810612762A CN 108809882 A CN108809882 A CN 108809882A
Authority
CN
China
Prior art keywords
soft decision
difference
magnitude
value
phase
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
CN201810612762.5A
Other languages
English (en)
Other versions
CN108809882B (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.)
Fujian Xinghai Communication Technology Co Ltd
Original Assignee
Fujian Xinghai Communication Technology 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 Fujian Xinghai Communication Technology Co Ltd filed Critical Fujian Xinghai Communication Technology Co Ltd
Priority to CN201810612762.5A priority Critical patent/CN108809882B/zh
Priority to CN202110255538.7A priority patent/CN113055330B/zh
Publication of CN108809882A publication Critical patent/CN108809882A/zh
Application granted granted Critical
Publication of CN108809882B publication Critical patent/CN108809882B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00Modulated-carrier systems
    • H04L27/26Systems using multi-frequency codes
    • H04L27/2601Multicarrier modulation systems
    • H04L27/2647Arrangements specific to the receiver only
    • H04L27/2649Demodulators
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00Modulated-carrier systems
    • H04L27/18Phase-modulated carrier systems, i.e. using phase-shift keying
    • H04L27/22Demodulator circuits; Receiver circuits
    • H04L27/233Demodulator circuits; Receiver circuits using non-coherent demodulation
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02DCLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
    • Y02D30/00Reducing energy consumption in communication networks
    • Y02D30/70Reducing energy consumption in communication networks in wireless communication networks

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Digital Transmission Methods That Use Modulated Carrier Waves (AREA)

Abstract

本发明公开了一种基于相位幅度特性的软判决方法,包括以下步骤:S1、将星座图分割为16份,根据解调出的子载波数据实部和虚部确定所述子载波数据在所述星座图中的相位块;划分规则为:最高位与星座点实部的符号位相同;第二位为:实部与虚部符号位的异或值;第三位为:若在1、3象限中实部绝对值大于虚部绝对值则为1,否则,为0;在2、4象限中相反;第四位为:根据星座点是否落在象限中间两块得到的值与第三位的异或,其中,若星座点落在象限中间两块,则值为1,否则,值为0。S2、通过比较前后两帧同一个子载波数据所处的相位块,得到对应所述子载波数据的软判决等级。该方法具有复杂程度低及节约FPGA资源等优点。

Description

一种基于相位幅度特性的软判决方法
技术领域
本发明涉及通信技术领域,具体涉及一种基于相位幅度的软判决方法。
背景技术
在正交频分复用技术(Orthogonal Frequency Division Multiplexing,OFDM)中,通常采用相移键控(Phase-Shift Keying,PSK)的方式对子载波进行调制,而在接收端需要通过快速傅立叶变换恢复出各子载波上调制的数据;而对于非相干解调,则通常采用差分相移键控(Differential Phase Shift Keying,DPSK)的调制方式。在PSK或DPSK调制中,由于噪声、多径、相偏以及频偏的影响,均会出现相位的偏转(如图1所示),使得解调容易产生误码。若是精确计算出所有子载波的相位,需要用到的坐标旋转数字计算方法(Coordinate Rotation Digital Computer,CORDIC算法)对于可编程门阵列(FieldProgrammable Gate Array,FPGA)来说实现起来数据处理量又太大。
发明内容
本发明所要解决的技术问题是:提供一种低复杂度且能节约现场FPGA资源的基于相位幅度特性的软判决方法。
为了解决上述技术问题,本发明采用的技术方案为:一种基于相位幅度特性的软判决方法,包括以下步骤:
S1、将星座图分割为16份,根据解调出的子载波数据的实部和虚部确定所述子载波数据在所述星座图中的相位块;
S2、通过比较前后两帧同一个子载波数据所处的相位块,得到对应所述子载波数据的软判决等级。
从上述描述可知,本发明的有益效果在于:本发明采用软判决法,不需要精确的相位信息,本方法在不降低后级维特比译码性能的情况下,大大降低了算法复杂度,减少了FPGA资源的消耗,提高了算法实时性。本发明方法能够很好的支持各种调制模式的星座解调,可以适用于对任何格雷星座映射进行软判决,极大的提高了软判决的速度和对各种模式的星座映射的通用性。
附图说明
图1为现有技术中由于各方面的影响发生偏移后的星座图;
图2为采用本发明方法划分的相位分块示意图;
图3为本发明实施例1~2的软判决逻辑流程图。
具体实施方式
为详细说明本发明的技术内容、所实现目的及效果,以下结合实施方式并配合附图予以说明。
一种基于相位幅度特性的软判决方法,包括以下步骤:
S1、将星座图分割为16份,根据解调出的子载波数据的实部和虚部确定所述子载波数据在所述星座图中的相位块;
S2、通过比较前后两帧同一个子载波数据所处的相位块,得到对应所述子载波数据的软判决等级。
从上述描述可知,本发明的有益效果在于:本发明采用软判决法,不需要精确的相位信息,本方法在不降低后级维特比译码性能的情况下,大大降低了算法复杂度,减少了FPGA资源的消耗,提高了算法实时性。本发明方法能够很好的支持各种调制模式的星座解调,可以适用于对任何格雷星座映射进行软判决,极大的提高了软判决的速度和对各种模式的星座映射的通用性。
进一步地,所述步骤S1中将星座图按以下方式分为十六份,划分规则为:
R1、最高位与星座点实部的符号位相同;
R2、第二位为:实部与虚部符号位的异或值;
R3、第三位为:若在1、3象限中实部绝对值大于虚部绝对值则为1,否则,为0;在2、4象限中相反;
R4、第四位为:根据星座点是否落在象限中间两块得到的值与第三位的异或,其中,若星座点落在象限中间两块,则值为1,否则,值为0。
从上述描述可知,本发明的有益效果在于:该划分规则使得象限分块值是连续增加的,方便后续计算。
进一步地,所述步骤S1中采用非相干解调方式解调子载波数据的实部与虚部。
进一步地,所述步骤S2的具体操作为:
S201、接收到第一帧时,以第一帧实部的符号位作为输出,并储存其相位块的值;
S202、当接收到第二帧时,将第二帧与储存值相减取绝对值,得到相位块的差值;
S203、当差值为0、1、15或7、8、9时,则根据差值输出不同的软判决等级并将新的相位块值储存;当差值为2、3、4、12、13、14或5、6、10、11时,不储存相位块值,并根据差值输出不同的软判决等级。
进一步地,所述步骤S203中根据差值输出的软判决等级具体标准为:当差值为0、1或15时,软判决等级为111;
当差值为2或14时,软判决等级为110;
当差值为3或13时,软判决等级为101;
当差值为4时,软判决等级为100;
当差值为12时,软判决等级为011;
当差值为5或11时,软判决等级为010;
当差值为6或10时,软判决等级为001;
当差值为7、8或9时,软判决等级为000。
进一步地,所述步骤S2还包括对子载波能量值的判断步骤,当实部和虚部都小于一阈值时,软判决等级减小一级。
进一步地,所述方法还包括S3、将所述步骤S2得到的软判决等级通过维特比译码后解码出实际数据。
进一步地,所述步骤S3为循环对每帧每个子载波使用步骤S2进行软判决,得到所有子载波上每比特的软判决等级,并将得到的软判决等级作为后续维特比译码的输入进行译码。
请参照图2~3所示,本发明的实施例一为:一种基于相位幅度特性的软判决方法,包括以下步骤:
S1、将星座图分割为16份,即每个象限4份,根据解调出的子载波数据的实部和虚部确定所述子载波数据在所述星座图中的相位块;
S2、通过比较前后两帧同一个子载波数据所处的相位块,得到对应所述子载波数据的软判决等级;
S3、将所述步骤S2得到的软判决等级通过维特比译码后解码出实际数据。
其中,所述步骤S1中将星座图按以下方式分为十六份,划分规则为:R1、最高位与星座点实部的符号位相同,即1、4象限为0,2、3象限为1;
R2、第二位为:实部与虚部符号位的异或值,即1、3象限为0,2、4象限为1;
R3、第三位为:若在1、3象限中实部绝对值大于虚部绝对值则为1,实部绝对值小于或等于虚部绝对值则为0;基在2、4象限中若实部绝对值小于或等于虚部绝对值则为1,否则为0;
R4、第四位为:根据星座点是否落在象限中间两块得到的值与第三位的异或,其中,若星座点落在象限中间两块,则值为1,否则,值为0。
划分后的结果如图2所示。
具体地,所述星座点落在象限中间两块的判决标准如下:如第一象限中星座点在1和2相位块时则为1,在0和3时为则为0。
以单个子载波为例,采用2DPSK的调制方式,在接收端进行非相干解调。所述步骤S2的具体操作为:
S201、接收到第一帧时,以第一帧实部的符号位作为输出,并储存其相位块的值;
S202、当接收到第二帧时,将第二帧与储存值相减取绝对值,得到相位块的差值;
S203、当差值在0、1、15或7、8、9时,说明相位差在0°和180°附近,则将新的相位块值储存,并根据其相位差输出111或000的软判决等级,以此循环;当差值在2、3、4、12、13、14或5、6、10、11时,由于偏移太大,不储存相位块值,并根据差值大小输出从110至100或001至011中的软判决等级。具体地,所述步骤S203中根据差值输出的软判决等级具体标准为:当差值为0、1或15时,软判决等级为111;
当差值为2或14时,软判决等级为110;
当差值为3或13时,软判决等级为101;
当差值为4时,软判决等级为100;
当差值为12时,软判决等级为011;
当差值为5或11时,软判决等级为010;
当差值为6或10时,软判决等级为001;
当差值为7、8或9时,软判决等级为000;
其中,判决值1的软判决等级从111、110、101到100依次递减,判决值0的软判决等级从000、001、010到011递减。当差值为0/1/15或7/8/9时候相当于两帧数据的相位差在0度或180度左右,此时的判决数据是比较可靠的,因此,判决等级最高为111或000。当差值越接近4或12时,说明相位差越接近90度,则表示判决数据是不可靠的,因此判决等级为最低的100或011。
本发明的实施例二为:一种基于相位幅度特性的软判决方法,包括以下步骤:
S1、将星座图分割为16份,即每个象限4份,根据解调出的子载波数据的实部和虚部确定所述子载波数据在所述星座图中的相位块;
S2、通过比较前后两帧同一个子载波数据所处的相位块,得到对应所述子载波数据的软判决等级;
S3、将所述步骤S2得到的软判决等级通过维特比译码后解码出实际数据。
其中,所述步骤S1中将星座图按以下方式分为十六份,划分规则为:
R1、最高位与星座点实部的符号位相同,即1、4象限为0,2、3象限为1;
R2、第二位为:实部与虚部符号位的异或值,即1、3象限为0,2、4象限为1;
R3、第三位为:若在1、3象限中实部绝对值大于虚部绝对值则为1,实部绝对值小于或等于虚部绝对值则为0;基在2、4象限中若实部绝对值小于或等于虚部绝对值则为1,否则为0;
R4、第四位为:根据星座点是否落在象限中间两块得到的值与第三位的异或,其中,若星座点落在象限中间两块,则值为1,否则,值为0。
划分后的结果如图2所示。
所述步骤S2的具体操作为:
S201、接收到第一帧时,采用非相干解调方式解调实部与虚部,以第一帧实部的符号位作为输出,并储存其相位块的值;
S202、当接收到第二帧时,将第二帧与储存值相减取绝对值,得到相位块的差值;
S203、当差值在0、1、15或7、8、9时,说明相位差在0°和180°附近,则将新的相位块值储存,并根据其相位差输出111或000的软判决等级,以此循环;当差值在2、3、4、12、13、14或5、6、10、11时,由于偏移太大,不储存相位块值,并根据差值大小输出从110至100或001至011中的软判决等级。具体地,进一步地,所述步骤S203中根据差值输出的软判决等级具体标准为:当差值为0、1或15时,软判决等级为111;
当差值为2或14时,软判决等级为110;
当差值为3或13时,软判决等级为101;
当差值为4时,软判决等级为100;
当差值为12时,软判决等级为011;
当差值为5或11时,软判决等级为010;
当差值为6或10时,软判决等级为001;
当差值为7、8或9时,软判决等级为000;
其中,判决值1的软判决等级从111到100递减,判决值0的软判决等级从000到011递减。
与此同时,所述步骤S2还包括对子载波能量值的判断步骤,当实部和虚部都小于一定阈值时,软判决等级减小一级。
所述步骤S3为:循环对每帧每个子载波使用步骤S2进行软判决,得到所有子载波上每比特的软判决等级,并将得到的软判决等级作为后续维特比译码的输入进行译码。
综上所述,本发明提供的一种基于幅度特性的软判决方法,该方法具有适用范围广、复杂程度低及节约FPGA资源等优点。
以上所述仅为本发明的实施例,并非因此限制本发明的专利范围,凡是利用本发明说明书及附图内容所作的等同变换,或直接或间接运用在相关的技术领域,均同理包括在本发明的专利保护范围内。

Claims (8)

1.一种基于相位幅度特性的软判决方法,其特征在于:包括以下步骤:
S1、将星座图分割为16份,根据解调出的子载波数据的实部和虚部确定所述子载波数据在所述星座图中的相位块;
S2、通过比较前后两帧同一个子载波数据所处的相位块,得到对应所述子载波数据的软判决等级。
2.根据权利要求1所述的基于相位幅度特性的软判决方法,其特征在于:所述步骤S1中将星座图按以下方式分为十六份,划分规则为:
R1、最高位与星座点实部的符号位相同;
R2、第二位为:实部与虚部符号位的异或值;
R3、第三位为:若在1、3象限中实部绝对值大于虚部绝对值则为1,否则,为0;在2、4象限中相反;
R4、第四位为:根据星座点是否落在象限中间两块得到的值与第三位的异或,其中,若星座点落在象限中间两块,则值为1,否则,值为0。
3.根据权利要求1所述的基于相位幅度特性的软判决方法,其特征在于:所述步骤S1中采用非相干解调方式解调子载波数据的实部与虚部。
4.根据权利要求1所述的基于相位幅度特性的软判决方法,其特征在于:所述步骤S2的具体操作为:
S201、接收到第一帧时,以第一帧实部的符号位作为输出,并储存其相位块的值;
S202、当接收到第二帧时,将第二帧与储存值相减取绝对值,得到相位块的差值;
S203、当差值为0、1、15或7、8、9时,则根据差值输出不同的软判决等级并将新的相位块值储存;当差值为2、3、4、12、13、14或5、6、10、11时,不储存相位块值,并根据差值输出不同的软判决等级。
5.根据权利要求4所述的基于相位幅度特性的软判决方法,其特征在于:所述步骤S203中根据差值输出的软判决等级具体标准为:当差值为0、1或15时,软判决等级为111;
当差值为2或14时,软判决等级为110;
当差值为3或13时,软判决等级为101;
当差值为4时,软判决等级为100;
当差值为12时,软判决等级为011;
当差值为5或11时,软判决等级为010;
当差值为6或10时,软判决等级为001;
当差值为7、8或9时,软判决等级为000;
其中,判决值1的软判决等级从111到100递减,判决值0的软判决等级从000到011递减。
6.根据权利要求3所述的基于相位幅度特性的软判决方法,其特征在于:所述步骤S2还包括对子载波能量值的判断步骤,当实部和虚部都小于一阈值时,软判决等级减小一级。
7.根据权利要求1所述的基于相位幅度特性的软判决方法,其特征在于:所述方法还包括步骤:
S3、将所述步骤S2得到的软判决等级通过维特比译码后解码出实际数据。
8.根据权利要求7所述的基于相位幅度特性的软判决方法,其特征在于:所述步骤S3为循环对每帧每个子载波使用步骤S2进行软判决,得到所有子载波上每比特的软判决等级,并将得到的软判决等级作为后续维特比译码的输入进行译码。
CN201810612762.5A 2018-06-14 2018-06-14 一种基于相位幅度特性的软判决方法 Active CN108809882B (zh)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN201810612762.5A CN108809882B (zh) 2018-06-14 2018-06-14 一种基于相位幅度特性的软判决方法
CN202110255538.7A CN113055330B (zh) 2018-06-14 2018-06-14 一种低复杂度的基于相位幅度特性的软判决方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810612762.5A CN108809882B (zh) 2018-06-14 2018-06-14 一种基于相位幅度特性的软判决方法

Related Child Applications (1)

Application Number Title Priority Date Filing Date
CN202110255538.7A Division CN113055330B (zh) 2018-06-14 2018-06-14 一种低复杂度的基于相位幅度特性的软判决方法

Publications (2)

Publication Number Publication Date
CN108809882A true CN108809882A (zh) 2018-11-13
CN108809882B CN108809882B (zh) 2021-04-02

Family

ID=64085970

Family Applications (2)

Application Number Title Priority Date Filing Date
CN202110255538.7A Active CN113055330B (zh) 2018-06-14 2018-06-14 一种低复杂度的基于相位幅度特性的软判决方法
CN201810612762.5A Active CN108809882B (zh) 2018-06-14 2018-06-14 一种基于相位幅度特性的软判决方法

Family Applications Before (1)

Application Number Title Priority Date Filing Date
CN202110255538.7A Active CN113055330B (zh) 2018-06-14 2018-06-14 一种低复杂度的基于相位幅度特性的软判决方法

Country Status (1)

Country Link
CN (2) CN113055330B (zh)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1614964A (zh) * 2004-09-03 2005-05-11 清华大学 全球数字广播用基于差分幅度和相位联合调制的方法
CN101136898A (zh) * 2006-08-28 2008-03-05 大唐移动通信设备有限公司 一种正交调幅软判决方法及装置
CN103236906A (zh) * 2013-05-10 2013-08-07 湖南大学 Ami系统载波通信模块基带映射及解映射方法
US8718205B1 (en) * 2013-04-30 2014-05-06 Douglas Howard Morais Hard and soft bit demapping for QAM non-square constellations

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101562499B (zh) * 2009-05-27 2011-12-07 西安电子科技大学 16apsk编码调制中软判决信息获取方法
KR101252439B1 (ko) * 2011-08-26 2013-04-08 한양대학교 산학협력단 진폭 위상 변조 시스템을 위한 연판정 디매핑 장치 및 방법
CN102307175A (zh) * 2011-10-08 2012-01-04 四川虹微技术有限公司 多载波系统的软判决方法
CN102710581B (zh) * 2012-05-07 2014-09-03 清华大学 一种软解映射方法
CN104283835B (zh) * 2013-07-02 2018-05-18 大唐联诚信息系统技术有限公司 一种正交幅度调制软比特解调方法及其装置
KR102021314B1 (ko) * 2015-05-15 2019-09-16 한국전자통신연구원 비정방 직교진폭변조 방식의 연판정 복조 장치 및 방법

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1614964A (zh) * 2004-09-03 2005-05-11 清华大学 全球数字广播用基于差分幅度和相位联合调制的方法
CN101136898A (zh) * 2006-08-28 2008-03-05 大唐移动通信设备有限公司 一种正交调幅软判决方法及装置
US8718205B1 (en) * 2013-04-30 2014-05-06 Douglas Howard Morais Hard and soft bit demapping for QAM non-square constellations
CN103236906A (zh) * 2013-05-10 2013-08-07 湖南大学 Ami系统载波通信模块基带映射及解映射方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
CHAO XU 等: "Soft-Decision Multiple-Symbol Differential Sphere Detection and Decision-Feedback Differential Detection for Differential QAM Dispensing with Channel Estimation in the Face of Rapidly Fading Channels", 《IEEE TRANSACTIONS ON WIRELESS COMMUNICATIONS》 *
张力 等: "OFDM系统中改进的16QAM软判决解调算法", 《现代电子技术》 *

Also Published As

Publication number Publication date
CN113055330A (zh) 2021-06-29
CN113055330B (zh) 2023-10-20
CN108809882B (zh) 2021-04-02

Similar Documents

Publication Publication Date Title
JP6438045B2 (ja) カバレージ改善のための低papr変調
CN107948113B (zh) 基于三维信号插入降低ofdm系统峰均功率比的方法及系统
CN103124250B (zh) 接收设备、信号处理设备以及信号处理方法
JP2002330188A (ja) Qam信号のビット対数尤度比演算の方法と装置
CN102307175A (zh) 多载波系统的软判决方法
ES2362759B1 (es) Procedimiento y dispositivo de comunicaciones digitales para la recepción de datos usando simbolos qam.
US8718205B1 (en) Hard and soft bit demapping for QAM non-square constellations
CN101540752B (zh) Qam软解调方法
WO2009115795A2 (en) Method and apparatus for performing log- likelihood calculations
US8023578B2 (en) Two-stage low-complexity max-log bit-level LLR calculator and method
CN103209151A (zh) 通用星座解调方法及系统
US8811539B2 (en) Demapper for rotated QAM constellations
CN101404564A (zh) 8psk格雷映射的一种软解调方法
CN103905370B (zh) 正交振幅调制信号的软解调方法和装置、移动终端
CN108809882A (zh) 一种基于相位幅度特性的软判决方法
CN103152303B (zh) 一种信号差分解调方法和装置
CN105530215B (zh) 一种基于cordic的8psk软差分译码方法
CN102045293A (zh) 一种ofdm系统接收端星座图解调功率归一化方法
CN101938333A (zh) 格雷编码π/M-MPSK调制软比特信息计算方法
CN115208725B (zh) 一种联合ofdm同步与信息调制的方法、装置及介质
CN104468453A (zh) 一种适用于光正交频分复用传输系统的相位盲估计方法
CN102769600A (zh) 一种ofdm系统接收机中缩放因子解调制的方法及装置
Sridhar et al. Constellation design for dual-mode OFDM-IM
CN112968858B (zh) 一种qam软判决解映射的fpga实现方法
CN114389749A (zh) 对数似然比获取方法、装置、设备及存储介质

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