CN113824488B - 基于判决反馈自适应对消的卫星通信非恶意干扰抑制方法 - Google Patents

基于判决反馈自适应对消的卫星通信非恶意干扰抑制方法 Download PDF

Info

Publication number
CN113824488B
CN113824488B CN202111055822.6A CN202111055822A CN113824488B CN 113824488 B CN113824488 B CN 113824488B CN 202111055822 A CN202111055822 A CN 202111055822A CN 113824488 B CN113824488 B CN 113824488B
Authority
CN
China
Prior art keywords
interference
signal
satellite communication
frequency
bandwidth
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
CN202111055822.6A
Other languages
English (en)
Other versions
CN113824488A (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.)
CETC 54 Research Institute
Original Assignee
CETC 54 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 CETC 54 Research Institute filed Critical CETC 54 Research Institute
Priority to CN202111055822.6A priority Critical patent/CN113824488B/zh
Publication of CN113824488A publication Critical patent/CN113824488A/zh
Application granted granted Critical
Publication of CN113824488B publication Critical patent/CN113824488B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B7/00Radio transmission systems, i.e. using radiation field
    • H04B7/14Relay systems
    • H04B7/15Active relay systems
    • H04B7/185Space-based or airborne stations; Stations for satellite systems
    • H04B7/1851Systems using a satellite or space-based relay
    • H04B7/18513Transmission in a satellite or space-based system
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02DCLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
    • Y02D30/00Reducing energy consumption in communication networks
    • Y02D30/70Reducing energy consumption in communication networks in wireless communication networks

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Astronomy & Astrophysics (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Noise Elimination (AREA)
  • Radio Relay Systems (AREA)

Abstract

本发明公开了一种基于判决反馈自适应对消的卫星通信非恶意干扰抑制方法,属于卫星通信技术领域。该方法利用非恶意干扰易获取频点、带宽等先验知识的优势,能够有效解决卫星通信中频域重叠非恶意干扰的抑制问题,适用于不同调制方式的扩频和非扩频卫星通信传输链路。该干扰抑制方法对卫通通信信号的调制方式和扩频处理无特殊要求,能够适用于不同的卫星通信传输链路,具有强的适应能力。

Description

基于判决反馈自适应对消的卫星通信非恶意干扰抑制方法
技术领域
本发明属于卫星通信技术领域,特别涉及一种基于判决反馈自适应对消的卫星通信非恶意干扰抑制方法,可用于不同调制方式的扩频和非扩频卫星通信传输系统。
背景技术
目前,卫星通信因其通信距离远和不受地理地貌限制等优势,得到了广泛应用。同时,卫星通信覆盖了从UHF、S、X、Ku和Ka等多个频段,频率范围可达几十个GHz。由于卫星通信工作频段广和传输信道的开放性,导致其传输链路不可避免地面临来自其他系统或用户的非恶意强干扰,将恶化链路传输性能,其中以UHF/S频段最为严重,比如地面雷达、移动4G、民航通信、地面集群等系统发射的同频或邻道信号等。因此,对于卫星通信传输系统需要进一步研究如何有效抑制频域重叠的非恶意强干扰。
卫星通信系统传输链路中的非恶意干扰相比于对抗环境中的恶意干扰,通常具有固定不变和部分参数已知的特点。因此,可充分利用非恶意干扰的先验知识,实现对卫星通信系统传输链路中非恶意干扰的有效抑制。在研究UHF频段卫星通信系统上行和下行链路干扰特性的基础上,利用干扰检测技术和动态信道分配相结合,能够一定程度上避免非恶意干扰,进而提高系统可用性,具体可参见文献:刘爱军,王杭先“UHF频段卫星通信系统干扰分析及检测技术研究”,军事通信技术,2007年,12卷,第4期,页码:21-25。对于直接序列扩频传输链路中的窄带干扰,可利用条件约束非线性内插滤波方法实现抑制,但对信号带宽有较为严格的限制,具体可参见文献:Leslie A R,Vincent P H“NarrowbandInterference Suppression in CDMA Spread Spectrum Communications”,IEEETransactions on Communications,1994年,42卷,第4期,页码:1969-1979。对于扩频传输系统中的窄带干扰抑制,还可采用小波变换和IIR格型滤波等抑制方法,不过小波变换实现复杂,而格型滤波只能适用于近似单音的极窄带干扰,具体可参见文献:王桁,吕智勇,杨龙“DSSS卫星通信中基于小波包变换的干扰抑制方法”,系统工程与电子技术,2016年,38卷,第6期,页码:1417-1422;龚文飞,吴嗣亮,李加琪“直扩系统中IIR格型滤波器抑制窄带干扰新方法与性能分析”,电子与信息学报,2010年,32卷,第10期,页码:2473-2478。在能够恢复干扰信号情况下,可采用干扰对消与信道失真补偿的方法对接收的干扰信号进行抑制,通过干扰信道状态估计和相应评判准则下的信号失真补偿实现干扰的有效抑制,具体可参见文献:Zhang S,Liew S,Wang H“Blind Known Interference Cancellation”,IEEEJournal on SelectedAreas in Communications,2013年,31卷,第8期,页码:1572-1582。
但是,对于卫星通信系统应用,不可避免面临资源冲突或者其他系统信号的干扰。因此,如何有效抑制传输链路的非恶意干扰问题,成为提升卫星通信系统传输性能的关键。
发明内容
本发明针对卫星通信系统,为实现对非恶意干扰的有效抑制,进而显著提升系统的链路传输整体性能,提供一种基于判决反馈自适应对消的卫星通信非恶意干扰抑制方法。该方法对干扰信号带宽无特殊要求,能够适用于不同调制方式的扩频和非扩频卫星通信传输链路。
为实现上述目的,本发明采用的技术方案为:
一种基于判决反馈自适应对消的卫星通信非恶意干扰抑制方法,其应用于卫星通信系统,用于对卫星通信信号进行非恶意干扰抑制处理,包括如下步骤:
(1)对模数转换器输出的数字卫星通信信号进行处理,以期望信号频率为中心,设定需要进行干扰检测的接收带宽BΔ,在该带宽内进行干扰信号检测;
(2)当接收信号中存在干扰时,依据包括估计或已知的干扰频点、带宽和调制方式在内的先验知识,从接收信号中恢复干扰信号,同时以多项式模型构造干扰信号的非线性分量;
(3)以重构的干扰信号作为自适应对消的参考信号,通过快速收敛自适应算法实现接收信号中的干扰抑制,自适应对消误差信号即为干扰抑制输出的待解调期望信号。
进一步的,步骤(1)包括以下步骤:
(1-1)依据卫星通信信号传输性能要求,确定模数转换器采样接收信号带宽BΔ,确保距离接收中心频点fi间隔不小于BΔ/2的信号将不对期望信号产生影响;
(1-2)对接收信号的数据序列r(n)进行加窗离散傅里叶变换处理:
Figure BDA0003254566040000031
式中,N表示傅里叶变换点数,X(k)表示变换输出,DFT(x)表示对x进行离散傅里叶变换计算,WN表示窗函数,w(n)表示窗函数对应的系数;
(1-3)对序列X(k)修正,得到对称的频谱序列Xi(k);计算Xi(k)的模值,与设定检测门限γ比较,若存在大于门限的模值,则表示存在干扰,此时,依据大于门限的模值所对应的位置,估计干扰信号中心频率。
进一步的,步骤(2)包括以下步骤:
(2-1)根据干扰信号检测结果和预知的先验知识,确定干扰信号的中心频率fJ、带宽BJ和调制方式;
(2-2)对接收信号进行数字混频,将干扰信号变至零频,在完成抽取和匹配滤波处理之后,进行定时和解调映射,得到干扰信号的0、1序列;
(2-3)根据解调映射输出的干扰0、1序列,进行原干扰信号调制方式的调制映射、内插和成形滤波处理,得到重构的基带干扰信号;
(2-4)利用多项式模型和重构的基带干扰信号构建其非线性分量,对重构干扰信号的线性和非线性分量进行混频,将频谱中心频率搬移至fJ,完成干扰信号重构。
进一步的,步骤(3)包括以下步骤:
(3-1)利用迭代变步长最小均方算法构造自适应对消器,确定最小均方算法的阶数;
(3-2)将干扰通道重构的干扰信号线性和非线性分量作为自适应对消参考信号,调整期望信号接收通道的延时与参考信号处理时延相匹配,将延时后的接收信号与加权的参考信号相减;
(3-3)依据接收信号与参考信号相减后的误差,反馈调整自适应对消器的加权系数,保证输出的误差信号均方值最小,完成自适应对消,将收敛状态下的稳态误差作为干扰抑制的输出信号。
本发明的有益效果在于:
1、本发明能够对卫星通信传输链路的干扰进行有无检测,同时实现非恶意干扰的有效抑制。
2、本发明对卫通通信信号的调制方式和扩频处理无特殊要求,能够适用于不同调制方式的扩频和非扩频卫星通信传输链路,具有较强的适应能力。
3、本发明利用非恶意干扰易获取频点、带宽等先验知识的优势,能够有效解决卫星通信中频域重叠非恶意干扰的抑制问题。
4、相比于频域陷波等干扰抑制方法,本发明能够显著降低对期望信号的损伤。
总之,本发明在进行干扰有无检测基础上,充分利用非恶意干扰易获取的先验知识,对接收信号中的干扰进行解调判决,再进行调制映射重构,以重构的干扰信号作为参考,利用快速收敛算法实现与接收信号的自适应对消,进而有效抑制接收信号中的非恶意干扰。
附图说明
图1为本发明实施例方法的原理框图;
图2为本发明实施例中非恶意干扰检测的实现原理图;
图3为本发明实施例中在不同信噪比条件下的干扰抑制性能图。
具体实施方式
下面结合附图和具体实施方式,对本发明的技术方案做进一步说明。
一种基于判决反馈自适应对消的非恶意干扰源抑制方法,包括:
(1)该干扰抑制方法以ADC采样后的数字处理为基础,对接收的宽带信号实现干扰检测;
(2)通过“解调+调制”架构实现干扰信号重构,重构采用多项式模型构造干扰的非线性分量;
(3)以重构的干扰信号为参考,与延迟后的接收混合信号相减,利用迭代变步长最小均方算法进行自适应对消;
(4)在自适应对消收敛状态下,对消输出的误差信号即可作为待解调的干扰抑制后信号。
其采用如下具体步骤实现:
首先,进行卫星通信系统的非恶意干扰检测以期望信号接收频率为中心,在BΔ频带范围内检测是否有强的其他信号存在,干扰通道提供干扰信号的频率、带宽和调制方式等特征参数。
对带宽BΔ的采样信号r(n),截取长度为N的序列进行N点加窗离散傅里叶变换,得到相应的频谱函数X(k)。因X(k)用以检测干扰并估计频率、带宽和调制方式等参数,所以其频率分辨率需要满足干扰参数估计的需求。
利用窗函数WN对基带信号r(n)进行加窗和N点DFT变换:
Figure BDA0003254566040000061
由于DFT变换得到序列包含正负频率分量,对序列X(k)修正,得到对称的频谱序列Xi(k):
Figure BDA0003254566040000062
计算Xi(k)模值,并与门限γ比较,得到最终检测结果。
存在干扰情况下,对接收信号中的干扰进行解调和再调制,同时进行频谱搬移,得到重构的干扰信号
Figure BDA0003254566040000063
再根据期望信号接收通道的L阶匹配滤波响应h(n)和多项式非线性模型计算线性和非线性分量,过程如下:
Figure BDA0003254566040000064
Figure BDA0003254566040000071
Figure BDA0003254566040000072
其中,k表示设置的非线性参数,通常取值为1,k1=2(k-1)+1,k2=2k+1,k3=2(k+1)+1表示非线性分量阶数,取值为1时表示线性分量。
设定最小均方算法阶数M,将重构干扰信号的各阶分量作为最小均方算法的参考信号输入向量
Figure BDA0003254566040000073
过程描述为:
Figure BDA0003254566040000074
其中,M=M1+M2+M3
采用最小均方算法,以信号接收通道处理和延时调整的信号rd(n)作为期望信号,重构的干扰信号
Figure BDA0003254566040000075
作为参考信号,自适应对消过程描述如下:
Figure BDA0003254566040000076
Figure BDA0003254566040000077
Figure BDA0003254566040000078
Figure BDA0003254566040000079
其中,e(n)为对消输出的误差信号,wM为最小均方算法的加权系数向量,
Figure BDA00032545660400000710
Figure BDA00032545660400000711
分别为M1×M1、M2×M2和M3×M3维的单位矩阵,μ1(n)、μ2(n)和μ3(n)分别为1阶、2阶和3阶参考信号所对应的不同时刻步长因子。在最小均方算法处于收敛状态后,误差信号e(n)即为判决反馈自适应对消干扰抑制的输出结果。
该方法可对卫星通信系统中的非恶意干扰进行有效抑制,能够消除频域重叠非恶意干扰对链路传输性能的影响,显著降低其传输误码率,进而提升系统的整体传输性能。
以下为一个更具体的例子:
设定某卫星通信系统采用8倍扩频的BPSK调制方式,扩频后的码片速率为9.6Mcps,其他系统产生的非恶意干扰同样为BPSK调制信号,其符号速率为960ksps,即干扰与通信信号的带宽占比为10%,采用的最小均方算法阶数为60,多项式非线性模型设为3阶,干扰信号的邻道功率比约为30dB。
在上述实例参数条件下,采用高斯白噪声信道对本发明的干扰抑制方法进行性能仿真,统计平均1000次蒙特卡罗仿真结果。图3中为不同信干比情况下,本方法在收敛状态时对非恶意干扰的抑制比随接收信号信噪比的变化曲线。可见,本方法在传输信道信噪比不低于10dB时,对干信比大于等于0dB的干扰,均能够取得20dB以上的干扰抑制比,实现对非恶意干扰的有效抑制。
总之,本发明方法在进行干扰有无检测基础上,充分利用非恶意干扰易获取的先验知识,实现干扰信号的时域重构和自适应对消,完成接收信号中干扰分量的有效抑制,相比于常用频域陷波方法能够减小对期望信号频谱的损伤,具有更高的处理增益。本发明方法可适用于不同调制方式的扩频和非扩频传输链路,能够有效抑制卫星通信系统面临的非恶意干扰,为改善系统传输性能提供有力支撑。

Claims (1)

1.一种基于判决反馈自适应对消的卫星通信非恶意干扰抑制方法,其特征在于,应用于卫星通信系统,用于对卫星通信信号进行非恶意干扰抑制处理,包括如下步骤:
(1)对模数转换器输出的数字卫星通信信号进行处理,以期望信号频率为中心,设定需 要进行干扰检测的接收带宽
Figure DEST_PATH_IMAGE001
,在该带宽内进行干扰信号检测;具体方式如下:
(1-1)依据卫星通信信号传输性能要求,确定模数转换器采样接收信号带宽
Figure 937374DEST_PATH_IMAGE001
,确保距 离接收中心频点
Figure 961174DEST_PATH_IMAGE002
间隔不小于
Figure DEST_PATH_IMAGE003
的信号将不对期望信号产生影响;
(1-2)对接收信号的数据序列
Figure 209752DEST_PATH_IMAGE004
进行加窗离散傅里叶变换处理:
Figure DEST_PATH_IMAGE005
式中,N表示傅里叶变换点数,
Figure 932858DEST_PATH_IMAGE006
表示变换输出,
Figure DEST_PATH_IMAGE007
表示对x进行离散傅里叶变换 计算,
Figure 216072DEST_PATH_IMAGE008
表示窗函数,
Figure DEST_PATH_IMAGE009
表示窗函数对应的系数;
(1-3)对序列X(k)修正,得到对称的频谱序列
Figure 366430DEST_PATH_IMAGE010
;计算
Figure 367884DEST_PATH_IMAGE010
的模值,与设定检测门限
Figure DEST_PATH_IMAGE011
比较,若存在大于门限的模值,则表示存在干扰,此时,依据大于门限的模值所对应的位 置,估计干扰信号中心频率;
(2)当接收信号中存在干扰时,依据包括估计或已知的干扰频点、带宽和调制方式在内的先验知识,从接收信号中恢复干扰信号,同时以多项式模型构造干扰信号的非线性分量;具体方式如下:
(2-1)根据干扰信号检测结果和预知的先验知识,确定干扰信号的中心频率
Figure 160260DEST_PATH_IMAGE012
、带宽
Figure DEST_PATH_IMAGE013
和调制方式;
(2-2)对接收信号进行数字混频,将干扰信号变至零频,在完成抽取和匹配滤波处理之后,进行定时和解调映射,得到干扰信号的0、1序列;
(2-3)根据解调映射输出的干扰0、1序列,进行原干扰信号调制方式的调制映射、内插和成形滤波处理,得到重构的基带干扰信号;
(2-4)利用多项式模型和重构的基带干扰信号构建其非线性分量,对重构干扰信号的 线性和非线性分量进行混频,将频谱中心频率搬移至
Figure 658499DEST_PATH_IMAGE012
,完成干扰信号重构;
(3)以重构的干扰信号作为自适应对消的参考信号,通过快速收敛自适应算法实现接收信号中的干扰抑制,自适应对消误差信号即为干扰抑制输出的待解调期望信号;具体方式如下:
(3-1)利用迭代变步长最小均方算法构造自适应对消器,确定最小均方算法的阶数;
(3-2)将干扰通道重构的干扰信号线性和非线性分量作为自适应对消参考信号,调整期望信号接收通道的延时与参考信号处理时延相匹配,将延时后的接收信号与加权的参考信号相减;
(3-3)依据接收信号与参考信号相减后的误差,反馈调整自适应对消器的加权系数,保证输出的误差信号均方值最小,完成自适应对消,将收敛状态下的稳态误差作为干扰抑制的输出信号。
CN202111055822.6A 2021-09-09 2021-09-09 基于判决反馈自适应对消的卫星通信非恶意干扰抑制方法 Active CN113824488B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111055822.6A CN113824488B (zh) 2021-09-09 2021-09-09 基于判决反馈自适应对消的卫星通信非恶意干扰抑制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111055822.6A CN113824488B (zh) 2021-09-09 2021-09-09 基于判决反馈自适应对消的卫星通信非恶意干扰抑制方法

Publications (2)

Publication Number Publication Date
CN113824488A CN113824488A (zh) 2021-12-21
CN113824488B true CN113824488B (zh) 2022-07-08

Family

ID=78914297

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111055822.6A Active CN113824488B (zh) 2021-09-09 2021-09-09 基于判决反馈自适应对消的卫星通信非恶意干扰抑制方法

Country Status (1)

Country Link
CN (1) CN113824488B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114866128B (zh) * 2022-04-07 2023-09-05 中国人民解放军战略支援部队信息工程大学 基于扩频信号的卫星通信临界干扰功率阈值估计方法及系统

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103067066A (zh) * 2012-12-18 2013-04-24 哈尔滨工业大学 双天线卫星通信系统下行链路干扰抑制方法
CN105049067A (zh) * 2015-08-31 2015-11-11 北京理工大学 一种基于导频信号的无源互调干扰对消方法
CN107872268A (zh) * 2016-09-26 2018-04-03 北京大学(天津滨海)新代信息技术研究院 一种用于卫星通信系统消除干扰的方法
WO2020134855A1 (zh) * 2018-12-28 2020-07-02 长沙天仪空间科技研究院有限公司 一种卫星通信系统
CN111628816A (zh) * 2019-01-03 2020-09-04 长沙天仪空间科技研究院有限公司 一种卫星通信系统的窄带干扰抑制方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103067066A (zh) * 2012-12-18 2013-04-24 哈尔滨工业大学 双天线卫星通信系统下行链路干扰抑制方法
CN105049067A (zh) * 2015-08-31 2015-11-11 北京理工大学 一种基于导频信号的无源互调干扰对消方法
CN107872268A (zh) * 2016-09-26 2018-04-03 北京大学(天津滨海)新代信息技术研究院 一种用于卫星通信系统消除干扰的方法
WO2020134855A1 (zh) * 2018-12-28 2020-07-02 长沙天仪空间科技研究院有限公司 一种卫星通信系统
CN111628816A (zh) * 2019-01-03 2020-09-04 长沙天仪空间科技研究院有限公司 一种卫星通信系统的窄带干扰抑制方法

Also Published As

Publication number Publication date
CN113824488A (zh) 2021-12-21

Similar Documents

Publication Publication Date Title
CN109639303B (zh) 一种基于加窗处理的干扰检测和抑制方法
US8699448B2 (en) Method and arrangement for improved interference rejection combining
CN102904604B (zh) 一种窄带干扰抑制方法和装置
CN106772457B (zh) 一种基于空时频架构的卫星导航抗干扰方法
KR100843253B1 (ko) 반복 수신 방법 및 반복 수신기
EP2928140B1 (en) Method and a device for cancelling a narrow band interference in a single carrier signal
CN113824488B (zh) 基于判决反馈自适应对消的卫星通信非恶意干扰抑制方法
CN110933007B (zh) 一种削除针对ofdm宽频信号干扰的装置及方法
WO2012055344A1 (zh) 估计窄带干扰的方法、装置及接收设备
KR101156131B1 (ko) 위성 중계기에서 간섭제거 방법 및 장치
CN101154983B (zh) 一种基于单倍采样的卫星通信抗干扰技术的实现方法
KR100748642B1 (ko) 이동 통신 중계기의 간섭 신호 제거 방법
CN113541706B (zh) 一种基于变换域处理的窄带干扰抑制方法
US9048917B1 (en) Method and apparatus for estimating noise covariance in a wireless network
CN104954304A (zh) 一种迭代干扰抑制方法及装置
CN113406672A (zh) 一种基于自适应加窗的分段线性调频干扰消除方法
JP6552753B2 (ja) チャネル推定を実行する方法及びデバイス
CN110798178A (zh) 一种时域自适应滤波方法
Chen et al. Research on anti-jamming algorithm of BDS signal in WiFi6 environment
CN100438337C (zh) 一种基于Laguerre结构的自适应非线性滤波方法
CN113630355B (zh) 基于空时功率倒置阵的宽带干扰抑制装置及抑制方法
Shirai et al. Frequency Utilization Efficiency Improvement by Using MLD in Spectrum Suppressed Transmission
Sun et al. Channel identification and interference compensation for OFDM system in long multipath environment
CN116318246A (zh) 一种基于窄带干扰抑制的直接序列扩频抗干扰方法
CN118233260A (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