CN112468215A - 适应低轨卫星高动态环境的下行频率盲补偿方法及系统 - Google Patents

适应低轨卫星高动态环境的下行频率盲补偿方法及系统 Download PDF

Info

Publication number
CN112468215A
CN112468215A CN202011281662.2A CN202011281662A CN112468215A CN 112468215 A CN112468215 A CN 112468215A CN 202011281662 A CN202011281662 A CN 202011281662A CN 112468215 A CN112468215 A CN 112468215A
Authority
CN
China
Prior art keywords
satellite
downlink frequency
blind compensation
orbit
type
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
CN202011281662.2A
Other languages
English (en)
Other versions
CN112468215B (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.)
Xian Institute of Space Radio Technology
Original Assignee
Xian Institute of Space Radio Technology
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 Xian Institute of Space Radio Technology filed Critical Xian Institute of Space Radio Technology
Priority to CN202011281662.2A priority Critical patent/CN112468215B/zh
Publication of CN112468215A publication Critical patent/CN112468215A/zh
Application granted granted Critical
Publication of CN112468215B publication Critical patent/CN112468215B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04B—TRANSMISSION
    • H04B7/00—Radio transmission systems, i.e. using radiation field
    • H04B7/14—Relay systems
    • H04B7/15—Active relay systems
    • H04B7/185—Space-based or airborne stations; Stations for satellite systems
    • H04B7/1851—Systems using a satellite or space-based relay
    • H04B7/18519—Operations control, administration or maintenance
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04B—TRANSMISSION
    • H04B7/00—Radio transmission systems, i.e. using radiation field
    • H04B7/14—Relay systems
    • H04B7/15—Active relay systems
    • H04B7/185—Space-based or airborne stations; Stations for satellite systems
    • H04B7/1851—Systems using a satellite or space-based relay
    • H04B7/18513—Transmission in a satellite or space-based system

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Astronomy & Astrophysics (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Radio Relay Systems (AREA)

Abstract

本发明公开了一种适应低轨卫星高动态环境的下行频率盲补偿方法及系统，其中，该方法如下步骤：(1)基于星历信息中的平均移动速度，计算出卫星的轨道半径；(2)依据波束空间位置分布和步骤(1)中卫星的轨道半径，计算各类型波束几何中心点与卫星的距离，同时得到各类型波束几何中心与星地连线的夹角；(3)计算各类型波束的下行频率盲补偿频率；(4)在星上处理时，对每个波束的均增加下行频率盲补偿频率。本发明提高了地面终端与卫星进行频率同步的速度，降低了地面终端的处理复杂度，提升了低轨卫星星座系统的接入效率。

Description

适应低轨卫星高动态环境的下行频率盲补偿方法及系统
技术领域
本发明属于卫星通信技术领域，尤其涉及一种适应低轨卫星高动态环境的下行频率盲补偿方法及系统。
背景技术
由于LEO卫星的高动态特性，因此地面终端通过下行信道发现并同步卫星信号是星地通信建立的首要条件，是地面终端开机时首先要完成的任务。同样还是由于LEO卫星的高动态特性，用户链路下行信道的信号到达地面终端时，会产生较大的多普勒频移，使得地面终端首次同步卫星信号的过程变得漫长且复杂，极大的影响了用户体验。
现有的技术针对低轨移动卫星星座系统，已有成果中的下行多普勒频率补偿算法主要是基于接收端的后处理，而不是在发射端进行预补偿或盲补偿。
发明内容
本发明解决的技术问题是：克服现有技术的不足，提供了一种适应低轨卫星高动态环境的下行频率盲补偿方法及系统，提高了地面终端与卫星进行频率同步的速度，降低了地面终端的处理复杂度，提升了低轨卫星星座系统的接入效率。
本发明目的通过以下技术方案予以实现：一种适应低轨卫星高动态环境的下行频率盲补偿方法，所述方法包括如下步骤：(1)基于星历信息中的平均移动速度，计算出卫星的轨道半径；(2)依据波束空间位置分布和步骤(1)中卫星的轨道半径，计算各类型波束几何中心点与卫星的距离，同时得到各类型波束几何中心与星地连线的夹角；(3)计算各类型波束的下行频率盲补偿频率；(4)在星上处理时，对每个波束的均增加下行频率盲补偿频率。
上述适应低轨卫星高动态环境的下行频率盲补偿方法中，在步骤(1)中，卫星的轨道半径为：
Figure BDA0002781009550000021
其中，R为卫星轨道半径，μ＝3.986005×1014m3/s2为以地球为中心的引力常数，n0星历信息中的平均移动速度。
上述适应低轨卫星高动态环境的下行频率盲补偿方法中，在步骤(2)中，波束类型的划分是依据卫星波束排布图中每个波束几何中心点到星下点之间距离，距离相同的波束划分为同一类波束。
上述适应低轨卫星高动态环境的下行频率盲补偿方法中，波束类型共3类，将这3类波束定义为中心波束、内圈波束、外圈波束。
上述适应低轨卫星高动态环境的下行频率盲补偿方法中，在步骤(3)中，各类型波束的下行频率盲补偿频率为：
Figure BDA0002781009550000022
其中，G＝6.67259×10-11Nm2/kg2为万有引力常数，M＝5.965×1024kg为地球质量，R为卫星的轨道半径，c＝3×108m/s为光速，f为通信频率，fi为各类型波束的下行频率盲补偿频率，i为波束的类型。
上述适应低轨卫星高动态环境的下行频率盲补偿方法中，中心波束几何中心距离卫星1114.5km，内圈波束几何中心距离卫星1491.3km，外圈波束结合中心距离卫星2231.7km。
上述适应低轨卫星高动态环境的下行频率盲补偿方法中，卫星的轨道半径为7448km。
一种适应低轨卫星高动态环境的下行频率盲补偿系统，包括：第一模块，用于基于星历信息中的平均移动速度，计算出卫星的轨道半径；第二模块，用于依据波束空间位置分布和卫星的轨道半径，计算各类型波束几何中心点与卫星的距离，同时得到各类型波束几何中心与星地连线的夹角；第三模块，用于计算各类型波束的下行频率盲补偿频率；第四模块，用于在星上处理时，对每个波束的均增加下行频率盲补偿频率。
上述适应低轨卫星高动态环境的下行频率盲补偿系统中，卫星的轨道半径为：
Figure BDA0002781009550000031
其中，R为卫星轨道半径，μ＝3.986005×1014m3/s2为以地球为中心的引力常数，n0星历信息中的平均移动速度。
上述适应低轨卫星高动态环境的下行频率盲补偿系统中，波束类型的划分是依据卫星波束排布图中每个波束几何中心点到星下点之间距离，距离相同的波束划分为同一类波束。
本发明与现有技术相比具有如下有益效果：
本发明基于每个类型的波束的空间位置分布信息、星历信息中的平均移动速度信息，分析星上每个类型的波束应该盲补偿的频率大小，然后在星上处理中进行补偿的方法，该方法可以有效的提高地面终端与卫星进行频率同步的速度，降低地面终端的处理复杂度，提升低轨卫星星座系统的接入效率。
附图说明
通过阅读下文优选实施方式的详细描述，各种其他的优点和益处对于本领域普通技术人员将变得清楚明了。附图仅用于示出优选实施方式的目的，而并不认为是对本发明的限制。而且在整个附图中，用相同的参考符号表示相同的部件。在附图中：
图1为本发明实施例下行多波束低轨卫星星座系统的典型链路的示意图；
图2为本发明实施例卫星波束排布的示意图；
图3为本发明实施例波束几何中心与卫星的距离的示意图；
图4为本发明实施例卫星处于不同位置时波束几何中心与星地连线夹角示意图；
图5为本发明实施例下行频率盲补偿处理框图。
具体实施方式
下面将参照附图更详细地描述本公开的示例性实施例。虽然附图中显示了本公开的示例性实施例，然而应当理解，可以以各种形式实现本公开而不应被这里阐述的实施例所限制。相反，提供这些实施例是为了能够更透彻地理解本公开，并且能够将本公开的范围完整的传达给本领域的技术人员。需要说明的是，在不冲突的情况下，本发明中的实施例及实施例中的特征可以相互组合。下面将参考附图并结合实施例来详细说明本发明。
具有下行多波束的低轨卫星星座系统的典型链路如图1所示，采用星上处理转发的方式，其中用户上下行链路均为多波束信号，是一种点对多点通信。该系统的一般特点是没有下行信标信号，也没有全覆盖的广播信号，仅有每个波束内的广播信号，因此简化了卫星星上的处理，但相应的增加了星地频率同步的难度，使地面终端变得复杂。
为了进一步提高地面终端与卫星进行频率同步的速度，降低地面终端的处理复杂度，提升低轨卫星星座系统的接入效率，本专利中设计的下行频率盲补偿方法，包括四个步骤：
(1)基于星历信息中的平均移动速度，计算出卫星的轨道半径；
(2)依据波束空间位置分布，和步骤(1)中计算出的轨道半径，计算不同类型的波束几何中心点与卫星的距离，同时可以得到各类型波束几何中心与星地连线的夹角。所述波束类型的划分是依据卫星波束排布图中每个波束几何中心点到星下点之间距离，距离相同的波束划分为同一类波束。
(3)计算各类波束的下行频率盲补偿大小。
(4)根据频率大小，在星上处理时，对每个波束的处理中均增加频偏盲补偿模块。
以下内容中，各个参数都给出了具体数值作为示例。
(1)轨道半径计算
星历信息中的平均移动速度：以n0＝13.5转/天为例，则依据公式1计算卫星轨道半径：
Figure BDA0002781009550000051
其中，R为卫星轨道半径，μ＝3.986005×1014m3/s2为以地球为中心的引力常数，n0星历信息中的平均移动速度。可以计算得R＝7448km，则卫星轨道据地面高度为1070km。
(2)波束几何中心点与卫星的距离计算
以单星52个波束为例，如图2所示，分为3圈，因此该例中波束类型共3类，现将这3类波束定义为中心波束、内圈波束、外圈波束；
卫星的中心波束、内圈波束、外圈波束的几何中心与卫星的距离，可用简单的几何方法计算出来，这里直接给出结果，如图3所示，中心波束几何中心距离卫星1114.5km，内圈波束几何中心距离卫星1491.3km，外圈波束结合中心距离卫星2231.7km；
各类型波束几何中心与星地连线的夹角也可用简单的几何方法计算，为：θ1＝15°、θ2＝40°、θ3＝53.5°，如图4所示。
(3)下行频率盲补偿大小计算
各个波束下行频率盲补偿大小可以按照公式2计算：
Figure BDA0002781009550000052
其中，G＝6.67259×10-11Nm2/kg2为万有引力常数，M＝5.965×1024kg为地球质量，R为卫星轨道半径，c＝3×108m/s为光速，f为通信频率，这里取L频段频率，i＝1,2,3。
则可计算得到：f1＝10.56kHz、f2＝26.24kHz、f3＝32.81kHz，
盲补偿后地面终端需处理的频偏对比如下表：
表1下行频率盲补偿前后的对比
Figure BDA0002781009550000053
Figure BDA0002781009550000061
(4)星上盲补偿处理
在按照以上方法获得了所有波束的下行频率盲补偿值后，将其实现于星上处理器中，进行补偿频率生成，就可以实现下行频率盲补偿处理，处理过程如图5所示，原始波束信号经过与产生的频率生成信号进行相乘，得到补偿后的信号。
本实施例还提供了一种适应低轨卫星高动态环境的下行频率盲补偿系统，包括：第一模块，用于基于星历信息中的平均移动速度，计算出卫星的轨道半径；第二模块，用于依据波束空间位置分布和卫星的轨道半径，计算各类型波束几何中心点与卫星的距离，同时得到各类型波束几何中心与星地连线的夹角；第三模块，用于计算各类型波束的下行频率盲补偿频率；第四模块，用于在星上处理时，对每个波束的均增加下行频率盲补偿频率。
上述实施中，卫星的轨道半径为：
Figure BDA0002781009550000062
其中，R为卫星轨道半径，μ＝3.986005×1014m3/s2为以地球为中心的引力常数，n0星历信息中的平均移动速度。
上述实施中，波束类型的划分是依据卫星波束排布图中每个波束几何中心点到星下点之间距离，距离相同的波束划分为同一类波束。
本发明基于每个类型的波束的空间位置分布信息、星历信息中的平均移动速度信息，分析星上每个类型的波束应该盲补偿的频率大小，然后在星上处理中进行补偿的方法，该方法可以有效的提高地面终端与卫星进行频率同步的速度，降低地面终端的处理复杂度，提升低轨卫星星座系统的接入效率。
本发明虽然已以较佳实施例公开如上，但其并不是用来限定本发明，任何本领域技术人员在不脱离本发明的精神和范围内，都可以利用上述揭示的方法和技术内容对本发明技术方案做出可能的变动和修改，因此，凡是未脱离本发明技术方案的内容，依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化及修饰，均属于本发明技术方案的保护范围。

Claims (10)

1.一种适应低轨卫星高动态环境的下行频率盲补偿方法，其特征在于，所述方法包括如下步骤：
(1)基于星历信息中的平均移动速度，计算出卫星的轨道半径；
(2)依据波束空间位置分布和步骤(1)中卫星的轨道半径，计算各类型波束几何中心点与卫星的距离，同时得到各类型波束几何中心与星地连线的夹角；
(3)计算各类型波束的下行频率盲补偿频率；
(4)在星上处理时，对每个波束的均增加下行频率盲补偿频率。
2.根据权利要求1所述的适应低轨卫星高动态环境的下行频率盲补偿方法，其特征在于：在步骤(1)中，卫星的轨道半径为：
Figure FDA0002781009540000011
其中，R为卫星轨道半径，μ＝3.986005×1014m3/s2为以地球为中心的引力常数，n0星历信息中的平均移动速度。
3.根据权利要求1所述的适应低轨卫星高动态环境的下行频率盲补偿方法，其特征在于：在步骤(2)中，波束类型的划分是依据卫星波束排布图中每个波束几何中心点到星下点之间距离，距离相同的波束划分为同一类波束。
4.根据权利要求3所述的适应低轨卫星高动态环境的下行频率盲补偿方法，其特征在于：波束类型共3类，将这3类波束定义为中心波束、内圈波束、外圈波束。
5.根据权利要求1所述的适应低轨卫星高动态环境的下行频率盲补偿方法，其特征在于：在步骤(3)中，各类型波束的下行频率盲补偿频率为：
Figure FDA0002781009540000012
其中，G＝6.67259×10-11Nm2/kg2为万有引力常数，M＝5.965×1024kg为地球质量，R为卫星的轨道半径，c＝3×108m/s为光速，f为通信频率，fi为各类型波束的下行频率盲补偿频率，i为波束的类型。
6.根据权利要求4所述的适应低轨卫星高动态环境的下行频率盲补偿方法，其特征在于：中心波束几何中心距离卫星1114.5km，内圈波束几何中心距离卫星1491.3km，外圈波束结合中心距离卫星2231.7km。
7.根据权利要求2所述的适应低轨卫星高动态环境的下行频率盲补偿方法，其特征在于：卫星的轨道半径为7448km。
8.一种适应低轨卫星高动态环境的下行频率盲补偿系统，其特征在于包括：
第一模块，用于基于星历信息中的平均移动速度，计算出卫星的轨道半径；
第二模块，用于依据波束空间位置分布和卫星的轨道半径，计算各类型波束几何中心点与卫星的距离，同时得到各类型波束几何中心与星地连线的夹角；
第三模块，用于计算各类型波束的下行频率盲补偿频率；
第四模块，用于在星上处理时，对每个波束的均增加下行频率盲补偿频率。
9.根据权利要求8所述的适应低轨卫星高动态环境的下行频率盲补偿系统，其特征在于：卫星的轨道半径为：
Figure FDA0002781009540000021
其中，R为卫星轨道半径，μ＝3.986005×1014m3/s2为以地球为中心的引力常数，n0星历信息中的平均移动速度。
10.根据权利要求8所述的适应低轨卫星高动态环境的下行频率盲补偿系统，其特征在于：波束类型的划分是依据卫星波束排布图中每个波束几何中心点到星下点之间距离，距离相同的波束划分为同一类波束。
CN202011281662.2A 2020-11-16 2020-11-16 适应低轨卫星高动态环境的下行频率盲补偿方法及系统 Active CN112468215B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011281662.2A CN112468215B (zh) 2020-11-16 2020-11-16 适应低轨卫星高动态环境的下行频率盲补偿方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011281662.2A CN112468215B (zh) 2020-11-16 2020-11-16 适应低轨卫星高动态环境的下行频率盲补偿方法及系统

Publications (2)

Publication Number Publication Date
CN112468215A true CN112468215A (zh) 2021-03-09
CN112468215B CN112468215B (zh) 2022-10-14

Family

ID=74837856

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011281662.2A Active CN112468215B (zh) 2020-11-16 2020-11-16 适应低轨卫星高动态环境的下行频率盲补偿方法及系统

Country Status (1)

Country Link
CN (1) CN112468215B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114189419A (zh) * 2021-11-29 2022-03-15 西安空间无线电技术研究所 一种低轨卫星ofdm体制下行同步方法及系统
CN116015411A (zh) * 2022-12-16 2023-04-25 中电防务科技有限公司 一种卫星移动通信系统高动态终端及切换方法

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CA2166366A1 (fr) * 1995-01-02 1996-07-03 Luigi Bella Procede et systeme de transmission de signaux radioelectriques via un reseau de satellites entre une station terrestre fixe et des terminaux mobiles d'usagers
CN102158457A (zh) * 2011-03-29 2011-08-17 西安电子科技大学 多点多用户ofdma系统的频偏补偿方法
CN102798845A (zh) * 2012-07-20 2012-11-28 西安空间无线电技术研究所 一种微波散射计多普勒频率补偿方法
CN103220773A (zh) * 2012-01-20 2013-07-24 电信科学技术研究院 卫星通信系统中实现初始同步的方法及装置
CN103630887A (zh) * 2013-11-11 2014-03-12 南京航空航天大学 一种多波束形成方法及使用该方法的多波束声纳
CN103975539A (zh) * 2011-12-06 2014-08-06 爱立信(中国)通信有限公司 多普勒频移补偿装置和方法
CN110545138A (zh) * 2019-09-29 2019-12-06 中兴通讯股份有限公司 一种信息指示方法、装置及计算机可读存储介质
CN110798422A (zh) * 2019-11-25 2020-02-14 成都天奥集团有限公司 低轨卫星多载波通信系统下行链路采样频偏估计与补偿方法
CN111385013A (zh) * 2018-12-29 2020-07-07 华为技术有限公司 广播数据的方法和装置

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CA2166366A1 (fr) * 1995-01-02 1996-07-03 Luigi Bella Procede et systeme de transmission de signaux radioelectriques via un reseau de satellites entre une station terrestre fixe et des terminaux mobiles d'usagers
CN102158457A (zh) * 2011-03-29 2011-08-17 西安电子科技大学 多点多用户ofdma系统的频偏补偿方法
CN103975539A (zh) * 2011-12-06 2014-08-06 爱立信(中国)通信有限公司 多普勒频移补偿装置和方法
CN103220773A (zh) * 2012-01-20 2013-07-24 电信科学技术研究院 卫星通信系统中实现初始同步的方法及装置
CN102798845A (zh) * 2012-07-20 2012-11-28 西安空间无线电技术研究所 一种微波散射计多普勒频率补偿方法
CN103630887A (zh) * 2013-11-11 2014-03-12 南京航空航天大学 一种多波束形成方法及使用该方法的多波束声纳
CN111385013A (zh) * 2018-12-29 2020-07-07 华为技术有限公司 广播数据的方法和装置
CN110545138A (zh) * 2019-09-29 2019-12-06 中兴通讯股份有限公司 一种信息指示方法、装置及计算机可读存储介质
CN110798422A (zh) * 2019-11-25 2020-02-14 成都天奥集团有限公司 低轨卫星多载波通信系统下行链路采样频偏估计与补偿方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
KEITH A. WEAR: "Directivity and Frequency-Dependent Effective Sensitive Element Size of Membrane Hydrophones: Theory Versus Experiment", 《IEEE TRANSACTIONS ON ULTRASONICS, FERROELECTRICS, AND FREQUENCY CONTROL ( VOLUME: 66, ISSUE: 11, NOVEMBER 2019)》 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114189419A (zh) * 2021-11-29 2022-03-15 西安空间无线电技术研究所 一种低轨卫星ofdm体制下行同步方法及系统
CN114189419B (zh) * 2021-11-29 2023-08-29 西安空间无线电技术研究所 一种低轨卫星ofdm体制下行同步方法及系统
CN116015411A (zh) * 2022-12-16 2023-04-25 中电防务科技有限公司 一种卫星移动通信系统高动态终端及切换方法

Also Published As

Publication number Publication date
CN112468215B (zh) 2022-10-14

Similar Documents

Publication Publication Date Title
CN112468215B (zh) 适应低轨卫星高动态环境的下行频率盲补偿方法及系统
US5906337A (en) Multiple altitude satellite relay system and method
US5971324A (en) Multiple altitude satellite relay system and method
CN109120561B (zh) 低轨卫星移动通信系统中频偏校正方法
CN106027138B (zh) 规避与同步卫星共线干扰的地面站系统及方法
US9461733B2 (en) Device and method for optimizing the ground coverage of a hybrid space system
CN107333241A (zh) 基于lte体制的卫星移动通信上行发射端定时调整方法
JP2001506465A (ja) 高緯度静止衛星システム
US10763967B2 (en) Communications relay satellite with a single-axis gimbal
US6990314B1 (en) Multi-node point-to-point satellite communication system employing multiple geo satellites
CN111650613A (zh) 一种分布式星历计算方法
CA2187833A1 (en) System and method for a scalable and reliable transmission of electronic software distribution
US20150326304A1 (en) World-wide, wide-band, low-latency, mobile internet and system therefor
EP1208660A2 (en) Multi-node point-to-point satellite communication system employing multiple geo satellites
JP2002323552A (ja) 測位システム及び測位装置
CN111817769A (zh) 一种星地数据链传播时延补偿的时隙设计方法
US11863289B2 (en) Satellite communications system with non-geosynchronous orbits
EP1076942B1 (en) Signal acquisition method and apparatus in a satellite communication system
CN117856876B (zh) 一种高低轨星间分布式协同通信系统及方法
CN220650885U (zh) 天空地一体化无线导航系统
CN114726470A (zh) 一种非地面网络的终端同步系统及方法
Cui et al. A Supplementary Antenna Alignment Method Based on Terahertz Communication between LEO-LEO
CN115604851A (zh) 一种基于固定截距的大规模卫星星座频率干扰规避方法
JPH0496528A (ja) 電波中継空中基地
US20200162151A1 (en) Constellation design for martian synchronous orbit

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