主权项 |
一种基于有限随机矩阵的频谱检测方法,包括如下步骤:(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>·</mo></mtd><mtd><mo>·</mo></mtd><mtd><mo>·</mo></mtd><mtd><mo>·</mo></mtd></mtr><mtr><mtd><mo>·</mo></mtd><mtd><mo>·</mo></mtd><mtd><mo>·</mo></mtd><mtd><mo>·</mo></mtd></mtr><mtr><mtd><mo>·</mo></mtd><mtd><mo>·</mo></mtd><mtd><mo>·</mo></mtd><mtd><mo>·</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>·</mo></mtd><mtd><mo>·</mo></mtd><mtd><mo>·</mo></mtd><mtd><mo>·</mo></mtd></mtr><mtr><mtd><mo>·</mo></mtd><mtd><mo>·</mo></mtd><mtd><mo>·</mo></mtd><mtd><mo>·</mo></mtd></mtr><mtr><mtd><mo>·</mo></mtd><mtd><mo>·</mo></mtd><mtd><mo>·</mo></mtd><mtd><mo>·</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>ξ</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>γ</mi><mi>ξ</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>时,判决为主用户存在,即当前频段频谱已被某用户占用,否则,判决为主用户不存在,即当前频段频谱为空闲状态,允许检测用户利用。 |