一种最小均方误差稀疏水声信道估计方法技术

本发明专利技术涉及一种最小均方误差稀疏水声信道估计方法，首先，针对稀疏水声信道进行建模，以期获得稀疏水声信道冲激响应函数的表达框架，在此基础上。采用周期训练模式，设计一种基于最小均方误差变步长算法迭代寻优得到稀疏水声信道冲激响应函数，即估计得到时域水声信道信息。有益效果体现在：本发明专利技术基于滑动窗函数控制步长变化，同时结合稀疏特性进行似零伪范数约束，有效地提高了算法的收敛速度，也能使本发明专利技术有更低的稳态误差。

【技术实现步骤摘要】
一种最小均方误差稀疏水声信道估计方法
本专利技术属于水声通信和水声信号处理等领域，涉及一种最小均方误差稀疏水声信道估计方法，冲激干扰下接收信号与训练序列对稀疏水声信道的估计方法，该估计方法将提高在冲激干扰下对稀疏水声信道的估计性能。适用于水声信道估计、水声通信等。
技术介绍
稀疏水声信道估计和水声通信等问题都可归结为冲激响应函数的估计优化问题，基于训练序列和接收到的信号用以对水声信道表达进行估计。目前，对稀疏水声信道的估计方法包括有限冲激响应框架和时域的逐块估计框架。有限冲激响应框架的算法细节具体参见《NewsparseadaptivealgorithmsbasedonthenaturalgradientandtheL0-norm》，该文2013年发表于《IEEEJournalofOceanicEngineering》第38期，起始页码为323。时域的逐块估计框架详见《Estimationofrapidlytime-varyingsparsechannels》该文2007年发表于《IEEEJournalofOceanicEngineering》第32期，起始页码为927。由于水声信道的扩展和时变特性，导致对水声信道的冲激响应函数估计极为困难，此外，接收信号可能受到冲激信号干扰，传统的自适应估计方法在大冲激干扰下将失去对信道的估计作用。因此，考虑到水声信道的扩展和时变特性可由水声信道的时域进行刻画。本专利技术将建立在该模型基础上，对水声信道进行估计。然而自适应估计方法的固定步长无法同时保证算法的收敛速度和收敛精度，因此，设计一种新的步长更新策略可以既保证收敛速度，也同时有望提高估计精度。考虑到实际中的水声信道冲激响应函数的时变性和稀疏性，本专利技术将采用移动窗函数控制步长的变化以及采用似零伪范数进行稀疏约束，以期提高时变稀疏水声信道的估计跟踪能力。
技术实现思路
要解决的技术问题为了避免现有技术的不足之处，本专利技术提出一种最小均方误差稀疏水声信道估计方法，利用基于似零伪范数的最小均方误差算法的迭代框架，采用移动窗函数进行变步长运算，从而实现迭代的估计误差最小化，最终精确地恢复出稀疏水声信道的冲激响应函数。技术方案一种最小均方误差稀疏水声信道估计方法，其特征在于步骤如下：步骤1：设参数L为水声信道冲激响应长度，设置控制步长的窗函数更新参数η＝0.99，控制步长的窗函数加权因子C＝1.483[1+5/(L-1)]，初始化水声信道冲激响应的梯度因子为初始化水声信道冲激响应函数为h0＝0；步骤2：设置似零伪范数加权因子γ＝10-4，似零伪范数近似控制项β＝20；步骤3：给定输入训练信号x和输出信号y，重复以下第i次迭代：计算估计误差：其中yi为第i时刻的接收信号离散值，为第i时刻的用于训练的发射信号的转置，hi-1为第i-1时刻的信道冲激响应函数；计算梯度因子绝对值：其中||xi||表示对输入信号求欧式范数；将过去L个梯度因子绝对值构成一个窗：Bi＝[bi,bi-1,...,bi-L+1]；梯度因子更新：其中η表示更新参数，C表示加权因子，min(Bi)表示从Bi中选取最小值；生成步长更新表达式：计算稀疏约束梯度向量：fi＝βmax(0,1-β|hi-1|)sign(hi-1)；信道估计迭代：hi＝hi-1+μibixi-μiγfi。所述迭代次数等于处理数据的长度和信道的长度之差。有益效果本专利技术提出的一种最小均方误差稀疏水声信道估计方法，首先，针对稀疏水声信道进行建模，以期获得稀疏水声信道冲激响应函数的表达框架，在此基础上。采用周期训练模式，设计一种基于最小均方误差变步长算法迭代寻优得到稀疏水声信道冲激响应函数，即估计得到时域水声信道信息。有益效果体现在：本专利技术基于滑动窗函数控制步长变化，同时结合稀疏特性进行似零伪范数约束，有效地提高了算法的收敛速度，也能使本专利技术有更低的稳态误差。附图说明图1是本专利技术方法(Least-mean-squaresfractionwithL0norm,LMSF2-L0)与传统的正则化最小均方误差1类算法(Normalizedleast-mean-squares,NLMS1)及2类算法(NLMS2)，最小均方误差分式1类算法(Least-mean-squaresfraction,LMSF1，主要是采用传统步长更新方式)，最小均方误差分式2类算法(LMSF2，采用滑动窗的步长更新方式)，基于最小均方误差分式2类算法和L1范数约束的改进算法(LMSF2-L1)等方法的对比结果图，稀疏信道的稀疏度设置分别为K＝10，K＝20，K＝30，K＝40。可以看出当稀疏度K＝10，本专利技术方法比其他方法相比，具有最低的稳态误差，且收敛速度保持最快，当然，这一优势的差距随着稀疏度增加而趋于减小，这是因为本专利技术是基于稀疏特性的水声信道而设计的。图2是针对传统的NLMS1，NLMS2，LMSF1方法设置不同的形状控制参数，以确保NLMS1，NLMS2方法在收敛速度与稳态误差之间做出平衡，并与LMSF2、LMSF2-L1以及本专利技术方法(LMSF2-L0)对稀疏水声信道的估计结果对比图，从中可以看出，尽管传统算法在收敛速度与稳态误差之间各有侧重，与本专利技术方法相比仍显得不足。图3是测试NLMS1，NLMS2，LMSF1，LMSF2，LMSF2-L1方法与本专利技术方法(LMSF2-L0)在不同高斯白噪声干扰下的收敛速率的学习效果对比图，同样验证了本专利技术方法的优越性。图4是测试NLMS1，NLMS2，LMSF1，LMSF2，LMSF2-L1方法与本专利技术方法(LMSF2-L0)在迭代中途对稀疏信道进行跳变设置的学习效果对比图，可以看出，尽管在迭代中途，由于稀疏信道彻底变换，所有算法的学习曲线都产生了跳变，然而本专利技术方法在保证收敛最快的同时，显示了最低的稳态误差。具体实施方式现结合实施例、附图对本专利技术作进一步描述：为解决稀疏水声信道的估计问题，提出的具体实现步骤为：(1)设参数L为水声信道冲激响应长度，设置控制步长的窗函数更新参数η＝0.99，控制步长的窗函数加权因子C＝1.483[1+5/(L-1)]，初始化水声信道冲激响应的梯度因子为初始化水声信道冲激响应函数为h0＝0。(2)设置似零伪范数加权因子γ＝10-4，似零伪范数近似控制项β＝20。(3)给定输入训练信号x和输出信号y，重复以下第i次迭代：计算估计误差：其中yi为第i时刻的接收信号离散值，为第i时刻的用于训练的发射信号的转置，hi-1为第i-1时刻的信道冲激响应函数；梯度因子绝对值：其中||xi||表示对输入信号求欧式范数；将过去L个梯度因子绝对值构成一个窗：Bi＝[bi,bi-1,...,bi-L+1]；梯度因子更新：其中η表示更新参数，C表示加权因子，min(Bi)表示从Bi中选取最小值；生成步长更新本文档来自技高网...

【技术保护点】
1.一种最小均方误差稀疏水声信道估计方法，其特征在于步骤如下：/n步骤1：设参数L为水声信道冲激响应长度，设置控制步长的窗函数更新参数η＝0.99，控制步长的窗函数加权因子C＝1.483[1+5/(L-1)]，初始化水声信道冲激响应的梯度因子为

【技术特征摘要】
1.一种最小均方误差稀疏水声信道估计方法，其特征在于步骤如下：
步骤1：设参数L为水声信道冲激响应长度，设置控制步长的窗函数更新参数η＝0.99，控制步长的窗函数加权因子C＝1.483[1+5/(L-1)]，初始化水声信道冲激响应的梯度因子为初始化水声信道冲激响应函数为h0＝0；
步骤2：设置似零伪范数加权因子γ＝10-4，似零伪范数近似控制项β＝20；
步骤3：给定输入训练信号x和输出信号y，重复以下第i次迭代：
计算估计误差：
其中yi为第i时刻的接收信号离散值，为第i时刻的用于训练的发射信号的转置，hi-1为第i-1时刻的信道冲激响应函数；

【专利技术属性】
技术研发人员：伍飞云，杨坤德，
申请(专利权)人：西北工业大学，
类型：发明
国别省市：陕西;61