发明名称 一种分布式天线系统及其功率分配方法
摘要 本发明提供了一种分布式天线系统的功率分配方法,根据每个用户的信道向量,组成一个等效的信道矩阵,并计算出这个等效的信道矩阵的伪逆,并进行能量归一化,即得到ZF波束形成矩阵,根据已有信息和求得的波束形成矩阵,进行初始化,采用迭代方法求解功率分配矩阵,最后波束形成矩阵和对角化的功率分配矩阵的开方相乘求得ZF发射预编码矩阵。本发明还提供了采用上述功率分配方法的分布式天线系统。与最优功率分配方法相比,本发明所提出的功率分配方法能够以很小的性能损失为代价有效地降低计算复杂度。
申请公布号 CN102468876B 申请公布日期 2015.03.18
申请号 CN201010532651.7 申请日期 2010.11.01
申请人 鼎桥通信技术有限公司 发明人 侯雪颖;王大飞;沈东栋
分类号 H04B7/04(2006.01)I;H04L1/00(2006.01)I;H04W52/18(2009.01)I 主分类号 H04B7/04(2006.01)I
代理机构 北京德琦知识产权代理有限公司 11018 代理人 牛峥;王丽琴
主权项 一种分布式天线系统的功率分配方法,其特征在于,包括如下步骤:A、将分布式天线系统中每个用户的信道向量合并,组成一个等效的信道矩阵H=[h<sub>1</sub>,…,h<sub>K</sub>]<sup>T</sup>,<img file="FDA0000590235730000011.GIF" wi="207" he="85" />为第k个用户与基站的多天线间信道向量,M为分布式天线系统中分布式天线的总数;计算出这个等效的信道矩阵的伪逆,并进行能量归一化,得到迫零波束形成矩阵<img file="FDA0000590235730000012.GIF" wi="456" he="166" />迫零波束形成矩阵G的第m行第k列的元素为g<sub>m</sub>,<sub>k</sub>;B、定义矩阵A的第m行第k列的元素为<img file="FDA0000590235730000013.GIF" wi="352" he="156" />定义S<sub>j</sub>为已选天线序号的集合,Τ<sub>j</sub>为剩余天线序号的集合,P为每个分布式天线的最大发射功率;设置<img file="FDA0000590235730000014.GIF" wi="185" he="74" /><img file="FDA0000590235730000015.GIF" wi="61" he="61" />为空集,Τ<sub>0</sub>={1,…,M},<img file="FDA0000590235730000016.GIF" wi="533" he="127" />j=0,根据<img file="FDA0000590235730000017.GIF" wi="837" he="193" />以及w<sub>j</sub>Ap<sup>T</sup>(w<sub>j</sub>)=1求出功率分配的初始值p(w<sub>0</sub>)=[p<sub>1</sub>(w<sub>0</sub>),…,p<sub>K</sub>(w<sub>0</sub>)],并进一步求出q<sub>m</sub>(w<sub>0</sub>)=a<sub>m</sub>p<sup>T</sup>(w<sub>0</sub>)‑1,m=1,…,M,(x)<sup>+</sup>=max(x,0);a<sub>m</sub>为矩阵A的第m行,<img file="FDA0000590235730000018.GIF" wi="146" he="69" />表示其类型为1×M的向量,σ<sub>k</sub><sup>2</sup>为第k个用户的接收机热噪声功率,w<sub>0</sub>为初始加权矩阵,j为迭代次数,w<sub>j</sub>为第j次迭代的加权矩阵,p<sub>k</sub>(w<sub>j</sub>)为第j次迭代分配给用户k的功率,w<sub>j,m</sub>为w<sub>j</sub>的第m个元素,p(w<sub>j</sub>)为第j次迭代分配给用户的功率向量,γ通过由wAp<sup>T</sup>(w)=1求得,其中p(w)=[p<sub>1</sub>(w),…,p<sub>K</sub>(w)];K为用户数目;C、设j=j+1,<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>s</mi><mi>j</mi></msub><mo>=</mo><mi>arg</mi><munder><mi>max</mi><mrow><mi>m</mi><mo>&Element;</mo><msub><mi>&Gamma;</mi><mrow><mi>j</mi><mo>-</mo><mn>1</mn></mrow></msub></mrow></munder><msub><mi>q</mi><mi>m</mi></msub><mrow><mo>(</mo><msub><mi>w</mi><mrow><mi>j</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>)</mo></mrow><mo>,</mo></mrow>]]></math><img file="FDA0000590235730000019.GIF" wi="494" he="105" /></maths>更新S<sub>j</sub>=S<sub>j‑1</sub>∪{s<sub>j</sub>},<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>&Gamma;</mi><mi>j</mi></msub><mo>=</mo><mo>{</mo><mi>m</mi><mo>&Element;</mo><msub><mi>&Gamma;</mi><mrow><mi>j</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>,</mo><mi>m</mi><mo>&NotElement;</mo><msub><mi>S</mi><mi>j</mi></msub><mo>}</mo><mo>;</mo></mrow>]]></math><img file="FDA00005902357300000110.GIF" wi="486" he="81" /></maths>根据公式<maths num="0003" id="cmaths0003"><math><![CDATA[<mfenced open='{' close=''><mtable><mtr><mtd><msub><mi>w</mi><mrow><mi>j</mi><mo>,</mo><mi>m</mi></mrow></msub><mo>=</mo><msub><mi>&mu;w</mi><mrow><mi>j</mi><mo>-</mo><mn>1</mn><mo>,</mo><mi>m</mi></mrow></msub><mo>+</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mi>&mu;</mi><mo>)</mo></mrow><mo>,</mo></mtd><mtd><mi>m</mi><mo>=</mo><msub><mi>s</mi><mi>j</mi></msub></mtd></mtr><mtr><mtd><msub><mi>w</mi><mrow><mi>j</mi><mo>,</mo><mi>m</mi></mrow></msub><mo>=</mo><msub><mi>&mu;w</mi><mrow><mi>j</mi><mo>-</mo><mn>1</mn><mo>,</mo><mi>m</mi></mrow></msub><mo>,</mo></mtd><mtd><mi>m</mi><mo>&NotEqual;</mo><msub><mi>s</mi><mi>j</mi></msub></mtd></mtr></mtable></mfenced>]]></math><img file="FDA00005902357300000111.GIF" wi="670" he="149" /></maths>更新w<sub>j</sub>;D、根据更新后的w<sub>j</sub>再计算出q<sub>m</sub>(w<sub>j</sub>),m={1,…,M},并定义<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><mi>q</mi><msub><mi>s</mi><mi>j</mi></msub></msub><mo>=</mo><msub><mi>max</mi><mrow><mi>m</mi><mo>&Element;</mo><msub><mi>S</mi><mi>j</mi></msub></mrow></msub><msub><mi>q</mi><mi>m</mi></msub><mrow><mo>(</mo><msub><mi>w</mi><mi>j</mi></msub><mo>)</mo></mrow><mo>,</mo><msub><mi>q</mi><msub><mi>&Gamma;</mi><mi>j</mi></msub></msub><mo>=</mo><msub><mi>max</mi><mrow><mi>m</mi><mo>&Element;</mo><msub><mi>&Gamma;</mi><mi>j</mi></msub></mrow></msub><msub><mi>q</mi><mi>m</mi></msub><mrow><mo>(</mo><msub><mi>w</mi><mi>j</mi></msub><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA00005902357300000112.GIF" wi="969" he="91" /></maths>E、判断是否<img file="FDA0000590235730000021.GIF" wi="236" he="76" />若是,返回步骤C;否则,由<img file="FDA0000590235730000022.GIF" wi="830" he="196" />及w<sub>j</sub>Ap<sup>T</sup>(w<sub>j</sub>)=1计算出p(w<sub>j</sub>),并得到最终的功率分配结果为:<img file="FDA0000590235730000023.GIF" wi="475" he="157" />则功率分配矩阵P=diag{p<sub>k</sub>};F、计算出迫零预编码矩阵为<img file="FDA0000590235730000024.GIF" wi="238" he="82" />并以该迫零预编码矩阵对待发射信号进行预编码;步骤C中,当j=1的时候,μ=0;当j>1的时候,由下述流程计算得到μ:C1、定义μ<sub>min</sub>=0和μ<sub>max</sub>=1;C2、计算μ=(μ<sub>min</sub>+μ<sub>max</sub>)/2,更新w<sub>j</sub>,并且由<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msub><mi>p</mi><mi>k</mi></msub><mrow><mo>(</mo><msub><mi>w</mi><mi>j</mi></msub><mo>)</mo></mrow><mo>=</mo><msup><mrow><mo>(</mo><mfrac><mn>1</mn><mrow><mi>&gamma;</mi><munderover><mi>&Sigma;</mi><mrow><mi>m</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><msub><mi>w</mi><mrow><mi>j</mi><mo>,</mo><mi>m</mi></mrow></msub><msub><mi>a</mi><mrow><mi>m</mi><mo>,</mo><mi>k</mi></mrow></msub></mrow></mfrac><mo>-</mo><msubsup><mi>&sigma;</mi><mi>k</mi><mn>2</mn></msubsup><mo>)</mo></mrow><mo>+</mo></msup><mo>,</mo><mi>k</mi><mo>=</mo><mn>1</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>K</mi><mo>,</mo></mrow>]]></math><img file="FDA0000590235730000025.GIF" wi="865" he="191" /></maths>及w<sub>j</sub>Ap<sup>T</sup>(w<sub>j</sub>)=1,来求得p(w<sub>j</sub>)和q<sub>m</sub>(w<sub>j</sub>),m∈S<sub>j</sub>;C3、如果<img file="FDA0000590235730000026.GIF" wi="471" he="91" />更新μ<sub>max</sub>=μ,否则,更新μ<sub>min</sub>=μ;C4、判断是否μ<sub>max</sub>‑μ<sub>min</sub><ε,若是,则完成对μ的计算,否则,返回步骤C2,ε为特定的门限值。
地址 100102 北京市朝阳区望京北路9号叶青大厦D座15层
您可能感兴趣的专利