发明名称 基于博弈论的正交频分多址接入中继系统资源分配方法
摘要 基于博弈论的正交频分多址接入中继系统资源分配的解决方法,设计了一种综合考虑中继节点选择、子载波配对和功率控制的资源分配方法。中继节点将计算得到的等效信道增益划分为量化区间,并反馈所属量化区间序列号,使得经过较少比特的信息交互可确定协同节点的选择。基站先分配第二跳时的各中继节点的子载波,各中继节点再选出用户节点。以“强强联合”为宗旨匹配两跳的子载波,利用最大流最小割原理,确定中继节点转发时在相应子载波上的发送功率。特别地,建立基于单位功率上的速率的效用函数,并引入基于链路质量与发射功率的代价机制,且代价因子可根据不同业务需求进行实时调整,实现用户节点在各被选中的子载波上的发送功率的优化配置。
申请公布号 CN101730109A 申请公布日期 2010.06.09
申请号 CN200910185438.0 申请日期 2009.11.09
申请人 中国人民解放军理工大学;东南大学 发明人 蔡跃明;吴丹;徐友云
分类号 H04W16/10(2009.01)I;H04W72/08(2009.01)I 主分类号 H04W16/10(2009.01)I
代理机构 南京经纬专利商标代理有限公司 32200 代理人 黄雪兰
主权项 1.一种基于博弈论的正交频分多址接入中继系统资源分配方法,其特征在于:第一步:在单小区OFDMA中继系统中,包括拥有N个子载波的基站d,在基站d的周围均匀分布有M<sub>r</sub>个固定中继节点<img file="F2009101854380C0000011.GIF" wi="413" he="55" />构成固定中继节点集<img file="F2009101854380C0000012.GIF" wi="568" he="91" />其中,M<sub>r</sub>为正整数,用以表示系统中的固定中继节点数,j为[1,M<sub>r</sub>]范围内任一正整数,且每个固定中继节点r<sub>j</sub>拥有N个子载波,单小区OFDMA中继系统还包括随机分布的M<sub>s</sub>个用户节点<img file="F2009101854380C0000013.GIF" wi="423" he="56" />构成用户节点集<img file="F2009101854380C0000014.GIF" wi="580" he="92" />其中,M<sub>s</sub>为正整数,用以表示系统所拥有的用户节点数,i为[1,M<sub>s</sub>]范围内任一正整数,且每个用户节点s<sub>i</sub>拥有N个子载波,每个用户节点s<sub>i</sub>的最大发射功率为<img file="F2009101854380C0000015.GIF" wi="156" he="53" />N为正整数;第二步:用户节点s<sub>i</sub>向所有固定中继节点r<sub>j</sub>和基站d广播信息,各固定中继节点和基站分别对其进行译码,并判断是否译码正确,当基站d译码正确时,基站d向用户节点s<sub>i</sub>和各固定中继节点r<sub>j</sub>发送信息“1”,用户节点s<sub>i</sub>在接收到信息“1”后,用户节点s<sub>i</sub>将确定采用直接传输方式传输,采用直接传输方式的用户节点s<sub>i</sub>在单小区中形成采用直接传输方式传输的用户节点集<img file="F2009101854380C0000016.GIF" wi="81" he="47" />当基站译码不正确时,基站d向用户节点s<sub>i</sub>和各固定中继节点发送信息“0”,当各固定中继节点接收到信息“0”后,能正确译码的固定中继节点r<sub>j</sub>计算自身的等效信道增益<img file="F2009101854380C0000017.GIF" wi="595" he="243" />其中,<img file="F2009101854380C0000018.GIF" wi="82" he="88" />和<img file="F2009101854380C0000019.GIF" wi="129" he="110" />分别为用户节点s<sub>i</sub>与可正确译码的固定中继节点r<sub>j</sub>之间在第n个子载波上的信道状态信息和信道增益,<img file="F2009101854380C00000110.GIF" wi="82" he="89" />和<img file="F2009101854380C00000111.GIF" wi="124" he="110" />分别为可正确译码的固定中继节点r<sub>j</sub>与基站d之间在第n个子载波上的信道状态信息和信道增益,其中,n为子载波序列号,且为[1,N]范围内任一正整数,为减小反馈量,各固定中继节点预先统一设定好等效信道增益的量化区间和量化区间序列号,量化区间预先统一设定的方法是选取等效信道增益范围,这一范围的上界为0,下界为最优等效信道增益<img file="F2009101854380C0000021.GIF" wi="559" he="153" />为利用已有的信道统计模型估算出固定中继节点与基站间的理想信道增益,再将这一范围等间隔划分为m份,m为正整数,从而形成m个量化区间,且将量化区间按照从小到大的顺序依次用自然数1,2,...,m作为序列号对量化区间进行标注,然后,可正确译码的固定中继节点r<sub>j</sub>根据自身的等效信道增益与量化区间的对应关系,确定可正确译码的固定中继节点r<sub>j</sub>的等效信道增益所属的量化区间,接着将对应的量化区间序列号作为反馈值反馈给用户节点s<sub>i</sub>,每个用户节点s<sub>i</sub>选取反馈值最大的固定中继节点r<sub>j</sub>作为帮助其实现协同传输的协同节点,则<img file="F2009101854380C0000022.GIF" wi="50" he="45" />中任一固定中继节点r<sub>j</sub>会同时为若干个用户节点实现协同传输,这些用户节点构成由共享固定中继节点r<sub>j</sub>来实现协同传输的用户节点集<img file="F2009101854380C0000023.GIF" wi="79" he="71" />第三步:各固定中继节点r<sub>j</sub>向基站d发送信息,基站d从中获知其自身与各固定中继节点在各子载波上的信道状态信息,经过比较,由基站d选出对各子载波而言最好的固定中继节点,即<img file="F2009101854380C0000024.GIF" wi="467" he="122" />n=1,2,...,N,基站d将各固定中继节点被分配的子载波数<img file="F2009101854380C0000025.GIF" wi="393" he="92" />范围内的正整数,以及被分配的子载波的序列号n告知各固定中继节点,各固定中继节点以此对其进行标记,即若固定中继节点<img file="F2009101854380C0000026.GIF" wi="273" he="85" />的第n个子载波被选中,则记为<img file="F2009101854380C0000027.GIF" wi="196" he="87" />若未被选中则记为<img file="F2009101854380C0000028.GIF" wi="209" he="89" />第四步:属于由共享固定中继节点r<sub>j</sub>来实现协同传输的用户节点集<img file="F2009101854380C0000029.GIF" wi="53" he="73" />中的用户节点,向对应的作为协同节点的固定中继节点r<sub>j</sub>发送确认信息,作为其协同节点的固定中继节点r<sub>j</sub>从接收的信息中提取与<img file="F2009101854380C00000210.GIF" wi="53" he="72" />中的各用户节点在各子载波上的信道状态信息,经过比较,作为协同节点的固定中继节点r<sub>j</sub>选出对各子载波而言能在该子载波上受到最好帮助的由共享固定中继节点r<sub>j</sub>来实现协同传输的用户节点,即<img file="F2009101854380C0000031.GIF" wi="797" he="128" />并依据已获知的由基站d分配给每个固定中继节点r<sub>j</sub>的子载波数<img file="F2009101854380C0000032.GIF" wi="95" he="70" />在信道增益从大至小的排列下,保留处于前<img file="F2009101854380C0000033.GIF" wi="67" he="69" />位的子载波分配情况,然后,固定中继节点r<sub>j</sub>将各用户节点被分配的子载波数<img file="F2009101854380C0000034.GIF" wi="199" he="66" />为<img file="F2009101854380C0000035.GIF" wi="166" he="99" />范围内的正整数,以及被分配的子载波的序列号n告知<img file="F2009101854380C0000036.GIF" wi="52" he="72" />中的各用户节点,各用户节点以此对其进行标记,即若用户节点<img file="F2009101854380C0000037.GIF" wi="285" he="98" />的第n个子载波被选中,则记为<img file="F2009101854380C0000038.GIF" wi="197" he="88" />若未被选中则记为<img file="F2009101854380C0000039.GIF" wi="211" he="89" />同时,属于采用直接传输方式传输的用户节点集<img file="F2009101854380C00000310.GIF" wi="52" he="42" />中的用户节点,向基站d发送确认信息,基站d从接收的确认信息中提取各用户节点与基站d在各子载波上的信道状态信息,确定用户节点的子载波分配情况,即<img file="F2009101854380C00000311.GIF" wi="465" he="118" />n=1,2,...,N,然后,基站d将各用户节点被分配的子载波数<img file="F2009101854380C00000312.GIF" wi="65" he="66" />以及具体的子载波的序列号n告知采用直接传输方式传输的用户节点集<img file="F2009101854380C00000313.GIF" wi="52" he="43" />中的各用户节点,各用户节点以此对其进行标注,即若用户节点<img file="F2009101854380C00000314.GIF" wi="267" he="71" />的第n个子载波被选中,则记为<img file="F2009101854380C00000315.GIF" wi="193" he="84" />若未被选中则记为<img file="F2009101854380C00000316.GIF" wi="202" he="83" />第五步:建立非合作博弈模型来确定各用户节点被选中的各子载波上的发送功率,每个用户节点以其单位功率上的速率为效用函数,并在效用函数中加入代价机制,以抑制用户节点盲目地增加发送功率,若属于由共享固定中继节点r<sub>j</sub>来实现协同传输的用户节点集<img file="F2009101854380C00000317.GIF" wi="54" he="72" />中的用户节点s<sub>i</sub>在子载波n上的发送功率为<img file="F2009101854380C00000318.GIF" wi="108" he="84" />且在各接收端都有相同的噪声方差σ<sup>2</sup>,则该用户节点s<sub>i</sub>的效用函数可表示为:<img file="F2009101854380C00000319.GIF" wi="800" he="164" />其中,<img file="F2009101854380C0000041.GIF" wi="445" he="151" />为基于用户节点发送功率和与固定中继节点间的信道增益的代价机制,<img file="F2009101854380C0000042.GIF" wi="45" he="50" />为该用户节点s<sub>i</sub>的代价因子,一般而言即为<img file="F2009101854380C0000043.GIF" wi="225" he="75" />具体的发送功率的迭代更新过程如下:(1)该用户节点s<sub>i</sub>初始化被选中的各子载波上的发送功率;(2)在第k次迭代开始时,设置剩余功率值<img file="F2009101854380C0000044.GIF" wi="278" he="69" />发送功率约束范围为<img file="F2009101854380C0000045.GIF" wi="189" he="89" />k=1,2,3,...;(3)该用户节点s<sub>i</sub>按照子载波的序列号依次计算当前迭代中各被选中的子载波上应分配的功率值,具体而言,第n个子载波上分配的功率值应该满足<img file="F2009101854380C0000046.GIF" wi="224" he="178" />即对该用户节点s<sub>i</sub>的效用函数<img file="F2009101854380C0000047.GIF" wi="70" he="84" />求关于<img file="F2009101854380C0000048.GIF" wi="82" he="84" />的导函数,并令其为零,从而得到当前迭代中第n个子载波上分配的功率值<img file="F2009101854380C0000049.GIF" wi="109" he="84" />(4)判断上步骤中获得的当前迭代中第n个子载波上分配的功率值是否满足功率约束条件,当<img file="F2009101854380C00000410.GIF" wi="308" he="91" />时,则子载波n以被分配的功率值作为更新后的发送功率值,否则,第n个子载波以当前剩余功率<img file="F2009101854380C00000411.GIF" wi="56" he="70" />作为更新后的发送功率值,且此时第n个子载波以后的子载波以0值作为迭代更新后的发送功率值,一起进入新一次的迭代;(5)更新剩余功率值<img file="F2009101854380C00000412.GIF" wi="351" he="84" />且更新功率约束范围,更新后的功率约束范围上界为0,下界为更新后的剩余功率值;(6)判断此次迭代第n个子载波上功率更新值与上一次迭代获得的发送功率值之差的二范数是否小于ε,其中ε为预先设定的较小数,它的取值依据系统对收敛速度和精度的要求而定,一般地,ε≤10<sup>-2</sup>,若满足二范数小于ε,则跳出迭代循环,此时所获得的发射功率为博弈算法的纳什均衡,该用户节点s<sub>i</sub>在第n个子载波上最终将以收敛时候的功率值<img file="F2009101854380C0000051.GIF" wi="82" he="83" />进行发射,若不满足,再进入下一次迭代更新过程;同样地,对于采用直接传输方式传输的用户节点集<img file="F2009101854380C0000052.GIF" wi="52" he="41" />中的用户节点s<sub>i</sub>,也以其单位功率上的速率为效用函数,并加入代价机制,具体可表示为:<img file="F2009101854380C0000053.GIF" wi="800" he="129" />其中,<img file="F2009101854380C0000054.GIF" wi="442" he="140" />为基于用户节点发送功率和与基站d之间的信道增益的代价机制,<img file="F2009101854380C0000055.GIF" wi="45" he="52" />为该用户节点s<sub>i</sub>的代价因子,一般而言即为<img file="F2009101854380C0000056.GIF" wi="227" he="75" />具体的发送功率的迭代更新过程如下:(1)该用户节点s<sub>i</sub>初始化被选中的各子载波上的发送功率;(2)在第k次迭代开始时,设置剩余功率值<img file="F2009101854380C0000057.GIF" wi="279" he="69" />发送功率约束范围为<img file="F2009101854380C0000058.GIF" wi="193" he="90" />k=1,2,3,...;(3)该用户节点s<sub>i</sub>按照子载波的序列号依次计算当前迭代中各被选中的子载波上应分配的功率值,具体而言,第n个子载波上分配的功率值应该满足<img file="F2009101854380C0000059.GIF" wi="229" he="179" />即对该用户节点s<sub>i</sub>的效用函数<img file="F2009101854380C00000510.GIF" wi="74" he="84" />求关于<img file="F2009101854380C00000511.GIF" wi="82" he="84" />的导函数,并令其为零,从而得到当前迭代中第n个子载波上分配的功率值<img file="F2009101854380C00000512.GIF" wi="109" he="85" />(4)判断上步骤中获得的当前迭代中第n个子载波上分配的功率值是否满足功率约束条件,当<img file="F2009101854380C00000513.GIF" wi="304" he="91" />时,则子载波n以被分配的功率值作为更新后的发送功率值,否则,第n个子载波以当前剩余功率<img file="F2009101854380C00000514.GIF" wi="56" he="69" />作为更新后的发送功率值,且此时第n个子载波以后的子载波以0值作为迭代更新后的发送功率值,一起进入新一次的迭代;(5)更新剩余功率值<img file="F2009101854380C00000515.GIF" wi="348" he="85" />且更新功率约束范围,更新后的功率约束范围上界为0,下界为更新后的剩余功率值;(6)判断此次迭代第n个子载波上功率更新值与上一次迭代获得的发送功率值之差的二范数是否小于ε,其中ε为预先设定的较小数,它的取值依据系统对收敛速度和精度的要求而定,一般地,ε≤10<sup>-2</sup>,若满足二范数小于ε,则跳出迭代循环,此时所获得的发射功率为博弈算法的纳什均衡,该用户节点s<sub>i</sub>在第n个子载波上最终将以收敛时候的功率值<img file="F2009101854380C0000061.GIF" wi="82" he="84" />进行发射,若不满足,再进入下一次迭代更新过程;第六步:对于需为<img file="F2009101854380C0000062.GIF" wi="53" he="72" />中用户节点转发的固定中继节点r<sub>j</sub>,结合先前已获知的将分配到的子载波序列号,完成转发前的子载波配对,以提高系统性能,所述子载波配对方法是:一方面,按信道增益<img file="F2009101854380C0000063.GIF" wi="128" he="110" />由大到小顺序,对固定中继节点r<sub>j</sub>中收到来自<img file="F2009101854380C0000064.GIF" wi="52" he="72" />中所有用户节点<img file="F2009101854380C0000065.GIF" wi="283" he="98" />信息的子载波进行重新排列,并对重新排列后的子载波的序列号进行重新编制,依次标记为<img file="F2009101854380C0000066.GIF" wi="470" he="75" />且重新排序后的子载波上的接收信噪比仍保持未重新排序时的该子载波上的接收信噪比,例如,若第n个子载波经重新排序后标记为n′,则重新排列后的第n′个子载波上的接收信噪比<img file="F2009101854380C0000067.GIF" wi="446" he="100" />另一方面,按信道增益<img file="F2009101854380C0000068.GIF" wi="125" he="126" />由大到小顺序,对固定中继节点r<sub>j</sub>向基站d发送信息的子载波进行重新排列,并对重新排列后的子载波的序列号进行重新编制,依次标记为<img file="F2009101854380C0000069.GIF" wi="503" he="76" />最后,将两次排列好的固定中继节点r<sub>j</sub>的子载波一一配对,即1′与1″配对、2′与2″配对,…,n′与n″配对,…,<img file="F2009101854380C00000610.GIF" wi="68" he="74" />与<img file="F2009101854380C00000611.GIF" wi="66" he="73" />配对,使得固定中继节点r<sub>j</sub>中收到来自<img file="F2009101854380C00000612.GIF" wi="54" he="72" />中用户节点信息的子载波上的信息将由与之配对的子载波转发至基站d对应的子载波上;此后,确定转发时所需的发送功率,依据最大流最小割原理,则有:<maths num="0001"><![CDATA[<math><mrow><mi>log</mi><mrow><mo>(</mo><mn>1</mn><mo>+</mo><msubsup><mi>&gamma;</mi><msub><mi>r</mi><mi>j</mi></msub><mrow><mo>(</mo><msup><mi>n</mi><mo>&prime;</mo></msup><mo>)</mo></mrow></msubsup><mo>)</mo></mrow><mo>=</mo><mi>log</mi><mrow><mo>(</mo><mn>1</mn><mo>+</mo><mrow><mo>(</mo><msubsup><mi>p</mi><msub><mi>r</mi><mi>j</mi></msub><mrow><mo>(</mo><msup><mi>n</mi><mrow><mo>&prime;</mo><mo>&prime;</mo></mrow></msup><mo>)</mo></mrow></msubsup><msubsup><mi>h</mi><mrow><msub><mi>r</mi><mi>j</mi></msub><mo>,</mo><mi>d</mi></mrow><mrow><mo>(</mo><msup><mi>n</mi><mrow><mo>&prime;</mo><mo>&prime;</mo></mrow></msup><mo>)</mo></mrow></msubsup><mo>)</mo></mrow><mo>/</mo><msup><mi>&sigma;</mi><mn>2</mn></msup><mo>)</mo></mrow></mrow></math>]]></maths>从而确定固定中继节点转发时在重编序列号后的子载波n″上的发送功率为:<maths num="0002"><![CDATA[<math><mrow><msubsup><mi>p</mi><msub><mi>r</mi><mi>j</mi></msub><mrow><mo>(</mo><msup><mi>n</mi><mrow><mo>&prime;</mo><mo>&prime;</mo></mrow></msup><mo>)</mo></mrow></msubsup><mo>=</mo><mrow><mo>(</mo><msubsup><mi>&gamma;</mi><mrow><msub><mi>s</mi><mi>i</mi></msub><mo>,</mo><msub><mi>r</mi><mi>j</mi></msub></mrow><mrow><mo>(</mo><msup><mi>n</mi><mo>&prime;</mo></msup><mo>)</mo></mrow></msubsup><msup><mi>&sigma;</mi><mn>2</mn></msup><mo>)</mo></mrow><mo>/</mo><msubsup><mi>h</mi><mrow><msub><mi>r</mi><mi>j</mi></msub><mo>,</mo><mi>d</mi></mrow><mrow><mo>(</mo><msup><mi>n</mi><mrow><mo>&prime;</mo><mo>&prime;</mo></mrow></msup><mo>)</mo></mrow></msubsup><mo>.</mo></mrow></math>]]></maths>
地址 210007 江苏省南京市白下区御道街标营2号