申请/专利权人：东南大学

申请日：2017-07-10

公开（公告）日：2020-01-07

公开（公告）号：CN107171986B

主分类号：H04L25/02(20060101)

分类号：H04L25/02(20060101)

优先权：

专利状态码：有效-授权

法律状态：2020.01.07#授权;2017.10.17#实质审查的生效;2017.09.15#公开

摘要：本发明公开了一种适用于多普勒失真水声信道的估计方法，该方法在鱼群算法的基础上，以迭代的方式分离多径分量，并在子迭代中自适应的调整人工鱼的位置和步长。本发明提出的改进鱼群算法利用水声信道固有的稀疏特性，显著降低了多扩展多时延信道的估计复杂度；仿真结果表明，IAFSA可以准确估计出每一条路径的参数对，在估计精度和计算复杂度上较正交匹配追踪算法OMP算法均有显著提升。

主权项：1.一种适用于多普勒失真水声信道的估计方法，其特征在于，包括如下步骤：1在问题空间中初始化鱼群位置，计算相应的适应度值，并将群体中最优适应度值和对应的位置记录在公告板上，进入子迭代过程；问题空间即为路径参数可能的取值空间，包括时延和多普勒扩展因子的取值范围，最大时延扩展为训练序列的时间长度，最大多普勒扩展为收发端最大相对运动速度与声波在海水中的速度的比值；2每一条人工鱼在其视野范围内执行聚群和追尾行为或觅食行为，更新自身位置和适应度值并更新公告板；3当子迭代次数大于设定值的一半时，若公告板中最优适应度值大于设定阈值且不发生变化，则将一半的人工鱼位置设置为最优适应度值对应的位置；4循环执行子迭代过程并不断调整步长，直至达到最大子迭代次数；5从公告板中得到最优位置，作为一条路径的参数，得到相应的信号分量，用以更新残余信号，进入下一次迭代。

全文数据：一种适用于多普勒失真水声信道的估计方法技术领域[0001]本发明涉及水声通信信道估计算法技术领域，尤其是一种适用于多普勒失真水声信道的估计方法。背景技术[0002]水声信道中显著的多普勒效应和严重的多径扩展给高速稳定的通信带来了很大的挑战。在水声系统中，声波的传输速度为1500ms，远远低于陆地无线通信中电磁波的传播速度。因此，收发端移动引起的多普勒效应十分显著，表现为引起时域上信号的压缩或扩展。因此，多普勒效应被处理为多普勒扩展因子。另一方面，严重的多径效应是由水下环境中大量的反射导致。由于声波传播速度慢，多径延时大，造成了严重的符号间干扰。为充分了解水声信道特点并克服其带来的挑战，对水声信道精确的建模和估计十分重要。[0003]如很多实验观察到的一样，不同路径的信号经历不同的多普勒扩展，在不同的时间点到达且具有不同的能量，接收信号是这些不同路径信号的叠加。所以多扩展多时延Multi-scalemulti-lag，MSML信道模型能够较好的描述水声信道的特点，为很多文献采用。根据MSML信道模型，每一条路径可以被参数化为多普勒扩展因子、时延和幅度三个参数。然而，严重的多径效应使得MSML信道的估计过于复杂。为了克服这个困难，很多研究者提出利用水声信道的稀疏特性，即大部分信道能量集中在较小的范围内。所以，MSML信道模型中，只有较少的抽头系数是非零的，需要被估计出来。因此，计算复杂度可以显著降低，并且很多利用信道稀疏特性的压缩感知算法得到了应用。[0004]基于压缩感知的算法主要分为两类：动态规划方法，如匹配追踪（Matchingpursuit，MP;线性规划方法，如基追踪Basispursuit，BPC3BP算法较高的计算复杂度限制了其应用，而MP算法得到了较为广泛的应用且出现了很多改进算法。[0005]MP算法通过迭代选取字典中与接收信号相关性最大的列来进行信道估计，并且在每次迭代结束时，从接收信号中减去相应的估计分量。在此基础上，通过使剩余信号与已选出的每一列正交，提出了正交匹配追踪（Orthogonalmatchingpursuit，0MP算法，OMP算法具有更优的估计精度和收敛速度。同时，也有一些算法提出自适应估计路径数，如稀疏自适应匹配追踪（Sparsityadaptivematchingpursuit，SaMP算法和自适应步长SaMP算法。进一步，为了降低计算量，有文献提出使用快速傅里叶变换简化OMP算法，但该方法降低的计算量有限因为其并没有改变字典本身的大小。另一种降低计算量的方法是分步估计时延和多普勒扩展，该方法仅适用于各条路径的多普勒扩展相差较小且经过了粗补偿的情况。[0006]因此，MP算法及其改进算法的不足之处在于其估计精度依赖于字典的大小，估计精度越高则字典的列数越多，因而计算量也就更大。对于时延-多普勒扩展较大的水声信道，MP算法的计算复杂度限制了高精度的参数估计，因而本发明提出了一种可以降低复杂度同时有较高的估计精度的估计算法。发明内容[0007]本发明所要解决的技术问题在于，提供一种适用于多普勒失真水声信道的估计方法，能够有较快的收敛速度和较高的估计精度。[0008]为解决上述技术问题，本发明提供一种适用于多普勒失真水声信道的估计方法，包括如下步骤：[0009]1在问题空间中初始化鱼群位置，计算相应的适应度值，并将群体中最优适应度值和对应的位置记录在公告板上，进入子迭代过程；[0010]2每一条人工鱼在其视野范围内执行聚群和追尾行为或觅食行为，更新自身位置和适应度值并更新公告板；[0011]3当子迭代次数大于设定值的一半时，若公告板中最优适应度值大于设定阈值且不发生变化，则将一半的人工鱼位置设置为最优适应度值对应的位置；[0012]⑷循环执行子迭代过程并不断调整步长，直至达到最大子迭代次数；[0013]5从公告板中得到最优位置，作为一条路径的参数，得到相应的信号分量，用以更新残余信号，进入下一次迭代。[0014]优选的，步骤⑴中，问题空间即为路径参数可能的取值空间，包括时延和多普勒扩展因子的取值范围，一般认为最大时延扩展为训练序列的时间长度，最大多普勒扩展为收发端最大相对运动速度与声波在海水中的速度的比值。[0015]优选的，步骤⑴中，人工鱼p的适应度值的计算公式为：[0017]其中r⑴为接收信号，s⑴为训练序列，Xp为人工鱼p的位置，为以Xp为时延-多普勒参数得到的训练序列。[0018]优选的，步骤2中，觅食行为是:人工鱼P在其视野范围内随机选取一个位置，若该位置的适应度值大于当前位置的适应度值，则向该位置移动一步;否则继续尝试，若尝试次数大于设定的最大值仍未成功，则随机移动一步。[0019]优选的，步骤2中，聚群行为是:人工鱼P在其视野范围内有Q个同伴，若Q0,计算Q个同伴的中心位置X。和相应的适应度值y。，若yc7QHyP，其中λ为拥挤度因子，则p向X。移动一步;若ycQ彡Xyp或Q=O，则执行觅食行为。[0020]优选的，步骤2中，追尾行为是:人工鱼p在其视野范围内内有Q个同伴，若Q0，找到具有最优适应度值的同伴Xq，若其适应度值yq满足yqQHyP，则p向乂0，计算Q个同伴的中心位置X。和相应的适应度值y。，若yjQHyp，其中λ为拥挤度因子，则ρ向X。移动一步;若ycQ彡Ayp或Q=O，则执行觅食行为。[0047]步骤2中，追尾行为是:人工鱼ρ在其视野范围内内有Q个同伴，若Q0，找到具有最优适应度值的同伴Xq，若其适应度值yq满足yqQHyP，则ρ向Xt^动一步，若yqQyP，则该鱼将向Xv，移动一步，即：[0073]其中Δ是步长，这个过程将重复I次直到有一个Xv满足要求；否则，该人工鱼将在视野范围内随机选取一点。[0074]人工鱼的聚群行为：[0075]令^为人工鱼P的当前位置，其视野范围内有Q个同伴，如果Q0,计算这Q个同伴的中心位置：[0077]定义λ为拥挤度因子，如果ycQHyp，则人工鱼p将会向Xc移动一步;否则，将执行觅食行为。若Q=O人工鱼也将执行觅食行为。[0078]人工鱼的追尾行为：[0079]人工鱼p的视野范围内有Q个同伴，如果Q0,找到具有最大适应度值yq的同伴Xqc3SyqQHyP，人工鱼p将向乂1移动一步，若yqQkmax2时，如果公告板保持不变且ycipte，将一半的鱼位置调整为X〇pt。[0091]⑶设置k=k+l，并调整步长为.，跳转到步骤3循环执行，直至kkmax。[0092]7从公告板选择最优位置1^作为路径1的时延和多普勒因子估计值，得到相应的时延-多普勒训练序列si，和最优适应度值7_作为路径1的幅度估计值i,，更新剩余信号：[0094]⑶如果I=L，停止迭代;否则，1=1+1，跳至步骤1.[0095]输出：[0096]估计参数对匕。[0097]注:路径数L可以在信号同步阶段获得；阈值ε根据接收端所能够检测到的信号能量值来设定。[0098]图2—图7给出了不同信道条件下，多普勒扩展因子估计的归一化均方误差、时延估计误差和剩余信号能量比随信噪比的变化而变化的仿真曲线，并与OMP算法做了比较。其中，信道1的参数设置为:路径数L=10,各路径信号的到达时间随机分布在0〜25ms，且将最小路径时延设为0。归一化路径幅值均匀分布，多普勒扩展因子随机分布在[1，1.02]，精确至IJ4位小数。采用长度为511的伪随机序列作训练序列，且用二进制相移键控调制。载波频率为10kHz，采样率为20kHz。对于OMP算法，所构造的字典多普勒因子分辨率为IXΚΓ4,时延分辨率为0.Ims，多普勒扩展为0.02，时延扩展为25ms，这也是IAFSA的问题空间。[0099]IAFSA的参数设置为：鱼群大小为50，拥挤度因子为0.3，视野范围为[0.005，I.Oms]，初始步长为0.2，最大子迭代次数等于10，最大尝试次数等于10，阈值ε=0.2。[0100]信道2采用BELLHOP产生:水深为100m，收发端水平距离2000m，发射端固定在80m深处，接收端位于50m深处，以15ms的水平速度向发射端靠近，声速设定为1500ms。海面和海底的反射系数分别为-0.9和0.7,声线图如图1所示。IAFSA的参数与信道1中相同，仅将问题空间中多普勒扩展改为0.01。[0101]从仿真图可见，在所有的实施例中本发明的性能都明显优于OMP算法。在计算复杂度上:设训练序列长度为Kl，对于OMP算法，字典中的列数为N=H，为时延和多普勒网格数之积。因此，一次迭代的乘积运算为P=NKl。对于信道1，队=250，1=200，因而~=5\104;对于信道2，Ντ=250，Na=100，因而N=2·5X104。[0102]而对于IAFSA，信道1和信道2中，每一次迭代包含的子迭代过程中，人工鱼分别执行聚群和追尾行为，最差的情况下需要搜索21次，因而一次迭代的乘积运算为P=KLPkmax2I，即P=IX104。可见，本发明的计算复杂度优于OMP算法。[0103]尽管本发明就优选实施方式进行了示意和描述，但本领域的技术人员应当理解，只要不超出本发明的权利要求所限定的范围，可以对本发明进行各种变化和修改。

权利要求：1.一种适用于多普勒失真水声信道的估计方法，其特征在于，包括如下步骤：1在问题空间中初始化鱼群位置，计算相应的适应度值，并将群体中最优适应度值和对应的位置记录在公告板上，进入子迭代过程；2每一条人工鱼在其视野范围内执行聚群和追尾行为或觅食行为，更新自身位置和适应度值并更新公告板；3当子迭代次数大于设定值的一半时，若公告板中最优适应度值大于设定阈值且不发生变化，则将一半的人工鱼位置设置为最优适应度值对应的位置；⑷循环执行子迭代过程并不断调整步长，直至达到最大子迭代次数；5从公告板中得到最优位置，作为一条路径的参数，得到相应的信号分量，用以更新残余信号，进入下一次迭代。2.如权利要求1所述的适用于多普勒失真水声信道的估计方法，其特征在于，步骤（1中，问题空间即为路径参数可能的取值空间，包括时延和多普勒扩展因子的取值范围，一般认为最大时延扩展为训练序列的时间长度，最大多普勒扩展为收发端最大相对运动速度与声波在海水中的速度的比值。3.如权利要求1所述的适用于多普勒失真水声信道的估计方法，其特征在于，步骤（1中，人工鱼P的适应度值的计算公式为：其中r⑴为接收信号，s⑴为训练序列，XA人工鱼p的位置，〃为以Xp为时延-多普勒参数得到的训练序列。4.如权利要求1所述的适用于多普勒失真水声信道的估计方法，其特征在于，步骤2中，觅食行为是:人工鱼P在其视野范围内随机选取一个位置，若该位置的适应度值大于当前位置的适应度值，则向该位置移动一步;否则继续尝试，若尝试次数大于设定的最大值仍未成功，则随机移动一步。5.如权利要求1所述的适用于多普勒失真水声信道的估计方法，其特征在于，步骤2中，聚群行为是:人工鱼P在其视野范围内有Q个同伴，若Q〇，计算Q个同伴的中心位置X。和相应的适应度值y。，若ycQHyp，其中λ为拥挤度因子，则P向X。移动一步;若ycQ彡Ayp或Q=O，则执行觅食行为。6.如权利要求1所述的适用于多普勒失真水声信道的估计方法，其特征在于，步骤2中，追尾行为是:人工鱼P在其视野范围内内有Q个同伴，若Q〇，找到具有最优适应度值的同伴Xq，若其适应度值yq满足yqQAyp，则P向Xt^动一步，若yqQ彡Ayp或Q=〇,则执行觅食行为。7.如权利要求1所述的适用于多普勒失真水声信道的估计方法，其特征在于，步骤4中，第k次子迭代步长的调整方法是：其中，△为初始步长，k为第k次子迭代，kmax为子迭代最大次数。8.如权利要求1所述的适用于多普勒失真水声信道的估计方法，其特征在于，步骤5中，更新剩余信号的方法是：其中，Si和為分别为估计出的第1条路径的时延-多普勒信号和路径幅值。

百度查询： 东南大学 一种适用于多普勒失真水声信道的估计方法

免责声明
1、本报告根据公开、合法渠道获得相关数据和信息，力求客观、公正，但并不保证数据的最终完整性和准确性。
2、报告中的分析和结论仅反映本公司于发布本报告当日的职业理解，仅供参考使用，不能作为本公司承担任何法律责任的依据或者凭证。

阅读全文 双屏查看 官方信息 专利公告 收藏专利 下载PDF 下载WORD