一种适用于快时变稀疏估计的最优导频图样的设计方法
【专利摘要】本发明提出一种适用于快时变稀疏估计的最优导频图样的设计方法，具体为：越来越多的测试实验表明无线信道具有稀疏性，而目前导频图样设计主要基于密集多径的传统信道估计方法，需要大量的导频信号来获取准确的信道状态信息，造成频谱资源利用率低；另一方面，目前稀疏信道估计方法大多适用于慢时变或时不变信道，对信道稀疏性的利用只局限于时延域而没有扩展到多普勒域，并且忽略了由多普勒频移产生的载波间干扰对信道造成的影响，因此在快时变信道中采用稀疏信道估计方法的精确性和可靠性会大幅下降。基于此，提出一种适用于快时变稀疏信道估计的最优导频图样设计方法，所提方法具有更好的鲁棒性和更高的信道估计精度。
【专利说明】
-种适用于快时变稀疏估计的最优导频图样的设计方法
技术领域
[0001] 本发明主要设及快时变信道中的信道估计领域，尤其设及基于导频图样设计的稀 疏信道估计方法领域。
【背景技术】
[0002] 随着无线通信系统的发展，特别是移动通信技术的不断提高，对数据传输速率、通 信系统性能和频带宽度提出了更高的要求，由于传统基于时分或频分的调制技术空间可用 频谱资源有限，高速传输时系统的通信质量会出现波动，因此迫切需要更先进的技术来进 一步提高频率资源利用率，W满足大容量、高速率的业务需求，同时还要克服高速数据在无 线信道下的多径衰落，降低多径干扰和多普勒频移造成对系统的影响，从而改进系统性能。 在运种需求推动下，正交频分复用（〇：rthogonal Frequenc^y Division Multiplexing, (FDM)和多天线多入多出(Multiple I吨ut Multiple Output,ΜΙΜΟ)逐渐成为研究热点，并 在移动通信中得到了越来越广泛的应用。但要想通过OFDM技术和ΜΙΜΟ技术实现系统性能的 提高，相关核屯、技术的开展必不可少，而信道估计技术就是其中的难点之一。
[0003] 于是，本发明主要基于信道估计领域展开研究。根据调研发现，为了获得移动信道 的频率响应，最常采用的方法是基于导频辅助的信道估计方法，而其中导频的时、频域间隔 和导频位置的不同成为影响整个系统性能的重要因素，所W有关导频图案的设计至关重 要。根据已有的研究基础发现，目前大部分的导频图样设计均是针对传统的信道估计方法 进行的，而针对稀疏信道估计的导频图样设计方法相对较少。但是，随着研究的深入，越来 越多的测试实验表明无线信道具有稀疏特性，即约少于10%的多径信道占据着信道85% W 上的能量，而传统的信道估计方法均假设无线信道是密集多径的，从而利用大量的导频信 号来获取准确的信道状态信息（化annel State Information，CSI)，导致频谱资源利用率 低。所W研究针对稀疏信道估计的导频图样设计方法具有重大的意义。本发明也旨在通过 设计最优的导频图样、利用更少的导频数估计出更优的信道频响特性，从而提高信道估计 精度和频谱资源利用率。
[0004] 另一方面，虽然国内外已经逐步在开展基于稀疏信道估计方法的研究和应用，并 也取得了一些实质性的成果。但是，对于该领域的研究仍然存在着很多问题仍未解决。根据 已有的研究基础发现，目前的研究大多是针对慢时变或者时不变信道，而运些方法对信道 稀疏性的利用只局限于时延域而没有扩展到多普勒域。同时，由于忽略了由多普勒频移产 生的载波间干扰（Inter Carrier Interference，ICI)对信道的影响，所W在快时变信道中 采用稀疏信道估计方法的精确性和可靠性会大幅下降。然而，随着高速铁路和高速公路的 开通和应用，未来移动通信系统往往面临高速移动环境，而在高速移动环境下，无线通信信 道往往是快时变信道。若不能适应运种变化，那么通信系统性能将会受到严重影响，从而极 大地降低了信息传输的速度和质量。所W研究适用于快时变信道的稀疏信道估计方法意义 重大。
[0005] 基于此，本发明针对快时变信道，提出了一种适用于快时变稀疏估计的最优导频 图样的设计方法，可w具有更好的鲁棒性和更高的信道估计精度。
【发明内容】
[0006] 发明目的:在快时变信道中，针对稀疏信道估计方法，提出一种最优的导频图样设 计方案。
[0007] 本发明的技术方案：
[000引为方便阐述，W单输入单输出-(FDM(Single I吨ut Single Output-0抑M，SIS0- OFDM)系统为例阐述问题，其中关于MIMO-0抑Μ的系统同理。在快时变信道中，假设有N个子 载波，通信系统的传输模型为：
[0009]
扣
[0010] 其中，x(t)是发射端信号，y(t)是接收端信号，z(t)是高斯白噪声，h(t，T)是信道 的时域冲击响应，可表示为：
[0011]
(2)
[001^ 其中，S为不为0的抽头数，as，Ts，fs分别是第s(s = l，...，S)条传播路径的衰落系 数、时延和多普勒频移。对于快时变信道来说，ICI的影响不能被忽略，所W其对应的离散模 型为：
[0013]
㈱
[0014] 其中，T为采样时间，L为信道最大的多径时延，1=0,. . .，kl，n = 0,. . .，N-1，Ts/T <L，S<L。
[0015] 则式(1)对应的离散形式如下所示：
[0016]
(4)
[0017] 可写为矩阵形式，得到：
[001 引 y 二 hx+Z (5)
[0019] 其中，y=(y[0],... ,y[N-]J)T，x=(x[0], ,x[N-：L])T，z=(z[0], ,Z[N-1] )T，h为信道矩阵，
[0020]
[0021 ] 令X表示发射端的数据信号，经快速傅里叶逆变换（Inverse Fast Fourier Transform，IFFT)后转换为发射信号X，而在接收端，接收信号y经快速傅里叶变换（Fast Fourier Transform，FFT)得到接收端的数据信号Y。矩阵Q表示离散傅里叶变换（Discrete 化11心1· T'ransform'DFT)矩阵，其中0(尸,(?）二1/^ -心,Ρ和q分别代表DFT矩阵的行 和列，l《p，q《N，则有:x = QHx，Y = Qy，高斯白噪声矩阵z =化，系统矩阵H =化班。假设系统 的循环前缀长度大于最大多径时延L，那么符号间干扰(Inter Symbol Inte计erenceJSI) 干扰可W近似忽略。则可得到完整的系统频域传输模型如下：
[0022] Υ =册+Z (6)
[0023] 为了充分的利用快时变信道模型的时频稀疏特性，将式（3)转换到延迟-多普勒 域，得到信道的基扩展模型函数，即：
[0029] 其中Φ3(1)和恥（d)分别描述了信道的时间延迟和多普勒频移，而(l)s(l)和扣(d) 的能量分别集中于W散射点Ts/T和fsNT为中点的邻域。在常见的移动通信系统中，多普勒频 移通常只局限于一个很小的范围，即|fsNT|<<l，所W对于s=l，...，S，托(d)的能量始终 集中于零点附近，即每个子载波的能量泄露主要只影响其相邻的几个子载波，则可W采用 带状近似矩阵。假设每个子载波的频率偏移最大的距离为D，则信道矩阵Η被近似为：
[0030]
Π 0)
[0031] 该带状近似矩阵只保留了 Η的对角线W及对角线两边距离D之间的元素。若ICI的 带宽D远小于载波数Ν，则近似模型的复杂度将大大降低。即假设恥(d)的能量主要集中于[- D，D]区间，则信道基扩展模型函数Us[d，l]W及u[d，l]的支撑集由[0，L-l]X[0，N-l]缩小 为[0，レリX[0，D]U[N-D，N-l]，其中[N-D，N-リ由[-D，-リ循环移位所得。因此，需要估计 的信道参数数量由化减少为(2D+1) L。
[0032] 而关于近似模型ICI带宽D的取值至关重要，其取值大小会影响整个估计方法的性 能。当D取值较大时，近似模型更精确，更接近于原始的信道矩阵，但同时也会带来更多的待 估计量，即（2D+1化越大，则大大增加了估计的复杂度，从而降低了后续稀疏重建算法的性 能。当D取值较小时，由近似操作带来的误差将会增大，但待估计量能够有效降低。
[0033] 为了更方便地估计信道参数，将式(6)转换为W下形式：
[0034] Y = Ah+Z (11)其中，信道冲击响应向 量h定义为：
[0038] 其中，k和q分别代表矩阵的行和列，n=(q-l)mod(N)，l = (q-1-n)/N。将待估计目 标h转换为信道的基扩展模型，而又由于u[d，l]的支撑集为[0，L-1]X[0，D]U[N-D，N-1]， 因此定义(2D+〇L维向量U有：
[0039]
[0040] 由于U具有稀疏特性，所W结合式(7)(12)(14),可得到如下关系：
[0041]
0巧
[0042] 其中，II为L阶单位矩阵，F为(2D+1)XN维的单位选择矩阵，表示从[1，N]抽取子集 [1，0川^-0+1，則的采样操作，9为0。1'矩阵。用11代替}1进行求解，则将式（11)转换到延迟- 多普勒域，得到：
[0043] Υ=巫 U+Z (16)
[0044] 其中，Φ称为测量矩阵，并有如下形式：
[0045]
(17)
[0046] 式（16)可W看作是求解稀疏信号的线性方程。根据式（13)(17)可推导出矩阵Φ中 的元素具有如下表达式：
[0047]
狂 8)
[004引其中，n=(mq-l)mod(N)，1 = (mq-1-n)/N，运里（2D+OL维向量m定义为：
[0049] 、0，，… (19) 三巧化 + \，...，Nk + 0斗1，械[+ N - D+l,.-.，Nk + N)
[0050] 由式（18)可知，在矩阵φ的第k行，〇k,q的组成成分包含发送符号X(k-n-l)mod(N)+l。当 Π _ 0时，X(k-n-l)mod(N)+l _Xk ;而当Π 声 0时，X(k-n-l)mod(N)+l为与Xk相邻的·f守亏，即X(k-n-l)mod(N)+l E (Xk-D,. . .，Xk-i，Xk+i,. .，Xk+D}。因此，第k个观测值化的取值不仅取决于Xk，还受到Xk的长度为 (2D+1)的邻域内其他符号的影响，而运些相邻符号对化的作用即等价于ICI干扰。
[0051] 假设子载波N中有P个导频符号，即Xn(i), . . . ,Χν(ρ)，则有：
[0 化 2]
(20)
[0化3] 可简写为：
[0化4] 否=遣11 +爱 (21)
[0化5]由于在Φ的第P行中，由于观测值Yn(p)受到与导频符号Xn(p)相邻的数据符号的干 扰，矩阵仍然存在部分未知元素，无法直接得到。于是，假设将Φ中包含的未知元素全部置 零，得到Φ的近似矩阵Φ'。设导频子载波集合为Ω =化(l)，k(2)，. . .，k(P)}，则史'可表示 为：
[0化6]
(22)
[0057]其中，n= (mq-l)mod(N)，且k(p) e Ω。由于只存在由导频符号构成的点，于是壶'完 全已知。假设W =(憂-Φ')ιι +玄，用W替代噪声Z，得到近似的测量方程：
[005引 ?=Φ'υ 十W (23)
[0化9] 因为Φ.·满足有限等距性质（Restricted IsometiT P;rope;rty，RIP)性质，之后即可 采用压缩感知（Compressive Sensing，CS)重构方法求解待估计量U，转而通过式（15)求解 信道参数h。
[0060] 采用运种近似处理方法会引入两种误差影响估计方法的性能：一是由式(22)的近 似处理引来的误差(憂-Φ')?ι，等价于ICI干扰;二是CS的重构误差。
[0061] 针对上述第一种误差，应降低ICI干扰对估计方法的影响，即降低导频符号所受到 的邻近的数据符号的干扰。因此，最优的导频排列方案是将所有导频符号集中放置在一起， 即连续导频图案。对于矩阵Φ而言，可近似认为其所有元素都只包含导频符号。因此， Φ?添，则由ICI造成的误差可W被忽略。
[0062] 针对上述第二种误差，应提高稀疏信号重建的精度和稳定性。而CS算法重建性能 主要由测量矩阵的RIP决定。与其他矩阵相比，随机矩阵具有最高的概率满足RIP性质。因 此，现有的稀疏信道估计方法，一般采用完全随机的导频图案，即将所有导频符号随机地插 入到数据符号中，W保证CS重建算法具有相对稳定的性能。但是，当导频码随机排列时，其 很可能与数据符号相邻，导致误差项(Φ.-Φ)Η增大。另一方面，连续导频排列方案也无法保 证感知矩阵Φ能满足RIP性质，可能导致重构误差增大。
[0063] 因此，为了能同时有效地降低ICI干扰和重建误差，提出采用随机分组导频图案。 即将导频符号分为若干个相同长度的组，再将运些导频组W随机的位置地插入到数据符号 中。当导频组的长度化足够大时，可W认为每个组中只有首尾的导频符号会受到数据符号 的干扰。因此，尽管该方案不能完全剔除误差项(Φ--Φ)α，却可W大大降低该误差对估计结 果的影响，而且对ICI的鲁棒性强于随机导频图案。同时，对比连续导频图案，该随机分组导 频图案由于具有更强的随机性，能够保证感知矩阵能W更高的概率满足RIP性质，使CS重建 算法的性能更优。
[0064] 由于目前导频所在的位置W及分组长度化具有不确定性，直接证明其RIP性质W 及计算RIP常数难度较大。因此，可采用矩阵的相干系数替代RIP常数来评价导频图案对于 CS重建算法的重建性能的影响。
[0065] 矩阵Φ的相干系数μ(Φ)定义为：
[0066]
(24)
[0067] 其中，Φ i和Φ j都是Φ的列向量，而μ( Φ )取值最小时，估计结果具有最低的误差， 而CS算法具有最好的重建性能。
[0068] 综上，本发明基于W上背景，提出了一种适用于快时变稀疏估计的最优导频图样 的设计方法。具体流程如下：
[0069] 步骤 1、初始化d = 〇，Pg=1;
[0070] 步骤2、建立信道模型方程式，根据(24)计算μ( Φ );
[0071] 步骤3、令Pg值保持不变，D = D+1，在0，^区间内逐步增大D值，寻找能使其μ(Φ) 最小的D值并保存D值；
[0072] 步骤4、保持寻找到最优的0值不变，令口(； = 口〇+1，在化^区间内逐步增大口〇值，寻 找能使其μ( Φ)最小的Pc值，并保存Pc值；
[0073] 步骤5、判断P/化是否为整数，若是，进入步骤6;若否，则令P/化=[P/Pg]職，其中 [·]取整表示按照四舍五入法进行取整，此时则存在其中某一个导频组比其他导频组多或者 少一个导频数的情况，进入步骤6;
[0074] 步骤6、判断P/Pe是否小于阔值ξ，若是，进入步骤7,若否，跳转步骤10;
[0075] 步骤7、随机的选择P/Pe个导频组信道得到原始的导频信道图案；
[0076] 步骤8、随机改变其中一个导频组位置(不再选择已最优处理过的导频组），重新计 算μ( Φ )，直至寻找到拥有最小的μ( Φ )的导频组位置，保存运个导频组位置；
[0077] 步骤9、重复步骤8第Ρ/化次，直到寻找到所有P/Pg个导频组的位置，保存，即寻找到 最优的导频图样；
[007引步骤10、将N个子信道均匀的分为[P/Pg/K]職个块空间，其中P/化是导频组信道个 数，而K为算法复杂度的一个参数，一般设为1-3;
[0079] 步骤11、在每个块区间中，随机的选择K个子信道作为导频组信道，从而得到原始 的导频信道图案；
[0080] 步骤12、随机选择一个块区间；
[0081] 步骤13、在区间内，随机改变其中一个导频组位置(不再选择已最优处理过的导频 组），重新计算μ( Φ )，直至寻找到拥有最小的μ( Φ )的导频组位置，保存运个导频组位置；
[0082] 步骤14、重复步骤13第Κ次，直到寻找到运个区间内所有Κ个导频组的最优位置，保 存；
[0083] 步骤15、随机选择一个新的块区间，重新执行步骤13、步骤14,寻找并保存运个区 间内Κ个导频组的最优位置；
[0084] 步骤16、重复步骤15第[P/Pg/K]??次，直至寻找到所有导频组的位置，保存并整 合，即寻找到最优的导频图样。
[0085] 关于判断导频图样是否采用分块处理，通过阔值ξ来决定，而一般设定阔值ξ>10。
[0086] 如上所述，本发明结合快时变信道，创新性地研究了针对稀疏信道估计的导频图 样设计方法，通过利用随机导频分组方案设计出最优的导频图样，可得到更优的信道频响 特性，从而提高了信道估计精度和频谱资源利用率。此外，本方法不仅可W适用于SISO- OFDM系统模型，也可适用于MIMO-OFDM等系统模型。
【附图说明】
[0087] 本发明的上述和/或附加的方面和优点结合下面附图对实施例的描述中将变得明 显和容易理解，其中：
[0088] 图1适用于快时变稀疏估计的最优导频图样的设计流程。
【具体实施方式】
[0089] 下面详细描述本发明的实施例，所述实施例的示例在附图中示出，其中自始至终 相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附 图描述的实施例是示例性的，仅用于解释本发明，而不能理解为对本发明的限制。
[0090] 在本发明的描述中，需要理解的是，术语"纵向"、"横向"、"上"、"下"、"前"、"后'、 "左"、"右"、"竖直"、"水平"、"顶"、"底"、"内"、"外"等指示的方位或位置关系为基于附图所 示的方位或位置关系，仅是为了便于描述本发明和简化描述，而不是指示或暗示所指的装 置或元件必须具有特定的方位、W特定的方位构造和操作，因此不能理解为对本发明的限 制。
[0091] 在本发明的描述中，除非另有规定和限定，需要说明的是，术语"安装"、"相连"、 "连接"应做广义理解，例如，可W是机械连接或电连接，也可W是两个元件内部的连通，可 W是直接相连，也可W通过中间媒介间接相连，对于本领域的普通技术人员而言，可W根据 具体情况理解上述术语的具体含义。
[0092] 下面结合附图1对本发明做进一步描述。
[0093] 参照附图1，一种适用于快时变稀疏估计的最优导频图样的设计方法，具体实施步 骤如下：
[0094] 步骤100,开始；
[0095] 步骤110，初始化D = 0，化=1，存储D值的变量a = 0，存储化值的变量b = 1，阔值ξ = 10;
[0096] 步骤120,建立信道模型方程式，根据式(24)计算μ(Φ);
[0097] 步骤130,令Pg值保持不变，D = D+1，逐步增大D值；
[009引步骤140,寻找能使其μ(φ)最小的D值，判断是否满足μ(Φ)ο+ι<μ(Φ)ο?若满足， 跳转步骤150,若不满足，跳转步骤160;
[0099] 步骤 150，a = D+l，进入步骤 170;
[0100] 步骤 160，a = D;
[0101] 步骤170，判断是否满足"<^?若满足，跳转步骤130，若不满足，将最终的3值赋 给D值，即D = a，保存D值后跳转步骤180;
[0102] 步骤180,令D值保持不变，& = Pg+1，逐步增大化值；
[010：3]步骤190,寻找能使μ(Φ)最小的Pe值，判断是否满足片(<?),。,1<从。\?若满足， 跳转步骤200，若不满足，则跳转步骤210;
[0104] 步骤200，b = PG+l，进入步骤 220;
[0105] 步骤 210，b = Pc;
[0106] 步骤220,判断是否满足Λ 若满足，跳转步骤180,若不满足，贝峭巧专步骤230;
[0107] 步骤230，将最终的b值赋给Pg值，即& = b，保存化值；
[010引步骤240,令P/Pg=[P/Pg]賴，其中[·]賴表示按照四舍五入法进行取整；
[0109] 步骤250，判断是否满足^ < I()?若满足，进入步骤260，若不满足，跳转步骤330;
[0110] 步骤260,随机的选择P/Pe个导频组信道得到原始的导频信道图案；
[0111] 步骤270,选定其中一个导频组(不再选择已最优处理过的导频组），其他导频组保 存不变；
[0112] 步骤280，改变选定导频组的位置；
[0113] 步骤290,重新计算μ(巫）；
[0114] 步骤300,判断此时的μ(Φ)是否最小?若是，跳转步骤310;若否，跳转步骤280;
[0115] 步骤310,保存此时运个导频组的位置；
[0116] 步骤320,判断是否所有的导频组都进行了最优处理，即寻找最优导频组的次数i 是否进行了 P/化次，写成数学表达式即：是否满足若满足，跳转步骤440,若不满足， 尸保 则跳转步骤270;
[0117] 步骤330,将N个子信道近似均匀的分为[P/Pg/K]賴个块空间，其中P/化是导频信道 个数，而K为算法复杂度的一个参数，一般设为1-3;
[0118] 步骤340，在每个块区间中，随机的选择K个子信道作为导频组信道，从而得到原始 的导频信道图案；
[0119] 步骤350,随机选择一个块区间（不再选择已最优处理过的块区间）；
[0120] 步骤360,在区间内，选定其中一个导频组位置(不再选择已最优处理过的导频组 位置），其他导频组保存不变
[0121 ]步骤370，改变选定导频组的位置；
[0122] 步骤380,计算此时的μ(巫）；
[0123] 步骤390,判断此时的μ(Φ)是否最小?若是，进入步骤400;若否，跳转步骤370;
[0124] 步骤400,保存此导频组位置；
[01巧]步骤410,判断是否所有的导频组都进行了最优处理，即若是，进入步骤 420,若否，跳转步骤370;
[01 %] 步骤420，保存此块区间Κ个导频组的图样；
[0127]步骤430,判断是否所有块区间都进行了最优处理，即是否寻找块区间最优导频图 样变量j是否执行了 [P/Pg/K]5?次，写成数学表达式即是否满足i> ?若满足，贝峭b -心取营 转步骤440，若不满足，则跳转步骤350;
[0128] 步骤440,保存并整合导频图样，即寻找到最优的导频图样；
[0129] 步骤450,结束。
[0130] 在本说明书的描述中，参考术语"一个实施例"、"一些实施例"、"示例"、"具体示 例"、或"一些示例"等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特 点包含于本发明的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不 一定指的是相同的实施例或示例。而且，描述的具体特征、结构、材料或者特点可W在一个 或多个实施例或示例中W合适的方式结合。
[0131] 尽管已经示出和描述了本发明的实施例，本领域的普通技术人员可W理解:在不 脱离本发明的原理和宗旨的情况下可W对运些实施例进行多种变化、修改、替换和变型，本 发明的范围由权利要求及其等同物限定。
【主权项】
1. 本发明提出一种适用于快时变稀疏估计的最优导频图样的设计方法，具体为： 针对快时变信道模型，并结合稀疏信道估计，提出一种最优化导频图样的设计方法； 该方法包含的创新有： S1，目前稀疏信道估计方法大多适用于慢时变或时不变信道，对信道稀疏性的利用只 局限于时延域而没有扩展到多普勒域，并且忽略了由多普勒频移产生的载波间干扰（Inter Carrier Interference，ICI)对信道造成的影响，因此针对在快时变信道中采用稀疏信道 估计方法的精确性和可靠性会大幅下降的问题，提出一种适用于快时变信道的稀疏信道估 计方法； S2,针对目前大部分导频图样设计均基于密集多径的传统信道估计方法进行，需要大 量的导频信号来获取准确的信道状态信息的问题，提出一种适用于快时变稀疏信道估计的 最优导频图样设计方法。2. 根据权利要求1所述的一种适用于快时变稀疏估计的最优导频图样的设计方法，其 特征在于，所述权利要求1中S1提出的一种适用于快时变信道的稀疏信道估计方法，包括： 为方便阐述，以单输入单输出-正交频分复用系统（3;1叫16 1叩1^3;[叫16〇1^口111：-Orthogonal Frequency Division Multiplexing，SIS〇-〇FDM)系统为例阐述问题，其中关 于多输入多输出_〇FDM(Multiple Input Multiple 0utput-0FDM，]\OM0-0FDM)的系统同理； 在快时变信道中，假设有N个子载波，通信系统的传输模型为：其中，x(t)是发射端信号，y(t)是接收端信号，z(t)是高斯白噪声，h(t，T)是信道的时 域冲击响应，可表不为：其中，S为不为0的抽头数，<^，^，^分别是第8(8=1，...，5)条传播路径的衰落系数、时 延和多普勒频移;对于快时变信道来说，ICI的影响不能被忽略，所以其对应的离散模型为：其中，T为采样时间，L为信道最大的多径时延，1=0,. . .，L-l，n = 0,. . .，N-1，ts/T<L，S <L；可写为矩阵形式，如下： y = hx+z 其中，y=(y[〇]，· · ·，y[N-l])T，x=(x[0]，· · ·，χ[Ν-1])τ，ζ = (ζ[0]，· · ·，z[N-l])T，h为 信道矩阵，有如下形式：令X表示发射端的数据信号，经快速傅里叶逆变换（Inverse Fast Fourier Transform，IFFT)后转换为发射信号x，在接收端，接收信号y经快速傅里叶变换（Fast Fourier Transform，FFT)后得到接收端的数据信号Y ;矩阵Q表示离散傅里叶变换 (Discrete Fourier Transform，DFT)矩阵P和q分别代表 DFT矩阵的行和列，1彡p，q彡N，则有:x = QHX，Y = Qy，高斯白噪声Z = Qz，系统矩阵H=QhQH;假 设系统的循环前缀长度大于最大的多径时延L，那么符号间干扰（Inter Symbo 1 Interference，I SI)可以近似忽略，则可得到完整的系统频域传输模型如下：Y = HX+Z 为了充分的利用快时变信道模型的时频稀疏特性， 转换到延迟-多普勒域，得到信道的基扩展模型函数，即：而1^[11，1]表示单个抽头的基扩展模型;分别设： Φ s(l) = sinc(3r(l-Ts/T))其中Φ8(1)和也⑷分别描述了信道的时间延迟和多普勒频移，而Φ8(1)和Wd)的能量 分别集中于以散射点is/I^PfsNT为中点的邻域，即每个子载波的能量泄露只影响其相邻的 几个子载波，所以可以采用带状近似矩阵;假设每个子载波的频率偏移最大的距离为D，则 矩阵Η被近似为：该带状近似矩阵只保留了 Η的对角线以及对角线两边距离D之间的元素;假设φ8((1)的能 量主要集中于[_D，D]区间，则信道基扩展模型函数us[d，l]以及u[d，l]的支撑集由[0，L-1] \[0少-1]缩小为[0丄-1]\[0,0]1^10，^1]，其中[^0，^1]由[-0，-1]循环移位所得； 因此，需要估计的信道参数数量由NL减少为(2D+1) L; 为了更方便地估计信道参数，将Y = HX+Z转换为以下形式： Y=Ah+Z 其中，信道冲击响应向量h定义为： h = (ho, . . . , hL-i), hi=(h[0,l], . . . ,h[N-l,l]) 贝1JN X NL维矩阵A中的元素定义为：其中，k和q分别代表行和列，n = (q-1 )mod(N)，1 = (q-l_n)/N;将待估计目标h转换为信 道的基扩展模型，而又由于u[d，l]的支撑集为[0丄-1]\[0,0川[^0，^1]，因此定义(20+ 1)L维向量u有： U= (u〇, . . . ,UL-l)T, ui=(u[l,0],...,u[l,D],u[l,N-D],...,u[l,N-1]) 因为u具有稀疏特性，所以结合:可得到如下关系： 其中，II为L阶单位矩阵，F为(2D+1) X N维的单位选择矩阵，表示从[1，N]抽取子集[1，D] U [N-D+l，N]的采样操作，Q为DFT矩阵；用u代替h进行求解，则将Y = Ah+Z转换到延迟-多普 勒域，得到： Υ=Φιι+Ζ 其中，Φ称为测量矩阵，且Φ = A(I, ?(FQ)H); 而γ= Φι1+Ζ可以看作是求解稀疏信号的线性方程；Φ = A(I￡ 0(FQ)H)可推导出矩阵Φ中的元素具有如下表达式：其中，n = (mq-l)mod(N)，1 = (mq-1-n)/N，这里（2D+1)L维向量m定义为： rn= (mo, . . . , niL-i), mk=(Nk+1,...,Nk+D+l,Nk+N-D+l,...,Nk+N)可知，在矩阵Φ的第k行，Φ k,q的组成成分包含发 η-1 )mod(N)+l ; ξΙξΙΠ - 0日寸，X(k-η-1 )mod(N)+l - Xk ； ΠΤ] ξΙξΙΠ 辛 0日寸，X(k-η-1 )mod (N) + l G { Xk-D，· · ·， Xk-hXk+i，. .，Xk+D};因此，第k个观测值Yk的取值不仅取决于Xk，还受到Xk的长度为（2D+1)的 邻域内其他符号的影响，而这些相邻符号对Yk的作用即等价于ICI干扰； 假设在子载波N中有P个导频符号，即XN(1)，. . .，Xn(p)，则有：在Φ的第P行中，由于观测值YN(P)受到与导频符号XN(P)相邻的数据符号的干扰，矩阵仍 然存在部分未知元素，无法直接得到；于是，假设将Φ中包含的未知元素全部置零，得到Φ 的近似矩阵Φ ;设导频子载波集合为Ω ={k(l)，k(2)，. . .，k(P)}，则φ可表示为：其中，n=(mq-l)mod(N)，且k(p) e Ω ;由于只存在由导频符号构成的点，于是.φΓ完全已 知;假设w = (Φ - Φ )n + Ζ，用w替代噪声之，得到近似的测量方程： Y = Φ'ιι..+ w 因为备满足有限等距性质（Restricted Isometry Property，RIP)，之后即可采用压缩 感知((3〇11^68 8；^36118；[1^，05)重构方法求解待估计量11，转而通过11_=.(1^@([(3严)11求得 信道参数h; 但是，采用这种近似处理方法会带来两种误差： 一是将Φ中包含的未知元素全部置零的近似处理引入的误差由）u,等价于ICI干 扰;二是CS的重构误差； 针对上述第一种误差，应降低ICI干扰对估计方法的影响，即降低导频符号所受到的邻 近的数据符号的干扰，因此最优的导频排列方案是采用连续导频图案，即将所有导频符号 集中放置在一起;而针对上述第二种误差，应提高稀疏信号重建的精度和稳定性;其中CS算 法重建性能主要由测量矩阵的RIP决定，所以在现有的稀疏信道估计方法中，最优的导频排 列方案是采用完全随机的导频图案，即将所有导频符号随机地插入到数据符号中，以保证 CS重建算法具有相对稳定的性能。 因此，为了能同时有效地降低ICI干扰和重建误差，提出采用随机分组导频图案：即将 导频符号分为若干个组，再将这些导频组以随机的位置地插入到数据符号中，Pc为导频的 分组长度； 由于直接证明RIP性质以及计算RIP常数难度大，因此本发明采用矩阵的相干系数替代 RIP常数来评价导频图案对于CS重建算法的重建性能的影响； 矩阵Φ的相干系数μ( Φ )定义为：其中，Φ i和Φ j都是Φ的列向量。3.根据权利要求1所述的一种适用于快时变稀疏估计的最优导频图样的设计方法，其 特征在于，所述权利要求1中S2提出的一种适用于快时变稀疏信道估计的最优导频图样设 计方法，包括： 基于权利要求2所述的信道估计方法，提出一种适用于快时变稀疏估计的最优导频图 样的设计方法，具体流程如下：步骤1、初始化D = 0，Pg=1; 步骤2、建立信道模型方程式 步骤3、令PC值保持不变，D = D+1:区间内逐步增大D值，寻找能使其μ( Φ )最小的D值并保存D值；步骤4、保持寻找到最优的D值不变，令PC = PC+1， 区间内逐步增大Pc值，寻找能 使其μ( Φ )最小的Pg值，并保存Pg值； 步骤5、判断P/Pe是否为整数，若是，进入步骤6;若否，则令P/Pe=[P/Pe]?e，其中[· ]?β 表示按照四舍五入法进行取整，此时则存在其中某一个导频组比其他导频组多或者少一个 导频数的情况，进入步骤6; 步骤6、判断P/Pc是否小于阈值ξ，若是，进入步骤7，若否，跳转步骤10; 步骤7、随机的选择P/PC个导频组信道得到原始的导频信道图案； 步骤8、随机改变其中一个导频组位置(不再选择已最优处理过的导频组），重新计算μ (Φ )，直至寻找到拥有最小的μ( Φ )的导频组位置，保存这个导频组位置； 步骤9、重复步骤8第P/PG次，直到寻找到所有P/PG个导频组的位置，保存，即寻找到最优 的导频图样； 步骤10、将N个子信道均勾的分为[P/PC/K]職个块空间，其中P/PC是导频组信道个数，而 K为算法复杂度的一个参数，一般设为1-3; 步骤11、在每个块区间中，随机的选择K个子信道作为导频组信道，从而得到原始的导 频信道图案； 步骤12、随机选择一个块区间； 步骤13、在区间内，随机改变其中一个导频组位置(不再选择已最优处理过的导频组）， 重新计算μ( Φ )，直至寻找到拥有最小的μ( Φ )的导频组位置，保存这个导频组位置； 步骤14、重复步骤13第Κ次，直到寻找到这个区间内所有Κ个导频组的最优位置，保存； 步骤15、随机选择一个新的块区间，重新执行步骤13、步骤14,寻找并保存这个区间内Κ 个导频组的最优位置； 步骤16、重复步骤15第[PAVK]職次，直至寻找到所有导频组的位置，保存并整合，即寻 找到最优的导频图样； 关于判断导频图样是否采用分块处理，通过阈值ξ来决定，而一般设定阈值ξ多10。
【文档编号】H04L5/00GK106059731SQ201610338017
【公开日】2016年10月26日
【申请日】2016年5月19日
【发明人】廖勇, 周昕
【申请人】重庆大学