一种基于时空数据的频谱预测方法

摘要

本发明公开了一种基于时空数据的频谱预测方法，方法包括：获取真实频谱数据矩阵；对所述真实频谱数据矩阵进行预处理，得到待测数据；通过空间性建模和时间性建模，搭建谱图卷积递归神经网络；将所述待测数据输入所述谱图卷积递归神经网络，得到预测结果；根据频谱熵和所述预测结果，确定待接入信道。本发明能够快速分析高精度大规模频谱历史数据,根据频谱信息时间上和空间上的相关性建模，实现对未来频谱的预测，可广泛应用于无线电技术领域。

权利要求

1.一种基于时空数据的频谱预测方法，其特征在于，包括：获取真实频谱数据矩阵；

对所述真实频谱数据矩阵进行截断奇异值分解，得到待压缩数据；

对所述待压缩数据进行压缩处理，得到压缩数据 所述压缩数据 的计算公式为：

其中，Dm×n代表真实频谱数据矩阵；Vnxk0 代表最优右奇异向量；m代表测量数量；n代表频率采样点的数量；k0代表压缩后保留维数；

根据所述真实频谱数据中的测量点构建拓扑图，并将所述压缩数据处理为时间序列，所述拓扑图和所述时间序列为待测数据；

通过空间性建模和时间性建模，搭建谱图卷积递归神经网络；

将所述待测数据输入所述谱图卷积递归神经网络，得到预测结果；

根据频谱熵和所述预测结果，确定待接入信道。

2.根据权利要求1所述的一种基于时空数据的频谱预测方法，其特征在于，所述获取真实频谱数据矩阵，包括：

根据预设的扫描频率和扫描次数，通过频谱分析仪获取真实频谱数据矩阵。

3.根据权利要求1所述的一种基于时空数据的频谱预测方法，其特征在于，所述通过空间性建模和时间性建模，搭建谱图卷积递归神经网络，包括：利用谱图卷积对频谱数据的空间性建模；

以及，

利用门控递归循环单元对频谱数据的时间性建模。

4.根据权利要求3所述的一种基于时空数据的频谱预测方法，其特征在于，所述利用谱图卷积对频谱数据的空间性建模，包括：在图信号和滤波器上进行谱图卷积运算，得到谱图卷积层；

将所述谱图卷积层上的多维特征映射到多维输出；

所述利用门控递归循环单元对频谱数据的时间性建模，包括：通过谱图卷积代替所述门控递归循环单元中的矩阵惩罚，得到谱图卷积门控递归单元；

根据所述谱图卷积门控递归单元，获取时间性建模的输出结果。

5.根据权利要求1所述的一种基于时空数据的频谱预测方法，其特征在于，所述根据频谱熵和所述预测结果，确定待接入信道，包括：将所述预测结果解压到原始数据空间，得到功率谱密度；

获取历史频谱数据；

根据所述历史频谱数据和所述功率谱密度，计算频谱熵；

从未来信道利用率最低的信道中选择频谱熵，确定待接入信道。

6.根据权利要求5所述的一种基于时空数据的频谱预测方法，其特征在于，所述根据所述历史频谱数据和所述功率谱密度，计算频谱熵中，所述频谱熵的计算公式为：其中，hi代表第i个信道的频谱熵；nk代表采样点处于第k个量化子区间的数量；N代采样点的总数；k代表第k个量化子区间。

说明书

一种基于时空数据的频谱预测方法

技术领域

[0001] 本发明涉及无线电技术领域，尤其是一种基于时空数据的频谱预测方法。

背景技术

[0002] 随着移动宽带服务日益增长的需求，无线电频谱资源日益短缺。认知无线电技术是解决频谱资源与频谱利用矛盾的有效途径，其核心思想是认知无线电的闲置资源可以在
不影响授权用户的情况下被认知用户使用，使频谱资源分配从静态变为动态，性能从低效
变为高效。认知无线电需要充分了解频谱特性，更加充分有效地利用频谱资源，其中最为关
键的就是频谱预测技术。

[0003] 先前关于频谱预测的技术，按照方法主要分为基于回归分析，基于马尔可夫和基于机器学习三大类。

[0004] 基回归分析的技术，主要包括移动平均预测模型和自回归(AR)预测模型。但是，这类方法大多是一步预测，而且在回归模型中更新系数的高复杂度问题没有被解决。

[0005] 基于马尔可夫的预测方法有则有着更多的发展。但是该方法是基于信道占用时间和空闲时间服从指数分布的假设，即信道先验知识是完全已知的，但真实环境往往是部分
已知或完全未知的，这限制了隐马尔可夫预测方法的应用范围。

[0006] 基于机器学习的预测方法主要分为支持向量机和神经网络。知识向量机在处理大规模数据集时，学习速度和预测精度明显下降，模型参数不好选择。而之前基于机器学习的
认知无线电网络预测技术的研究往往仅仅考虑在空域或时域上的相关性，即固定地点在一
定时间内预测未来的频谱状态，或者未考虑时间，在固定空间由已知点的频谱数据预测未
知点的频谱数据，而忽略了分析跨越空间和时间域的相关性。

[0007] 从实际测量的频谱数据之间不难发现，数据在时域、频域、空域和能量域等多个维度上密切相关。为了能够从多个维度建立频谱数据模型，分析、挖掘和利用数据之间的潜在
相关性，需要在频率、空间和时间域中获得精确的频谱数据。采集这样一组不同维度的频谱
测量数据会产生巨大的数据量。并且，传统的数据挖掘方法难以处理这么大的数据量。

发明内容

[0008] 有鉴于此，本发明实施例提供一种基于时空数据的频谱预测方法。

[0009] 本发明实施例提供了一种基于时空数据的频谱预测方法，包括：

[0010] 获取真实频谱数据矩阵；

[0011] 对所述真实频谱数据矩阵进行预处理，得到待测数据；

[0012] 通过空间性建模和时间性建模，搭建谱图卷积递归神经网络；

[0013] 将所述待测数据输入所述谱图卷积递归神经网络，得到预测结果；

[0014] 根据频谱熵和所述预测结果，确定待接入信道。

[0015] 优选地，所述获取真实频谱数据矩阵，包括：

[0016] 根据预设的扫描频率和扫描次数，通过频谱分析仪获取真实频谱数据矩阵。

[0017] 优选地，所述对所述真实频谱数据矩阵进行预处理，得到待测数据，包括：

[0018] 对所述真实频谱数据矩阵进行截断奇异值分解，得到待压缩数据；

[0019] 对所述待压缩数据进行压缩处理，得到压缩数据；

[0020] 根据所述真实频谱数据中的测量点构建拓扑图，并将所述压缩数据处理为时间序列。

[0021] 优选地，所述通过空间性建模和时间性建模，搭建谱图卷积递归神经网络，包括：

[0022] 利用谱图卷积对频谱数据的空间性建模；

[0023] 以及，

[0024] 利用门控递归循环单元对频谱数据的时间性建模。

[0025] 优选地，所述利用谱图卷积对频谱数据的空间性建模，包括：

[0026] 在图信号和滤波器上进行谱图卷积运算，得到谱图卷积层；

[0027] 将所述谱图卷积层上的多维特征映射到多维输出；

[0028] 所述利用门控递归循环单元对频谱数据的时间性建模，包括：

[0029] 通过谱图卷积代替所述门控递归循环单元中的矩阵惩罚，得到谱图卷积门控递归单元；

[0030] 根据所述谱图卷积门控递归单元，获取时间性建模的输出结果。

[0031] 优选地，所述根据频谱熵和所述预测结果，确定待接入信道，包括：

[0032] 将所述预测结果解压到原始数据空间，得到功率谱密度；

[0033] 获取历史频谱数据；

[0034] 根据所述历史频谱数据和所述功率谱密度，计算频谱熵；

[0035] 从未来信道利用率最低的信道中选择频谱熵，确定待接入信道。

[0036] 优选地，所述根据所述历史频谱数据和所述功率谱密度，计算频谱熵中，所述频谱熵的计算公式为：

[0037]

[0038] 其中，hi代表第i个信道的频谱熵；nk代表采样点处于第k个量化子区间的数量；N代采样点的总数；k代表第k个量化子区间。

[0039] 本发明的实施例获取真实频谱数据矩阵；对所述真实频谱数据矩阵进行预处理，得到待测数据；通过空间性建模和时间性建模，搭建谱图卷积递归神经网络；将所述待测数
据输入所述谱图卷积递归神经网络，得到预测结果；根据频谱熵和所述预测结果，确定待接
入信道。本发明能够快速分析高精度大规模频谱历史数据,根据频谱信息时间上和空间上
的相关性建模，实现对未来频谱的预测。

附图说明

[0040] 为了更清楚地说明本申请实施例中的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本申请的一些实施例，对于
本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他
的附图。

[0041] 图1为本发明实施例的步骤流程图。

具体实施方式

[0042] 为了使本申请的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本申请进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本申请，并
不用于限定本申请。

[0043] 针对现有技术存在的问题，首先，需要对大面积的频谱测量记录一个非常宽的波段，具有高分辨率，并记录了时间戳信息和地理位置信息，涵盖了时间和空间变化。然后为
了实现快速分析，需要对收集到的频谱数据做预处理。将收集频谱数据按波段每100Mh处理
为真实数据矩阵Dm×n,并用截断奇异值分解作降维处理。真实数据矩阵Dm×n经过奇异值分解
后，通过逐一检验校正，排除掉不包含有用频谱信息的主成分，确定最终保留的维数k0(k0＜
＜m)，得到包含频谱使用的信息右奇异值矩阵 然后通过计算真实数据矩阵Dm×n在每
个主成分上的投影得到压缩数据 然后，在压缩数据上执行预测算法，可以快速计
算。紧接着，将压缩数据按照时序处理为时间序列，并将一个个测量点建模为拓扑图上的
点，边的权重为测量点之间的距离，并引入谱图卷积递归神经网络，这是一种融合频谱信息
中空间和时间相关性的深度学习框架。谱图卷积递归神经网络能够捕捉拓扑图上的空间依
赖性，而时间依赖则使用具有计划采样的编解码器体系结构，实现对频谱的多步时空预测。
通过频谱预测结果判断未来的信道占用情况，结合该点的历史各个信道的频谱熵，可以更
加准确地辅助认知用户选择合适信道。本发明实施例提供了一种基于时空数据的频谱预测
方法，如图1所示，包括：

[0044] 获取真实频谱数据矩阵；

[0045] 对所述真实频谱数据矩阵进行预处理，得到待测数据；

[0046] 通过空间性建模和时间性建模，搭建谱图卷积递归神经网络；

[0047] 将所述待测数据输入所述谱图卷积递归神经网络，得到预测结果；

[0048] 根据频谱熵和所述预测结果，确定待接入信道。

[0049] 优选地，所述获取真实频谱数据矩阵，包括：

[0050] 根据预设的扫描频率和扫描次数，通过频谱分析仪获取真实频谱数据矩阵。

[0051] 优选地，所述对所述真实频谱数据矩阵进行预处理，得到待测数据，包括：

[0052] 对所述真实频谱数据矩阵进行截断奇异值分解，得到待压缩数据；

[0053] 对所述待压缩数据进行压缩处理，得到压缩数据；

[0054] 根据所述真实频谱数据中的测量点构建拓扑图，并将所述压缩数据处理为时间序列。

[0055] 优选地，所述通过空间性建模和时间性建模，搭建谱图卷积递归神经网络，包括：

[0056] 利用谱图卷积对频谱数据的空间性建模；

[0057] 以及，

[0058] 利用门控递归循环单元对频谱数据的时间性建模。

[0059] 优选地，所述利用谱图卷积对频谱数据的空间性建模，包括：

[0060] 在图信号和滤波器上进行谱图卷积运算，得到谱图卷积层；

[0061] 将所述谱图卷积层上的多维特征映射到多维输出；

[0062] 所述利用门控递归循环单元对频谱数据的时间性建模，包括：

[0063] 通过谱图卷积代替所述门控递归循环单元中的矩阵惩罚，得到谱图卷积门控递归单元；

[0064] 根据所述谱图卷积门控递归单元，获取时间性建模的输出结果。

[0065] 优选地，所述根据频谱熵和所述预测结果，确定待接入信道，包括：

[0066] 将所述预测结果解压到原始数据空间，得到功率谱密度；

[0067] 获取历史频谱数据；

[0068] 根据所述历史频谱数据和所述功率谱密度，计算频谱熵；

[0069] 从未来信道利用率最低的信道中选择频谱熵，确定待接入信道。

[0070] 优选地，所述根据所述历史频谱数据和所述功率谱密度，计算频谱熵中，所述频谱熵的计算公式为：

[0071]

[0072] 其中，hi代表第i个信道的频谱熵；nk代表采样点处于第k个量化子区间的数量；N代表采样点的总数；k代表第k个量化子区间。

[0073] 下面结合说明书附图，对本发明的具体实现过程进行详细描述：

[0074] 设在一个几十平方区域部署N台带有含全球定位系统(北斗或者GPS)模块的频谱分析仪作为数据的收集设备。选用ESMD频谱分析仪每100M赫兹分辨率为26215，每5分钟扫
描一次300Mhz‑4Ghz频段，共分为37个100M赫兹频带。数据收集时长为10个月，所有时长每
100M赫兹共扫描次数为60k次，将扫描频段每100M赫兹整合为真实频谱数据矩阵Dm×n。

[0075] 第一步，对Dm×n作截断奇异值分解并压缩数据。

[0076] 本步骤为对真实频谱数据的功率谱密度部分作数据预处理，但需要完整保留每次扫描的时间戳和地理位置信息。

[0077] 将每个100MHz波段的PSD数据，处理为具有m行和n列的矩阵Dm×n，其中m是测量数量，n是100Mhz中频率采样点的数量。

[0078] 将Dm×n作保留k维的截断奇异值分解公式如下：

[0079] Dm×n≈Um×k·Sk×k·(Vn×k)T

[0080] 其中，Vn×k＝[v1,v2,v3,...,vk]最优右奇异向量，包含关于频谱使用的有用信息。这k列向量称为Dm×n前k个主成分。

[0081] 其中，通过逐一校验校正，得到最佳保留维数。通过对数数据的分析，发现数据的绝大部分有用信息在1000维以下，从1000维从大到小开始，逐一计算每一主成分与高斯分
布的夏皮罗‑威尔克校验可信度，直到找到一个k0使得其为高斯分布的可信度低于阈值
0.99，确定保留的维数为700。

[0082] 在频域奇异值分解后，取每100Mhz的最大最优维数700。对所有每个100MHz频段作截断奇异值分解并得到主成分矩阵 然后通过计算在Dm×n每个主成分上的投影来压缩
数据集并减小数据维数，得到压缩数据 至此，将所有的数据降为k0维，公式如下：

[0083]

[0084] 第二步，输入谱图卷积递归神经网络前的数据预处理。

[0085] 本步骤根据保留的地理位置信息和时间信息，分别将测量点建成拓扑图和将压缩数据处理为时间序列。

[0086] 将测量点组成的网络表示为加权拓扑图G＝(v,ε,w)，其中v是测量节点集合，|v|N×N
＝N,ε是一组边，W∈R 是表示节点邻近性的加权邻接矩阵，边的权重由数据中的经纬度信
息根据半正矢公式测算距离，如下：

[0087]

[0088] 其中， 和 表示两点的纬度，λ1和λ2表示两点的经度。

[0089] 在G上观察到频谱数据表示为图信号 其中，N为测量点的数量,k0上述为确定保留的维数。

[0090] 设X(t)表示在时间t处观察到的图形信号，频谱预测问题旨在学习一个h(·)函数，给定一个图G，函数h(·)将历史拓扑图信号映射到未来的拓扑图信号，

[0091]

[0092] 第三步，谱图卷积递归神经网络搭建及训练。

[0093] 本步骤主要为利用谱图卷积对频谱数据的空间性建模，和利用门控递归循环单元对频谱数据的时间性建模。以扩散卷积门控循环单元搭建序列到序列体系结构，实现多步
预测。

[0094] 谱图卷积学习图结构数据的表示，可以使用基于随机梯度的方法来训练，其在图信号 和滤波器fθ上得到的谱图卷积运算定义如下，

[0095]

[0096] 其中，谱图卷积利用归一化图拉普拉斯的概念有 参数fθ为Λ的K阶多项式，并利用稳定的切比雪夫多项式基计算。归一化处理，用
表示将特征值从[0，λmax]映射到[‑1,1]，其中λmax表示为特征值矩阵L最大的特征值。
为关于 的K阶多项式， 为归一化后的待训练参数，X:,p为图信号.

[0097] 利用上式子中定义的卷积运算，我们可以建立一个谱图卷积层，将k0维特征映射到Q维输出,如下式，

[0098]N×Q

[0099] 其中， 为输入，H∈R 为输出，fΘ为滤波器，a为ReLU激活函数。

[0100] 然后，使用递归神经网络(RNNs)来捕捉时间依赖性。这里使用门控递归单元(GRU)，这是RNNs的一个简单而强大的变体。用谱图卷积代替GRU中的矩阵乘法形成谱图卷
积门控递归单元(GCGRU),其公式如下，

[0101]

[0102]

[0103] 其中，其中X(t)、H(t)分别表示时间t的输入和输出，r(t)、u(t)分别是时间t的复位门和更新门。表示上述定义的谱图卷积，Θu，ΘC是相应滤波器的参数。用GCGRU来建立递归
神经网络层，并通过时间反向传播进行训练。为了实现多步预测中，基于GCGRU搭建序列到
序列体系结构，编码器和解码器都包含两个递归层，每个循环层有64个单位。

[0104] 在训练过程中，将历史时间一个小时序列输入编码器，并使用其最终状态初始化‑2
解码器，得到未来一个小时的频谱预测，反复训练。初始学习速率为1e ，从第20个epoch开
始，每10个epoch减少 整个网络是通过最大限度地利用时间反向传播生成目标未来时
间序列的来训练。在测试时，输入的标签被编码器本身产生的预测所取代。在收集到的数据
集中，将频谱数据的70％用于培训，20％用于测试，其余10％用于验证。训练和测试的输入
分布之间的差异会导致性能下降。为了缓解这一问题，将计划抽样集成到模型中，通过将输
入解码器的真实标签概率逐渐减少到0，以允许模型学习测试分布。

[0105] 第四步，利用谱图卷积递归神经网络得到预测结果，并引入频谱熵辅助认知用户选择信道。

[0106] 本步骤通过网络预测得到的未来一个小时预测结果向量，将估计的投影解压回原始数据空间得到未来一个小时测量点处高精度的功率谱密度，并且结合过去24小时历史频
谱数据，通过计算每小时的信道利用率，进一步计算得到衡量各个信道稳定的频谱熵(在不
同的频段，对信道的划分不同)，结合未来的信道状态，辅助认知用户更好选择接入信道。

[0107] 通过由历史和预测得到的未来功率谱密度与功率阈值(‑100dbm)进行比较，计算信道的利用率。并且结合给定信道历史利用率{ui(1),ui(2),…,ui(N)}，可以根据以下公式
计算频谱熵信道衡量信道的稳定性：

[0108]

[0109] 最后在未来信道利用率最低的信道中选择频谱熵最低的，选择最稳定的信道。

[0110] 综上所述，本发明首先对原始频谱数据采用了特别设计的截断式奇异值分解，能够有效地保留原始数据特征的同时，降低数据的维数，使得后续在压缩空间进行运算，使得
对数据的快速分析成为可能。其次，本发明提出了对时空数据类型的频谱信息的建模方式，
使用调试的深度学习框架，有效提取频谱数据时间上和空间上的相关性，区别于以往的频
谱预测仅仅依赖时间或空间上的相关性进行预测，并且预测的精度更高，结合频谱熵有利
于认知用户选择更合适的信道接入。第三，本发明频谱数据对应的地域较广，对频谱数据的
预测和收集有利于当地的无线电管理局查明哪些频带在大面积长期无活跃信号，并可以考
虑为二次用户开放频带，提高频谱利用率。

[0111] 在本说明书的描述中，参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特
点包含于本发明的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不
一定指的是相同的实施例或示例。而且，描述的具体特征、结构、材料或者特点可以在任何
的一个或多个实施例或示例中以合适的方式结合。

[0112] 尽管已经示出和描述了本发明的实施例，本领域的普通技术人员可以理解：在不脱离本发明的原理和宗旨的情况下可以对这些实施例进行多种变化、修改、替换和变型，本
发明的范围由权利要求及其等同物限定。

[0113] 以上是对本发明的较佳实施进行了具体说明，但本发明并不限于所述实施例，熟悉本领域的技术人员在不违背本发明精神的前提下还可做作出种种的等同变形或替换，这
些等同的变形或替换均包含在本申请权利要求所限定的范围内。