CN111679269B - 一种基于变分的多雷达融合航迹状态估计方法 - Google Patents

一种基于变分的多雷达融合航迹状态估计方法 Download PDF

Info

Publication number
CN111679269B
CN111679269B CN202010560351.3A CN202010560351A CN111679269B CN 111679269 B CN111679269 B CN 111679269B CN 202010560351 A CN202010560351 A CN 202010560351A CN 111679269 B CN111679269 B CN 111679269B
Authority
CN
China
Prior art keywords
track
target
fusion
variation
radar
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
Application number
CN202010560351.3A
Other languages
English (en)
Other versions
CN111679269A (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.)
Strategic Early Warning Research Institute Of People's Liberation Army Air Force Research Institute
Original Assignee
Strategic Early Warning Research Institute Of People's Liberation Army Air Force Research Institute
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 Strategic Early Warning Research Institute Of People's Liberation Army Air Force Research Institute filed Critical Strategic Early Warning Research Institute Of People's Liberation Army Air Force Research Institute
Priority to CN202010560351.3A priority Critical patent/CN111679269B/zh
Publication of CN111679269A publication Critical patent/CN111679269A/zh
Application granted granted Critical
Publication of CN111679269B publication Critical patent/CN111679269B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S13/00Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/02Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
    • G01S13/50Systems of measurement based on relative movement of target
    • G01S13/58Velocity or trajectory determination systems; Sense-of-movement determination systems
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S13/00Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/66Radar-tracking systems; Analogous systems
    • G01S13/72Radar-tracking systems; Analogous systems for two-dimensional tracking, e.g. combination of angle and range tracking, track-while-scan radar
    • G01S13/723Radar-tracking systems; Analogous systems for two-dimensional tracking, e.g. combination of angle and range tracking, track-while-scan radar by using numerical data
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S13/00Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/87Combinations of radar systems, e.g. primary radar and secondary radar
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Mathematical Physics (AREA)
  • Theoretical Computer Science (AREA)
  • Data Mining & Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Pure & Applied Mathematics (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Computational Mathematics (AREA)
  • Algebra (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明涉及一种基于变分的多雷达融合航迹状态估计方法,属于多传感器目标跟踪与数据融合领域。本发明为了克服雷达情报处理系统在融合航迹时常在多个原始航迹间来回摇摆,呈锯齿状,影响航迹质量,以及常规方法易出现转弯“鼓包”甚至跟踪失败现象的问题,通过对多雷达的观测信息进行统一管理,构建包含目标运动速度变换曲线的变分模型,并给出了变分模型的快速求解方式,求解的运动速度曲线用于所有观测雷达的航迹外推,最后对目标状态估计值进行估计更新。本发明最大程度上考虑了观测信息的相关性,降低航迹相对偏差和目标突然机动对目标估计的不利影响,提高了目标运动速度和融合位置估计的精度,提高了航迹质量,获得了对目标的可靠跟踪。

Description

一种基于变分的多雷达融合航迹状态估计方法
技术领域
本发明属于多传感器目标跟踪与数据融合技术领域,具体涉及一种基于变分的多雷达融合航迹状态估计方法。
背景技术
目前各级雷达情报处理系统普遍采用分布式处理架构,即各雷达装备分别进行目标探测和航迹跟踪,产生单雷达原始航迹,然后送至上级融合处理中心,由中心完成各原始航迹的关联和融合,并最终生成组网融合航迹。这种架构中上下级系统松耦合,通信量少,便于工程实现。
因探测周期不同、开机时刻不一致等,各雷达对目标的探测时刻难以同步,融合中心需对各原始航迹进行时间配准(航迹外推至同一时刻),再加权平均得到融合航迹点。但由于雷达探测精度不一,网络时延不同,探测环境易受电磁、气象、季节等因素影响,以及装备使用过程中逐渐累积的标校误差等因素,探测数据中系统偏差和随机误差并存,原始航迹间存在较大差异,时间配准步骤中难以有效利用其他航迹的观测信息,融合航迹常在多个原始航迹间来回摇摆,呈锯齿状,严重影响航迹质量。此外,常规方法也难以对强机动目标进行稳健跟踪,易出现转弯“鼓包”甚至跟踪失败现象。
发明内容
(一)要解决的技术问题
本发明要解决的技术问题是如何提供一种基于变分的多雷达融合航迹状态估计方法,以克服现有的雷达情报处理系统在融合航迹时常在多个原始航迹间来回摇摆,呈锯齿状,影响航迹质量,以及常规方法易出现转弯“鼓包”甚至跟踪失败现象的问题。
(二)技术方案
为了解决上述技术问题,本发明提出一种基于变分的多雷达融合航迹状态估计方法,该方法包括如下步骤:
步骤1:创建观测数据队列;
步骤2:向所述队列中添加实时观测数据;
步骤3:建立并求解目标X轴上的运动速度的变分模型;
步骤4:X轴上的原始航迹位置外推;
步骤5:X轴上的融合航迹位置估计;
步骤6:按照与步骤3-5相同的方式,计算Y、Z轴上的位置估计。
进一步地,对每批融合航迹,在起批之初都创建一个动态更新的观测队列,用以保存其最新N秒时间内各原始雷达的观测移动记录。
进一步地,N为20-60。
进一步地,所述步骤2具体包括:一旦相应的原始航迹进行了位置更新,都向队列中添加记录(tk,1,tk,2,dk,x,dk,y,dk,z),其中,k>0表示在其队列中的索引编号,tk,2、tk,1分别是雷达最新航迹点和上一个航迹点的记录时刻,dk,x,dk,y,dk,z分别为该段时间内目标在X、Y、Z轴上的移动距离。
进一步地,所述步骤3中建立目标X轴上的运动速度的变分模型具体包括:具体包括:设当前融合航迹的状态估计时刻为tc,为求目标X轴上的运动速度曲线Vx,可选取最近N秒时间内的队列数据,建立变分模型:
式(1)中,等号右侧第一项为正则项,使Vx尽量平滑,等号右侧第二项为保真项,使得Vx趋同于观测到的数据记录,使得上式代价最小的Vx即为最终解;表示差分,λ∈[0,1]为正则项加权因子,αk为第k条观测记录的加权系数,与雷达装备的探测性能密切相关,
进一步地,所述λ值越小,表示模型对目标的机动适应能力越强,但取值过小,可能会导致过拟合现象。
进一步地,所述步骤3中求解目标X轴上的运动速度的变分模型具体包括:使用迭代法求解所述变分模型,先生成初始值然后逐渐更新生成/> 使得求解结果越来越精确,具体过程如下:
(1)取为观测记录的加权:
(2)根据多项式分解,计算下一步的理论调整步长
式中,n为大于0的整数;
(3)更新
式中,β为防止过拟合而设置的更新步长因子,0<β<1;
(4)不断重复步骤(2)和步骤(3),直到趋于稳定,则最终的/>即为所求的Vx(t)。
进一步地,所述步骤4具体包括:对相应的M路原始航迹进行外推,预测其在当前时刻tc的位置:
式中,1≤m≤M,tlast分别是第m路原始航迹中最后的观测时刻和观测位置。
进一步地,所述步骤5具体包括:根据所述步骤4的各原始航迹外推结果,加权得到融合航迹的位置估计:
式中,∑mγm=1,γm为加权系数,与雷达装备的探测性能密切相关。
进一步地,所述步骤6具体包括:按照与所述步骤3-5相同的方式,计算融合航迹在Y轴、Z轴上的位置估计
(三)有益效果
本发明提出一种基于变分的多雷达融合航迹状态估计方法,通过对多雷达的观测信息进行统一管理,构建包含目标运动速度变换曲线的变分模型,并给出了变分模型的快速求解方式,求解的运动速度曲线用于所有观测雷达的航迹外推,最后对目标状态估计值进行估计更新。本发明最大程度上考虑了观测信息的相关性,降低航迹相对偏差和目标突然机动对目标估计的不利影响,提高了目标运动速度和融合位置估计的精度,提高了航迹质量,获得了对目标的可靠跟踪。
附图说明
图1为本发明的基于变分的多雷达融合航迹状态估计方法流程图。
具体实施方式
为使本发明的目的、内容和优点更加清楚,下面结合附图和实施例,对本发明的具体实施方式作进一步详细描述。
本发明提出一种基于变分的多雷达融合航迹状态估计方法,包括如下步骤:
步骤1:创建观测数据队列
对每批融合航迹,在起批之初都创建一个动态更新的观测队列,用以保存其最新N秒时间内各原始雷达的观测移动记录。N一般可取20—60;
步骤2:向队列中添加实时观测数据
一旦相应的原始航迹进行了位置更新,都向队列中添加记录(tk,1,tk,2,dk,x,dk,y,dk,z),其中,k>0表示在其队列中的索引编号,tk,2、tk,1分别是雷达最新航迹点和上一个航迹点的记录时刻,dk,x,dk,y,dk,z分别为该段时间内目标在X、Y、Z轴上的移动距离,坐标系可采用直角坐标或经纬度坐标。
需要说明的是,雷达测量误差主要包含技术限制所致的固有随机误差和因装备老化、天线失校等因素引起的系统偏差,系统偏差均值不为零,但在不同方位距离上可能存在较大差异,在目标飞行过程中,目标航迹的系统偏差呈缓慢变化,本发明中队列数据采用移动距离,而非位置数据本身,能有效消除系统偏差对融合估计的不利干扰。
步骤3:建立并求解目标运动速度的变分模型
设当前融合航迹的状态估计时刻为tc,为求目标X轴上的运动速度曲线Vx,可选取最近N秒时间内的队列数据,建立变分模型:
式中,式中第一项为正则项,使Vx尽量平滑,第二项为保真项,使得Vx趋同于观测到的数据记录,使得上式代价最小的Vx即为最终解。表示差分,λ∈[0,1]为正则项加权因子,值越小,表示模型对目标的机动适应能力越强,但取值过小,可能会导致过拟合现象。对一般飞行目标,可取值默认0.8,αk为第k条观测记录的加权系数,与雷达装备的探测性能密切相关,/>
本发明使用迭代法求解,先生成初始值然后逐渐更新生成/> 等,使得求解结果越来越精确,过程取下:
(1)取为观测记录的加权:
(2)根据多项式分解,计算下一步的理论调整步长
式中,n为大于0的整数;
(3)更新
式中,β为防止过拟合而设置的更新步长因子,0<β<1,本发明默认取值0.4;
(4)不断重复步骤(2)和步骤(3),直到趋于稳定,则最终的/>即为所求的Vx(t)。
步骤4:原始航迹位置外推
对相应的M路原始航迹进行外推,预测其在当前时刻tc的位置:
式中,1≤m≤M,tlast分别是第m路原始航迹中最后的观测时刻和观测位置。
步骤5:融合航迹位置估计
根据各原始航迹外推结果,加权得到融合航迹的位置估计:
式中,∑mγm=1,γm为加权系数,与雷达装备的探测性能密切相关,可参考公式1中的αk取值。
步骤6:计算Y、Z轴上的位置估计
按照与步骤3—5相同的方式,计算融合航迹在Y轴、Z轴上的位置估计
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明技术原理的前提下,还可以做出若干改进和变形,这些改进和变形也应视为本发明的保护范围。

Claims (7)

1.一种基于变分的多雷达融合航迹状态估计方法,其特征在于,该方法包括如下步骤:
步骤1:创建观测数据队列;
步骤2:向所述队列中添加实时观测数据;
步骤3:建立并求解目标X轴上的运动速度的变分模型;
步骤4:X轴上的原始航迹位置外推;
步骤5:X轴上的融合航迹位置估计;
步骤6:按照与步骤3-5相同的方式,计算Y、Z轴上的位置估计;
其中,
所述步骤1具体包括:对每批融合航迹,在起批之初都创建一个动态更新的观测队列,用以保存其最新N秒时间内各原始雷达的观测移动记录;
所述步骤2具体包括:一旦相应的原始航迹进行了位置更新,都向队列中添加记录(tk,1,tk,2,dk,x,dk,y,dk,z),其中,k>0表示在其队列中的索引编号,tk,2、tk,1分别是雷达最新航迹点和上一个航迹点的记录时刻,dk,x,dk,y,dk,z分别为该段时间内目标在X、Y、Z轴上的移动距离;
所述步骤3中建立目标X轴上的运动速度的变分模型具体包括:设当前融合航迹的状态估计时刻为tc,为求目标X轴上的运动速度曲线Vx,可选取最近N秒时间内的队列数据,建立变分模型:
式(1)中,等号右侧第一项为正则项,使Vx尽量平滑,等号右侧第二项为保真项,使得Vx趋同于观测到的数据记录,使得上式代价最小的Vx即为最终解;表示差分,λ∈[0,1]为正则项加权因子,αk为第k条观测记录的加权系数,与雷达装备的探测性能密切相关,
2.如权利要求1所述的基于变分的多雷达融合航迹状态估计方法,其特征在于,N为20-60。
3.如权利要求1所述的基于变分的多雷达融合航迹状态估计方法,其特征在于,所述λ值越小,表示模型对目标的机动适应能力越强。
4.如权利要求1所述的基于变分的多雷达融合航迹状态估计方法,其特征在于,所述步骤3中求解目标X轴上的运动速度的变分模型具体包括:使用迭代法求解所述变分模型,先生成初始值然后逐渐更新生成/> 使得求解结果越来越精确,具体过程如下:
(1)取为观测记录的加权:
(2)根据多项式分解,计算下一步的理论调整步长
式中,n是迭代次数,为大于0的整数;
(3)更新
式中,β为防止过拟合而设置的更新步长因子,0<β<1;
(4)不断重复步骤(2)和步骤(3),直到趋于稳定,则最终的/>即为所求的Vx(t)。
5.如权利要求3或4所述的基于变分的多雷达融合航迹状态估计方法,其特征在于,所述步骤4具体包括:对相应的M路原始航迹进行外推,预测其在当前时刻tc的位置:
式中,1≤m≤M,tlast分别是第m路原始航迹中最后的观测时刻和观测位置。
6.如权利要求5所述的基于变分的多雷达融合航迹状态估计方法,其特征在于,所述步骤5具体包括:根据所述步骤4的各原始航迹外推结果,加权得到融合航迹的位置估计:
式中,∑mγm=1,γm为加权系数,与雷达装备的探测性能密切相关。
7.如权利要求6所述的基于变分的多雷达融合航迹状态估计方法,其特征在于,所述步骤6具体包括:按照与所述步骤3-5相同的方式,计算融合航迹在Y轴、Z轴上的位置估计
CN202010560351.3A 2020-06-18 2020-06-18 一种基于变分的多雷达融合航迹状态估计方法 Active CN111679269B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010560351.3A CN111679269B (zh) 2020-06-18 2020-06-18 一种基于变分的多雷达融合航迹状态估计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010560351.3A CN111679269B (zh) 2020-06-18 2020-06-18 一种基于变分的多雷达融合航迹状态估计方法

Publications (2)

Publication Number Publication Date
CN111679269A CN111679269A (zh) 2020-09-18
CN111679269B true CN111679269B (zh) 2023-09-29

Family

ID=72436508

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010560351.3A Active CN111679269B (zh) 2020-06-18 2020-06-18 一种基于变分的多雷达融合航迹状态估计方法

Country Status (1)

Country Link
CN (1) CN111679269B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112433849B (zh) * 2020-11-11 2022-09-02 成都民航空管科技发展有限公司 一种集成塔台系统多监视源数据分布式融合的方法及系统
CN116520307B (zh) * 2023-02-20 2023-10-27 扬州宇安电子科技有限公司 一种基于目标追踪的雷达协同检测系统及方法

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105788252A (zh) * 2016-03-22 2016-07-20 连云港杰瑞电子有限公司 基于定点检测器和信号配时数据融合的城市干道车辆轨迹重构方法
CN106371091A (zh) * 2016-08-24 2017-02-01 四川九洲空管科技有限责任公司 Ads‑b与二次雷达监视信息的数据融合方法及装置
CN106680806A (zh) * 2016-11-24 2017-05-17 清华大学 一种多雷达点迹融合方法
CN108646732A (zh) * 2018-04-20 2018-10-12 华东交通大学 基于驾驶人操控意图的车辆轨迹预测方法、装置及系统
CN109655792A (zh) * 2018-12-12 2019-04-19 北京无线电测量研究所 一种雷达数据并行处理方法、装置及存储介质
CN110068805A (zh) * 2019-05-05 2019-07-30 中国人民解放军国防科技大学 基于变分贝叶斯推论的高速目标hrrp重构方法
CN111178385A (zh) * 2019-12-02 2020-05-19 江苏大学 一种鲁棒在线多传感器融合的目标跟踪方法
WO2020116156A1 (ja) * 2018-12-03 2020-06-11 日立オートモティブシステムズ株式会社 車両制御装置

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10310068B2 (en) * 2014-12-08 2019-06-04 Northrop Grumman Systems Corporation Variational track management

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105788252A (zh) * 2016-03-22 2016-07-20 连云港杰瑞电子有限公司 基于定点检测器和信号配时数据融合的城市干道车辆轨迹重构方法
CN106371091A (zh) * 2016-08-24 2017-02-01 四川九洲空管科技有限责任公司 Ads‑b与二次雷达监视信息的数据融合方法及装置
CN106680806A (zh) * 2016-11-24 2017-05-17 清华大学 一种多雷达点迹融合方法
CN108646732A (zh) * 2018-04-20 2018-10-12 华东交通大学 基于驾驶人操控意图的车辆轨迹预测方法、装置及系统
WO2020116156A1 (ja) * 2018-12-03 2020-06-11 日立オートモティブシステムズ株式会社 車両制御装置
CN109655792A (zh) * 2018-12-12 2019-04-19 北京无线电测量研究所 一种雷达数据并行处理方法、装置及存储介质
CN110068805A (zh) * 2019-05-05 2019-07-30 中国人民解放军国防科技大学 基于变分贝叶斯推论的高速目标hrrp重构方法
CN111178385A (zh) * 2019-12-02 2020-05-19 江苏大学 一种鲁棒在线多传感器融合的目标跟踪方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
一种模板匹配的快速实现方法;李超 等;航天返回与遥感;第37卷(第01期);第63-70页 *
多雷达系统数据融合与航迹预测;吴显亮 等;数学的实践与认识;第40卷(第15期);第160-173页 *

Also Published As

Publication number Publication date
CN111679269A (zh) 2020-09-18

Similar Documents

Publication Publication Date Title
CN109633590B (zh) 基于gp-vsmm-jpda的扩展目标跟踪方法
CN111679269B (zh) 一种基于变分的多雷达融合航迹状态估计方法
CN107315171B (zh) 一种雷达组网目标状态与系统误差联合估计算法
CN108255791B (zh) 基于分布式传感器一致性的机动目标跟踪方法
WO2008038889A1 (en) Radar tracking device and method thereof
WO2024001629A1 (zh) 一种面向智能驾驶车辆的多传感器融合方法及系统
CN115808682B (zh) 基于车载毫米波雷达的稳定多目标跟踪方法及系统
CN112034445B (zh) 基于毫米波雷达的车辆运动轨迹跟踪方法和系统
CN109917373B (zh) 运动补偿搜索的动平台雷达的动态规划检测前跟踪方法
CN112230195B (zh) 基于非线性最优平滑层策略的平滑变结构滤波方法及系统
CN110133586A (zh) 基于线性校正的toa联合同步与定位方法
CN114139109A (zh) 一种目标跟踪方法、系统、设备、介质及数据处理终端
CN112067007B (zh) 地图生成方法、计算机存储介质及电子设备
CN116878510A (zh) 基于因子图优化的多源融合室内定位方法、装置及系统
CN110007298B (zh) 一种目标超前预测跟踪方法
CN116576868A (zh) 一种多传感器融合精确定位及自主导航方法
CN115884368A (zh) 一种基于高斯滤波器信念传播的协作式车辆定位方法
CN110031797B (zh) 用于被动传感系统对具有非连续特性目标的检测跟踪方法
CN112285697A (zh) 一种多传感器多目标空时偏差校准与融合方法
Jianxing et al. A Particle Filtering Method for Radar Target Tracking
CN116930949A (zh) 一种星载雷达网协同目标探测系统及方法
CN117250606B (zh) 轨迹跟踪方法、装置、设备及存储介质
CN115392117B (zh) 一种水下高速机动平台高帧率无模糊声学导航方法
CN112966720B (zh) 一种基于blue的雷达与红外观测数据融合方法、系统
CN112526450B (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