发明名称 移动场景中的基于时间估计频谱感知系统用户容量的方法
摘要 移动场景中的基于时间估计频谱感知系统用户容量的方法,涉及信息与通信技术领域。本发明是为了获得频谱感知系统中感知用户的用户容量的问题。本发明将频谱感知分为两个阶段,引入了感知用户容量的概念,首先判断在网络覆盖区域,感知用户数是否达到了容限值。如果没有,则下一步进行正常的能量检测工作,研究感知用户和主用户及其他的感知用户之间是否有干扰,如果没有,则可以使用频谱。由于引入了感知用户容量,提高了频谱检测的效率。本发明适用于移动场景中的基于时间估计频谱感知系统用户容量。
申请公布号 CN104796898A 申请公布日期 2015.07.22
申请号 CN201510114700.8 申请日期 2015.03.16
申请人 哈尔滨工业大学 发明人 贾敏;张光宇;顾学迈;郭庆;王欣玉;王世龙;王雪;王振永
分类号 H04W16/10(2009.01)I;H04B17/382(2015.01)I 主分类号 H04W16/10(2009.01)I
代理机构 哈尔滨市松花江专利商标事务所 23109 代理人 张宏威
主权项 移动场景中的基于时间估计频谱感知系统用户容量的方法,其特征是:它由以下步骤实现:步骤一、在网络域A中,设主用户在匀速直线运动长度L=l<sub>i</sub>时,进入感知用户(CR)的干扰域;则:对于k个匀速运动的时间段主用户没有进入干扰域的概率为:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>P</mi><mrow><mo>(</mo><mi>M</mi><mo>></mo><mi>k</mi><mo>|</mo><msub><mi>x</mi><mi>CR</mi></msub><mo>)</mo></mrow><mo>=</mo><munder><mo>&Integral;</mo><msub><mi>L</mi><mn>1</mn></msub></munder><mo>.</mo><mo>.</mo><mo>.</mo><munder><mo>&Integral;</mo><msub><mi>L</mi><mi>k</mi></msub></munder><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msub><mi>P</mi><msub><mi>C</mi><mrow><mo>(</mo><msub><mi>x</mi><mi>CR</mi></msub><mo>,</mo><msub><mi>l</mi><mn>1</mn></msub><mo>)</mo></mrow></msub></msub><mo>)</mo></mrow><mo>.</mo><mo>.</mo><mo>.</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msub><mi>P</mi><msub><mi>C</mi><mrow><mo>(</mo><msub><mi>x</mi><mi>CR</mi></msub><mo>,</mo><msub><mi>l</mi><mi>k</mi></msub><mo>)</mo></mrow></msub></msub><mo>)</mo></mrow><msub><mi>dl</mi><mn>1</mn></msub><mo>.</mo><mo>.</mo><mo>.</mo><msub><mi>d</mi><msub><mi>l</mi><mi>k</mi></msub></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000682728470000011.GIF" wi="1671" he="181" /></maths>式中:M表示从开始到进入干扰域的匀速直线运动的时间段数目,<img file="FDA0000682728470000012.GIF" wi="189" he="121" />是第i次运动主用户进入感知用户干扰域的概率;i为正整数,且i小于k;k为正整数;x<sub>CR</sub>表示感知用户的位置;步骤二、由于最后的匀速运动使到达感知用户(CR)的干扰域的位置是均匀分布的,因此得到:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>P</mi><mrow><mo>(</mo><mi>M</mi><mo>></mo><mi>k</mi><mo>|</mo><msub><mi>x</mi><mi>CR</mi></msub><mo>)</mo></mrow><mo>=</mo><msup><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msub><mi>P</mi><mi>g</mi></msub><mrow><mo>(</mo><msub><mi>x</mi><mi>CR</mi></msub><mo>)</mo></mrow><mo>)</mo></mrow><mi>k</mi></msup><mo>=</mo><msup><mrow><mo>(</mo><mn>1</mn><mo>-</mo><munder><mo>&Integral;</mo><mi>L</mi></munder><msub><mi>P</mi><msub><mi>C</mi><mrow><mo>(</mo><msub><mi>x</mi><mi>CR</mi></msub><mo>,</mo><mi>l</mi><mo>)</mo></mrow></msub></msub><msub><mi>f</mi><mi>L</mi></msub><mrow><mo>(</mo><mi>l</mi><mo>)</mo></mrow><mi>dl</mi><mo>)</mo></mrow><mi>k</mi></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000682728470000013.GIF" wi="1690" he="219" /></maths>式中:f<sub>L</sub>(l)是每步运动大小的概率密度函数;步骤三、将关于x<sub>CR</sub>的M的概率分布函数表示为:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>P</mi><mi>g</mi></msub><mrow><mo>(</mo><msub><mi>x</mi><mi>CR</mi></msub><mo>)</mo></mrow><mo>=</mo><msub><mo>&Integral;</mo><mi>L</mi></msub><msub><mi>P</mi><mrow><mi>C</mi><mrow><mo>(</mo><msub><mi>x</mi><mi>CR</mi></msub><mo>,</mo><mi>l</mi><mo>)</mo></mrow></mrow></msub><msub><mi>f</mi><mi>L</mi></msub><mrow><mo>(</mo><mi>l</mi><mo>)</mo></mrow><mi>dl</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000682728470000014.GIF" wi="1137" he="149" /></maths>步骤四、则P<sub>g</sub>(x<sub>CR</sub>)相当于是每一步主用户成功进入干扰域的概率,而E[M|x<sub>CR</sub>]相当于是M的平均值,即:平均经过E[M|x<sub>CR</sub>]个运动步骤主用户进入任意个CR用户的干扰域,于是有:E[M|x<sub>CR</sub>]=(P<sub>g</sub>(x<sub>CR</sub>))<sup>‑1</sup>    (4)步骤五、每一步主用户成功进入干扰域的概率P<sub>g</sub>(x<sub>CR</sub>)则表示为:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mfenced open='' close=''><mtable><mtr><mtd><msub><mi>P</mi><mi>g</mi></msub><mrow><mo>(</mo><msub><mi>x</mi><mi>CR</mi></msub><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><munderover><mo>&Integral;</mo><mi>R</mi><mrow><mi>R</mi><mo>+</mo><mi>L</mi></mrow></munderover><mn>2</mn><mi>&pi;x</mi><mo>&CenterDot;</mo><mfrac><mrow><mn>2</mn><mi>arccos</mi><mrow><mo>(</mo><mfrac><mrow><msup><mi>x</mi><mn>2</mn></msup><mo>+</mo><msup><mi>L</mi><mn>2</mn></msup><mo>-</mo><msup><mi>R</mi><mn>2</mn></msup></mrow><mrow><mn>2</mn><mi>xL</mi></mrow></mfrac><mo>)</mo></mrow></mrow><mrow><mn>2</mn><mi>&pi;</mi></mrow></mfrac><mi>dx</mi></mrow><msup><mi>a</mi><mn>2</mn></msup></mfrac></mtd></mtr><mtr><mtd><mo>=</mo><mfrac><mrow><munderover><mo>&Integral;</mo><mi>R</mi><mrow><mi>R</mi><mo>+</mo><mi>L</mi></mrow></munderover><mn>2</mn><mi>x</mi><mo>&CenterDot;</mo><mi>arccos</mi><mrow><mo>(</mo><mfrac><mrow><msup><mi>x</mi><mn>2</mn></msup><mo>+</mo><msup><mi>L</mi><mn>2</mn></msup><mo>-</mo><msup><mi>R</mi><mn>2</mn></msup></mrow><mrow><mn>2</mn><mi>xL</mi></mrow></mfrac><mo>)</mo></mrow><mi>dx</mi></mrow><msup><mi>a</mi><mn>2</mn></msup></mfrac></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000682728470000021.GIF" wi="1318" he="514" /></maths>其中:L表示主用户此次运动的每段匀速直线运动的平均距离;R表示主用户的保护域半径;步骤六、判断感知的用户是否为静止的,如果判断结果为是,则执行步骤六一;如果判断结果为否,则执行步骤八;步骤六一、由于主用户在每个感知用户干扰域外的时间需要大于自己通信所需的最小时间,因此得:<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><mfrac><mi>t</mi><mrow><mi>K</mi><mo>&CenterDot;</mo><msub><mi>P</mi><mi>g</mi></msub><mrow><mo>(</mo><msub><mi>x</mi><mi>CR</mi></msub><mo>)</mo></mrow></mrow></mfrac><mo>&GreaterEqual;</mo><mi>T</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000682728470000022.GIF" wi="1040" he="143" /></maths>式中:t为主用户做每一段匀速直线运动的时间,T为保证主用户通信所需要的最短的时间间隔,K为感知用户的个数;步骤七、根据步骤六一算出的结果,则得到感知用户的容量为:<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><msub><mi>K</mi><mi>max</mi></msub><mo>=</mo><mfrac><mi>t</mi><mrow><mi>T</mi><mo>&CenterDot;</mo><msub><mi>P</mi><mi>g</mi></msub><mrow><mo>(</mo><msub><mi>x</mi><mi>CR</mi></msub><mo>)</mo></mrow></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0000682728470000023.GIF" wi="1108" he="141" /></maths>步骤八、设感知同一频谱的感知用户当前的状态是有K<sub>1</sub>个用户在运动,K<sub>2</sub>个感知用户处于静止状态,并且正在做运动的感知用户的运动也可抽象成几段的匀速直线运动,主用户在任何时刻选择方向,角度上都服从在[0,2π]之间的均匀分布的,则把运动的感知用户,抽象成是使用同一频率的主用户,得到:<maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><mfrac><msup><mi>t</mi><mo>&prime;</mo></msup><mrow><mrow><mo>(</mo><msub><mi>K</mi><mn>1</mn></msub><mo>+</mo><msub><mi>K</mi><mn>2</mn></msub><mo>)</mo></mrow><mo>&CenterDot;</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msub><mi>P</mi><mi>g</mi></msub><mrow><mo>(</mo><msub><mi>x</mi><mi>CR</mi></msub><mo>)</mo></mrow><mo>)</mo></mrow><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msub><mi>K</mi><mn>1</mn></msub><msub><mi>P</mi><mi>g</mi></msub><mrow><mo>(</mo><msub><mi>x</mi><msub><mi>CR</mi><msub><mi>K</mi><mn>2</mn></msub></msub></msub><mo>)</mo></mrow><mo>)</mo></mrow><mo>)</mo></mrow></mrow></mfrac><mo>&GreaterEqual;</mo><mi>T</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>8</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000682728470000024.GIF" wi="1608" he="189" /></maths>其中:t′表示保证主用户和运动中感知用户通信的最短的时间,<img file="FDA0000682728470000025.GIF" wi="255" he="110" />表示运动中的感知用户每运动一步进入静止的感知用户和主用户干扰域的概率;步骤九、判断运动中的感知用户数是否已知,如果判断结果为是,则执行步骤九一;如果判断结果为否,则执行步骤十;步骤九一、根据公式:<maths num="0008" id="cmaths0008"><math><![CDATA[<mrow><msub><mi>K</mi><mi>max</mi></msub><mo>=</mo><mfrac><mi>t</mi><mrow><mi>T</mi><mo>&CenterDot;</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msub><mi>P</mi><mi>g</mi></msub><mrow><mo>(</mo><msub><mi>x</mi><mi>CR</mi></msub><mo>)</mo></mrow><mo>)</mo></mrow><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msub><mi>K</mi><mn>1</mn></msub><msub><mi>P</mi><mi>g</mi></msub><mrow><mo>(</mo><msub><mi>x</mi><msub><mi>CR</mi><msub><mi>K</mi><mn>2</mn></msub></msub></msub><mo>)</mo></mrow><mo>)</mo></mrow><mo>)</mo></mrow></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>9</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000682728470000031.GIF" wi="1351" he="186" /></maths>获得容纳感知用户的容量K<sub>max</sub>;步骤十、根据公式(9)作为已知感知用户的容量是一个随K<sub>1</sub>单调递减的函数,且取K<sub>2</sub>=0,则根据公式:<maths num="0009" id="cmaths0009"><math><![CDATA[<mrow><msub><mi>K</mi><mi>max</mi></msub><mo>=</mo><msub><mi>K</mi><mn>1</mn></msub><mo>=</mo><mfrac><msup><mi>t</mi><mo>&prime;</mo></msup><mi>T</mi></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>10</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000682728470000032.GIF" wi="1130" he="132" /></maths>获得容纳感知用户的容量K<sub>max</sub>。
地址 150001 黑龙江省哈尔滨市南岗区西大直街92号