CN106093645B - 确定微波部件宽功率变化范围高阶无源互调电平的方法 - Google Patents

确定微波部件宽功率变化范围高阶无源互调电平的方法 Download PDF

Info

Publication number
CN106093645B
CN106093645B CN201610437928.5A CN201610437928A CN106093645B CN 106093645 B CN106093645 B CN 106093645B CN 201610437928 A CN201610437928 A CN 201610437928A CN 106093645 B CN106093645 B CN 106093645B
Authority
CN
China
Prior art keywords
passive intermodulation
power
value
rank
higher order
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
CN201610437928.5A
Other languages
English (en)
Other versions
CN106093645A (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.)
Xian Institute of Space Radio Technology
Original Assignee
Xian Institute of Space Radio Technology
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 Xian Institute of Space Radio Technology filed Critical Xian Institute of Space Radio Technology
Priority to CN201610437928.5A priority Critical patent/CN106093645B/zh
Publication of CN106093645A publication Critical patent/CN106093645A/zh
Application granted granted Critical
Publication of CN106093645B publication Critical patent/CN106093645B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R31/00Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
    • G01R31/001Measuring interference from external sources to, or emission from, the device under test, e.g. EMC, EMI, EMP or ESD testing

Abstract

本发明涉及确定微波部件宽功率变化范围高阶无源互调电平的方法,首先对微波部件3阶无源互调电平的测量值对载波功率之和进行多项式拟合,基于拟合后的3阶无源互调电平随载波功率之和的变化曲线,以载波功率之和的最小值向更小取整为起始功率,以载波功率之和的最大值向更大取整为终点功率,以待预测高阶无源互调的载波功率为中心,依次采用奇数个功率及互调电平测量值对3阶无源互调电平随功率变化关系的幂级数展开式进行参数估计,从而实现对高阶无源互调电平的预测,该方法实现微波部件宽功率变化范围高阶无源互调的准确预测,为后续型号在高阶无源互调测试系统不具备条件下微波部件高阶无源互调电平的评估提供了有效手段。

Description

确定微波部件宽功率变化范围高阶无源互调电平的方法
技术领域
本发明涉及一种确定微波部件宽功率变化范围高阶无源互调电平的方法,属于微波部件无源互调领域。
背景技术
无源互调(Passive-Intermodulation,简称PIM)是指在大功率条件下,当输入两个或者两个以上载波时,由于微波无源部件的非线性导致载波信号相互调制,产生载波频率的组合产物落入接收通带内造成干扰的现象,当互调电平较低时,会使接收信号底噪抬高,使接收机信噪比降低,误码率升高;当互调电平进一步增高时,会影响整个通信系统的正常工作,被迫降低功率使用,或分通道使用;严重时互调产物将淹没接收信号,导致通道阻塞,通信中断,使整个系统处于瘫痪状态。因此无源互调效应是影响航天器载荷性能、高可靠性的重要因素。
在空间应用中,受频率资源的限制,发射信号和接收带宽通常相距较远,无源互调的阶数通常较高。而高阶的无源互调的功率电平更小,对小信号的观测较困难导致高阶的无源互调测量困难较大,而测试低阶的无源互调相对容易。因此在进行微波部件无源互调分析评价时,通常需要根据低阶无源互调来预测高阶无源互调电平,从而实现对微波部件无源互调特性的评价。
而低阶无源互调电平测量值随功率增长在较宽功率范围内并非线性关系,而传统的预测高阶无源互调电平的方法要么假设低阶无源互调电平测量值随功率增长为线性关系,要么采用所有预测值进行高阶多项式拟合,存在电平奇异值,预测误差大,无法满足微波部件无源互调性能评估的要求,因此需要实现宽功率变化范围微波部件高阶无源互调电平的准确预测。
发明内容
本发明的目的在于克服现有技术的上述缺陷,提供一种确定微波部件宽功率变化范围高阶无源互调电平的方法,该方法能够实现微波部件高阶无源互调的准确预测,为后续型号在高阶无源互调测试系统不具备条件下微波部件无源互调水平的评估提供了有效手段。
本发明的上述目的主要是通过如下技术方案予以实现的:
确定微波部件宽功率变化范围高阶无源互调电平的方法,包括如下步骤:
步骤(1)、预先给定N个功率的3阶无源互调电平测量值Pmea3(i),i=1…N,其中无源互调测量时输入的两路载波信号功率相等,Ps(i)为第i个无源互调测量输入两路载波信号的总功率,;
步骤(2)、设定功率范围内的无源互调功率点总数为NN,第j个无源互调测量输入的单路载波信号功率为P'c(j),则:
NN=Pmax-Pmin+1;
P'c(j)=P's(j)-3,j=1…NN;
其中:Pmin为min{Ps(i)}以10向小取整;
Pmax为max{Ps(i)}以10向大取整;
P's(j)=Pmin+j-1;
步骤(3)、采用最小二乘法获得最小时对应的系数p1,p2,p3,其中:P3f(i)的表达式如下:
P3f(i)=p1+p2*Ps(i)+p3*(Ps(i))2*log(Ps(i))i=1……N;
步骤(4)、根据系数p1,p2,p3得到3阶无源互调电平测量值Pmea3(i)的修正值P'mea3(j);
步骤(5)、从NN个功率点中第1个功率点开始选取奇数M个点,获得M个功率点的中心点的5阶无源互调值PIM5(j)、7阶无源互调值PIM7(j)和9阶无源互调值PIM9(j);从NN个功率点中第2个功率点开始选取奇数M个点,获得M个功率点的中心点的5阶无源互调值PIM5(j)、7阶无源互调值PIM7(j)和9阶无源互调值PIM9(j);依次类推,共获得NN-M+1个功率点的5阶无源互调值PIM5(j)、7阶无源互调值PIM7(j)和9阶无源互调值PIM9(j)。
在上述确定微波部件宽功率变化范围高阶无源互调电平的方法中,所述步骤(4)中根据系数p1,p2,p3,通过如下公式得到3阶无源互调电平测量值Pmea3(i)的修正值P'mea3(j):
P'mea3(j)=p1+p2*Ps'(j)+p3*(Ps'(j))2*log(Ps'(j))j=1…NN。
在上述确定微波部件宽功率变化范围高阶无源互调电平的方法中,所述步骤(5)中获得NN-M+1个功率点的5阶无源互调值PIM5(j)、7阶无源互调值PIM7(j)和9阶无源互调值PIM9(j)的具体方法如下:
(1)、令ii=0;
(2)、令jj=ii+1,取M个功率点P'c(jj),P'c(jj+1),……P'c(jj+M-1)及其对应的无源互调值,采用最小二乘法获得最小时对应的系数a3,a5,a7,a9,a11,其中PIM3(j)的表达式如下:
(3)、采用步骤(2)得到的系数,根据如下公式获得M个功率点的中心点的5阶无源互调值PIM5(j)、7阶无源互调值PIM7(j)和9阶无源互调值PIM9(j):
(4)、判断jj>NN-M+1,若否,令ii加1,返回步骤(2);若是,进入步骤(5);
(5)、结束。
在上述确定微波部件宽功率变化范围高阶无源互调电平的方法中,所述获得的NN-M+1个功率点的5阶无源互调值PIM5(j)、7阶无源互调值PIM7(j)和9阶无源互调值PIM9(j),其对应的载波信号的总功率为
本发明与现有技术相比具有如下有益效果:
(1)、本发明给出一种新型的确定微波部件宽功率变化范围高阶无源互调电平的方法,首先对微波部件3阶无源互调电平的测量值对载波功率之和进行多项式拟合,基于拟合后的3阶无源互调电平随载波功率之和的变化曲线,以载波功率之和的最小值向更小取整为起始功率,以载波功率之和的最大值向更大取整为终点功率,以待预测高阶无源互调的载波功率为中心,依次采用奇数个功率及互调电平测量值对3阶无源互调电平随功率变化关系的幂级数展开式进行参数估计,采用获得的参数结合高阶无源互调表达式进行高阶无源互调的计算,从而实现对高阶无源互调电平的预测,;
(2)、本发明方法针对无源互调在宽功率变化范围内随载波功率呈非线性关系的问题,很好的解决了微波部件高阶无源互调预测的难题,计算结果与高阶测量结果吻合良好,预测方法具有较高的准确性;
(3)、本发明方法对微波部件宽功率变化范围高阶无源互调进行了准确预测,可以代替实际测试得到准确的高阶无源互调值,省去了测试的繁琐程序,降低成本,为后续型号在高阶无源互调测试系统不具备条件下微波部件无源互调水平的评估提供了有效手段。
附图说明
图1为本发明确定微波部件宽功率变化范围高阶无源互调电平的方法流程图;
图2为本发明实施例中的基于3阶无源互调电平预测的5、7、9阶无源互调的结果与测量值的对比图。
具体实施方式
下面结合附图和具体实施例对本发明作进一步详细的描述:
如图1所示为本发明确定微波部件宽功率变化范围高阶无源互调电平的方法流程图,本发明微波部件宽功率变化范围高阶无源互调电平的方法具体包括如下步骤:
步骤(1)、预先给定N个功率的3阶无源互调电平测量值Pmea3(i),i=1…N,其中无源互调测量时输入的两路载波信号功率相等,Ps(i)为第i个无源互调测量输入两路载波信号的总功率;功率单位为dBm。
步骤(2)、设定功率变化范围内的无源互调功率点总数为NN,第j个无源互调测量输入的单路载波信号功率为P'c(j),则:
NN=Pmax-Pmin+1;
P'c(j)=P's(j)-3,j=1…NN;功率单位为dBm;
其中:Pmin为min{Ps(i)}以10向小取整,
Pmax为max{Ps(i)}以10向大取整;
P's(j)=Pmin+j-1;
例如10~19以10向小取整为10;10~19以10向大取整为20。再例如20~29以10向小取整为20;20~29以10向大取整为30,依次类推。
步骤(3)、采用最小二乘法获得最小时对应的系数p1,p2,p3,其中:P3f(i)的表达式如下:
P3f(i)=p1+p2*Ps(i)+p3*(Ps(i))2*log(Ps(i))i=1……N;
步骤(4)、根据系数p1,p2,p3,通过如下公式得到3阶无源互调电平测量值Pmea3(i)的修正值P'mea3(j):
P'mea3(j)=p1+p2*Ps'(j)+p3*(Ps'(j))2*log(Ps'(j))j=1…NN。
步骤(5)、从NN个功率点中第1个功率点开始选取奇数M个点,获得M个功率点的中心点的5阶无源互调值PIM5(j)、7阶无源互调值PIM7(j)和9阶无源互调值PIM9(j);从NN个功率点中第2个功率点开始选取奇数M个点,获得M个功率点的中心点的5阶无源互调值PIM5(j)、7阶无源互调值PIM7(j)和9阶无源互调值PIM9(j);依次类推,通过循环的方法,共获得NN-M+1个功率点的5阶无源互调值PIM5(j)、7阶无源互调值PIM7(j)和9阶无源互调值PIM9(j)。
获得NN-M+1个功率点的5阶无源互调PIM5(j)、7阶无源互调PIM7(j)和9阶无源互调PIM9(j)的具体方法如下:
(5.1)、令ii=0;
(5.2)、令jj=ii+1,取M个功率点P'c(jj),P'c(jj+1),……P'c(jj+M-1)及其对应的无源互调值,采用最小二乘法获得最小时对应的系数a3,a5,a7,a9,a11,其中PIM3(j)的表达式如下:
其中:PIM3(j)为3阶无源互调随功率变化关系的幂级数展开式。
(5.3)、采用步骤(2)得到的系数,根据如下公式获得M个功率点的中心点的5阶无源互调值PIM5(j)、7阶无源互调值PIM7(j)和9阶无源互调值PIM9(j):
(5.4)、判断jj>NN-M+1,若否,令ii加1(即令ii的值增加1),返回步骤(5.2);若是,进入步骤(5.5);
(5.5)、结束。
获得的NN-M+1个功率点的5阶无源互调值PIM5(j)、7阶无源互调值PIM7(j)和9阶无源互调值PIM9(j),即第个功率点的无源互调值,其对应的载波信号的总功率为P's(j),
实施例
(1)、N=47;给定3阶的无源互调测量值,第一列为Ps(i),第二列为对应的3阶无源互调测量电平Pmea3(i),i=1…N。
(2)、设定功率变化范围内的无源互调功率点总数NN=31,Pmin=20,Pmax=50,P’s(j)=j+19,P’c(j)=j+16,j=1…NN。
(3)、采用最小二乘法获得最小时对应的系数p1,p2,p3,其中:P3f(i)的表达式如下:
P3f(i)=p1+p2*Ps(i)+p3*(Ps(i))2*log(Ps(i))i=1……N;
得到:
p1=-186.097323058506;
p2=3.66221275703211;
p3=-0.00675616299861527;
(4)、根据系数p1,p2,p3,通过如下公式得到3阶无源互调电平测量值Pmea3(i)的修正值P'mea3(j):
P'mea3(j)=p1+p2*Ps'(j)+p3*(Ps'(j))2*log(Ps'(j))j=1…NN。
得到:
(5)、M=7,采取上述步骤(5.1)~(5.5)的方法,获得共25个(NN-M+1个)5阶、7阶和9阶的无源互调预测值。
首先第1个循环从31个功率点数中选取序号为第1~第7的7个点,获得7个功率点的中心点,即第4个点的5阶、7阶和9阶的无源互调预测值;接着第2个循环从31个功率点数中选取序号为第2~第8的7个点,获得7个功率点的中心点,即第5个点的5阶、7阶和9阶的无源互调预测值;接着第3个循环从31个功率点数中选取序号为第3~第9的7个点,获得7个功率点的中心点,即第6个点的5阶、7阶和9阶的无源互调预测值;依次类推,直至获得最后1个循环中7个功率点的中心点,即第28个点的5阶、7阶和9阶的无源互调预测值。
(6)、25个5阶、7阶和9阶的无源互调预测值无源互调预测值如下:
如图2所示为本发明实施例中的基于3阶无源互调电平预测的5、7、9阶无源互调的结果与测量值的对比图,由图可知,本发明方法的计算结果与高阶测量结果吻合良好,本发明预测方法具有较高的准确性。
以上所述,仅为本发明最佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。
本发明说明书中未作详细描述的内容属于本领域专业技术人员的公知技术。

Claims (4)

1.确定微波部件宽功率变化范围高阶无源互调电平的方法,其特征在于:包括如下步骤:
步骤(1)、预先给定N个功率的3阶无源互调电平测量值Pmea3(i),i=1…N,其中无源互调测量时输入的两路载波信号功率相等,Ps(i)为第i个无源互调测量输入两路载波信号的总功率;
步骤(2)、设定功率范围内的无源互调功率点总数为NN,第j个无源互调测量输入的单路载波信号功率为P'c(j),则:
NN=Pmax-Pmin+1;
P'c(j)=P's(j)-3,j=1…NN;
其中:Pmin为min{Ps(i)}以10向小取整;
Pmax为max{Ps(i)}以10向大取整;
P's(j)=Pmin+j-1;
步骤(3)、采用最小二乘法获得最小时对应的系数p1,p2,p3,其中:P3f(i)的表达式如下:
P3f(i)=p1+p2*Ps(i)+p3*(Ps(i))2*log(Ps(i))i=1……N;
步骤(4)、根据系数p1,p2,p3得到3阶无源互调电平测量值Pmea3(i)的修正值P'mea3(j);
步骤(5)、从NN个功率点中第1个功率点开始选取奇数M个点,获得M个功率点的中心点的5阶无源互调值PIM5(j)、7阶无源互调值PIM7(j)和9阶无源互调值PIM9(j);从NN个功率点中第2个功率点开始选取奇数M个点,获得M个功率点的中心点的5阶无源互调值PIM5(j)、7阶无源互调值PIM7(j)和9阶无源互调值PIM9(j);依次类推,共获得NN-M+1个功率点的5阶无源互调值PIM5(j)、7阶无源互调值PIM7(j)和9阶无源互调值PIM9(j)。
2.根据权利要求1所述的确定微波部件宽功率变化范围高阶无源互调电平的方法,其特征在于:所述步骤(4)中根据系数p1,p2,p3,通过如下公式得到3阶无源互调电平测量值Pmea3(i)的修正值P'mea3(j):
P'mea3(j)=p1+p2*Ps'(j)+p3*(Ps'(j))2*log(Ps'(j))j=1…NN。
3.根据权利要求1所述的确定微波部件宽功率变化范围高阶无源互调电平的方法,其特征在于:所述步骤(5)中获得NN-M+1个功率点的5阶无源互调值PIM5(j)、7阶无源互调值PIM7(j)和9阶无源互调值PIM9(j)的具体方法如下:
(1)、令ii=0;
(2)、令jj=ii+1,取M个功率点P'c(jj),P'c(jj+1),……P'c(jj+M-1)及其对应的无源互调值,采用最小二乘法获得最小时对应的系数a3,a5,a7,a9,a11,其中PIM3(j)的表达式如下:
其中:PIM3(j)为3阶无源互调值;
(3)、采用步骤(2)得到的系数,根据如下公式获得M个功率点的中心点的5阶无源互调值PIM5(j)、7阶无源互调值PIM7(j)和9阶无源互调值PIM9(j):
(4)、判断jj>NN-M+1,若否,令ii加1,返回步骤(2);若是,进入步骤(5);
(5)、结束。
4.根据权利要求1或3所述的确定微波部件宽功率变化范围高阶无源互调电平的方法,其特征在于:所述获得的NN-M+1个功率点的5阶无源互调值PIM5(j)、7阶无源互调值PIM7(j)和9阶无源互调值PIM9(j),其对应的载波信号的总功率为P's(j),
CN201610437928.5A 2016-06-17 2016-06-17 确定微波部件宽功率变化范围高阶无源互调电平的方法 Active CN106093645B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610437928.5A CN106093645B (zh) 2016-06-17 2016-06-17 确定微波部件宽功率变化范围高阶无源互调电平的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610437928.5A CN106093645B (zh) 2016-06-17 2016-06-17 确定微波部件宽功率变化范围高阶无源互调电平的方法

Publications (2)

Publication Number Publication Date
CN106093645A CN106093645A (zh) 2016-11-09
CN106093645B true CN106093645B (zh) 2019-01-11

Family

ID=57236332

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610437928.5A Active CN106093645B (zh) 2016-06-17 2016-06-17 确定微波部件宽功率变化范围高阶无源互调电平的方法

Country Status (1)

Country Link
CN (1) CN106093645B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107154826B (zh) * 2017-04-28 2020-11-10 西安空间无线电技术研究所 一种确定大功率微波部件高阶无源互调电平的方法
CN109359423B (zh) * 2018-11-16 2020-10-23 北京理工大学 一种快速确定矩形波导法兰3阶无源互调产物电平的方法
CN110716095A (zh) * 2019-09-05 2020-01-21 北京航空航天大学 无源互调非线性特性参数的多状态测量方法和装置
CN115604744A (zh) * 2021-06-28 2023-01-13 华为技术有限公司(Cn) 一种无源互调源数目确定方法及相关设备

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102025432A (zh) * 2010-12-17 2011-04-20 镇江市澳华测控技术有限公司 隔离器正反向无源互调测试装置及方法
CN102544680A (zh) * 2012-01-19 2012-07-04 福建三元达通讯股份有限公司 一种改善微波器件无源互调的方法
CN103368665A (zh) * 2012-03-30 2013-10-23 中国联合网络通信有限公司广东省分公司 多信号无源互调测试方法、设备及系统
CN104104455A (zh) * 2013-04-09 2014-10-15 中兴通讯股份有限公司 无源互调位置检测方法及装置
CN105069206A (zh) * 2015-07-24 2015-11-18 西安空间无线电技术研究所 一种确定微波部件的无源互调产物的方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0486573A (ja) * 1990-07-30 1992-03-19 Uchu Tsushin Kiso Gijutsu Kenkyusho:Kk 簡易多周波混変調測定装置

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102025432A (zh) * 2010-12-17 2011-04-20 镇江市澳华测控技术有限公司 隔离器正反向无源互调测试装置及方法
CN102544680A (zh) * 2012-01-19 2012-07-04 福建三元达通讯股份有限公司 一种改善微波器件无源互调的方法
CN103368665A (zh) * 2012-03-30 2013-10-23 中国联合网络通信有限公司广东省分公司 多信号无源互调测试方法、设备及系统
CN104104455A (zh) * 2013-04-09 2014-10-15 中兴通讯股份有限公司 无源互调位置检测方法及装置
CN105069206A (zh) * 2015-07-24 2015-11-18 西安空间无线电技术研究所 一种确定微波部件的无源互调产物的方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
《Prediction of passive intermodulation between rough waveguide flanges based on fractal theory》;Chunjiang Bai;《Radio Science Meeting (Joint with AP-S Symposium), 2015 USNC-URSI》;20151026;全文
《微波与射频频段无源互调干扰研究》;张世全;《中国博士学位论文全文数据库 基础科学辑》;20050615;全文

Also Published As

Publication number Publication date
CN106093645A (zh) 2016-11-09

Similar Documents

Publication Publication Date Title
CN106093645B (zh) 确定微波部件宽功率变化范围高阶无源互调电平的方法
CN104698273A (zh) 谐波责任划分方法和系统
CN111478737A (zh) 射频测试方法、装置、电子装置及存储介质
US20150036729A1 (en) Method for testing sensitivity of a data packet signal transceiver
JP2016519376A5 (zh)
CN104333492A (zh) 一种评估通信网络结构合理性的方法及装置
CN103837767A (zh) 一种利用矢网对互易混频器进行特性表征的方法
CN103472333A (zh) 风电并网电能质量综合性能检测方法
WO2019205003A1 (en) Apparatus and method for transmission system
CN105429719A (zh) 基于功率谱和多尺度小波变换分析强干扰信号检测方法
CN103475431A (zh) 一种低信噪比条件下的基于最佳随机共振的频谱感知方法
CN106597392B (zh) 一种sar系统相位误差补偿方法
CN109959825A (zh) 一种基于拉格朗日插值法的非线性压控衰减器拟合校准方法
CN108494508A (zh) 基于mwc相关支撑集恢复的高效频谱检测方法
CN107154826A (zh) 一种确定大功率微波部件高阶无源互调电平的方法
CN102932077B (zh) 一种频谱占用度的测量方法
US20150033844A1 (en) Frequency modulation continuous wave radar level meter and measuring method for the same
CN106162681B (zh) 一种gsm900二次谐波干扰的检测方法及装置
CN113805156B (zh) 一种低信噪比的信号还原方法及系统
CN105445549B (zh) 一种微弱电磁信号频谱测量方法
CN200994119Y (zh) 一种栅压自适应补偿装置
CN106487461A (zh) 一种驻波峰值确定方法和设备
KR20170079842A (ko) 중계기 및 이의 신호 감쇄 방법
CN109525333B (zh) 天线辐射性能的测试方法和测试系统
CN103384374B (zh) 一种语音业务质量的评估方法和设备

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant