发明名称 分布式MIMO系统基站侧天线位置优化方法
摘要 本发明公开了一种分布式MIMO系统基站侧天线位置优化方法,包括步骤:A.初始阶段分布式天线端口DAU向移动台MS发送一个测试信号;B. 移动台MS响应测试信号后将目前信道状态信息反馈给分布式天线端口DAU;C. 智能天线通过对天线阵列接收的反馈信号在整个分布式天线端口DAU功率受限的条件下通过加权处理形成天线波束;D. 天线波束通过加权矢量调整天线方向图使其主瓣指向需要进行通信的移动台MS方向。进一步的运用天线选择策略和合理的天线功率分配方案,实现节能环保的绿色通信,并在此基础上通过采用粒子群算法对系统不同目标函数进行优化,获取天线在不同区域的最优位置,实现其最优覆盖。
申请公布号 CN103746729B 申请公布日期 2017.02.01
申请号 CN201310728444.2 申请日期 2013.12.25
申请人 河海大学 发明人 李岳衡;崔磊;彭文杰;王莉;羌佳琳;田月萍;蔡俊;付明浩;居美艳;蒋德富;黄平
分类号 H04B7/04(2006.01)I 主分类号 H04B7/04(2006.01)I
代理机构 南京经纬专利商标代理有限公司 32200 代理人 朱小兵
主权项 一种分布式MIMO系统基站侧天线位置优化方法,其特征在于:所述基站侧天线包括:配置于分布式天线端口DAU中,由两根及两根以上权重天线组成的天线阵列;以及配置于移动台MS中的单根天线;优化方法包括以下步骤:A.在系统正常工作时,由分布式天线端口DAU向移动台MS发送一个测试信号;B.移动台MS响应测试信号后将目前信道状态信息反馈给分布式天线端口DAU;C.通过对天线阵列接收的反馈信号在整个分布式天线端口DAU功率受限的条件下进行加权处理形成天线波束;D.通过加权矢量调整天线方向图使天线波束主瓣指向需要进行通信的移动台MS方向;所述的加权矢量调整的方法为:设第i个分布式天线端口DAU<sub>i</sub>与特定位置上的移动台MS之间的信道衰落系数为h<sub>i</sub>,各个分布式天线端口DAU<sub>i</sub>的加权矩阵为W=(w<sub>1</sub>,w<sub>2</sub>,w<sub>3</sub>,w<sub>4</sub>,w<sub>5</sub>),i=1,2,…,5;移动台MS的接收信号表示为:<maths num="0001"><math><![CDATA[<mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><msub><mi>y</mi><mn>1</mn></msub><mo>=</mo><msubsup><mi>h</mi><mn>1</mn><mi>H</mi></msubsup><mrow><mo>(</mo><msub><mi>w</mi><mn>1</mn></msub><msub><mi>s</mi><mn>1</mn></msub><mo>)</mo></mrow><mo>+</mo><mi>n</mi></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>y</mi><mn>2</mn></msub><mo>=</mo><msubsup><mi>h</mi><mn>2</mn><mi>H</mi></msubsup><mrow><mo>(</mo><msub><mi>w</mi><mn>2</mn></msub><msub><mi>s</mi><mn>2</mn></msub><mo>)</mo></mrow><mo>+</mo><mi>n</mi></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>y</mi><mn>3</mn></msub><mo>=</mo><msubsup><mi>h</mi><mn>3</mn><mi>H</mi></msubsup><mrow><mo>(</mo><msub><mi>w</mi><mn>3</mn></msub><msub><mi>s</mi><mn>3</mn></msub><mo>)</mo></mrow><mo>+</mo><mi>n</mi></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>y</mi><mn>4</mn></msub><mo>=</mo><msubsup><mi>h</mi><mn>4</mn><mi>H</mi></msubsup><mrow><mo>(</mo><msub><mi>w</mi><mn>4</mn></msub><msub><mi>s</mi><mn>4</mn></msub><mo>)</mo></mrow><mo>+</mo><mi>n</mi></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>y</mi><mn>5</mn></msub><mo>=</mo><msubsup><mi>h</mi><mn>5</mn><mi>H</mi></msubsup><mrow><mo>(</mo><msub><mi>w</mi><mn>5</mn></msub><msub><mi>s</mi><mn>5</mn></msub><mo>)</mo></mrow><mo>+</mo><mi>n</mi></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0001024983160000011.GIF" wi="406" he="454" /></maths>其中,y<sub>i</sub>是移动台MS的接收信号标量;h<sub>i</sub>是M<sub>i</sub>×1的信道向量,M<sub>i</sub>是分布式天线端口DAU<sub>i</sub>配置的天线数目;上标()<sup>H</sup>表示共轭转置;w<sub>i</sub>是各个发射天线端口上的加权列向量;s<sub>i</sub>是各发射天线端口上的发送符号;n是加性复高斯白噪声;E.对分布式天线端口DAU进行天线功率分配:通过移动台MS反馈的信道信息,分布式天线端口DAU在天线阵列中根据信道增益系数的大小为其分配相应的功率,信道的增益系数越小分配的功率越少;在功率受限的条件下,天线功率分配是按照各信道增益系数的比例来分配的,从而在特定方向上与移动台建立通信链路;F.移动台MS与接收信噪比最大的分布式天线端口DAU进行通信;所述的移动台MS接收端的信噪比SNR表示为:<img file="FDA0001024983160000012.GIF" wi="334" he="141" />其中,P<sub>i</sub>是天线端口的发送功率,σ是加性复高斯白噪声的方差;G.在分布式MIMO系统不同的区域内,通过粒子群算法对分布式MIMO系统的接收信噪比SNR、系统容量C、中断概率P<sub>out</sub>三个系统性能指标作为目标函数进行优化,获取天线在相应区域的最优位置;目标优化函数基于移动台MS特定位置下的数学表达式表示为:SNR(ρ,θ)=E<sub>h</sub>{r<sub>i</sub>}C(ρ,θ)=E<sub>h</sub>{log<sub>2</sub>(1+r<sub>i</sub>)}P<sub>out</sub>(ρ,θ)=E<sub>h</sub>{P<sub>r</sub>(r<sub>i</sub>≤r<sub>th</sub>)}其中,<img file="FDA0001024983160000021.GIF" wi="315" he="141" />为移动台MS接收端信噪比,E为数学期望,下标h表示信道系数,P<sub>r</sub>是r<sub>i</sub>≤r<sub>th</sub>的概率,r<sub>th</sub>表示中断门限;将在近似蜂窝的圆形小区内移动用户的分布概率密度函数记作f(ρ,θ),ρ为天线极半径,θ为天线极角,为了反映天线位置μ对系统性能的影响,将上述表达式对移动台MS位置进行统计平均,具体表达式如下<maths num="0002"><math><![CDATA[<mfenced open = "" close = ""><mtable><mtr><mtd><mrow><mover><mrow><mi>S</mi><mi>N</mi><mi>R</mi></mrow><mo>&OverBar;</mo></mover><mrow><mo>(</mo><mi>&mu;</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>E</mi><mrow><mi>&rho;</mi><mo>,</mo><mi>&theta;</mi></mrow></msub><mo>{</mo><mi>S</mi><mi>N</mi><mi>R</mi><mrow><mo>(</mo><mi>&rho;</mi><mo>,</mo><mi>&theta;</mi><mo>)</mo></mrow><mo>}</mo></mrow></mtd></mtr><mtr><mtd><mrow><mo>=</mo><munderover><mo>&Integral;</mo><mn>0</mn><mi>R</mi></munderover><munderover><mo>&Integral;</mo><mrow><mn>2</mn><mi>&pi;</mi></mrow><mn>0</mn></munderover><mi>S</mi><mi>N</mi><mi>R</mi><mrow><mo>(</mo><mi>&rho;</mi><mo>,</mo><mi>&theta;</mi><mo>)</mo></mrow><mi>f</mi><mrow><mo>(</mo><mi>&rho;</mi><mo>,</mo><mi>&theta;</mi><mo>)</mo></mrow><mi>d</mi><mi>&theta;</mi><mi>d</mi><mi>&rho;</mi></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0001024983160000022.GIF" wi="773" he="255" /></maths><maths num="0003"><math><![CDATA[<mfenced open = "" close = ""><mtable><mtr><mtd><mrow><mover><mi>C</mi><mo>&OverBar;</mo></mover><mrow><mo>(</mo><mi>&mu;</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>E</mi><mrow><mi>&rho;</mi><mo>,</mo><mi>&theta;</mi></mrow></msub><mo>{</mo><mi>C</mi><mrow><mo>(</mo><mi>&rho;</mi><mo>,</mo><mi>&theta;</mi><mo>)</mo></mrow><mo>}</mo></mrow></mtd></mtr><mtr><mtd><mrow><mo>=</mo><munderover><mo>&Integral;</mo><mn>0</mn><mi>R</mi></munderover><munderover><mo>&Integral;</mo><mrow><mn>2</mn><mi>&pi;</mi></mrow><mn>0</mn></munderover><mi>C</mi><mrow><mo>(</mo><mi>&rho;</mi><mo>,</mo><mi>&theta;</mi><mo>)</mo></mrow><mi>f</mi><mrow><mo>(</mo><mi>&rho;</mi><mo>,</mo><mi>&theta;</mi><mo>)</mo></mrow><mi>d</mi><mi>&theta;</mi><mi>d</mi><mi>&rho;</mi></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0001024983160000023.GIF" wi="661" he="253" /></maths><maths num="0004"><math><![CDATA[<mfenced open = "" close = ""><mtable><mtr><mtd><mrow><mover><msub><mi>P</mi><mrow><mi>o</mi><mi>u</mi><mi>t</mi></mrow></msub><mo>&OverBar;</mo></mover><mrow><mo>(</mo><mi>&mu;</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>E</mi><mrow><mi>&rho;</mi><mo>,</mo><mi>&theta;</mi></mrow></msub><mo>{</mo><msub><mi>P</mi><mrow><mi>o</mi><mi>u</mi><mi>t</mi></mrow></msub><mrow><mo>(</mo><mi>&rho;</mi><mo>,</mo><mi>&theta;</mi><mo>)</mo></mrow><mo>}</mo></mrow></mtd></mtr><mtr><mtd><mrow><mo>=</mo><munderover><mo>&Integral;</mo><mn>0</mn><mi>R</mi></munderover><munderover><mo>&Integral;</mo><mrow><mn>2</mn><mi>&pi;</mi></mrow><mn>0</mn></munderover><msub><mi>P</mi><mrow><mi>o</mi><mi>u</mi><mi>t</mi></mrow></msub><mrow><mo>(</mo><mi>&rho;</mi><mo>,</mo><mi>&theta;</mi><mo>)</mo></mrow><mi>f</mi><mrow><mo>(</mo><mi>&rho;</mi><mo>,</mo><mi>&theta;</mi><mo>)</mo></mrow><mi>d</mi><mi>&theta;</mi><mi>d</mi><mi>&rho;</mi></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0001024983160000024.GIF" wi="758" he="255" /></maths>将上述不同的目标函数统一用f<sub>obj</sub>表示,采用复化Simpson积分公式将目标函数表达式近似为:<maths num="0005"><math><![CDATA[<mrow><mover><msub><mi>f</mi><mrow><mi>o</mi><mi>b</mi><mi>j</mi></mrow></msub><mo>&OverBar;</mo></mover><mrow><mo>(</mo><mi>&mu;</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><mn>2</mn><mi>&pi;</mi><mi>R</mi></mrow><mrow><mn>9</mn><mi>P</mi><mi>Q</mi></mrow></mfrac><munderover><mo>&Sigma;</mo><mrow><mi>p</mi><mo>=</mo><mn>0</mn></mrow><mi>P</mi></munderover><munderover><mo>&Sigma;</mo><mrow><mi>q</mi><mo>=</mo><mn>0</mn></mrow><mi>Q</mi></munderover><msub><mi>&delta;</mi><mrow><mi>p</mi><mo>,</mo><mi>q</mi></mrow></msub><msub><mi>&rho;</mi><mi>p</mi></msub><msub><mi>f</mi><mrow><mi>o</mi><mi>b</mi><mi>j</mi></mrow></msub><mrow><mo>(</mo><msub><mi>&rho;</mi><mi>p</mi></msub><mo>,</mo><msub><mi>&theta;</mi><mi>q</mi></msub><mo>)</mo></mrow><mi>f</mi><mrow><mo>(</mo><msub><mi>&rho;</mi><mi>p</mi></msub><mo>,</mo><msub><mi>&theta;</mi><mi>q</mi></msub><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001024983160000025.GIF" wi="1035" he="151" /></maths>其中,R是圆形小区半径,P和Q分别是划分极半径ρ和极角θ的等距节点数,δ<sub>p,q</sub>是权值,为矩阵A的第p+1行,第q+1列的元素,p=0,1,2,…,P;q=0,1,2,…,Q;矩阵<img file="FDA0001024983160000026.GIF" wi="718" he="71" />其中[·]<sup>T</sup>表示矩阵转置,<img file="FDA0001024983160000027.GIF" wi="38" he="39" />为克罗内克积。
地址 210098 江苏省南京市鼓楼区西康路1号