发明名称 认知无线网络中基于H无穷滤波理论的协作频谱感知方法
摘要 一种认知无线网络中基于H无穷滤波理论的协作频谱感知方法,步骤如下:针对认知用户进行本地频谱检测,然后将感知数据及本地感知结果发送给数据融合中心;数据融合中心根据每个认知用户发送来的感知数据构建信道增益的状态空间模型,并使用H无穷滤波方法估计每个认知用户对应的信道阴影衰落增益的变化情况;利用每个认知用户对应的信道阴影衰落增益估计值计算测试统计量;将该测试统计量与给定的阈值比较,确定哪些认知用户存在拜占庭式网络攻击行为;将存在攻击行为的认知用户数据和本地感知结果在数据融合中心进行删除;数据融合中心利用或(OR)融合准则对认知用户的本地频谱感知结果进行决策融合,得到协作频谱感知结果。本方法能够更好地实现对授权用户的协作频谱感知。
申请公布号 CN104348563A 申请公布日期 2015.02.11
申请号 CN201410682946.0 申请日期 2014.11.23
申请人 辽宁大学 发明人 景源;马利;马技;李鹏;杨晓峰;李皓宇;牛斌
分类号 H04B17/30(2015.01)I 主分类号 H04B17/30(2015.01)I
代理机构 沈阳杰克知识产权代理有限公司 21207 代理人 罗莹
主权项 认知无线网络中基于H无穷滤波理论的协作频谱感知方法,其特征在于,具体步骤如下:(1)认知无线网络中R个认知用户进行本地频谱检测:认知用户通过天线接收某一频段内的无线信号。因此,第r个认知用户对应的基带等效离散时间接收信号可以表示为:x<sub>r</sub>[n]=s[n]+v<sub>r</sub>[n],n=1,2,...,M式中s[n]表示授权用户经过信道传输的无线信号采样值,v<sub>r</sub>[n]表示相应的噪声采样值,M为接收信号的采样点数;对于每一个认知用户,计算其接收信号的对应的功率谱估计值X<sub>r</sub>[n]:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>X</mi><mi>r</mi></msub><mo>[</mo><mi>n</mi><mo>]</mo><mfrac><mn>1</mn><mi>N</mi></mfrac><msup><mrow><mo>|</mo><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><msub><mi>x</mi><mi>r</mi></msub><mo>[</mo><mi>k</mi><mo>]</mo><msup><mi>e</mi><mrow><mo>-</mo><mi>j</mi><mn>2</mn><mi>&pi;nk</mi></mrow></msup><mo>|</mo></mrow><mn>2</mn></msup><mo>,</mo><mi>n</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><msub><mi>N</mi><mi>f</mi></msub></mrow>]]></math><img file="FDA0000614240390000011.GIF" wi="779" he="140" /></maths>然后,计算X<sub>r</sub>[n]的峰度统计量S<sub>X</sub>,即<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>S</mi><mi>X</mi></msub><mo>=</mo><mfrac><mrow><mn>1</mn><mo>/</mo><msub><mi>N</mi><mi>f</mi></msub><munderover><mi>&Sigma;</mi><mrow><mi>n</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>N</mi><mi>f</mi></msub></munderover><msup><mrow><mo>(</mo><msub><mi>X</mi><mi>r</mi></msub><mo>[</mo><mi>n</mi><mo>]</mo><mo>-</mo><mi>&mu;</mi><mo>)</mo></mrow><mn>3</mn></msup></mrow><msup><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mo>/</mo><msub><mi>N</mi><mi>f</mi></msub><munderover><mi>&Sigma;</mi><mrow><mi>n</mi><mo>=</mo><mn>1</mn></mrow><msub><mi>N</mi><mi>f</mi></msub></munderover><msup><mrow><mo>(</mo><msub><mi>X</mi><mi>r</mi></msub><mo>[</mo><mi>n</mi><mo>]</mo><mo>-</mo><mi>&mu;</mi><mo>)</mo></mrow><mn>2</mn></msup><mo>)</mo></mrow><mrow><mn>3</mn><mo>/</mo><mn>2</mn></mrow></msup></mfrac></mrow>]]></math><img file="FDA0000614240390000012.GIF" wi="559" he="274" /></maths>式中,μ为X<sub>r</sub>[n]的均值。进而,可以计算得到非高斯性测度测试统计量<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>T</mi><mi>X</mi></msub><mo>=</mo><msubsup><mi>S</mi><mi>X</mi><mn>2</mn></msubsup><mo>;</mo></mrow>]]></math><img file="FDA0000614240390000013.GIF" wi="178" he="70" /></maths>根据给定的虚警概率P<sub>f</sub>,通过反复试验可以得到频谱检测阈值τ。将测试统计量T<sub>X</sub>与检测阈值τ进行比较,如果T<sub>X</sub>大于τ,说明授权用户信号存在,反之,表明授权用户不存在。至此,每个认知用户都可以计算得到相应的本地频谱检测结果{D<sub>r</sub>,r=1,2,...,R}。D<sub>r</sub>=1表示第r个认知用户检测到授权用户存在;D<sub>r</sub>=0表示第r个认知用户检测到授权用户不存在。然后,每个认知用户将感知数据(即,接收信号强度)和本地检测结果发送给数据融合中心;(2)数据融合中心根据每个认知用户发送来的感知数据,建立信道增益的状态空间模型:假设第r个认知用户发送来的N个时刻的感知信号强度表示为:P<sub>r</sub>[n]=P<sub>0</sub>+γ10log(d<sub>0</sub>)‑γ10log(d<sub>r</sub>[n])+h<sub>r</sub>[n]+w<sub>r</sub>[n],n=1,2,...,N;式中,γ表示无线信道的路径衰落系数,d<sub>0</sub>为参考距离,P<sub>0</sub>为参考距离d<sub>0</sub>对应的接收信号强度,d<sub>r</sub>[n]为时刻n第r个认知用户与授权用户之间的距离,h<sub>r</sub>[n]表示时刻n第r个认知用户与授权用户之间的信道阴影衰落增益,w<sub>r</sub>[n]表示环境噪声。该方程也可以称为是信道阴影衰落增益的观测方程;假设信道阴影衰落增益的状态变化可以用一阶自回归模型来表示,即h<sub>r</sub>[n]=ah<sub>r</sub>[n‑1]+v<sub>r</sub>[n],n=1,2,...,N;式中,a为模型系数,v<sub>r</sub>[n]为模型驱动噪声。至此,以上两个公式给出了第r个认知用户与授权用户之间的信道阴影衰落增益的状态空间模型;(3)根据步骤(2)中建立的信道阴影衰落增益的状态空间模型,利用如下H无穷滤波准则,对h<sub>r</sub>[n]进行迭代估计<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mi>sup</mi><mo>{</mo><mi>J</mi><mo>=</mo><mfrac><mrow><munderover><mi>&Sigma;</mi><mrow><mi>n</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msup><mrow><mo>|</mo><msub><mi>h</mi><mi>r</mi></msub><mo>[</mo><mi>n</mi><mo>]</mo><mo>-</mo><msub><mover><mi>h</mi><mo>^</mo></mover><mi>r</mi></msub><mo>[</mo><mi>n</mi><mo>]</mo><mo>|</mo></mrow><mn>2</mn></msup></mrow><mrow><msubsup><mrow><mo>|</mo><msub><mi>h</mi><mi>r</mi></msub><mo>[</mo><mn>1</mn><mo>]</mo><mo>-</mo><msub><mover><mi>h</mi><mo>^</mo></mover><mi>r</mi></msub><mo>[</mo><mn>1</mn><mo>]</mo><mo>|</mo></mrow><msubsup><mi>Q</mi><mn>0</mn><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><mn>2</mn></msubsup><mo>+</mo><munderover><mi>&Sigma;</mi><mrow><mi>n</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><mo>{</mo><msub><mi>w</mi><mi>r</mi></msub><mo>[</mo><mi>n</mi><mo>]</mo><msubsup><mo>|</mo><msup><mi>W</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mn>2</mn></msubsup><mo>+</mo><msubsup><mrow><mo>|</mo><msub><mi>v</mi><mi>r</mi></msub><mo>[</mo><mi>n</mi><mo>]</mo><mo>|</mo></mrow><msup><mi>V</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mn>2</mn></msubsup></mrow></mfrac><mo>}</mo><mo>&le;</mo><msup><mi>&gamma;</mi><mn>2</mn></msup><mo>;</mo></mrow>]]></math><img file="FDA0000614240390000021.GIF" wi="1084" he="268" /></maths>式中,<img file="FDA0000614240390000022.GIF" wi="110" he="79" />为信道阴影衰落增益h<sub>r</sub>[n]的估计值,<img file="FDA0000614240390000023.GIF" wi="179" he="74" />W&gt;0,V&gt;0,γ为给定的噪声水平界限。上述H无穷滤波等价于如下所示极小极大优化问题<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><munder><mi>min</mi><mrow><msub><mover><mi>h</mi><mo>^</mo></mover><mi>r</mi></msub><mo>[</mo><mi>n</mi><mo>]</mo></mrow></munder><munder><mi>max</mi><mrow><msub><mi>w</mi><mi>r</mi></msub><mo>[</mo><mi>n</mi><mo>]</mo><mo>,</mo><msub><mi>v</mi><mi>r</mi></msub><mo>[</mo><mi>n</mi><mo>]</mo><mo>,</mo><msub><mover><mi>h</mi><mo>^</mo></mover><mi>r</mi></msub><mo>[</mo><mn>0</mn><mo>]</mo></mrow></munder><mi>J</mi><mo>=</mo><mo>-</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><msup><mi>&gamma;</mi><mn>2</mn></msup><msubsup><mrow><mo>|</mo><msub><mi>h</mi><mi>r</mi></msub><mo>[</mo><mn>1</mn><mo>]</mo><mo>-</mo><msub><mover><mi>h</mi><mo>^</mo></mover><mi>r</mi></msub><mo>|</mo></mrow><msubsup><mi>Q</mi><mn>0</mn><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><mn>2</mn></msubsup><mo>+</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mo>[</mo><munderover><mi>&Sigma;</mi><mrow><mi>n</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msup><mrow><mo>|</mo><msub><mi>h</mi><mi>r</mi></msub><mo>[</mo><mi>n</mi><mo>]</mo><mo>-</mo><msub><mover><mi>h</mi><mo>^</mo></mover><mi>r</mi></msub><mo>[</mo><mi>n</mi><mo>]</mo><mo>|</mo></mrow><mn>2</mn></msup><mo>-</mo><msup><mi>&gamma;</mi><mn>2</mn></msup><mrow><mo>(</mo><mo>{</mo><mo>|</mo><msub><mi>w</mi><mi>r</mi></msub><mo>[</mo><mi>n</mi><mo>]</mo><msubsup><mo>|</mo><msup><mi>W</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mn>2</mn></msubsup><mo>+</mo><msub><mrow><mo>|</mo><mi>v</mi></mrow><mi>r</mi></msub><mo>[</mo><mi>n</mi><mo>]</mo><msubsup><mo>|</mo><msup><mi>V</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mn>2</mn></msubsup><mo>}</mo><mo>)</mo></mrow><mo>]</mo></mrow>]]></math><img file="FDA0000614240390000024.GIF" wi="1703" he="143" /></maths>通过求解该优化问题,可以得到信道阴影衰落增益的迭代估计为:<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><msub><mover><mi>h</mi><mo>^</mo></mover><mi>r</mi></msub><mo>[</mo><mi>n</mi><mo>]</mo><mo>=</mo><mi>a</mi><msub><mover><mi>h</mi><mo>^</mo></mover><mi>r</mi></msub><mo>[</mo><mi>n</mi><mo>-</mo><mn>1</mn><mo>]</mo><mo>+</mo><mi>H</mi><mo>[</mo><mi>n</mi><mo>]</mo><mo>[</mo><msub><mi>P</mi><mi>r</mi></msub><mo>[</mo><mi>n</mi><mo>]</mo><mo>-</mo><msub><mi>P</mi><mn>0</mn></msub><mo>-</mo><mi>&lambda;</mi><mn>10</mn><mi>log</mi><mrow><mo>(</mo><msub><mi>d</mi><mn>0</mn></msub><mo>)</mo></mrow><mo>-</mo><mi>&lambda;</mi><mn>10</mn><mi>log</mi><mrow><mo>(</mo><msub><mi>d</mi><mi>r</mi></msub><mo>[</mo><mi>n</mi><mo>]</mo><mo>)</mo></mrow><mo>-</mo><msub><mover><mi>h</mi><mo>^</mo></mover><mi>r</mi></msub><mo>[</mo><mi>n</mi><mo>-</mo><mn>1</mn><mo>]</mo><mo>]</mo></mrow>]]></math><img file="FDA0000614240390000025.GIF" wi="1414" he="104" /></maths>式中H[n]=aP[n](V+P[n])<sup>‑1</sup>,P<sup>‑1</sup>[n]=[a(P<sup>‑1</sup>[n]+V<sup>‑1</sup>)<sup>‑1</sup>a+W]<sup>‑1</sup>‑γ<sup>‑2</sup>;(4)对于认知无线网络中的R个认知用户,数据融合中心通过估计得到的每个认知用户对应的信道阴影衰落增益<img file="FDA0000614240390000031.GIF" wi="108" he="80" />在N个时刻内的变化情况,检测存在拜占庭攻击行为的认知用户。该检测准则对应的信道阴影衰落增益变化测试量为<maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><msub><mi>T</mi><mi>A</mi></msub><mo>=</mo><msup><mrow><mo>(</mo><mo>|</mo><msub><mover><mi>h</mi><mo>^</mo></mover><mi>r</mi></msub><mo>[</mo><mi>n</mi><mo>]</mo><mo>-</mo><msub><mover><mi>h</mi><mo>^</mo></mover><mi>r</mi></msub><mo>[</mo><mi>n</mi><mo>-</mo><mn>1</mn><mo>]</mo><msup><mo>|</mo><mn>2</mn></msup><mo>-</mo><mo>|</mo><msub><mover><mi>h</mi><mo>^</mo></mover><mn>1</mn></msub><mo>[</mo><mi>n</mi><mo>-</mo><mn>1</mn><mo>]</mo><mo>-</mo><msub><mover><mi>h</mi><mo>^</mo></mover><mi>r</mi></msub><mo>[</mo><mi>n</mi><mo>-</mo><mn>2</mn><mo>]</mo><msup><mo>|</mo><mn>2</mn></msup><mo>)</mo></mrow><mn>2</mn></msup></mrow>]]></math><img file="FDA0000614240390000032.GIF" wi="901" he="132" /></maths>通过该准则可以看出,T<sub>A</sub>反应了信道阴影衰落增益在短时间内的变化程度。由于实际的信道在短时间内阴影衰落增益不会产生较大的突变,因此当T<sub>A</sub>大于给定阈值时,则认定第r个认知用户存在拜占庭攻击行为,该阈值可以通过反复试验获得;(5)在数据融合中心,将存在拜占庭式攻击行为的认知用户对应的本地频谱感知结果进行删除;(6)在数据融合中心,根据步骤(3)中的H无穷滤波结果以及步骤(4)中的认知用户攻击行为的识别结果,数据融合中心对剩余的认知用户本地频谱感知结果进行或决策融合,得到最终的授权用户协作频谱感知结果。
地址 110000 辽宁省沈阳市沈北新区道义南大街58号
您可能感兴趣的专利