CN113765607B - 一种频谱检测方法及装置 - Google Patents

一种频谱检测方法及装置 Download PDF

Info

Publication number
CN113765607B
CN113765607B CN202110745871.6A CN202110745871A CN113765607B CN 113765607 B CN113765607 B CN 113765607B CN 202110745871 A CN202110745871 A CN 202110745871A CN 113765607 B CN113765607 B CN 113765607B
Authority
CN
China
Prior art keywords
result
spectrum
detection
fuzzy
preset
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
CN202110745871.6A
Other languages
English (en)
Other versions
CN113765607A (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.)
State Grid Corp of China SGCC
State Grid Jiangsu Electric Power Co Ltd
Global Energy Interconnection Research Institute
State Grid Fujian Electric Power Co Ltd
Information and Telecommunication Branch of State Grid Fujian Electric Power Co Ltd
Original Assignee
State Grid Corp of China SGCC
State Grid Jiangsu Electric Power Co Ltd
Global Energy Interconnection Research Institute
State Grid Fujian Electric Power Co Ltd
Information and Telecommunication Branch of State Grid Fujian Electric Power Co Ltd
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 State Grid Corp of China SGCC, State Grid Jiangsu Electric Power Co Ltd, Global Energy Interconnection Research Institute, State Grid Fujian Electric Power Co Ltd, Information and Telecommunication Branch of State Grid Fujian Electric Power Co Ltd filed Critical State Grid Corp of China SGCC
Publication of CN113765607A publication Critical patent/CN113765607A/zh
Application granted granted Critical
Publication of CN113765607B publication Critical patent/CN113765607B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B17/00Monitoring; Testing
    • H04B17/30Monitoring; Testing of propagation channels
    • H04B17/391Modelling the propagation channel
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B17/00Monitoring; Testing
    • H04B17/30Monitoring; Testing of propagation channels
    • H04B17/309Measuring or estimating channel quality parameters
    • H04B17/336Signal-to-interference ratio [SIR] or carrier-to-interference ratio [CIR]
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B17/00Monitoring; Testing
    • H04B17/30Monitoring; Testing of propagation channels
    • H04B17/391Modelling the propagation channel
    • H04B17/3913Predictive models, e.g. based on neural network models
    • 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)
  • Electromagnetism (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Quality & Reliability (AREA)
  • Artificial Intelligence (AREA)
  • Evolutionary Computation (AREA)
  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明提供了一种频谱检测方法及装置,其中,该方法包括:获取多个初始频谱感知结果;若多个初始频谱感知结果中不包含一位判决结果,或,一位判决结果满足预设条件,根据初始频谱感知结果中的两位判决结果和预设模糊判决模型得到模糊逻辑输出结果;若模糊逻辑输出结果属于第一预设模糊级别,获取历史感知结果序列;对历史感知结果序列进行指数平滑预测,得到指数平滑预测结果;根据指数平滑预测结果得到最终频谱检测结果。在本发明提供的频谱检测方法及装置中,当无法直接根据模糊逻辑输出结果确定最终频谱检测结果时,结合多个时期的感知结果对当前时期的感知结果进行指数平滑预测,从而使得到的当前时期的最终频谱检测结果更准确。

Description

一种频谱检测方法及装置
技术领域
本发明涉及频谱感知技术领域,具体涉及一种频谱检测方法及装置。
背景技术
随着泛在电力物联网建设的推进,越来越多的用户和各类服务都将应用到电网中,这将导致电网对频谱资源的需求不断增加。然而,由于电力无线专用网的频谱资源有限,鼓励采用动态频谱共享技术来利用频谱资源。因此,认知无线电技术在电网中的应用具有重要意义。频谱感知技术是认知无线电的前提,因此提高频谱感知技术的检测性能成为研究热点。目前常用的频谱感知技术为能量检测算法,但是,在低信噪比环境下,如果噪声存在不确定性,能量检测算法的检测性能会下降。
为了提高检测性能,现有技术中采用传统的能量检测算法得到能量检测值、第一门限值、第二门限值,当能量检测值位于第一门限值与第二门限值之间时,将第一门限值与第二门限值之间的区域划分为多个模糊区域,并根据能量检测值所属的模糊区域和模糊规则集得到模糊输出,模糊输出包括三类:Low、Medium、High,当模糊输出为Low时,确定能量检测结果为不存在频谱信号,当模糊输出为High时,确定能量检测结果为存在频谱信号,当模糊输出为Medium时,表示无判决结果,以上一时期的检测结果作为模糊输出为Medium时的检测结果。
通过以上方法确定检测结果时,若模糊输出为Medium,则将上一时期的检测结果作为当前时期的检测结果,但是,无法保证上一时期的检测结果是准确的,因此,当上一时期的检测结果存在错误时,当前时期的检测结果一定是错误的。
由此可见,通过现有技术进行频谱感知时,得到的感知结果可靠性较差。
发明内容
因此,本发明要解决的技术问题在于克服现有技术中得到的感知结果可靠性较差的缺陷,从而提供一种频谱检测方法及装置。
本发明第一方面提供了一种频谱检测方法,包括:获取多个初始频谱感知结果;若多个初始频谱感知结果中不包含一位判决结果,或,一位判决结果满足预设条件,根据初始频谱感知结果中的两位判决结果和预设模糊判决模型得到模糊逻辑输出结果;若模糊逻辑输出结果属于第一预设模糊级别,获取历史感知结果序列;对历史感知结果序列进行指数平滑预测,得到指数平滑预测结果;根据指数平滑预测结果得到最终频谱检测结果。
可选地,在本发明提供的频谱检测方法中,历史感知结果序列中包括多个时期的指最终频谱检测结果,对历史感知结果序列进行指数平滑预测,得到指数平滑预测结果的步骤,包括:根据历史感知结果序列获取上一时期的最终频谱检测结果,以及上一时期的平滑值;将上一时期的最终频谱检测结果和上一时期的平滑值输入到预设平滑预测模型中,得到当前时期的指数平滑预测结果。
可选地,在本发明提供的频谱检测方法中,根据历史感知结果序列获取上一时期的平滑值,包括:采用预设迭代公式和各时期的最终频谱检测结果对各时期的平滑值进行迭代计算,得到上一时期的平滑值,预设迭代公式为:Ψ'(l)=β*Ψ(l)+(1-β)*Ψ'(l-1),其中,Ψ'(l)表示l时期的平滑值,Ψ'(l-1)表示l-1时期的平滑值,Ψ(l)表示l时期的最终频谱检测结果,β表示静态平滑指数。
可选地,在本发明提供的频谱检测方法中,预设平滑预测模型为:
Figure BDA0003144345810000031
其中,
Figure BDA0003144345810000032
表示当前时期的指数平滑预测结果,Ψ(Lh)表示上一时期的最终频谱检测结果,Ψ'(Lh)表示上一时期的平滑值,β表示静态平滑指数。
可选地,在本发明提供的频谱检测方法中,根据指数平滑预测结果得到最终频谱检测结果,包括:若指数平滑预测结果小于或等于预设阈值,将最终频谱检测结果确定为不存在频谱信号;若指数平滑预测结果大于预设阈值,将最终频谱检测结果确定为存在频谱信号。
可选地,在本发明提供的频谱检测方法中,获取初始频谱感知结果,包括:采用基于能量检测的双门限检测模型得到检测统计量、第一门限值、第二门限值,将第一门限值与第二门限值之间的区域划分为多个模糊级别;当检测统计量小于第一门限值,或,检测统计量大于第二门限值时,根据检测统计量与第一门限值和第二门限值的比较结果得到一位判决结果;当检测统计量大于或等于第一门限值,且,检测统计量小于或等于第二门限值时,根据检测量统计量所属模糊级别得到两位判决结果。
可选地,在本发明提供的频谱检测方法中,预设模糊判决模型为:
Figure BDA0003144345810000041
其中,Θ表示模糊逻辑输出结果,Medium表示第一预设模糊级别,Low表示第二预设模糊级别,High表示第三预设模糊级别,M2表示两位判决结果的数量,
Figure BDA0003144345810000042
表示向上取整,θi为利用可信度的量化隶属度函数所表示的第i个初始频谱感知结果对应的模糊层次,当检测统计量属于不同的模糊级别时,
Figure BDA0003144345810000043
取不同的值。
可选地,本发明提供的频谱检测方法还包括:若多个初始频谱感知结果中包含一位判决结果,且,一位判决结果不满足预设条件,将最终频谱检测结果确定为存在频谱信号。
可选地,在本发明提供的频谱检测方法中,第二预设模糊级别中对应的检测统计量的值大于或等于第一门限值,第三预设模糊级别中对应的检测统计量的值小于或等于第二门限值;第二预设模糊级别、第一预设模糊级别、第三模糊级别所对应的检测统计量的值依次增大;本发明提供的频谱检测方法还包括:若模糊逻辑输出结果属于第二预设模糊级别,将最终频谱检测结果确定为不存在频谱信号;若模糊逻辑输出结果属于第三预设模糊级别,将最终频谱检测结果确定为存在频谱信号。
本发明第二方面提供了一种频谱检测装置,包括:初始频谱感知结果获取模块,用于获取多个初始频谱感知结果;模糊逻辑输出结果获取模块,若多个初始频谱感知结果中不包含一位判决结果,或,一位判决结果满足预设条件,模糊逻辑输出结果获取模块用于根据初始频谱感知结果中的两位判决结果和预设模糊判决模型得到模糊逻辑输出结果;历史感知结果获取模块,若模糊逻辑输出结果属于第一预设模糊级别,历史感知结果获取模块用于获取历史感知结果序列;指数平滑预测模块,用于对历史感知结果序列进行指数平滑预测,得到指数平滑预测结果;最终频谱检测结果获取模块,用于根据指数平滑预测结果得到最终频谱检测结果。
本发明第三方面提供了一种计算机设备,包括:至少一个处理器;以及与至少一个处理器通信连接的存储器;其中,存储器存储有可被至少一个处理器执行的指令,指令被至少一个处理器执行,从而执行如本发明第一方面提供的频谱检测方法。
本发明第四方面提供了一种计算机可读存储介质,其特征在于,计算机可读存储介质存储有计算机指令,计算机指令用于使计算机执行如本发明第一方面提供的频谱检测方法。
本发明技术方案,具有如下优点:
本发明提供的频谱检测方法及装置,在获取多个初始频谱感知结果后,当初始频谱感知结果中不包含一位判决结果,或,一位判决结果满足预设条件时,根据两位判决结果所属的模糊级别得到模糊逻辑输出结果,当模糊逻辑输出结果属于第一预设模糊级别时,获取历史感知结果序列,并对历史感知结果序列进行指数平滑预测,得到指数平滑预测结果,根据指数平滑预测结果得到最终频谱检测结果,在本发明提供的频谱检测方法及装置中,当无法直接根据模糊逻辑输出结果确定最终频谱检测结果时,结合多个时期的感知结果对当前时期的感知结果进行指数平滑预测,从而得到当前时期的最终频谱检测结果,即使上一时期的最终频谱检测结果存在错误,结合多个时期的感知结果仍然能够得到可靠性更高的当前时期的最终频谱检测结果,因此,本发明提供的频谱检测方法及装置具有较强的鲁棒性,通过实施本发明能够得到更准确的最终频谱检测结果。
附图说明
为了更清楚地说明本发明具体实施方式或现有技术中的技术方案,下面将对具体实施方式或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图是本发明的一些实施方式,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例中协作检测系统模型示意图
图2为本发明实施例中频谱检测方法的一个具体示例的流程图;
图3为本发明实施例中第一门限值和第二门限值之间的不确定区域的模糊级别划分示意图;
图4为本发明实施例中Δi=2dB,β=0.2,Lh=3时,次用户数量M取不同值时得到的检测概率与信噪比曲线;
图5为本发明实施例中M=5,Lh=5时,不同静态平滑指数和不同噪声不确定性下,本发明实施例提供的频谱检测方法对信噪比的检测性能曲线;
图6为本发明实施例中Δi=2dB,β=0.1,M=5时,不同Lh下的检测概率与信噪比曲线;
图7(a)为本发明实施例中Lh=3,M=5,Δi=5dB,β=0.1时,不同噪声不确定度下三种检测算法的检测性能与信噪比的关系;
图7(b)为本发明实施例中Lh=3,M=5,Δi=1dB,β=0.1时,不同噪声不确定度下三种检测算法的检测性能与信噪比的关系;
图7(c)为本发明实施例中Lh=3,M=5,Δi=3dB,β=0.3时,不同噪声不确定度下三种检测算法的检测性能与信噪比的关系;
图7(d)为本发明实施例中Lh=3,M=5,Δi=5dB,β=0.3时,不同噪声不确定度下三种检测算法的检测性能与信噪比的关系;
图8为本发明实施例中频谱检测装置的一个具体示例的原理框图;
图9为本发明实施例中计算机设备的一个具体示例的原理框图。
具体实施方式
下面将结合附图对本发明的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
在本发明的描述中,需要说明的是,术语“第一”、“第二”、“第三”仅用于描述目的,而不能理解为指示或暗示相对重要性。
此外,下面所描述的本发明不同实施方式中所涉及的技术特征只要彼此之间未构成冲突就可以相互结合。
本发明实施例提供了一种频谱检测方法,应用于如图1所示的协作系统模型中的融合中心,在图1所示的协作系统模型中,包括一个主用户(PU),M个次用户(SU)和一个融合中心(FC),在本发明实施例中,每个SU都是相互独立的,感知过程中每个SU将其本地感知结果发送给FC,由FC得到最终的检测结果,确定是否存在PU信号。
本发明实施例提供的频谱检测方法,如图2所示,包括:
步骤S10:获取多个初始频谱感知结果。
初始频谱感知结果是如图1所示的各次用户发送至融合中心的,次用户采用传统基于能量检测的双门限检测模型得到初始频谱感知结果。各次用户在采用传统基于能量检测的双门限检测模型计算过程中,会得到检测统计量Ri、第一门限值λl_i、第二门限值λh_i,当Ril_i,或,Rih_i时,通过检测统计量Ri得到的初始频谱感知结果Di为一位判决结果D1i,当λl_i≤Ri≤λh_i时,通过检测统计量Ri得到的初始频谱感知结果Di为两位判决结果D2i,两位判决结果采用模糊逻辑FL的方式得到。
步骤S20:判断多个初始频谱感知结果中是否包含一位判决结果,当初始频谱感知结果中包含一位判决结果时,判断一位判决结果是否满足预设条件。若初始频谱感知结果中不包含一位判决结果,或,一位判决结果满足预设条件,执行步骤S30。
在一可选实施例中,在获取到多个初始频谱感知结果后,若多个初始频谱感知结果中包含一位判决结果,会先根据一位判决结果进行第一步判决,得到第一步判决结果,第一步判决结果存在两种情况:
Figure BDA0003144345810000081
或,
Figure BDA0003144345810000082
Figure BDA0003144345810000083
表示存在频谱信号,
Figure BDA0003144345810000084
表示不存在频谱信号,当
Figure BDA0003144345810000085
时,可直接将最终频谱检测结果确定为存在频谱信号,但是,当
Figure BDA0003144345810000086
时,需要进一步结合两位判决结果进行第二步判决。
因此,判断一位判决结果是否满足预设条件,实际是判断根据一位判决结果得到的检测结果是否为不存在频谱信号,当一位判决结果满足预设条件时,表示根据一位判决结果得到的检测结果为不存在频谱信号,需要执行步骤S30。当一位判决结果不满足预设条件时,表示根据一位判决结果得到的检测结果为存在频谱信号,此时将最终频谱检测结果确定为存在频谱信号。
步骤S30:根据频谱感知结果中的两位判决结果和预设模糊判决模型得到模糊逻辑输出结果。
如上述步骤S10中记载,当λl_i≤Ri≤λh_i时,得到两位判决结果D2i,两位判决结果位于两门限之间的不确定区域中,不确定区域被划分为多个模糊级别,检测统计量的值不同时,两位判决结果所属的模糊级别也有所不同。
不同的次用户得到的两位判决结果可能属于不同的模糊级别,并且,双门限之间本身就是不确定区域,因此为了保证足够低的虚警概率,在本发明实施例中通过预设模糊判决模型得到各两位判决结果综合的模糊逻辑输出结果,结合模糊逻辑输出结果,得到第二步判决结果
Figure BDA0003144345810000091
模糊逻辑输出结果同样分为多个模糊级别,可以通过一些模糊级别直接得到最终频谱检测结果,但是,仍然存在一些模糊级别,无法直接确定最终频谱检测结果,需要进一步执行第三步判决。
在本发明实施例中,将无法直接得到最终频谱检测结果的模糊级别统称为第一预设模糊级别,执行步骤S30得到模糊逻辑输出结果后,还需执行如下步骤:
步骤S40:判断模糊逻辑输出结果是否属于第一预设模糊级别。
若模糊逻辑输出结果属于第一预设模糊级别,执行第三步判决,即步骤S50-步骤S70;
若模糊逻辑输出结果属于第二预设模糊级别,将最终频谱检测结果确定为不存在频谱信号;若模糊逻辑输出结果属于第三预设模糊级别,将最终频谱检测结果确定为存在频谱信号。
第二预设模糊级别中对应的检测统计量的值大于或等于第一门限值,第三预设模糊级别中对应的检测统计量的值小于或等于所述第二门限值;第二预设模糊级别、第一预设模糊级别、第三模糊级别所对应的检测统计量的值依次增大;。
步骤S50:获取历史感知结果序列。
历史感知结果序列的长度,即,历史感知结果序列中所包含的最终检测结果的数量可以根据实际需求进行设定,示例性地,历史感知结果序列中可以包括距离当前时期最近的5个时期的最终频谱检测结果,也可以包括距离当前时期最近的7个时期的最终频谱检测结果。
历史感知结果序列是实时更新的,在每次获得当前时期最终频谱检测结果后,会在历史感知结果序列中加入当前时期的最终频谱检测结果,从历史感知结果序列中删除距离当前时期最远的时期的最终频谱检测结果,从而完成对历史感知结果序列的更新。
步骤S60:对历史感知结果序列进行指数平滑预测,得到指数平滑预测结果。
在一可选实施例中,对历史感知结果序列进行指数平滑预测时,先确定平滑初始值,结合平滑初始值和历史感知结果序列中各时期的最终频谱检测结果进行迭代计算,得到上一时期的平滑值,然后将上一时期的平滑值和上一时期的最终频谱检测结果输入到平滑预测模型中计算得到指数平滑预测结果。
步骤S70:根据指数平滑预测结果得到最终频谱检测结果。
在一可选实施例中,若指数平滑预测结果小于或等于预设阈值,将最终频谱检测结果确定为不存在频谱信号;若指数平滑预测结果大于预设阈值,将最终频谱检测结果确定为存在频谱信号。
示例性地,可以将预设阈值设置为0.5,即,当指数平滑预测结果小于或等于0.5时,将最终频谱检测结果确定为不存在频谱信号;当指数平滑预测结果大于0.5时,将最终频谱检测结果确定为存在频谱信号:
Figure BDA0003144345810000111
本发明实施例提供的频谱检测方法,在获取多个初始频谱感知结果后,当初始频谱感知结果中不包含一位判决结果,或,一位判决结果满足预设条件时,根据两位判决结果所属的模糊级别得到模糊逻辑输出结果,当模糊逻辑输出结果属于第一预设模糊级别时,获取历史感知结果序列,并对历史感知结果序列进行指数平滑预测,得到指数平滑预测结果,根据指数平滑预测结果得到最终频谱检测结果。在本发明实施例提供的频谱检测方法中,当无法直接根据模糊逻辑输出结果确定最终频谱检测结果时,结合多个时期的感知结果对当前时期的感知结果进行指数平滑预测,从而得到当前时期的最终频谱检测结果,即使上一时期的最终频谱检测结果存在错误,结合多个时期的感知结果仍然能够得到可靠性更高的当前时期的最终频谱检测结果,因此,本发明实施例提供的频谱检测方法具有较强的鲁棒性,通过实施本发明实施例能够得到更准确的最终频谱检测结果。
在一可选实施例中,在上述步骤S60中,计算上一时期的平滑值的迭代公式为:
Ψ'(l)=β*Ψ(l)+(1-β)*Ψ'(l-1),
其中,Ψ'(l)表示l时期的平滑值,Ψ'(l-1)表示l-1时期的平滑值,Ψ(l)表示l时期的最终频谱检测结果,β表示静态平滑指数,在一可选实施例中,可以设置平滑初始值为Ψ'(0)=(∑Ψ)/Lh,Lh表示历史感知结果序列长度。
在一可选实施例中,在上述步骤S60中,用于计算指数平滑预测结果的预设平滑预测模型为:
Figure BDA0003144345810000121
其中,
Figure BDA0003144345810000122
表示当前时期的指数平滑预测结果,Ψ(Lh)表示上一时期的最终频谱检测结果,Ψ'(Lh)表示上一时期的平滑值,β表示静态平滑指数。
在一可选实施例中,次用户(SU)获取初始频谱感知结果时,检测过程为:
Figure BDA0003144345810000123
其中,H0和H1分别表示信道内不存在频谱信号和存在频谱信号,ni(t)是服从正态分布的加性高斯白噪声,均值为0,方差为
Figure BDA0003144345810000124
si(t)是主用户信号,均值为μ方差为
Figure BDA0003144345810000125
di(t)是加入的强度为ξ的直流噪声。
次用户根据接收到的K个样本的能量值确定检测统计量:
Figure BDA0003144345810000126
在K值足够大的情况下,Ri近似服从高斯分布:
Figure BDA0003144345810000127
其中,
Figure BDA0003144345810000128
第i个次用户的检测概率Pdi与虚警概率Pfi可分别表示为:
Figure BDA0003144345810000129
其中,
Figure BDA00031443458100001210
当虚警概率一定时,可以计算得到单门限λi
Figure BDA0003144345810000131
其中,ξ的大小为
Figure BDA0003144345810000132
假定噪声的不确定度为ΔidB,则
Figure BDA0003144345810000133
则实际噪声方差的范围为
Figure BDA0003144345810000134
Figure BDA0003144345810000135
基于噪声不确定度的第一门限值λl_i和第二门限值λh_i可表示为:
Figure BDA0003144345810000136
则传统的能量检测的双门限检测结果可表示为:
Figure BDA0003144345810000137
其中,ND表示无判决结果。
当Ril_i,或,Rih_i时,得到一位判决结果D1i
Figure BDA0003144345810000138
当检测量Ri位于两门限之间时,采用模糊逻辑FL的方式得到两位判决结果D2i。模糊系统使用一组语言规则来表示模糊情况,利用模糊隶属度函数将这些规则转化为数学表达式。第一门限值和第二门限值之间的不确定区域的划分方式可根据实际需求进行调整,在如图3所示的实施例中,将第一门限值和第二门限值之间的不确定区域划分为三个模糊级别:Low、Medium、High,其中,λ1_i=(2λl_ih_i)/3,λ2_i=(2λh_il_i)/3,当λl_i≤Ri1_i时,表示两位判决结果属于模糊级别Low,当λ1_i≤Ri≤λ2_i时,表示两位判决结果属于模糊级别Medium,当λ2_i<Ri≤λh_i时,表示两位判决结果属于模糊级别High。
利用可信度的量化隶属度函数来表示各个模糊层次,用θ表示:
Figure BDA0003144345810000141
θi=-1表示主用户有很大可能不存在,θi=1表示主用户很有可能存在,θi=0表示无法判断的情况。
在一可选实施例中,三种不同的模糊级别可以用两个bit进行表示,从而得到D2i
Figure BDA0003144345810000142
综上,次用户获取的初始频谱感知结果可表示为:
Figure BDA0003144345810000143
在一可选实施例中,如上述步骤S10中记载,在获取到多个初始频谱感知结果后,若多个初始频谱感知结果中包含一位判决结果,会先根据一位判决结果进行第一步判决,示例性地,可以使用OR准则结合进行第一步判决,得到第一步判决检测结果
Figure BDA0003144345810000144
Figure BDA0003144345810000145
其中,M1表示一位判决结果的数量,D1i表示第i个一位判决结果。
当通过本发明实施例计算第一步判决检测结果
Figure BDA0003144345810000146
时,步骤S20中所称的预设条件为
Figure BDA0003144345810000147
在一可选实施例中,在上述步骤S30中,获取模糊逻辑输出结果时所使用的预设模糊判决模型可以根据模糊判决规则确定,示例性地,模糊逻辑规则可以为:
表1
SU 1 SU 2 SU 3 Out
Low Low Low Low
Low Low Medium Low
Low Low High Low
Low Medium Medium Low
Low Medium High Medium
Low High High Medium
Medium Medium Medium Medium
Medium Medium High Medium
Medium High High High
High High High High
根据表1的模糊逻辑规则得到的预设模糊判决模型为:
Figure BDA0003144345810000151
在一可选实施例中,上述步骤S30中所称的结合模糊逻辑输出结果得到的第二步判决结果
Figure BDA0003144345810000152
可以为:
Figure BDA0003144345810000153
由上式可知,当模糊逻辑结果所属的模糊级别为Low时,可以得到最终频谱检测结果为不存在频谱信号,当模糊逻辑结果所属的模糊级别为High时,可以得到最终频谱检测结果为存在频谱信号,但是,当模糊逻辑结果所属的模糊级别为Medium时,无法直接得到最终频谱检测结果
Figure BDA0003144345810000154
需要进行三步判决。在执行本发明实施例得到第二步判决结果时,上述步骤S40中的第一预设模糊级别为Medium,第二预设模糊级别为Low,第三预设模糊级别为High。
在本发明实施例提供的技术方案中,当模糊逻辑结果所属的模糊级别为Medium时,执行上述步骤S50-步骤S70,对历史感知结果序列进行指数平滑预测,基于预测结果得到当前时期的最终频谱检测结果,但是,在背景技术所记载的对比技术中,当模糊逻辑结果所属的模糊级别为Medium时,直接将上一时期的最终频谱检测结果确定为当前时期的最终频谱检测结果,相较于背景技术中记载的对比技术,通过实施本发明实施例提供的技术方案得到的当前时期的最终频谱检测结果更准确。
本发明实施例提供的频谱检测方法,采用基于双阈值的模糊逻辑对局部检测中的不确定区域进行量化,提高了在噪声不确定环境下的检测性能,并且实施本发明可以在不增加传输开销的基础上提高检测性能。
在一可选实施例中,对提出的算法进行蒙特卡洛仿真,假设在传输过程中没有比特错误,每一个SU的仿真参数如下表2所示,
表2
Figure BDA0003144345810000161
如图4所示为Δi=2dB,β=0.2,Lh=3时,次用户数量M取不同值时得到的检测概率与信噪比曲线。
从图4中可以看出,随着M的增大,检测概率逐渐增大。但是,当M增加到一定程度时,检测性能的效率会逐渐降低,随之,协同网络的系统开销仍然会线性增加。为了平衡检测性能和系统开销,取M的值为5。
如图5所示为M=5,Lh=5时,不同静态平滑指数和不同噪声不确定性下,本发明实施例提供的频谱检测方法对信噪比的检测性能。从图5可以看出,随着Δ和β的变化,本发明实施例提供的频谱检测方法的虚警概率(PF)基本保持不变,而检测概率(PD)的变化则不同。
在低噪声不确定度下,随着的β增加,检测概率略有降低。当Δi较大时,在低信噪比环境下,β=0.3的检测概率优于β=0.1,当信噪比增加到一定值时,β=0.1的检测概率优于β=0.3。这是因为当β小于0.3时,一阶指数平滑预测适用于预测数据变化平缓的场景,而当β大于0.3时,一阶指数平滑预测适用于预测数据变化剧烈的场景。在低信噪比环境下,当噪声不确定性较大时,检测结果具有较大的波动性,因此β=0.3时的检测概率优于β=0.1。当信噪比较好或噪声不确定度较小时,检测结果的精度较高。此时数据变化平稳,所以β=0.1时的检测概率优于β=0.3,与理论结果一致。当Δi大于3dB时,设置β=0.3,否则设置β=0.1。
如图6所示为Δi=2dB,β=0.1,M=5时,不同Lh下的检测概率与信噪比曲线。从图6可以看出,随着Lh值的增加,检测概率增加不大。这是因为一阶指数平滑预测是短期的预测,所以增加历史传感数据的长度并不能提高检测性能。从图6可以看出,与Lh=2相比,当Lh=3时,检测性能更平稳,因此选择Lh=3。
为了验证本发明实施例提供的频谱检测方法的有效性,本发明实施例中将本发明实施例提供的频谱检测算法与背景技术中提出的对比技术以及传统的能量检测算法进行了比较。为便于描述,将本发明实施例提供的频谱检测方法标注为“proposed”,将背景技术中提出的对比技术标注为“reference”,将传统的能量检测算法标注为“ED”。
图7(a)、图7(b)、图7(c)、图7(d)为不同噪声不确定度下三种检测算法的检测性能与信噪比的关系,图7(a)、图7(b)、图7(c)、图7(d)中,Lh的值都为3,M的值都为5,在图7(a)中,Δi=5dB,β=0.1,在图7(b)中,Δi=1dB,β=0.1,在图7(c)中,Δi=3dB,β=0.3,在图7(d)中,Δi=5dB,β=0.3。
由图7(a)、图7(b)、图7(c)、图7(d)对比可知,在噪声不确定度较低的情况下,在不增加虚警概率的情况下,本发明实施例提供的频谱检测方法的检测性能(detectionprefprmance)优于其他算法,进而验证了通过执行本发明实施例,可以预测得到准确性更高的最终频谱检测结果。在噪声不确定性较大的情况下,该算法还可以牺牲微弱的虚警概率来提高检测性能。
通过分析可以看出,本发明实施例提供的频谱检测方法的系统传输消耗与背景技术中提出的对比技术相同,在双阈值内均为2比特。
通过计算可知,本发明实施例提供的频谱检测方法的计算量主要包括M*K+2Lh个乘法运算和M*K+M+Lh个加法运算,背景技术中提出的对比技术的计算量主要包括M*K+M2个乘法运算和M*K+M个加法运算。然而,由于Lh=3,所以本发明实施例提供的频谱检测方法的计算复杂度近似等于背景技术中提出的对比技术的计算复杂度,且本发明实施例提供的频谱检测方法的存储消耗可以忽略不计。因此,与其他两种算法相比,本发明实施例提供的频谱检测方法具有一定的优势。
本发明实施例提供了一种频谱检测装置,如图8所示,包括:
初始频谱感知结果获取模块10,用于获取多个初始频谱感知结果,详细内容参见上述实施例中对步骤S10的描述,在此不再赘述。
模糊逻辑输出结果获取模块20,若多个初始频谱感知结果中不包含一位判决结果,或,一位判决结果满足预设条件,模糊逻辑输出结果获取模块20用于根据初始频谱感知结果中的两位判决结果和预设模糊判决模型得到模糊逻辑输出结果,详细内容参见上述实施例中对步骤S20、步骤S30的描述,在此不再赘述。
历史感知结果获取模块30,若模糊逻辑输出结果属于第一预设模糊级别,历史感知结果获取模块30用于获取历史感知结果序列,详细内容参见上述实施例中对步骤S40、步骤S50的描述,在此不再赘述。
指数平滑预测模块40,用于对历史感知结果序列进行指数平滑预测,得到指数平滑预测结果,详细内容参见上述实施例中对步骤S60的描述,在此不再赘述。
最终频谱检测结果获取模块50,用于根据指数平滑预测结果得到最终频谱检测结果,详细内容参见上述实施例中对步骤S70的描述,在此不再赘述。
本发明实施例提供的频谱检测装置,在获取多个初始频谱感知结果后,当初始频谱感知结果中不包含一位判决结果,或,一位判决结果满足预设条件时,根据两位判决结果所属的模糊级别得到模糊逻辑输出结果,当模糊逻辑输出结果属于第一预设模糊级别时,获取历史感知结果序列,并对历史感知结果序列进行指数平滑预测,得到指数平滑预测结果,根据指数平滑预测结果得到最终频谱检测结果,在本发明实施例提供的频谱检测装置中,当无法直接根据模糊逻辑输出结果确定最终频谱检测结果时,结合多个时期的感知结果对当前时期的感知结果进行指数平滑预测,从而得到当前时期的最终频谱检测结果,即使上一时期的最终频谱检测结果存在错误,结合多个时期的感知结果仍然能够得到可靠性更高的当前时期的最终频谱检测结果,因此,本发明提供的频谱检测装置具有较强的鲁棒性,通过实施本发明能够得到更准确的最终频谱检测结果。
本发明实施例提供了一种计算机设备,如图9所示,该计算机设备主要包括一个或多个处理器61以及存储器62,图9中以一个处理器61为例。
该计算机设备还可以包括:输入装置63和输出装置64。
处理器61、存储器62、输入装置63和输出装置64可以通过总线或者其他方式连接,图9中以通过总线连接为例。
处理器61可以为中央处理器(Central Processing Unit,CPU)。处理器61还可以为其他通用处理器、数字信号处理器(Digital Signal Processor,DSP)、专用集成电路(Application Specific Integrated Circuit,ASIC)、现场可编程门阵列(Field-Programmable Gate Array,FPGA)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件等芯片,或者上述各类芯片的组合。通用处理器可以是微处理器或者该处理器也可以是任何常规的处理器等。存储器62可以包括存储程序区和存储数据区,其中,存储程序区可存储操作系统、至少一个功能所需要的应用程序;存储数据区可存储根据频谱检测装置的使用所创建的数据等。此外,存储器62可以包括高速随机存取存储器,还可以包括非暂态存储器,例如至少一个磁盘存储器件、闪存器件、或其他非暂态固态存储器件。在一些实施例中,存储器62可选包括相对于处理器61远程设置的存储器,这些远程存储器可以通过网络连接至频谱检测装置。输入装置63可接收用户输入的计算请求(或其他数字或字符信息),以及产生与频谱检测装置有关的键信号输入。输出装置64可包括显示屏等显示设备,用以输出计算结果。
本发明实施例提供了一种计算机可读存储介质,该计算机可读存储介质存储计算机指令,计算机存储介质存储有计算机可执行指令,该计算机可执行指令可执行上述任意方法实施例中的频谱检测方法。其中,存储介质可为磁碟、光盘、只读存储记忆体(Read-Only Memory,ROM)、随机存储记忆体(Random Access Memory,RAM)、快闪存储器(FlashMemory)、硬盘(Hard Disk Drive,缩写:HDD)或固态硬盘(Solid-State Drive,SSD)等;存储介质还可以包括上述种类的存储器的组合。
显然,上述实施例仅仅是为清楚地说明所作的举例,而并非对实施方式的限定。对于所属领域的普通技术人员来说,在上述说明的基础上还可以做出其它不同形式的变化或变动。这里无需也无法对所有的实施方式予以穷举。而由此所引伸出的显而易见的变化或变动仍处于本发明创造的保护范围之中。

Claims (12)

1.一种频谱检测方法,其特征在于,包括:
获取多个初始频谱感知结果;
若多个初始频谱感知结果中不包含一位判决结果,或,所述一位判决结果满足预设条件,根据所述初始频谱感知结果中的两位判决结果和预设模糊判决模型得到模糊逻辑输出结果;
若所述模糊逻辑输出结果属于第一预设模糊级别,
获取历史感知结果序列;
对所述历史感知结果序列进行指数平滑预测,得到指数平滑预测结果;
根据指数平滑预测结果得到最终频谱检测结果。
2.根据权利要求1所述的频谱检测方法,其特征在于,所述历史感知结果序列中包括多个时期的最终频谱检测结果,
对所述历史感知结果序列进行指数平滑预测,得到指数平滑预测结果的步骤,包括:
根据所述历史感知结果序列获取上一时期的最终频谱检测结果,以及上一时期的平滑值;
将所述上一时期的最终频谱检测结果和所述上一时期的平滑值输入到预设平滑预测模型中,得到当前时期的指数平滑预测结果。
3.根据权利要求2所述的频谱检测方法,其特征在于,根据所述历史感知结果序列获取上一时期的平滑值,包括:
采用预设迭代公式和各时期的最终频谱检测结果对各时期的平滑值进行迭代计算,得到所述上一时期的平滑值,所述预设迭代公式为:
Ψ'(l)=β*Ψ(l)+(1-β)*Ψ'(l-1),
其中,Ψ'(l)表示l时期的平滑值,Ψ'(l-1)表示l-1时期的平滑值,Ψ(l)表示l时期的最终频谱检测结果,β表示静态平滑指数。
4.根据权利要求2或3所述的频谱检测方法,其特征在于,所述预设平滑预测模型为:
Figure FDA0003144345800000021
其中,
Figure FDA0003144345800000022
表示当前时期的指数平滑预测结果,Ψ(Lh)表示上一时期的最终频谱检测结果,Ψ'(Lh)表示上一时期的平滑值,β表示静态平滑指数。
5.根据权利要求1-4中任一项所述的频谱检测方法,其特征在于,根据指数平滑预测结果得到最终频谱检测结果,包括:
若所述指数平滑预测结果小于或等于预设阈值,将所述最终频谱检测结果确定为不存在频谱信号;
若所述指数平滑预测结果大于所述预设阈值,将所述最终频谱检测结
果确定为存在所述频谱信号。
6.根据权利要求1所述的频谱检测方法,其特征在于,获取初始频谱感知结果,包括:
采用基于能量检测的双门限检测模型得到检测统计量、第一门限值、第二门限值,将第一门限值与第二门限值之间的区域划分为多个模糊级别;
当所述检测统计量小于第一门限值,或,所述检测统计量大于所述第二门限值时,根据所述检测统计量与第一门限值和第二门限值的比较结果得到一位判决结果;
当所述检测统计量大于或等于所述第一门限值,且,所述检测统计量小于或等于所述第二门限值时,根据所述检测量统计量所属模糊级别得到两位判决结果。
7.根据权利要求6所述的频谱检测方法,其特征在于,预设模糊判决模型为:
Figure FDA0003144345800000031
其中,Θ表示模糊逻辑输出结果,Medium表示第一预设模糊级别,Low表示第二预设模糊级别,High表示第三预设模糊级别,M2表示两位判决结果的数量,
Figure FDA0003144345800000032
表示向上取整,θi为利用可信度的量化隶属度函数所表示的第i个初始频谱感知结果对应的模糊层次,当检测统计量属于不同的模糊级别时,θi取不同的值。
8.根据权利要求6所述的频谱检测方法,其特征在于,还包括:
若所述多个初始频谱感知结果中包含一位判决结果,且,所述一位判决结果不满足预设条件,将所述最终频谱检测结果确定为存在频谱信号。
9.根据权利要求7所述的频谱检测方法,其特征在于,
所述第二预设模糊级别中对应的检测统计量的值大于或等于所述第一门限值,所述第三预设模糊级别中对应的检测统计量的值小于或等于所述第二门限值;第二预设模糊级别、第一预设模糊级别、第三模糊级别所对应的检测统计量的值依次增大;
所述方法还包括:
若所述模糊逻辑输出结果属于第二预设模糊级别,将所述最终频谱检测结果确定为不存在频谱信号;
若所述模糊逻辑输出结果属于第三预设模糊级别,将所述最终频谱检测结果确定为存在频谱信号。
10.一种频谱检测装置,其特征在于,包括:
初始频谱感知结果获取模块,用于获取多个初始频谱感知结果;
模糊逻辑输出结果获取模块,若多个初始频谱感知结果中不包含一位判决结果,或,所述一位判决结果满足预设条件,所述模糊逻辑输出结果获取模块用于根据所述初始频谱感知结果中的两位判决结果和预设模糊判决模型得到模糊逻辑输出结果;
历史感知结果获取模块,若所述模糊逻辑输出结果属于第一预设模糊级别,所述历史感知结果获取模块用于获取历史感知结果序列;
指数平滑预测模块,用于对所述历史感知结果序列进行指数平滑预测,得到指数平滑预测结果;
最终频谱检测结果获取模块,用于根据指数平滑预测结果得到最终频谱检测结果。
11.一种计算机设备,其特征在于,包括:
至少一个处理器;以及与所述至少一个处理器通信连接的存储器;其中,所述存储器存储有可被所述至少一个处理器执行的指令,所述指令被所述至少一个处理器执行,从而执行如权利要求1-9中任一项所述的频谱检测方法。
12.一种计算机可读存储介质,其特征在于,所述计算机可读存储介质存储有计算机指令,所述计算机指令用于使所述计算机执行如权利要求1-9中任一项所述的频谱检测方法。
CN202110745871.6A 2021-03-15 2021-07-01 一种频谱检测方法及装置 Active CN113765607B (zh)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
CN202110273690 2021-03-15
CN2021102736908 2021-03-15

Publications (2)

Publication Number Publication Date
CN113765607A CN113765607A (zh) 2021-12-07
CN113765607B true CN113765607B (zh) 2022-08-30

Family

ID=78787531

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110745871.6A Active CN113765607B (zh) 2021-03-15 2021-07-01 一种频谱检测方法及装置

Country Status (1)

Country Link
CN (1) CN113765607B (zh)

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101710848A (zh) * 2009-11-13 2010-05-19 北京邮电大学 基于模糊积分和最优化理论的合作式频谱感知方法
CN103763706A (zh) * 2013-12-30 2014-04-30 河海大学 基于分簇的协作频谱感知模糊融合方法
CN103856946A (zh) * 2014-01-13 2014-06-11 南京邮电大学 基于差分能量检测的双门限协作频谱感知方法
CN105578480A (zh) * 2015-12-15 2016-05-11 西安电子科技大学 面向宽带调制变换器的欠采样频谱感知预判决方法
CN107592174A (zh) * 2017-09-12 2018-01-16 河南工业大学 一种智能电网通信中的高效频谱感知方法
CN107623554A (zh) * 2017-08-24 2018-01-23 西安电子科技大学 一种基于信任值的协作频谱感知方法、认知无线网络
CN109272142A (zh) * 2018-08-14 2019-01-25 苏州朗动网络科技有限公司 基于事件的数据预测方法、装置、计算机设备和存储介质
CN109347584A (zh) * 2018-11-02 2019-02-15 西安电子科技大学昆山创新研究院 基于指数平滑预测的频谱感知方法
CN111342921A (zh) * 2020-03-06 2020-06-26 西京学院 基于认知无线电的跳频通信中频谱共享方法、系统及应用
CN112073138A (zh) * 2020-09-17 2020-12-11 重庆大学 一种基于量化的双门限协作频谱感知方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10651955B2 (en) * 2018-01-19 2020-05-12 Arizona Board Of Regents On Behalf Of Arizona State University Active sequential xampling receiver for spectrum sensing

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101710848A (zh) * 2009-11-13 2010-05-19 北京邮电大学 基于模糊积分和最优化理论的合作式频谱感知方法
CN103763706A (zh) * 2013-12-30 2014-04-30 河海大学 基于分簇的协作频谱感知模糊融合方法
CN103856946A (zh) * 2014-01-13 2014-06-11 南京邮电大学 基于差分能量检测的双门限协作频谱感知方法
CN105578480A (zh) * 2015-12-15 2016-05-11 西安电子科技大学 面向宽带调制变换器的欠采样频谱感知预判决方法
CN107623554A (zh) * 2017-08-24 2018-01-23 西安电子科技大学 一种基于信任值的协作频谱感知方法、认知无线网络
CN107592174A (zh) * 2017-09-12 2018-01-16 河南工业大学 一种智能电网通信中的高效频谱感知方法
CN109272142A (zh) * 2018-08-14 2019-01-25 苏州朗动网络科技有限公司 基于事件的数据预测方法、装置、计算机设备和存储介质
CN109347584A (zh) * 2018-11-02 2019-02-15 西安电子科技大学昆山创新研究院 基于指数平滑预测的频谱感知方法
CN111342921A (zh) * 2020-03-06 2020-06-26 西京学院 基于认知无线电的跳频通信中频谱共享方法、系统及应用
CN112073138A (zh) * 2020-09-17 2020-12-11 重庆大学 一种基于量化的双门限协作频谱感知方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
基于信道认知在线可定义的电力线载波通信方法;李建岐,陆阳,高鸿坚;《中国电机工程学报》;20151020;第35卷(第20期);5235-5242 *

Also Published As

Publication number Publication date
CN113765607A (zh) 2021-12-07

Similar Documents

Publication Publication Date Title
US8918375B2 (en) Content aware chunking for achieving an improved chunk size distribution
CN112287986B (zh) 一种图像处理方法、装置、设备及可读存储介质
US11010451B2 (en) Techniques for automated Bayesian posterior sampling using Markov Chain Monte Carlo and related schemes
CN111326169B (zh) 一种语音质量的评价方法及装置
CN113014340A (zh) 一种基于神经网络的卫星频谱资源动态分配方法
CN113270107A (zh) 音频信号中噪声响度的获取方法、装置和电子设备
CN114676825A (zh) 一种神经网络模型量化方法、系统、设备及介质
CN113657249B (zh) 训练方法、预测方法、装置、电子设备以及存储介质
CN112101543A (zh) 神经网络模型确定方法、装置、电子设备及可读存储介质
CN109219055B (zh) 一种主用户占空比估计方法
CN113765607B (zh) 一种频谱检测方法及装置
CN112419312A (zh) 相似房源信息检测方法、装置、电子设备和可读介质
CN112073138B (zh) 一种基于量化的双门限协作频谱感知方法
CN115984084A (zh) 一种基于动态可拆分网络的遥感分布式数据处理方法
CN115941084A (zh) 基于时频图模板匹配的水声通信前导信号检测方法及装置
CN112003662B (zh) 认知网络中的基于降维和聚类的协作频谱感知方法、装置
Wang et al. Estimation and confidence bands of a conditional survival function with censoring indicators missing at random
CN108833138B (zh) 一种错误累积敏感的增量式动态社区发现方法及系统
CN112148902A (zh) 数据处理方法、装置、服务器及存储介质
CN112668832A (zh) 基于指标管理体系的风险量化评估方法、装置以及电子设备
CN106685549B (zh) 一种主用户频谱感知方法与装置
CN109951237B (zh) 一种信号检测的方法、装置及计算机存储介质
WO2018203510A1 (ja) 質問推定装置
WO2020037629A1 (en) Co-operative sensing
KR20190031982A (ko) 인지 무선 통신 방법

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