CN109725648B - 一种相对导航卫星伴飞机动窗口计算方法、装置及介质 - Google Patents
一种相对导航卫星伴飞机动窗口计算方法、装置及介质 Download PDFInfo
- Publication number
- CN109725648B CN109725648B CN201811495677.1A CN201811495677A CN109725648B CN 109725648 B CN109725648 B CN 109725648B CN 201811495677 A CN201811495677 A CN 201811495677A CN 109725648 B CN109725648 B CN 109725648B
- Authority
- CN
- China
- Prior art keywords
- moment
- satellite
- orbit
- main satellite
- main
- 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
Landscapes
- Navigation (AREA)
Abstract
一种相对导航卫星伴飞机动窗口计算方法,首先根据主星和伴星的位置和速度参数,计算给定主星点火时刻和轨道转移时间情况下的变轨速度增量和与伴星交会时的相对速度;其次,给定主星点火时刻范围和轨道转移时间范围,以点火时刻为横坐标、轨道转移时间为纵坐标,获取主星加速的速度增量和交会时相对伴星的速度等高线图;最后,根据速度等高线图,得到满足主星变轨速度增量约束的点火时刻和轨道转移时间。此外,根据速度等高线图,还可得到主星变轨速度增量最优和轨道转移时间最优的机动窗口。该发明还可对给定的点火时刻和轨道转移时间,得到主星点火的方向和速度增量需求,以及交会时相对伴星的速度大小和方向。
Description
技术领域
本发明涉及一种相对导航卫星伴飞机动窗口计算方法、装置及介质,属于卫星控制技术领域。
背景技术
空间已成为维护国家安全和利益至关重要的战略制高点,开展空间相对导航任务研究,是我国空间技术向智能化、自主化迈出的重要一步。相对导航任务轨道机动涉及主星和伴星两个目标的相互运动,由主星择机点火进行轨道机动,在规定时间实现与伴星的伴飞任务。
由于相对导航任务双星均在空间环绕地球运行,满足轨道转移条件(变轨速度增量、点火方向、轨道转移行时间等约束)的窗口并非时刻存在,需对轨道机动分布、速度增量分布情况进行研究分析,为双星伴飞构型构建提供基础。传统方法主要是基于点火诸元解算,对给定的点火时刻和速度增量,计算点火方向,实现对伴星的交会,但无法对伴飞任务轨道机动窗口和变轨速度的整体分布进行分析。
发明内容
本发明要解决的技术问题是:克服现有技术的不足,提供了一种相对导航卫星伴飞机动窗口计算方法、装置及介质,在给定主星轨道机动时刻范围和轨道转移时间范围,即可根据主星和伴星轨道参数,获取主星变轨速度和交会时相对伴星的速度等高线图,解决了主星轨道机动窗口和变轨速度增量整体分布的问题,为空间相对导航任务伴飞构型的构建提供依据。
本发明目的通过以下技术方案予以实现:
一种相对导航卫星伴飞机动窗口计算方法,包括如下步骤:
步骤一、确定主星点火时刻TL(i)和轨道转移时间TOF(j),其中第一序数i=1且i≤N1,第二序数j=1;
步骤二、根据步骤一中所述的主星点火时刻TL(i)和轨道转移时间TOF(j),获得主星变轨时刻TL(i,j)和交会时刻TF(i,j);
步骤三、根据T0时刻主星轨道参数和步骤二所述TL(i,j)时刻,获得TL(i,j)时刻主星的轨道参数;根据T1时刻伴星轨道参数和步骤二所述TF(i,j)时刻获得TF(i,j)时刻伴星的轨道参数;
步骤四、根据步骤三中所述TL(i,j)时刻主星的轨道参数和TF(i,j)时刻伴星的轨道参数,采用求解兰伯特问题的方法,获得TL(i,j)时刻和TF(i,j)时刻的主星变轨速度增量大小和方向;根据TL(i,j)时刻和TF(i,j)时刻的主星变轨速度增量大小和方向,获得TL(i,j)时刻和TF(i,j)时刻的主星速度大小和方向;
步骤五、判断第二序数j≤N2是否成立,若成立,则令j的值增加1,然后转入步骤二;否则转入步骤六;
步骤六、判断第一序数i≤N1是否成立,若成立,则令i的值增加1且j等于1,然后转入步骤二;否则转入步骤七;
步骤七、基于步骤四获得的所有时刻主星变轨速度增量大小和主星速度大小,即获得相对导航卫星伴飞机动窗口。
上述相对导航卫星伴飞机动窗口计算方法,所述第一序数i取值范围的最大值N1为点火时刻分布区间内的采样点数量。
上述相对导航卫星伴飞机动窗口计算方法,所述第二序数j取值范围的最大值N2为主星飞行时间区间内的采样点数量。
一种相对导航卫星伴飞机动窗口计算装置,包括参数输入模块、数据计算模块、循环迭代模块;
所述参数输入模块用于接收外部参数,所述外部参数包括T0时刻主星轨道参数、T1时刻伴星轨道参数、主星点火时刻TL(i)、轨道转移时间TOF(j),然后参数输入模块将外部参数输出给所述循环迭代模块;
所述循环迭代模块用于将接收到的外部参数输出给所述数据计算模块,其中,循环迭代模块将不同的主星点火时刻TL(i)和不同的轨道转移时间TOF(j)进行排列组合;
所述数据计算模块根据所述的主星点火时刻TL(i)和轨道转移时间TOF(j),获得主星变轨时刻TL(i,j)和交会时刻TF(i,j);然后根据T0时刻主星轨道参数和TL(i,j)时刻,获得TL(i,j)时刻主星的轨道参数,根据T1时刻伴星轨道参数和TF(i,j)时刻获得TF(i,j)时刻伴星的轨道参数;最后根据所述TL(i,j)时刻主星的轨道参数和TF(i,j)时刻伴星的轨道参数,采用求解兰伯特问题的方法,获得TL(i,j)时刻和TF(i,j)时刻的主星变轨速度增量大小和方向,根据TL(i,j)时刻和TF(i,j)时刻的主星变轨速度增量大小和方向,获得TL(i,j)时刻和TF(i,j)时刻的主星速度大小和方向。
上述相对导航卫星伴飞机动窗口计算装置,还包括显示模块,所述显示模块用于将所述数据计算模块获得的所有时刻主星变轨速度增量大小和主星速度大小进行显示。
上述相对导航卫星伴飞机动窗口计算装置,所述第一序数i取值范围的最大值N1为点火时刻分布区间内的采样点数量。
上述相对导航卫星伴飞机动窗口计算装置,所述第二序数j取值范围的最大值N2为主星飞行时间区间内的采样点数量。
一种计算机存储介质,其上存储有计算机程序,该程序被处理器执行时实现上述相对导航卫星伴飞机动窗口计算方法的步骤。
本发明相比于现有技术具有如下有益效果:
(1)本发明方法与现有技术相比,通过绘制速度等高线图,可从图中直观获取变轨窗口,得到主星变轨最佳时机;
(2)本发明方法与现有技术相比,可以根据主星变轨时刻、轨道转移时间即交会时刻,得到主星变轨速度方向和交会时刻相对伴星的速度需求;
(3)本发明方法与现有技术相比,可以根据变轨时刻、轨道转移时间和变轨速度约束,得到主星的变轨方向;
(4)本发明方法面向相对导航领域,对主星变轨窗口、变轨速度增量、交会时刻相对伴星速度、轨道转移时间的整体分布清晰明了,对相对导航任务具有极强参考价值。
附图说明
图1为本发明方法的步骤流程图;
图2为本发明实施例的步骤流程图;
图3为本发明实施例的主星变轨速度等高线示意图;
图4为本发明实施例的交会时刻主星相对伴飞的速度等高线示意图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,下面将结合附图对本发明的实施方式作进一步详细描述。
一种相对导航卫星伴飞机动窗口计算方法,如图1所示,包括如下步骤:
步骤101、确定主星点火时刻TL(i)和轨道转移时间TOF(j),其中第一序数i=1且i≤N1,第二序数j=1;第一序数i取值范围的最大值N1为点火时刻分布区间内的采样点数量;第二序数j取值范围的最大值N2为主星飞行时间区间内的采样点数量;
步骤102、根据步骤101中所述的主星点火时刻TL(i)和轨道转移时间TOF(j),获得主星变轨时刻TL(i,j)和交会时刻TF(i,j);
步骤103、根据T0时刻主星轨道参数和步骤102所述TL(i,j)时刻,获得TL(i,j)时刻主星的轨道参数;根据T1时刻伴星轨道参数和步骤102所述TF(i,j)时刻获得TF(i,j)时刻伴星的轨道参数;
步骤104、根据步骤103中所述TL(i,j)时刻主星的轨道参数和TF(i,j)时刻伴星的轨道参数,采用求解兰伯特问题的方法,获得TL(i,j)时刻和TF(i,j)时刻的主星变轨速度增量大小和方向;根据TL(i,j)时刻和TF(i,j)时刻的主星变轨速度增量大小和方向,获得TL(i,j)时刻和TF(i,j)时刻的主星速度大小和方向;
步骤105、判断第二序数j≤N2是否成立,若成立,则令j的值增加1,然后转入步骤102;否则转入步骤106;
步骤106、判断第一序数i≤N1是否成立,若成立,则令i的值增加1且j等于1,然后转入步骤102;否则转入步骤107;
步骤107、基于步骤104获得的所有时刻主星变轨速度增量大小和主星速度大小,即获得相对导航卫星伴飞机动窗口。
一种相对导航卫星伴飞机动窗口计算装置,包括参数输入模块、数据计算模块、循环迭代模块、显示模块;
所述参数输入模块用于接收外部参数,所述外部参数包括T0时刻主星轨道参数、T1时刻伴星轨道参数、主星点火时刻TL(i)、轨道转移时间TOF(j),然后参数输入模块将外部参数输出给所述循环迭代模块;第一序数i取值范围的最大值N1为点火时刻分布区间内的采样点数量;第二序数j取值范围的最大值N2为主星飞行时间区间内的采样点数量;
所述循环迭代模块用于将接收到的外部参数输出给所述数据计算模块,其中,循环迭代模块将不同的主星点火时刻TL(i)和不同的轨道转移时间TOF(j)进行排列组合;
所述数据计算模块根据所述的主星点火时刻TL(i)和轨道转移时间TOF(j),获得主星变轨时刻TL(i,j)和交会时刻TF(i,j);然后根据T0时刻主星轨道参数和TL(i,j)时刻,获得TL(i,j)时刻主星的轨道参数,根据T1时刻伴星轨道参数和TF(i,j)时刻获得TF(i,j)时刻伴星的轨道参数;最后根据所述TL(i,j)时刻主星的轨道参数和TF(i,j)时刻伴星的轨道参数,采用求解兰伯特问题的方法,获得TL(i,j)时刻和TF(i,j)时刻的主星变轨速度增量大小和方向,根据TL(i,j)时刻和TF(i,j)时刻的主星变轨速度增量大小和方向,获得TL(i,j)时刻和TF(i,j)时刻的主星速度大小和方向;
所述显示模块用于将所述数据计算模块获得的所有时刻主星变轨速度增量大小和主星速度大小进行显示。
一种计算机存储介质,其上存储有计算机程序,该程序被处理器执行时实现上述相对导航卫星伴飞机动窗口计算方法的步骤。
实施例:
一种相对导航卫星伴飞机动窗口计算方法,首先根据主星和伴星的位置和速度参数,计算给定主星点火时刻和轨道转移时间情况下的变轨速度增量和与伴星交会时的相对速度;其次,给定主星点火时刻范围和轨道转移时间范围,以点火时刻为横坐标、轨道转移时间为纵坐标,获取主星加速的速度增量和交会时相对伴星的速度等高线图;最后,根据速度等高线图,得到满足主星变轨速度增量约束的点火时刻和轨道转移时间。此外,根据速度等高线图,还可得到主星变轨速度增量最优和轨道转移时间最优的机动窗口。该发明还可对给定的点火时刻和轨道转移时间(与伴星交会时刻),得到主星点火的方向和速度增量需求,以及交会时相对伴星的速度大小和方向。
具体的,在给定主星轨道机动时刻范围和轨道转移时间范围,即可根据主星和伴星轨道参数,获取主星变轨速度和交会时相对伴星的速度等高线图,解决了主星轨道机动窗口和变轨速度增量整体分布的问题,为空间相对导航任务伴飞构型的构建提供依据。本发明方法包括如下步骤,如图2所示:
(1)确定T0时刻主星轨道参数和T1时刻伴星轨道参数,以及主星点火时刻TL(i)和轨道转移时间TOF(j)。其中:
i=1,2,…,N1
j=1,2,…,N2
(2)根据主星点火时刻TL(i)和轨道转移时间TOF(j),得到主星变轨时刻TL(i,j)和交会时刻TF(i,j);
(3)根据T0时刻主星轨道参数和T1时刻伴星轨道参数,得到TL(i,j)时刻主星的轨道参数和TF(i,j)时刻伴星轨道参数;
(4)根据TL(i,j)时刻主星的轨道参数和TF(i,j)时刻伴星轨道参数,通过求解兰伯特问题,得到TL(i,j)时刻和TF(i,j)时刻的主星变轨速度增量大小和方向,据此得到主星TL(i,j)时刻的速度大小和方向,以及TF(i,j)时刻与伴星交会的速度大小和方向;
(5)判断是否满足j小于等于N2。若是,则令j=j+1,返回第(2)步循环计算;若否,则进入第(6)步;
(6)判断是否满足i小于等于N1。若是,则令i=i+1,j=1,返回第(2)步循环计算;若否,则进入第(7)步;
(7)根据计算得到的变轨时刻矩阵TL和飞行时间分布矩阵TOF,绘制变轨速度和交会速度等高线图,得到变轨窗口分布情况。
图3为本发明实施例的主星变轨速度等高线图。图中每个点均代表一个变轨时机,图中不同曲线代表不同的变轨速度增量需求,同一曲线上的点代表变轨速度增量大小相同。
图4为本发明实施例的交会时刻主星相对伴星的速度等高线图。图中每个点均代表一个变轨时机,图中不同曲线代表不同的相对伴星的速度,同一曲线的点代表相对速度大小相同。
一种相对导航卫星伴飞机动窗口计算方法,通过多次求解轨道中的兰伯特问题,获取点火时刻和交会时刻的速度等高线图,得到轨道机动窗口和速度增量的整体分布,为我国空间相对导航任务的伴飞构型构建窗口选择提供直观依据。
本发明说明书中未作详细描述的内容属本领域技术人员的公知技术。
Claims (4)
1.一种相对导航卫星伴飞机动窗口计算方法,其特征在于:包括如下步骤:
步骤一、确定主星点火时刻TL(i)和轨道转移时间TOF(j),其中第一序数i的初始值为1,第二序数j的初始值为1;
步骤二、根据步骤一中所述的主星点火时刻TL(i)和轨道转移时间TOF(j),获得主星变轨时刻TL(i,j)和交会时刻TF(i,j);
步骤三、根据T0时刻主星轨道参数和步骤二所述TL(i,j)时刻,获得TL(i,j)时刻主星的轨道参数;根据T1时刻伴星轨道参数和步骤二所述TF(i,j)时刻获得TF(i,j)时刻伴星的轨道参数;
步骤四、根据步骤三中所述TL(i,j)时刻主星的轨道参数和TF(i,j)时刻伴星的轨道参数,采用求解兰伯特问题的方法,获得TL(i,j)时刻和TF(i,j)时刻的主星变轨速度增量大小和方向;根据TL(i,j)时刻和TF(i,j)时刻的主星变轨速度增量大小和方向,获得TL(i,j)时刻和TF(i,j)时刻的主星速度大小和方向;
步骤五、判断第二序数j≤N2是否成立,若成立,则令j的值增加1,然后转入步骤二;否则转入步骤六;
步骤六、判断第一序数i≤N1是否成立,若成立,则令i的值增加1且j等于1,然后转入步骤二;否则转入步骤七;
步骤七、基于步骤四获得的所有时刻主星变轨速度增量大小和主星速度大小,即获得相对导航卫星伴飞机动窗口;
所述第一序数i取值范围的最大值N1为点火时刻分布区间内的采样点数量;所述第二序数j取值范围的最大值N2为主星飞行时间区间内的采样点数量。
2.一种相对导航卫星伴飞机动窗口计算装置,其特征在于:包括参数输入模块、数据计算模块、循环迭代模块;
所述参数输入模块用于接收外部参数,所述外部参数包括T0时刻主星轨道参数、T1时刻伴星轨道参数、主星点火时刻TL(i)、轨道转移时间TOF(j),然后参数输入模块将外部参数输出给所述循环迭代模块;
所述循环迭代模块用于将接收到的外部参数输出给所述数据计算模块,其中,循环迭代模块将不同的主星点火时刻TL(i)和不同的轨道转移时间TOF(j)进行排列组合;
所述数据计算模块根据所述的主星点火时刻TL(i)和轨道转移时间TOF(j),获得主星变轨时刻TL(i,j)和交会时刻TF(i,j);然后根据T0时刻主星轨道参数和TL(i,j)时刻,获得TL(i,j)时刻主星的轨道参数,根据T1时刻伴星轨道参数和TF(i,j)时刻获得TF(i,j)时刻伴星的轨道参数;最后根据所述TL(i,j)时刻主星的轨道参数和TF(i,j)时刻伴星的轨道参数,采用求解兰伯特问题的方法,获得TL(i,j)时刻和TF(i,j)时刻的主星变轨速度增量大小和方向,根据TL(i,j)时刻和TF(i,j)时刻的主星变轨速度增量大小和方向,获得TL(i,j)时刻和TF(i,j)时刻的主星速度大小和方向;
第一序数i取值范围的最大值N1为点火时刻分布区间内的采样点数量;
第二序数j取值范围的最大值N2为主星飞行时间区间内的采样点数量。
3.根据权利要求2所述的一种相对导航卫星伴飞机动窗口计算装置,其特征在于:还包括显示模块,所述显示模块用于将所述数据计算模块获得的所有时刻主星变轨速度增量大小和主星速度大小进行显示。
4.一种计算机存储介质,其上存储有计算机程序,其特征在于:该程序被处理器执行时实现权利要求1所述方法的步骤。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811495677.1A CN109725648B (zh) | 2018-12-07 | 2018-12-07 | 一种相对导航卫星伴飞机动窗口计算方法、装置及介质 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811495677.1A CN109725648B (zh) | 2018-12-07 | 2018-12-07 | 一种相对导航卫星伴飞机动窗口计算方法、装置及介质 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109725648A CN109725648A (zh) | 2019-05-07 |
CN109725648B true CN109725648B (zh) | 2020-09-18 |
Family
ID=66295550
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811495677.1A Active CN109725648B (zh) | 2018-12-07 | 2018-12-07 | 一种相对导航卫星伴飞机动窗口计算方法、装置及介质 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109725648B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112000121B (zh) * | 2020-07-14 | 2022-07-12 | 哈尔滨工业大学 | 一种多服务飞行器空间在轨服务燃料最优轨道的设计方法 |
CN114537714B (zh) * | 2022-02-28 | 2022-11-25 | 中国人民解放军63921部队 | 一种高轨卫星变轨调控方法及系统 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107589756A (zh) * | 2017-09-12 | 2018-01-16 | 北京理工大学 | 一种奔月卫星编队初始化方法 |
CN108216687A (zh) * | 2017-12-25 | 2018-06-29 | 中国空间技术研究院 | 基于粒子群算法的geo卫星变轨策略计算方法、系统及介质 |
CN108897023A (zh) * | 2018-04-26 | 2018-11-27 | 北京空间飞行器总体设计部 | 一种星上自主的非合作机动目标跟踪保持变轨方法 |
CN108957499A (zh) * | 2018-05-04 | 2018-12-07 | 北京空间飞行器总体设计部 | 基于观测量频谱分析与最优估计的伴飞目标相对导航方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CA3000108A1 (en) * | 2015-10-02 | 2017-04-06 | Axis Labs Inc. | External motor drive system for window covering system with continuous cord loop |
-
2018
- 2018-12-07 CN CN201811495677.1A patent/CN109725648B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107589756A (zh) * | 2017-09-12 | 2018-01-16 | 北京理工大学 | 一种奔月卫星编队初始化方法 |
CN108216687A (zh) * | 2017-12-25 | 2018-06-29 | 中国空间技术研究院 | 基于粒子群算法的geo卫星变轨策略计算方法、系统及介质 |
CN108897023A (zh) * | 2018-04-26 | 2018-11-27 | 北京空间飞行器总体设计部 | 一种星上自主的非合作机动目标跟踪保持变轨方法 |
CN108957499A (zh) * | 2018-05-04 | 2018-12-07 | 北京空间飞行器总体设计部 | 基于观测量频谱分析与最优估计的伴飞目标相对导航方法 |
Non-Patent Citations (4)
Title |
---|
伴随卫星伴随运动运动学分析;王志刚 等;《空军工程大学学报(自然科学版)》;20010228;全文 * |
小天体探测轨道设计与仿真;李旦伟;《CNKI》;20091231;全文 * |
服务型航天器在轨机动窗口;蔡远文 等;《北京航空航天大学学报》;20101231;全文 * |
近地小行星探测转移轨道方案设计与优化;乔栋 等;《哈尔滨工业大学学报》;20060630;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN109725648A (zh) | 2019-05-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103884237B (zh) | 基于目标概率分布信息的多对一协同制导方法 | |
CN103557867B (zh) | 一种基于稀疏a*搜索的三维多uav协同航迹规划方法 | |
CN103112600B (zh) | 一种星际转移轨道设计方法 | |
CN103136444B (zh) | 一种导航卫星不同控制模式下太阳帆板的受照分析方法 | |
CN109725648B (zh) | 一种相对导航卫星伴飞机动窗口计算方法、装置及介质 | |
CN104374388A (zh) | 一种基于偏振光传感器的航姿测定方法 | |
CN103414451B (zh) | 一种应用于飞行器姿态估计的扩展卡尔曼滤波方法 | |
CN103869701A (zh) | 基于姿态序列解算的飞行器新型实时制导方法 | |
CN103970142A (zh) | 在轨拖曳的组合体航天器姿轨复合控制方法 | |
Wolf et al. | Toward improved landing precision on Mars | |
CN104503471A (zh) | 一种机动飞行器多终端约束反演滑模末制导方法 | |
CN103453907A (zh) | 基于分层大气模型的行星进入段导航滤波方法 | |
US20150207327A1 (en) | Systems and methods for estimating net solar energy production for airborne photovoltaic systems | |
CN105069311A (zh) | 一种远程火箭发射初态误差传播估计方法 | |
CN102508492B (zh) | 一种飞行器在等高航路点间的定高度大圆飞行实现方法 | |
CN104729510A (zh) | 一种空间目标相对伴飞轨道确定方法 | |
Zhang et al. | Receding horizon trajectory optimization with terminal impact specifications | |
CN105138808A (zh) | 基于摄动理论的滑翔弹道误差传播分析方法 | |
CN104677359A (zh) | 一种动态载体运行轨迹预测预报的方法 | |
Way et al. | EDL simulation results for the Mars 2020 landing site safety assessment | |
Martinson | Obstacle avoidance guidance and control algorithm for spacecraft maneuvers | |
CN103440358B (zh) | 一种基于dem数据的坡度拟合方法 | |
CN104820228A (zh) | 一种全球卫星导航系统接收机星座优选方法 | |
Jiang et al. | Multi-target tracking algorithm of vehicle-mounted millimeter-wave radar based on IMM-JPDA-STUKF | |
CN111090830B (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 |