CN113237492B - 一种初相稳定性判定方法、存储介质及电子设备 - Google Patents
一种初相稳定性判定方法、存储介质及电子设备 Download PDFInfo
- Publication number
- CN113237492B CN113237492B CN202110456445.0A CN202110456445A CN113237492B CN 113237492 B CN113237492 B CN 113237492B CN 202110456445 A CN202110456445 A CN 202110456445A CN 113237492 B CN113237492 B CN 113237492B
- Authority
- CN
- China
- Prior art keywords
- ref
- phase
- phi
- vlbi
- calculating
- 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 31
- 238000000819 phase cycle Methods 0.000 claims abstract description 14
- 238000005070 sampling Methods 0.000 claims description 3
- 238000005259 measurement Methods 0.000 description 8
- 238000005305 interferometry Methods 0.000 description 6
- 230000009286 beneficial effect Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C25/00—Manufacturing, calibrating, cleaning, or repairing instruments or devices referred to in the other groups of this subclass
Landscapes
- Engineering & Computer Science (AREA)
- Manufacturing & Machinery (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Measuring Phase Differences (AREA)
Abstract
本发明涉及信息获取与处理,特别涉及一种初相稳定性判定方法、存储介质及电子设备。该方法包括:步骤1,生成基准信号Sref;步骤2,按预设时间段连续获取VLBI记录系统的记录数据Sn(t);步骤3,对所述基准信号Sref以及所述记录数据Sn(t)进行互相关计算;步骤4,根据互相关计算结果Scorn(t)进行残余相位序列计算;步骤5,根据残余相位序列的计算结果Φn(t)计算相位变化率均值φref;步骤6,根据相位变化率均值φref计算预设时间段端点时刻的差分相位以及基准差分相位Φref;步骤7,根据所述差分相位以及所述基准差分相位Φref判断VLBI记录系统的初相是否稳定。本发明能够解决对于初相稳定性的判定方法准确度较低的问题。
Description
技术领域
本发明涉及信息获取与处理技术领域,尤其涉及一种初相稳定性判定方法、存储介质及电子设备。
背景技术
VLBI技术是支撑深空探测器轨道测定的重要技术手段,通常采用分时采集记录模式,为有效抑制系统误差,通常采用时域差分的工作模式,即要求系统时延、相位在差分弧段内保持稳定,采集记录设备的初相稳定性是影响VLBI测量精度的重要因素,直接影响系统测量精度,但现有技术中对于初相稳定性的判定方法准确度较低。
发明内容
本发明所要解决的技术问题是提供一种初相稳定性判定方法、存储介质及电子设备。
本发明解决上述技术问题的技术方案如下:一种基于VLBI记录系统的初相稳定性判定方法,包括:
步骤1,生成基准信号Sref;
步骤2,按预设时间段连续获取VLBI记录系统的记录数据Sn(t);
步骤3,对所述基准信号Sref以及所述记录数据Sn(t)进行互相关计算;
步骤4,根据互相关计算结果Scorn(t)进行残余相位序列计算;
步骤5,根据残余相位序列的计算结果Φn(t)计算相位变化率均值φref;
本发明的有益效果是:通过利用高稳频率源产生基准信号Sref,VLBI记录系统分时段、变参数采集记录数据Sn(t),对记录数据Sn(t)与基准信号Sref进行互相关计算,获取全部采集序列的残余相位信息,通过拟合残余相位,判断系统初相稳定性指标,可以实现对单站记录设备的初相绝对稳定性的评估,还可以有效抑制系统误差,并在深空测控干涉测量领域具有显著的应用前景。
在上述技术方案的基础上,本发明还可以做如下改进。
进一步,步骤1具体为:
通过高稳频率源生成基准信号Sref,基准信号Sref的计算公式为:
Sref=exp(j·2π·fref·t)
其中,fref为基准信号频率,t为信号时间,j为虚数单位。
进一步,步骤2具体为:
获取VLBI记录系统的记录数据Sn(t)的公式为:
Sn(t)=exp[j·2π·(fref-fLn)·tM]
其中,n表示第n个时间段,fLn为n个时间段的本振频率,
tM为该时间段内的采样点序列,Fn为第n个时间段的采样频率。
进一步,步骤3具体为:
对所述基准信号Sref以及所述记录数据进行互相关计算,得到互相关计算结果Scorn(t),公式为:
其中,conj()为复共轭计算,M为进行互相关计算的数据点数,Sn(ti)为第n个预设时间段的VLBI记录系统的记录数据。
进一步,步骤4具体为:
根据互相关计算结果Scorn(t)进行残余相位序列计算,得到残余相位序列的计算结果Φn(t),计算公式为:
Φn(t)=angle{Scorn(t)}
其中,angle{}为取幅角计算。
进一步,步骤5具体为:
根据残余相位序列的计算结果Φn(t)计算相位变化率均值φref,公式为:
其中,N为全部N个时间段,R为每个时间段内互相关计算得到的残余相位数,ti为第i个时间段的起始时刻,Φn(ti+1)为第n个预设时间内的第i+1个残余相位,Φn(ti)为第n个预设时间内的第i个残余相位。
进一步,步骤6具体为:
Φref=φref·(tn+1,1-tn,M)
其中,tn+1,1为第n+1个时间段的第一个点的时间,tn,M为第n个时间段的最后一个点的时间,Φn+1(t1)为第n+1个预设时间内的第一个残余相位,Φn(tM)为第n个预设时间内的第M个残余相位。
进一步,步骤7具体为:
本发明解决上述技术问题的另一种技术方案如下:一种存储介质,所述存储介质中存储有指令,当计算机读取所述指令时,使所述计算机执行如上述任一项所述的一种基于VLBI记录系统的初相稳定性判定方法。
本发明的有益效果是:通过利用高稳频率源产生基准信号Sref,VLBI记录系统分时段、变参数采集记录数据Sn(t),对记录数据Sn(t)与基准信号Sref进行互相关计算,获取全部采集序列的残余相位信息,通过拟合残余相位,判断系统初相稳定性指标,可以实现对单站记录设备的初相绝对稳定性的评估,还可以有效抑制系统误差,并在深空测控干涉测量领域具有显著的应用前景。
本发明解决上述技术问题的另一种技术方案如下:一种电子设备,包括存储器、处理器及存储在所述存储器上并在所述处理器上运行的程序,所述处理器执行所述程序时实现如上述任一项所述的一种基于VLBI记录系统的初相稳定性判定方法的步骤。
本发明的有益效果是:通过利用高稳频率源产生基准信号Sref,VLBI记录系统分时段、变参数采集记录数据Sn(t),对记录数据Sn(t)与基准信号Sref进行互相关计算,获取全部采集序列的残余相位信息,通过拟合残余相位,判断系统初相稳定性指标,可以实现对单站记录设备的初相绝对稳定性的评估,还可以有效抑制系统误差,并在深空测控干涉测量领域具有显著的应用前景。
附图说明
图1为本发明一种初相稳定性判定方法实施例提供的流程示意图;
图2为本发明一种初相稳定性判定方法实施例提供的实际采集数据获取的残余相位结果示意图。
具体实施方式
以下对本发明的原理和特征进行描述,所举实例只用于解释本发明,并非用于限定本发明的范围。
如图1所示,一种基于VLBI记录系统的初相稳定性判定方法,包括:
步骤1,生成基准信号Sref;
步骤2,按预设时间段连续获取VLBI记录系统的记录数据Sn(t);
步骤3,对基准信号Sref以及记录数据Sn(t)进行互相关计算;
步骤4,根据互相关计算结果Scorn(t)进行残余相位序列计算;
步骤5,根据残余相位序列的计算结果Φn(t)计算相位变化率均值φref;
在一些可能的实施方式中,通过利用高稳频率源产生基准信号Sref,VLBI记录系统分时段、变参数采集记录数据Sn(t),对记录数据Sn(t)与基准信号Sref进行互相关计算,获取全部采集序列的残余相位信息,通过拟合残余相位,判断系统初相稳定性指标,可以实现对单站记录设备的初相绝对稳定性的评估,还可以有效抑制系统误差,并在深空测控干涉测量领域具有显著的应用前景,获取的残余相位可参考图2。
需要注意的是,全部数据包含N个预设时间段,每个预设时间段内连续记录数据,每个预设时间段内分成若干个子数据段,每个子数据段的原始数据点数为M,每个子数据段对应一个相位值。通常采用高稳频率源生成基准信号Sref,为了便于理解,预设时间段以1s为例,按预设时间段连续获取VLBI记录系统的记录数据Sn(t)可以理解为:按照分时以及变参数的形式采集记录数据,按工作模式,系统是一段时间一段时间的记录,每段时间称为SCAN,相邻两段时间有一定间隔,所以是分时;每段时间的采集参数可能变化也可能不变,所以是变参数。
在上述技术方案的基础上,本发明还可以做如下改进。
优选地,在上述任意实施例中,步骤1具体为:
通过高稳频率源生成基准信号Sref,基准信号Sref的计算公式为:
Sref=exp(j·2π·fref·t)
其中,fref为基准信号频率,t为信号时间,j为虚数单位。
优选地,在上述任意实施例中,步骤2具体为:
获取VLBI记录系统的记录数据Sn(t)的公式为:
Sn(t)=exp[j·2π·(fref-fLn)·tM]
其中,n表示第n个时间段,fLn为n个时间段的本振频率,
tM为该时间段内的采样点序列,Fn为第n个时间段的采样频率。
优选地,在上述任意实施例中,步骤3具体为:
对基准信号Sref以及所述记录数据进行互相关计算,得到互相关计算结果Scorn(t),公式为:
其中,conj()为复共轭计算,M为进行互相关计算的数据的点数,Sn(ti)为第n个预设时间段的VLBI记录系统的记录数据。
需要说明的是,每个预设时间段内可能包含多个Scorn(t),若包含多个Scorn(t)则需要对每个Scorn(t)进行计算得出结果。
优选地,在上述任意实施例中,步骤4具体为:
根据互相关计算结果Scorn(t)进行残余相位序列计算,得到残余相位序列的计算结果Φn(t),计算公式为:
Φn(t)=angle{Scorn(t)}
其中,angle{}为取幅角计算。
优选地,在上述任意实施例中,步骤5具体为:
根据残余相位序列的计算结果Φn(t)计算相位变化率均值φref,公式为:
其中,N为全部N个时间段,R为每个时间段内互相关计算得到的残余相位数,ti为第i个时间段的起始时刻,Φn(ti+1)为第n个预设时间内的第i+1个残余相位,Φn(ti)为第n个预设时间内的第i个残余相位。
优选地,在上述任意实施例中,步骤6具体为:
Φref=φref·(tn+1,1-tn,M)
其中,tn+1,1为第n+1个时间段的第一个点的时间,tn,M为第n个时间段的最后一个点的时间,Φn+1(t1)为第n+1个预设时间内的第一个残余相位,Φn(tM)为第n个预设时间内的第M个残余相位。
优选地,在上述任意实施例中,步骤7具体为:
一种存储介质,所述存储介质中存储有指令,当计算机读取所述指令时,使所述计算机执行如上述任一项所述的一种基于VLBI记录系统的初相稳定性判定方法。
在一些可能的实施方式中,通过利用高稳频率源产生基准信号Sref,VLBI记录系统分时段、变参数采集记录数据Sn(t),对记录数据Sn(t)与基准信号Sref进行互相关计算,获取全部采集序列的残余相位信息,通过拟合残余相位,判断系统初相稳定性指标,可以实现对单站记录设备的初相绝对稳定性的评估,还可以有效抑制系统误差,并在深空测控干涉测量领域具有显著的应用前景。
一种电子设备,包括存储器、处理器及存储在所述存储器上并在所述处理器上运行的程序,所述处理器执行所述程序时实现如上述任一项所述的一种基于VLBI记录系统的初相稳定性判定方法的步骤。
在一些可能的实施方式中,通过利用高稳频率源产生基准信号Sref,VLBI记录系统分时段、变参数采集记录数据Sn(t),对记录数据Sn(t)与基准信号Sref进行互相关计算,获取全部采集序列的残余相位信息,通过拟合残余相位,判断系统初相稳定性指标,可以实现对单站记录设备的初相绝对稳定性的评估,还可以有效抑制系统误差,并在深空测控干涉测量领域具有显著的应用前景。。
读者应理解,在本说明书的描述中,参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中,对上述术语的示意性表述不必针对的是相同的实施例或示例。而且,描述的具体特征、结构、材料或者特点可以在任一个或多个实施例或示例中以合适的方式结合。此外,在不相互矛盾的情况下,本领域的技术人员可以将本说明书中描述的不同实施例或示例以及不同实施例或示例的特征进行结合和组合。
在本申请所提供的几个实施例中,应该理解到,所揭露的装置和方法,可以通过其它的方式实现。例如,以上所描述的方法实施例仅仅是示意性的,例如,步骤的划分,仅仅为一种逻辑功能划分,实际实现时可以有另外的划分方式,例如多个步骤可以结合或者可以集成到另一个步骤,或一些特征可以忽略,或不执行。
上述方法如果以软件功能单元的形式实现并作为独立的产品销售或使用时,可以存储在一个计算机可读取存储介质中。基于这样的理解,本发明的技术方案本质上或者说对现有技术做出贡献的部分,或者该技术方案的全部或部分可以以软件产品的形式体现出来,该计算机软件产品存储在一个存储介质中,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行本发明各个实施例方法的全部或部分步骤。而前述的存储介质包括:U盘、移动硬盘、只读存储器(ROM,Read-OnlyMemory)、随机存取存储器(RAM,RandomAccessMemory)、磁碟或者光盘等各种可以存储程序代码的介质。
以上,仅为本发明的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到各种等效的修改或替换,这些修改或替换都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应以权利要求的保护范围为准。
Claims (3)
1.一种基于VLBI记录系统的初相稳定性判定方法,其特征在于,包括:
步骤1,生成基准信号Sref;
步骤2,按预设时间段连续获取VLBI记录系统的记录数据Sn(t);
步骤3,对所述基准信号Sref以及所述记录数据Sn(t)进行互相关计算;
步骤4,根据互相关计算结果Scorn(t)进行残余相位序列计算;
步骤5,根据残余相位序列的计算结果Φn(t)计算相位变化率均值φref;
其中,步骤1具体为:
通过高稳频率源生成基准信号Sref,基准信号Sref的计算公式为:
Sref=exp(j·2π·fref·t)
其中,fref为基准信号频率,t为信号时间,j为虚数单位;
步骤2具体为:
获取VLBI记录系统的记录数据Sn(t)的公式为:
Sn(t)=exp[j·2π·(fref-fLn)·tM]
其中,n表示第n个时间段,fLn为n个时间段的本振频率,
tM为该时间段内的采样点序列,Fn为第n个时间段的采样频率;
步骤3具体为:
对所述基准信号Sref以及所述记录数据进行互相关计算,得到互相关计算结果Scorn(t),公式为:
其中,conj()为复共轭计算,M为进行互相关计算的数据点数,Sn(ti)为第n个预设时间段的VLBI记录系统的记录数据;
步骤4具体为:
根据互相关计算结果Scorn(t)进行残余相位序列计算,得到残余相位序列的计算结果Φn(t),计算公式为:
Φn(t)=angle{Scorn(t)}
其中,angle{}为取幅角计算;
步骤5具体为:
根据残余相位序列的计算结果Φn(t)计算相位变化率均值φref,公式为:
其中,N为全部N个时间段,R为每个时间段内互相关计算得到的残余相位数,ti为第i个时间段的起始时刻,Φn(ti+1)为第n个预设时间内的第i+1个残余相位,Φn(ti)为第n个预设时间内的第i个残余相位;
步骤6具体为:
Φref=φref·(tn+1,1-tn,M)
其中,tn+1,1为第n+1个时间段的第一个点的时间,tn,M为第n个时间段的最后一个点的时间,Φn+1(t1)为第n+1个预设时间内的第一个残余相位,Φn(tM)为第n个预设时间内的第M个残余相位;步骤7具体为:
2.一种存储介质,其特征在于,所述存储介质中存储有指令,当计算机读取所述指令时,使所述计算机执行如权利要求1所述的一种基于VLBI记录系统的初相稳定性判定方法。
3.一种电子设备,包括存储器、处理器及存储在所述存储器上并在所述处理器上运行的程序,其特征在于,所述处理器执行所述程序时实现如权利要求1所述的一种基于VLBI记录系统的初相稳定性判定方法的步骤。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110456445.0A CN113237492B (zh) | 2021-04-27 | 2021-04-27 | 一种初相稳定性判定方法、存储介质及电子设备 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110456445.0A CN113237492B (zh) | 2021-04-27 | 2021-04-27 | 一种初相稳定性判定方法、存储介质及电子设备 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113237492A CN113237492A (zh) | 2021-08-10 |
CN113237492B true CN113237492B (zh) | 2022-05-31 |
Family
ID=77129427
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110456445.0A Active CN113237492B (zh) | 2021-04-27 | 2021-04-27 | 一种初相稳定性判定方法、存储介质及电子设备 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113237492B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115166613B (zh) * | 2022-06-10 | 2023-09-29 | 北京航天飞行控制中心 | 一种测站设备时延稳定性的检测方法、系统、介质及设备 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2004287734A (ja) * | 2003-03-20 | 2004-10-14 | Japan Research Institute Ltd | パラメータ推定方法、データ予測方法、パラメータ推定装置、データ予測装置、コンピュータプログラム、及び記録媒体 |
JP2010108402A (ja) * | 2008-10-31 | 2010-05-13 | Suwa Trust:Kk | 2系列の時系列データの解析装置及び2系列の時系列データの解析プログラムを記録したコンピュータ読み取り可能な記録媒体 |
CN103490856A (zh) * | 2013-09-16 | 2014-01-01 | 中国科学院上海天文台 | 用于深空探测的Mark5B格式VLBI数据接收解码与纠错系统及方法 |
CN105659869B (zh) * | 2011-09-27 | 2014-06-25 | 北京航天飞行控制中心 | 一种基于数传信号的甚长基线干涉测量处理方法 |
CN107766291A (zh) * | 2017-09-15 | 2018-03-06 | 中国人民解放军63920部队 | 一种获取甚长基线干涉测量中残余时延的方法 |
CN110823257A (zh) * | 2019-12-16 | 2020-02-21 | 中国人民解放军63921部队 | 深空多普勒测量数据时标偏差的自校准方法 |
CN111625413A (zh) * | 2020-04-23 | 2020-09-04 | 平安科技(深圳)有限公司 | 指标异常分析方法、装置及存储介质 |
CN112187262A (zh) * | 2020-06-16 | 2021-01-05 | 中国人民解放军战略支援部队航天工程大学 | 一种dbbc的同步采集时延控制方法 |
-
2021
- 2021-04-27 CN CN202110456445.0A patent/CN113237492B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2004287734A (ja) * | 2003-03-20 | 2004-10-14 | Japan Research Institute Ltd | パラメータ推定方法、データ予測方法、パラメータ推定装置、データ予測装置、コンピュータプログラム、及び記録媒体 |
JP2010108402A (ja) * | 2008-10-31 | 2010-05-13 | Suwa Trust:Kk | 2系列の時系列データの解析装置及び2系列の時系列データの解析プログラムを記録したコンピュータ読み取り可能な記録媒体 |
CN105659869B (zh) * | 2011-09-27 | 2014-06-25 | 北京航天飞行控制中心 | 一种基于数传信号的甚长基线干涉测量处理方法 |
CN103490856A (zh) * | 2013-09-16 | 2014-01-01 | 中国科学院上海天文台 | 用于深空探测的Mark5B格式VLBI数据接收解码与纠错系统及方法 |
CN107766291A (zh) * | 2017-09-15 | 2018-03-06 | 中国人民解放军63920部队 | 一种获取甚长基线干涉测量中残余时延的方法 |
CN110823257A (zh) * | 2019-12-16 | 2020-02-21 | 中国人民解放军63921部队 | 深空多普勒测量数据时标偏差的自校准方法 |
CN111625413A (zh) * | 2020-04-23 | 2020-09-04 | 平安科技(深圳)有限公司 | 指标异常分析方法、装置及存储介质 |
CN112187262A (zh) * | 2020-06-16 | 2021-01-05 | 中国人民解放军战略支援部队航天工程大学 | 一种dbbc的同步采集时延控制方法 |
Non-Patent Citations (1)
Title |
---|
中国VLBI网软件相关处理机测地应用精度分析;刘磊等;《测绘学报》;20170715(第07期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN113237492A (zh) | 2021-08-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US5792062A (en) | Method and apparatus for detecting nonlinearity in an electrocardiographic signal | |
KR100590561B1 (ko) | 신호의 피치를 평가하는 방법 및 장치 | |
Pham et al. | Identifying low variance pathways for free energy calculations of molecular transformations in solution phase | |
Rosenstein et al. | Reconstruction expansion as a geometry-based framework for choosing proper delay times | |
US8332163B2 (en) | Methods and apparatus for estimating rotor slots | |
CN113237492B (zh) | 一种初相稳定性判定方法、存储介质及电子设备 | |
JP2011507137A (ja) | サーボデータを用いたフライングハイト制御のためのシステム及び方法 | |
CN109946518B (zh) | 基于贝叶斯方法的电力谐波信号分析方法与分析设备 | |
US20100271724A1 (en) | Systems and Methods for Nyquist Tone Harmonic Measurements | |
US20110255700A1 (en) | Detecting Musical Structures | |
Henry | Spectral analysis techniques using Prism signal processing | |
Henry | An ultra-precise fast Fourier transform | |
WO2007072394A2 (en) | Audio structure analysis | |
US20080278837A1 (en) | Frequency, phase, and gain estimation technique for use in a read channel receiver of a hard disk drive | |
US11593360B2 (en) | Empirically providing data privacy with reduced noise | |
Martin et al. | Estimating covariant Lyapunov vectors from data | |
KR20030045068A (ko) | 산란 시그네춰를 이용한 구조 확인 개선 | |
Labuda | Daniell method for power spectral density estimation in atomic force microscopy | |
US20100191790A1 (en) | System and method for correlation scoring of signals | |
CN114343675B (zh) | 一种脑电成分提取方法 | |
Henry | An ultra-precise Fast Fourier Transform | |
JP6286933B2 (ja) | 小節間隔推定およびその推定のための特徴量抽出を行う装置、方法、およびプログラム | |
EP3768895A1 (en) | Asphalt density estimation system, and related method of reducing signal noise | |
CN111832148B (zh) | 波形间隙填充 | |
CN115452010A (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 |