基于动态阈值和残差卷积网络的电磁频谱占用度二维预测方法

1.本发明涉及一种电磁频谱占用度二维预测方法，属于电磁环境频谱占用预测技术领域。


背景技术：

2.随着信息和通信领域的迅猛发展，对复杂电磁环境的研究在进一步推动了人类社会和技术进步，与此同时，众多用频设备的快速扩展使得受限的频谱资源日益稀缺和拥塞。特别是在原有的静态频谱分配策略中，主要用户占用特定频段工作引起的频谱空洞会造成频谱资源的浪费和低效，这导致了频谱资源的可持续风险增加。因此，频谱预测技术吸引了很多关注，其通过感知频谱空洞的发生，使次用户能够动态的接入频谱空洞，在不对主用户的活动产生干扰的情况下，使更多的用户可以在有限的频谱范围内有序高效地工作，从而提升频谱利用率，降低频谱资源的可持续风险。
3.近年来深度学习技术被引入到电磁环境分析挖掘领域，rnn、seq2seq和lstm等深度学习网络已被广泛应用于频谱预测，并取得了较好的预测性能。但是现有研究中普遍存在两大劣势，一是忽视了电磁频谱先验信息的引入，电磁频谱由于其不同于一般时间序列的特性，往往具有时间近邻性、周期性和趋势性特征，而本方法的多维张量数据在频谱预测中融入了特征先验信息，提升了网络预测性能；二是预测维度单一，现有研究中很少有频谱预测技术能够同时预测时频二维的频谱占用情况，单频点或单步长的预测应用价值较低，而本方法实现了多尺度频谱占用度图像预测，具有较高的用频指导意义。


技术实现要素：

4.本发明为解决目前对于电磁频谱占用度预测任务上先验信息少、预测尺度单一、预测性能差的问题，进而提出一种基于动态阈值和残差卷积网络的电磁频谱占用度二维预测方法。
5.本发明为解决上述问题采取的技术方案是：本发明的步骤如下：
6.步骤一、对原始频谱数据进行数据缺失和异常处理，获得预处理后的频谱数据；
7.步骤二、对预处理后的电磁频谱数据，采用动态阈值算法提取电磁频谱底噪数据并建立电磁频谱底噪拟合模型；
8.步骤三、在电磁频谱底噪模型上加上3db作为电磁频谱占用状态阈值，得到电磁频谱占用状态数据，并计算每15分钟的频谱占用度；
9.步骤四、对频谱占用度数据采用三个多尺度滑动窗口获得时频二维的频谱占用度数据，并以多维张量的形式构建时频块；
10.步骤五、将tfblock输入到结合注意机制的残差卷积网络和lstm网络中，实现多尺度电磁频谱占用度图像预测。
11.进一步的，步骤二中对电磁频谱数据进行动态阈值的底噪数据提取，电磁频谱底
噪的自适应阈值模型表达式如下：
[0012][0013]
公式(1)中nk是在对频谱数据进行排序之后，原始频谱数据行中具有k1比率的时隙的数量，m(j)是k1比率时隙的功率谱密度值之和，k2是功率谱密度数据范围的权重，f
max
和f
min
表示最大和最小功率谱密度值；
[0014]
在确定电磁环境的动态噪声阈值之后，将低于阈值的数据拟合为电磁频谱底噪，选择的拟合模型为正交级数耦合模型，电磁环境噪声模型的统一表达式为：
[0015][0016]
公式(2)中c0为常数，ξj为常数，ζj为常数，n表示级数阶数，j表示虚数单位，ρ表示频点序列，w表示级数角频率。
[0017]
进一步的，步骤三中将建模后的底噪阈值加上3db作为频谱占用状态阈值，得到频谱占用数据，电磁频谱占用状态阈值表达式如下：
[0018][0019][0020]
公式(3)和(4)中，表示电磁频谱占用状态阈值，s
i,j
表示第i个时隙第j个频点的功率值，x
i,j
表示频谱占用状态。
[0021]
进一步的，步骤四中采用三个多尺度滑动窗口获取时频二维的频谱占用度数据，分别建模电磁频谱占用度的时间近邻性、周期性和趋势性，构建多维张量时频块。
[0022]
进一步的，步骤五中将时频块的三个分量输入到具有se注意机制的卷积残差网络中获取二维特征信息，将二维特征展平为一维特征输入到lstm后，基于矩阵对三个分量进行融合，最后通过tanh激活函数,实现多尺度电磁频谱占用度图像预测。
[0023]
本发明的有益效果是：
[0024]
1、本发明采用多尺度滑动窗口充分提取频谱占用度的时间近邻性、周期性和趋势性，包含丰富的先验信息，并采用基于矩阵的多特征融合方法，实现了从长期尺度上挖掘深层电磁频谱相关信息提升预测性能；
[0025]
2、本发明设计的网络模型结合了深度残差网络的深层频谱数据挖掘特征和lstm网络的时间序列记忆特征，充分挖掘了电磁频谱的图像特征和时间序列特征，预测性能好，工程可实现性强；
[0026]
3、本发明的核心在于设计了动态阈值的底噪和占用度提取算法，采用基于电磁频谱占用度时间近邻特征、周期特征和趋势特征作为先验信息，通过tfbrl网络融合残差卷积网络和lstm时序预测网络，实现电磁频谱占用度预测，有效提升预测性能。
附图说明
[0027]
图1是本发明的流程框图；
[0028]
图2是多尺度滑动窗口构建多维张量数据示意图；
[0029]
图3是tfbrl网络结构图；
[0030]
图4是预测误差值mse随预测尺度的变化示意图；
[0031]
图5是预测误差值rmse随预测尺度的变化示意图；
[0032]
图6是预测误差值mae随预测尺度的变化示意图。
具体实施方式
[0033]
具体实施方式一：结合图1说明本实施方式，本实施方式所述基于动态阈值和残差卷积网络的电磁频谱占用度二维预测方法的步骤如下：
[0034]
步骤一、对原始频谱数据进行数据缺失和异常处理，获得预处理后的频谱数据；
[0035]
步骤二、对预处理后的电磁频谱数据，采用动态阈值算法提取电磁频谱底噪数据并建立电磁频谱底噪拟合模型；
[0036]
步骤三、在电磁频谱底噪模型上加上3db作为电磁频谱占用状态阈值，得到电磁频谱占用状态数据，并计算每15分钟的频谱占用度；
[0037]
步骤四、对频谱占用度数据采用三个多尺度滑动窗口获得时频二维的频谱占用度数据，并以多维张量的形式构建时频块；
[0038]
步骤五、将tfblock输入到结合注意机制的残差卷积网络和lstm网络中，实现多尺度电磁频谱占用度图像预测。
[0039]
本实施方式中，所述对公开的电磁频谱监测数据芬兰图尔库数据集进行预处理，电磁频谱数据可以表示为：
[0040]
s＝(s1,s2,s3,s4,...,sn),s∈rm×n[0041]
此处m为监测时长，共计m个时隙，n是监测频率点的数量，s1表示m个时隙中的频点1获得的功率谱值向量。
[0042]
具体实施方式二：结合图1说明本实施方式，本实施方式所述基于动态阈值和残差卷积网络的电磁频谱占用度二维预测方法的步骤二中对电磁频谱数据进行动态阈值的底噪数据提取，电磁频谱底噪的自适应阈值模型表达式如下：
[0043][0044]
公式(1)中nk是在对频谱数据进行排序之后，原始频谱数据行中具有k1比率的时隙的数量，m(j)是k1比率时隙的功率谱密度值之和，k2是功率谱密度数据范围的权重，f
max
和f
min
表示最大和最小功率谱密度值；
[0045]
在确定电磁环境的动态噪声阈值之后，将低于阈值的数据拟合为电磁频谱底噪，选择的拟合模型为正交级数耦合模型，电磁环境噪声模型的统一表达式为：
[0046][0047]
公式(2)中c0为常数，ξj为常数，ζj为常数，n表示级数阶数，j表示虚数单位，ρ表示频点序列，w表示级数角频率。
[0048]
具体实施方式三：结合图1说明本实施方式，本实施方式所述基于动态阈值和残差卷积网络的电磁频谱占用度二维预测方法的步骤三中将建模后的底噪阈值加上3db作为频
谱占用状态阈值，得到频谱占用数据，电磁频谱占用状态阈值表达式如下：
[0049][0050][0051]
公式(3)和(4)中，表示电磁频谱占用状态阈值，s
i,j
表示第i个时隙第j个频点的功率值，x
i,j
表示频谱占用状态。
[0052]
具体实施方式四：结合图1说明本实施方式，本实施方式所述基于动态阈值和残差卷积网络的电磁频谱占用度二维预测方法的步骤四中采用三个多尺度滑动窗口获取时频二维的频谱占用度数据，分别建模电磁频谱占用度的时间近邻性、周期性和趋势性，构建多维张量时频块。
[0053]
具体实施方式五：结合图1说明本实施方式，本实施方式所述基于动态阈值和残差卷积网络的电磁频谱占用度二维预测方法的步骤五中将时频块的三个分量输入到具有se注意机制的卷积残差网络中获取二维特征信息，将二维特征展平为一维特征输入到lstm后，基于矩阵对三个分量进行融合，最后通过tanh激活函数,实现多尺度电磁频谱占用度图像预测。
[0054]
实施例
[0055]
基于动态阈值和残差卷积网络的电磁频谱占用度二维预测方法是根据时频特征对频谱数据占用进行预测，以芬兰图尔库数据集为例，其频域分辨率为78.125khz，时域分辨率为3s，时间总长为7天。
[0056]
对原始频谱进行预处理。由于原始电磁频谱数据中没有数据缺失，因此无需填补缺失的电磁频谱数据。对于与整体数据严重偏离的异常数据值，可使用其周围3
×
3区域的平均功率谱密度值进行替换。原始电磁频谱数据具有非常精细的分辨率，这使得数据值稀疏、计算量大和挖掘效率低。另一方面，由于原始频谱数据在采集过程中容易与电磁环境噪声混合，因此每隔20个时隙对原始时域分辨率为3s的数据进行平均，以在一定程度上去除设备噪声。
[0057]
对预处理后的电磁频谱数据，采用动态阈值算法提取电磁频谱底噪数据，本方法中电磁频谱底噪的自适应阈值模型表达式如下：
[0058][0059]
其中nk是在对频谱数据进行排序之后，原始频谱数据行中具有k1比率的时隙的数量，m(j)是k1比率时隙的功率谱密度值之和，k2是功率谱密度数据范围的权重，f
max
和f
min
表示最大和最小功率谱密度值。
[0060]
在确定电磁环境的动态噪声阈值之后，将低于阈值的数据拟合为电磁频谱底噪。本方法选择的拟合模型为正交级数耦合模型。电磁环境噪声模型的统一表达式如下：
[0061][0062]
此处，c0、ξj和ζj为常数。
[0063]
对于所述s3将建模后的底噪阈值加上3db作为频谱占用状态阈值，得到频谱占用数据，电磁频谱占用状态阈值表达式如下：
[0064][0065][0066]
此处的表示电磁频谱占用状态阈值，s
i,j
表示第i个时隙，第j个频点的功率值，x
i,j
为频谱占用状态。
[0067]
如图2所示，采用三个多尺度滑动窗口获取时频二维的频谱占用度数据，分别建模电磁频谱占用度的时间近邻性、周期性和趋势性，构建多维张量时频块(tfblock)。
[0068]
构建好电磁频谱预测数据集后，将tfblock的三个分量输入到tfbrl网络中，该网络所具有的se注意机制能有效提升残差卷积网络的预测性能，其中的“挤压”和“激励”核心运算环节可以分别表示为：
[0069][0070]
此处是网络的特征值。
[0071]
δ＝f
ex
(η,w)＝σ(g(η,w))＝σ(w2δ(w1η))
[0072]
这里，通过两个完全连接的层w1和w2来处理在前一步骤中获得的向量η，以获得所需的信道权重值δ。在两个完全连接的层之后，δ中的不同值表示不同通道的权重信息，为通道赋予不同的权重。然后将计算的权重矩阵与输入特征张量相乘，并将权重分配给输入特征张量。
[0073]
然后基于矩阵对三个分量进行融合，最后通过tanh激活函数输出预测二维频谱占用图像。tfbrl网络模型如图3所示。
[0074]
如图4所示，将tfbrl网络模型预测效果与其他常用模型进行指标对比。
[0075]
以上所述，仅是本发明的较佳实施例而已，并非对本发明作任何形式上的限制，虽然本发明已以较佳实施例揭露如上，然而并非用以限定本发明，任何熟悉本专业的技术人员，在不脱离本发明技术方案范围内，当可利用上述揭示的技术内容做出些许更动或修饰为等同变化的等效实施例，但凡是未脱离本发明技术方案内容，依据本发明的技术实质，在本发明的精神和原则之内，对以上实施例所作的任何简单的修改、等同替换与改进等，均仍属于本发明技术方案的保护范围之内。