CN109358651A - 一种直升机动态航线规划方法 - Google Patents
一种直升机动态航线规划方法 Download PDFInfo
- Publication number
- CN109358651A CN109358651A CN201811568504.8A CN201811568504A CN109358651A CN 109358651 A CN109358651 A CN 109358651A CN 201811568504 A CN201811568504 A CN 201811568504A CN 109358651 A CN109358651 A CN 109358651A
- Authority
- CN
- China
- Prior art keywords
- threat
- course line
- coordinate
- track points
- new
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 9
- 238000004891 communication Methods 0.000 claims abstract description 27
- 238000013439 planning Methods 0.000 claims abstract description 17
- 238000005516 engineering process Methods 0.000 claims abstract description 9
- 238000013210 evaluation model Methods 0.000 claims abstract description 4
- 230000035772 mutation Effects 0.000 claims abstract description 4
- 230000001174 ascending effect Effects 0.000 claims description 5
- 238000006073 displacement reaction Methods 0.000 claims description 4
- 235000007164 Oryza sativa Nutrition 0.000 claims description 3
- 238000001514 detection method Methods 0.000 claims description 3
- 235000009566 rice Nutrition 0.000 claims description 3
- 240000007594 Oryza sativa Species 0.000 claims 1
- 238000012546 transfer Methods 0.000 abstract description 4
- 230000006835 compression Effects 0.000 abstract description 2
- 238000007906 compression Methods 0.000 abstract description 2
- 238000005457 optimization Methods 0.000 abstract description 2
- 230000005540 biological transmission Effects 0.000 abstract 1
- 241000209094 Oryza Species 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 1
- 238000005336 cracking Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05D—SYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
- G05D1/00—Control of position, course, altitude or attitude of land, water, air or space vehicles, e.g. using automatic pilots
- G05D1/10—Simultaneous control of position or course in three dimensions
- G05D1/101—Simultaneous control of position or course in three dimensions specially adapted for aircraft
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C21/00—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
- G01C21/20—Instruments for performing navigational calculations
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Automation & Control Theory (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Aviation & Aerospace Engineering (AREA)
- Traffic Control Systems (AREA)
Abstract
本发明涉及一种直升机动态航线规划方法,其特征在于采用以下步骤:构建通讯信号威胁模型f1和地形威胁模型f2,并由通讯信号威胁模型f1和地形表面威胁模型f2加权叠加,构成综合威胁模型F;构建动态任务序列,实现新威胁信息的隐蔽链动态传递,由同组直升机获取新威胁参数、进行飞行路径规划,包括生成初始航线集合,对航线进行交叉、复制、变异操作,根据航线评价模型fv计算飞行代价,最终实现航线可飞性优化。本发明考虑到同组直升机信息传递问题进行航线协同规划,采用航迹点数据快速压缩技术实现信息传递的及时性、隐蔽性及安全性,同时去中心化的隐蔽传递技术也为实现直升机航线实时规划提供了有效的技术保障。
Description
技术领域
本发明属于直升机航线规划领域,涉及到使用直升机性能和环境威胁等信息规划出最优飞行航线方法。
背景技术
到目前为止,国内外在复杂环境下航线规划方面的研究还比较初步,尚没有形成完整而系统的直升机动态飞行航线规划模型与方法。复杂地理环境的直升机动态航线规划目前还没有得到很好的解决,如何能应用隐蔽传递技术来规划安全性能最大的飞行航线,以及实现动态航空救援协同规划是亟待解决的问题,此问题的解决无论对军用还是民用都具有极其重要的应用价值。
发明内容
本发明的目的是提供一种在复杂地理环境下、能应用隐蔽传递技术来规划安全性能最大的飞行航线、以实现动态航空救援协同规划的直升机动态航线规划方法。包括以下步骤:
1)构建通讯信号威胁模型f1:
nxi=(xti+xti+1)/2
其中,公式(1)中n是航路中的航迹点总数,xn=Xe,yn=Ye,yhn=He,(xti,yti,hti)和(Xe,Ye,He)分别表示第i个通讯威胁的航迹点、目标点坐标,TN是通讯威胁的总数,(xck,yck)和Rk分别表示第k个通讯威胁的中心点平面坐标和通讯威胁影响半径,(nxi,nyi)表示航迹点回避威胁目标后调整的新平面坐标;
(2)构建地形威胁模型f2,其数学描述如下:
公式(2)中,(xd,yd,hd)、(xoi,yoi,Hmaxi)分别表示地形威胁点的三维坐标和第i个山峰峰顶的三维坐标,(xsi,ysi)表示山峰随着x、y方向下降的调节参数,M表示山峰的总数;
(3)由通讯信号威胁模型f1和地形表面威胁模型f2加权叠加,构成综合威胁模型F,其数学描述如下:
p1+p2=1
由此可得:
其中,F为航线综合威胁影响,(xt,yt,ht)表示通讯威胁的航迹点坐标,p1、p2分别为地形威胁权重和航程约束权重,f1、f2分别为通讯威胁和地形威胁;
(4)构建动态任务序列,新威胁信息的隐蔽链动态传递
根据机载雷达识别的新威胁构建新的任务信息序列,并加入隐蔽链传递末端,新威胁信息采用下式表达:
Tm(N)={xc,yc,hc,flag} (4)
Tm={Tm(1),Tm(2),…,Tm(N)}
公式(4)中,Tm(N)表示第N个任务信息,(xc,yc,hc)表示新威胁的三维坐标,纵横坐标(xc,yc)单位为公里,高程hc单位为米;flag表示新威胁解决情况,flag取值为0时表示新威胁得到解决,取值为1时表示新威胁尚未解决;通过hash计算实现新威胁信息加密;
(5)获取新威胁参数
邻机采用隐蔽传递的信息探测技术,并对hash加密的威胁信息进行解密,获取动态规划任务参数:新威胁(xc,yc,hc)的任务完成情况为flag,当flag值为0时,整个航线规划完毕;当flag值为1时,进入步骤(6);
(6)生成初始航线集合,根据航线评价模型fv计算飞行代价,若飞行代价满足飞行安全的要求并小于最大航程,执行步骤(7),否则对航线进行交叉、复制、变异操作,并重复步骤(6):
fpbest={(xi,yi,hi)HI},i=[1,2,i,…,n] (6)
公式(5)~(6)中,H(j)i为第i条航线的第j个航迹点的高程,最小航迹点高程HJmin为最小航迹点高程差之和,m为航线集内航线的总数,HI为最小航迹点高程差之和对应的航线序号,fpbest为最优航线集,(xi,yi,hi)HI指对应于序号HI航线集第i个航迹点的三维坐标。
(7)航线可飞性优化
首先对航迹点按照x坐标由小到大进行排序,此时y和h坐标随着x坐标进行相同的移动;或者按照y坐标由小到大进行排序,此时x和h坐标随着x坐标进行相同的移动:
公式(7)中排序之前航线的表示为P,第i个航迹点坐标为(xi,yi,hi),排序由sort表示,排序后的航线为Ps,排序后的第i个航迹点坐标为(xi′,yi′,hi′);
然后去除航线的冗余航迹点:
公式(8)~(9)中li表示航迹点P(xi-1,yi-1,hi-1)和P(xi,yi,hi)之间航线的长度,LS和lmin分别表示最小航迹段集合、长度值,lCi-1、lCi+1表示航迹段向前和向后合并的航迹段长度值;
最后结合直升机的性能对上述航线针对转弯角、俯仰角、乘坐品质及B-Spline航线光滑性进行优化,即得最优飞行航线。
本发明与现有技术相比,优点在于:
a.实现隐蔽环境下动态航线规划。军机在隐蔽环境下飞行为避免通信信号被敌方捕获带来极大的安全隐患,或民用飞机在飞行中发现突发威胁提醒队友及时避让,都要求实现信息的及时性、隐蔽性传递,隐蔽传递技术的去中心化特点为解决该问题提供了一种有效的方法。
b.利用航迹点数据快速压缩技术实现航线优化。通过航线规划算法选择的航迹点存在大量冗余及不满足可飞性的问题,通过对航迹点坐标按照x或y坐标值的大小进行排序,对减少航线往复绕行,航迹段不满足航迹段长度约束进行改进,能够提高航线的可飞性,缩短航程有重要作用。
具体实施方式
下面结合实施例对本发明技术方案做进一步说明。实施例采用如下假设:有2架米-28武装直升机飞行,其中直升机1在飞行中探测到威胁的位置,且对危险并未处理,而是将危险信息上传至隐蔽传递,直升机2从隐蔽传递中获取威胁信息参数并进行处理。在(20,50)、(15,25)、(20,25)、(10,15)、(35,35)、(20,40)、(5,17)、(13,76)、(36,48)、(29,41)、(51,19)、(48,27)、(19,82)和(24,41)处共有14个地形威胁;通信威胁(55,55)、(10,50)、(15,25)、(30,25)、(10,15)、(40,15)和(40,10)处7个;航迹段的最小长度为5km,最大航程470千米。航线规划的参数为:待选航线集个数100,每条航线的航迹点个数为12,航线起始点坐标、新威胁的坐标分别为(1km,1km,150m)和(60km,60km,150m),新威胁未处理:flag=0。具体步骤为:
(1)根据假设构建通讯信号威胁模型f1:
nxi=(xti+xti+1)/2
(2)构建地形威胁模型f2,其数学描述如下:
f(hdi)=100*exp(-(xdi-10)^2./5^2-(ydi-50)^2./5^2)+100*exp(-(xdi-15)^2./7^2-(ydi-25)^2./6^2)+200*exp(-(xdi-20)^2./4^2-(ydi-25)^2./5^2)+250*exp(-(xdi-10)^2./7^2-(ydi-15)^2./6^2)+300*exp(-(xdi-35)^2/3^2-(ydi-35)^2/4^2)+240*exp(-(xdi-20)^2/3^2-(ydi-40)^2/6^2)+360*exp(-(xdi-5)^2/8^2-(ydi-17)^2/9^2)+410*exp(-(xdi-13)^2/2^2-(ydi-76)^2/4^2)+190*exp(-(xdi-36)^2/5^2-(ydi-48)^2/6^2)+360*exp(-(xdi-29)^2/9^2-(ydi-41)^2/3^2)+380*exp(-(xdi-51)^2/8^2-(ydi-19)^2/8^2)+250*exp(-(xdi-48)^2/4^2-(ydi-27)^2/6^2)+180*exp(-(xdi-19)^2/4^2-(ydi-82)^2/7^2)+170*exp(-(xdi-24)^2/5^2-(ydi-41)^2/4^2) (2)
(3)由通讯信号威胁模型f1和地形表面威胁模型f2加权叠加,构成综合威胁模型F,其数学描述如下:
其中,F为航线综合威胁影响,(xt,yt,ht)表示通讯威胁的航迹点坐标,p1、p2分别为地形威胁权重和航程约束权重,f1、f2分别为通讯威胁和地形威胁;
(4)构建动态任务序列,新威胁信息的隐蔽链动态传递
根据机载雷达识别的新威胁构建新的任务信息序列,并加入隐蔽链传递末端,新威胁信息采用下式表达:
Tm(1)={100,100,200,1},Tm(2)={60,60,150,0} (4)
Tm={Tm(1),Tm(2)}
(5)获取新威胁参数
邻机采用隐蔽传递的信息探测技术,并对hash加密的威胁信息进行解密,获取动态规划任务参数:新威胁(60km,60km,150m)的任务完成情况为flag,flag值为0,进入步骤(6);
(6)生成初始航线集合,根据航线评价模型fv计算飞行代价,对航线进行交叉、复制、变异操作,从能够回避地形威胁和通讯威胁的航线中选出代价最小的最优的初始航线fpbest:
fpbest={(xi,yi,hi)HI},i=[1,2,i,…,12] (6)
(7)航线可飞性优化
首先对航迹点按照x坐标由小到大进行排序,此时y和h坐标随着x坐标进行相同的移动:
Ps{(x′1,y′1,h′1),(x′2,y′2,h′2),…(x′12,y′12,h′12)}=sort(P{(x1,y1,h1),(x2,y2,h2)…(x12,y12,h12)},x) (7)
公式(7)中排序之前航线的表示为P,第i个航迹点坐标为(xi,yi,hi),排序由sort表示,排序后的航线为Ps,排序后的第i个航迹点坐标为(xi′,yi′,hi′);
然后去除航线的冗余航迹点:
公式(8)~(9)中li表示航迹点P(xi-1,yi-1,hi-1)和P(xi,yi,hi)之间航线的长度,LS表示最小航迹段集合,lCi-1、lCi+1表示航迹段向前和向后合并的航迹段长度值;
最后结合直升机的性能对上述航线针对转弯角、俯仰角、乘坐品质及B-Spline航线光滑性进行优化,即得最优飞行航线。
表1实施例的通讯和地形威胁参数
表2直升机动态航迹计算和总航程评价结果
表2显示了:实验中设置地形威胁和通信信号威胁具有相同的权重0.5,规划航程为112.55公里,小于最大航程,规划时间短,仅为14.2秒,满足实时航线规划的需要,迭代代数为1,说明航线具有很快的寻优能力。
Claims (1)
1.一种直升机动态航线规划方法,其特征在于采用以下步骤:
(1)构建通讯信号威胁模型f1:
nxi=(xti+xti+1)/2
其中,公式(1)中n是航路中的航迹点总数,xn=Xe,yn=Ye,yhn=He,(xti,yti,hti)和(Xe,Ye,He)分别表示第i个通讯威胁的航迹点、目标点坐标,TN是通讯威胁的总数,(xck,yck)和Rk分别表示第k个通讯威胁的中心点平面坐标和通讯威胁影响半径,(nxi,nyi)表示航迹点回避威胁目标后调整的新平面坐标;
(2)构建地形威胁模型f2,其数学描述如下:
公式(2)中,(xd,yd,hd)、(xoi,yoi,Hmaxi)分别表示地形威胁点的三维坐标和第i个山峰峰顶的三维坐标,(xsi,ysi)表示山峰随着x、y方向下降的调节参数,M表示山峰的总数;
(3)由通讯信号威胁模型f1和地形表面威胁模型f2加权叠加,构成综合威胁模型F,其数学描述如下:
p1+p2=1
由此可得:
其中,F为航线综合威胁影响,(xt,yt,ht)表示通讯威胁的航迹点坐标,p1、p2分别为地形威胁权重和航程约束权重,f1、f2分别为通讯威胁和地形威胁;
(4)构建动态任务序列,新威胁信息的隐蔽链动态传递
根据机载雷达识别的新威胁构建新的任务信息序列,并加入隐蔽链传递末端,新威胁信息采用下式表达:
Tm(N)={xc,yc,hc,flag} (4)
Tm={Tm(1),Tm(2),…,Tm(N)}
公式(4)中,Tm(N)表示第N个任务信息,(xc,yc,hc)表示新威胁的三维坐标,纵横坐标(xc,yc)单位为公里,高程hc单位为米;flag表示新威胁解决情况,flag取值为0时表示新威胁得到解决,取值为1时表示新威胁尚未解决;通过hash计算实现新威胁信息加密;
(5)获取新威胁参数
邻机采用隐蔽传递的信息探测技术,并对hash加密的威胁信息进行解密,获取动态规划任务参数:新威胁(xc,yc,hc)的任务完成情况为flag,当flag值为0时,整个航线规划完毕;当flag值为1时,进入步骤(6);
(6)生成初始航线集合,根据航线评价模型fv计算飞行代价,若飞行代价满足飞行安全的要求并小于最大航程,执行步骤(7),否则对航线进行交叉、复制、变异操作,并重复步骤(6):
fpbest={(xi,yi,hi)HI},i=[1,2,i,…,n] (6)
公式(5)~(6)中,H(j)i为第i条航线的第j个航迹点的高程,最小航迹点高程HJmin为最小航迹点高程差之和,m为航线集内航线的总数,HI为最小航迹点高程差之和对应的航线序号,fpbest为最优航线集,(xi,yi,hi)HI指对应于序号HI航线集第i个航迹点的三维坐标。
(7)航线可飞性优化
首先对航迹点按照x坐标由小到大进行排序,此时y和h坐标随着x坐标进行相同的移动;或者按照y坐标由小到大进行排序,此时x和h坐标随着x坐标进行相同的移动:
公式(7)中排序之前航线的表示为P,第i个航迹点坐标为(xi,yi,hi),排序由sort表示,排序后的航线为Ps,排序后的第i个航迹点坐标为(xi′,yi′,hi′);
然后去除航线的冗余航迹点:
公式(8)~(9)中li表示航迹点P(xi-1,yi-1,hi-1)和P(xi,yi,hi)之间航线的长度,LS和lmin分别表示最小航迹段集合、长度值,lCi-1、lCi+1表示航迹段向前和向后合并的航迹段长度值;
最后结合直升机的性能对上述航线针对转弯角、俯仰角、乘坐品质及B-Spline航线光滑性进行优化,即得最优飞行航线。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811568504.8A CN109358651B (zh) | 2018-12-21 | 2018-12-21 | 一种直升机动态航线规划方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811568504.8A CN109358651B (zh) | 2018-12-21 | 2018-12-21 | 一种直升机动态航线规划方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109358651A true CN109358651A (zh) | 2019-02-19 |
CN109358651B CN109358651B (zh) | 2021-07-30 |
Family
ID=65329351
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811568504.8A Expired - Fee Related CN109358651B (zh) | 2018-12-21 | 2018-12-21 | 一种直升机动态航线规划方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109358651B (zh) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105700553A (zh) * | 2016-01-28 | 2016-06-22 | 中国科学院自动化研究所 | 多无人机自主协同决策快速集成系统 |
CN105807266A (zh) * | 2016-05-19 | 2016-07-27 | 中国人民解放军军械工程学院 | 预警雷达航迹数据传输压缩方法 |
US9500488B2 (en) * | 2012-03-14 | 2016-11-22 | Lockheed Martin Corporation | Route re-planning using enemy force lethality projection |
CN108353081A (zh) * | 2015-09-28 | 2018-07-31 | 13部门有限公司 | 无人机入侵检测和对抗 |
CN108592921A (zh) * | 2018-05-02 | 2018-09-28 | 山东理工大学 | 一种分段最速下降混合航线规划方法 |
-
2018
- 2018-12-21 CN CN201811568504.8A patent/CN109358651B/zh not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9500488B2 (en) * | 2012-03-14 | 2016-11-22 | Lockheed Martin Corporation | Route re-planning using enemy force lethality projection |
CN108353081A (zh) * | 2015-09-28 | 2018-07-31 | 13部门有限公司 | 无人机入侵检测和对抗 |
CN105700553A (zh) * | 2016-01-28 | 2016-06-22 | 中国科学院自动化研究所 | 多无人机自主协同决策快速集成系统 |
CN105807266A (zh) * | 2016-05-19 | 2016-07-27 | 中国人民解放军军械工程学院 | 预警雷达航迹数据传输压缩方法 |
CN108592921A (zh) * | 2018-05-02 | 2018-09-28 | 山东理工大学 | 一种分段最速下降混合航线规划方法 |
Non-Patent Citations (6)
Title |
---|
GUANGHUI WANG,XUEFENG SUN,LIPING ZHANG,CHAO LV: "Saturation attack based route planning and threat avoidance algorithm for cruise missiles", 《JOURNAL OF SYSTEMS ENGINEERING AND ELECTRONICS》 * |
XIANGYIN ZHANG, HAIBIN DUAN: "An improved constrained differential evolution algorithm for unmanned aerial vehicle global route planning", 《APPLIED SOFT COMPUTING》 * |
ZHOU SHUDAO, WANG JUN, JIN YONGQI: "Route Planning for Unmanned Aircraft Based on Ant Colony Optimization and Voronoi Diagram", 《11 SECOND INTERNATIONAL CONFERENCE ON INTELLIGENT SYSTEM DESIGN AND ENGINEERING APPLICATION》 * |
刘丽峰,张树清: "利用改进威胁模型的电势理论的三维飞行路径规划", 《中国科学院研究生院学报》 * |
田疆: "无人机三维约束多目标航迹规划", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 * |
郭琳,周德云,汪凌霄: "低可探测性飞机威胁建模及航迹规划仿真研究", 《计算机仿真》 * |
Also Published As
Publication number | Publication date |
---|---|
CN109358651B (zh) | 2021-07-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110609552B (zh) | 一种水下无人航行器编队平面航迹规划方法 | |
CN103913172B (zh) | 一种适用于复杂低空下飞行器的路径规划方法 | |
CN102175245B (zh) | 一种基于海流历史统计信息的水下潜器路径规划方法 | |
CN101477169B (zh) | 巡检飞行机器人对电力线路的检测方法 | |
CN1947032B (zh) | 海上地震测量方法和系统 | |
CN102759357B (zh) | 通信延迟下多无人机协同实时航路规划方法 | |
CN101286071B (zh) | 基于微粒群优化和遗传算法的多无人机三维编队重构方法 | |
CN110031004A (zh) | 基于数字地图的无人机静态和动态路径规划方法 | |
CN104991895A (zh) | 一种基于三维空域网格的低空救援航空器航迹规划方法 | |
CN102591358A (zh) | 一种多无人机的动态编队控制方法 | |
CN102749080A (zh) | 一种基于流体力学的无人机三维航路生成方法 | |
CN104317305B (zh) | 一种面向复杂战场威胁的射前航迹确定方法 | |
CN103557867A (zh) | 一种基于稀疏a*搜索的三维多uav协同航迹规划方法 | |
Wang et al. | Ship route planning based on double-cycling genetic algorithm considering ship maneuverability constraint | |
CN109445463A (zh) | 一种无人机动态航线规划方法 | |
Liu et al. | Three-dimensional mountain complex terrain and heterogeneous multi-UAV cooperative combat mission planning | |
CN108762280A (zh) | 一种考虑海洋环流影响的uuv基于能量消耗优化的远程航海路径规划方法 | |
Malaek et al. | Dynamic based cost functions for TF/TA flights | |
CN113359737A (zh) | 一种考虑队形伸缩的船舶编队自适应事件触发控制方法 | |
CN101825723B (zh) | 海上地震测量方法和系统 | |
Shan | Study on submarine path planning based on modified ant colony optimization algorithm | |
Chen et al. | Three dimensional path planning of UAV based on adaptive particle swarm optimization algorithm | |
CN111762341A (zh) | 小天体形貌特征在线提取的安全附着多模切换制导方法 | |
CN108592921B (zh) | 一种分段最速下降混合航线规划方法 | |
CN109358651A (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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20210730 Termination date: 20211221 |
|
CF01 | Termination of patent right due to non-payment of annual fee |