发明名称 基于多进程D-S证据理论的合作频谱感知方法
摘要 本发明基于多进程D-S证据理论的合作频谱感知方法,通过双门限检测方案的两个门限值找出不确定的信息,认知用户再根据D-S证据理论处理不确定信息,最后由簇头节点收集簇内认知用户的感知结果并传给融合中心,融合中心通过多进程D-S证据理论进行融合并判断,这样除去过时感知信息对融合中心正确判断频谱状态的影响,保证判断结果的正确性,同时缩短了判断时间提高了频谱的利用效率。
申请公布号 CN103581922A 申请公布日期 2014.02.12
申请号 CN201310467363.1 申请日期 2013.10.09
申请人 北京科技大学 发明人 刘健;齐爱想;隆克平
分类号 H04W16/14(2009.01)I;H04B17/00(2006.01)I 主分类号 H04W16/14(2009.01)I
代理机构 成都行之专利代理事务所(普通合伙) 51220 代理人 温利平
主权项 1.基于多进程D-S证据理论的合作频谱感知方法,其特征在于,包括以下步骤:(1)、生成多进程D-S证据理论合作频谱感知的帧结构:将第(i-1)个认知用户SU<sub>i-1</sub>的感知结果作为报告与第i,i=1,2,…N个认知用户SU<sub>i</sub>的感知组成新的频谱感知字段,新频谱感知字段的数据传输周期不变,每个认知用户的感知和报告的周期相同均为T,则N个认知用户完成感知和报告的周期为(N+1)T,新频谱感知字段和数据接收、传输字段构成多进程D-S证据理论合作频谱感知的帧结构;(2)、认知用户根据D-S证据理论计算不确定信息的信度函数:(2.1)、计算第i个认知用户接受到的能量<img file="FDA0000392961410000011.GIF" wi="114" he="89" />表达式为:<maths num="0001"><![CDATA[<math><mrow><msub><mi>x</mi><msub><mi>E</mi><mi>i</mi></msub></msub><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>t</mi><mo>=</mo><mn>1</mn></mrow><mrow><mn>2</mn><mi>v</mi></mrow></munderover><msup><mrow><mo>|</mo><msub><mi>y</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>|</mo></mrow><mn>2</mn></msup></mrow></math>]]></maths>其中,v为时间带宽,y<sub>i</sub>(t)为第i个认知用户接收到的信号,其表达式为:<maths num="0002"><![CDATA[<math><mrow><msub><mi>y</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mfenced open='{' close=''><mtable><mtr><mtd><msub><mi>h</mi><mi>i</mi></msub><msub><mi>s</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>n</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><msub><mi>H</mi><mn>1</mn></msub></mtd></mtr><mtr><mtd><msub><mi>n</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><msub><mi>H</mi><mn>0</mn></msub></mtd></mtr></mtable></mfenced></mrow></math>]]></maths>其中,h<sub>i</sub>为信道衰落系数,s<sub>i</sub>(t)为主用户信号,n<sub>i</sub>(t)为加性高斯白噪声,H<sub>0</sub>为主用户信号不存在,H<sub>1</sub>为主用户信号存在;根据中心极限定理,当v足够大时,<img file="FDA0000392961410000014.GIF" wi="83" he="90" />为高斯随机变量,其均值和方差的为:主用户信号不存在时的均值μ<sub>0i</sub>和方差<img file="FDA0000392961410000015.GIF" wi="108" he="82" /><maths num="0003"><![CDATA[<math><mfenced open='{' close='' separators=''><mtable><mtr><mtd><msub><mi>&mu;</mi><mrow><mn>0</mn><mi>i</mi></mrow></msub><mo>=</mo><msub><mi>N</mi><mn>0</mn></msub><mi>v</mi></mtd></mtr><mtr><mtd><msubsup><mi>&sigma;</mi><mrow><mn>0</mn><mi>i</mi></mrow><mn>2</mn></msubsup><mo>=</mo><msubsup><mi>N</mi><mn>0</mn><mn>2</mn></msubsup><mi>v</mi></mtd></mtr></mtable><mo>;</mo></mfenced></math>]]></maths>主用户信号存在时的均值μ<sub>1i</sub>和方差<img file="FDA0000392961410000017.GIF" wi="100" he="87" /><maths num="0004"><![CDATA[<math><mfenced open='{' close='' separators=''><mtable><mtr><mtd><msub><mi>&mu;</mi><mrow><mn>1</mn><mi>i</mi></mrow></msub><mo>=</mo><msub><mi>N</mi><mn>0</mn></msub><mi>v</mi><mrow><mo>(</mo><mi>SNR</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow></mtd></mtr><mtr><mtd><msubsup><mi>&sigma;</mi><mrow><mn>1</mn><mi>i</mi></mrow><mn>2</mn></msubsup><mo>=</mo><msubsup><mi>N</mi><mn>0</mn><mn>2</mn></msubsup><mi>v</mi><mrow><mo>(</mo><mn>2</mn><mi>SNR</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow></mtd></mtr></mtable><mo>;</mo></mfenced></math>]]></maths>其中,N<sub>0</sub>为独立复高斯随机变量的单边功率谱密度,SNR为信噪比;(2.2)、通过双门限检测方案找出不确定信息:计算双门限检测方案的门限值λ<sub>0</sub>和λ<sub>1</sub>;在h<sub>i</sub>确定条件下,第i个认知用户的虚警概率Q<sub>f,i</sub>:<maths num="0005"><![CDATA[<math><mrow><msub><mi>Q</mi><mrow><mi>f</mi><mo>,</mo><mi>i</mi></mrow></msub><mo>=</mo><mi>P</mi><mo>{</mo><msub><mi>x</mi><msub><mi>E</mi><mi>i</mi></msub></msub><mo>></mo><msub><mi>&lambda;</mi><mn>1</mn></msub><mo>|</mo><msub><mi>H</mi><mn>0</mn></msub><mo>}</mo><mo>=</mo><mfrac><mrow><mi>&Gamma;</mi><mrow><mo>(</mo><mi>v</mi><mo>,</mo><msub><mi>&lambda;</mi><mn>1</mn></msub><mo>/</mo><mn>2</mn><mo>)</mo></mrow></mrow><mrow><mi>&Gamma;</mi><mrow><mo>(</mo><mi>v</mi><mo>)</mo></mrow></mrow></mfrac></mrow></math>]]></maths>其中,Γ(ν)为完全伽玛函数;Γ(ν,λ<sub>1</sub>/2)为非完全伽玛函数;Γ(ν,λ<sub>1</sub>/2)=Γ(ν)-Γ(ν)·P(ν,λ<sub>1</sub>/2);<maths num="0006"><![CDATA[<math><mrow><mi>P</mi><mrow><mo>(</mo><mi>v</mi><mo>,</mo><msub><mi>&lambda;</mi><mn>1</mn></msub><mo>/</mo><mn>2</mn><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mrow><mi>&Gamma;</mi><mrow><mo>(</mo><mi>v</mi><mo>)</mo></mrow></mrow></mfrac><munderover><mo>&Integral;</mo><mn>0</mn><mrow><msub><mi>&lambda;</mi><mn>1</mn></msub><mo>/</mo><mn>2</mn></mrow></munderover><msup><mi>t</mi><mrow><mi>v</mi><mo>-</mo><mn>1</mn></mrow></msup><mo>&CenterDot;</mo><msup><mi>e</mi><mrow><mo>-</mo><mi>t</mi></mrow></msup><mi>dt</mi><mo>;</mo></mrow></math>]]></maths>其中P(ν,λ<sub>1</sub>/2)表示的是低阶非完全伽玛函数;由此我们可以得到λ<sub>1</sub>=2P<sup>-1</sup>(ν,1-Q<sub>f,i</sub>);λ<sub>0</sub>=c·λ<sub>1</sub>式中c为常数,可根据实际需要确定该值;根据λ<sub>0</sub>和λ<sub>1</sub>找出不确定信息;第i个认知用户接受到的能量值<img file="FDA0000392961410000023.GIF" wi="174" he="89" />为H<sub>0</sub>,<img file="FDA0000392961410000024.GIF" wi="173" he="89" />为H<sub>1</sub>,<img file="FDA0000392961410000025.GIF" wi="227" he="89" />为不确定信息;(2.3)、通过D-S证据理论计算不确定信息的信度函数:根据D-S证据理论,通过双门限计算第i个感知用户的信度函数m<sub>i</sub>(H<sub>0</sub>)和m<sub>i</sub>(H<sub>1</sub>),<maths num="0007"><![CDATA[<math><mrow><mfenced open='' close=''><mtable><mtr><mtd><msub><mi>m</mi><mi>i</mi></msub><mrow><mo>(</mo><msub><mi>H</mi><mn>0</mn></msub><mo>)</mo></mrow><mo>=</mo><mi>P</mi><mo>{</mo><msub><mi>x</mi><msub><mi>E</mi><mi>i</mi></msub></msub><mo>&lt;</mo><msub><mi>&lambda;</mi><mn>0</mn></msub><mo>|</mo><msub><mi>H</mi><mn>0</mn></msub><mo>}</mo><mo>+</mo><mi>P</mi><mo>{</mo><msub><mi>&lambda;</mi><mn>0</mn></msub><mo>&lt;</mo><msub><mi>x</mi><msub><mi>E</mi><mi>i</mi></msub></msub><mo>&lt;</mo><msub><mi>&lambda;</mi><mn>1</mn></msub><mo>|</mo><msub><mi>H</mi><mn>0</mn></msub><mo>}</mo></mtd></mtr><mtr><mtd><mo>=</mo><munderover><mo>&Integral;</mo><mrow><mo>-</mo><mo>&infin;</mo></mrow><msub><mi>&lambda;</mi><mn>1</mn></msub></munderover><mfrac><mn>1</mn><mrow><msqrt><mn>2</mn><mi>&pi;</mi></msqrt><msub><mi>&sigma;</mi><mrow><mn>0</mn><mi>i</mi></mrow></msub></mrow></mfrac><mi>exp</mi><mo>[</mo><mo>-</mo><msup><mrow><mo>(</mo><mi>x</mi><mo>-</mo><msub><mi>&mu;</mi><mrow><mn>0</mn><mi>i</mi></mrow></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>/</mo><msubsup><mi>&sigma;</mi><mrow><mn>0</mn><mi>i</mi></mrow><mn>2</mn></msubsup><mo>]</mo><mi>dx</mi></mtd></mtr></mtable></mfenced><mo>;</mo></mrow></math>]]></maths><maths num="0008"><![CDATA[<math><mrow><mfenced open='' close=''><mtable><mtr><mtd><msub><mi>m</mi><mi>i</mi></msub><mrow><mo>(</mo><msub><mi>H</mi><mn>1</mn></msub><mo>)</mo></mrow><mo>=</mo><mi>P</mi><mo>{</mo><msub><mi>x</mi><msub><mi>E</mi><mi>i</mi></msub></msub><mo>&lt;</mo><msub><mi>&lambda;</mi><mn>1</mn></msub><mo>|</mo><msub><mi>H</mi><mn>1</mn></msub><mo>}</mo><mo>+</mo><mi>P</mi><mo>{</mo><msub><mi>&lambda;</mi><mn>0</mn></msub><mo>&lt;</mo><msub><mi>x</mi><msub><mi>E</mi><mi>i</mi></msub></msub><mo>&lt;</mo><msub><mi>&lambda;</mi><mn>1</mn></msub><mo>|</mo><msub><mi>H</mi><mn>1</mn></msub><mo>}</mo></mtd></mtr><mtr><mtd><mo>=</mo><munderover><mo>&Integral;</mo><msub><mi>&lambda;</mi><mn>0</mn></msub><mrow><mo>+</mo><mo>&infin;</mo></mrow></munderover><mfrac><mn>1</mn><mrow><msqrt><mn>2</mn><mi>&pi;</mi></msqrt><msub><mi>&sigma;</mi><mrow><mn>1</mn><mi>i</mi></mrow></msub></mrow></mfrac><mi>exp</mi><mo>[</mo><mo>-</mo><msup><mrow><mo>(</mo><mi>x</mi><mo>-</mo><msub><mi>&mu;</mi><mrow><mn>1</mn><mi>i</mi></mrow></msub><mi></mi><mo>)</mo></mrow><mn>2</mn></msup><mo>/</mo><msubsup><mi>&sigma;</mi><mrow><mn>1</mn><mi>i</mi></mrow><mn>2</mn></msubsup><mo>]</mo><mi>dx</mi></mtd></mtr></mtable></mfenced><mo>;</mo></mrow></math>]]></maths>(2.4)、簇头节点收集簇内认知用户的感知结果并传给融合中心:将步骤(2.3)中计算得到的信度函数传给簇头节点CH<sub>h</sub>(h=1,2…,M)作为最终感知结果R<sub>Ch</sub>传给融合中心;(3)、融合中心基于多进程D-S证据理论的合作频谱感知方法进行融合并判决:(3.1)、来自簇头节点的新感知结果到达融合中心后,融合中心便开启一个新的D-S进程;(3.2)、开启的D-S进程对融合数据进行计算并作出判决:(3.2.1)、定义基本概率指派函数m:基本概率指派为2<sup>Ω</sup>到[0,1]的函数m,m满足:<maths num="0009"><![CDATA[<math><mfenced open='{' close=''><mtable><mtr><mtd><mi>m</mi><mrow><mo>(</mo><mi>&Phi;</mi><mo>)</mo></mrow><mo>=</mo><mn>0</mn></mtd></mtr><mtr><mtd><munder><mi>&Sigma;</mi><mrow><mi>A</mi><mo>&SubsetEqual;</mo><mi>&Omega;</mi></mrow></munder><mi>m</mi><mrow><mo>(</mo><mi>A</mi><mo>)</mo></mrow><mo>=</mo><mn>1</mn></mtd></mtr></mtable></mfenced></math>]]></maths>式中m(A)表示的是A的可信度,也叫做证据,Ω表示一个互斥且完整的假设集,称为识别域,分别定义了信度函数bel和似然函数pl,对于所有<img file="FDA0000392961410000038.GIF" wi="169" he="62" /><maths num="0010"><![CDATA[<math><mfenced open='{' close=''><mtable><mtr><mtd><mi>bel</mi><mrow><mo>(</mo><mi>A</mi><mo>)</mo></mrow><mo>=</mo><munder><mi>&Sigma;</mi><mrow><mi>B</mi><mo>&SubsetEqual;</mo><mi>A</mi></mrow></munder><mi>m</mi><mrow><mo>(</mo><mi>B</mi><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mi>pl</mi><mrow><mo>(</mo><mi>A</mi><mo>)</mo></mrow><mo>=</mo><munder><mi>&Sigma;</mi><mrow><mi>A</mi><mo>&cap;</mo><mi>B</mi><mo>&NotEqual;</mo><mi>&phi;</mi></mrow></munder><mi>m</mi><mrow><mo>(</mo><mi>B</mi><mo>)</mo></mrow></mtd></mtr></mtable></mfenced></math>]]></maths>因此[bel(A),pl(A)]就表示了m的不确定区间;(3.2.2)、根据多进程D-S证据理论的融合法则,对任意K个基本概率指派函数进行正交和,构成新的基本概率指派函数m<sub>K</sub>(A<sub>K</sub>),K属于[1,M],M是网络中簇的个数:<maths num="0011"><![CDATA[<math><mfenced open='' close=''><mtable><mtr><mtd><msub><mi>m</mi><mi>K</mi></msub><mrow><mo>(</mo><msub><mi>A</mi><mi>K</mi></msub><mo>)</mo></mrow><mo>=</mo><msub><mi>m</mi><mn>1</mn></msub><mrow><mo>(</mo><msub><mi>A</mi><mn>1</mn></msub><mo>)</mo></mrow><mo>&CirclePlus;</mo><msub><mi>m</mi><mn>2</mn></msub><msub><mi>A</mi><mn>2</mn></msub><mo>&CirclePlus;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><msub><mi>m</mi><mi>M</mi></msub><mrow><mo>(</mo><msub><mi>A</mi><mi>M</mi></msub><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mo>=</mo><mfrac><mrow><munder><mi>&Sigma;</mi><mrow><msub><mi>A</mi><mn>1</mn></msub><mo>&cap;</mo><msub><mi>A</mi><mn>2</mn></msub><mo>&cap;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><msub><mi>A</mi><mi>M</mi></msub><mo>=</mo><mi>A</mi></mrow></munder><msub><mi>m</mi><mn>1</mn></msub><mrow><mo>(</mo><msub><mi>A</mi><mn>1</mn></msub><mo>)</mo></mrow><msub><mi>m</mi><mn>2</mn></msub><mrow><mo>(</mo><msub><mi>A</mi><mn>2</mn></msub><mo>)</mo></mrow><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><msub><mi>m</mi><mi>M</mi></msub><mrow><mo>(</mo><msub><mi>A</mi><mi>M</mi></msub><mo>)</mo></mrow></mrow><mrow><munder><mi>&Sigma;</mi><mrow><msub><mi>A</mi><mn>1</mn></msub><mo>&cap;</mo><msub><mi>A</mi><mn>2</mn></msub><mo>&cap;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><msub><mi>A</mi><mi>M</mi></msub><mo>&NotEqual;</mo><mi>A</mi></mrow></munder><msub><mi>m</mi><mn>1</mn></msub><mrow><mo>(</mo><msub><mi>A</mi><mn>1</mn></msub><mo>)</mo></mrow><msub><mi>m</mi><mn>2</mn></msub><mrow><mo>(</mo><msub><mi>A</mi><mn>2</mn></msub><mo>)</mo></mrow><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><msub><mi>m</mi><mi>M</mi></msub><mrow><mo>(</mo><msub><mi>A</mi><mi>M</mi></msub><mo>)</mo></mrow></mrow></mfrac></mtd></mtr></mtable></mfenced></math>]]></maths>(3.2.3)、对融合后的基本概率指派函数计算:融合中心接收到K个感知结果,对应K个D-S证据理论进程的最终结果<img file="FDA0000392961410000034.GIF" wi="84" he="66" />(K)和<img file="FDA0000392961410000035.GIF" wi="78" he="66" />(K)为:<maths num="0012"><![CDATA[<math><mrow><msub><mi>m</mi><msub><mi>H</mi><mn>0</mn></msub></msub><mrow><mo>(</mo><mi>K</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><munder><mi>&Sigma;</mi><mrow><msub><mi>A</mi><mn>1</mn></msub><mo>&cap;</mo><msub><mi>A</mi><mn>2</mn></msub><mo>&cap;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><msub><mi>A</mi><mi>M</mi></msub><mo>=</mo><msub><mi>H</mi><mn>0</mn></msub></mrow></munder><munderover><mi>&Pi;</mi><mrow><mi>K</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><msub><mi>m</mi><mi>K</mi></msub><mrow><mo>(</mo><msub><mi>A</mi><mi>K</mi></msub><mo>)</mo></mrow></mrow><mrow><mn>1</mn><mo>-</mo><munder><mi>&Sigma;</mi><mrow><msub><mi>A</mi><mn>1</mn></msub><mo>&cap;</mo><msub><mi>A</mi><mn>2</mn></msub><mo>&cap;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><msub><mi>A</mi><mi>M</mi></msub><mo>=</mo><mi>&Theta;</mi></mrow></munder><munderover><mi>&Pi;</mi><mrow><mi>K</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><msub><mi>m</mi><mi>K</mi></msub><mrow><mo>(</mo><msub><mi>A</mi><mi>K</mi></msub><mo>)</mo></mrow></mrow></mfrac></mrow></math>]]></maths><maths num="0013"><![CDATA[<math><mrow><msub><mi>m</mi><msub><mi>H</mi><mn>1</mn></msub></msub><mrow><mo>(</mo><mi>K</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><munder><mi>&Sigma;</mi><mrow><msub><mi>A</mi><mn>1</mn></msub><mo>&cap;</mo><msub><mi>A</mi><mn>2</mn></msub><mo>&cap;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><msub><mi>A</mi><mi>M</mi></msub><mo>=</mo><msub><mi>H</mi><mn>1</mn></msub></mrow></munder><munderover><mi>&Pi;</mi><mrow><mi>K</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><msub><mi>m</mi><mi>K</mi></msub><mrow><mo>(</mo><msub><mi>A</mi><mi>K</mi></msub><mo>)</mo></mrow></mrow><mrow><mn>1</mn><mo>-</mo><munder><mi>&Sigma;</mi><mrow><msub><mi>A</mi><mn>1</mn></msub><mo>&cap;</mo><msub><mi>A</mi><mn>2</mn></msub><mo>&cap;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><mo>&CenterDot;</mo><msub><mi>A</mi><mi>M</mi></msub><mo>=</mo><mi>&Theta;</mi></mrow></munder><munderover><mi>&Pi;</mi><mrow><mi>K</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><msub><mi>m</mi><mi>K</mi></msub><mrow><mo>(</mo><msub><mi>A</mi><mi>K</mi></msub><mo>)</mo></mrow></mrow></mfrac></mrow></math>]]></maths>其中,Θ表示空集;(3.2.4)、将最终的融合结果进行比较作出判断,判断准则为:<maths num="0014"><![CDATA[<math><mrow><msub><mi>H</mi><mn>0</mn></msub><mo>:</mo><msub><mi>m</mi><msub><mi>H</mi><mn>0</mn></msub></msub><mrow><mo>(</mo><mi>K</mi><mo>)</mo></mrow><mo>></mo><msub><mi>m</mi><msub><mi>H</mi><mn>1</mn></msub></msub><mrow><mo>(</mo><mi>K</mi><mo>)</mo></mrow></mrow></math>]]></maths><maths num="0015"><![CDATA[<math><mrow><msub><mi>H</mi><mn>1</mn></msub><mo>:</mo><msub><mi>m</mi><msub><mi>H</mi><mn>1</mn></msub></msub><mrow><mo>(</mo><mi>K</mi><mo>)</mo></mrow><mo>&GreaterEqual;</mo><msub><mi>m</mi><msub><mi>H</mi><mn>0</mn></msub></msub><mrow><mo>(</mo><mi>K</mi><mo>)</mo></mrow></mrow></math>]]></maths>
地址 100083 北京市海淀区学院路30号