一种水声块结构稀疏特性的信道估计方法
【专利摘要】本发明公开一种水声块结构稀疏特性的信道估计方法，包括三大模块：OFDM系统整体工作模块；块正交匹配追踪算法流程模块；基于块稀疏压缩感知信道估计模块。此种信道估计方法可有效降低信号稀疏度，进而降低估计所需的导频数量，大幅降低算法估计所需要的时间，且提高重构的精度。
【专利说明】
一种水声块结构稀疏特性的信道估计方法
技术领域
[0001] 本发明属于水声通信领域，涉及一种0FDM水声通信系统中块稀疏结构的信道估计 方法，特别适合于稀疏多径呈块状分布的水声信道。
【背景技术】
[0002] 随着在海洋中的活动越来越频繁，人们对水声通信速率的要求越来越高。正交频 分复用技术(CFDM)以其较高的频带利用率、低复杂性成为高速率水声通信的首选技术。然 而，(FDM对定时和频偏都比较敏感，多普勒频移会产生频率同步误差，引起载波间干扰。在 水声0FDM系统接收机设计中，往往需要已知信道状态信息(Channel State Information, CSI)，只有准确地估计出信道的时域或者频域响应，才能对接收到的数据进行校正与恢复， 以获得相干检测的性能增益。然而，实际通讯中无法获取准确的信道状态信息，因此信道估 计起到尤为重要的作用。
[0003 ]目前，水声通信中主要采取以下两类方法进行信道估计：
[0004] 1、盲信道估计。该类方法主要利用接收端接收的未知信号，通过复杂的数据统计 特性得到信道冲激响应。虽然不需要导频，但是一般都比较复杂，运算量大，不适于硬件实 现，目前还只停留在理论研究阶段。
[0005] 2、基于导频的信道估计。该类方法通过估计混合在发送数据中的导频处的信道响 应，然后利用一些插值算法拟合出其他数据处的信道响应。这种方法因其复杂度低，易实现 等优点在实际工程中得到了广泛应用。传统的方法是使用最小二乘(LS)和最小均方误差 (MMSE)等方法进行信道估计，该类方法需要在发送的数据流中插入大量的导频才能取得准 确的信道状态信息，增加了系统的开销，降低了频谱资源利用率，并且该类方法是在无线信 道路径丰富时估计性能才较好。通过物理信道测量发现，水声信道往往呈现稀疏性，即信道 的能量主要集中在少数的几条多径时延扩展和多普勒扩展上，而其它大部分的多径延迟和 多普勒扩展上的能量几乎为零。由于稀疏信道只有少数非零抽头，采样到信道零抽头时传 统的信道估计方法无法准确地估计出信道响应。压缩感知理论指出，只要信号是可压缩的 或在某个变换域下可被稀疏表示，则可以用一个与变换基不相关的观测矩阵将变换所得到 的高维信号投影到一个低维空间上并且可从这些少量的测量数据中以高概率重构出原始 信号。水声信道本身所具有的稀疏特性是使用压缩感知理论的前提，使用压缩感知进行信 道估计，不仅降低了估计所需的导频数量，而且提高了重构的精度。目前，国内外已有许多 学者使用压缩感知理论进行稀疏信道估计，并且取得了大量的研究成果。
【发明内容】
[0006] 本发明的目的，在于提供一种水声块结构稀疏特性的信道估计方法，其可有效降 低信号稀疏度，进而降低估计所需的导频数量，大幅降低算法估计所需要的时间，且提高重 构的精度。
[0007] 为了达成上述目的，本发明的解决方案是：
[0008] 一种水声块结构稀疏特性的信道估计方法，包括如下三大模块：
[0009] (1)0FDM系统整体工作模块；
[0010] (2)块正交匹配追踪算法流程模块；
[0011] (3)基于块稀疏压缩感知信道估计模块。
[0012] 上述0FDM系统整体工作模块的步骤如下：
[0013]步骤S1-1:将输入的比特序列经星座调制映射成星座点上Xk的复数点，插入导频， 得到频域信号X=[X0，Xi，.. .，Χν-ι]τ，（ ·）?示转置，其中，ke [0，Ν-1-ΝΡ]，Ν为0FDM子载波 个数，NP为插入的导频数量；
[0014] 步骤S1-2:经IFFT变换后频域信号X转变为时域信号X=[XQ，X1，. . .，XN-
[0015] 步骤S卜3:经信道并去除循环前缀后的接收信号为7 = [70,71，. . .，yN-0τ;
[0016] 步骤S1-4:接收端，信号y经FFT变换后恢复到频域信号ΥζΚο,Υι . . .
[0017] 步骤S1 -5:频域信号Υ经解调、并串变化后得到二进制的输出比特流。
[0018]上述块正交匹配追踪算法流程模块的步骤如下：
[0019] 步骤S2-1:输入:Mm X 1〇|测量矩阵Φ，测量向量y，±夬稀疏度Κ;
[0020] 输出：原始块稀疏信号的一个估计；
[0021 ]步骤S2-2:初始化:令块索引集Λ =# :，迭代次数1 = 1，初始化迭代余量r〇 = y;
[0022] 步骤S2-3:迭代：
[0023] 选择出迭代余量与测量矩阵Φ每列内积最大的块索引：
[0024]更新块索引集：八^八卜⑴匕其中山为上一步计算所得的索引号；
[0025] 得到估计信号：= ，其中，Φ、..表示从测量矩阵Φ中选择块索引集Λ :指示的 列所组成的子矩阵，Φ1,表示为~的广义逆矩阵，即Φ1, =(Φ?,Φ~)-％1,;
[0026] 更新迭代余量?,ν-Φ.、,ν，直到| |ri| |2<ε，ε为设定的迭代误差。
[0027] 上述步骤S2-3中，ε取1〇-5。
[0028] 上述基于块稀疏压缩感知信道估计模块的步骤如下：
[0029]步骤S3-1:使用选择矩阵记录导频插入位置，从发送信号中选出导频处的信号ΧΡ， 以及对应导频处的接收信号yP;
[0030] 步骤S3-2:将插入的导频符号ΧΡ与经DFT变换矩阵生成的F矩阵相乘，得到测量矩 阵T;
[0031] 步骤S3-3:将已知的测量矩阵T和接收信号yP传入块正交匹配追踪算法，求得信道 冲击响应估计向量h
[0032] 步骤S3-4:由导频处的|向量，通过H=Fh得到信道频域响应采样值Η。
[0033] 上述步骤S3-2中，测量矩阵Φ采用傅里叶矩阵、伯努利举证或部分哈达玛矩阵。
[0034] 采用上述方案后，本发明采用基于导频的信道估计方法，由于目前使用的标准压 缩感知理论是建立在稀疏模型的基础上，没有进一步考虑稀疏信号的内部结构关系，但是 若忽略信号的内在结构，运用普通的重构算法进行重构，则会造成一定的重构误差且降低 重构的效率，因此，本发明针对目前忽略的内部结构，在第二类方法中的标准压缩感知理论 的基础上，进一步考虑信号内在的结构，常用的信号结构有:聚类结构稀疏信号（又称块稀 疏），这类信号的非零元素(或者零元素)呈现聚集成块的形式;树结构稀疏信号，这一类稀 疏信号的非零元素呈树状分布，常见于小波变换系数中；联合稀疏信号:某个信号的非零元 素的位置在多次测量里均保持不变。考虑到海洋信道的水声传播过程，水声信道介质的非 均匀性等特性将导致声线以簇的形式传播。因此，本发明将稀疏信道进行分块，使用块稀疏 压缩感知理论的相关算法进行信道估计，可有效降低信道的稀疏度、降低恢复时间、提高信 道估计的性能。
【附图说明】
[0035]图1是本发明中0FDM系统整体工作模块示意图；
[0036]图2是Β0ΜΡ算法的流程图；
[0037]图3是基于块稀疏压缩感知信道估计示意图；
[0038]图4是块稀疏信道模型；
[0039]图5是稀疏信道模型。
【具体实施方式】
[0040] 以下将结合附图，对本发明的技术方案及有益效果进行详细说明。
[0041] 本发明提供一种水声块结构稀疏特性的信道估计方法，其针对水声信道的特点， 在基于标准压缩感知信道估计的基础之上，进一步考虑稀疏信号的内部结构，将稀疏信号 分块有效降低信号的稀疏度，特别适用于水声等块稀疏结构的信道中。
[0042]下面首先介绍一下0FDM和水声信道的数学模型及接收信号的相关信息。
[0043] 1.0FDM数学模型
[0044] 设0FDM子载波个数为N，第η个子载波的频率为fn = fc+n Δ f，其中Δ f = 1/T，为子载 波间隔，Τ为OFDM符号持续时间，f。为载波频率，则发送端的发射信号可表示为：
[0046]其中，x(t)为发送端要传输的数据符号，里面包含要传输的信息符号，也有用于信 道估计的导频符号；0FDM符号前缀的长度大于最大可能的路径时延，NXN矩阵X = diag(x (0)，χ(1)，...，χ(Ν-1))表示OFDM符号内的数据，包含用户数据处理映射后的信号的导频信 号。
[0047] 2.水声信道模型
[0048]水声相干多径信道的冲激响应函数可表示为
[0050]其中，L为信道多径的数目，lu(t)为t时刻第1条路径的复增益，n(t)为t时刻第1条 路径的时延。
[00511假设信道的相干时间远大于0FDM的符号周期，则在一个0FDM符号里，信道冲击响 应可以认为是时不变的，即为h⑴。对h⑴用0FDM系统的采样周期Ts?行采样，得到水声信 道的离散时间模型为：
[0053]上式(3)中，L为离散时间系统模型中抽头时延线的总个数，即信道长度。信道的稀 疏性体现在[h^lu，...，匕-^中大部分的信道抽头是0或接近于0。
[0054]块稀疏信号的定义：
[0055]假设一维离散实值信号XeRDxl由给定长度为dP(l彡ρ<Ρ)的若干块组成，用χ[ρ] 用来表示第ρ个长度SdP的子块，那么X可以被
[0056]
，其中I( I |χ[ρ] I h>0)是指标函数，当 | |x[P] I h>〇 时值为1
是标准的Edulide范数。如果信号XeR Dxl在1 = (6, . . .，dP} 上至多有kz个x[p]有非零范数，即| |x| |q,i彡kz，就称X为kz阶块稀疏的。当dP=l时，块稀疏 就和标准意义上的稀疏完全等价了。
[0057]假设上式（3)中的向量hecLxl是由C个块级联而成且每个块有d个信道抽头，则上 式(3)的块稀疏多径信道h可以重写为：
[0059] 其中，L = Cd，在C个块中，只有m(m<<C)个块内元素是非零的，块稀疏信道模型如 图4所示，图5为稀疏信道模型。
[0060] 3.接收信号
[0061] 映射后的包含导频的发射信号经过水声块稀疏信道后，则接收端得到NX 1样值向 量：
[0062] y=XH+n=XFh+n (5)
[0063] y表不NX 1向量的接收信号，X表不NX N矩阵的发送信号，η表不均值为零的复加性 高斯白噪声，F为Ν X NDFT变换矩阵的前L列：
[0066]所述信道估计方法对输入信号没有特别要求，包括如下三大模块：
[0067] (1)0FDM系统整体工作模块
[0068] 配合图1所示，包括如下步骤：
[0069] 步骤S1-1:将输入的比特序列经星座调制（16QAM)映射成星座点上Xk的复数点，k e[0，N-1-NP]，N P为插入的导频数量，插入导频，得到频域信号ΧΚΧοΑ,.,.,Χν-小，（·）τ 表示转置；
[0070] 步骤S1-2:经IFFT变换后频域信号X转变为时域信号x=[XQ，X1，. . .，χν-Τ;
[0071] 步骤S1-3:经信道并去除循环前缀后的接收信号为7 = [7()，71，. . . jh]1;
[0072] 步骤S1-4:接收端，信号y经FFT变换后恢复到频域信号ΥζΚο,Υι . . .
[0073] 步骤S1 -5:频域信号Υ经解调、并串变化后得到二进制的输出比特流。
[0074]图1为本发明的OFDM系统整体工作模块示意图，在一个OFDM系统中，二进制数据经 过16QAM调制后均匀插ANP个导频，经过步骤S1-2的IFFT变换则可得到系统的时域发射信
经0FDM系统的采样周期1^采样，得到离散发 射信号。经过水声信道后，在接收端经步骤S1-4后系统接收到的信号S:y = XH+n。
[0075] (2)块正交匹配追踪(Β0ΜΡ)算法流程模块
[0076] 在信道重构中，本发明将采用块结构稀疏压缩感知贪婪算法里的Β0ΜΡ算法，因为 该算法复杂度低，实现简单。
[0077]压缩采样的测量过程为:y=〇x
[0079] 其中，I |x| |2,〇表示非零块的个数。
[0080] 配合图2所示，该模块具体包括如下步骤：
[0081] 步骤S2-1:输入:MmXMn维测量矩阵Φ，测量向量y，块稀疏度K;
[0082]输出：原始块稀疏信号的一个估计；
[0083] 步骤S2-2:初始化:令块索引集A = p 迭代次数1 = 1，初始化迭代余量r〇 = y;
[0084] 步骤S2-3:迭代：
[0085] 选择出迭代余量与测量矩阵Φ每列内积最大的块索引
[0086] 更新块索引集，其中^为上一步计算所得的索引号；
[0087] 得到估计信号：= Φ.Ι，，其中，ΦΑ?表示从测量矩阵Φ中选择块索引集Λ :指示的 列所组成的子矩阵，Φ〗,表示为的广义逆矩阵，即Φ、. = (Φ:^ΦΛ?)
[0088] 更新迭代余量-Φν.ν，直到I |ri| |2<ε(ε为设定的迭代误差，通常取ΚΓ5)。
[0089] (3)基于块稀疏压缩感知信道估计模块
[0090] 将0FDM系统形成的信号，利用块稀疏压缩感知的Β0ΜΡ算法进行信道估计，得到信 道响应。
[0091] 信号的稀疏表示是压缩感知理论的一个重要前提，直接关系到原始信号的重构精 度。只有选择一个合适的正交基函数对信号进行稀疏分解，才能保证原始信号的精确恢复。 在压缩感知理论中，若信号稀疏度为κ，φ是1 XMn的测量矩阵，将X投影到Φ上得y = Φ X，要 准确地重构信号χ，Φ必须满足受限等距特性(RIP)准则，即存在一个很小的常数Ske(〇，i) 使得对任意k阶稀疏量X有式(7)成立：
[0093]常用的测量矩阵有傅里叶矩阵、伯努利举证、部分哈达玛矩阵等。
[0094]本发明采用CS进行信道估计，是考虑到基于CS的信道估计方法所需导频数较少， 且估计精确，适合本发明所述的信道估计方法。基于导频的信道估计问题就相当于已知发 送导频信号构成的矩阵X和相应的导频接收信号矩阵Y，根据某种准则求解信道冲激响应H。 本发明采用时域方法进行信道估计，首先估计出信道冲击响应，再进行傅里叶变换得到信 道频域响应。
[0095]由上式(5)看出，接收端的接收信号为y=XH+n=XFh+n。
[0096] 设S为NPXN的选择矩阵，从N个子载波中选择出NPf导频所插入的位置，S矩阵是从 NXN单位矩阵中选择与导频位置对应的^行，接收端对应导频位置的接收信号可以表示 为：
[0097] yp = XpFph+np = Th+nP (8)
[0098] 其中，yP = Sy为接收的导频信号构成的NP X 1的向量;XP = SXS '为NP X NP的导频信 号;FP = SF为NP X N的DFT矩阵;nP = Sn为导频位置对应的噪声组成的NP X 1向量。
[0099] 由于h具有块结构稀疏性，将信道估计问题转化为块稀疏信号的重构问题。经步骤 S3-2,即可以把^看成是测量矩阵T = XPF对信道h的测量矩阵T也按h的方法进行分块），其中 yP，XP，F P，办对于信号接收端均为已知信号，只有h是未知的，接收端经步骤S3-3采用块正交 匹配追踪(Β0ΜΡ)重构算法恢复h向量，信道频域响应Η可以由H=Fh得到。
[0100] 上述过程的具体步骤如下：
[0101] 步骤S3-1:使用选择矩阵S记录NP导频插入位置，从发送信号中选出导频处的信号 Xp，以及对应导频处的接收信号yp;
[0102] 步骤S3-2:将插入的导频符号ΧΡ与经DFT变换矩阵生成的F矩阵相乘，得到测量矩 阵T = XPF;
[0103] 步骤S3-3:将已知的测量矩阵T和接收信号yP传入块正交匹配追踪(Β0ΜΡ)算法，求 得信道冲击响应估计向量L
[0104] 步骤S3-4:由导频处的S向量，通过H=Fh得到信道频域响应采样值H。
[0105] Β0ΜΡ算法是标准压缩感知理论的0ΜΡ算法的推广，0ΜΡ算法在迭代时每次选出与迭 代余量最匹配的一列的索引，而Β0ΜΡ算法每次选出与迭代余量最匹配的一块的索引。假设 块稀疏信号每块的抽头数为d，非零抽头块的个数为p(即非零抽头数为dp)，Β0ΜΡ算法需要 迭代P次可得到估计信号，0ΜΡ算法至少需要迭代dp次才能得到估计的信号。因此在考虑了 信号的块稀疏先验信息后可有效降低稀疏度，降低算法的运行时间，同时可以进一步减少 导频的使用。
[0106]综合上述，本发明一种水声块结构稀疏特性的信道估计方法，针对目前的水声信 道估计方法主要是基于信道的稀疏结构模型的思想，本发明在稀疏模型基础上进一步考虑 水声稀疏信号内部的块状分布结构，即稀疏信号的非零元素(或者零元素)不是随机分布而 是以块的形式呈现，将稀疏信号分块可有效降低信号稀疏度，进而降低估计所需的导频数 量，大幅降低算法估计所需要的时间，且提高重构的精度。
[0107]以上实施例仅为说明本发明的技术思想，不能以此限定本发明的保护范围，凡是 按照本发明提出的技术思想，在技术方案基础上所做的任何改动，均落入本发明保护范围 之内。
【主权项】
1. 一种水声块结构稀疏特性的信道估计方法，其特征在于包括如下三大模块： (I )OFDM系统整体工作模块； (2) 块正交匹配追踪算法流程模块； (3) 基于块稀疏压缩感知信道估计模块。2. 如权利要求1所述的一种水声块结构稀疏特性的信道估计方法，其特征在于：所述 OFDM系统整体工作模块的步骤如下： 步骤Sl-I:将输入的比特序列经星座调制映射成星座点上Xk的复数点，插入导频，得到 频域信号X= [Xo，Xi，...，Xn-I]τ，（ ·广表示转置，其中，k e [〇，N-I-Np]，N为OFDM子载波个数， Np为插入的导频数量； 步骤S1-2:经IFFT变换后频域信号X转变为时域信号X= [XQ，X1，...，Xn-Ot; 步骤S1-3:经信道并去除循环前缀后的接收信号为7=[列，71，. . .，ynf; 步骤S1-4:接收端，信号y经FFT变换后恢复到频域信号Y= [YoJ1, ... ,Yn-Ot; 步骤S1-5:频域信号Y经解调、并串变化后得到二进制的输出比特流。3. 如权利要求1所述的一种水声块结构稀疏特性的信道估计方法，其特征在于:所述块 正交匹配追踪算法流程模块的步骤如下： 步骤S2-1:输入:Mm X 1〇|测量矩阵Φ，测量向量y，块稀疏度K; 输出：原始块稀疏信号的一个估计； 步骤S2-2:初始化:令块索引集A =供，迭代次数I = 1，初始化迭代余量r〇 = y;步骤S2-3:迭代： 选择出迭代余量与测量矩阵Φ每列内积最大的块索引： 更新块索引集：A1= Ap1U ii，其中，ii为上一步计算所得的索引号； 得到估计信号3 = 其中，φΛ,表示从测量矩阵Φ中选择块索引集A1指示的列所 组成的子矩阵，表示为ΦΑ?的广义逆矩阵，艮更新迭代余量Κ=)·-ΦΛι，直到I Ir1I |2<ε，ε为设定的迭代误差。4. 如权利要求3所述的一种水声块结构稀疏特性的信道估计方法，其特征在于:所述步 骤 S2-3中，ε取10-5。5. 如权利要求1所述的一种水声块结构稀疏特性的信道估计方法，其特征在于:所述基 于块稀疏压缩感知信道估计模块的步骤如下： 步骤S3-1:使用选择矩阵S记录导频插入位置，从发送信号中选出导频处的信号ΧΡ，以及 对应导频处的接收信号yP; 步骤S3-2:将插入的导频符号Xp与经DFT变换矩阵生成的F矩阵相乘，得到测量矩阵T; 步骤S3-3:将已知的测量矩阵T和接收信号yP传入块正交匹配追踪算法，求得信道冲击 响应估计向量 步骤S3-4:由导频处的I向量，通过H=Fh得到信道频域响应采样值H。6. 如权利要求5所述的一种水声块结构稀疏特性的信道估计方法，其特征在于:所述步 骤S3-2中，测量矩阵Φ采用傅里叶矩阵、伯努利举证或部分哈达玛矩阵。
【文档编号】H04L27/26GK106027445SQ201610318947
【公开日】2016年10月12日
【申请日】2016年5月13日
【发明人】王彪, 朱芹
【申请人】江苏科技大学