发明名称 用预条件共轭梯度法的码分多址通讯系统多用户检测方法
摘要 一种预条件共轭梯度算法的码分多址通讯系统的多用户检测方法,其包含以下步骤:1.检测所有用户所占全部码道的扩频信息、加扰信息、信道冲击响应信息;2.产生系统的自相关矩阵R,该矩阵R反映上述用户的扩频信息、加扰信息、信道冲击响应信息;3.生成接收到的用户数据矢量d<SUB>maf</SUB>:并以一定的排列规律排列,形成矢量d<SUB>maf</SUB>;4.由上述矩阵R、用户数据矢量d<SUB>maf</SUB>,得到要求解的计算各个码道发送数据的系统方程Rd=d<SUB>maf</SUB>;5.预条件,对系统相关矩阵R和用户数据矢量d<SUB>maf</SUB>进行预处理,得到矩阵R′和矢量d<SUB>maf</SUB>′,得出方程的方程R′d′=d<SUB>maf</SUB>′;6.通过上述解出的方程的解d′,通过反预条件d=(P<SUP>H</SUP>)<SUP>-1</SUP>d′得到解出d。
申请公布号 CN100365943C 申请公布日期 2008.01.30
申请号 CN200410052844.7 申请日期 2004.07.14
申请人 凯明信息科技股份有限公司 发明人 王炳立;李兴伟
分类号 H04B1/707(2006.01);H04J13/02(2006.01) 主分类号 H04B1/707(2006.01)
代理机构 上海新天专利代理有限公司 代理人 衷诚宣
主权项 1.一种利用预条件共轭梯度算法的码分多址通讯系统的多用户检测方法,其包含以下步骤:步骤1.在码分多址的通讯系统中,接收装置检测所有用户所占全部码道的扩频信息、加扰信息、信道冲击响应信息;步骤2.接收装置产生系统的自相关矩阵R,该矩阵R反映上述用户的扩频信息、加扰信息、信道冲击响应信息;步骤3.生成接收到的用户数据矢量dmaf:接收装置使用匹配滤波器或者瑞克接收机,检测得到接收信号经过解扩后所有激活的用户数据,并按一定的排列规律排列,形成矢量dmaf;步骤4.由上述系统相关矩阵R、用户数据矢量dmaf,得到需要求解的计算各个码道发送数据的系统方程<maths num="0001"><math>&lt;mrow&gt;&lt;mi&gt;R&lt;/mi&gt;&lt;mover&gt;&lt;mi&gt;d&lt;/mi&gt;&lt;mo&gt;^&lt;/mo&gt;&lt;/mover&gt;&lt;mo&gt;=&lt;/mo&gt;&lt;msub&gt;&lt;mi&gt;d&lt;/mi&gt;&lt;mi&gt;maf&lt;/mi&gt;&lt;/msub&gt;&lt;/mrow&gt;</math></maths>步骤5.多用户检测模块进行预条件,对系统相关矩阵R和用户数据矢量d<sub>maf</sub>进行预处理,得到矩阵R′和矢量d′<sub>maf</sub>,得出方程的方程R′d′=d′<sub>maf</sub>;其中,预条件,对系统相关矩阵R和接受数据矢量进行预处理是指:对Rd=d<sub>maf</sub>进行预条件,其中,<img file="C2004100528440002C2.GIF" wi="673" he="307" />式中:N是系统所规定的每个用户的符号数,K是步骤1所检测到的所有用户所占全部码道数目;5.1:定义P为预条件矩阵;5.2:使用以下变换:<maths num="0002"><math>&lt;mrow&gt;&lt;msup&gt;&lt;mi&gt;P&lt;/mi&gt;&lt;mrow&gt;&lt;mo&gt;-&lt;/mo&gt;&lt;mn&gt;1&lt;/mn&gt;&lt;/mrow&gt;&lt;/msup&gt;&lt;mi&gt;R&lt;/mi&gt;&lt;mrow&gt;&lt;mo&gt;(&lt;/mo&gt;&lt;msup&gt;&lt;mi&gt;P&lt;/mi&gt;&lt;msup&gt;&lt;mi&gt;H&lt;/mi&gt;&lt;mrow&gt;&lt;mo&gt;-&lt;/mo&gt;&lt;mn&gt;1&lt;/mn&gt;&lt;/mrow&gt;&lt;/msup&gt;&lt;/msup&gt;&lt;msup&gt;&lt;mi&gt;P&lt;/mi&gt;&lt;mi&gt;H&lt;/mi&gt;&lt;/msup&gt;&lt;mo&gt;)&lt;/mo&gt;&lt;/mrow&gt;&lt;mi&gt;d&lt;/mi&gt;&lt;mo&gt;=&lt;/mo&gt;&lt;mrow&gt;&lt;mo&gt;(&lt;/mo&gt;&lt;msup&gt;&lt;mi&gt;P&lt;/mi&gt;&lt;mrow&gt;&lt;mo&gt;-&lt;/mo&gt;&lt;mn&gt;1&lt;/mn&gt;&lt;/mrow&gt;&lt;/msup&gt;&lt;msub&gt;&lt;mi&gt;d&lt;/mi&gt;&lt;mi&gt;maf&lt;/mi&gt;&lt;/msub&gt;&lt;mo&gt;)&lt;/mo&gt;&lt;/mrow&gt;&lt;/mrow&gt;</math></maths><maths num="0003"><math>&lt;mrow&gt;&lt;mrow&gt;&lt;mo&gt;(&lt;/mo&gt;&lt;msup&gt;&lt;mi&gt;P&lt;/mi&gt;&lt;mrow&gt;&lt;mo&gt;-&lt;/mo&gt;&lt;mn&gt;1&lt;/mn&gt;&lt;/mrow&gt;&lt;/msup&gt;&lt;mi&gt;R&lt;/mi&gt;&lt;msup&gt;&lt;mi&gt;P&lt;/mi&gt;&lt;msup&gt;&lt;mi&gt;H&lt;/mi&gt;&lt;mrow&gt;&lt;mo&gt;-&lt;/mo&gt;&lt;mn&gt;1&lt;/mn&gt;&lt;/mrow&gt;&lt;/msup&gt;&lt;/msup&gt;&lt;mo&gt;)&lt;/mo&gt;&lt;/mrow&gt;&lt;mo&gt;&amp;CenterDot;&lt;/mo&gt;&lt;mrow&gt;&lt;mo&gt;(&lt;/mo&gt;&lt;msup&gt;&lt;mi&gt;P&lt;/mi&gt;&lt;mi&gt;H&lt;/mi&gt;&lt;/msup&gt;&lt;mi&gt;d&lt;/mi&gt;&lt;mo&gt;)&lt;/mo&gt;&lt;/mrow&gt;&lt;mo&gt;=&lt;/mo&gt;&lt;mrow&gt;&lt;mo&gt;(&lt;/mo&gt;&lt;msup&gt;&lt;mi&gt;P&lt;/mi&gt;&lt;mrow&gt;&lt;mo&gt;-&lt;/mo&gt;&lt;mn&gt;1&lt;/mn&gt;&lt;/mrow&gt;&lt;/msup&gt;&lt;msub&gt;&lt;mi&gt;d&lt;/mi&gt;&lt;mi&gt;maf&lt;/mi&gt;&lt;/msub&gt;&lt;mo&gt;)&lt;/mo&gt;&lt;/mrow&gt;&lt;/mrow&gt;</math></maths>上式中:<img file="C2004100528440003C1.GIF" wi="653" he="308" />P<sup>-1</sup>是一个NK×NK的矩阵;矩阵L为下三角矩阵,它是矩阵R<sub>0</sub>的乔莱斯基分解后的下三角矩阵:即满足R<sub>0</sub>=L×L<sup>H</sup>,其中L矩阵具有以下的结构:<img file="C2004100528440003C2.GIF" wi="586" he="305" />5.3:重写线性方程<maths num="0004"><math>&lt;mrow&gt;&lt;mi&gt;R&lt;/mi&gt;&lt;mover&gt;&lt;mi&gt;d&lt;/mi&gt;&lt;mo&gt;^&lt;/mo&gt;&lt;/mover&gt;&lt;mo&gt;=&lt;/mo&gt;&lt;msub&gt;&lt;mi&gt;d&lt;/mi&gt;&lt;mi&gt;maf&lt;/mi&gt;&lt;/msub&gt;&lt;/mrow&gt;</math></maths>预条件过程如下:<img file="C2004100528440003C4.GIF" wi="1814" he="302" />令:d′<sub>maf</sub>=P<sup>-1</sup>d<sub>maf</sub>,即有:<img file="C2004100528440003C5.GIF" wi="864" he="312" />则,继续简化得到:<img file="C2004100528440003C6.GIF" wi="1478" he="325" />5.4:由于R<sub>0</sub>=L×L<sup>H</sup>,上式可得:<img file="C2004100528440004C1.GIF" wi="1491" he="320" />其中I为单位矩阵;重写方程如下:<img file="C2004100528440004C2.GIF" wi="845" he="308" />该方程简写为:R′d′=d′<sub>maf</sub>式中令:<maths num="0005"><math>&lt;mrow&gt;&lt;msub&gt;&lt;mi&gt;V&lt;/mi&gt;&lt;mn&gt;1&lt;/mn&gt;&lt;/msub&gt;&lt;mo&gt;=&lt;/mo&gt;&lt;msup&gt;&lt;mi&gt;L&lt;/mi&gt;&lt;mrow&gt;&lt;mo&gt;-&lt;/mo&gt;&lt;mn&gt;1&lt;/mn&gt;&lt;/mrow&gt;&lt;/msup&gt;&lt;msub&gt;&lt;mi&gt;R&lt;/mi&gt;&lt;mn&gt;1&lt;/mn&gt;&lt;/msub&gt;&lt;msup&gt;&lt;mi&gt;L&lt;/mi&gt;&lt;msup&gt;&lt;mi&gt;H&lt;/mi&gt;&lt;mrow&gt;&lt;mo&gt;-&lt;/mo&gt;&lt;mn&gt;1&lt;/mn&gt;&lt;/mrow&gt;&lt;/msup&gt;&lt;/msup&gt;&lt;mo&gt;=&lt;/mo&gt;&lt;msup&gt;&lt;mrow&gt;&lt;mo&gt;(&lt;/mo&gt;&lt;msup&gt;&lt;mi&gt;L&lt;/mi&gt;&lt;mrow&gt;&lt;mo&gt;-&lt;/mo&gt;&lt;mn&gt;1&lt;/mn&gt;&lt;/mrow&gt;&lt;/msup&gt;&lt;msubsup&gt;&lt;mi&gt;R&lt;/mi&gt;&lt;mn&gt;1&lt;/mn&gt;&lt;mi&gt;H&lt;/mi&gt;&lt;/msubsup&gt;&lt;msup&gt;&lt;mi&gt;L&lt;/mi&gt;&lt;msup&gt;&lt;mi&gt;H&lt;/mi&gt;&lt;mrow&gt;&lt;mo&gt;-&lt;/mo&gt;&lt;mn&gt;1&lt;/mn&gt;&lt;/mrow&gt;&lt;/msup&gt;&lt;/msup&gt;&lt;mo&gt;)&lt;/mo&gt;&lt;/mrow&gt;&lt;mi&gt;H&lt;/mi&gt;&lt;/msup&gt;&lt;mo&gt;,&lt;/mo&gt;&lt;/mrow&gt;</math></maths><maths num="0006"><math>&lt;mrow&gt;&lt;msubsup&gt;&lt;mi&gt;V&lt;/mi&gt;&lt;mn&gt;1&lt;/mn&gt;&lt;mi&gt;H&lt;/mi&gt;&lt;/msubsup&gt;&lt;mo&gt;=&lt;/mo&gt;&lt;msup&gt;&lt;mi&gt;L&lt;/mi&gt;&lt;mrow&gt;&lt;mo&gt;-&lt;/mo&gt;&lt;mn&gt;1&lt;/mn&gt;&lt;/mrow&gt;&lt;/msup&gt;&lt;msubsup&gt;&lt;mi&gt;R&lt;/mi&gt;&lt;mn&gt;1&lt;/mn&gt;&lt;mi&gt;H&lt;/mi&gt;&lt;/msubsup&gt;&lt;msup&gt;&lt;mi&gt;L&lt;/mi&gt;&lt;msup&gt;&lt;mi&gt;H&lt;/mi&gt;&lt;mrow&gt;&lt;mo&gt;-&lt;/mo&gt;&lt;mn&gt;1&lt;/mn&gt;&lt;/mrow&gt;&lt;/msup&gt;&lt;/msup&gt;&lt;/mrow&gt;</math></maths>5.5上式中,V<sub>1</sub>与V<sub>1</sub><sup>H</sup>中,只需求解其中的一个矩阵即可;步骤6.多用户检测模块通过上述方程求解d′;所述的步骤6包含以下步骤:6.1:多用户检测模块进行初始参量的选择6.1.1:初始迭代值选取d<sub>0</sub>′;6.1.2:选择初始残量:r<sub>0</sub>=Rd′<sub>0</sub>-d′<sub>maf</sub>;6.1.3:收敛性检测,如果满足设定的收敛条件,则d′<sub>0</sub>就是所要求的解d′,退出运算过程;6.1.4:初始搜索方向的选择:方向q<sub>0</sub>=r<sub>0</sub>;6.2:进入迭代计算过程,以下角标k表示第k次迭代;k=1,2,3,...6.2.1:计算第k次迭代时的矩阵与矢量的乘积:Rq<sub>k-1</sub>;6.2.2:计算前进步长:α<sub>k-1</sub>=-(r<sub>k-1</sub>,r<sub>k-1</sub>)/(Rq<sub>k-1</sub>,q<sub>k-1</sub>),该步骤包括两个内积运算和一个除法运算;6.2.3:计算第k次迭代值:d′<sub>k</sub>=d′<sub>k-1</sub>+α<sub>k-1</sub>q<sub>k-1</sub>;6.2.4:更新残量:r<sub>k</sub>=r<sub>k-1</sub>+α<sub>k-1</sub>Rq<sub>k-1</sub>;6.2.5:收敛性检测,如果满足一定的收敛条件则停止计算,上述则d′<sub>k</sub>就是所要求的解d′,退出运算过程;跳出循环;6.2.6:计算参量,令λ<sub>k-1</sub>=(r<sub>k</sub>,r<sub>k</sub>)/(r<sub>k-1</sub>,r<sub>k-1</sub>);6.2.7:更新搜索方向:q<sub>k</sub>=r<sub>k-1</sub>+λ<sub>k-1</sub>q<sub>k-1</sub>;6.2.8:循环执行步骤6.2.1至6.2.7,直到步骤6.2.5的条件满足而停止计算跳出循环;步骤7.反预条件,对上述估计数据矢量d′进行反预处理,得到用户数据<img file="C2004100528440005C1.GIF" wi="50" he="65" />其中,反预条件是指:从d′解出<img file="C2004100528440005C2.GIF" wi="50" he="62" />根据<maths num="0007"><math>&lt;mrow&gt;&lt;msup&gt;&lt;mi&gt;P&lt;/mi&gt;&lt;mi&gt;H&lt;/mi&gt;&lt;/msup&gt;&lt;mover&gt;&lt;mi&gt;d&lt;/mi&gt;&lt;mo&gt;^&lt;/mo&gt;&lt;/mover&gt;&lt;mo&gt;=&lt;/mo&gt;&lt;msup&gt;&lt;mi&gt;d&lt;/mi&gt;&lt;mo&gt;&amp;prime;&lt;/mo&gt;&lt;/msup&gt;&lt;/mrow&gt;</math></maths>得到:<maths num="0008"><math>&lt;mrow&gt;&lt;mover&gt;&lt;mi&gt;d&lt;/mi&gt;&lt;mo&gt;^&lt;/mo&gt;&lt;/mover&gt;&lt;mo&gt;=&lt;/mo&gt;&lt;msup&gt;&lt;mrow&gt;&lt;mo&gt;(&lt;/mo&gt;&lt;msup&gt;&lt;mi&gt;P&lt;/mi&gt;&lt;mi&gt;H&lt;/mi&gt;&lt;/msup&gt;&lt;mo&gt;)&lt;/mo&gt;&lt;/mrow&gt;&lt;mrow&gt;&lt;mo&gt;-&lt;/mo&gt;&lt;mn&gt;1&lt;/mn&gt;&lt;/mrow&gt;&lt;/msup&gt;&lt;msup&gt;&lt;mi&gt;d&lt;/mi&gt;&lt;mo&gt;&amp;prime;&lt;/mo&gt;&lt;/msup&gt;&lt;/mrow&gt;</math></maths>即:<img file="C2004100528440005C5.GIF" wi="852" he="304" />
地址 201108上海市莘庄工业区申旺路18号