发明名称 一种认知无线电系统中基于效用函数的分布式功率控制方法
摘要 一种认知无线电(Cognitive Radio)系统中基于效用函数的分布式功率控制方法,属于无线通信技术领域,特别涉及应用在认知无线电系统中的功率控制。本发明是一种既满足了对主用户的保护又支持了认知用户的服务质量要求的折中方案。此方案基于非合作式的博弈论,即,它只需要各个CRU自行调节其发射功率,因此避免了传统分布式功率控制需要大量信息交互的弊端。通过分析效用函数的几何特性,巧妙地设置了效用函数的参数,同时结合接入控制,保证了功率的收敛。同时,通过调整参数,能够适应现有网络业务的需要。此外,我们设计了与功率控制方法相适应的帧结构和协议流程,它们可以适用于任意分布式认知无线电系统下,因而具有较强的实用性。
申请公布号 CN101626260A 申请公布日期 2010.01.13
申请号 CN200810045492.0 申请日期 2008.07.08
申请人 电子科技大学 发明人 车越岭;陈劼;王军;唐万斌
分类号 H04B7/005(2006.01)I;H04B1/69(2006.01)I 主分类号 H04B7/005(2006.01)I
代理机构 代理人
主权项 1、一种认知无线电(Cognitive Radio)系统中基于效用函数的分布式功率控制方法,其特征在于:认知无线电功率控制的两个最基本的要求是实现对PUs的保护(保证CRUs对PUs的干扰不会超过PUs能承受的干扰范围)和提供对CRUs的QoS支持。一般来讲,QoS的好坏与SINR的满足与否有较大的关系。本方案的增益函数G(i)表示CRUi的SINR要求——CRUi接收端获得的信干噪比(SINR)γ<sub>i</sub>要达到QoS要求的信干噪比γ<sub>d</sub>:<maths num="0001"><![CDATA[<math><mrow><msub><mi>&gamma;</mi><mi>i</mi></msub><mo>=</mo><mfrac><mrow><msub><mi>p</mi><mi>i</mi></msub><mo>&CenterDot;</mo><msub><mi>h</mi><mi>ii</mi></msub></mrow><mrow><msubsup><mi>&Sigma;</mi><munder><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mrow><mi>j</mi><mo>&NotEqual;</mo><mi>i</mi></mrow></munder><mi>N</mi></msubsup><msub><mi>p</mi><mi>j</mi></msub><mo>&CenterDot;</mo><msub><mi>h</mi><mi>ji</mi></msub><mo>+</mo><msub><mi>n</mi><mi>o</mi></msub></mrow></mfrac><mo>&GreaterEqual;</mo><msub><mi>&gamma;</mi><mi>d</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow></math>]]></maths>此外,利用代价函数C(i)表示对PUs的保护——本发明假设在进行功率控制之前,我们已经得到了PUs能承受的最小干扰,P<sub>T</sub>,则对PUs的保护可以表示为:<maths num="0002"><![CDATA[<math><mrow><msubsup><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></msubsup><msub><mi>p</mi><mi>i</mi></msub><mo>&CenterDot;</mo><msub><mi>g</mi><mi>i</mi></msub><mo>&le;</mo><msub><mi>P</mi><mi>T</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow></math>]]></maths>这样,每个想接入本系统的CRU都拥有了一个效用函数U(i)=G(i)-C(i)。直接由公式(1)的SINR要求,增益函数表示如下:<maths num="0003"><![CDATA[<math><mrow><mi>G</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>&lambda;</mi><mi>i</mi></msub><msqrt><msub><mi>&gamma;</mi><mi>i</mi></msub><mo>-</mo><msub><mi>&gamma;</mi><mi>d</mi></msub></msqrt><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow></math>]]></maths>上式要求γ<sub>i</sub>≥γ<sub>d</sub>以使函数有意义。其中,λ<sub>i</sub>是一个可调参数。代价函数可表示为:<maths num="0004"><![CDATA[<math><mrow><mi>C</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>=</mo><msup><mrow><mo>(</mo><mfrac><mrow><msub><mi>p</mi><mi>i</mi></msub><mo>&CenterDot;</mo><msub><mi>g</mi><mi>i</mi></msub></mrow><mrow><msub><mi>P</mi><mi>T</mi></msub><mo>-</mo><msubsup><mi>&Sigma;</mi><munder><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mrow><mi>j</mi><mo>&NotEqual;</mo><mi>i</mi></mrow></munder><mi>N</mi></msubsup><msub><mi>p</mi><mi>j</mi></msub><mo>&CenterDot;</mo><msub><mi>g</mi><mi>j</mi></msub></mrow></mfrac><mo>)</mo></mrow><msub><mi>a</mi><mi>i</mi></msub></msup><mo>&CenterDot;</mo><msub><mi>p</mi><mi>i</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow></math>]]></maths>p<sub>i</sub>可以表示为γ<sub>i</sub>的函数,结合公式(3)和公式(4),本发明的效用函数如下:<maths num="0005"><![CDATA[<math><mrow><mi>U</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>&lambda;</mi><mi>i</mi></msub><msqrt><msub><mi>&gamma;</mi><mi>i</mi></msub><mo>-</mo><msub><mi>&gamma;</mi><mi>d</mi></msub></msqrt><mo>-</mo><msup><mrow><mo>(</mo><mfrac><msub><mi>g</mi><mi>i</mi></msub><mrow><msub><mi>P</mi><mi>T</mi></msub><mo>-</mo><msubsup><mi>&Sigma;</mi><munder><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mrow><mi>j</mi><mo>&NotEqual;</mo><mi>i</mi></mrow></munder><mi>N</mi></msubsup><msub><mi>p</mi><mi>j</mi></msub><mo>&CenterDot;</mo><msub><mi>g</mi><mi>j</mi></msub></mrow></mfrac><mo>)</mo></mrow><msub><mi>a</mi><mi>i</mi></msub></msup><mo>&CenterDot;</mo><msup><mrow><mo>(</mo><mfrac><mrow><msub><mi>&gamma;</mi><mi>i</mi></msub><mo>&CenterDot;</mo><msub><mi>Intf</mi><mi>i</mi></msub></mrow><msub><mi>h</mi><mi>ii</mi></msub></mfrac><mo>)</mo></mrow><mrow><msub><mi>a</mi><mi>i</mi></msub><mo>+</mo><mn>1</mn></mrow></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow></math>]]></maths>则功率控制问题就可被详细描述为:寻找到一组功率向量P=[p<sub>1</sub>,p<sub>2</sub>…p<sub>N</sub>],使得对于所有合适的CRUs有下面的问题成立:<maths num="0006"><![CDATA[<math><mrow><mfenced open='{' close=''><mtable><mtr><mtd><munder><mi>max</mi><mrow><msub><mi>p</mi><mi>i</mi></msub><mo>&GreaterEqual;</mo><mn>0</mn></mrow></munder><mi>U</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow></mtd></mtr><mtr><mtd><msub><mi>&gamma;</mi><mi>i</mi></msub><mo>&GreaterEqual;</mo><msub><mi>&gamma;</mi><mi>d</mi></msub></mtd></mtr><mtr><mtd><mi>U</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>&GreaterEqual;</mo><mn>0</mn></mtd></mtr><mtr><mtd><msub><mi>p</mi><mi>i</mi></msub><mo>&CenterDot;</mo><msub><mi>g</mi><mi>i</mi></msub><mo>&le;</mo><msub><mi>P</mi><mi>T</mi></msub><mo>-</mo><msubsup><mi>&Sigma;</mi><munder><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mrow><mi>j</mi><mo>&NotEqual;</mo><mi>i</mi></mrow></munder><mi>N</mi></msubsup><msub><mi>p</mi><mi>j</mi></msub><mo>&CenterDot;</mo><msub><mi>g</mi><mi>j</mi></msub></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow></math>]]></maths>令p<sub>i</sub><sup>*</sup>是对CRUi而言的局部最优值,则:<maths num="0007"><![CDATA[<math><mrow><mfrac><mrow><mo>&PartialD;</mo><mi>U</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow></mrow><msubsup><mrow><mo>&PartialD;</mo><mi>p</mi></mrow><mi>i</mi><mo>*</mo></msubsup></mfrac><mo>=</mo><mn>0</mn><mo>&DoubleRightArrow;</mo><msubsup><mi>&gamma;</mi><mi>i</mi><msub><mi>a</mi><mi>i</mi></msub></msubsup><mo>&CenterDot;</mo><msqrt><msub><mi>&gamma;</mi><mi>i</mi></msub><mo>-</mo><msub><mi>&gamma;</mi><mi>d</mi></msub></msqrt><mo>=</mo><mfrac><msub><mi>&lambda;</mi><mi>i</mi></msub><mn>2</mn></mfrac><mo>&CenterDot;</mo><mfrac><mn>1</mn><mrow><mrow><mo>(</mo><msub><mi>a</mi><mi>i</mi></msub><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>&CenterDot;</mo><msub><mi>&beta;</mi><mi>i</mi></msub></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow></mrow></math>]]></maths>其中,<maths num="0008"><![CDATA[<math><mrow><msub><mi>&beta;</mi><mi>i</mi></msub><mo>=</mo><msup><mrow><mo>(</mo><mfrac><msub><mi>g</mi><mi>i</mi></msub><mrow><msub><mi>P</mi><mi>T</mi></msub><mo>-</mo><msubsup><mi>&Sigma;</mi><munder><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mrow><mi>j</mi><mo>&NotEqual;</mo><mi>i</mi></mrow></munder><mi>N</mi></msubsup><msub><mi>p</mi><mi>j</mi></msub><mo>&CenterDot;</mo><msub><mi>g</mi><mi>j</mi></msub></mrow></mfrac><mo>)</mo></mrow><msub><mi>a</mi><mi>i</mi></msub></msup><mo>&CenterDot;</mo><msup><mrow><mo>(</mo><mfrac><msub><mi>Intf</mi><mi>i</mi></msub><msub><mi>h</mi><mi>ii</mi></msub></mfrac><mo>)</mo></mrow><mrow><msub><mi>a</mi><mi>i</mi></msub><mo>+</mo><mn>1</mn></mrow></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>8</mn><mo>)</mo></mrow></mrow></math>]]></maths>令<maths num="0009"><![CDATA[<math><mrow><mi>f</mi><mrow><mo>(</mo><msub><mi>&gamma;</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>=</mo><msubsup><mi>&gamma;</mi><mi>i</mi><msub><mi>a</mi><mi>i</mi></msub></msubsup><mo>&CenterDot;</mo><msqrt><msub><mi>&gamma;</mi><mi>i</mi></msub><mo>-</mo><msub><mi>&gamma;</mi><mi>d</mi></msub></msqrt><mo>,</mo></mrow></math>]]></maths>则最优解为:<maths num="0010"><![CDATA[<math><mrow><msubsup><mi>&gamma;</mi><mi>i</mi><mo>*</mo></msubsup><mo>=</mo><msup><mi>f</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mrow><mo>(</mo><mfrac><msub><mi>&lambda;</mi><mi>i</mi></msub><mn>2</mn></mfrac><mo>&CenterDot;</mo><mfrac><mn>1</mn><mrow><mrow><mo>(</mo><msub><mi>a</mi><mi>i</mi></msub><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>&CenterDot;</mo><msub><mi>&beta;</mi><mi>i</mi></msub></mrow></mfrac><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>9</mn><mo>)</mo></mrow></mrow></math>]]></maths>相应的功率分配为:<maths num="0011"><![CDATA[<math><mrow><msubsup><mi>p</mi><mi>i</mi><mo>*</mo></msubsup><mo>=</mo><msubsup><mi>&gamma;</mi><mi>i</mi><mo>*</mo></msubsup><mo>&CenterDot;</mo><mfrac><msub><mi>Intf</mi><mi>i</mi></msub><msub><mi>h</mi><mi>ii</mi></msub></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>10</mn><mo>)</mo></mrow></mrow></math>]]></maths>因此,CRUi根据以下原则更新其在第k+1步的发射功率p<sub>i</sub><sup>*</sup>(k+1):<maths num="0012"><![CDATA[<math><mrow><msubsup><mi>p</mi><mi>i</mi><mo>*</mo></msubsup><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>=</mo><msubsup><mi>&gamma;</mi><mi>i</mi><mo>*</mo></msubsup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>&CenterDot;</mo><mfrac><mrow><msub><mi>Intf</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mrow><mrow><msub><mi>h</mi><mi>ii</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mrow></mfrac><mo>=</mo><mfrac><mrow><msubsup><mi>&gamma;</mi><mi>i</mi><mo>*</mo></msubsup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mrow><mrow><msub><mi>&gamma;</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></mrow></mfrac><mo>&CenterDot;</mo><msub><mi>p</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>11</mn><mo>)</mo></mrow></mrow></math>]]></maths>然而,一些研究结果表明公式(11)的形式表明,当存在一个合适的功率分配使得任意接入用户i的SINR满足γ<sub>i</sub>≥γ<sub>d</sub>时,式(11)能保证功率收敛。但若上述条件不成立,公式(11)的方法会导致功率不收敛。因此,本发明结合接入控制来拒绝不受欢迎的CRUs以解决收敛问题。现在的问题是以什么样的标准判断一个CRUi是否是受欢迎的或是允许接入的,换句话说,如何找到一个合适的γ<sub>i</sub><sup>*</sup>以满足式(6)中的所有条件。通过调整参数β<sub>i</sub>,会有三种不同的C(i)的曲线,它们分别和G(i)相离,相切和相交。当<maths num="0013"><![CDATA[<math><mrow><msub><mi>p</mi><mi>i</mi></msub><mo>&CenterDot;</mo><msub><mi>g</mi><mi>i</mi></msub><mo>=</mo><msub><mi>P</mi><mi>T</mi></msub><mo>-</mo><msub><mi>&Sigma;</mi><munder><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mrow><mi>j</mi><mo>&NotEqual;</mo><mi>i</mi></mrow></munder></msub><msub><mi>p</mi><mi>j</mi></msub><mo>&CenterDot;</mo><msub><mi>g</mi><mi>j</mi></msub></mrow></math>]]></maths>时,如果G(i)和C(i)相切,则令相应的<maths num="0014"><![CDATA[<math><mrow><msub><mi>&beta;</mi><mi>i</mi></msub><mo>=</mo><msubsup><mi>&beta;</mi><mrow><mi>i</mi><mo>,</mo></mrow><mi>t</mi></msubsup></mrow></math>]]></maths><maths num="0015"><![CDATA[<math><mrow><msub><mi>&gamma;</mi><mi>i</mi></msub><mo>=</mo><msubsup><mi>&gamma;</mi><mi>i</mi><mi>t</mi></msubsup><mo>.</mo></mrow></math>]]></maths>[性质1]对任一想接入本系统的CRUi而言,若其<maths num="0016"><![CDATA[<math><mrow><msub><mi>&beta;</mi><mi>i</mi></msub><mo>&le;</mo><msubsup><mi>&beta;</mi><mi>i</mi><mi>t</mi></msubsup><mo>,</mo></mrow></math>]]></maths>则认为此认知用户是受欢迎的。[性质2]当<maths num="0017"><![CDATA[<math><mrow><msub><mi>&beta;</mi><mi>i</mi></msub><mo>=</mo><msubsup><mi>&beta;</mi><mi>i</mi><mi>t</mi></msubsup></mrow></math>]]></maths>时,令<maths num="0018"><![CDATA[<math><mrow><msub><mi>&theta;</mi><mi>i</mi></msub><mo>=</mo><mfrac><mrow><msub><mi>&gamma;</mi><mi>d</mi></msub><mo>&CenterDot;</mo><msub><mi>Intf</mi><mi>i</mi></msub><mo>&CenterDot;</mo><msub><mi>g</mi><mi>i</mi></msub></mrow><mrow><msub><mi>remainder</mi><mi>i</mi></msub><mo>&CenterDot;</mo><msub><mi>h</mi><mi>ii</mi></msub></mrow></mfrac><mo>,</mo></mrow></math>]]></maths>有:<maths num="0019"><![CDATA[<math><mrow><msub><mi>a</mi><mi>i</mi></msub><mo>=</mo><mfrac><mn>1</mn><mrow><mn>2</mn><mo>&CenterDot;</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msub><mi>&theta;</mi><mi>i</mi></msub><mo>)</mo></mrow></mrow></mfrac><mo>-</mo><mn>1</mn><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>12</mn><mo>)</mo></mrow></mrow></math>]]></maths>令p<sub>i,ph</sub><sup>max</sup>表示发射机允许的物理最大发射功率,其相应的SINR为γ<sub>i,ph</sub><sup>max</sup>,则γ<sub>i</sub>的最大值γ<sub>i</sub><sup>max</sup>为:<maths num="0020"><![CDATA[<math><mrow><msubsup><mi>&gamma;</mi><mi>i</mi><mi>max</mi></msubsup><mo>=</mo><mi>max</mi><mrow><mo>(</mo><msubsup><mi>&gamma;</mi><mrow><mi>i</mi><mo>,</mo><mi>ph</mi></mrow><mi>max</mi></msubsup><mo>,</mo><msubsup><mi>&gamma;</mi><mi>i</mi><mi>t</mi></msubsup><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>13</mn><mo>)</mo></mrow></mrow></math>]]></maths>有了以上两个性质作为理论基础,下面我们给出本发明中的算法的步骤——基于效用函数的分布式功率控制(UPCAC,Utility based Power Control jointing with Admission Control):步骤1更新Intf<sub>i</sub>(k),g<sub>i</sub>(k),<img file="A2008100454920004C8.GIF" wi="340" he="90" />和h<sub>ii</sub>(k),计算α<sub>i</sub>(k),β<sub>i</sub>(k)和β<sub>i</sub><sup>t</sup>(k)。如果<maths num="0021"><![CDATA[<math><mrow><msub><mi>&beta;</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>&le;</mo><msubsup><mi>&beta;</mi><mi>i</mi><mi>t</mi></msubsup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>,</mo></mrow></math>]]></maths>至步骤2;否则,至步骤4;步骤2令与γ<sub>i</sub><sup>max</sup>(k)相应的功率为p<sub>i</sub><sup>max</sup>(k+1),与γ<sub>d</sub>对应的功率为p<sub>i</sub><sup>min</sup>(k+1)。若<maths num="0022"><![CDATA[<math><mrow><msubsup><mi>p</mi><mi>i</mi><mi>min</mi></msubsup><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>&le;</mo><msubsup><mi>p</mi><mi>i</mi><mi>max</mi></msubsup><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>,</mo></mrow></math>]]></maths>至步骤3;否则,至步骤4;步骤3在闭区间[p<sub>i</sub><sup>min</sup>(k+1),p<sub>i</sub><sup>max</sup>(k+1)]上找到最优点p<sub>i</sub>*(k+1)以最大化U(i),至步骤5;步骤4p<sub>i</sub>*(k+1)=0,U(i)(k+1)=0;步骤5令k←k+1,至步骤1。若一个想接入本网络的CRU的最终的功率非零,则认为此CRU是受欢迎的并且允许接入本网络。由于在UPCAC中,我们保证了0<inf(i)<1,因此我们确保了对主用户的保护。
地址 610054四川省成都市建设北路二段四号通信大楼525