用于高相对带宽通信系统时变信道的多采样率自适应均衡技术的制作方法
【专利摘要】名称为用于高相对带宽通信系统时变信道的多采样率自适应均衡技术，所属技术领域为通信工程。本项技术是提供一种用来处理高相对带宽通信系统时变信道影响的自适应滤波技术。通过反复“刷式”更新过程，在Turbo软解码器的辅助下，对相位和滤波权重进行调整。重采样(Resampling)结构对于信道矩阵的特征值(Eigenvalue)具有向“1”处收敛的作用，明显改善了信道矩阵的条件数(condition number)。经过滤波后，观察滤波符号的星座图可以看出对于BPSK的发送符号，带有重采样结构的结构可以较为清楚的恢复出BPSK的星座分布，其中重复刷式更新可以进一步优化滤波效果。可以用于推动海底石油勘探、潜艇水下通信、水下传感器远程采集、海洋污染监测、海岸水下监测网、水下无人机器人控制等领域的发展。
【专利说明】用于高相对带宽通信系统时变信道的多采样率自适应均衡 技术 所属技术领域
[0001] 通信工程。
【背景技术】
[0002] 无线通信系统大多采用低相对带宽系统，即通信系统的频率绝对带宽与载波的中 心频率的比值较低（例如< 10% )。比如，当前广泛采用的基于IEEE 802. 11标准的WiFi 无线局域网络设备的中心频率为2. 4GHz，而其标准基带带宽只有W = 20MHz，其相对带 宽只有大概λ =W/f；= 0.83%。与之相对的，高相对带宽的通信系统的绝对带宽与载波 的中心频率的比值一般超过25%，比如追求高分辨率的超宽带雷达系统。其他高相对带宽 的通信系统还包括军事领域的超/极低频无线电通信系统，和远距离水下声波通信系统。 例如采用4kHz~8kHz频段的水声通信系统，W = 4kHz，fe= 6kHz，则相对带宽高达λ = 66%，是一个典型的高相对带宽通信系统。
[0003] 高速运动的通信终端将对通信系统都带来多普勒效应，从而带来信道的时变性。 在低相对带宽通信系统中，多普勒效应表现为频率偏移，对于这种频率偏移的补偿技术已 经得到了充分的研究和成功案例，例如对于基于DVB-Η技术的移动电视广播系统、4G移动 手机通讯系统等。然而，高相对带宽通信系统如何支持高速运动的通信终端还是一个亟待 解决的课题。已经得到证实，在高相对带宽通信系统中，多普勒效应反应的是波形伸缩因 子。
[0004] 首先，使用数学方式描述这种伸缩因子的存在。任何时变通信系统都可以描述 为：
[0005]
[0006] 其中为接收信号，只〇为发送信号，河〇为背景噪声，而时变信道的冲击响应 可以具体表述为
[0007]
[0008] 这里的^是第1条路径的衰减，而τ i (t)代表着伴随第1条路径的随着时间变化 的延迟信息。
[0009] 如果我们假设，发送端沿着第1条路径为径向速度分量为Vl(T)，而接收端沿着第1 条路径为径向速度分量为 Vl(R)，那么我们可以具体表述τια)为
[0010]
[0011] 这里与时间无关的τ i表示第1条路径的初始距离决定的延迟常量，而c则是传 播速度，例如水声传播速度大概为c = 1500m/s。
[0012] 按照多普勒现象的物理性质，我们定义与该第1条路径对应的多普勒伸缩因子为
[0013]
[0014] 例如，对于静止发射段（即Vl(T)= 0)和20节径向时速的接收端（即大约v /R) = 10. 28m/s)，将会最大产生伸缩因子约α 1= 1. 0069。至此我们可以修改时变的时延信息为
[0015]
[0016] 进而将它代入信道的冲击响应h(t，τ)中，我们可以修改一般性的信道模型，重 新表达下述的信道模型，表达如下
[0017]
[0018] 因此可以看出伸缩因子（^对于信号s(t)产生了波形伸缩变形作用。同时由于多 条通信路径，这个伸缩因子在不同路径上反映为不同的数值，从而变现为"多伸缩多时延" 特点。对于这种时变信道，现在的均衡技术主要可以分为三类。
[0019] (1)第一类将"多伸缩"时变性使用频率偏移的近似。但实际上，使用频率偏移对 多普勒效应进行描述，是低相对带宽通信系统的一种合理近似。但是对于高相对带宽通信 系统，这种做法只是一种不合理的近似方式。
[0020] 首先，从通信原理的数学基础上阐明这种做法的不合理性。在忽略加性噪声时，上 面带有"多伸缩多时延"信道模型的等价频率域的表达式如下：
[0021]
[0022] 这里α1= l+β i，用/)为F⑷的傅里叶变换，歹(/)为印)的傅里叶变换。如果 定义相对径向速度％ = 4R) -<Τ)，根据传播速度c 一般远大于通信终端的速度，我们有 那么根据
[0026] 注意，这里的频率f在一定带宽内，即f e [flOT，fhigh]，而带宽W = fhigh-flOT，中心
[0023]
[0024]
[0025] 频率/ = /high + L，则相对带宽为λ =w/f。。这样频率f的动态范围限制为 c 2
[0027]
[0028] 注意，对于低相对带宽系统（即λ = W/f；<< 1)，变量β f~β f。，因此才允许 将通信系统公式近似地表述为
[0029]
[0030] 其中频率偏移为ω1= β?·。。由上式可以看出，原有的多普勒伸缩因子对于低相 对带宽系统的信号s(t)不再具备伸缩作用，而是表现为频率偏移。
[0031] 然而，对于高相对带宽系统，不可能将β?·近似为β?·。，无法获得上面频率偏移的 近似公式。在直观上，我们假设一个只有两个路径的时变信道，并将具有高相对带宽的信号 波形通过时变信道，其中第1信道带有伸缩因子 α :=0.99而第1信道带有伸缩因子a i =1. 01。构造该信道的频域信道矩阵如下图1 (a)所示，多普勒效应作用于信号后体现为明 显的变形。但是，使用低相对带宽系统的频率偏移来近似的话，如图1(b)显示的是频率的 平行移动。因此，两者的区别很明显，即在高相对带宽通信系统中使用频率偏移来近似多普 勒效应是不合理的。
[0032] 第二类是使用"单伸缩因子"近似"多伸缩"时变性。这类方法使用单一的伸缩 因子和载波频率偏移（CF0)对"多伸缩"时变信道进行近似。例如近几年A.-B. Salberg 等在 IEEE Transactions on Wireless Communications 发表的学术论文〈〈Doppler and frequency-offset synchronization in wideband OFDM〉〉，B. Li 等在 IEEE Journal of Oceanic Engineering 发表的学术论文〈〈Multicarrier communication over underwater acoustic channels with nonuniform Doppler shifts)), S. Yerramalli 等在 IEEE Transaction on Wireless Communications 发表的〈〈Optimal resampling of OFDM signals for multiscale-multilag underwater acoustic channels〉〉，这些前沿的石开究工 作的共同点是使用单一的伸缩因子的处理方式对高相对带宽通信系统的时变信道进行处 理。
[0033] 第三类并不对"多伸缩"时变性进行近似，而是正视这种时变性从而提出处理方 式，但是这些工作还需要假设信道信息的已知性，例如Y. Jiang等在IEEE Transactions on signal processing 上发表的学术论文〈〈Discrete time-scale characterization of wideband time-varying systems》和 A. R. Margetts 等在 IEEE Transactions on Wireless Communications 发表的学术论文〈〈Joint scale-lag diversity in wideband mobile direct sequence spread spectrum systems》在理论上对使用"多伸缩"可以带来的优 势做出了阐述。
【申请人】在IEEE Transactions on signal processing上发表的学术论文 ?Multi-Rate Block Transmission Over Wideband Multi-Scale Multi-Lag Channels)) 和在EURASIP Journal on Wireless Communications and Networking上发表的学术论文 《Iterative equalization for OFDM systems over wideband Multi-Scale Multi-Lag channels》利用"多伸缩"的优势对信道均衡技术做出了设计，但这些工作都需要假设信道 信息的已知性，或者对于这种多伸缩信道进行预先估计（这种估计十分困难），因此实用性 不强。
[0034] 针对以上的问题，本发明设计一种用于高相对带宽通信系统时变信道的多采样率 自适应均衡技术。首先，多伸缩多时延的时变性被正视，不再使用频率偏移或单一伸缩性进 行近似，同时回避对于信道信息的预知性，这样既能够提高信道均衡的性能，又能提高信道 均衡的实用性。
【发明内容】
[0035] 首先介绍发射信号符号具体表述为
[0036]
[0037] 其中N(如N= 256)为数据块长度,f。为中心频率，而sn为数据序列的符号或者训 练序列的符号。若乂为数据序列的序号集合，4为训练序列的序号集合，那么当n e仏时Sn 为数据序列符号，当neArp时snS训练序列符号。这里所有符号采用binary phase-shift keying(BPSK)，即
[0038]
η 、 '
[0039] 数据段以较长的训练序列片段为起始，之后由训练序列和数据序列交叉而成。在 上面的例子中，数据序列符号的数目为N D= 256。
[0040] 将高相对带宽通信的接收信号按照"多伸缩多时延"系统的方式进行表示如下
[0041]
[0042] 其中Q为具有不同伸缩因子（即，当1古Γ时&1古&1,)的通信路径数目，而L 为具有不同伸缩因子（即，当m古m'时^古τ 的通信路径数目。
[0043] 其中p (t)是与发射信号J(〇具有同等带宽（即1/TS)的低通滤波器，而;是在第 q层（q = 0,1，. . .，Q-1)上通过采样速率17&|11采集的第η个采样数据。
[0044] 下面对某一层的滤波过程进行描述（省略层数上标号（q)，这里的q = 0, 1，...， Q-1)如下：
[0045]
[0046]
[0047]
[0048]
[0049]
[0050]
[0051] 首先按照RLS滤波器的阶数L将第η个基带数据符号yn前后邻近的L个符号合 并为一个数组，即
[0052]
[0053] 这里上角标T表示矩阵转置。而η的计数是η = 0, 1, ···, N-1,这里的N是一个 数据块长度，锁相环PLL与滤波器更新的运算如下
[0054]
[0055] 其中\ = 。而这里的u对应在第p次扫描过程中得到的对第η个发送基 带数据符号的估计，而θ η, Ρ是第Ρ次扫描过程中PLL的相位输出，滤波器的L阶权重组合 为cn,p=c。Cl…CuT。注意此处提到的扫描过程，指的是滤波器在Ν长的数据块上进行向 前和向后的滤波过程，比如P = 1，2, ... P。即当ρ为奇数，滤波的顺序沿着η = 0,1，...， Ν-2的向前扫描，则土中取-号；而当ρ为偶数时，滤波的顺序沿着η = Ν-1，Ν-2，. . .，1的 向后扫描，则土中取+号。这样的扫描的目的主要通过增加滤波在训练序列上的作用次数 从而提高滤波稳定性。比如Ρ = 3,则扫描为三次。因此，这里的更新过程是一种反复"刷 式"更新过程。
[0056] 判决器（decision mode)的工作方式如下：首先根据符号位置η是否属于训练序 列和t的情况，决定参考符号#并计算误差、=d，P。具体的讲，
[0057] 第一步是根据训练序列计算出滤波后的符号均值μ和均方差σ 2,即
[0058]
[0059] 并通过他们计算出的对数似然比（Log-Likelihood Ratio, LLR)输入数为
[0060]
[0061] 而LLR运算模块输出的llr，被SIS0 decoder (软输入软输出解码器）用来计算 比特概率值如下
[0062]
[0063]
[0064] 其中前者指7",p l6"=Q和7",p L"=1分别代表估计发射机发送的当前位置的比特是0或 者1的概率。
[0065]
[0066] 其中Γ为概率门限值（取Γ =80%)，而sgn{}指取数值的符号（+1或-1)。这 样可以计算出误差^ = C k
[0067] 第二步，根据误差κ n,p来更新滤波器的增益权重数组c n,p，具体如下
[0068]
[0069]
[0070]
[0071]
[0072] cn>p= cn±1>p+K n>pgn>p
[0073] 其中土按照扫描过程的向前或向后的滤波过程来选择-号或+号，遗忘因子λ =0. 99。
[0074] 第二步，对于锁相环DPPL相位Θ ^的更新采取以下的过程：
[0075]
[0076] ηη±1>ρ= η η,ρ+(-ι)ρΘη,ρ
[0077] θη±1>ρ= θ η ρ+ΚιΘη>ρ+(-1)ρΚ2ηη±1>ρ
[0078] 其中1和Κ 2分别指分式和积分相位跟踪（proportional and integral phase-tracking)常数（如 1= 0. 02 和 K 2= 0. 04)，而 η n,p初始值设为零。
[0079] W下而描沭的高相对带甯的诵#r系统为例，
[0080]
[0081] 其中伸缩因子aq均匀的分布在[0.99,1.01)范围上。对于水声通信而言，这样 的伸缩因子对应最大约15m/s的径向相对速度，而时延T q均匀的分布在[0,200ms)范围； 信道增益\服从独立分布正态随机变量，其分布均值为0,而分布均方差为 6_%#。信道多径 数目取Q = 5。背景加性噪声坷0为白噪声。
[0082] 首先，对上述信道进行分析发现，重采样（Resampling)结构对于信道矩阵的特征 值（Eigenvalue)具有向"1"处收敛的作用，与没有Resampling前端处理的处理机相比，如 图2(a)所示，明显改善了信道矩阵的条件数（condition number)。图2(b)所示，对于这样 应用，采用了 Resampling方式的明显提高了后端自适应滤波的收敛速度，大概只需要不采 用重采样（Resampling)方式的一半循环次数，从而大幅度减少了处理时间。可以说明，采 用Resampling方式，是高速运动终端进行水声通信的必要处理结构。
[0083] 经过滤波后，观察在不同方式下滤波符号的星座图如图3。可以看出对于BPSK的 发送符号，没有任何重采样结构的滤波过程无法恢复出BPSK的星座分布，如图3(a)所示； 带有重采样结构的结构可以较为清楚的恢复出BPSK的星座分布，其中重复刷式更新可以 进一步优化滤波效果，例如比较下图3(b)和（c)。
[0084] 本发明提出的自适应滤波的效果还可以从图4进一步得到认证。从图4显示测试 结果中可以看出采用SIS0(软输入软输出）的软（Soft)解码方式，性能优于硬（hard)解 码方式，而且滤波的性能随着刷式更新次数的提高而有所优化。
[0085] 对本发明的自适应滤波结构做出下面简要说明。
[0086] 其中输入为经过"多伸缩多时延"信道的被接收的高相对带宽信号印)。对于此接 收信号，首先在第q层上根据aqf。的频率和p (a qt)(这里q = 0, 1，. . .，Q-1)进行下变频和 低通滤波，并根据Ts/aq的采样率进行采样获得第q层的基带信号采样值。
[0087] 对于采样值2^首先通过6-?进行相位的调整，而后面的滤波器（Filter)使用权 重沙)按照递归式最小均方（RLS)算法进行滤波，得到妒之后取和运算获得i 。 η,ρ η,ρ η,ρ / -J η,ρ 注意，这里滤波器（Filter)的更新和锁相环（PLL)相位的更新都需要根据判决器 (decision mode)而进行的。将U安照ND(如256)的长度形成向量对数似然比 (Log-Likelihood Ratio, LLR)运算结构根据计算出并且根据软解码器（SIS0 decoder)获得概率Yn,p|b_=。和γ n,p|b_=1。该概率将参与判决器（decision mode)。此 外，软解码器（SIS0 decoder)还会解算处最终的比特估计值。
【附图说明】
[0088] 图1是高、低相对宽带使用频率偏移来近似多普勒效应结果示意图，图l_(a)为高 相对带宽使用频率偏移近似多普勒效应结果，图l-(b)为低相对带宽使用频率偏移近似 多普勒效应结果； 图2是采用Resampling方式，是高速运动终端进行水声通信的必要处理结构示意图， 图2-(a)是水声信道矩阵的特征值，图2-(b)是自适应滤波的收敛； 图3是经过滤波后，观察在不同方式下滤波符号的星座图示意图，图3-(a)是无重采样 结构，图3-(b)是带有重采样结构、自适应滤波、无刷式更新，图3-(c)是带有重采样结构、 自适应滤波、P= 10刷式更新； 图4是采用SIS0 (软输入软输出）的软（Soft)解码方式，性能优于硬（hard)解码方 式，而且滤波的性能随着刷式更新次数的提高而有所优化的示意图； 图5是对本发明的自适应滤波结构示意图； 图6是bn是发送的二进制比特数值（0或者1)，发送数据块的结构示例图； 图7是高相对带宽通信的接收信号对于该信号的接收前端"重采样"构架图； 图8是一种带数字锁相环DPLL的Turbo式解调器，它具有Q层处理结构，每层都是单 独的进行自适应滤波，而对每层滤波得到的结构进行整合后进行统一的判决； 图9是反复"刷式"更新过程示意图。 【具体实施方式】 水声通信信号是典型的高相对带宽通信时变系统。特别是具有高速运动通信终端的水 声通信系统将面临的是典型的"多伸缩多时延"的时变信道。因此本发明可以用于推动海 底石油勘探、潜艇水下通信、水下传感器远程采集、海洋污染监测、海岸水下监测网、水下无 人机器人控制等众多重要领域的发展。
【主权项】
1. 一种时变信道的均衡技术，其特征是，对于接收到的信号进行多速率（如Q= 3,即 三种速率a。，ai，a2)的均衡运算，即首先进行对其进行不同速率的"重新采样"，然后按照 镂相环原理分别进行相位调整，再根据递归最小二次方（RLS)算法分别进行自适应滤波运 算，并且这个自适应滤波运算的迭代更新是依据判决反馈的方式进行控制的。2. 根据权利要求1所述的时变信道的均衡技术，其特征是，所述的"重新采样结构"包 括对射频信号的中间频率的变速处理以抵消载频频偏的问题，然后通过变速后的P(a qt)匹 配滤波和Ts/aq采样周期实现对基带信号的重新采样。3. 根据权利要求1和2所述的时变信道的均衡技术，其特征是，所述的"相位调整"和 "自适应滤波器"是针对每路重新采样结构都单独配置一套，它们的更新过程由后端的判决 器（decision mode)进行分别控制。4. 根据权利要求1，2和3所述的时变信道的均衡技术，其特征是，所述的自适应迭代更 新的判决器是依据Turbo解码器计算出的"软信息"进行判决的，即似然度大于某个门限时 进行更新，反而不进行跟新。5. 根据权利要求1，2, 3和4所述的时变信道的均衡技术，其特征是，所述的多速率的均 衡运算，在只存在一种时变因子速率（如Q= 1，即一种速率a。）时，该技术成为单速率的均 衡运算，则均衡技术简化为单路运算，而后续结构不变。6. 根据权利要求1，2, 3,4和5所述的时变信道的均衡技术，其特征是，该技术更加适用 于高相对带宽通信系统，但对于相对带宽较小的传统窄带通信系统，只不过区别在于所估 计的a q都很小。
【文档编号】H04L25/03GK105897625SQ201410201994
【公开日】2016年8月24日
【申请日】2014年5月10日
【发明人】许涛
【申请人】许涛, 曹义铭