CN114332418A - 一种目标轨迹模拟方法及装置 - Google Patents
一种目标轨迹模拟方法及装置 Download PDFInfo
- Publication number
- CN114332418A CN114332418A CN202111565184.2A CN202111565184A CN114332418A CN 114332418 A CN114332418 A CN 114332418A CN 202111565184 A CN202111565184 A CN 202111565184A CN 114332418 A CN114332418 A CN 114332418A
- Authority
- CN
- China
- Prior art keywords
- target
- current
- time step
- path point
- speed
- 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
Images
Landscapes
- Navigation (AREA)
Abstract
本发明公开了一种目标轨迹模拟方法及装置,通过设置时间步长,对于当前时间步长,计算当前时间步长结束时目标速度以及对应的加速度和角速度,并据此计算当前时间步长内平均速度和平均角速度,从而预测出本时间步长剩余时间内最大运动距离和最大转向角,并计算下一个路径点相对目标的方位角、俯仰角和斜距,确定出目标运动朝向的方位角,然后运动到下一个位置,更新当前目标状态,在运动到最后一个路径点时结束运动模拟。本发明对于目标运动的模拟更加精确。
Description
技术领域
本申请属于目标轨迹模拟技术领域,尤其涉及一种目标轨迹模拟方法及装置。
背景技术
在侦察***的装备功能性能试验中,往往需要利用合作目标进行试验。合作目标可采用真实目标和数字仿真目标。使用真实目标(如无人机)测试更接近实际的作战场景,但测试成本较高,测试过程也较为复杂。数字仿真目标轨迹模拟会建立目标的运动模型,根据目标的固有属性和运动趋势计算运动姿态,并发送给需要测试的装备或模拟器。
现有技术对目标运动轨迹的模拟大多基于wgs-84或类似的地理坐标系,均为三维直角坐标系。在距离较近时,动目标在地面上按最短路径到达某点的运动轨迹接近直线,在直角坐标系下对目标轨迹进行模拟可以得出近似解。但侦察***区域范围较大,考虑到地球是一个椭球体,地面目标与低慢小飞行器的运动轨迹并非直线,仍在直角坐标系下模拟目标轨迹会产生较大误差。
发明内容
本申请的目的是提供一种目标轨迹模拟方法及装置,以克服在直角坐标系下模拟目标轨迹会产生较大误差的问题。
为了实现上述目的,本申请技术方案如下:
一种目标轨迹模拟方法,包括:
步骤S1、设置时间步长,输入目标的属性及经过的路径点,初始化目标的当前状态;
步骤S2、对于当前时间步长,计算当前时间步长结束时目标速度以及对应的加速度和角速度,并据此计算当前时间步长内平均速度和平均角速度,从而预测出本时间步长剩余时间内最大运动距离和最大转向角;
步骤S3、将目标当前的位置和下一个路径点的经纬度代入Vincenty逆解公式,计算下一个路径点相对目标的方位角、俯仰角和斜距;
步骤S4、比较下一个路径点相对目标的方位角与本时间步长剩余时间内最大转向角,确定出目标运动朝向的方位角;
步骤S5、判断目标当前位置与下一个路径点之间的斜距是否大于本时间步长剩余时间内最大运动距离,如果大于,转入步骤S6,否则,转入步骤7;
步骤S6、目标按照下一个路径点相对目标的俯仰角和目标运动朝向的方位角运动,运动距离为所述最大运动距离,并通过Vincenty正解公式计算出当前时间步长结束时的结束方位和位置,输出目标当前的运动状态,重置剩余时间转入步骤S2进行下一个时间步长的计算;
步骤S7、目标按照下一个路径点相对目标的俯仰角和目标运动朝向的方位角运动,运动距离为所述斜距,并通过Vincenty正解公式计算出到达下一路径点时的结束方位和位置,判断是否还存在下一个路径点,如果存在,则计算剩余时间后返回步骤S3,否则输出目标当前的运动状态,目标轨迹模拟完成。
进一步的,所述对于当前时间步长,计算当前时间步长结束时目标速度以及对应的加速度和角速度,并据此计算当前时间步长内平均速度和平均角速度,从而预测出本时间步长剩余时间内最大运动距离和最大转向角,包括:
根据目标的当前速度vc和加速度ac计算出当前步长可加速到的最快速度vm=vc+act,若vm>vmax,则当前时间步长结束时目标的速度vg=vmax;若vm≤vmax,则当前时间步长结束时目标的速度vg=vm;vmax表示的目标的最大速度;
进一步的,所述比较下一个路径点相对目标的方位角与本时间步长剩余时间内最大转向角,确定出目标运动朝向的方位角,包括:
进一步的,所述目标在一次运动的过程中,加速度不变,俯仰角不变。
进一步的,所述目标在一次运动结束后,用结束方位和位置更新目标当前的方位和位置。
本申请还提出了一种目标轨迹模拟装置,包括处理器以及存储有若干计算机指令的存储器,所述计算机指令被处理器执行时实现所述目标轨迹模拟方法的步骤。
本申请提出的一种目标轨迹模拟方法及装置,对于目标运动的模拟更加精确,更接近真实目标在地球上沿一系列路径点按最短路径运动的实际运动轨迹,能更好的验证对应装备或模拟器的功能性能。
附图说明
图1为本申请目标轨迹模拟方法流程图;
图2为本申请实施例目标位置关系示意图。
具体实施方式
为了使本申请的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本申请进行进一步详细说明。应当理解,此处描述的具体实施例仅用以解释本申请,并不用于限定本申请。
大地主题解算是大地主题正算和反算问题的统称,是椭球面大地测量学研究的核心问题之一,也是作战仿真领域兵力实***置推算和弹道计算的核心算法。其中,Vincenty公式是ThaddeusVincenty于1975年提出的一种以椭球体来计算两点距离的方法,它是基于地球椭球面上经纬度和两点间距离方位转换计算求解大地解算问题,其采用贝塞尔球作为过渡,先将椭球面各大地元素按特定要求转换到辅助球面,并在球面上解算,然后通过椭球改正项把球面解算的结果再归算到椭球面上。
Vincenty公式以贝塞尔公式为基础,推导出嵌套系数公式,可适用于任何距离的大地主题解算,包括正解公式和反解公式。Vincenty正解公式,即已知一点的经纬度坐标及它与另一点的距离和方位,求另一点的经纬度坐标;Vincenty反解公式,即已知两点的经纬度坐标,求它们之间的距离和方位。Vincenty公式的数学表达式已经是公开的技术,这里不再赘述。本申请将Vincenty公式应用于目标轨迹的模拟,提高目标轨迹模拟的精度。
目标轨迹的模拟过程基于时间步长,用户可设置最短为1ms的仿真步长,***时钟按照此时间步长依次推进推演进程,在每个时间步长内,可以忽略目标及态势的状态变化,近似的认为目标及态势未发生变化,在每个时间步长结束时,更新目标状态。
在一个实施例中,如图1所示,提供了一种目标轨迹模拟方法,包括:
步骤S1、设置时间步长,输入目标的属性及经过的路径点,初始化目标的当前状态。
本申请技术方案中,将所有目标视作质点,每个步长目标的运动视为在地面(或近地)上沿某方位前进的匀(加)速运动。
首先输入目标的属性(类型、初始方位角、加速度a、最大速度vmax、转向半径r)及经过的路径点p1、p2…pi…pm(经纬高),有m个路径点。目标类型包括人员、车辆、动物及低慢小物体,路径点可在gis地图上点选或手动输入经纬高坐标。
不同类型的目标属性上会有区别,如人的最大速度会比车小很多,使用者也可以自行修改属性。路径点是目标需要经过的一组坐标点,包含经度、纬度和高度信息,至少要有2个路径点。
初始化目标的当前状态(当前方位角当前速度vc=0、当前角速度ωc=0、当前俯仰角Θc=0、当前加速度ac=a、当前位置pc=p1、下一路径点pi=p2)与当前步长剩余运动时间tc=t,t为设置的时间步长。目标初始位置在第一个路径点处,朝第二个路径点运动。
步骤S2、对于当前时间步长,计算当前时间步长结束时目标速度以及对应的加速度和角速度,并据此计算当前时间步长内平均速度和平均角速度,从而预测出本时间步长剩余时间内最大运动距离和最大转向角。
在达到最高速度vmax前,目标会尽量加速。
在第一个时间步长时,记录目标当前速度vc=0,当前加速度ac=a。根据目标的当前速度vc和加速度ac计算出当前步长可加速到的最快速度vm=vc+at,若vm>vmax,则在当前时间步长结束时目标加速到最快速度,即当前时间步长结束时目标的速度vg=vmax。若vm≤vmax,则当前时间步长结束时目标的速度vg=vm。
在其他时间步长时,记录目标当前速度vc,加速度为ac。根据目标的当前速度vc和加速度ac计算出当前步长可加速到的最快速度vm=vc+act,若vm>vmax,则在当前时间步长结束时目标加速到最快速度,即当前时间步长结束时目标的速度vg=vmax。若vm≤vmax,则当前时间步长结束时目标的速度vg=vm。
步骤S3、将目标当前的位置和下一个路径点的经纬度代入Vincenty逆解公式,计算下一个路径点相对目标的方位角、俯仰角和斜距。
然后根据下一个路径点pi与目标当前位置pc的高度差h计算出俯仰角Θp以及斜距stgt:
Θp=tan-1(h/sh)
步骤S4、比较下一个路径点相对目标的方位角与本时间步长剩余时间内最大转向角,确定出目标运动朝向的方位角。
目标会尽量转向pi相对pc的方位若说明目标可以转到方位则目标运动朝向的方位角若则转动最大转向角转向目标运动朝向的方位角为当前方位角减去最大转向角,即(逆时针),或为当前方位角加上最大转向角即(顺时针)。
步骤S5、判断目标当前位置与下一个路径点之间的斜距是否大于本时间步长剩余时间内最大运动距离,如果大于,转入步骤S6,否则,转入步骤7。
本步骤为判断步骤,判断在当前时间步长内能否到达下一路径点。即比较目标与下一路径点pi的斜距stgt和当前步长的最大运动距离sm,若stgt>sm,说明剩余运动时间tc内目标无法到达下一路径点pi,否则说明剩余运动时间tc内目标可以到达下一路径点pi。
步骤S6、目标按照下一个路径点相对目标的俯仰角和目标运动朝向的方位角运动,运动距离为所述最大运动距离,并通过Vincenty正解公式计算出当前时间步长结束时的结束方位和位置,输出目标当前的运动状态,重置剩余时间转入步骤S2进行下一个时间步长的计算。
方位角和俯仰角确定之后,则目标就开始当前时间步长的运动,运动距离为最大运动距离sm。在本次运动就是沿该方位角和俯仰角运动的,目标的每次运动都认为是在球面上沿某一方位俯仰近似直线的匀变速运动。加速度、俯仰角在一次运动内都认为是不变的,方位角随着角速度计算会有微小的变化。
根据俯仰角Θp计算目标在地面上的运动距离sl=sm×cosΘp。将目标当前位置pc的经纬度、当前方位角地面上的运动距离sl代入Vincenty正解公式,计算目标在当前时间步长结束时到达位置的结束方位与位置pe,更新当前朝向当前位置pc=pe。
此时当前时间步长计算完成,更新vc=vg,vc=ωg,ac=ag,Θc=Θp。输出目标当前的运动状态(当前朝向当前速度vc、当前角速度ωc、当前俯仰Θc、当前加速度ac、当前位置pc),下一时间步长重置剩余时间tc=t,返回步骤S2。
步骤S7、目标按照下一个路径点相对目标的俯仰角和目标运动朝向的方位角运动,运动距离为所述斜距,并通过Vincenty正解公式计算出到达下一路径点时的结束方位和位置,判断是否还存在下一个路径点,如果存在,则计算剩余时间后返回步骤S3,否则输出目标当前的运动状态,目标轨迹模拟完成。
根据俯仰Θp计算目标在地面上的运动距离sl=stgt×cosΘp。将目标当前位置pc的经纬度、当前朝向地面上的运动距离sl代入Vincenty正解公式,计算目标在到达路径点pi时的结束方位与坐标pe,更新当前朝向当前位置p=pe。
判断路径点pc后是否还存在下一个路径点pi+1。若存在pi+1,则下一路径点pi更新为pi+1,更新当前步长剩余运动时间更新vc=vg,ωc=ωg,ac=ag,Θc=Θp回到步骤S3,继续当前时间步长内的计算。需要说明的是,返回步骤S3后,还在当前时间步长内,继续当前时间步长内的运动,即在当前时间步长内有两次运动。
在生成运动距离跨度较大的运动目标时,本申请目标轨迹模拟更接近实际目标在地球上的运动轨迹。
为了验证本申请技术方案的计算精度,可以连续进行Vincenty公式的正解与反解计算,将正解计算的结果作为反解公式的初值进行反算,比较反解公式结果与正算初值之间的偏差,就可以反映计算精度。取地球表面面上距离不同的两点,计算偏差结果如下表所示:
距离(km) | 方位角(°) | 距离偏差(m) | 角度偏差(°) |
5 | 30 | 0.0091 | 0.00016 |
10 | 30 | 0.018 | 0.00016 |
20 | 30 | 0.037 | 0.00015 |
50 | 30 | 0.091 | 0.00015 |
表1
考虑性能原因,Vincenty公式在实现中迭代次数较少,误差比理论较大,但在50公里内,最大距离偏差仅为0.091m,角度偏差极小,足以满足***模拟需要。
上述Vincenty计算过程中椭球(地球模型)各参数定义如下:a为椭球长半轴,b为椭球短半轴,f为扁率f=1-b/a,在WGS-84坐标系中,取a=6378137,b=6356752.3142,f=1/298.257223563。
在另一个实施例中,还提出了一种目标轨迹模拟装置,包括处理器以及存储有若干计算机指令的存储器,所述计算机指令被处理器执行时实现所述目标轨迹模拟方法的步骤。
关于目标轨迹模拟装置的具体限定可以参见上文中对于目标轨迹模拟方法的限定,在此不再赘述。上述目标轨迹模拟装置可全部或部分通过软件、硬件及其组合来实现。可以硬件形式内嵌于或独立于计算机设备中的处理器中,也可以以软件形式存储于计算机设备中的存储器中,以便于处理器调用执行以上对应的操作。
存储器和处理器之间直接或间接地电性连接,以实现数据的传输或交互。例如,这些元件互相之间可以通过一条或多条通讯总线或信号线实现电性连接。存储器中存储有可在处理器上运行的计算机程序,所述处理器通过运行存储在存储器内的计算机程序,从而实现本发明实施例中的网络拓扑布局方法。
其中,所述存储器可以是,但不限于,随机存取存储器(Random Access Memory,RAM),只读存储器(Read Only Memory,ROM),可编程只读存储器(Programmable Read-OnlyMemory,PROM),可擦除只读存储器(Erasable Programmable Read-Only Memory,EPROM),电可擦除只读存储器(Electric Erasable Programmable Read-Only Memory,EEPROM)等。其中,存储器用于存储程序,所述处理器在接收到执行指令后,执行所述程序。
所述处理器可能是一种集成电路芯片,具有数据的处理能力。上述的处理器可以是通用处理器,包括中央处理器(Central Processing Unit,CPU)、网络处理器(NetworkProcessor,NP)等。可以实现或者执行本发明实施例中公开的各方法、步骤及逻辑框图。通用处理器可以是微处理器或者该处理器也可以是任何常规的处理器等。
以上所述实施例仅表达了本申请的几种实施方式,其描述较为具体和详细,但并不能因此而理解为对发明专利范围的限制。应当指出的是,对于本领域的普通技术人员来说,在不脱离本申请构思的前提下,还可以做出若干变形和改进,这些都属于本申请的保护范围。因此,本申请专利的保护范围应以所附权利要求为准。
Claims (6)
1.一种目标轨迹模拟方法,其特征在于,所述目标轨迹模拟方法,包括:
步骤S1、设置时间步长,输入目标的属性及经过的路径点,初始化目标的当前状态;
步骤S2、对于当前时间步长,计算当前时间步长结束时目标速度以及对应的加速度和角速度,并据此计算当前时间步长内平均速度和平均角速度,从而预测出本时间步长剩余时间内最大运动距离和最大转向角;
步骤S3、将目标当前的位置和下一个路径点的经纬度代入Vincenty逆解公式,计算下一个路径点相对目标的方位角、俯仰角和斜距;
步骤S4、比较下一个路径点相对目标的方位角与本时间步长剩余时间内最大转向角,确定出目标运动朝向的方位角;
步骤S5、判断目标当前位置与下一个路径点之间的斜距是否大于本时间步长剩余时间内最大运动距离,如果大于,转入步骤S6,否则,转入步骤7;
步骤S6、目标按照下一个路径点相对目标的俯仰角和目标运动朝向的方位角运动,运动距离为所述最大运动距离,并通过Vincenty正解公式计算出当前时间步长结束时的结束方位和位置,输出目标当前的运动状态,重置剩余时间转入步骤S2进行下一个时间步长的计算;
步骤S7、目标按照下一个路径点相对目标的俯仰角和目标运动朝向的方位角运动,运动距离为所述斜距,并通过Vincenty正解公式计算出到达下一路径点时的结束方位和位置,判断是否还存在下一个路径点,如果存在,则计算剩余时间后返回步骤S3,否则输出目标当前的运动状态,目标轨迹模拟完成。
2.根据权利要求1所述的目标轨迹模拟方法,其特征在于,所述对于当前时间步长,计算当前时间步长结束时目标速度以及对应的加速度和角速度,并据此计算当前时间步长内平均速度和平均角速度,从而预测出本时间步长剩余时间内最大运动距离和最大转向角,包括:
根据目标的当前速度vc和加速度ac计算出当前步长可加速到的最快速度vm=vc+act,若vm>vmax,则当前时间步长结束时目标的速度vg=vmax;若vm≤vmax,则当前时间步长结束时目标的速度vg=vm;vmax表示的目标的最大速度;
4.根据权利要求1所述的目标轨迹模拟方法,其特征在于,所述目标在一次运动的过程中,加速度不变,俯仰角不变。
5.根据权利要求1所述的目标轨迹模拟方法,其特征在于,所述目标在一次运动结束后,用结束方位和位置更新目标当前的方位和位置。
6.一种目标轨迹模拟装置,包括处理器以及存储有若干计算机指令的存储器,其特征在于,所述计算机指令被处理器执行时实现权利要求1至权利要求5中任意一项所述方法的步骤。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111565184.2A CN114332418B (zh) | 2021-12-20 | 2021-12-20 | 一种目标轨迹模拟方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111565184.2A CN114332418B (zh) | 2021-12-20 | 2021-12-20 | 一种目标轨迹模拟方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114332418A true CN114332418A (zh) | 2022-04-12 |
CN114332418B CN114332418B (zh) | 2023-01-03 |
Family
ID=81052152
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111565184.2A Active CN114332418B (zh) | 2021-12-20 | 2021-12-20 | 一种目标轨迹模拟方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114332418B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116883633A (zh) * | 2023-09-07 | 2023-10-13 | 之江实验室 | 一种基于gis模拟智能体移动及可调速的可视化方法和装置 |
CN117407319A (zh) * | 2023-12-15 | 2024-01-16 | 航天中认软件测评科技(北京)有限责任公司 | 用于雷达***软件测试的目标生成方法、装置及电子设备 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101615301A (zh) * | 2009-07-29 | 2009-12-30 | 腾讯科技(深圳)有限公司 | 计算机虚拟环境中目标的路径控制方法和*** |
CN103674028A (zh) * | 2013-12-27 | 2014-03-26 | 上海大唐移动通信设备有限公司 | 一种室内行进轨迹的定位测试方法及装置 |
CN113581213A (zh) * | 2021-08-31 | 2021-11-02 | 北京轻舟智航科技有限公司 | 一种预测车辆运动轨迹的处理方法 |
-
2021
- 2021-12-20 CN CN202111565184.2A patent/CN114332418B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101615301A (zh) * | 2009-07-29 | 2009-12-30 | 腾讯科技(深圳)有限公司 | 计算机虚拟环境中目标的路径控制方法和*** |
CN103674028A (zh) * | 2013-12-27 | 2014-03-26 | 上海大唐移动通信设备有限公司 | 一种室内行进轨迹的定位测试方法及装置 |
CN113581213A (zh) * | 2021-08-31 | 2021-11-02 | 北京轻舟智航科技有限公司 | 一种预测车辆运动轨迹的处理方法 |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116883633A (zh) * | 2023-09-07 | 2023-10-13 | 之江实验室 | 一种基于gis模拟智能体移动及可调速的可视化方法和装置 |
CN116883633B (zh) * | 2023-09-07 | 2023-12-12 | 之江实验室 | 一种基于gis模拟智能体移动及可调速的可视化方法和装置 |
CN117407319A (zh) * | 2023-12-15 | 2024-01-16 | 航天中认软件测评科技(北京)有限责任公司 | 用于雷达***软件测试的目标生成方法、装置及电子设备 |
CN117407319B (zh) * | 2023-12-15 | 2024-04-02 | 航天中认软件测评科技(北京)有限责任公司 | 用于雷达***软件测试的目标生成方法、装置及电子设备 |
Also Published As
Publication number | Publication date |
---|---|
CN114332418B (zh) | 2023-01-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107544067B (zh) | 一种基于高斯混合近似的高超声速再入飞行器跟踪方法 | |
CN114332418B (zh) | 一种目标轨迹模拟方法及装置 | |
CN110146909A (zh) | 一种定位数据处理方法 | |
CN109059907B (zh) | 轨迹数据处理方法、装置、计算机设备和存储介质 | |
CN110609290A (zh) | 激光雷达匹配定位方法及装置 | |
CN104019828A (zh) | 高动态环境下惯性导航***杆臂效应误差在线标定方法 | |
US10656259B2 (en) | Method for determining trajectories of moving physical objects in a space on the basis of sensor data of a plurality of sensors | |
CN112441260B (zh) | 用于航天器自主脱轨的***和方法 | |
RU2571567C2 (ru) | Способ управления движением аэробаллистического летательного аппарата по заданной пространственной траектории | |
CN109045699B (zh) | 虚拟比赛的方法、装置、计算机设备和存储介质 | |
CN113984044A (zh) | 一种基于车载多感知融合的车辆位姿获取方法及装置 | |
US20210192777A1 (en) | Method, device and storage medium for positioning object | |
CN115493601A (zh) | 基于路网匹配的车辆自主定位方法、装置、存储介质 | |
CN105242252A (zh) | 基于图像匹配的下降轨聚束sar雷达定位方法 | |
CN111812669A (zh) | 绕机检查装置及其定位方法、存储介质 | |
CN109857127A (zh) | 训练神经网络模型以及飞行器姿态解算的方法、装置 | |
CN106885567B (zh) | 一种惯导协作定位方法及定位设备 | |
CN110794434B (zh) | 一种位姿的确定方法、装置、设备及存储介质 | |
CN110411449B (zh) | 一种航空侦察载荷目标定位方法、***及终端设备 | |
CN116805047A (zh) | 多传感器融合定位的不确定性表达方法、装置及电子设备 | |
KR102489865B1 (ko) | 다중 필터 및 센서 퓨전 기반의 차량의 위치 추정 방법 | |
CN114763998A (zh) | 基于微型雷达阵列的未知环境并行导航方法和*** | |
CN114660587A (zh) | 基于Jerk模型的跳跃滑翔弹道目标跟踪方法及*** | |
CN113761662A (zh) | 一种滑翔类目标的轨迹预测管道的生成方法 | |
CN116227237B (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 |