主权项 |
一种分布式天线系统的功率分配方法,其特征在于,包括如下步骤: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>∈</mo><msub><mi>Γ</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>Γ</mi><mi>j</mi></msub><mo>=</mo><mo>{</mo><mi>m</mi><mo>∈</mo><msub><mi>Γ</mi><mrow><mi>j</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>,</mo><mi>m</mi><mo>∉</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>μ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>μ</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>μ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>≠</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>∈</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>Γ</mi><mi>j</mi></msub></msub><mo>=</mo><msub><mi>max</mi><mrow><mi>m</mi><mo>∈</mo><msub><mi>Γ</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>γ</mi><munderover><mi>Σ</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>σ</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,ε为特定的门限值。 |