专利名称：盲目适应波束形成技术的制作方法
技术领域：
本发明涉及一种智能型天线系统，特别地涉及一种用于诸如WCDMA和CDMA-2000的无线通信系统的智能型天线系统中的增强的盲自适应波束形成技术。
智能天线系统通常与基站共址，其将天线阵与数字信号处理能力结合起来，以便能以自适应的、空间敏感的方式发送和接收。换言之，该系统能够根据所在的信号环境自动地变换它的辐射图的定向。这样就可以明显地提高无线系统的性能特征(如容量)。
自适应波束形成天线系统是所述智能型天线系统中的一种。
自适应天线技术代表了目前最先进的智能型天线方案。自适应系统使用了各种新的信号处理算法，它的优势在于其能够有效地定位和跟踪各种类型的信号，以动态地使干扰达到最小，使预定的信号接收达到最大。
在自适应天线系统中，通过用于控制各种复杂运行条件的组合的先进的技术(算法)提高了基本的信号处理能力。典型地，将来自空间分布的天线单元的接收的信号乘以一个加权矢量，对幅度和相位进行复合的调整。将这些信号组合，产生天线阵的输出。自适应算法根据预先确定的目标控制加权矢量。这些动态的计算使系统能够改变其辐射方向图，以使信号的接收最佳化。
自适应波束形成算法分为基于定方的波束成形、基于训练序列的波束成形和盲波束形成算法。
所述基于定方的波束形成算法检测来自移动终端的信号的各种到达角度，以形成最佳的波束方向图。但是，该算法要求阵列校准具有高的精确度，因此在无线通信系统中有限制地使用。
所述基于训练序列的波束形成算法使用基准信号(训练信号)来估算加权矢量，以便根据加权矢量形成波束方向图。
最后，所述盲波束形成算法的特征在于，在加权矢量的估算中无需使用基准(训练)信号。
本发明的波束形成算法属于盲波束形成算法。
已研究和开发了很多的用于形成最佳波束方向图的波束形成算法。Seungwon Chui et.al.开发了一种自适应的波束形成算法(Seungwon Choi和Donghee Shim的“一种用于CDMA移动通信环境中的智能型天线系统的新的自适应波束形成算法”，IEEE Trans.OnVehicular Technology，Vol.49，No.5 pp1793-1806，2000年9月)，它的特征在于将信噪比最大(SNR)，而不是使信号与干扰加上噪音的比(SINR)最大。但是，该波束成形算法在高能量干扰的环境中不能很好地工作，因为该波束成形算法使SNR最大。
Donghee Shim et.al.开发了另一个波束成形算法(Donghee Shim，Fakhrul Alam和Jeongho Kim的“采用盲自适应算法的智能型天线系统的性能分析”，VT MPRG workshop，2001)，它的特征在于，根据自相关矩阵的倒置使SINR最大。但是，由于该波束形成算法在更新加权矢量的等式中使用近似值，因此可对系统的性能造成不利影响。
本发明的目的是提供一种自适应的波束形成方法，在宽带码分多址(WCDMA)、CDMA-2000等的智能型天线系统中使用梯度搜索算法将SINR作为性能数据，使SINR最大。
为实现上述目的，本发明的自适应波束形成方法包括以下步骤设置初始加权矢量，根据新的信号的矢量计算各种矢量和/或标量的值，使用新信号的矢量和标量的值在梯度搜索方法中更新加权矢量，使加权矢量归一化，并基于更新的加权矢量形成最佳波束的方向图。
根据以下cost函数的最大，将所述加权矢量设置为使SINR最大J(wk,pλ)=wk,pHRsswk,p+λ(1-wk,pHRuuwk,p)]]>其中H是复共轭转置矩阵(厄密共轭算子)，λ是拉格朗日算子，Wk，p是第k个用户终端的第p个多路径信号的加权矢量，Rss是所需信号的自协方差矩阵，以及Ruu是不需要信号的自协方差矩阵，即干扰信号和噪声的和(u)。
所述梯度搜索方法由下列等式构成wk,p(m+1)=wk,p(m)+12μ▿w]]>=wk,p(m)+μ(Rss(m)wk,p(m)-λ(m)Ruu(m)wk,p(m))]]>这里，w(m)是函数J(Wk，p，λ)在瞬态(snapshot)指数m的梯度矢量，μ是步进尺寸。
根据以下等式，分别计算所需的和不需要的信号的自相关矩阵Rss，k，p(m)和Ruu，k，p(m)。
Rss，h，p(m)=fRss，k，p(m-1)+s(m)sH(m)和Ruu，k，p(m)＝fRuu，k，p(m-1)+μ(m)μH(m)其中f是忘记因子(forgetting factor)(0≤f≤1)。
可用以下等式表示λ(m)λ(m)=wm,pH(m)Rss,k,p(m)wk,p(m)wk,pH(m)Ruu,k,p(m)wk,p(m)]]>用以下等式更新加权矢量wk,p(m+1)=wk,p(m)+μ[y~k,p(m)-λ(m)q~k,p(m)]]]>通过以下等式使更新的加权矢量归一化wk,p(m+1)=wk,p(m+1)||wk,p(m+1)||]]>
图1是一个表示了根据本发明的优选实施例的盲自适应波束形成方法的流程图。
在本发明所述的自适应波束形成方法中，可以用采用拉格朗日算子的cost函数将按如下方式计算加权矢量，以使SINR最大J(wk,p·λ)=wk,pHRsswk,p+λ(1-wk,pHRuuwk,p)---(1)]]>其中H是复共轭转置矩阵(厄密共轭算子)，λ是拉格朗日算子，Wk，p是第k个用户终端的第p个多路径信号的加权矢量，Rss是所需信号的自协方差矩阵，以及Ruu是不需要的信号的自协方差矩阵，即干扰信号和噪声的和(u)。
可用以下公式表示所述加权矢量Wk，pwk,p=[wk,p(0)wk,p1···wk,pN-1]H---(2)]]>其中N是天线的数量。从等式(2)中，使函数J(Wk，p，λ)最大的加权矢量Wk，p可通过约束wk,pHRuuwk,p10=1]]>条件得出。
Rsswk，p＝λRuuwk，p(3)从等式(1)中，可通过如下等式计算函数J(Wk，p，λ)的与wk，p*有关的梯度矢量w(m)w(m)=2(Rss(m)Wk，p(m)-λ(m)Ruu(m)Wm，p(m))(4)其中m是瞬态指数，*是复共轭。
根据梯度搜索方法，按如下方式更新加权矢量wk,p(m+1)=wk,p(m)+12μ▿w---(5)]]>=wk,p(m)+μ(Rss(m)wk,p(m)-λ(m)Ruu(m)wk,p(m))]]>其中，μ为步进尺寸。可通过如下等式获得协方差矩阵Rss和RuuRss，k，p(m)=fRss，k，p(m-1)+s(m)sH(m) (6)Ruu，k，p(m)＝fRuu，k，p(m-1)+u(m)uH(m) (7)其中f是忘记因子(0≤f≤1)。从等式(5)中，λ(m)可表示如下λ(m)=wm,pH(m)Rss,k,p(m)wk,p(m)wk,pH(m)Ruu,k,p(m)wk,p(m)---(8)]]>同样，按如下方式将所述加权矢量wk，p归一化wk,p(m+1)=wk,p(m+1)||wk,p(m+1)||---(9)]]>使用协方差矩阵，可以推出能够减少盲梯度搜索方法(BGSA)的计算复杂程度的线性化的盲自适应算法。
可按如下方式定义瞬态指数m中的拉格朗日算子λλ(m)=λnum(m)λden(m)---(10)]]>由公式(10)，将λnum(m)表示如下λnum(m)=wk,pH(m)Rss(m)wk,p(m)]]>=wm,pH(m)≥[fRss(m-1)+s(m)sH(m)]wk,p(m)---(11)]]>=fwk,pH(m)Rss(m-1)wk,p(m)+wk,pH(m)s(m)sH(m)wk,p(m)]]>≈fwk,pH(m-1)Rss(m-1)wk,p(m-1)+wk,pH(m)s(m)sH(m)wk,p(m)]]>同样，从等式(11)中，可将λnum(m)表示如下λnum(m)＝fλnum(m-1)+|yk，p(m)|2(12)从等式(12)中，可yk，p(m)将表示如下yk,p(m)=sk,pH(m)wk,p(m)---(13)]]>以同样的方式，将λden表示如下λden(m)＝fλden(m-1)+|qk，p(m)|2(14)从等式(14)中，可将qk，p(m)表示如下qk，p(m)＝uH(m)wk，p(m) (15)同样，可将加权矢量更新等式(5)写成如下wk，p(m+1)＝wk，p(m)+μzk，p(m) (16)在等式(16)中，可将zk，p(m)表示如下zk，p(m)＝[Rss(m)-λ(m)Ruu(m)]wk，p(m) (17)在等式(17)中，可将zk，p(m)表示为zk，p(m)＝［fRss(m-1)+s(m)sH(m)-λ(m)(fRuu(m-1)+u(m)uH(m))]wk，p(m)＝[(fRss(m-1)wk，p(m)+s(m)sH(m)wk，p(m))-fλ(m)Ruu(m-1)wk，p(m)-λ(m)u(m)wk，p(m)](18)从等式(13)和(18)中，可将 表示为y~k,p(m)=fRss(m-1)wk,p(m)+s(m)sH(m)wk,p(m)]]>≈fRss(m-1)wk,p(m-1)+s(m)sH(m)wk,p(m)---(19)]]>=fy~k,p(m-1)+yk,p(m)s(m)]]>其中可将初始值表示为y~k,p(0)=yk,p(0)s(0).]]>以相同的方式，从公式(15)和(18)可将 表示如下q~k,p(m)=fRuu(m-1)wk,p(m)+u(m)uH(m)wk,p(m)---(20)]]>=fq~k,p(m-1)+qk,p(m)u(m)]]>其中，可将初始值表示为q~k,p(0)=qk,p(0)u(0).]]>相应地，可将zk，p(m)表示为zk,p(m)=[y~k,p(m)-λ(m)qk,p(m)]---(21)]]>将参照图1对本发明所述的盲自适应波束形成方法进行描述。
当本发明所述的盲自适应波束形成算法用于接收的信号时，利用如下(22)-(27)的等式计算的初始矢量和/或标量，在步骤S101来设置初始加权矢量wk，p(0)wk，p(0)＝s(0) (22)yk，p(0)＝sH(0)wk，p(0)(23)qk，p(0)＝uH(0)wk，p(0)(24)λnum(0)＝|yk，p(0)|2(25)λden(0)=|qk,p(0)|2---(26)]]>λ(0)=λnum(0)λden(0)---(27)]]>接收到新的信号时，就在步骤S102按如下方法使用新信号的信号矢量更新以前的矢量和/或标量yk，p(m)＝sH(m)wk，p(m) (28)qk，p(m)＝uH(m)wk，p(m) (29)λnum(m)＝fλnum(m-1)+|yk，p(m)|2(30)λden(m)＝fλden(m-1)+|qk，p(m)|2(31)λ(m)=λnum(m)λden(m)---(32)]]>y~(m)=fy~k,p(m-1)+yk,p(m)sk,p(m)---(33)]]>q~k,p(m)=fq~k,p(m-1)+qk,p(m)u(m)---(34)]]>上述值更新后，在步骤S103中，采用梯度搜索方法中的更新的矢量和/或标量值来更新加权矢量如下wk,p(m+1)=wk,p(m)+μ[y~k,p(m)-λ(m)q~k,p(m)]---(35)]]>更新后，在步骤S104，按如下方法对加权矢量进行归一化wk,p(m+1)=wk,p(m+1)||wk,p(m+1)||---(36)]]>如果继续进行自适应的更新程序，则自适应波束成形算法重复地回到步骤S102。
一旦计算出加权矢量，在步骤S105，根据更新的加权矢量计算天线阵的输出。
当结合目前认为最为实用的和优选的实施例，对本发明作了描述后，可以理解，本发明不仅限于公开的实施例，相反，本发明包括在权利要求中所述的精神和范围之内的各种修改和等效设置。
如上所述，在本发明的盲自适应波束形成算法中，由于形成最佳波束方向图的加权矢量是在梯度搜索方法中计算的，因此有可能使同信道的干扰和噪声最小，结果，使系统容量、覆盖区域和可靠的高速数据通信服务达到最佳。
同样，本发明所述的盲自适应波束形成算法可用于第四代无线通信系统中，也可用于诸如WCDMA、CDMA-2000的其他无线系统中，以使系统性能最佳。
权利要求
1.一种自适应波束形成方法，包括步骤设置初始加权矢量；根据新的信号的矢量计算各种矢量和/或标量的值；使用新信号的矢量和标量的值在梯度搜索方法更新加权矢量；根据更新的加权矢量形成最佳波束方向图。
2.根据权利要求1所述的自适应波束形成方法，其中所述加权矢量设置为使SINR最大。
3.根据权利要求2所述的自适应波束形成方法，其中根据以下cost函数的最大使SINR最大J(wk,p·λ)=wk,pHRsswk,p+λ(1-wk,pHRuuwk,p)]]>其中H是复共轭转置矩阵(厄密共轭算子)，λ是拉格朗日算子，wk，p是第k个用户终端的第p个多径信号的加权矢量，Rss是所需信号的自协方差矩阵，以及Ruu是不需要的信号的自协方差矩阵，即干扰信号和噪声的和(u)。
4.根据权利要求2所述的自适应波束形成方法，其中使用拉格朗日算子使SINR最大。
5.根据权利要求3所述的自适应波束形成方法，其中梯度搜索方法用以下等式表示wk,p(m+1)=wk,p(m)+12μ▿w]]>=wk,p(m)+μ(Rss(m)wk,p(m)-λ(m)Ruu(m)wk,p(m))]]>其中，w(m)是函数J(Wk，p，λ)瞬态指数m的梯度矢量，以及μ是步进尺寸。
6.根据权利要求5所述的自适应波束形成方法，其中根据以下等式，分别计算所需的和不需要的信号的自协方差矩阵Rss，k，p(m)和Ruu，k，p(m)Rss，k，p(m)＝fRss，k，p(m-1)+s(m)sH(m)和Ruu，k，p(m)＝fRuu，k，p(m-1)+u(m)uH(m)其中f是忘记因子(0≤f≤1)。
7.根据权利要求5所述的自适应波束形成方法，其中用以下等式表示λ(m)λ(m)=wm,pH(m)Rss,k,p(m)wk,p(m)wk,pH(m)Ruu,k,p(m)wk,p(m)]]>
8.根据权利要求6所述的自适应波束形成方法，其中通过下列等式更新所述加权矢量wk,p(m+1)=wk,p(m)+μ[y~k,p(m)-λ(m)q~k,p(m)]]]>
9.一种自适应波束形成方法，包括步骤设置初始加权矢量；根据新的信号的矢量计算各种矢量和/或标量的值；使用新信号的矢量和标量的值在梯度搜索方法中更新加权矢量；使所述加权矢量归一化；和根据更新的加权矢量形成最佳波束方向图。
10.根据权利要求9所述的自适应波束形成方法，其中所述加权矢量设置为使SINR最大。
11.根据权利要求10所述的自适应波束形成方法，其中根据以下cost函数的最大使SINR最大J(wk,p·λ)=wk,pHRsswk,p+λ(1-wk,pHRuuwk,p)]]>其中H是复共轭转置矩阵(厄密共轭算子)，λ是拉格朗日算子，wk，p是第k个用户终端的第p个多径信号的加权矢量，Rss是所需信号的自协方差矩阵，以及Ruu是不需要的信号的自协方差矩阵，即干扰信号和噪声的和(u)。
12.根据权利要求10所述的自适应波束形成方法，其中使用拉格朗日算子使SINR最大。
13.根据权利要求11所述的自适应波束形成方法，其中梯度搜索方法用以下等式表示wk,p(m+1)=wk,p(m)+12μ▿w]]>=wk,p(m)+μ(Rss(m)wk,p(m)-λ(m)Ruu(m)wk,p(m))]]>其中，w(m)是函数J(Wk，p，λ)瞬态图指数m的梯度矢量，以及μ是步进尺寸。
14.根据权利要求13所述的自适应波束形成方法，其中根据以下等式，分别计算所需的和不需要的信号的自协方差矩阵Rss，k，p(m)和Ruu，k，p(m)Rss，k，p(m)＝fRss，k，p(m-1)+s(m)sH(m)和Ruu，k，p(m)＝fRuu，k，p(m-1)+u(m)uH(m)其中f是忘记因子(0≤f≤1)。
15.根据权利要求14所述的自适应波束形成方法，其中用以下等式表示λ(m)λ(m)=wm,pH(m)Rss,k,p(m)wk,p(m)wk,pH(m)Ruu,k,p(m)wk,p(m)]]>
16.根据权利要求15所述的自适应波束形成方法，其中通过下列等式更新所述加权矢量wk,p(m+1)=wk,p(m)+μ[y~k,p(m)-λ(m)q~k,p(m)]]]>
17.根据权利要求15所述的自适应波束形成方法，其中通过以下等式将所述加权矢量归一化wk,p(m+1)=wk,p(m+1)||wk,p(m+1)||]]>
全文摘要
本发明公开了一种自适应波束形成方法，它包括以下步骤设置初始加权矢量，根据新的信号的矢量计算各种矢量和/或标量的值，使用新信号的矢量和标量的值在梯度搜索方法中更新所述加权矢量，使所述加权矢量归一化，并基于更新的加权矢量形成最佳波束方向图。由于本发明的形成最佳波束方向图的所述加权矢量是按照盲自适应波束形成算法中的梯度搜索方法计算的，因此可使同信道的干扰和噪声最小，结果，使系统容量、覆盖区域和可靠的高速数据通信业务达到最佳。
文档编号H04J13/00GK1467941SQ0310845
公开日2004年1月14日 申请日期2003年4月11日 优先权日2002年6月17日
发明者金相俊 申请人:Lg电子株式会社