CN113708873B - 低轨卫星tdma静中通系统的返向链路多普勒频移估计方法 - Google Patents
低轨卫星tdma静中通系统的返向链路多普勒频移估计方法 Download PDFInfo
- Publication number
- CN113708873B CN113708873B CN202110973070.5A CN202110973070A CN113708873B CN 113708873 B CN113708873 B CN 113708873B CN 202110973070 A CN202110973070 A CN 202110973070A CN 113708873 B CN113708873 B CN 113708873B
- Authority
- CN
- China
- Prior art keywords
- satellite
- sat
- return
- hub
- time
- 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
- 238000000034 method Methods 0.000 title claims abstract description 25
- 230000003068 static effect Effects 0.000 title claims abstract description 11
- 230000005540 biological transmission Effects 0.000 claims abstract description 20
- 238000004891 communication Methods 0.000 claims description 10
- HOWHQWFXSLOJEF-MGZLOUMQSA-N systemin Chemical compound NCCCC[C@H](N)C(=O)N[C@@H](CCSC)C(=O)N[C@@H](CCC(N)=O)C(=O)N[C@@H]([C@@H](C)O)C(=O)N[C@@H](CC(O)=O)C(=O)OC(=O)[C@@H]1CCCN1C(=O)[C@H]1N(C(=O)[C@H](CC(O)=O)NC(=O)[C@H](CCCN=C(N)N)NC(=O)[C@H](CCCCN)NC(=O)[C@H](CO)NC(=O)[C@H]2N(CCC2)C(=O)[C@H]2N(CCC2)C(=O)[C@H](CCCCN)NC(=O)[C@H](CO)NC(=O)[C@H](CCC(N)=O)NC(=O)[C@@H](NC(=O)[C@H](C)N)C(C)C)CCC1 HOWHQWFXSLOJEF-MGZLOUMQSA-N 0.000 claims description 3
- 108010050014 systemin Proteins 0.000 claims description 3
- 230000033001 locomotion Effects 0.000 description 6
- 238000010586 diagram Methods 0.000 description 5
- 230000000694 effects Effects 0.000 description 2
- 238000013213 extrapolation Methods 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 1
- 230000001934 delay Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
Images
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04J—MULTIPLEX COMMUNICATION
- H04J3/00—Time-division multiplex systems
- H04J3/02—Details
- H04J3/06—Synchronising arrangements
- H04J3/0602—Systems characterised by the synchronising information used
- H04J3/0617—Systems characterised by the synchronising information used the synchronising signal being characterised by the frequency or phase
-
- 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
-
- Y—GENERAL 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02D—CLIMATE 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/00—Reducing energy consumption in communication networks
- Y02D30/70—Reducing energy consumption in communication networks in wireless communication networks
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Physics & Mathematics (AREA)
- Astronomy & Astrophysics (AREA)
- Aviation & Aerospace Engineering (AREA)
- General Physics & Mathematics (AREA)
- Radio Relay Systems (AREA)
Abstract
本申请提供一种透明转发低轨卫星TDMA静中通系统的返向链路多普勒频移估计方法,包括步骤:计算出卫星可视窗口内卫星/主站传输时延与星历时间的函数、卫星对主站的径向速度与星历时间的函数以及卫星对端站的径向速度与星历时间的函数;根据返向突发到达主站目标时隙的期望星历时间,估算出返向突发到达卫星的星历时间;根据返向突发到达卫星的星历时间以及卫星对主站的径向速度与星历时间的函数和卫星对端站的径向速度与星历时间的函数估算出返向突发到达卫星时,卫星对主站的径向速度和卫星对端站的径向速度;估算出返向上行链路和返向下行链路的多普勒频移,并相加,得到整个返向链路的多普勒频移。
Description
技术领域
本发明属于数字通信技术领域,涉及TDMA卫星通信系统的频率同步技术,尤其涉及一种透明转发低轨卫星TDMA静中通系统的返向链路多普勒频移估计方法。
背景技术
在无线通信中,当信源和信宿在径向方向上存在相对运动时,信宿接收信号的频率与信源发射信号的频率不一致,即信宿接收信号的频率发生了偏移,这种现象称为多普勒效应。收发两端信号频率的偏差称为多普勒频移。如图1所示,信源静止,信宿以速度v延箭头指向的方向移动,并在移动中与信源进行通信。信源发射的信号的频率是fc。当信宿在Δt时间内从位置A移动到位置A′时,其接收信号产生的多普勒频移是
其中,θ是位置A′处信宿移动方向与信号传播方向的夹角,c是电磁波的传播速度,v·cos(θ)表示位置A′处信宿对信源的径向速度。
如图2所示,在透明转发低轨卫星TDMA通信系统中,返向链路是一条从端站到卫星,再到主站的传输链路,由返向上行链路和返向下行链路两部分组成。其中,返向上行链路是指从端站到卫星的传输链路,返向下行链路是指从卫星到主站的传输链路。
对于静中通场景,主站和端站的位置固定不变,而卫星处于运动状态。这等价于在返向上行链路中,信源静止,而信宿移动;在返向下行链路中,信源移动,而信宿静止。因此,在返向上行链路和返向下行链路中均存在信源和信宿在径向方向上的相对运动,即在返向上行链路和返向下行链路中均存在多普勒频移。分别估计出返向上行链路和返向下行链路的多普勒频移,然后将二者相加,就可以得到整个返向链路的多普勒频移。
在低轨卫星通信系统中,由于卫星的位置是随时间变化的,返向链路不是一条固定不变的传输链路,而是随时间变化而变化的。准确估计返向链路多普勒频移的首要前提是正确找到待估计的返向链路是哪一条。
对于静中通场景,因为主站和端站的位置都是固定不变的,所以返向链路由卫星的位置唯一确定。而卫星的位置是时间的函数,因此,只要确定了返向突发到达卫星的时间,就能够确定相应时刻卫星的位置,进而唯一地确定传输该返向突发的返向链路。举例,如图2所示,返向突发BF[n]在t″n时刻从端站发出,在t′n时刻到达卫星,随后被卫星转发给主站。在t′n时刻,卫星的位置是P(t′n),因此,传输BF[n]的返向链路就是R→P(t′n)→H。综上所述,确定返向突发到达卫星的时间是准确估计返向链路多普勒频移的关键。
发明内容
针对上述现状,本发明提供一种透明转发低轨卫星TDMA静中通系统的返向链路多普勒频移估计方法,该方法采用了基于线性近似的外推估计法,计算简便;在卫星星历信息准确的条件下,可以获得良好的估计性能。
为了实现本发明的目的,本发明拟通过以下技术方案实现:
一种透明转发低轨卫星TDMA静中通系统的返向链路多普勒频移估计方法,包括以下步骤:
S2、在直角坐标系中,建立两个线性方程:一个是反恒等直线方程l1:τ′SAT/HUB(t)=-t,另一个是卫星/主站传输时延与星历时间的函数τSAT/HUB(t)=g1(t)的坐标平移版本τ′SAT/HUB(t)=τSAT/HUB(t+tn)=g1(t+tn)在区间[-T,0]上的局部近似直线段方程l2:T是一个时间增量,且 表示返向下行链路自由空间传输时延的最大值;
其中,直线段方程l2建立方法是:在曲线τ′SAT/HUB(t)=g1(t+tn),t∈[-T,0]上选取两个相邻的点(0,g1(tn))和(-T,g1(tμ)),其中,tμ=tn-T;然后,根据上述两点的坐标建立直线段方程,得
S3、求解直线l1和直线段l2的交点坐标tintersect:
S4、估算返向突发BF[n]到达卫星的星历时间,得
t′n=tn+tintersect;
S5、根据预先计算的两组参数:卫星对主站的径向速度与星历时间的函数vSAT/HUB(t)=h1(t)以及卫星对端站的径向速度与星历时间的函数vSAT/RCST(t)=h2(t),估算返向突发BF[n]到达卫星时,卫星对主站的径向速度vSAT/HUB(t′n)和卫星对端站的径向速度vSAT/RCST(t′n),得
vSAT/HUB(t′n)=h1(t′n)
vSAT/RCST(t′n)=h2(t′n);
S6、根据下式
本发明有益效果在于,提供一种透明转发低轨卫星TDMA静中通系统的返向链路多普勒频移估计方法。该方法采用了基于线性近似的外推估计法,计算简便;在卫星星历信息准确的条件下,可以获得良好的估计性能。
附图说明
图1是分析多普勒效应机理及计算多普勒频移的示意图。
图2是透明转发低轨卫星TDMA通信系统的返向链路示意图。
图3是本申请实施例中估算返向突发BF[n]到达卫星的星历时间t′n的方法的原理示意图。
图4是本申请实施例中步骤S1~S3的几何示意图。
具体实施方式
为了使本申请的目的、技术方案和具体实施方法更为清楚,结合附图实例对本申请进行进一步详细说明。
本申请实施例提供一种透明转发低轨卫星TDMA静中通系统的返向链路多普勒频移估计方法,设计思路如下:
首先,根据卫星的星历信息以及主站和端站的GNSS位置信息预先计算出卫星可视窗口内,卫星/主站传输时延与星历时间的函数τSAT/HUB(t)=g1(t)、卫星对主站的径向速度与星历时间的函数vSAT/HUB(t)=h1(t)以及卫星对端站的径向速度与星历时间的函数vSAT/RCST(t)=h2(t)。
然后,根据返向突发BF[n]到达主站目标时隙的期望星历时间tn(已知量)和函数关系τSAT/HUB(t)=g1(t),估算出返向突发BF[n]到达卫星的星历时间t′n。
然后,根据返向突发BF[n]到达卫星的星历时间t′n以及函数关系vSAT/HUB(t)=h1(t)和vSAT/RCST(t)=h2(t)估算出返向突发BF[n]到达卫星时,卫星对主站的径向速度vSAT/HUB(t′n)和卫星对端站的径向速度vSAT/RCST(t′n)。
然后,根据卫星对端站的径向速度和卫星对主站的径向速度分别估算出返向上行链路的多普勒频移和返向下行链路的多普勒频移。
最后,将返向上行链路和返向下行链路的多普勒频移相加,得到整个返向链路的多普勒频移。
为使用本申请实施例所述方法准备两个前提条件:
1、在主站侧,建立NCR时间与星历时间的映射关系。
2、端站完成前向链路时间同步,即NCR同步。
在上述前提条件下,为了简化返向链路多普勒频移估计的计算复杂度,需要预先准备三组参数:卫星/主站传输时延与星历时间的函数τSAT/HUB(t)=g1(t)、卫星对主站的径向速度与星历时间的函数vSAT/HUB(t)=h1(t)以及卫星对端站的径向速度与星历时间的函数vSAT/RCST(t)=h2(t)。
其中,计算卫星/主站传输时延与星历时间的函数的方法如下:
1)根据卫星的星历信息计算出卫星可视窗口内,卫星在地心地固(ECEF)坐标系下的位置与星历时间的函数。
2)根据主站的GNSS位置换算得到主站在ECEF坐标系下的位置。
3)在ECEF坐标系下,根据卫星和主站的位置信息,计算出卫星/主站距离与星历时间的函数dSAT/HUB(t)=f1(t)。
4)用dSAT/HUB(t)=f1(t)除以电磁波的传播速度c,得到卫星/主站传输时延与星历时间的函数τSAT/HUB(t)=dSAT/HUB(t)/c。
其中,计算卫星对主站的径向速度与星历时间的函数的方法如下:
1)根据卫星的星历信息计算出卫星可视窗口内,卫星在ECEF坐标系下的线速度与星历时间的函数v(t)=q(t)。
2)根据主站的GNSS位置换算得到主站在ECEF坐标系下的位置。
3)在ECEF坐标系下,根据卫星的星历信息和主站的位置信息,计算出卫星可视窗口内,卫星的线速度方向与卫星/主站连线的夹角θSAT/HUB(t)与星历时间的函数θSAT/HUB(t)=s1(t)。4)根据v(t)=q(t)和θSAT/HUB(t)=s1(t)计算出卫星对主站的径向速度与星历时间的函数,即vSAT/HUB(t)=v(t)·cos(θSAT/HUB(t))。
参照上述方法,可以计算出卫星对端站的径向速度与星历时间的函数。
在上述条件下,进行返向链路的多普勒频移估计:
S2、在直角坐标系中,建立两个线性方程:一个是反恒等直线方程l1:τ′SAT/HUB(t)=-t,另一个是卫星/主站传输时延与星历时间的函数τSAT/HUB(t)=g1(t)的坐标平移版本τ′SAT/HUB(t)=τSAT/HUB(t+tn)=g1(t+tn)在区间[-T,0]上的局部近似直线段方程l2:
其中,直线段方程l2建立方法是:
首先,在曲线τ′SAT/HUB(t)=g1(t+tn),t∈[-T,0]上选取两个相邻的点(0,g1(tn))和(-T,g1(tμ)),其中,tμ=tn-T,T是一个时间增量,且 表示返向下行链路自由空间传输时延的最大值;
然后,根据上述两点的坐标建立直线段方程,得
S3、求解直线l1和直线段l2的交点坐标tintersect:
S4、估算返向突发BF[n]到达卫星的星历时间,得
t′n=tn+tintersect。
具体的,估算返向突发BF[n]到达卫星的星历时间t′n的方法如下:如图2所示,在t″n时刻,返向突发BF[n]从端站发出,并在t′n时刻到达卫星。此时,卫星的位置是P(t′n)。随后,BF[n]被卫星转发给主站,并在tn时刻到达主站。而在tn-t′n的时间内,卫星也从位置P(t′n)移动到了新的位置P(tn)。卫星在轨道P(t′n)→P(tn)上的运动时间与BF[n]在返向下行链路P(t′n)→H上的传输时延恰好相等,即直线τSAT/HUB(t)=-t+tn与曲线τSAT/HUB(t)=g1(t)在tn之前必然相交,且交点就是t′n(如图3所示,其代表的场景是卫星与主站的距离逐渐减小,因此返向下行链路的传输时延也逐渐减小。在本图中,g1(t′n)=tn-t′n)。因此,求解非线性方程组
就可以得到返向突发BF[n]到达卫星的星历时间t′n。
通常,在低轨卫星通信系统中,返向下行链路的自由空间传输时延很小,卫星在这段时间的移动距离很短,其运动轨迹可以近似为直线段。因此,曲线τSAT/HUB(t)=g1(t)在tn附近的区域也可以近似为直线段。令此直线段的方程为其中,T是一个时间增量。进而,上述问题可以简化为直线τSAT/HUB(t)=-t+tn和直线段的相交问题。综上所述,求解线性方程组
就可以得到返向突发BF[n]到达卫星的星历时间t′n。
S5、估算返向突发BF[n]到达卫星时,卫星对主站的径向速度vSAT/HUB(t′n)和卫星对端站的径向速度vSAT/RCST(t′n),得
S6、根据下式
Claims (3)
1.一种透明转发低轨卫星TDMA静中通系统的返向链路多普勒频移估计方法,其特征在于,包括以下步骤:
S2、在直角坐标系中,建立两个线性方程:一个是反恒等直线方程l1:τ′SAT/HUB(t)=-t,另一个是卫星/主站传输时延与星历时间的函数τSAT/HUB(t)=g1(t)的坐标平移版本τ′SAT/HUB(t)=τSAT/HUB(t+tn)=g1(t+tn)在区间[-T,0]上的局部近似直线段方程l2:T是一个时间增量,且 表示返向下行链路自由空间传输时延的最大值;
S3、求解直线l1和直线段l2的交点坐标tintersect;
S4、估算返向突发BF[n]到达卫星的星历时间,得
t′n=tn+tintersect;
S5、根据预先计算的两组参数:卫星对主站的径向速度与星历时间的函数vSAT/HUB(t)=h1(t)以及卫星对端站的径向速度与星历时间的函数vSAT/RCST(t)=h2(t),估算返向突发BF[n]到达卫星时,卫星对主站的径向速度vSAT/HUB(t′n)和卫星对端站的径向速度vSAT/RCST(t′n),得
S6、根据下式
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110973070.5A CN113708873B (zh) | 2021-08-24 | 2021-08-24 | 低轨卫星tdma静中通系统的返向链路多普勒频移估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110973070.5A CN113708873B (zh) | 2021-08-24 | 2021-08-24 | 低轨卫星tdma静中通系统的返向链路多普勒频移估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113708873A CN113708873A (zh) | 2021-11-26 |
CN113708873B true CN113708873B (zh) | 2023-04-14 |
Family
ID=78654273
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110973070.5A Active CN113708873B (zh) | 2021-08-24 | 2021-08-24 | 低轨卫星tdma静中通系统的返向链路多普勒频移估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113708873B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114978829B (zh) * | 2022-05-12 | 2024-05-28 | 重庆金美通信有限责任公司 | 一种低轨卫星通信系统反向频偏预校正机制 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0901240A2 (en) * | 1997-09-06 | 1999-03-10 | ICO Services Ltd. | User terminal doppler compensation |
US6307503B1 (en) * | 1998-05-07 | 2001-10-23 | Hughes Electronics Corporation | Method and system of determining user terminal position using a measured signal propagation delay and Doppler shift of a communications link |
CN108027442A (zh) * | 2015-10-12 | 2018-05-11 | 迪尔公司 | 具有定点西格玛罗滤波器的卫星导航接收机 |
CN108259079A (zh) * | 2017-12-29 | 2018-07-06 | 中国电子科技集团公司第二十研究所 | 基于星历的高速移动平台tdma卫星通信同步控制方法 |
CN110324075A (zh) * | 2019-05-31 | 2019-10-11 | 西安空间无线电技术研究所 | 一种低轨星座上行信号时频校正周期确定方法 |
CN112019252A (zh) * | 2020-08-03 | 2020-12-01 | 航天科工空间工程发展有限公司 | 低轨卫星通信系统返向频偏校正方法和系统 |
-
2021
- 2021-08-24 CN CN202110973070.5A patent/CN113708873B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0901240A2 (en) * | 1997-09-06 | 1999-03-10 | ICO Services Ltd. | User terminal doppler compensation |
US6307503B1 (en) * | 1998-05-07 | 2001-10-23 | Hughes Electronics Corporation | Method and system of determining user terminal position using a measured signal propagation delay and Doppler shift of a communications link |
CN108027442A (zh) * | 2015-10-12 | 2018-05-11 | 迪尔公司 | 具有定点西格玛罗滤波器的卫星导航接收机 |
CN108259079A (zh) * | 2017-12-29 | 2018-07-06 | 中国电子科技集团公司第二十研究所 | 基于星历的高速移动平台tdma卫星通信同步控制方法 |
CN110324075A (zh) * | 2019-05-31 | 2019-10-11 | 西安空间无线电技术研究所 | 一种低轨星座上行信号时频校正周期确定方法 |
CN112019252A (zh) * | 2020-08-03 | 2020-12-01 | 航天科工空间工程发展有限公司 | 低轨卫星通信系统返向频偏校正方法和系统 |
Non-Patent Citations (4)
Title |
---|
Shervin Amiri ; Mohammad Mehdipour."Accurate Doppler Frequency Shift Estimation for any Satellite Orbit".《2007 3rd International Conference on Recent Advances in Space Technologies》.2007,全文. * |
张骞丹 ; 田红心 ; .GPS系统多普勒频移估算的研究.无线电工程.2007,(第04期),全文. * |
彭耿 ; 黄知涛 ; 李强 ; 姜文利 ; .中低轨卫星信号的多普勒频移估计与补偿.系统工程与电子技术.2009,(第02期),全文. * |
李杰."火星探测器行星际定轨研究".《中国博士学位论文全文数据库信息科技辑》.2011,全文. * |
Also Published As
Publication number | Publication date |
---|---|
CN113708873A (zh) | 2021-11-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110545136B (zh) | 低轨卫星多载波通信系统上行链路载波频偏估计与补偿方法 | |
CN107333241B (zh) | 基于lte体制的卫星移动通信上行发射端定时调整方法 | |
CN101933304B (zh) | 一种移动装置定位方法及设备 | |
JP4060368B2 (ja) | 低高度地球軌道衛星2機を用いた受動的位置決定 | |
EP2360967B1 (en) | A method for communicating in a wireless communication network | |
EP1110100B1 (en) | Accurate range and range rate determination in a satellite communications system | |
CN113708874B (zh) | 低轨卫星tdma动中通系统的返向链路多普勒频移估计方法 | |
CN110366095B (zh) | 一种水下节点分布式定位系统及方法 | |
KR20000062203A (ko) | 로컬 위치 결정 시스템 | |
CN113708873B (zh) | 低轨卫星tdma静中通系统的返向链路多普勒频移估计方法 | |
CN111448480B (zh) | 具有全球导航卫星系统信号生成装置和辐射电缆的定位系统 | |
US20210396832A1 (en) | Method and Apparatus for Determining the Angle of Departure | |
CN111262640A (zh) | 多普勒偏移的测试方法及装置、可读存储介质 | |
CN111856524B (zh) | 同频双向连续高精度测量方法及系统 | |
CN106054226B (zh) | 一种移动蜂窝网与卫星导航系统联合卫星信号盲区定位方法 | |
CN105487094A (zh) | 一种数据链与卫星导航协同定位方法及定位系统 | |
EP3879295A1 (en) | Cellular positioning based verification of an estimated vehicle position | |
KR20200041684A (ko) | 위치 측정 서비스를 제공하는 사용자 단말, 방법 및 컴퓨터 프로그램 | |
GB2332111A (en) | Compensating for motion in satellite telephone handset | |
CN113708872B (zh) | 低轨卫星tdma静中通系统的前向链路时延估计方法 | |
CN113708871B (zh) | 低轨卫星tdma静中通系统的返向链路时延估计方法 | |
CN113708875B (zh) | 低轨卫星tdma动中通系统的前向链路时延估计方法 | |
CN113708870B (zh) | 低轨卫星tdma动中通系统的返向突发发送时间估计方法 | |
KR20190129578A (ko) | 고정형 앵커 기반 멀티 홉 시각 동기 장치 및 이를 이용한 owr 측위 시스템 | |
US7965237B2 (en) | Mobile system and method for position estimation |
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 |