发明名称 基于有限随机矩阵的频谱检测方法
摘要 本发明公开了一种基于有限随机矩阵的频谱检测方法,主要解决现有基于特征值频谱检测方法获取检测门限困难和能量频谱检测方法受噪声不确定度影响大的问题。其实现步骤是:(1)检测用户根据所要检测的频段,采集该频段的数据x(n);(2)检测用户根据采集的数据构建协方差矩阵R<sub>x</sub>,并对该矩阵进行Cholesky分解;(3)计算检测统计量T<sub>ξ</sub>(4)分析检测统计量T<sub>ξ</sub>的概率分布,根据T<sub>ξ</sub>的概率分布计算目标虚警概率下的检测门限γ<sub>ξ</sub>;(5)检测用户对检测统计量与检测门限进行比较,判决主用户是否存在。本发明具抗噪声不确定度、检测门限精准、检测性能高的优点,可用于无线通信。
申请公布号 CN103973382B 申请公布日期 2016.04.27
申请号 CN201410211335.8 申请日期 2014.05.19
申请人 西安电子科技大学 发明人 周福辉;李赞;杨鼎;刘向丽;杜军朝;刘惠;高锐;黄海燕;关磊;熊天意;胡伟龙
分类号 H04B17/382(2015.01)I 主分类号 H04B17/382(2015.01)I
代理机构 陕西电子工业专利中心 61205 代理人 王品华;黎汉华
主权项 一种基于有限随机矩阵的频谱检测方法,包括如下步骤:(1)将占用当前频段的用户信号定义为主用户,将检测当前频段上主用户存在与否以试图占用该频段的用户信号定义为检测用户,检测用户只配备单根天线,检测用户根据所要观察的频段,采集该频段的数据x(n),其中n=1,…,N,N为采样的点数;(2)检测用户根据采集的数据,构建检测统计量T<sub>ξ</sub>:(2.1)检测用户根据采集到的数据,构建数据矩阵X和协方差矩阵R<sub>x</sub>,其中数据矩阵X为:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>X</mi><mo>=</mo><mfenced open = "[" close = "]"><mtable><mtr><mtd><mrow><mi>x</mi><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow></mtd><mtd><mrow><mi>x</mi><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow></mtd><mtd><mo>...</mo></mtd><mtd><mrow><mi>x</mi><mrow><mo>(</mo><msub><mi>N</mi><mi>s</mi></msub><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><mi>x</mi><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow></mtd><mtd><mrow><mi>x</mi><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow></mtd><mtd><mo>...</mo></mtd><mtd><mrow><mi>x</mi><mrow><mo>(</mo><msub><mi>N</mi><mi>s</mi></msub><mo>+</mo><mn>1</mn><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><mrow><mi>x</mi><mrow><mo>(</mo><mi>M</mi><mo>)</mo></mrow></mrow></mtd><mtd><mrow><mi>x</mi><mrow><mo>(</mo><mi>M</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow></mrow></mtd><mtd><mo>...</mo></mtd><mtd><mrow><mi>x</mi><mrow><mo>(</mo><mi>M</mi><mo>+</mo><msub><mi>N</mi><mi>s</mi></msub><mo>-</mo><mn>1</mn><mo>)</mo></mrow></mrow></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0000865954450000011.GIF" wi="1068" he="364" /></maths>协方差矩阵为:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>R</mi><mi>x</mi></msub><mo>=</mo><mfrac><mn>1</mn><msub><mi>N</mi><mi>s</mi></msub></mfrac><msup><mi>XX</mi><mi>H</mi></msup><mo>,</mo></mrow>]]></math><img file="FDA0000865954450000012.GIF" wi="357" he="158" /></maths>其中N<sub>s</sub>为每段数据的点数,M为数据段数,(·)<sup>H</sup>为Heimitian转置<sub>;</sub>(2.2)检测用户对协方差矩阵进行Cholesky分解,得到分解后的上三角矩阵,即:R<sub>x</sub>=L<sup>T</sup>L,其中,L为上三角矩阵,其表示为:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>L</mi><mo>=</mo><mfenced open = "[" close = "]"><mtable><mtr><mtd><msub><mi>l</mi><mn>11</mn></msub></mtd><mtd><msub><mi>l</mi><mn>12</mn></msub></mtd><mtd><mo>...</mo></mtd><mtd><msub><mi>l</mi><mrow><mn>1</mn><mi>M</mi></mrow></msub></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><msub><mi>l</mi><mn>22</mn></msub></mtd><mtd><mo>...</mo></mtd><mtd><msub><mi>l</mi><mrow><mn>2</mn><mi>M</mi></mrow></msub></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd><mtd><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mo>...</mo></mtd><mtd><msub><mi>l</mi><mrow><mi>M</mi><mi>M</mi></mrow></msub></mtd></mtr></mtable></mfenced><mo>,</mo></mrow>]]></math><img file="FDA0000865954450000013.GIF" wi="581" he="350" /></maths>其中,l<sub>ij</sub>为上三角矩阵L的第i行j列元素,i=1,…,M,j=1,…,M;(2.3)检测用户根据分解后的上三角矩阵L,构建检测统计量T<sub>ξ</sub>:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><mi>T</mi><mi>&xi;</mi></msub><mo>=</mo><mfrac><msubsup><mi>l</mi><mn>11</mn><mn>2</mn></msubsup><msubsup><mi>l</mi><mrow><mi>M</mi><mi>M</mi></mrow><mn>2</mn></msubsup></mfrac><mo>,</mo></mrow>]]></math><img file="FDA0000865954450000021.GIF" wi="214" he="143" /></maths>其中l<sub>11</sub>为上三角矩阵L的主对角线的第一个元素,l<sub>MM</sub>为上三角矩阵L的主对角线的第M个元素;(3)检测用户根据检测统计量T<sub>ξ</sub>的概率分布,计算检测门限γ<sub>ξ</sub>:<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msub><mi>&gamma;</mi><mi>&xi;</mi></msub><mo>=</mo><mfrac><mi>m</mi><mi>n</mi></mfrac><msubsup><mi>P</mi><mi>F</mi><mrow><mo>-</mo><mn>1</mn></mrow></msubsup><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msub><mi>P</mi><mrow><mi>f</mi><mi>a</mi></mrow></msub><mo>,</mo><mi>m</mi><mo>,</mo><mi>n</mi><mo>)</mo></mrow><mo>,</mo></mrow>]]></math><img file="FDA0000865954450000022.GIF" wi="517" he="118" /></maths>其中,<img file="FDA0000865954450000023.GIF" wi="134" he="78" />为F分布F(m,n)的累积量分布函数P<sub>F</sub>(p,m,n)的逆函数,m=N<sub>s</sub>,n=N<sub>s</sub>‑M+1,P<sub>fa</sub>为虚警概率,取值为(0,1),p为P<sub>F</sub>(p,m,n)的自变量,取值为(0,1);(4)将步骤(2.3)中得到的检测统计量T<sub>ξ</sub>与步骤(3)中得到的检测门限γ<sub>ξ</sub>进行比较,当T<sub>ξ</sub>≥γ<sub>ξ</sub>时,判决为主用户存在,即当前频段频谱已被某用户占用,否则,判决为主用户不存在,即当前频段频谱为空闲状态,允许检测用户利用。
地址 710071 陕西省西安市太白南路2号