CN102186175B - 基于振荡器相位同步的认知网络动态频谱分配方法 - Google Patents

Abstract

本发明提出了一种基于振荡器相位同步的认知网络动态频谱分配方法，克服现有技术中频谱使用效率较低且用户花费较高的问题。其实现步骤是：(1)绘制频谱干扰图；(2)取补图；(3)初始化；(4)求解相位值；(5)判定相位同步；(6)节点分组；(7)同步组排序；(8)分配频谱；(9)获得频谱分配方案。本发明提出的方法减少了频谱分配过程中所需的时间，有效提高了频谱的使用效率；同时，该方法在频谱分配的过程中充分考虑了用户购买频谱时的花费，给尽可能多的用户分配费用较低的频谱，因而降低了认知网络中用户的总花费。

Description

基于振荡器相位同步的认知网络动态频谱分配方法

技术领域

本发明属于通信技术领域，更进一步涉及无线通信网络中对认知网络进行动态频谱分配的方法，利用振荡器相位同步原理，提高了并行处理能力，可快速有效地得到认知网络的频谱分配方案。

背景技术

认知无线电技术被认为是解决无线电通信中频谱资源使用效率低下的最佳方案，该技术从频谱再利用的思想出发，在不影响授权用户正常通信的条件下，将授权用户所空闲的频谱分配给认知网络中的用户使用，以提高空闲频谱的使用效率。对认知网络进行动态的频谱分配是认知网络中的关键技术之一，人们也提出了很多方法来对认知网络进行动态的频谱分配，主要包括基于博弈理论的方法，基于图论的方法，和基于拍卖理论的方法。其中基于图论的方法是将认知网络的频谱分配问题转换成图着色模型，然后利用一些解决图着色问题的算法，如遗传算法，启发式算法等，来对该模型进行图着色，从而获得认知网络的频谱分配结果。现有的基于图论的频谱方法，大都没有考虑进行分配时所用的时间和认知用户在购买频谱时的花费问题，因此现有的方法有可能造成频谱资源在时间上的浪费和增加了用户通讯时的花费。

例如，东南大学申请的专利“认知无线电系统中一种基于图论的频谱分配方法”(专利申请号200910034227.7，公开号CN101662783A)。该专利申请所公开的基于图论的频谱分配方法是在计算分配优先级时，将信道可用度和信道连续度综合考虑，同时引入区域呼损率参数，减低了呼损率。该专利申请的方法存在不足之处是，不能实时地得到频谱分配结果，需用相对较长的时间来完成对整个网络进行频谱分配的过程，从而导致了频谱资源在时间上的浪费；同时，该方法在进行频谱分配的过程中侧重考虑了减低呼损率而忽视了认知用户的通讯花费，导致认知用户购买频谱时花费增多，影响了认知网络中频谱资源的优化性。

发明内容

本发明的目的在于克服上述已有技术的不足，提出一种基于振荡器相位同步的认知网络动态频谱分配方法，以有效地提高频谱的使用效率，同时降低认知网络中用户购买频谱时的总花费。

本发明的具体步骤如下：

(1)绘制频谱干扰图：将N个认知用户组成的认知网络拓扑结构抽象成图，图中每一个节点与每一个认知用户一一对应。

(2)取补图：对频谱干扰图取补图，写出该补图的邻接矩阵，若补图节点i与节点j之间有连接边，则其邻接矩阵中的对应元素a_ij＝1，否则a_ij＝0，矩阵对角线上的元素a_ii＝1。

(3)初始化：生成与认知网络中可使用的频谱数M相同个数的空同步组。

(4)求解相位值：补图中的每一个节点与一个振荡器一一对应，振荡器的相位变化值对应节点的相位变化值，在MATLAB环境下，通过求解以下方程获得每一个振荡器在0～t时刻内的相位变化值：

${\stackrel{\·}{\mathrm{\θ}}}_i=w_i+\frac{K}{K_{\max }}\underset{j=1}{\overset{N}{\mathrm{\Σ}}}{a}_{\mathrm{ij}}\sin ({\mathrm{\θ}}_j-{\mathrm{\θ}}_i),(i=\mathrm{1,2}\·\·\·N)$

其中，

为振荡器i的相位随时间的变化率，w_i为振荡器i的固有频率，该固有频率在[-0.1，0.1]之间随机产生且服从均匀分布，K为耦合强度，K_max为补图中所有节点度的最大值，a_ij为补图的邻接矩阵中的对应元素，θ_i为振荡器i的相位，振荡器的初始相位在[0，2π]之间随机产生且服从均匀分布。

(5)判定相位同步：将N个振荡器在t时刻的相位值θ_1t，θ_2t…θ_Nt两两进行比较，若相位θ_it与相位θ_jt满足关系式T＜cos(θ_it-θ_jt)，T(0＜T＜1)为阈值常数，则判定振荡器i与振荡器j相位同步，否则，判定振荡器i与振荡器j相位不同步。

(6)节点分组：在步骤(5)中得到的相位同步结果中，将每一组相位同步的振荡器所对应的节点放入一个同步组中，以此类推，完成将N个节点划分到M个同步组中的过程。

(7)同步组排序：将M个同步组按照其节点数从大到小的顺序依次排序。

(8)分配频谱：将M个频谱按照其费用从小到大的顺序依次分配给步骤(7)中排序后的同步组，每一个同步组都对应一个频谱。

(9)获得频谱分配方案：将同步组所对应的频谱分配给该同步组中所有节点所代表的用户共同使用，以此类推，完成给认知网络中的所有用户分配频谱的过程。

本发明与现有技术相比存在以下优点：

第一，由于本发明利用了振荡器的相位同步原理，可并行处理认知网络的频谱分配问题，克服了现有技术中处理较慢，耗费时间较长的问题。本发明的求解振荡器相位变化过程是一个并行处理过程，可大大减少对整个网络进行频谱分配时所需的时间，能够实时地对认知网络进行频谱分配，从而有效地提高频谱的使用效率。

第二，由于本发明在对认知网络进行频谱分配的过程中考虑了用户购买频谱时的花费，克服了现有技术中忽视认知用户通讯花费，导致认知用户购买频谱时花费增多的问题。本发明给尽可能多的用户分配费用较低的频谱，从而降低了认知网络中用户的总花费。

附图说明

图1为本发明的流程图；

图2为本发明实施例中认知网络的频谱干扰图；

图3为本发明实施例中认知网络的频谱干扰图的补图；

图4为本发明实施例中振荡器的相位变化图；

图5为本发明实施例中对认知网络用户的频谱分配结果图；

图6为本发明仿真实验中所采用的三组网络的频谱干扰图；

图7为本发明仿真实验中振荡器的相位变化图。

具体实施方式

参照图1，本发明的具体实施步骤如下：

步骤1.绘制频谱干扰图

本发明的实施例选择一个具体的认知网络进行频谱分配，绘制出的该认知网络的频谱干扰图如图2所示，图2中的12个节点1～12分别代表认知网络中12个不同的用户，有边连接的两个节点所代表的用户在通讯时会产生干扰，不能分配同一个频谱。

步骤2.取补图

对图2中的认知网络频谱干扰图取补图的步骤是，保持认知网络频谱干扰图中节点的网络拓扑结构不变，将每对节点间原频谱干扰图中有连接边的去掉，没有连接边的加上，就得到了该认知网络频谱干扰图的补图。图3中每一条边所连接的两个节点表示通讯时没有干扰的用户，即它们可共用频谱。若补图中节点i与节点j之间有连接边，则其邻接矩阵中的对应元素a_ij＝1，否则a_ij＝0，且a_ii＝1，那么可写出图3的邻接矩阵为：

$\left\{\begin{array}{ccccc}1& 0& \·\·\·& 1& 1\\ 0& 1& \·\·\·& 1& 1\\ \·& \·& & \·& \·\\ \·& \·& & \·& \·\\ \·& \·& & \·& \·\\ 1& 1& \·\·\·& 1& 0\\ 1& 1& \·\·\·& 0& 1\end{array}\right\}$

步骤3.初始化

生成与认知网络中可使用的频谱数相同个数的空同步组。本发明实施例中，该认知网络中有3个可使用的频谱，因此初始3个空同步组C₁，C₂，C₃。用户购买频谱P₁，P₂，P₃的花费分别为b₁，b₂，b₃，满足b₁＜b₂＜b₃。

步骤4.求解相位值

参照图3，补图中的每一个节点与一个振荡器一一对应，振荡器的相位变化值对应节点的相位变化值，此12个振荡器分别记为1～12。在MATLAB环境下，通过求解以下方程获得每一个振荡器在0～3秒时刻内的相位变化值，取值时间间隔为0.01秒：

${\stackrel{\·}{\mathrm{\θ}}}_i=w_i+\frac{K}{K_{\max }}\underset{j=1}{\overset{N}{\mathrm{\Σ}}}{a}_{\mathrm{ij}}\sin ({\mathrm{\θ}}_j-{\mathrm{\θ}}_i),(i=\mathrm{1,2}\·\·\·12)$

其中，

为振荡器i的相位随时间的变化率；w_i为振荡器i的固有频率，该固有频率在[-0.1，0.1]之间随机产生且服从均匀分布；耦合强度K＝6；补图中节点度指一个节点所拥有的连接边的个数，那么补图中节点度的最大值K_max＝4；a_ij为图3邻接矩阵中的对应元素；θ_i为振荡器i的相位，初始相位在[0，2π]之间随机产生且服从均匀分布。

通过求解振荡器的相位微分方程，可得到该12个振荡器的相位随时间的变化图，如图4所示。

步骤5.判定相位同步

将12个振荡器在第3秒时的相位值θ_1t，θ_2t…θ_Nt两两进行比较，若相位θ_it与相位θ_jt满足关系式T＜cos(θ_it-θ_jt)，T＝0.99为阈值常数，则判定振荡器i与振荡器j相位同步；否则，判定振荡器i与振荡器j相位不同步。判定后可知振荡器1，2，3，4，5之间两两相位同步，振荡器6，7，8，9之间两两相位同步，振荡器10，11，12之间两两相位同步。

步骤6.节点分组

在步骤5中得到的振荡器相位同步结果中，由于振荡器1，2，3，4，5之间两两相位同步，则将其对应的节点1，2，3，4，5放入第一个同步组C₁，即C₁{1，2，3，4，5}。以此类推，完成将12个节点划分到3个同步组中的过程，可得第二个同步组C₂{6，7，8，9}和第三个同步组C₃{10，11，12}。

步骤7.同步组排序

将同步组C₁，C₂，C₃按照其节点数从大到小的顺序依次排列。由于同步组C₁中有5个节点，同步组C₂中有4个节点，同步组C₃中有3个节点，因而排列好的同步组顺序为C₁，C₂，C₃。

步骤8.分配频谱

将3个频谱按照其费用从小到大的顺序依次分配给步骤7中排序后的同步组。因为用户购买频谱P₁，P₂，P₃时的花费分别为b₁，b₂，b₃，满足b₁＜b₂＜b₃，那么将频谱P₁，P₂，P₃依次分配给同步组C₁，C₂，C₃。

步骤9.获得频谱分配方案

将同步组所对应的频谱分配给该同步组中所有节点所代表的用户共同使用，即节点1，2，3，4，5所代表的用户共用频谱P₁，节点6，7，8，9所代表的用户共用频谱P₂，节点10，11，12所代表的用户共用频谱P₃，就得到了本实施例中认知网络的频谱分配结果，如图5所示。该频谱分配结果中，给尽可能多的用户分配费用较低的频谱，从而降低了认知网络中用户的总花费。

本发明的效果可以通过以下实验进一步说明：

1.仿真条件：

在CPU为core2 2.4GHZ、内存2G、WINDOWS XP系统上使用MATLAB进行了仿真。

2.仿真内容：

选取三组不同结构的网络作为实验对象，此三组网络的频谱干扰图如图6所示，其中图6(a)、图6(b)和图6(c)分别为网络1、网络2和网络3的频谱干扰图。分别用基于遗传算法的图着色方法和本发明中提出的方法求出这三组网络的频谱分配结果，并记录其运行的时间。

基于遗传算法的图着色方法是将认知网络的频谱分配问题转换成传统的图着色模型，且每一个着色方案对应一个频谱分配结果，并用遗传算法的原理来求出该图着色模型的着色方案。该方法基于遗传算法原理，即在整个解空间里，通过对初始种群不断进行交叉操作和变异操作来不断产生下一代种群，进而逐渐逼近着色结果。而本发明中提出的方法主要利用振荡器的相位同步原理来对认知网络进行动态频谱分配。

实验中，分别用以上两种方法对这三组认知网络进行频谱分配，求出其频谱分配结果，并记录其运行时间。在采用本发明中提出的方法进行频谱分配时，能够得到这三组网络所对应的振荡器的相位变化图，如图7所示，其中图7(a)、图7(b)和图7(c)分别为网络1、网络2和网络3在实验中所对应的相位变化图，根据该相位变化图可以得到这三个网络的频谱分配结果。同时，下表为实验中记录的数据结果。

表1两种方法下求解频谱分配结果所需的时间

由表1可知，采用基于遗传算法的图着色方法和本发明中提出的方法，都可以求出满足认知网络中所有用户正常通讯时所需要的最少频谱数。但是，在求同一个认知网络的频谱分配结果时，本发明中提出的方法所需的时间远远小于基于遗传算法的图着色方法。

从以上说明可以看出，基于振荡器相位同步的认知网络动态频谱分配方法可大大减少对整个网络进行频谱分配过程中所需的时间，能够实时地对认知网络进行频谱分配，从而有效地提高频谱的使用效率。

上述实施方式仅是本发明的一个实例，不构成对本发明的任何限制，例如用本发明的方法还可以对包含20个节点的认知网络，对包含50个节点认知网络，以及对包含更多节点的认知网络进行频谱分配。

Claims (2)

1.一种基于振荡器相位同步的认知网络动态频谱分配方法，具体步骤如下：

(1)绘制频谱干扰图：将N个认知用户组成的认知网络拓扑结构抽象成图，图中每一个节点与每一个认知用户一一对应；

(2)取补图：对频谱干扰图取补图，写出该补图的邻接矩阵，若补图节点i与节点j之间有连接边，则其邻接矩阵中的对应元素a_ij＝1，否则a_ij＝0，矩阵对角线上的元素a_ii＝1；

(3)初始化：生成与认知网络中可使用的频谱数M相同个数的空同步组；

(4)求解相位值：补图中的每一个节点与一个振荡器一一对应，振荡器的相位变化值对应节点的相位变化值，在MATLAB环境下，通过求解以下方程获得每一个振荡器在0～t时刻内的相位变化值：

${\stackrel{\·}{\mathrm{\θ}}}_i=w_i+\frac{K}{K_{\max }}\underset{j=1}{\overset{N}{\mathrm{\Σ}}}{a}_{\mathrm{ij}}\sin ({\mathrm{\θ}}_j-{\mathrm{\θ}}_i),(i=\mathrm{1,2}\·\·\·N)$

其中，

为振荡器i的相位随时间的变化率，w_i为振荡器i的固有频率，该固有频率在[-0.1，0.1]之间随机产生且服从均匀分布，K为耦合强度，K_max为补图中所有节点度的最大值，a_ij为补图的邻接矩阵中的对应元素，θ_i为振荡器i的相位，振荡器的初始相位在[0，2π]之间随机产生且服从均匀分布；

(5)判定相位同步：将N个振荡器在t时刻的相位值θ_1t，θ_2t…θ_Nt两两进行比较，若相位θ_it与相位θ_jt满足关系式T＜cos(θ_it-θ_jt)，T(0＜T＜1)为阈值常数，则判定振荡器i与振荡器j相位同步，否则，判定振荡器i与振荡器j相位不同步；

(6)节点分组：在步骤(5)中得到的相位同步结果中，将每一组相位同步的振荡器所对应的节点放入一个同步组中，以此类推，完成将N个节点划分到M个同步组中的过程；

(7)同步组排序：将M个同步组按照其节点数从大到小的顺序依次排序；

(8)分配频谱：将M个频谱按照其费用从小到大的顺序依次分配给步骤(7)中排序后的同步组，每一个同步组都对应一个频谱；

(9)获得频谱分配方案：将同步组所对应的频谱分配给该同步组中所有节点所代表的用户共同使用，以此类推，完成给认知网络中的所有用户分配频谱的过程。

2.根据权利要求1所述的基于振荡器相位同步的认知网络动态频谱分配方法，其特征在于，所述步骤(2)中对频谱干扰图取补图的步骤是，保持认知网络频谱干扰图中节点的网络拓扑结构不变，将每对节点间原频谱干扰图中有连接边的去掉这个连接边，没有连接边的加上连接边，就得到了该认知网络频谱干扰图的补图。

Images