第32卷第11期 电子与信息学报 Vo1.32NO.11 2010年¨月 Journal of Electronics&Information Technology NOV.2010 基于分布式空频码的协作通信系统的信道估计及优化 张晴川① 束锋① 孙锦涛 王东明 f南京理工大学电子工程与光电技术学院 南京210094) (东南大学移动通信国家重点实验室 南京 210096) 摘要:该文研究了协作通信中分布式空频码系统的信道估计问题。为使信道可辨识,在中继节点处加入了循环卷 积滤波a(CCF)。推导了训练模式中导频序列和中继节点处CCF的准最优设计方案,并由此得到了源节点与中继 节点之间功率分配的闭合表达式。计算机仿真验证了本文工作的有效性。 关键词:协作通信;正交频分复用;信道估计;分布式系统 中图分类号:TN92 文献标识码:A 文章编号:1009—5896(2010)11—2565—06 DOI:10.3724/SP.J.1146.2009.01532 Channel Estimation and Optimization for Cooperative Communication Systems Based on Distributed Space—frequency Coding Zhang Qing—chuan① Shu Feng① Sun Jin—tao① Wang Dong—ming② (School of Electronic Engineering and Optoelectronic Technique, Nanjing University of Science and Technology,Nanjing 210094,China) (National Mobile Communications Research Laboratory,Southeast University,Nanjing 210096,China) Abstract:In this paper,the issue of channel estimation is investigated for distributed space—frequency coding system in cooperative communications.To make the channel being identifiable,Cyclic Convolution Filters(CCF) are utilized at the relay nodes.In the training stage,the suboptimal design is derived for both the pilot sequence and CCFs adopted at relay nodes.Furthermore,the closed expression of power allocation between source node and relay nodes is obtained.Finally,computer simulations demonstrate the effectiveness of the proposed approach. Key words:Cooperative communications;OFDM:Channel estimation;Distributed systems 1 引言 统中进行信道估计与传统点对点系统相比有很大的 协作通信旨在通过多用户间互相配合,用协作 不同。首先,不仅需要设计源节点的训练序列,还 的方式传输信息,来获取额外的分集以对抗信道衰 需要设计中继节点的信号发送方案;其次,总的信 落I1,2】。用户只需配备单天线就可获取与多天线系统 道是由多个信道合成的,这就意味着需要更大的训 相近的性能,协作通信这一巨大的优势使这方面的 练开销。另外,用户之间存在着独立的功率约束, 研究成为近年来的热点。 同时也有随之而来的功率分配fPower Allocation, 随着无线系统追求的传输速率越来越高,无线 PA1问题。 信道会呈现出明显的频率选择性。如果能同时获取 文献『5,6]讨论了协作通信中的信道估计问题, 频率分集与协作分集,必将大大提升系统性能。为 但只适用于平坦衰落信道。文献f7,8]针对的是基于 此,已有文献提出在正交频分复用fOrthogonal OFDM的协作系统,但限于两个用户节点的情况。 Frequency Division Multiplexing,OFDM1的基础 文献f91采用时域的方法讨论了多用户中继时频率选 上,用所谓的分布式空频码fDistributed Space— 择性信道的信道估计问题。而本文则根据空频码系 Frequency Coding,DSFC)来同时获取这两种分 统的特点从频域角度出发,讨论了专门针对协作 集[314】。绝大部分的协作系统都假设已知信道信息, DSFC系统的信道估计问题。本文提出在中继节点 这通常通过信道估计获得。而在协作式的DSFC系 处运用循环卷积滤波器来处理接收信号,这样保证 了多中继协作信道的可辨识。不同于上述文献,本 文全面的研究了相应的导频序列、中继端信号发送 2009—12—01收到,2010 04—28改回 方案以及包括在中继、目标用户处于不同信噪比条 国家自然科学基金(60702028)资助课题 件时的功率分配在内的一系列优化设计问题,提出 通信作者:张晴川schillerchang@hotmail.com 2566 电子与信息学报 第32卷 了一种准最优训练设计方案和具有闭合表达式的用 户间功率分配方案,并在能同时获取协作和多径分 集的DSFC编码系统中验证了本文所提出方法的优 越性。 本文中, ,O和。分别表示kronecker积、 Hadamard积以及循环卷积; 和F分别表示单位 矩阵和DFT变换矩阵;0,1和e“)分别表示全零 向量、全1向量和单位向量;(・) ,(.) 以及(・)一分 其中g。全E: ,DU处的加性高斯噪声t,均值为零, 方差为 。发送信号 是原始调制信号 经线性星 座预编码[3,1o1后的信号,z=@毫。由于囊是对角阵, 式f3)可进一步表示为 Y=E 2:1 。 @意+tl, (4) DSFC包含 以及@的设计[3I4],旨在同时获取协作 分集与多径分集。但前提条件是DU已知信道 别表示转置、共轭转置以及矩阵的逆; 和diag{a} 表示对角向量为a的对角阵;E(・)表示以下标为对 象求期望;tr(A)表示矩阵的迹;to 表示下取整操 作;l Ia Il和I a 1分别表示向量a的2范数和模向量; A, 表示矩阵A的第1到L行组成的矩阵,f A1 ,表 示矩阵第i行第J列上的元素。 2系统模型 假设1个信息源用户fSource User,SU), 个 中继用户(Relay User,RU)和1个信息目标用户 (Destination User,DU)组成一个协作系统,用户间 信道为有频率选择性的瑞利衰落信道,每一个用户 终端都只配备单天线。采用子载波数为Ⅳ的OFDM 系统,假设用户问多径信道的在采样间隔上的可分 辨路径数都为 。SU到RU的信道冲击响应为 {元 ) ,…, ,RU到DU的信道为{ ) , ,它们 的元素均是方差为1均值为0的独立的高斯变量。分 布式空频码协作系统工作分为两个阶段。第1阶段, SU发送RU接收。将SU的频域发送信号用向量 来表示,则RU处的接收信号为 =√ + (1) 其中i=1,…,M;E(1 lz I12)=1,尸0表示SU发送 一个OFDM符号的能量;hi全 无 ,表示SU与 第i个RU之间的信道的子载波增益;高斯白噪声向 量 均值为0,方差为 。接着在第2阶段中,RU 发送DU接收。第i个中继节点在频域上的发送信号 为 = (2) 其中 是RU的频域线性变换矩阵,通常取其为对角 阵,O/是相应的功率常数。假设第i个RU的功率约 束为 ,则唾的对角向量可取【1,ej t/ ,…, e伽(N-1) /Ⅳ】【3],相应地, =√ /P0+Ⅳ ;或者, 对角向量取[1T+l 0 e 1 ] :Jv,J=lN/MLJI4],相 应有 =,/MP,/( +JMLa ̄)。DUG接收信号为 M Y=E Xi+t, i=1 M M =Ea p +E 唾 +t, i=1 0=1 { o ) ,M,而这通常需要通过信道估计来获 取。 3信道估计 本节假设系统工作于训练模式,此时SU发送 导频序列8,除此之外,为表达方便,其余变量均 沿用上节中的含义。正如本节下文中将指出的,为 使信道可辨识,每一个中继端将接收信号去CP后, 都需要通过一个循环卷积滤波器fCyclic Convolution Filter,CCF)。然后,将滤波后的信号 加上CP再发送。将NX1维的时域CCF向量用 {专) .一, 表示,对应的频域CCF向量毫=F毛。 中继端在频域上的发送信号为 =p Y 。类似于上 节的分析,系统工作在训练模式时的模型为 Y=E 。 s+∑ +t, (5) 为便于分析,在式f5)中,我们将用户节点的功率约 束包含在s以及岛之中。在训练模式中,SU的功率 约束为芦o。本文假设各中继节点处的噪声方差相等 (等于 ),因此鉴于信道间的这种对称性,假设M 个RU有相同的功率限制 。本文的目的就是估计 DSFC系统的合成信道{g o ) 一.M,同时满足以 下功率约束 Il 8 P0 I … ‰ ( , :1,…,M『 J 频域合成信道与时域信道响应间的关系为 g。ohi=F(雪。,0 0 0) (7) 其中雪 和元 分别表示 。和元 补零后的N X 1向 量。令于 为向量豆油④无 ,0去掉尾部的零元素后的 全2 一1维向量。现在问题转化为估计时域合成 信道{于 ) ,,M。由式(7),式(5)可化为 Y=I—D,,D ̄F1:Ⅳ… I— 一日: 1I+'112, (8) 垒于 合成噪声向量W的协方差矩阵为 第11期 张晴川等:基于分布式空频码的协作通信系统的信道估计及优化 2567 A垒 z} l( "H)=diag ∑f岛f。+oral f(9) \ i=1 , 束下,式(11)和式(13)不仅退化为同一个待优化的代 价函数,而且由于A与 成比例,其表达式也得以 简化,在此基础上可以推出最优设计条件并获得具 不失一般地,本文假设N>KM。考查式(8),如果 RU处不做循环卷积处理,则N×KM维矩阵P的秩 体的设计方案。将式(14)的约束条件代入式(9)和式 为 ,是一个奇异矩阵,信道向量于无法辨识。而 (11),再联合式(6)的功率约束,8和{岛) ,…, 的最 我们在RU处引入CCF,就可以使P列满秩,这就 解释了我们在中继端使用CCF的原因。对于式(8), 于的最小二乘(Least Square,LS)估计为 ,Ls=RLsY=(P“P) P“Y (10) 于LS的均方误差(Mean Square Error,MSE)为 E(1l于 s一于lI2)=tr(RLsARLHS) =tr{(P“P) P“AP(P P) )(11) 式(11)中,期望运算的对象与式(9)中的相同。考虑 到噪声协方差矩阵A通常并不与 成比例,采用加 权最小二乘(Weighted Least Square,WLS)估计可 以提高性能 。WLS估计的表达式为 于wLs=RwLsY=(P“ P) P“ Y (12) 而于w 的均方误差为 E(1 l于w 一于l 12)=tr(Rw AR昌 ) =tr{(P“ P) } (13) WLS的性能通常都优于LS,除非A与 成比例。 虽然WLS性能更好,但WLS估计需要知道DU与 RU处具体的噪声方差,而LS并不需要,因此LS 估计的实用性更佳。为了使在噪声统计特性未知的 情况下,使用LS估计较之WLS不会有性能的损失, 本文提出给CCF附加如下的约束条件: } tl (14) 其中t是一个可调参数。对CCF附加式(14)的约束 后,将使得DU在每个子载波上的接收信号在统计 平均意义上具有相同的增益幅度。显然,现在于的 LS估计与WLS估计两者等价。此时LS有与WLS 完全相同的性能。附加式(141的意义还不仅在于此, 在噪声方差已知时,我们希望优化信道估计的性能, 但该性能优化问题是非常复杂的非线性优化问题, 然而有了式(14),问题将得以简化。这将在下一节 中讨论。 4优化设计 在DSFC系统中,为使信道估计最优,需要优 化设计导频序列s和CCF向量{岛} M。信道估 计的性能度量是估计结果的MSE。对于LS和WLS, 其MSE分别是式(11)和式(13)。然而,在式(6)的约 束下,直接按式(11)或式(13)优化8和{专) 一,M是 非常困难的。但是引入了式(14),在这个更紧的约 优设计可归结为如下优化问题: {(!G2Rt ̄G D2)-I-DIHe) ) ( 5) .v s.t.∑l s。1。 o (16a) i=1 Ⅳ ∑ ,l。+ ,l。 户 ,i:1,…,M(16b) j=l M ∑ l 一t:0,J=1,…,N (16c) =1 其中将s与岛表示为s=[s¨…,sⅣ 以及岛= , …,毛.Ⅳ] ;同时须注意,待优化的变量还包括可调 参数t。本文假设协作系统中有一个中央控制单元, 控制着用户问的能量分配。在这种情况下,信道估 计有一个总功率约束户。式(16a)和式(16b)是分布式 系统自然形成的功率约束,而本文引入的式(16c)将 使得该优化问题具有形式简单的闭合解,从而大大 提升了本文提出的优化设计方案的实用性。将式(151 中的矩阵 展开表示为 B=( £+ 刍) 硪 l曩: …磁 『。p : 砝D H吖 E: …堪D H qs1 E. (17) 对于KM×KM维的正定矩阵B的逆矩阵的迹,有 如下关系 : tr(B ) ∑ (18) 式(18)中的等号仅当B是对角阵时成立。求解目标 函数式(15)在约束条件式(16)下的优化问题的思路 是:首先求出∑ KM[曰] 的最小值并得到达到此最 小值的条件(简称为条件1);再根据式(18)等号成立 的条件,得到条件2;如果能找到同时满足条件1 和条件2以及约束条件式(16)的设计方案f即当目标 函数的下界∑ KM[B坛 达到其最小值时,式(18)的 等号也能同时成立),那么该方案必是式(15),式(16) 原始优化问题的最优解。考查式(17),有 K∑ ==M1 ∑f=1 I ∑ i=1 u l uD 广cJf ()19, 电子与信息学报 第32卷 同时,将约束条件式(16b)的M个不等式相加,得到 M N ∑∑I毫,J l。I sj f。+盯 J l。 i=1 j=1 .Ⅳ =同时不难得出,{岛) ,…,M的所有分量也是等功率 的。于是,结合s和{毫) 一,M的恒模特性OAR条件 1的式(23)和式(24),可得到一种满足条件的导频序 列和CCF设计方案为 :.∑t I sj I。+Ⅳ j=l 户l (20) ,其中我们利用了约束条件式(16e)。于是可以得到 V厩 -v I- oj 丽ej2 ̄'(i-1 1一 ’ ∑∑ J I。I勺I。 i=1 j=1 .v 虽然式(26)中的.[ ) ..,Ⅳ可为任意值,但我们设置 它为一个预先确定的伪随机序列,使得SU的发送 =∑t I I。 min{tPo,M 一lv ̄t}(21) 进一步根据式(19),式(21)以及不等式∑ ÷ 。,‘/ 朋 0 ,可以推出 —— t=l。 KM 黠‰(22 式(22)右边仅仅是t的函数,显然,当tPo=MPl 一Ⅳ4t,即 t::丝 f231 0+ V 五 此时式(、 22)的右式达到最小值,这就是 =I fB … 的最小值。而达到此最小值需要式(21)和式(22)的等 号同时成立,这要求满足下面的条件 Ⅳ ∑l瓯l =po i=1 N N (24) ∑ I。I s』I。=∑ l。I l。,Vi,J∈ 一, f=1 f=1 将式(23),式(24)联合,即得到条件1。接着需要找 到条件2。而tr(B )要达到其下界,即式(18)中的 等号成立,则要求B是对角阵。同时,条件1中的 式(24)意味着B的对角元都相等,等于t /M。于 是,有 f£芦o, . . 堪O H。 曩:Ⅳ={【 0 (25) , i≠J ’式(25)即条件2。可以证明,式(25)等价于两个条 件[12]:各个序列的分量有相等的模;序列间满足相 位偏移正交性。 接下来,需要找到符合以上两个条件的设计方 案。根据文献[12】,满足条件2的相位偏移正交性的 一种设计方法为s 毛j=,√t /NMej2r(i-1) (j-1)/Ⅳ。同 时,条件2要求.[f 毛,J j都相等,必然有 ∑ I 8i l。I岛, J。=∑ l I。l岛,J 2 Vi,J。再根据式 f16c),可以得出s的分量是恒模的,也即等功率的, 信号保持较低的峰均功率比。容易验证,式(26)满 足原始约束条件式(16)。至此,我们得到了式(15)、 式(161中约束优化问题的最优化解。 分布式协作系统中,SU和RU用户间存在着总 功率约束下的功率分配问题。定义Co( o,芦1)全 min{tr(B ))。将式(23)代入式(22)的右边,可以得 到满足最优条件时,则tr(B )的最小值Co(Po,P1) 为 Co(Po, ): (27) 考虑总功率约束 0+MP1 ,令f(Po,芦1) 芦0+M芦1~ ,根据Karush-Kuhn-Tucker条件, 最优功率分配 和 满足 壶Co(P ,P 1)+恃 _0)i-0 }(28) Af(P0,P1)=0, 0 I限于篇幅,我们略去了上式的求解过程,直接给出 以下结果。 P~ P/2M,P~o:P/2,P=1 P~ l= , = —MP~ ;,p喾 (29) 其中Q全√ +N(o-2+ )芦+N ,P全 / 。有趣的是,根据式(29),当SU到RU间信道 和RU到DU间信道的噪声功率相同时,最佳功率 分配非常简单,就是SU分配总功率的一半,其余 平均分配给各个RU。需要指出的是,本节的结论 是在附加了式(14)的约束后得到的,因此严格意义 上不一定是全局最优的,本文称其为准最优的设计。 5仿真结果 本节中的仿真参数设定为:OFDM的子载波数 目Ⅳ:64,CP长度为10,RU节点数M=2,4, 各用户节点间的多径信道的抽头数L=8,其功率 延迟分布服从指数关系 ~,?=0,…,三一1,其中归 一化常数 =1/ ? e~。假设用户间信道在一个 第ll期 张晴川等:基于分布式 频码的协作通信系统的信道估计及优化 2569 OFDM符号周期内保持不变。每个OFDM符号时 间内,所有用户节点间具有总功率约束 。传输信 息数据时,每个用户节点分配相同的功率。信噪比 小,这是由于式(29)的取值逐渐接近于等功率分配。 需要指出的是,式f29)表明,除了P=1的特殊情况, SU和RU间的功率分配并不只取决于p,与噪声方 差的绝对大小也相关。 的定义为芦/N4。假定在每一帧中,导频序列占用 一个OFDM符号。为了比较的公平性,定义信道系 最后,图3是P=1时,不同训练方案问的误符 数的平均均方误差(Average MSE,AMSE)为 E(1I于.一.尹lf)/KM。选择如下3种方案与本文提 出的式(26)的设计方案进行对比:(S1)8为高斯序 号率(Symbol Error Ratio,SER)性能比较。我们取 DSFC对信息数据的预编码分组长度【 】为8,此时的 协作系统可获得的最大分集数为8。当M=2,在 列,{ ) ,… M是高斯随机向量,仅满足约束条件式 (6),不满足式(14),采用WLS信道估计器;(S2)s SER=10 时,使用最优功率分配式(29)的准最优 设计式(261比采用理想信道知识时的性能仅差1 dB 左右,而比随机训练至少有2 dB的性能提升。当 为高斯序列,{ ) 、,M是等功率的随机序列,满足 式(6)和式(14),即式(16);(s3)s和{岛) ,. 都是 等功率的随机序列,满足式(6)和式(14)。 首先,比较使用不同导频和CCF时的信道估计 性能。假设P一1,且在训练模式中,所有的SU和 4时,这一性能增益变得非常巨大。非最优的 训练方案的性能变得非常差,而准最优设计则始终 =表现优异。当P取其它值时,本文提出的设计方案 的SER性能也明显优于其它的对比方案,且与使用 RU节点都分配相同的功率。图1是AMSE的仿真 曲线。本文提出的式(26)的设计方案优于其它满足 式(161约束的随机训练方案,这符合第4节中的分 理想信道知识时较为接近,限于篇幅,这里略去了 仿真图。 6结论 本文讨论了分布式空频码系统的信道估计问 析。而且不考虑附加约束式(16)的(S1)性能仅相当于 (S21,明显比本文提出的设计方案差。图1说明了 式f261的设计具有较高的信道估计精度。尤其当 =4时,待估计参数的数量较大,该方案的优势 题。在中继端对接收信号做循环卷积滤波可以保证 信道可辨识。通过引入一个合理的附加约束,不仅 可以使得常用的LS信道估计器避免了性能损失。而 且使原本高度复杂的信道估计性能优化问题变得可 解,从而得到了准最优的训练设计准则以及具体的 实现方案。在此基础上,推导了功率分配的解析表 达。本文提出方法的估计精度高,性能稳健,在与 多个有代表性的随机训练方案的比较中,体现了明 显的优势,并且其性能与使用理想信道知识时较为 接近。 非常明显。 其次,在准最优设计式(26)的前提下,比较使 用式(29)的最优功率分配与使用等功率分 ̄(EPA1 时的性能差异。AMSE曲线如图2所示。其中包含 了P分别取0.1,1以及10三种情况。当P=0.1时, DU处的加性噪声较弱,最优功率分配有很大的 MSE性能增益。随着P的增大,性能的差异逐渐减 图1 M=2,4时,等功率分配下, 图2 P取不同值时,等功率分配与 图3 M=2,4时,不同的功率分配 不同训练设计方案的AMSE性能曲线 最优功率分配的AMSE性能比较 以及训练没计对系统SER性能的影响 参考文献 [1】 Sendonaris A,Erkip E,and Aazhang B.User cooperation diversity--part I,II[J].IEEE Transactions On Communications, 2003,51(11):1927—1948. [2】Laneman J N,Tse D N C,and Wornell G W.Cooperative diversity in wireless networks:eficient protocolfs and outage behavior[J].IEEE Transactions On Information Theory,2004, 50(12):3062—3080. [3】 Zhang W,Li Y,and Xia X,et a1..Distributed space- requency codifng for cooperative diversity in broadband 2570 电子与信息学报 第32卷 wireless Ad hoc networks[J].IEEE Transactions on Wireless Communications,2008,7(3):995—1003. Seddik K G and Liu K J R.Distributed space-frequency channel estimation method for OFDM—based cooperative diversity systems with multiple amplify-and-forward relays[J] EURASIP Joureat on Wireless Communication and coding over broadband relay channels[J].IEEE Transactions on Wireless Communications,2008,7(11):4748—4759. Patel C S and Stuber G L.Channel estimation for amplify Networking,http://downloads.hindawi.com/journals/wcn/ 2008/149803.pdf,2008,9 pages. u Z,Xin Y,and Giannakis G B.Linear constellation 10] Liprecoding for OFDM with maximum multipath diversity and and forward relay based cooperation diversity systems[J]. IEEE Transactions on Wielress Communications,2007,6(6): coding gains[J].IEEE Transactions on Communications, 2348—2356. Gao F,Cui T,and Nallanathan A.On channel estimation and optimal training design for amplify and forward relay networks[J1. IEEE Transactions on Wireless Communications,2008,7(5):1907—1916 Kim K,Kim H,and Park H.OFDM channel estimation for the amp and-forward cooperative channel[C].The 65th IEEE Vehicular Technology Conference,Dublin,2007: 1642—1646. Dogan H.Maximmn a posteriori channel estimation for cooperative diversity orthogonal frequency division multiplexing systems in amplify-and—forward mode[J[./ET Communications,2009,3(4):501—511. Yah K,Ding S,and Qiu Y,et a1..A low-complexity LMMSE 2003,51(3):416—427, 【n】 Kay S M.Fundumentals of Statistical Signal Processing: Estimation Theory[M].Englewood Clifs,NJ,Prentice-Hall, 1993:225—226. [121 Barhumi I,Leus G,and Marc M.Optimal training design for MIMO OFDM systems in mobile wireless channels[J].IEEE rTansactions on Signal Processing,2003,51(6):1615—1624. 张晴川: 男 1981年生,博士生,研究方向为无线通信中的信号 处理 、协作通信等. 柬锋 男。 1973年生,博士,副教授,研究方向为无线通信技 术、 卫星定位等. 孙锦涛 男, 1942年生,教授,博士生导师,研究方向为现代通 信理 论、通信信号处理、雷达信号处理等.