CN103957528A - 一种基于认知无线电的多策略动态频谱接入方法 - Google Patents

Abstract

本发明提供了一种基于认知无线电的多策略动态频谱接入方法，应用于认知无线电系统，系统中的主次用户有着不同的使用优先级，当一个次级用户到达，SAP首先感知cP信道，如果有空闲信道便随即挑选一个接入，否则SAP感知RP信道，再根据结果做出决定；主用户接入过程与次级用户相似但接入信道顺序相反，先接入预留信道，若预留信道没有空闲，则尝试RP信道，若RP信道也没有空闲，最后尝试cP信道。通过使用本发明的多策略动态频谱接入方法，主系统通过SLAB暂时将共享信道中的频谱使用权租赁给次级用户，并从次级系统获取回报，次系统根据本发明提出的MS—DSA方法为新到达的用户分配频谱，一旦获得SLAB的许可，次级系统便暂时获得接入许可。

Description

一种基于认知无线电的多策略动态频谱接入方法

技术领域

本发明涉及无线通信领域，尤其涉及一种多策略动态频谱接入方法。

背景技术

传统的无线通信系统采用的是静态(固定)频谱分配政策，即授权用户对其频谱具有独占性，其他用户不允许使用。虽然这种频谱分配方式有利于保证系统的服务质量，但由于通信业务在地域、时域、频域上的不均衡性，频谱资源并没有得到充分利用。以美国为例,美国联邦通信委员会(FCC)的大量研究报告说明频谱的利用情况极不平衡,一些非授权频段占用拥挤,而有些授权频段则经常空闲。而随着无线通信业务需求的快速增长，无线频谱资源显得日益缺乏。可用频谱资源稀缺及频谱利用效率低等问题使得一种全新的、优化使用频谱资源的无线通信模式变得十分必要。因此，基于认知无线电(Cognitive Radio,CR)的动态频谱接入(DynamicSpectrum Access，DSA)技术应运而生。根据IEEE1900.1草案标准的定义，DSA是一种可以在一定频谱使用权限范围内动态地选择工作频谱，并利用在一定时域、空域和频域上出现的频谱空洞或白区进行通信的无线电技术。

DSA是一个广义的概念，它包含了多种频谱改革方法。DSA策略可概括为三种模型：动态排他使用模型(Dynamic Exclusive Use Model)、开放共享模型(Open Sharing Model)和分层接入模型(Hierarchical AccessModel)。

1.动态排他使用模型

动态排他使用模型仍沿用了目前静态频谱分配政策的基本结构，其主要思想是引入灵活性以提高频谱利用效率。目前该模型主要包含两种技术：频谱产权(Spectrum Property Rights)和动态频谱分配(Dynamic SpectrumAllocation)。

频谱产权可解释为在指定频段上传输信号的权利，前提条件是在指定的时间段内以及指定的地域范围内使用频谱，并且信号功率不超过指定强度。频谱产权策略允许某些频段的授权用户自由地交易或者出租其频段。然而对于这种频谱改革技术，清楚地定义频谱的权利并且将其作为实物予以执行具有很大的技术难度。

动态频谱分配这种DSA策略利用通信业务在时域和空域上的统计特性，通过对共存的无线业务实现时间相关和空间相关的频谱共享，以达到提高频谱利用效率的目的。动态频谱分配策略的主要缺点在于它并不能完全消除和充分利用由于通信业务的突发性而产生的频谱空洞。

2.开放共享模型

这个模型向所有用户开放频谱使其共享，如工业、科学和医用(ISM)频段的开放共享方式。基于基础设施的集中式频谱共享和分布式频谱共享是应用于该模型的两类主要频谱共享技术。

3.分层接入模型

可看作是动态排他使用模型和开放共享模型二者的结合。其基本思想是在次要用户(Secondary User)不对授权用户产生任何有害干扰的前提下，向次要用户开放已授权的频段，从而实现授权用户和次要用户对频谱的共享。用来实现此共享的技术包括频谱衬垫(Spectrum Underlay)和频谱填充(Spectrum Overlay)。

与动态排他使用模型和开放共享模型相比，分层接入模型可能是与现行的频谱管理政策和传统的无线系统最兼容。而且，频谱衬垫和频谱填充方法可同时使用来进一步提高频谱效率。本发明所涉及的动态频谱接入方法就是基于该模型的一种改进。

发明内容

为了解决现有技术中问题，本发明提供了一种基于认知无线电的多策略动态频谱接入方法，能够实现效用最佳的动态频谱分配方案。

本发明通过如下技术方案实现：

一种基于认知无线电的多策略动态频谱接入(MS-DSA)方法，应用于认知无线电系统，采用随机抢先DSA(RP-DSA)和协作DSA(CP-DSA)策略，该方法包括以下步骤：

S1：通过频谱租赁与拍卖SLAB完成交易，获得频谱，初始化信道总数

N，预留信道个数γ_Res、CP信道个数γ_CP和RP信道个数γ_RP；

S2：判断新到达的用户的类别，若为主用户PU，则转到步骤S3，若为次级用户SU，则转到步骤S10；

S3：若PU占用预留信道的个数小于等于γ_Res，则转到S4，否则，转到步骤S5；

S4：PU随机接入一个空闲的预留信道；

S5：若PU占用RP信道的个数与PU和SU同时占用RP信道的个数之和小于等于γ_RP，则转到S6，否则，转到步骤S7；

S6：PU按照RP-DSA策略接入一个RP信道；

S7：若PU占用CP信道的个数与PU和SU同时占用CP信道的个数之和小于等于CP信道的总数，则转到步骤8，否则，转到步骤S9；

S8：按照RP-DSA策略接入一个RP信道；

S9：主用户PU被阻塞；

S10：若SU占用CP信道的个数与PU和SU同时占用CP信道的个数之和小于等于γ_CP，则转到S11；

S11：次级用户访问点SAP感知CP信道，若感知到有空闲的CP信道，则SU按照CP策略接入一个CP信道，否则，转到S12；

S12：若SU占用RP信道的个数与PU和SU同时占用RP信道的个数之和小于等于γ_RP，则转到S13；

S13：SAP感知RP信道，若感知到有空闲的RP信道，则SU按照RP策略接入一个RP信道；否则，转到S14；

S14：次级用户SU被阻塞；

S15：SAP感知被SU所占用的信道，若检测到PU在该信道上传输，则SU从正在传输的信道中离开。

本发明的另一方面，提出了一种认知无线电系统，所述系统包括主系统、次级系统和频谱租赁和拍卖SLAB；所述主系统将下行授权频谱资源划分为两个部分，一部分是仅为主用户传输留下的预留信道，另一部分是与次级用户共享的使用动态频谱分配技术的信道，主系统通过SLAB暂时将共享信道中的频谱使用权租赁给次级用户，并从次级系统获取回报；所述次级系统包括：DSA策略数据库、频谱获取模块、频谱感知模块、接入控制模块和性能评估模块；在一个协调器的管理下次级系统建立一个次级用户访问点SAP，SAP发送给所述SLAB一个包括频谱需求、支持的动态频谱接入DSA策略、频谱租用价格和对主用户干扰的惩罚的请求，所述认知无线电系统采用如权利要求1所述的多策略动态频谱接入方法为新到达的用户分配频谱，一旦获得SLAB的许可，次级系统便暂时获得接入许可。

本发明的有益效果是：本发明的基于认知无线电的多策略动态频谱接入方法，相比于现有技术，打破了认知网络只能运行单一的动态频谱接入策略的限制，通过与单一策略的仿真对比，以复杂度增加的代价获取了吞吐量及效用的提高。

附图说明

图1是本发明的多策略动态频谱接入系统的次级系统框架图；

图2是本发明的多策略动态频谱接入方法的信道分配示意图；

图3是本发明的多策略动态频谱接入方法流程图。

具体实施方式

下面结合附图说明及具体实施方式对本发明进一步说明。

认知无线电网络的基础结构由三个实体组成：主系统，次级系统和SLAB(频谱租赁和拍卖)。

对于主系统，它将下行授权频谱资源划分为两个部分，一部分是仅为主用户传输留下的预留信道，另一部分是可以与次级用户共享的使用动态频谱分配技术的信道。主系统通过SLAB暂时将共享信道中的频谱使用权租赁给次级用户，并从次级系统获取回报。SLAB介于主系统和次级系统之间具有管理和授权作用。此外为了保证安全，假设主系统在任何时刻都不给SLAB和次级系统发送信息。

对于次级系统，次级用户在一个协调器的管理下建立一个SAP(次级用户访问点)。首先，SAP发送给SLAB一个包括频谱需求、支持的DSA(动态频谱接入)策略、频谱租用价格和对主用户干扰的惩罚的请求。一旦获得SLAB的许可，次级系统便暂时获得接入许可。此外，假设有一个控制信道来控制SAP和次用户之间的联络，以便SAP对所有次用户的传输进行协调和避免次用户之间的冲突。此外，还假设感知和接入时间要远小于服务时间的持续时间。

如附图1所示，本发明提出的多策略动态频谱接入(MS-DSA)系统中的次级系统由5个基本要素组成：DSA策略数据库、频谱获取、频谱感知、接入控制和性能评估。

首先，DSA策略可以总结为三类：RP-DSA(随机抢先DSA)、CP-DSA(协作DSA)、NP-DSA(非抢占式DSA)。本发明所提出的方法主要用到RP-DSA和CP-DSA，并在主次系统共享信道中使用它们。

在RP-DSA策略中，主系统拥有绝对的优先级，次级系统对于主系统而言是完全透明的——即当一个主用户到达主系统将该主用户随机分配到任意信道而忽视次级用户的活动。在CP-DSA策略中，一个新到的主用户将不会被分配到正被次级用户所使用的信道，从而避免对次级用户传输的干扰。因此，只有在信道全忙的情况下次级用户才会产生掉话。为方便起见，在这里将使用RP-DSA和CP-DSA策略的信道分别称为“RP信道”和“CP信道”，而主用户所单独使用的信道称之为“预留信道”。在执行MS-DSA协议之前，假设光谱的定价已经通过SLAB完成。ω_RP和ω_CP表示每个RP和CP信道的租用价格，γ_RP、γ_CP和γ_Res分别表示RP、CP和预留信道的个数，N表示信道总数，γ_RP+γ_CP+γ_Res＝N。

在不同的DSA策略中，主次用户有着不同的使用优先级。当一个次级用户到达，SAP首先感知CP信道，如果有空闲信道便随即挑选一个接入，否则SAP感知RP信道，再根据结果做出决定。如附图2所示，主用户PU接入过程与次级用户SU相似但接入信道顺序相反，先接入预留信道然后RP最后CP信道。

因考虑按需的感知机理，所以当一个主用户接入一个正在被次级用户使用的信道时，次级用户不能立刻腾出信道。因此主次用户共享信道，有四种状态：闲置、主用户占用、次级用户占用和主次用户同时占用。用i_CP、i_RP和i_Res分别表示主用户占用CP、RP和预留信道的个数；用j_CP和j_RP分别表示次级用户占用CP和RP信道的个数；用k_CP和k_RP分别表示主次用户同时占用CP和RP信道的个数。用N_CP,IDLE＝γ_CP-i_CP-j_CP-k_CP，N_RP,IDLE＝γ_RP-i_RP-j_RP-k_RP和N_Res,IDLE＝γ_Res-i_Res分别表示CP、RP和预留信道的空闲信道个数。如附图3所示，本发明的MS-DSA协议的算法概述如下：

【性能评估】

对本发明的基于认知无线电的多策略动态频谱接入方法性能的评估，是评估次级系统中次级用户的掉话率、阻塞率、吞吐量、对主用户的干扰量以及效用函数，次级系统根据每个次级用户的服务质量表现(QOS)来进行系统性能评估。

下面利用多维CTMC模型来模拟本发明提出的MS-DSA方法，表示为{S(t),t≥0}...。假设主用户和次级用户的到达率分别为泊松分布λ_P、λ_S，主次用户的服务时间服从指数分布，分别为μ_P、μ_S。假设所有信道为同类信道，因此次级用户的吞吐量与次级用户传输数量和相应的占用时间成正比。

状态S(t)可以表示为一个7维向量s＝(i_CP,i_RP,i_Res,j_CP,j_RP,k_CP,k_RP)，状态空间为Ω＝{s}，

$s.t.\left\{\begin{array}{c}{i}_{\mathrm{CP}},i_{\mathrm{RP}},i_{\mathrm{Res}},j_{\mathrm{CP}},j_{\mathrm{RP}},k_{\mathrm{CP}},k_{\mathrm{RP}}\≥0\\ i_{\mathrm{CP}}+j_{\mathrm{CP}}+k_{\mathrm{CP}}\≤{\mathrm{\γ}}_{\mathrm{CP}}\\ i_{\mathrm{RP}}+j_{\mathrm{RP}}+k_{\mathrm{RP}}\≤{\mathrm{\γ}}_{\mathrm{RP}}\\ i_{\mathrm{Res}}\≤{\mathrm{\γ}}_{\mathrm{Res}}\end{array}\right.$

为了方便起见，定义s'＝f(s,d(s))∈Ω来反应状态s与状态s'的关系，其中d(s)表示s的状态变化， $d\left(s\right)\∈\{i_{\mathrm{CP}}^{\±},i_{\mathrm{RP}}^{\±},i_{\mathrm{Res}}^{\±},j_{\mathrm{CP}}^{\±},k_{\mathrm{CP}}^{\±},k_{\mathrm{RP}}^{\±}\}.$

用Re(s,s')表示在事件e下引起的状态s到状态s'的变化率,e∈{PUArrival(PUA),PU Departure(PUD),SU Arrival(SUA),SU Departure(SUD)},下面讨论这四种事件发生下的状态变化。

A.e＝PUA

当事件e＝PUA发生，一个到达的主用户会产生下列5中状态变化。

·第一种：当i_Res<γ_Res，到达的主用户直接接入一个空闲的预留信道，状态由s变为

·第二种：当i_Res＝γ_Res而且i_RP+k_RP<γ_RP，到达的主用户以概率

N_RP,IDLE/(γ_RP-i_RP-k_RP)随机接入一个空闲的RP信道，或者以概率j_RP/(γ_RP-i_RP-k_RP)接入一个被次级用户占用的RP信道，他们分别使状态由s变为和 $f(s,j_{\mathrm{RP}}^{-},k_{\mathrm{RP}}^{+});$

·第三种：当i_Res＝γ_Res,i_RP+k_RP＝γ_RP而且N_CP,IDLE>0，到达的主用户直接接入一个空闲的CP信道，使状态转变为

·第四种：当i_Res＝γ_Res,i_RP+k_RP＝γ_RP而且N_CP,IDLE＝0,j_CP>0，到达的主用户接入一个被次级用户所占用的信道，从而令状态变为

·第五种：当i_CP+i_RP+i_Res+k_CP+k_RP＝N，所有信道都被中用户占用，这时新到达的主用户被阻塞。

总结，由主用户到达引起的状态转移率可以写为：

$R_{\mathrm{PUA}}(s,s^{\&prime;})=\left\{\begin{array}{cc}\begin{array}{c}{\mathrm{\&lambda;}}_p\\ \frac{N_{\mathrm{RP},\mathrm{IDLE}}}{({\mathrm{\&gamma;}}_{\mathrm{RP}}-i_{\mathrm{RP}}-k_{\mathrm{Rp}})}{\mathrm{\&lambda;}}_p\\ \frac{j_{\mathrm{RP}}}{({\mathrm{\&gamma;}}_{\mathrm{RP}}-i_{\mathrm{RP}}-k_{\mathrm{RP}})}{\mathrm{\&lambda;}}_p\\ {\mathrm{\&lambda;}}_p\\ {\mathrm{\&lambda;}}_p\\ 0\end{array}& \begin{array}{c}{s}^{\&prime;}=f(s,i_{\mathrm{Res}}^{+}),i_{\mathrm{Res}}\&le;{\mathrm{\&gamma;}}_{\mathrm{Res}}\\ s^{\&prime;}=f(s,i_{\mathrm{RP}}^{+}),i_{\mathrm{Res}}={\mathrm{\&gamma;}}_{\mathrm{Res}},i_{\mathrm{RP}}+k_{\mathrm{RP}}\&le;{\mathrm{\&gamma;}}_{\mathrm{RP}}\\ s^{\&prime;}=f(s,j_{\mathrm{RP}}^{-},k_{\mathrm{RP}}^{+}),i_{\mathrm{Res}}={\mathrm{\&gamma;}}_{\mathrm{Res}},i_{\mathrm{RP}}+k_{\mathrm{RP}}\&le;{\mathrm{\&gamma;}}_{\mathrm{RP}}\\ s^{\&prime;}=f(s,i_{\mathrm{CP}}^{+}),i_{\mathrm{Res}}={\mathrm{\&gamma;}}_{\mathrm{Res}},\\ i_{\mathrm{RP}}+k_{\mathrm{RP}}={\mathrm{\&gamma;}}_{\mathrm{RP}},N_{\mathrm{CP},\mathrm{IDLE}}>0\\ s^{\&prime;}=f(s,j_{\mathrm{CP}}^{-},k_{\mathrm{CP}}^{+}),i_{\mathrm{Res}}={\mathrm{\&gamma;}}_{\mathrm{Res}},\\ i_{\mathrm{Rp}}+k_{\mathrm{Rp}}={\mathrm{\&gamma;}}_{\mathrm{RP}},N_{\mathrm{CP},\mathrm{IDLE}}=0,j_{\mathrm{CP}}>0\\ \mathrm{otherwise}\end{array}\end{array}\right.$

B.e＝SUA

次级用户到达，SAP首先感知当前频谱状态。但由于算法、信噪比、噪声的不确定性和样本空间，频谱感知还不够完善。不失一般性，所有这些因素可以反应为关于虚警概率P_f和漏检概率P_e的ROC曲线。为方便起见，用δ＝(P_f,P_e)表示感知性能。当SAP需要感知所有N条信道中有多少信道被主用户(i)、次级用户(j)或者主次共同(k)占用时，感知对m条空闲信道的概率和感知错n条空闲信道的概率可以写为：

$P_{\mathrm{IDLE}}^{\mathrm{right}}\left(m\right|\mathrm{\&delta;},\underset{\&OverBar;}{N},\underset{\&OverBar;}{i},\underset{\&OverBar;}{j},\underset{\&OverBar;}{k})=C_{\underset{\&OverBar;}{N_{\mathrm{IDLE}}}}^m{(1-P_f)}^m{P}_f^{\underset{\&OverBar;}{N_{\mathrm{IDLE}}}-m}$

$P_{\mathrm{IDLE}}^{\mathrm{wrong}}\left(n\right|\mathrm{\&delta;},\underset{\&OverBar;}{N},\underset{\&OverBar;}{i},\underset{\&OverBar;}{j},\underset{\&OverBar;}{k})=C_{\underset{\&OverBar;}{i}}^n{P}_e^n{(1-P_e)}^{\underset{\&OverBar;}{i}-n}$

其中N _IDLE ＝N-i-j-k,0≤m≤N _IDLE ,0≤n≤i。

令m＝0,n＝0，可以获得SAP确定所有信道忙的概率，表示为：

$\begin{array}{c}{P}_{\mathrm{BUSY}}(\mathrm{\&delta;},\underset{\&OverBar;}{N},\underset{\&OverBar;}{i},\underset{\&OverBar;}{j},\underset{\&OverBar;}{k})=P_{\mathrm{IDLE}}^{\mathrm{right}}\left(0\right|\mathrm{\&delta;},\underset{\&OverBar;}{N},\underset{\&OverBar;}{i},\underset{\&OverBar;}{j},\underset{\&OverBar;}{k})\&CenterDot;P_{\mathrm{IDLE}}^{\mathrm{wrong}}\left(0\right|\mathrm{\&delta;},\underset{\&OverBar;}{N},\underset{\&OverBar;}{i},\underset{\&OverBar;}{j},\underset{\&OverBar;}{k})\\ =P_f^{\underset{\&OverBar;}{N_{\mathrm{IDLE}}}}{(1-P_e)}^{\underset{\&OverBar;}{i}}\end{array}$

假设SAP以等概率将感知到的空闲信道分配给新到达的次级用户。因此，次级用户接入正确和错误信道的概率为：

$P_{\mathrm{access}}^{\mathrm{right}}(\mathrm{\&delta;},\underset{\&OverBar;}{N},\underset{\&OverBar;}{i},\underset{\&OverBar;}{j},\underset{\&OverBar;}{k})=\underset{m=1}{\overset{\underset{\&OverBar;}{N_{\mathrm{IDLE}}}}{\mathrm{\&Sigma;}}}{P}_{\mathrm{IDLE}}^{\mathrm{right}}\left(m\right|\mathrm{\&delta;},\underset{\&OverBar;}{N},\underset{\&OverBar;}{i},\underset{\&OverBar;}{j},\underset{\&OverBar;}{k})\underset{n=0}{\overset{\underset{\&OverBar;}{i}}{\mathrm{\&Sigma;}}}\frac{m}{n+m}{P}_{\mathrm{IDLE}}^{\mathrm{wrong}}\left(n\right|\mathrm{\&delta;},\underset{\&OverBar;}{N},\underset{\&OverBar;}{i},\underset{\&OverBar;}{j},\underset{\&OverBar;}{k})$

$P_{\mathrm{access}}^{\mathrm{right}}(\mathrm{\&delta;},\underset{\&OverBar;}{N},\underset{\&OverBar;}{i},\underset{\&OverBar;}{j},\underset{\&OverBar;}{k})=\underset{m=1}{\overset{\underset{\&OverBar;}{N_{\mathrm{IDLE}}}}{\mathrm{\&Sigma;}}}{P}_{\mathrm{IDLE}}^{\mathrm{right}}\left(m\right|\mathrm{\&delta;},\underset{\&OverBar;}{N},\underset{\&OverBar;}{i},\underset{\&OverBar;}{j},\underset{\&OverBar;}{k})\underset{n=0}{\overset{\underset{\&OverBar;}{i}}{\mathrm{\&Sigma;}}}\frac{m}{n+m}{P}_{\mathrm{IDLE}}^{\mathrm{wrong}}\left(n\right|\mathrm{\&delta;},\underset{\&OverBar;}{N},\underset{\&OverBar;}{i},\underset{\&OverBar;}{j},\underset{\&OverBar;}{k})$

基于此，事件e＝SUA的发生，会产生以下五种情况：

·第一种：以概率正确的接入一个CP信道，使状态变为

·第二种：以概率错误的接入一个CP信道，使状态变为

·第三种：以概率P_BUSY(δ,γ_CP,i_CP,j_CP,k_CP)正确的接入一个RP信道，使状态变为

·第四种：以概率P_BUSY(δ,γ_CP,i_CP,j_CP,k_CP)错误的接入一个RP信道，使状态变为

·第五种：当感知结果为信道全忙时，新到达次级用户被阻塞。

总结，由次级用户到达引起的状态转移率可以写为：

$R_{\mathrm{SUA}}^{\mathrm{Access}}(s,s^{\&prime;}|\mathrm{\&delta;})=\left\{\begin{array}{cc}\begin{array}{c}{P}_{\mathrm{access}}^{\mathrm{right}}(\mathrm{\&delta;},{\mathrm{\&gamma;}}_{\mathrm{CP}},i_{\mathrm{CP}},j_{\mathrm{CP}},k_{\mathrm{CP}}){\mathrm{\&lambda;}}_s\\ P_{\mathrm{access}}^{\mathrm{wrong}}(\mathrm{\&delta;},{\mathrm{\&gamma;}}_{\mathrm{CP}},i_{\mathrm{CP}},j_{\mathrm{CP}},k_{\mathrm{CP}}){\mathrm{\&lambda;}}_s\\ P_{\mathrm{BUSY}}(\mathrm{\&delta;},{\mathrm{\&gamma;}}_{\mathrm{CP}},i_{\mathrm{CP}},j_{\mathrm{CP}},k_{\mathrm{CP}})\&times;\\ P_{\mathrm{access}}^{\mathrm{right}}(\mathrm{\&delta;},{\mathrm{\&gamma;}}_{\mathrm{CP}},i_{\mathrm{CP}},j_{\mathrm{CP}},k_{\mathrm{CP}}){\mathrm{\&lambda;}}_s\\ P_{\mathrm{BUSY}}(\mathrm{\&delta;},{\mathrm{\&gamma;}}_{\mathrm{CP}},i_{\mathrm{CP}},j_{\mathrm{CP}},k_{\mathrm{CP}})\&times;\\ P_{\mathrm{access}}^{\mathrm{wrong}}(\mathrm{\&delta;},{\mathrm{\&gamma;}}_{\mathrm{CP}},i_{\mathrm{CP}},j_{\mathrm{CP}},k_{\mathrm{CP}}){\mathrm{\&lambda;}}_s\\ 0\end{array}& \begin{array}{c}{s}^{\&prime;}=f(s,j_{\mathrm{CP}}^{+}),\&ForAll;s\\ s^{\&prime;}=f(s,i_{\mathrm{CP}}^{-},k_{\mathrm{CP}}^{+}),\&ForAll;s\\ s^{\&prime;}=f(s,j_{\mathrm{RP}}^{+}),\&ForAll;s\\ s^{\&prime;}=f(s,i_{\mathrm{RP}}^{-},k_{\mathrm{RP}}^{+}),\&ForAll;s\\ \mathrm{otherwise}\end{array}\end{array}\right.$

此外，为了避免在连续的两次感知动作间对主用户产生干扰，SAP不得不感知被次级用户占用的信道。定义和为次级用户从CP和RP信道中错误的掉话数，和为次级用户从CP和RP信道中正确的掉话数。因此，根据参数次级用户被挤掉的概率可表示为：

$\begin{array}{c}{P}_{\mathrm{SUI}}(s,s^{\&prime;}|\mathrm{\&delta;})=C_{j_{\mathrm{CP}}}^{{\mathrm{\&zeta;}}_{\mathrm{jCP}}}{C}_{j_{\mathrm{RP}}}^{{\mathrm{\&zeta;}}_{\mathrm{jRP}}}{P}_f^{{\mathrm{\&zeta;}}_{\mathrm{jCP}}+{\mathrm{\&zeta;}}_{\mathrm{jRP}}}{(1-P_f)}^{(j_{\mathrm{CP}}+j_{\mathrm{RP}})-({\mathrm{\&zeta;}}_{\mathrm{iCP}}+{\mathrm{\&zeta;}}_{\mathrm{jRP}})}\&times;\\ C_{k_{\mathrm{CP}}}^{{\mathrm{\&zeta;}}_{k_{\mathrm{CP}}}}{C}_{k_{\mathrm{RP}}}^{{\mathrm{\&zeta;}}_{k_{\mathrm{RP}}}}{(1-P_e)}^{{\mathrm{\&zeta;}}_{k_{\mathrm{CP}}}+{\mathrm{\&zeta;}}_{k_{\mathrm{RP}}}}{P}_e^{(k_{\mathrm{CP}}+k_{\mathrm{RP}})-({\mathrm{\&zeta;}}_{k_{\mathrm{CP}}}+{\mathrm{\&zeta;}}_{k_{\mathrm{RP}}})}\end{array}$

因此，事件次级用户到达产生的状态转移率可推导为：

$R_{\mathrm{SUA}}(s,s^{\&prime;}|\mathrm{\&delta;})=\underset{s^{\&prime;\&prime;}\&Element;\mathrm{\&Omega;}}{\mathrm{\&Sigma;}}{R}_{\mathrm{SUA}}^{\mathrm{Access}}(s,s^{\&prime;\&prime;}|\mathrm{\&delta;})\&CenterDot;P_{\mathrm{SUI}}(s^{\&prime;\&prime;},s^{\&prime;}|\mathrm{\&delta;})$

C.e＝PUDande＝SUD

在没有主用户和次级用户到达的情况，由主用户和次级用户离开所产生的状态转移率可写为：

$R_{\mathrm{PUD}}(s,s^{\&prime;})=\left\{\begin{array}{cc}{i}_{\mathrm{CP}}{\mathrm{\&mu;}}_p& s^{\&prime;}=f(s,i_{\mathrm{CP}}^{-}),\&ForAll;s\\ i_{\mathrm{RP}}{\mathrm{\&mu;}}_p& s^{\&prime;}=f(s,i_{\mathrm{RP}}^{-}),\&ForAll;s\\ i_{\mathrm{Res}}{\mathrm{\&mu;}}_p& s^{\&prime;}=f(s,i_{\mathrm{Res}}^{-}),\&ForAll;s\\ i_{\mathrm{CP}}{\mathrm{\&mu;}}_p& s^{\&prime;}=f(s,j_{\mathrm{CP}}^{+},k_{\mathrm{CP}}^{-}),\&ForAll;s\\ i_{\mathrm{RP}}{\mathrm{\&mu;}}_p& s^{\&prime;}=f(s,j_{\mathrm{RP}}^{+},k_{\mathrm{RP}}^{-}),\&ForAll;s\\ 0& \mathrm{otherwise}\end{array}\right.R_{\mathrm{SUD}}(s,s^{\&prime;})=\left\{\begin{array}{cc}{j}_{\mathrm{CP}}{\mathrm{\&mu;}}_s& s^{\&prime;}=f(s,j_{\mathrm{CP}}^{-}),\&ForAll;s\\ j_{\mathrm{RP}}{\mathrm{\&mu;}}_s& s^{\&prime;}=f(s,j_{\mathrm{RP}}^{-}),\&ForAll;s\\ k_{\mathrm{CP}}{\mathrm{\&mu;}}_s& s^{\&prime;}=f(s,i_{\mathrm{CP}}^{+},k_{\mathrm{CP}}^{-}),\&ForAll;s\\ k_{\mathrm{RP}}{\mathrm{\&mu;}}_s& s^{\&prime;}=f(s,i_{\mathrm{RP}}^{+},k_{\mathrm{Rp}}^{-}),\&ForAll;s\\ 0& \mathrm{otherwise}\end{array}\right.$

总的来说，可以获得在CTMC下MS-DSA模型的无穷小生成矩阵，如下：

$Q=\left\{\begin{array}{c}\begin{array}{cc}{R}_{\mathrm{PUA}}+R_{\mathrm{PUD}}+R_{\mathrm{SUA}}+R_{\mathrm{SUD}}& s\&NotEqual;s^{\&prime;}\end{array}\\ Q(s,s)=-\underset{s^{\&prime;}\&NotEqual;s}{\mathrm{\&Sigma;}}Q(s,s^{\&prime;})\end{array}\right.$

这里Q(s,s')代表状态由s到s'的转移率，而且来确保Q的行和始终为0。

根据Rate Boundedness Assumption,RBA(转移率有界假设)和Positive Recurrent Assumption，PRA(正递归假设)。通过ΠQ＝0,andΠe＝1来求取稳态概率。其中Π＝{π(s),s∈Ω}是由稳态概率按序组成的行向量，同时e＝[1,1,...,1]^T。特别的，用e＝1代换Q的第一列得到Q₁，设行向量b＝[100...]。所以可以得出ΠQ₁＝b，这样稳态概率就可以表示为：

1)阻塞率：所有到达的次级用户中被阻塞用户占有的比例。表示为：

$\begin{array}{c}{P}_{\mathrm{Block}}=\underset{s\&Element;\mathrm{\&Omega;}}{\mathrm{\&Sigma;}}{P}_{\mathrm{BUSY}}(\mathrm{\&delta;},{\mathrm{\&gamma;}}_{\mathrm{CP}},i_{\mathrm{CP}},j_{\mathrm{CP}},k_{\mathrm{CP}})\&CenterDot;P_{\mathrm{BUSY}}(\mathrm{\&delta;},{\mathrm{\&gamma;}}_{\mathrm{RP}},i_{\mathrm{RP}},j_{\mathrm{Rp}},k_{\mathrm{Rp}})\&CenterDot;\mathrm{\&pi;}\left(s\right)\\ =\underset{s\&Element;\mathrm{\&Omega;}}{\mathrm{\&Sigma;}}{P}_f^{(N_{\mathrm{RP},\mathrm{IDLE}}+N_{\mathrm{CP},\mathrm{IDLE}})}\&CenterDot;{(1-P_e)}^{(i_{\mathrm{Rp}}+i_{\mathrm{CP}})}\&CenterDot;\mathrm{\&pi;}\left(s\right)\end{array}$

2)掉话率：所有到达的次级用户中被强制掉话的比例。表示为：

$P_{\mathrm{Drop}}=\frac{\underset{s\&Element;\mathrm{\&Omega;}}{\mathrm{\&Sigma;}}\underset{s^{\&prime;}\&Element;\mathrm{\&Omega;}}{\mathrm{\&Sigma;}}({\mathrm{\&zeta;}}_{j_{\mathrm{CP}}}+{\mathrm{\&zeta;}}_{j_{\mathrm{RP}}}+{\mathrm{\&zeta;}}_{k_{\mathrm{CP}}}+{\mathrm{\&zeta;}}_{k_{\mathrm{RP}}})P_{\mathrm{SUI}}(s,s^{\&prime;}|\mathrm{\&delta;})\&CenterDot;\mathrm{\&pi;}\left(s\right)}{1-P_{\mathrm{Block}}}$

其中 $s^{\&prime;}=(i_{\mathrm{CP}}+{\mathrm{\&zeta;}}_{k_{\mathrm{CP}}},i_{\mathrm{RP}}+{\mathrm{\&zeta;}}_{k_{\mathrm{RP}}},i_{\mathrm{Res}},j_{\mathrm{CP}}-{\mathrm{\&zeta;}}_{j_{\mathrm{CP}}},j_{\mathrm{RP}}-{\mathrm{\&zeta;}}_{j_{\mathrm{RP}}},k_{\mathrm{CP}}-{\mathrm{\&zeta;}}_{k_{\mathrm{CP}}},k_{\mathrm{RP}}-{\mathrm{\&zeta;}}_{k_{\mathrm{RP}}})$

3)有效吞吐量：每单位时间内成功完成传输的次级用户。可以通过掉话率和阻塞率得到，表示为：

η＝λ_s(1-P_Block)(1-P_Drop)

4)干扰量：每单位时间内次级用户对主用户的累积干扰时间。根据PASTA性得出，表示为：

${\mathrm{\&tau;}}_{\mathrm{INT}}=\underset{s\&Element;{\mathrm{\&Omega;}}_{\mathrm{MS}}}{\mathrm{\&Sigma;}}\underset{t\&RightArrow;\&infin;}{\lim }\frac{1}{t}{\&Integral;}_0^{\mathrm{\&tau;}}(k_{\mathrm{CP}}+k_{\mathrm{RP}})I_s\left(\mathrm{\&tau;}\right)\mathrm{d\&tau;}=\underset{s\&Element;\mathrm{\&Omega;}}{\mathrm{\&Sigma;}}(k_{\mathrm{CP}}+k_{\mathrm{RP}})\&CenterDot;\mathrm{\&pi;}\left(s\right)$

其中I_s(τ)是指标函数，当S(τ)＝s，则I_s(τ)＝1，其他情况I_s(τ)＝0备注：当频谱感知足够准确，如果s'＝(i_CP+k_CP,i_RP+k_RP,i_Res,j_CP,j_RP,0,0)则P_SUI(s,s')≈1，其他P_SUI(s,s')≈0。所以可以根据阻塞率和掉话率得出的τ_INT≈P_Drop(1-P_Block)来估计干扰量。因此，次级系统可以根据掉话率和阻塞率估计出干扰量。

5)效用函数：基于上述的性能指标，提出次级系统的效用函数，如下：

u＝η·ω_SU-τ_INT·ω_INT-γ_CP·ω_CP-γ_RP·ω_RP

其中ω_SU代表一次成功的次级用户服务的服务价格；ω_INT表示对主用户干扰的惩罚代价；ω_CP和ω_RP为每个RP和CP信道单位时间的租用价格。

【系统参数优化】

定义向量表示在MS-DSA系统中DSA策略的信道分配情况。相应地，最佳的信道分配方案会使效用函数最大化，即：

γ^*＝argmaxu

因为当每个DSA策略对应的信道数量变化时，整个状态空间就会产生变化，而且所有的性能指标也需要重新推导，所以通过优化方案很难获取γ*。因此，采用通过估计复杂度O(N^d)的穷举法来取得优化方案，其中N代表共享信道个数、d代表使用几种DSA策略。只有当d有限且足够小的时候，才能估算出其复杂度。而这里d≤2，可以使用穷举法求得。

就前面提及的，可以由ROC曲线来得知频谱感知的性能。鉴于感知算法、信噪比和漏检概率Pe是关于虚警概率Pf的单调递减函数。一方面，它带来了一种新的可以改善系统性能的方法。特别的，当P_f变大，频谱感知趋向保守，这样会使降低到达次级用户与主用户发生碰撞的风险，同时也会导致传输效率的降低。在另一方面，虚警概率P_f很小会使频谱感知更加积极，从而通过承担更大的干扰代价来获取更多的频谱资源。

尝试寻找最佳的频谱感知策略来使整个系统的效用函数最大化。因此，P_e和P_f的映射问题，这个问题也可以理解为寻找最佳P_f的问题。

$P_f^{*}=\arg \max u$

改写效用函数：

u＝∑π(s)r(s)-C＝ΠR-C

其中Π＝{π(s)},R＝{r(s)}T,这里C＝γ_CP·ω_CP+γ_RP·ω_RP与P_f相互独立，r(s)代表关于状态s的奖励函数。可以推导出r(s)的表达式：

$\begin{array}{c}r\left(s\right)=(1-P_f^{(N_{\mathrm{Rp},\mathrm{IDLE}}+N_{\mathrm{CP},\mathrm{IDLE}})}\&CenterDot;{(1-P_e)}^{(i_{\mathrm{RP}}+i_{\mathrm{CP}})})-\\ \underset{s^{\&prime;}\&Element;\mathrm{\&Omega;}}{\mathrm{\&Sigma;}}({\mathrm{\&zeta;}}_{j_{\mathrm{CP}}}+{\mathrm{\&zeta;}}_{j_{\mathrm{RP}}}+{\mathrm{\&zeta;}}_{k_{\mathrm{CP}}}+{\mathrm{\&zeta;}}_{k_{\mathrm{RP}}})P_{\mathrm{SUI}}(s,s^{\&prime;}|\mathrm{\&delta;}){\mathrm{\&omega;}}_{\mathrm{SU}}-(k_{\mathrm{CP}}+k_{\mathrm{RP}}){\mathrm{\&omega;}}_{\mathrm{INT}}\end{array}$

然后利用牛顿法来解决这个问题：

${\left(P_f\right)}^{t+1}={\left(P_f\right)}^t+{\left({\&dtri;}_{P_f}^{2}u\right)}^{-1}{\&dtri;}_{P_f}u$

其中t是迭代次数。在经过有限次数的迭代，获得P_f的近似最优解。

采用牛顿法，求出效用函数关于u的一阶和二阶导数是非常必要的。但是，因为Π和P_f之间没有明确关系，u关于P_f的导数无法获得。为了解决这个问题，用Q和R关于P_f的导数来表示Π的导数。

这样效用函数的导数可以表述为：

$\&dtri;u=\&prod;(\&dtri;Q(-{\left(Q\right)}^{\#}R)+\&dtri;R)$

${\&dtri;}^{2}u=2\&prod;\&dtri;{\mathrm{QQ}}^{\#}(\&dtri;Q{\left(Q\right)}^{\#}R-\&dtri;R)-\&prod;({\&dtri;}^2{\mathrm{QQ}}^{\#}R-{\&dtri;}^{2}R)$

其中Q^#是Q的群逆，满足Q^#＝(Q-eΠ)^-1+eΠ。此时利用牛顿法，就可求出P_f的近似最优解。

本发明的基于认知无线电的多策略动态频谱接入方法，相比于现有技术，打破了认知网络只能运行单一的动态频谱接入策略的限制，通过与单一策略的仿真对比，以复杂度增加的代价获取了吞吐量及效用的提高。

以上内容是结合具体的优选实施方式对本发明所作的进一步详细说明，不能认定本发明的具体实施只局限于这些说明。对于本发明所属技术领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干简单推演或替换，都应当视为属于本发明的保护范围。

Claims (2)

1.一种基于认知无线电的多策略动态频谱接入(MS-DSA)方法，应用于认知无线电系统，采用随机抢先DSA(RP-DSA)和协作DSA(CP-DSA)策略，所述方法包括以下步骤：

S1：通过频谱租赁与拍卖SLAB完成交易，获得频谱，初始化信道总数

N，预留信道个数γ_Res、CP信道个数γ_CP和RP信道个数γ_RP；

S2：判断新到达的用户的类别，若为主用户PU，则转到步骤S3，若为次级用户SU，则转到步骤S10；

S3：若PU占用预留信道的个数小于等于γ_Res，则转到S4，否则，转到步骤S5；

S4：PU随机接入一个空闲的预留信道；

S5：若PU占用RP信道的个数与PU和SU同时占用RP信道的个数之和小于等于γ_RP，则转到S6，否则，转到步骤S7；

S6：PU按照RP-DSA策略接入一个RP信道；

S7：若PU占用CP信道的个数与PU和SU同时占用CP信道的个数之和小于等于CP信道的总数，则转到步骤8，否则，转到步骤S9；

S8：按照RP-DSA策略接入一个RP信道；

S9：主用户PU被阻塞；

S10：若SU占用CP信道的个数与PU和SU同时占用CP信道的个数之和小于等于γ_CP，则转到S11；

S11：次级用户访问点SAP感知CP信道，若感知到有空闲的CP信道，则SU按照CP策略接入一个CP信道，否则，转到S12；

S12：若SU占用RP信道的个数与PU和SU同时占用RP信道的个数之和小于等于γ_RP，则转到S13；

S13：SAP感知RP信道，若感知到有空闲的RP信道，则SU按照RP策略接入一个RP信道；否则，转到S14；

S14：次级用户SU被阻塞；

S15：SAP感知被SU所占用的信道，若检测到PU在该信道上传输，则SU从正在传输的信道中离开。

2.一种认知无线电系统，其特征在于：所述认知无线电系统包括主系统、次级系统和频谱租赁和拍卖SLAB；所述主系统将下行授权频谱资源划分为两个部分，一部分是仅为主用户传输留下的预留信道，另一部分是与次级用户共享的使用动态频谱分配技术的信道，主系统通过SLAB暂时将共享信道中的频谱使用权租赁给次级用户，并从次级系统获取回报；所述次级系统包括：DSA策略数据库、频谱获取模块、频谱感知模块、接入控制模块和性能评估模块；在一个协调器的管理下次级系统建立一个次级用户访问点SAP，SAP发送给所述SLAB一个包括频谱需求、支持的动态频谱接入DSA策略、频谱租用价格和对主用户干扰的惩罚的请求，所述认知无线电系统采用如权利要求1所述的多策略动态频谱接入方法为新到达的用户分配频谱，一旦获得SLAB的许可，次级系统便暂时获得接入许可。