CN109787676B - 一种高动态下的零陷展宽方法 - Google Patents
一种高动态下的零陷展宽方法 Download PDFInfo
- Publication number
- CN109787676B CN109787676B CN201811607382.9A CN201811607382A CN109787676B CN 109787676 B CN109787676 B CN 109787676B CN 201811607382 A CN201811607382 A CN 201811607382A CN 109787676 B CN109787676 B CN 109787676B
- Authority
- CN
- China
- Prior art keywords
- covariance matrix
- updated
- data
- rotation
- time
- 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
Images
Landscapes
- Radio Transmission System (AREA)
Abstract
本发明一种高动态下的零陷展宽方法,包括:获取阵列天线接收信号采样数据,进行带通滤波、Hilbert变换以及幅相校正对所获取采样数据等间隔抽取后分两路输出,其中一路输入FIFO模块存储,及另一路进行等间隔抽取后送入QR分解单元进行协方差矩阵估计,输出得到更新的协方差矩阵;进行中间变量的前向迭代更新,根据更新的协方差矩阵和中间变量进行后向迭代计算,得到更新的最优权值;将FIFO模块所存储估计协方差矩阵的采样数据,并与更新的最优权值进行乘加运算后滤波输出。本发明等效扩展了干扰零陷,增强了滤波的稳定性,并进行了FPGA实现,增强了滤波的稳定性。
Description
技术领域
本发明属于卫星通信技术领域,涉及一种高动态下的零陷展宽方法。
背景技术
随着卫星通信技术的发展与应用,其接收机面临的窄带或宽带压制式干扰问题也日益突出。基于阵列信号处理技术的接收机可以通过调整各接收阵元的加权复系数来自适应地消除这些干扰。工程实现中,较常采用LMS算法及SMI算法进行自适应滤波,LMS算法实现简单且计算量小,较早的进行了工程实现,但算法的收敛性能受相关矩阵的特征值分布影响较大,因此抑制强干扰的能力有限。SMI算法是开环结构,并不受特征值分散影响,具有较好的干扰抑制效果,目前工程实现中更多的采用SMI算。
当卫星通信系统应用于高动态场景时,SMI的滤波性能会严重降低,基于高动态环境的抗干扰研究也逐渐受到重视。文献《卫星导航系统中平台运动对天线阵列性能的影响分析》分析了平台运动对自适应阵列抗干扰性能的影响,文献《高动态条件下统计空时零陷加宽方法》根据弹载导航接收机的特点,提出了一种基于统计的空时零陷加宽方法,但是该方法需要在原协方差矩阵的基础上点乘一个扩张矩阵,实现过程复杂。文献《一种小快拍时对角加载的快速算法》将对角加载与数据矩阵的QR分解结合起来,进行权矢量的实时更新。
发明内容
发明所要解决的课题是,克服现有技术的不足,在QRD-SMI算法基础上,结合高动态环境的特点,提出抽取式QRD-SMI滤波方法,等效扩展了干扰零陷,增强了滤波的稳定性,方法实现过程简单,利于FPGA实现。
QRD-SMI算法在传统SMI算法基础上通过对接收端的数据阵进行QR分解估计协方差矩阵,避免了直接矩阵求逆,进而简化了最优权求解过程,可对数据矩阵X进行QR分解,假设存在L×L的酉矩阵Q满足QXH=R,其中R为M×M维的非奇异上三角阵。由此可得:XXH=(RHQ)(QHR)=RHR,最优权估计则通过下式中的前后向代入得到:
本发明在QRD-SMI算法基础上,用等间隔抽取式数据块代替连续数据块估计协方差矩阵。
为了解决上述课题,本发明提出一种高动态下的零陷展宽方法,包括以下步骤:
步骤1、获取AD采集的阵列天线信号数据,对经阵列天线信号数据进行进行带通滤波、Hilbert变换以及幅相校正后进入步骤2;
步骤2、对采样数据等间隔抽取后分两路输出,其中一路输入FIFO模块存储,及另一路进行等间隔抽取后送入QR分解单元进行协方差矩阵估计,输出得到更新的协方差矩阵R;
步骤3、根据更新的协方差矩阵R代入式RHV=La(θd)中进行中间变量V的前向迭代更新,得到更新的中间变量V,其中a(θd)为角度θd对应的导向矢量;L为采样次数;
步骤4、根据更新的协方差矩阵R和中间变量V进行后向迭代计算,得到更新的最优权值Wopt;
步骤5、将FIFO模块所存储采样数据,并与更新的最优权值Wopt进行乘加运算后滤波输出。
进一步地,作为本发明的一种优选技术方案,具体为:
每间隔T时刻抽取一次采样快拍,并连续操作N次;
获取第i组NT时段快拍的入射角θi,得到接收信号模型为,及其对应的协方差矩,根据该协方差矩阵矩阵求逆引理,得到加权矢量的最优近似解。
进一步地,作为本发明的一种优选技术方案,所述步骤2中QR分解单元进行QR分解,由N次Givens旋转进行输入数据消元实现,抽取间隔Td大于一次Givens旋转时间,协方差估计样本数N取值大于2M,M为接收阵元数。
进一步地,作为本发明的一种优选技术方案,所述步骤2中QR分解单元进行QR分解具体为:
将所获取估计协方差矩阵的采样数据经过延迟处理分别依次延迟一个脉动间隔t,在第ti个脉动时刻,第一排数据x1(t1)进入脉动结构,经过边界单元计算后得出旋转因子c1/s1,该边界单元计算对角元素rii和产生内部单元所需要的Givens的旋转因子c和s,其中使用两次CORDIC旋转,第一次使用CORDIC旋转来实现输入数据的求模运算得到输入数据的模值|xin|和旋转角θ,第二次使用CORDIC旋转来实现对角元素rii和旋转角更新;计算旋转角θ和的正余弦值,更新得到旋转因子c和s;
在第ti+1个脉动时刻,第二排数据x2(t1)进入脉动结构,经过内部单元的一次计算,得出经过一次Givens的旋转因子c1/s1的上三角矩阵R的非对角元素r12和经过变换后待消元元素xout12;
在第ti+2个脉动时刻,第三排数据x3(t1)经过内部单元的一次计算,得出经过一次Givens的旋转因子c1/s1的上三角矩阵R的非对角元素r13和经过变换后待消元元素xout13,依次类推,每经过一次边界单元消元一个元素,采用流水线结构经过7个脉动周期便可以将新进的元素消元,完成更新后的上三角矩阵,及在t=ti+7时刻输出更新的协方差矩阵R。
进一步地,作为本发明的一种优选技术方案所述步骤4得到更新的最优权值Wopt采用公式:RWopt=V。
发明效果
本发明针对高动态环境下阵列信号接收建立了近似模型,为高动态情况下阵列信号处理提供了理论基础,并利用近似模型从理论上推导验证了本发明的正确性。
本发明通过等间隔抽取接收数据作为协方差矩阵的估计基本元素,相同的采样数据量包括更多的方位信息,高动态环境下具有更精确的协方差估计,可以达到零陷展宽的效果。本发明中以等间隔抽取代替传统连续抽取,抵消了协方差矩阵的更新延迟,可以实现信号的实时滤波。
本发明的零陷展宽方法进行FPGA实现时,只需在传统QRD-SMI结构基础上增加一个FIFO实现数据暂存,无需增加额外的硬件结构,相比其他零陷展宽方法实现简单,易于与高性能FPGA相结合。
因此,本发明在QRD-SMI算法基础上,结合高动态环境的特点,提出抽取式QRD-SMI滤波方法,等效扩展了干扰零陷,增强了滤波的稳定性,并进行了FPGA实现。
附图说明
图1为本发明方法的流程图;
图2为本发明等间隔抽取GR-TSA结构示意图;
图3为本发明边界单元FPGA实现结构图;
图4为本发明内部单元FPGA实现结构图;
图5为本发明中间变量V的FPGA实现结构图;
图6为本发明权值W的FPGA实现结构图。
具体实施方式
以下,基于附图针对本发明进行详细地说明。
如图1所示,本发明设计了一种高动态下的零陷展宽方法,该方法是等间隔抽取式的改进QRD-SMI方法,高动态环境下,可以有效增强滤波稳健性。所述改进QRD-SMI方法包括输入数据延迟寄存,数据抽取,QR分解和实时滤波,具体包括以下步骤:
步骤1、获取接收信号的估计协方差矩阵的采样数据,对估计协方差矩阵的采样数据进行带通滤波、Hilbert变换以及幅相校正后进入步骤2;
步骤2、对所获取估计协方差矩阵的采样数据等间隔抽取后分两路输出,其中一路输入FIFO模块存储,等待权值更新后进行滤波处理;及另一路进行等间隔抽取后送入QR分解单元进行协方差矩阵估计,输出得到更新的协方差矩阵R,其中QR分解由N次Givens旋转进行输入数据消元实现,抽取间隔Td需大于一次Givens旋转时间,协方差估计样本数N取值应大于2M,M为接收阵元数。
首先,本方法进行等间隔抽取过程为:
假设干扰入射角度随时间连续变化,对估计协方差矩阵Rxx的采样数据做如下处理,每间隔T时刻抽取一次采样快拍,并连续操作N次。当NT持续时间较短时,N次快拍的入射角可近似认为相等,这里用θ1表示。在下一个NT时段等间隔抽取N次快拍,其入射角用θ2表示。依次类推,第i组NT时段的入射角则看做θi。由此接收信号模型可近似修改为:其对应的协方差矩阵可表示为:其中σ2为高斯白噪声的功率;AJ=[a(θ1) a(θ2) … a(θL)];RJ=diag(δ1,δ1,…δL)为干扰信号功率。最终得到的加权矢量的最优解近似为:
然后,等间隔抽取送入QR分解单元进行协方差矩阵估计,进行QR分解滤波的具体操作步骤如下具体如下:
将所获取估计协方差矩阵的采样数据经过延迟处理分别依次延迟一个脉动间隔t,在第ti个脉动时刻,第一排数据x1(t1)进入脉动结构,经过边界单元计算后得出旋转因子c1/s1,该边界单元计算对角元素rii和产生内部单元所需要的Givens的旋转因子c和s,其中使用两次CORDIC旋转,第一次使用CORDIC旋转来实现输入数据的求模运算得到输入数据的模值|xin|和旋转角θ,第二次使用CORDIC旋转来实现对角元素rii和旋转角更新;计算旋转角θ和的正余弦值,更新得到旋转因子c和s;
在第ti+1个脉动时刻,第二排数据x2(t1)进入脉动结构,经过内部单元的一次计算,得出经过一次Givens的旋转因子c1/s1的上三角矩阵R的非对角元素r12和经过变换后待消元元素xout12;
在第ti+2个脉动时刻,第三排数据x3(t1)经过内部单元的一次计算,得出经过一次Givens的旋转因子c1/s1的上三角矩阵R的非对角元素r13和经过变换后待消元元素xout13,依次类推,每经过一次边界单元消元一个元素,采用流水线结构经过7个脉动周期便可以将新进的元素消元,完成更新后的上三角矩阵,及在t=ti+7时刻输出更新的协方差矩阵R。
步骤3、N次旋转后根据更新的协方差矩阵R代入式RHV=La(θd)中进行中间变量V的前向迭代更新,得到更新的中间变量V,其中a(θd)为角度θd对应的导向矢量;L为采样次数。
步骤4、根据更新的协方差矩阵R和中间变量V进行后向迭代计算,采用公式RWopt=V得到更新的最优权值Wopt;
步骤5、将FIFO模块所存储估计协方差矩阵的采样数据,并与更新的最优权值Wopt进行乘加运算后滤波输出,实现数据的实时滤波。
为了验证本发明方法能够在高动态环境下具有更精确的协方差估计,可以达到零陷展宽的效果,特举一个验证例进行说明。
验证例1
本验证例的方法,涉及高动态应用环境下阵列信号的稳健自适应滤波领域。针对高动态应用环境,提供了一种高动态下的零陷展宽方法,可以与高性能FPGA结合的稳健滤波,具有较强的可操作性。具体发明应用于由4阵元组成的圆形阵列,详细步骤如下:
(1)对采样的四路AD信号进行带通滤波,Hilbert变换以及幅相校正后进入步骤(2);
(2)对步骤(1)中的接收数据分两路输出,一路进入FIFO模块存储,等待权值更新后进行滤波处理,一路进行等间隔抽取后送入QR分解单元,进行协方差矩阵估计,然进入步骤(3);
QR分解的脉动阵结构GR-TSA如图2所示,四路数据为了满足脉动阵的时序结构经过延迟处理分别依次延迟一个脉动间隔t,在第ti个脉动时刻,第一排数据进入脉动结构,x1(t1)经过边界单元计算后得出旋转因子c1/s1,边界单元的过程如图3所示,主要用来计算对角元素rii和产生内部单元所需要的Givens旋转因子c和s,若xin=0,则c=1,s=0,rii=rii;若xin≠0,则为了避开除法运算,当xin≠0时,边界单元的运算可以改写为:其中使用两次CORDIC旋转,第一次使用CORDIC旋转来实现输入数据的求模运算得到输入数据的模值|xin|和旋转角θ,第二次使用CORDIC旋转来实现对角元素更新和旋转角根据更新后旋转角θ和调用FPGA中的直接数字频率合成器DDS IP核,计算出sinθ、cosθ,调用乘法器,根据公式更新旋转因子c和s,进入步骤(3);
(3)内部单元主要作用是根据边界单元传递过来的Givens变换因子进行旋转消元,消元求出非对角元素rij。内部单元中计算非对角元素r′ij=s*xin+crij,经过一次Givens变换后的待消元元素xout=cxin-srij,根据边界单元计算的c和s,结合非对角元素rij,xin为阵列天线信号输入数据(第二次Givens变换时输入的xin是第一次Givens变换后计算出的输出数据xout),由于变换因子s、非对角元素rij、信号输入数据xin均为复数,内部单元FPGA实现如图4所示,执行两次复数乘法和四次实数乘法求出一次Givens变换xout和r′ij。
在第ti+1个脉动时刻,第二排数据进入脉动结构,此时x2(t1)经过内部单元的一次计算,得出经过一次Givens的旋转因子c1/s1的上三角矩阵R的非对角元素r12和经过变换后待消元元素xout12,在第ti+2个脉动时刻,x3(t1)经过内部单元的一次计算,得出经过一次Givens的旋转因子c1/s1的上三角矩阵R的非对角元素r13和经过变换后待消元元素xout13,依次类推,每经过一次边界单元消元一个元素,采用流水线结构经过7个脉动周期便可以将新进的元素消元,完成更新后的上三角矩阵,在t=ti+7时刻输出R。进入步骤(4);
(6)FIFO模块将存储的数据输出,并与权值W乘加后滤波输出。
综上,本发明方法在QRD-SMI算法基础上,结合高动态环境的特点,提出抽取式QRD-SMI滤波方法,等效扩展了干扰零陷,增强了滤波的稳定性,并进行了FPGA实现。
需要说明的是,以上说明仅是本发明的优选实施方式,应当理解,对于本领域技术人员来说,在不脱离本发明技术构思的前提下还可以做出若干改变和改进,这些都包括在本发明的保护范围内。
Claims (4)
1.一种高动态下的零陷展宽方法,其特征在于,包括以下步骤:
步骤1、获取接收信号的估计协方差矩阵的采样数据,对估计协方差矩阵的采样数据进行带通滤波、Hilbert变换以及幅相校正后进入步骤2;
步骤2、对所获取估计协方差矩阵的采样数据等间隔抽取后分两路输出,其中一路输入FIFO模块存储,及另一路进行等间隔抽取后送入QR分解单元进行协方差矩阵估计,输出得到更新的协方差矩阵R;其中QR分解由N次Givens旋转进行输入数据消元实现,抽取间隔Td大于一次Givens旋转时间,协方差估计样本数N取值大于2M,M为接收阵元数;
步骤3、根据更新的协方差矩阵R代入式RHV=La(θd)中进行中间变量V的前向迭代更新,得到更新的中间变量V,其中a(θd)为角度θd对应的导向矢量;L为采样次数;
步骤4、根据更新的协方差矩阵R和中间变量V进行后向迭代计算,得到更新的最优权值Wopt;
步骤5、将FIFO模块所存储估计协方差矩阵的采样数据输出 ,并与更新的最优权值Wopt进行乘加运算后滤波输出。
2.根据权利要求1所述高动态下的零陷展宽方法,其特征在于,所述步骤2中对所获取估计协方差矩阵的采样数据等间隔抽取,具体为:
每间隔T时刻抽取一次采样快拍,并连续操作N次;
获取第i组NT时段快拍的入射角θi,得到接收信号模型及其对应的协方差矩阵,根据协方差矩阵求逆引理,得到加权矢量的最优近似解。
3.根据权利要求1所述高动态下的零陷展宽方法,其特征在于,所述步骤2中QR分解单元进行QR分解具体为:
将所获取估计协方差矩阵的采样数据经过延迟处理分别依次延迟一个脉动间隔t,在第ti个脉动时刻,第一排数据x1(t1)进入脉动结构,经过边界单元计算后得出旋转因子c1/s1,该边界单元计算对角元素rii和产生内部单元所需要的Givens的旋转因子c和s,其中使用两次CORDIC旋转,第一次使用CORDIC旋转来实现输入数据的求模运算得到输入数据的模值|xin|和旋转角θ,第二次使用CORDIC旋转来实现对角元素rii和旋转角更新;计算旋转角θ和的正余弦值,更新得到旋转因子c和s;
在第ti+1个脉动时刻,第二排数据x2(t1)进入脉动结构,经过内部单元的一次计算,得出经过一次Givens的旋转因子c1/s1的上三角矩阵R的非对角元素r12和经过变换后待消元元素xout12;
在第ti+2个脉动时刻,第三排数据x3(t1)经过内部单元的一次计算,得出经过一次Givens的旋转因子c1/s1的上三角矩阵R的非对角元素r13和经过变换后待消元元素xout13,依次类推,每经过一次边界单元消元一个元素,采用流水线结构经过7个脉动周期便可以将新进的元素消元,完成更新后的上三角矩阵,及在t=ti+7时刻输出更新的协方差矩阵R。
4.根据权利要求1所述高动态下的零陷展宽方法,其特征在于,所述步骤4得到更新的最优权值Wopt采用公式:RWopt=V。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811607382.9A CN109787676B (zh) | 2018-12-27 | 2018-12-27 | 一种高动态下的零陷展宽方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811607382.9A CN109787676B (zh) | 2018-12-27 | 2018-12-27 | 一种高动态下的零陷展宽方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109787676A CN109787676A (zh) | 2019-05-21 |
CN109787676B true CN109787676B (zh) | 2021-07-06 |
Family
ID=66497757
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811607382.9A Active CN109787676B (zh) | 2018-12-27 | 2018-12-27 | 一种高动态下的零陷展宽方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109787676B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115033843B (zh) * | 2022-08-09 | 2022-11-08 | 之江实验室 | 基于三角脉动阵列的协方差矩阵计算的电路实现方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5956624A (en) * | 1994-07-12 | 1999-09-21 | Usa Digital Radio Partners Lp | Method and system for simultaneously broadcasting and receiving digital and analog signals |
CN101340648A (zh) * | 2008-08-11 | 2009-01-07 | 中兴通讯股份有限公司 | 一种多用户波束赋形发射权值的生成方法 |
CN101482605A (zh) * | 2009-01-22 | 2009-07-15 | 中国民航大学 | 加宽零陷的高动态卫星导航系统干扰抑制方法 |
CN102508213A (zh) * | 2011-11-04 | 2012-06-20 | 西安电子科技大学 | 基于零陷加宽的波束域自适应干扰抑制方法 |
CN104702379A (zh) * | 2013-12-10 | 2015-06-10 | 北京信威通信技术股份有限公司 | 一种无线通信系统中的动态零陷方法 |
CN108663693A (zh) * | 2018-07-25 | 2018-10-16 | 电子科技大学 | 一种基于空时处理的高动态gnss零陷展宽干扰抑制方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7904243B2 (en) * | 2004-01-15 | 2011-03-08 | The Boeing Company | Real-time data aiding for enhanced GPS performance |
-
2018
- 2018-12-27 CN CN201811607382.9A patent/CN109787676B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5956624A (en) * | 1994-07-12 | 1999-09-21 | Usa Digital Radio Partners Lp | Method and system for simultaneously broadcasting and receiving digital and analog signals |
CN101340648A (zh) * | 2008-08-11 | 2009-01-07 | 中兴通讯股份有限公司 | 一种多用户波束赋形发射权值的生成方法 |
CN101482605A (zh) * | 2009-01-22 | 2009-07-15 | 中国民航大学 | 加宽零陷的高动态卫星导航系统干扰抑制方法 |
CN102508213A (zh) * | 2011-11-04 | 2012-06-20 | 西安电子科技大学 | 基于零陷加宽的波束域自适应干扰抑制方法 |
CN104702379A (zh) * | 2013-12-10 | 2015-06-10 | 北京信威通信技术股份有限公司 | 一种无线通信系统中的动态零陷方法 |
CN108663693A (zh) * | 2018-07-25 | 2018-10-16 | 电子科技大学 | 一种基于空时处理的高动态gnss零陷展宽干扰抑制方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109787676A (zh) | 2019-05-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Liu et al. | Sparsity-inducing direction finding for narrowband and wideband signals based on array covariance vectors | |
CN110113085B (zh) | 一种基于协方差矩阵重构的波束形成方法及系统 | |
CN103901395B (zh) | 一种冲击噪声环境下相干信号波达方向动态跟踪方法 | |
CN103954931A (zh) | 一种远场和近场混合信号源的定位方法 | |
CN105137454B (zh) | 一种基于协方差矩阵特征分解的抗干扰算法的fpga实现方法及实现装置 | |
CN112180320B (zh) | 一种无人机无源定位系统及方法 | |
CN110515038A (zh) | 一种基于无人机-阵列的自适应无源定位装置及实现方法 | |
CN109143154A (zh) | 一种应用于l型阵列的信号二维doa与频率联合估计方法 | |
CN109061686B (zh) | 基于递归广义最大互熵的自适应多径估计方法 | |
CN110850445A (zh) | 一种基于空时采样协方差求逆的脉冲干扰抑制方法 | |
CN112731273A (zh) | 一种基于稀疏贝叶斯的低复杂度信号波达方向估计方法 | |
CN106603036A (zh) | 一种基于低阶内插滤波器的自适应时延估计方法 | |
CN109787676B (zh) | 一种高动态下的零陷展宽方法 | |
CN111880198A (zh) | 基于交替极化敏感阵列的空时极化抗干扰方法 | |
CN113504549B (zh) | 基于广义旁瓣相消器的导航空时抗干扰方法 | |
CN105353338B (zh) | 宽带信号超分辨测向中的阵列通道幅相不一致性误差校正方法 | |
CN104950282B (zh) | 连续域内稀疏重构实现的宽带信号超分辨测向方法及装置 | |
CN114884841A (zh) | 基于高阶统计和非均匀阵列的欠定参数联合估计方法 | |
CN105242236B (zh) | 宽带信号超分辨测向中的阵元位置误差校正方法 | |
CN110749855A (zh) | 一种基于协方差域零化的均匀线阵波达方向估计方法 | |
CN105227227B (zh) | 一种基于小波的智能天线波束形成系统和方法 | |
CN105182279B (zh) | 基于空域稀疏优化的宽带信号超分辨测向误差校正方法 | |
CN110058192A (zh) | 基于辅助阵元的原子范数互耦doa估计方法 | |
CN105223541B (zh) | 宽带信号超分辨测向中的阵元间互耦误差校正方法 | |
Chen et al. | Error analysis of convolutional beamspace algorithms |
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 |