一种基于压缩感知技术的ofdm系统参数化信道估计及均衡方法
【专利摘要】本发明公开了一种基于压缩感知技术的OFDM系统参数化信道估计及均衡方法，信道估计方法包括以下步骤：1)由导频子载波处的OFDM频域接收数据，基于信道的参数化模型，对信道进行参数估计：首先，根据MDL准则，得到信道的多径数目估计值然后，利用ESPRIT技术估计每条路径的延时参数，得到多径延时估计值2)利用参数估计结果构建冗余字典矩阵，通过BPDN算法求解信道响应估计值。本发明基于参数化信道建模，将信号的参数估计技术与压缩感知框架相结合，提高了信道估计的精度和通信效率，降低了通信误码率。与传统信道估计方法相比，本发明具有明显的性能优势。
【专利说明】
一种基于压缩感知技术的OFDM系统参数化信道估计及均衡 方法
技术领域
[0001]本发明涉及通信技术领域，主要涉及一种基于压缩感知技术的OFDM系统参数化信 道估计及均衡方法。
【背景技术】
[0002] 无线通信系统性能受无线信道的影响，其通信环境不同于有线通信，具有较大的 不确定性和随机性。通常来说，阴影衰落、多径效应和多普勒效应等都会对无线信道产生影 响，使无线信道呈现时间选择性衰落和频率选择性衰落。无线通信的接收端需要对无线信 道进行补偿，这就需要进行信道估计，以获得无线信道中信道阶数、冲激响应、多普勒频移 等参数。信道估计的精度会直接影响整个系统的性能，也是进行相关检测、解调和均衡等技 术的基础。因此，信道估计是无线通信系统中的一项关键技术。
[0003] 正交频分复用（Orthogonal Frequency Division Multiplexing,0FDM)多载波的 出现，为信道估计技术的应用提供了更广阔的空间。CFDM系统的多个子载波在频率上相互 重叠而且正交，形成了多个子信道，调制信号在每个子信道上进行传输，从而能够传递高速 的数据流。(FDM系统具有抗多径干扰能力强，频带利用率高等优点，因此受到广泛关注。目 前，作为物理层的核心技术，OFDM已经被多个通信标准采纳。例如:欧洲数字音频广播(DAB) 标准采用了使用差分相位调制方式的OFDM技术，数字视频广播(DVB)标准使用了多幅度调 制方式的OFDM技术，无线局域网IEEE 802. Ila标准、无线城域网IEEE 802.16标准等都采用 了OFDM作为核心技术。
[0004]压缩感知是近几年在应用数学领域和信号处理领域受到广泛关注的一种数学方 法。压缩感知理论利用变换空间描述信号，通过直接采集得到少数"精挑细选"的线性观测 数据，将信号的采样转变成信息的采样，并通过信号重建算法恢复原始信号。现有的利用压 缩感知技术进行OFDM系统信道估计的方法有一个共同的缺陷，就是在使用过完备的字典来 提升重建精度的同时，由于字典项数的变大，计算量大大增加。
[0005] 另一类利用信道稀疏特性的信道估计方式是基于参数化模型的信道估计方法。这 类方法通过建立参数化(多径数和多径延时）的信道模型，降低信道估计问题的维度，从而 得到更好的信道估计效果。然而，与基于压缩感知的信道估计方法一样，利用参数化模型的 信道估计方法都有一个共同的问题就是计算复杂度太高，同时，与采用过完备、超分辨率字 典的基于压缩感知的信道估计方法相比，估计效果也没有明显的优势。
[0006] 因此，有必要设计一种兼顾算法精度与复杂度的OFDM系统信道估计方法。
【发明内容】
[0007] 本发明所解决的技术问题是，针对现有技术的不足，提出了一种基于压缩感知技 术的OFDM系统参数化信道估计及均衡方法，利用参数估计方法改进压缩感知技术进行信道 估计及均衡，在不增大字典矩阵的大小的情况下，提高OFDM系统信道估计精度，改善通信系 统的性能。
[0008] 为了解决上述技术问题，本发明所提出的技术方案是：
[0009] 一种基于压缩感知技术的OFDM系统参数化信道估计方法，包括以下步骤：
[0010] 1)由导频子载波处的OFDM频域接收数据，基于信道的参数化模型，对信道进行参 数估计:首先，根据MDL【Minimum Description Length，最小描述长度】准则，得到信道的多 径数目估计值?;然后，利用ESPRIT技术估计每条路径的延时参数，得到多径延时估计值
[0011] 2)利用参数估计结果构建冗余字典矩阵，使用BPDN算法作为压缩感知信号重建算 法，通过BPDN算法求解信道响应估计值。
[0012] 利用信道冲激响应的稀疏特性建立信道的参数化模型如下式：
[0013]
[0014] 其中，Μτ，τ)为时变冲激响应，t为时间，τ为延迟时间参数，hKt)表示第1条路径 在时亥叶的复增益，δ(τ-τι)表示位移^后的冲激函数(狄拉克函数）， τι表示第1个条路径的 时延，L表示多径数目；Ts是OFDM系统的基带采样间隔，有Ts =|，B为OFDM系统的基带带宽。 对于瑞利信道而言，lu(t)是一个瑞利随机过程，即广义窄带复高斯随机过程，其均值为0， 方差为1，并具有Jakas形式的功率谱密度。不同的路径增益之间不相关;为以基带采 样间隔为单位的延迟参数，为了简单起见，假设ξι为整数。
[0015] 所述计算信道的多径数目估计值￡的方法为：
[0016] 首先计算导频子载波处的最小二乘估计值：
[0017]
[0018] 其中，表示第i个OFDM符号中第m个导频子载波处的信道响应最小二乘估计 值;P(m)表示第m个导频所在的子载波编号，为第m个导频的发射数据，Y1,Ρω表示第i个 (FDM符号中第m个导频子载波的接收信号，N为OFDM系统子载波个数，m,P(m)表示第i个OFDM 符号中第m个导频子载波的加性高斯噪声;令,爲>:ρ(Μ_υ；|7表示第i 个OFDM符号中所有M个导频子载波处的信道响应最小二乘估计值，M〈N【子载波中有一部分 作为导频子载波，其余的子载波为数据子载波】，为了使用MDL准则估计L的值，将重 构成如下的快照矩阵(Snap Shot Array):
[0019]
[0020] 其中，K、M为可调参数，K的取值影响相关矩阵的精度【在仿真实验中，取M=导频子 载波数目，K = 5,可以获得较好的估计结果】，该参数需满足条f
[0021] 然后，计算第i个OFDM符号的相关矩阵：
[0022]
[0023]其中，Q(i)H表示矩阵Q(i)的共辄转置，J为(M-K+l) X (M-K+1)的单位矩阵，3;表示 矩阵A的共辄矩阵；
[0024]然后，对连续I个OFDM符号的相关矩阵求平均：
[0025]
[0026] 然后，对1进行特征值分解：
[0027]
[0028] 其中，之.1.2 S ··:_·: 2 是降序排列的特征值，_Μι，是对应的特征 向量；
[0029] 根据MDL准则，多径数目估计值￡由下式确定：
[0030]
[0031] 其中，arg min表示求使得MDL(P)达到最小值时的P;
[0032]
[0033] 其中，Π 为连乘符号。
[0034]采用ESPRIT技术计算多径延时估计值的方法为：
[0035] 首先将务最大的I：个特征值对应的特征向量排列成一个矩阵：
[0036]
[0037]
[0038]
[0039]
[0040] 其中，Ιμ-κ为(M-K) X (M-K)的单位矩阵；
[0041 ] ESPRIT计算矩阵为：
[0042]
[0043] 对Φ进行特征值分解，得到$的￡个特征值为·则第1条路径的路径延时估 计值为：
[0044]
[0045] 其中，argOf)表示复数特征值的共辄柯啲相角，Df表示导频间隔（即相邻导频之 间子载波的个数）。
[0046] 所述步骤2)具体为：
[0047] 2.1)利用步骤1)得到的多径数目估计值I：和多径延时估计值f，构建冗余字典矩阵 Ψ：
[0048] 首先计算出每条信道路径对应的字典原子项数目，第1条路径对应的字典原子项 数目2m_l【该值没有特定的取值限制，在本方法中可以简单取Ifj】，保证所有信道路径对应 的字典原子项数目之和为给定的字典规模【给定的字典规模可以是任意值，本方法实现中 选取为OFDM子载波个数，即N】；
[0049] 接着对每条信道路径，分别构造其过完备子字典Ψ i:
[0050]
[0051] q= [0,1，···，N-1]T^1条路径对应的过完备字典的规模为NX (2m-l)，所有的m满足分配 关系
[0052] 最后，将所有Ψι合并为新的过完备字I
[0053] 2.2)使用BPDN算法求得信道频率响应估计值M，其中第i个OFDM符号的信道频率响 应估计值记为見,计算方法如下：
[0054]
[0055]
[0056]
[0057] 其中，Cn表示N维复数空间；丫户出,[)(())丄,[)(1)，"土, [)0?-1)]7，民为第1个0?01符号 的信道冲激响应估计值，μ是调节算法精度与稀疏度的参数，μ越大得到的馬越稀疏，σ为信 噪比；Φ为一个MXN的矩阵，满足关系
该矩阵符合准托 普林兹矩阵结构，即每一行都是由上一行移位得到。相关研究证明，准托普林兹矩阵是符合 RIP原则的，因此可以作为观测矩阵。
[0058] 一种基于压缩感知技术的OFDM系统参数化信道均衡方法，采用上述的基 于压缩感知技术的OFDM系统参数化信道估计方法得到信道频率响应估计值見利用
对接收信号进行信道均衡，以抵消信道乘性噪声对无线 信号的影响：
[0059]
[0060]其中，Yi(k)为第i个OFDM符号中第k个子载波的接收信号；
[0061]
iI即为信道均衡后的接收信号。
[0062] 有益效果：
[0063]本发明对传统的OFDM系统信道估计技术进行优化，在基于压缩感知技术的无线信 道估计技术基础上，利用参数估计技术，对压缩感知中的冗余字典矩阵(基矩阵)进行动态 构建和优化，在不增大字典矩阵的大小的情况下，使信道估计的结果误差更小，提高系统的 通信效率，改善通信系统的性能。
[0064] 本发明基于参数化信道建模，将信号的参数估计技术与压缩感知框架相结合，提 高了信道估计的精度和通信效率，降低了通信误码率。与传统信道估计方法相比，本发明具 有明显的性能优势。
【附图说明】
[0065] 图1是使用本发明进行无线通信的方案流程图；
[0066] 图2是使用本发明进行无线通信的系统框图；
[0067]图3是本发明的算法实施示意图1;
[0068]图4是本发明的算法实施不意图2;
[0069]图5是压缩感知框架信号处理框图；
[0070] 图6是本发明的信道估计误差仿真结果；
[0071] 图7是本发明用于OFDM系统的通信误码率仿真结果。
【具体实施方式】
[0072]下面结合附图和特定的实例，对本发明的【具体实施方式】作进一步详细描述。
[0073]使用本发明提供的一种基于压缩感知技术的OFDM系统参数化信道估计方法进行 无线通信的流程如下：
[0074] S1、在数据发送端将待发送的数据插入梳状导频后进行OFDM调制并发送；
[0075] S2、在数据接收端，对OFDM频域接收数据，利用ESPRIT(estimation of signal parameters by rotational invariance techniques，基于旋转不变技术的信号参数估 计)技术和基于信道的参数化模型，获得信道的参数估计;包括无线信道的多径数目估计和 多径延时估计；
[0076] S3、利用S2获得的参数估计结果，包括信道多径数目和多径延时的估计参数，构建 冗余字典矩阵，使用BPDN算法作为压缩感知信号重建算法，通过BPDN算法求解信道响应估 计值【重建信道冲激响应和频率响应】；
[0077] S4、使用S3中求解得到的信道估计结果对接收数据进行信道均衡，并进行OFDM解 调输出。
[0078] 具体流程如图1所示，实施过程如下：
[0079] 第一步，如图2前9部分所示，在数据发送端，待发送的初始数据经过16QAM(即包含 16种符号的QAM调制方式)调制后得到I/Q通道的待发送数据，串/并转换步骤将待发送数据 分配到不同的数据子载波，并将用于信道估计的导频数据分配到导频子载波上，得到OFDM 频域发送数据;OFDM频域发送数据经过IFFT运算后得到时域发送数据，插入循环前缀，进行 并/串转换和数模转换后得到基带时域发送波形;将基带时域发送波形进行上变频处理调 制到射频频段，并经天线辐射到自由空间。
[0080] 如图2第10-15部分，无线信号经信道传递后到达数据接收端(接收机），下变频后 得到基带时域接收信号;将基带时域接收信号经模数转换处理进行采样量化，并根据OFDM 系统同步结果将采样量化后的基带时域接收信号分隔为单个OFDM符号帧;将OFDM符号帧进 行串/并转换，并去除循环前缀后，经过FFT得到OFDM频域接收数据。
[0081] 第二步，对OFDM频域接收数据按图2第16部分所示进行处理:算法实施流程如图3 和图4所示，利用ESPRIT技术和基于信道的参数化模型，对信道进行参数估计，得到信道的 多径数目估计L和多径延时估计f = , 亥步骤包括两部分:首先，由导频子载波 处的接收数据，根据MDUMinimum Description Length，最小描述长度)准则，估计出信道 的多径数目然后，利用ESPRIT技术估计每条路径的延时参数，得到多径延时估计 τ = Iri, TyjTiJl:
[0082] 信道的参数化模型如下式：
[0083]
[0084]其中，lu(t)表示第1条路径的复增益，T1表示第1个条路径的时延，L表示多径数目； Ts是OFDM系统的基带采样间隔，有7〉=|,B为OFDM系统的基带带宽。对于瑞利信道而言，hi (t)是一个瑞利随机过程，即广义窄带复高斯随机过程，其均值为0,方差为1，并具有Jakas 形式的功率谱密度。不同的路径增益之间不相关。& = g为以基带采样间隔为单位的延迟 参数，为了简单起见，假设ξι为整数。
[0085] 首先计算导频子载波处的最小二乘估计值：
[0086]
[0087] 其中，岛滅_表示第i个OFDM符号中第m个导频子载波处的信道响应估计值;p(m) 表示第m个导频所在的子载波编号，为第m个导频的发射数据，Y1,pW表示第i个OFDM符号 中第m个导频子载波的接收信号，N为OFDM系统子载波个数，m, P (m)表示第i个OFDM符号中第m 个导频子载波的加性高斯噪声。
[0088] -
良示第i个OFDM符号中所有M个导频子载 波处的信道响应最小二乘估计值，M〈N【子载波中有一部分作为导频子载波，其余的子载波 为数据子载波，如图4所示】;为了使用MDL准则估计L的值，首先将重构成如下的快照 矩阵（Snap Shot Array):
[0089]
[0090] 其中，K、M为可调参数，K的取值影响相关矩阵的精度。该参数需满足条件M_K+1> Cl ο
[0091] 然后，计算相关矩阵：
[0092]
[0093] 其中，Q(i)Η表示矩阵Q⑴的共辄转置，J为(M-K+l) X (Μ-Κ+1)的单位矩阵，互:表示 矩阵A的共辄矩阵。
[0094]然后，对连续I个OFDM符号的相关矩阵求平均：
[0095]
[0096] 然后，对嚴进行特征值分解：
[0097]
[0098] 其中，I1之I2之…2 是降序排列的特征值，是对应的特征 向量。MDL准0丨丨可以表示为：
[0099]
[0100] 中，Π 为连乘符号；
[0101]根据MDL准则，多径数目的估计值由下式确定：
[0102]
[0103] 采用ESPRIT计算多径延时估计值的方法如下，首先将免最大的I个特征值对应的 特征向量排列成一个矩阵：
[0104]
[0105]
[0106]
[0107]
[0108] 其中，Im-k为(M-K) X (M-K)的单位矩阵。
[0109] ESPRIT计算矩阵为：
[0110]
[0111] 对Φ进行特征值分解，$的￡个特征值为{^好=1，则第1条路径的路径延时估计值 为：
[0112]
[0113] 其中，argOf)表示复数特征值的共辄扣的相角，Df表示导频间隔（即相邻导频之间 子载波的个数)。
[0114] 第三步，如图2第17部分所示，利用S2步骤得到的多径数目估计￡和多径延时估计 元构建冗余字典Ψ:
[0115] 先计算出每个信道路径对应的字典原子项数目，保证所有信道多径对应的字典原 子项数目之和为给定的字典规模。
[0116] 接着对每个信道路径，分别构造其过完备子字典：
[0117]
[0118] q= [0，1，…，N_1]T。其中，第1条路径对应的过完备字典的规模为NX (2m-l)，所有的m满足 分配关J
[0119] 最后，将所有W1合并为新的过完备字典少=[^,?,…,;在此基础上，
[0120] 如图5所示，使用BPDN算法求得信道频率响应估计值身:，其中第i个OFDM符号的信道 频率响应估计值记为总,，计算方法如下：
[0121]
[0122]
[0123]
[0124] 其中，为第i个OFDM符号的信道冲激响应估计值，CN表示N维复数空间；Yi = [Yi,P(Q)，Yi,P(1)，…，Yi,P(M-υ] Τ，μ是调节算法精度与稀疏度的参数，μ越大得到的越稀疏，。 为信噪比；Φ为一个MXN的矩阵，满足关系
亥矩阵符合 准托普林兹矩阵结构，即每一行都是由上一行移位得到。相关研究证明，准托普林兹矩阵是 符合RIP原则的，因此可以作为观测矩阵。
[0125] 第四步，如图2第1 8，1 9，2 0，2 1部分所示，利用信道响应估计结果 瓦=[反(0),/?#)風(iV -Ι)]对接收信号进行信道均衡，以抵消信道乘性噪声对无线 信号的影响：
[0126]
[0127] 其中，Yi(k)为第i个OFDM符号中第k个子载波的接收信号。
[0128]
卩为信道均衡后的接收信号。均衡后的接收信号 经过并/串转换后得到接收调制数据序列，该序列进行16QAM解调步骤后，即得到OFDM通信 系统输出数据。
[0129] 本发明的性能仿真结果如图6、图7所示。由图6可以看出，与传统的基于压缩感知 技术的信道估计方法相比，本发明可以有效提高OFDM系统的信道估计精度；由图6可以看 出，与基于传统冗余字典的压缩感知信道估计方法相比，本发明可以实现较低的通信误码 率。因此，本发明与传统的基于压缩感知的信道估计方案相比，具有明显的性能提升。
[0130]由以上实施例可以看出，本发明能用来有效地提高OFDM系统的信道估计效果。与 传统的方法相比，在不降低系统数据传输效率的情况下，本发明具有更高的信道估计精度， 从而降低系统的通信误码率，提高通信可靠性。
【主权项】
1. 一种基于压缩感知技术的OFDM系统参数化信道估计方法，其特征在于，包括W下步 骤： 1) 由导频子载波处的OFDM频域接收数据，基于信道的参数化模型，对信道进行参数估 计:首先，根据最小描述长度MDL准则，得到信道的多径数目估计值￡;然后，利用ESPRIT技术 估计每条路径的延时参数，得到多径延时估计2) 利用参数估计结果构建冗余字典矩阵，使用BPDN算法作为压缩感知信号重建算法， 通过BPDN算法求解信道响应估计值。2. 根据权利要求1所述的基于压缩感知技术的OFDM系统参数化信道估计方法，其特征 在于，利用信道冲激响应的稀疏特性建立信道的参数化模型如下式：其中，h(T，t)为时变冲激响应，t为时间，T为延迟时间参数，hi(t)表示第1条路径在时刻 t的复增益，S(T-Tl)表示位移Tl后的冲激函数(狄拉克函数），Tl表示第1个条路径的时延，L 表示多径数目；Ts是OFDM系统的基带采样间隔，有为OFDM系统的基带带宽I W基带采样间隔为单位的延迟参数，为了简单起见，假设Cl为整数。3. 根据权利要求2所述的基于压缩感知技术的OFDM系统参数化信道估计方法，其特征 在于，所述计算信道的多径数目估计值I；的方法为： 首先计算导频子载波处的最小二乘估计值：其中，巧，表示第i个OFDM符号中第m个导频子载波处的信道响应最小二乘估计值;P (m)表示第m个导频所在的子载波编号，丫 m为第m个导频的发射数据，Yi,PW表示第i个OFDM 符号中第m个导频子载波的接收信号； 令成化；二[巧,>([))，巧,灿)，…，巧诚M-i)f表示第i个0抑M符号中所有M个导频子载波处的 信道响应最小二乘估计值，M<N，N为(FDM系统子载波个数;为了使用MDL准则估计L的值，将 所WV重构成如下的快照矩阵：其中，K、M为可调参数，K的取值影响相关矩阵的精度，该参数需满足条件M-K+1 > U， 巧 然后，计算第i个OFDM符号的相关矩阵：其中，Q(i)H表示矩阵W1;的巧狐巧置，」刃W-K+U X 的单位矩阵，a表示矩阵A 的共辆矩阵； 然后，对连续I个O抑M符号的相关矩阵求平均： 然后，对度进行特征值分解：其中，义1兰毛含"'至為w:-K+'i是降序排列的特征值，化1, ...,Wm-K+1.是对应的特征向量； 根据MDL准则，多径数目估计值I；由下式确定：其中，arg min表示求使得MDL(P)达到最小值时的P;其中，n为连乘符号。4. 根据权利要求3所述的基于压缩感知技术的OFDM系统参数化信道估计方法，其特征 在于，采用ESPR 口技术计算多径延时估计值的方法为：首先将巧最大的I：个特征值对应的特佈向量排列成一个巧阵： 然后，令： 其中，Im-K为(M-K) X (M-K)的J ESPR 口计算矩阵为： 对4进行特征值分解，得到参的￡个特祉但刃，则第1条路径的路径延时估计值 为：其中，a巧Op表示复数符化但的巧狐巧的和用，Df巧不导娜间隔。5. 根据权利要求4所述的基于压缩感知技术的OFDM系统参数化信道估计方法，其特征 在于，所述步骤2)具体为： 2.1)利用步骤1)得到的多径数目估计值￡和多径延时估计值t，构建冗余字典矩阵W : 首先计算出每条信道路径对应的字典原子项数目，第I条路径对应的字典原子项数目 化1-1，保证所有信道路径对应的字典原子项数目之和为给定的字典规模； 接着对每条信道路径，分别构造其过完备子字典W1: 其4- Q= [0,1，…，N-1]|;第1条路径对应的过完备字典的规模为NX (2ni-l)，所有的m满足分配 关署最后，将所有W1合并为新的过完备字I2.2)使用BPDN算法求得信道频率响应估计值化，其中第i个OFDM符号的信道频率响应估 计值记为蘇,心十算方法如下： 0 =巫W其中，砖为第i个OFDM符号的信道沖徽响应估计值，C"表示N维复数空间；Yi= [Yi,P(O), Yi,P(I),…，Yi,p(M-i)]T，y是调节算法精度与稀疏度的参数，y越大得到的馬越稀疏，O为信噪 比；O为一个M X N的矩阵，满足W下关系：6. -种基于压缩感知技术的OFDM系统参数化信道均衡方法，其特征在于，采用权利要 求1~5中任一项所述的基于压缩感知技术的OFDM系统参数化信道估计方法得到信道频率 响应估计值寐，利用I接收信号进行信道均衡，W抵消信 道乘性噪声对无线信号的影响：其中，Yi化)为第i个OFDM符号中第k个子载波的接收信号；为信道均衡后的接收信号。
【文档编号】H04L25/02GK105915473SQ201610357016
【公开日】2016年8月31日
【申请日】2016年5月26日
【发明人】邓宏贵, 杜捷
【申请人】中南大学