发明名称 一种大规模天线系统中基于遗传算法的导频分配方法
摘要 一种大规模天线系统中基于遗传算法的导频分配方法,以提高系统中用户的和速率为目标,采用遗传算法,得到最优的导频分配方案。首先,引入基因编码的概念,通过交换编码产生不同的导频分配方案,进而计算不同导频分配方案下的系统中用户的和速率,同时采用优胜劣汰的准则,保留对应的基因编码;其次,对保留的基因编码进行复制、交叉和变异操作,产生新的基因编码;最后,重复以上步骤,经过多轮的优胜劣汰和基因的交叉变异,得到使得系统中用户和速率最大的导频分配方案。本发明采用遗传算法,突破局部最优解的限制,从而缓解大规模天线系统的导频污染影响,并降低导频分配策略的复杂度。
申请公布号 CN105024793A 申请公布日期 2015.11.04
申请号 CN201510357846.5 申请日期 2015.06.25
申请人 山东大学 发明人 白智全;张标;孔凡堂;苏英彦;高鹏;孙秀凯
分类号 H04L5/00(2006.01)I;G06N3/12(2006.01)I 主分类号 H04L5/00(2006.01)I
代理机构 济南金迪知识产权代理有限公司 37219 代理人 许德山
主权项 一种大规模天线系统中基于遗传算法的导频分配方法,用于在L个小区中所设置的大规模天线系统,其中每个小区含有一个基站和K个单天线终端用户,基站天线数目为理想情况,即其趋近于无穷大,通信过程采用时分双工机制,并考虑信道间的互易性,同时,考虑有K个正交导频在L个小区中复用,以提高系统中各小区用户和速率为目标,通过采用遗传算法得到最优的导频分配方案,首先,通过引入基因编码的概念,利用交换编码产生不同的导频分配方案,进而计算不同导频分配方案下的系统中的用户和速率,同时采用优胜劣汰的准则,保留对应的基因编码;其次,对保留的基因编码进行复制、交叉和变异操作,产生新的基因编码;最后,重复以上步骤,经过多轮的优胜劣汰和基因的交叉变异,得到使系统中用户和速率最大的导频分配方案,该方法具体步骤如下:1)随机初始化一组导频分配方案为<img file="FDA0000745541160000011.GIF" wi="485" he="108" />且<img file="FDA0000745541160000012.GIF" wi="842" he="127" />将其用于初始化第l个小区导频分配方案;复制导频分配方案<img file="FDA0000745541160000013.GIF" wi="421" he="118" />为任意一组导频分配方案,记为M组,即<img file="FDA0000745541160000014.GIF" wi="838" he="116" />第l个小区的导频分配方案为<img file="FDA0000745541160000015.GIF" wi="138" he="99" />其序列<img file="FDA0000745541160000016.GIF" wi="372" he="102" />由整数1到K的一种有序排列构成,且<img file="FDA0000745541160000017.GIF" wi="105" he="103" />表示第l个小区中分配第k个导频的用户编号,因而<img file="FDA0000745541160000018.GIF" wi="96" he="99" />可以唯一确定第l个小区的导频分配方案;2)随机初始化C组基因编码为<img file="FDA0000745541160000019.GIF" wi="783" he="119" />这里每组基因编码有L条染色体组成,其中第c组基因编码中的第l条染色体为<img file="FDA00007455411600000110.GIF" wi="614" he="103" />其中<img file="FDA00007455411600000111.GIF" wi="592" he="95" />表示染色体<img file="FDA00007455411600000112.GIF" wi="105" he="95" />的第n个基因,N表示该染色体上的基因个数,并且满足以下约束条件,<img file="FDA00007455411600000113.GIF" wi="701" he="107" />这里设置C=M,基因编码与导频分配方案一一对应,具体地,第c组基因编码<img file="FDA00007455411600000114.GIF" wi="422" he="112" />用于第c组导频分配方案<img file="FDA00007455411600000115.GIF" wi="408" he="112" />的遗传;3)初始化遗传代数t=0,并设置最大遗传代数为t<sub>max</sub>;4)分别根据C组基因编码<img file="FDA00007455411600000116.GIF" wi="779" he="112" />对C组导频分配方案<img file="FDA00007455411600000117.GIF" wi="761" he="114" />进行交换编码,其中第c组基因编码中的染色体<img file="FDA00007455411600000118.GIF" wi="94" he="96" />用做第l个小区用户的导频分配方案<img file="FDA0000745541160000021.GIF" wi="95" he="99" />的交换编码,可得到C组新的导频分配方案<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mo>{</mo><msubsup><mi>P</mi><mrow><mi>c</mi><mn>1</mn></mrow><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></msubsup><mo>,</mo><msubsup><mi>P</mi><mrow><mi>c</mi><mn>2</mn></mrow><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></msubsup><mo>,</mo><mn>...</mn><mo>,</mo><msubsup><mi>P</mi><mrow><mi>c</mi><mi>L</mi></mrow><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></msubsup><mo>}</mo><mo>,</mo><mi>c</mi><mo>&Element;</mo><mo>{</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>,</mo><mn>...</mn><mo>,</mo><mi>C</mi><mo>}</mo><mo>,</mo></mrow>]]></math><img file="FDA0000745541160000022.GIF" wi="777" he="113" /></maths>对于<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mo>&ForAll;</mo><mi>c</mi><mo>&Element;</mo><mo>{</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>,</mo><mn>...</mn><mo>,</mo><mi>c</mi><mo>}</mo><mo>,</mo></mrow>]]></math><img file="FDA0000745541160000023.GIF" wi="375" he="76" /></maths>l∈{1,2,...,L},交换编码即是指根据<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msubsup><mi>T</mi><mrow><mi>c</mi><mi>l</mi></mrow><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></msubsup><mo>=</mo><mo>&lsqb;</mo><msubsup><mi>t</mi><mrow><mi>c</mi><mi>l</mi><mn>1</mn></mrow><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></msubsup><mo>,</mo><msubsup><mi>t</mi><mrow><mi>c</mi><mi>l</mi><mn>2</mn></mrow><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></msubsup><mo>,</mo><mn>...</mn><mo>,</mo><msubsup><mi>t</mi><mrow><mi>c</mi><mi>l</mi><mi>N</mi></mrow><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></msubsup><mo>&rsqb;</mo></mrow>]]></math><img file="FDA0000745541160000024.GIF" wi="528" he="103" /></maths>对<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msubsup><mi>p</mi><mrow><mi>c</mi><mi>l</mi></mrow><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></msubsup><mo>=</mo><mo>&lsqb;</mo><msubsup><mi>p</mi><mrow><mi>c</mi><mi>l</mi><mn>1</mn></mrow><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></msubsup><mo>,</mo><msubsup><mi>p</mi><mrow><mi>c</mi><mi>l</mi><mn>2</mn></mrow><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></msubsup><mo>,</mo><mn>...</mn><mo>,</mo><msubsup><mi>p</mi><mrow><mi>c</mi><mi>l</mi><mi>K</mi></mrow><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></msubsup><mo>&rsqb;</mo></mrow>]]></math><img file="FDA0000745541160000025.GIF" wi="547" he="104" /></maths>进行交换编码,其具体操作如下:(1)把<img file="FDA0000745541160000026.GIF" wi="217" he="94" />保存到临时变量p<sub>tmp</sub>中;(2)对于<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><mo>&ForAll;</mo><mi>n</mi><mo>&Element;</mo><mo>{</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>,</mo><mn>...</mn><mo>,</mo><mi>N</mi><mo>-</mo><mn>1</mn><mo>}</mo><mo>,</mo></mrow>]]></math><img file="FDA0000745541160000027.GIF" wi="460" he="82" /></maths>依次<img file="FDA0000745541160000028.GIF" wi="296" he="103" />赋值给<img file="FDA0000745541160000029.GIF" wi="242" he="118" />(3)把p<sub>tmp</sub>赋值给<img file="FDA00007455411600000210.GIF" wi="245" he="130" />这里对于<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><mo>&ForAll;</mo><mi>n</mi><mo>&Element;</mo><mo>{</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>,</mo><mn>...</mn><mo>,</mo><mi>N</mi><mo>}</mo><mo>,</mo></mrow>]]></math><img file="FDA00007455411600000211.GIF" wi="379" he="73" /></maths><img file="FDA00007455411600000212.GIF" wi="224" he="96" />表示<img file="FDA00007455411600000213.GIF" wi="99" he="99" />的第<img file="FDA00007455411600000214.GIF" wi="92" he="102" />个元素值;5)计算C个适应值,即分别采用C组导频分配方案<img file="FDA00007455411600000215.GIF" wi="761" he="111" />时,来得到系统中用户的和速率<img file="FDA00007455411600000216.GIF" wi="544" he="77" />用户的和速率<img file="FDA00007455411600000217.GIF" wi="118" he="76" />的获得,依赖于导频分配方案<img file="FDA00007455411600000218.GIF" wi="442" he="128" />且<img file="FDA00007455411600000219.GIF" wi="124" he="73" />由下式给出:<maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><msubsup><mi>R</mi><mi>c</mi><mrow><mi>s</mi><mi>u</mi><mi>m</mi></mrow></msubsup><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>l</mi><mo>=</mo><mn>1</mn></mrow><mi>L</mi></munderover><munderover><mo>&Sigma;</mo><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>K</mi></munderover><msub><mi>R</mi><mrow><mi>l</mi><mo>,</mo><mi>k</mi></mrow></msub><mo>.</mo></mrow>]]></math><img file="FDA00007455411600000220.GIF" wi="433" he="200" /></maths>其中,R<sub>l,k</sub>表示第l个小区中第k个用户的速率,并由下式给出,<maths num="0008" id="cmaths0008"><math><![CDATA[<mrow><msub><mi>R</mi><mrow><mi>l</mi><mo>,</mo><mi>k</mi></mrow></msub><mo>=</mo><msub><mi>log</mi><mn>2</mn></msub><mrow><mo>(</mo><mn>1</mn><mo>+</mo><mfrac><msubsup><mi>&beta;</mi><mi>llk</mi><mn>2</mn></msubsup><mrow><msub><mi>&Sigma;</mi><mrow><mi>j</mi><mo>&NotEqual;</mo><mn>1</mn></mrow></msub><msubsup><mi>&beta;</mi><mi>jlk</mi><mn>2</mn></msubsup></mrow></mfrac><mo>)</mo></mrow><mo>.</mo></mrow>]]></math><img file="FDA00007455411600000221.GIF" wi="580" he="177" /></maths>其中,β<sub>llk</sub>表示第l个小区中第k个用户到第l个基站的大尺度衰落,β<sub>jlk</sub>(j≠l)表示由导频分配方案<img file="FDA00007455411600000222.GIF" wi="425" he="113" />确定的第j个小区和第l个小区中第k个使用相同导频的用户到第l个基站的大尺度衰落;6)选择操作,对适应值<img file="FDA00007455411600000223.GIF" wi="519" he="85" />进行排序,采用蒙特卡洛选择方法,保留<img file="FDA00007455411600000315.GIF" wi="57" he="40" />组基因编码<img file="FDA0000745541160000031.GIF" wi="888" he="128" />所述蒙特卡洛选择方法具体步骤如下:(1)首先归一化适应值<img file="FDA0000745541160000032.GIF" wi="566" he="77" />得到对应的归一化值{p<sub>c</sub>},c∈{1,2,...,C}其中p<sub>c</sub>∈[0,1]作为第c组基因编码<img file="FDA0000745541160000033.GIF" wi="408" he="111" />能保留到下一代的概率,即适应度;(2)第t代C组基因编码<img file="FDA0000745541160000034.GIF" wi="756" he="119" />分别按照适应度{p<sub>c</sub>},c∈{1,2,...,C}保留到下一代,表示为<img file="FDA0000745541160000035.GIF" wi="888" he="132" />这里集合<img file="FDA00007455411600000316.GIF" wi="63" he="59" />包含了被保留下的基因编码组;7)对保留下来的<img file="FDA00007455411600000317.GIF" wi="63" he="57" />组基因编码<img file="FDA0000745541160000036.GIF" wi="852" he="114" />依次进行复制、交叉和变异操作,生成新的C组基因编码<img file="FDA0000745541160000037.GIF" wi="928" he="130" />(1)复制操作:第t+1代的基因编码复制第t代保留下来的<img file="FDA00007455411600000314.GIF" wi="53" he="55" />组基因编码<img file="FDA0000745541160000038.GIF" wi="943" he="113" />并扩充至C组基因编码<maths num="0009" id="cmaths0009"><math><![CDATA[<mrow><mo>{</mo><msubsup><mi>T</mi><mrow><mi>c</mi><mn>1</mn></mrow><mrow><mo>(</mo><mi>t</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow></msubsup><mo>,</mo><msubsup><mi>T</mi><mrow><mi>c</mi><mn>2</mn></mrow><mrow><mo>(</mo><mi>t</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow></msubsup><mo>,</mo><mn>...</mn><mo>,</mo><msubsup><mi>T</mi><mrow><mi>c</mi><mi>L</mi></mrow><mrow><mo>(</mo><mi>t</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow></msubsup><mo>}</mo><mo>,</mo><mi>c</mi><mo>&Element;</mo><mo>{</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>,</mo><mn>...</mn><mo>,</mo><mi>C</mi><mo>}</mo><mo>;</mo></mrow>]]></math><img file="FDA0000745541160000039.GIF" wi="929" he="131" /></maths>(2)交叉操作:第t+1代的C组基因编码<img file="FDA00007455411600000310.GIF" wi="896" he="112" />之间进行交叉,即随机交换不同组基因编码中相同位置的某个基因,例如,将<img file="FDA00007455411600000311.GIF" wi="128" he="96" />和<img file="FDA00007455411600000312.GIF" wi="129" he="101" />进行交换,这里c≠c′;(3)变异操作:第t+1代的某一组基因编码<img file="FDA00007455411600000313.GIF" wi="545" he="112" />中某一条染色体上的某一个基因随机选择变异;8)判断t是否达到最大遗传代数t<sub>max</sub>,如果t=t<sub>max</sub>,执行步骤9;否则,遗传代数增加1,即t=t+1,转入步骤4);9)最终导频分配方案由<img file="FDA0000745541160000041.GIF" wi="565" he="98" />给出。
地址 250199 山东省济南市历城区山大南路27号