CN111431626A - 一种基于训练的自适应水声信道跟踪方法 - Google Patents
一种基于训练的自适应水声信道跟踪方法 Download PDFInfo
- Publication number
- CN111431626A CN111431626A CN202010138589.7A CN202010138589A CN111431626A CN 111431626 A CN111431626 A CN 111431626A CN 202010138589 A CN202010138589 A CN 202010138589A CN 111431626 A CN111431626 A CN 111431626A
- Authority
- CN
- China
- Prior art keywords
- channel
- tracking
- kalman
- training
- adaptive
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 86
- 230000003044 adaptive effect Effects 0.000 title claims abstract description 55
- 238000012549 training Methods 0.000 title claims abstract description 54
- 238000012937 correction Methods 0.000 claims abstract description 16
- 230000002457 bidirectional effect Effects 0.000 claims abstract description 11
- 238000005259 measurement Methods 0.000 claims abstract description 7
- 239000011159 matrix material Substances 0.000 claims description 63
- 230000007704 transition Effects 0.000 claims description 43
- 238000001914 filtration Methods 0.000 claims description 15
- 230000008569 process Effects 0.000 claims description 14
- 230000009467 reduction Effects 0.000 claims description 10
- 230000000694 effects Effects 0.000 claims description 8
- 238000004364 calculation method Methods 0.000 claims description 6
- 230000008859 change Effects 0.000 claims description 6
- 238000000354 decomposition reaction Methods 0.000 claims description 6
- 230000006835 compression Effects 0.000 claims description 3
- 238000007906 compression Methods 0.000 claims description 3
- 230000009286 beneficial effect Effects 0.000 abstract description 2
- 230000006872 improvement Effects 0.000 description 5
- 238000013524 data verification Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B13/00—Transmission systems characterised by the medium used for transmission, not provided for in groups H04B3/00 - H04B11/00
- H04B13/02—Transmission systems in which the medium consists of the earth or a large mass of water thereon, e.g. earth telegraphy
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B17/00—Monitoring; Testing
- H04B17/30—Monitoring; Testing of propagation channels
- H04B17/391—Modelling the propagation channel
- H04B17/3913—Predictive models, e.g. based on neural network models
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L25/00—Baseband systems
- H04L25/02—Details ; arrangements for supplying electrical power along data transmission lines
- H04L25/0202—Channel estimation
- H04L25/024—Channel estimation channel estimation algorithms
- H04L25/0256—Channel estimation using minimum mean square error criteria
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Power Engineering (AREA)
- Artificial Intelligence (AREA)
- Evolutionary Computation (AREA)
- Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
- Filters That Use Time-Delay Elements (AREA)
Abstract
本发明提供了一种基于训练的自适应水声信道跟踪方法,通过基于训练的自适应ASRMAE信道跟踪方法得到跟踪的信道,使用自适应卡尔曼跟踪器自适应的校正预测误差来改善跟踪性能,并通过训练数据获得自适应卡尔曼跟踪器的参数,通过基于双向卡尔曼ASRMAE信道跟踪方法,通过使用前向和后向时刻的测量值来改善信道主成分的跟踪。本发明的有益效果是:跟踪性能大幅提升。
Description
技术领域
本发明涉及水声信道跟踪方法,尤其涉及一种基于训练的自适应水声信道跟踪方法。
背景技术
现有研究表明,许多水声信道中的多路径间是具有一定相关性的。利用多路径的互相关性,将高维信道分解成具有低维的信道主成分和相应的信道子空间,提出了一类有效的水声信道跟踪方法,在提升跟踪性能的同时降低计算复杂度。为了跟踪信道,将信道主成分建模为自回归 (Autoregressive,AR)过程,在此模型的基础上使用卡尔曼滤波器对其进行跟踪。信道子空间也可以通过递归算法进行跟踪。但是,这类方法,如基于子空间跟踪的降维的模基估计(Adaptive Subspace-tracking with Reduced-rank Model-basedAmplitude Estimation,ASRMAE)方法留下了许多不确定的参数,不同的参数导致方法具有不同的跟踪性能;对于复杂多变的水声信道,先验模型的失配是不可避免的。这些都会导致跟踪性能的下降。
发明内容
为了解决现有技术中的问题,本发明提供了一种基于训练的自适应水声信道跟踪方法。
本发明提供了一种基于训练的自适应水声信道跟踪方法,通过基于训练的自适应ASRMAE(基于子空间跟踪的降维的模基估计,Adaptive Subspace-tracking withReduced-rank Model-based Amplitude Estimation)信道跟踪方法得到跟踪的信道,使用自适应卡尔曼跟踪器自适应的校正预测误差来改善跟踪性能,并通过训练数据获得自适应卡尔曼跟踪器的参数,通过基于双向卡尔曼ASRMAE信道跟踪方法,通过使用前向和后向时刻的测量值来改善信道主成分的跟踪。
作为本发明的进一步改进,假设有M条路径的水声信道 h(n)=[h1(n),h2(n),…,hM(n)]T,根据信道路径间的相关性,将高维度信道分解为低维度的信道主成分zr(n)和相应的信道子空间Qr(n),其中r为降秩后的维度,并有r<<M,通过分别跟踪信道主成分zr(n)和相应的信道子空间Qr(n)实现对信道的跟踪:
h(n)=Qr(n)zr(n)
对于信道主成分,首先对其建立自回归模型,简称为AR模型,然后在此基础上采用自适应卡尔曼跟踪器实现对其跟踪。
作为本发明的进一步改进,基于训练的自适应的ASRMAE信道跟踪方法的步骤如下:首先通过LMS算法获得信道的粗估计通过训练过程,信道子空间Qr通过信道自相关矩阵的 EVD(特征值分解,eigenvalue decompositio)分解获得,然后通过粗估计的信道主成分获得AR模型阶数ptr、训练序列的模型参数状态转移矩阵Φz和噪声协方差矩阵RN,这时,AR模型得以建立,通过训练确定信道降维处理后的最优维度rtr、可变遗忘因子α(n)的最优步长μtr和最小值在每一时刻,通过PASTd方法跟踪信道子空间 Qr(n),并通过自适应卡尔曼跟踪器获得信道主成分状态向量z*(n),再从中提取出信道主成分zr(n),最后得到跟踪的信道
作为本发明的进一步改进,基于双向卡尔曼的ASRMAE信道跟踪方法步骤如下:在对信道主成分进行跟踪前,和基于训练的自适应ASRMAE 的信道跟踪方法一样,得到AR模型参数和每一时刻的信道子空间Qr(n),然后使用双向卡尔曼滤波算法分别获得前向估计值z*f(n|n)和后向估计值通过两者得到信道主成分的最优估计从中提取出信道主成分zr(n),最后就可以得到信道的跟踪结果
作为本发明的进一步改进,基于训练的自适应的ASRMAE信道跟踪方法包括:
A1、信道主成分的状态转移模型;
(1)状态转移模型的建立
对于一个p阶马尔科夫离散随机过程,AR模型表示为如下形式:
式中Φ(n,l)——代表状态转移矩阵;
η(n)——代表高斯过程噪声,并有η(n)=[η1(n),η2(n),…ηr(n)]T;
为了能够得到卡尔曼的状态转移方程,将(1)改写为如下形式:
假设传输的信号向量为d(n)=[d(n),d(n-1),…,d(n-M)]T,接收信号为 r(n)=d(n)Hh(n)+v(n),v(n)假设为服从均值为零,方差为σ2=E[v(n)v*(n)]的高斯白噪声;
定义一个扩张数据向量D(n)=[dH(n)Qr(n)01×r(p-1)],信道主成分状态向量 z*(n)=[zr T(n),…,zr T(n-p+1)]T,因此获得状态转移方程和预测方程:
z*(n)=Φz(n)z*(n-1)+N(n) (3)
r(n)=D(n)z*(n)+v(n) (4)
式中Φz(n)——代表(2)中的状态转移矩阵;
N(n)——代表过程噪声,N(n)=[ηT(n),01×r(p-1)]T,并有协方差矩阵
RN=E[N(n)NH(n)];
通过训练序列估计出状态转移矩阵Φz(n)和噪声协方差矩阵RN,并将两个参数带入到卡尔曼滤波算法中就可以实现对信道的实时跟踪;
(2)AR模型的定阶及参数估计;
为了寻找最佳的AR模型阶数,通过训练序列使用AIC准则(赤池信息量准则,Akaike information criterion),计算得到不同阶数的AIC值,来进行AR 模型的定阶,p阶AR模型的AIC定义为:
AIC(p)=2p+Ntrln(S/Ntr) (5)
式中S——代表残差平方和;
Ntr——代表训练数据序列的个数;
状态转移矩阵假设为时不变,即Φz(n)=Φz,因此,采用最小二乘法对训练序列的模型参数状态转移矩阵Φz和噪声协方差矩阵RN进行估计,并用于后续时刻的信道跟踪过程中;
A2、基于训练的自适应ASRMAE信道跟踪方法;
(1)自适应卡尔曼跟踪器;
采用自适应卡尔曼跟踪器对信道主成分进行跟踪,通过可变遗忘因子α对卡尔曼增益和误差协方差进行自适应更新,自适应的校正由于模型失配带来的误差;
使用可变遗忘因子的卡尔曼滤波增益K(n)和误差协方差矩阵P(n+1|n) 的计算公式如下所示:
定义预测误差方差其中跟踪的估计误差 e(n|n-1)e(n|n-1)=r(n)-hH(n|n-1)d(n),h(n|n-1)为由h(n-1)预测出的n时刻信道状态,当有明显的增大时,则自适应的增加卡尔曼增益K(n)和误差协方差矩阵P(n+1|n),突出新的测量值对预测结果的校正作用,α(n)的值越小就表示对预测结果的作用越大,为了实现α(n)对误差的自适应校正作用,令:
α(n)=α(n-1)+μ·β (8)
式中μ——代表步长;
β用来决定可变遗忘因子的增大或减小,其值由前后两个时刻的预测误差方差的比值来决定,在跟踪过程中,若预测误差方差有较大的增加趋势,并且前后两次的比值预测误差方差达到门限值Gh时,取β=-1,否则取β=1,即
另外,为了防止滤波发散,当α(n)超出限定范围时,需要对α(n)的值进行修正,当α(n)大于1时就修正为1,当α(n)小于最小值αmin时则修正为αmin,即
自适应卡尔曼滤波器流程包括:首先建立系统模型,输入发送信号和接收信号,由h(n-1)预测出n时刻的信道状态h(n|n-1),并计算信号预测误差e(n|n-1)以及误差方差然后将与上一时刻的误差方差进行比较,判断增长是否达到门限值,若达到门限值,则根据式(8)和(10)更新遗忘因子α(n)并进行修正,然后利用新的α(n)进行校正更新得到h(n|n)和 P(n+1|n),如此重复直到达到跟踪次数;
(2)自适应卡尔曼跟踪器参数的设置;
对于参数μ和αmin,通过使用网格搜索法来寻优,对于参数μ和αmin,用均方误差来评价算法的跟踪性能,均方误差简称为MSE:
同样的,对于降秩后的维度r,也使用MSE来描述自适应卡尔曼跟踪器的性能,通过训练序列获取性能最佳的参数设置值。
作为本发明的进一步改进,基于双向卡尔曼的ASRMAE信道跟踪方法包括:
B1、双向卡尔曼跟踪器;
前向卡尔曼滤波从时刻1向时刻N顺序进行跟踪,即通过时刻n的状态值预测出时刻n+1的状态值,对于后向卡尔曼,则是从时刻N向时刻1顺序进行跟踪,即通过n时刻的状态值预测出n-1的状态值,已知前向卡尔曼的状态转移方程:
z*(n)=ΦZz*(n-1)+N(n) (12)
由于ΦZ是一个非零的对角矩阵,所以在等式两边同乘ΦZ -1有:
ΦZ -1z*(n)=z*(n-1)+ΦZ -1N(n) (13)
将上式进行整理,得到后向状态转移方程:
z*(n-1)=Φbz*(n)+Nb(n) (14)
其中Φb=ΦZ -1,Nb(n)=-ΦZ -1N(n)并且其方差为ΦbRfNΦb H;
B2、信道主成分的最优估计;
假设有两个如下的线性系统:
y1=A1w+q1 (15a)
y2=A2w+q2 (15b)
式中y1——代表系统一的观测向量;
y2——代表系统二的观测向量;
A1——代表系统一的状态转移矩阵;
A2——代表系统二的状态转移矩阵;
q1——代表系统一的高斯观测噪声;
q2——代表系统二的高斯观测噪声,与q1互不相关;
线性最小均方误差简称为LMMSE,每个系统的w的LMMSE估计由下式给出:
Rw——代表w的协方差矩阵;
Rq1——代表q1的协方差矩阵;
同理有
Rq2——代表q2的协方差矩阵。
计算同时给定y1和y2时w的估计值,有级联线性系统如下:
y=Aw+q (20)
通过求解式(20)得到:
在对w没有其自相关的先验信息的时候,去掉Rw,得到信道主成分z*(n) 的最优估计:
M(n|n)=(Pf(n|n)-1+Pb(n|n)-1)-1 (24)
其中Pf(n|n)和Pb(n|n)分别由前向卡尔曼滤波算法和后向卡尔曼滤波算法计算得到。
本发明的有益效果是:提出了基于训练的自适应ASRMAE信道跟踪方法,通过使用自适应卡尔曼滤波器自适应的校正预测误差,改善跟踪性能;并通过训练数据获得跟踪器的参数,还提出了基于双向卡尔曼ASRMAE信道跟踪方法,通过使用前向和后向时刻的测量值来改善信道主成分的跟踪,从另一个角度弥补模型失配带来的跟踪性能损失,通过海洋数据验证,提出的上述两种信道跟踪方法较现有方法有较低的信号预测误差,跟踪性能大幅提升。
附图说明
图1是本发明一种基于训练的自适应水声信道跟踪方法的自适应卡尔曼滤波器流程图。
具体实施方式
下面结合附图说明及具体实施方式对本发明作进一步说明。
一种基于训练的自适应水声信道跟踪方法,过使用自适应卡尔曼滤波器自适应的校正预测误差,改善跟踪性能。并通过训练数据获得跟踪器的参数。我们还提出了基于双向卡尔曼ASRMAE(基于子空间跟踪的降维的模基估计,Adaptive Subspace-tracking withReduced-rank Model-based Amplitude Estimation)信道跟踪方法,通过使用前向和后向时刻的测量值来改善信道主成分的跟踪,从另一个角度弥补模型失配带来的跟踪性能损失。通过海洋数据验证,我们提出的上述两种信道跟踪方法较现有方法有较低的信号预测误差,跟踪性能大幅提升。
为了改善现有水声信道跟踪方法由于不定参数和模型失配导致跟踪性能下降的问题,本发明提供了一种基于训练的自适应水声信道跟踪方法,包括基于训练的自适应ASRMAE信道跟踪方法和基于双向卡尔曼的 ASRMAE信道跟踪方法。
假设有M条路径的水声信道h(n)=[h1(n),h2(n),…,hM(n)]T,根据信道路径间的相关性,将高维度信道分解为低维度的信道主成分zr(n)和相应的信道子空间Qr(n),其中r为降秩后的维度,并有r<<M,通过分别跟踪两者实现对信道的跟踪:
h(n)=Qr(n)zr(n)
首先使用LMS方法获得信道的粗估计然后采用基于压缩的投影近似子空间跟踪(Projection Approximation Subspace Tracking based on the deflationtechnique,PASTd)方法对信道子空间进行跟踪。
对于信道主成分,首先对其建立AR模型,然后在此基础上采用卡尔曼滤波器实现对其跟踪。
1.信道主成分的状态转移模型
(1)状态转移模型的建立
对于一个p阶马尔科夫离散随机过程,AR模型可以表示为如下形式:
式中Φ(n,l)——代表状态转移矩阵;
η(n)——代表高斯过程噪声,并有η(n)=[η1(n),η2(n),…ηr(n)]T。
为了能够得到卡尔曼的状态转移方程,将(1)改写为如下形式:
假设传输的信号向量为d(n)=[d(n),d(n-1),…,d(n-M)]T,接收信号为 r(n)=d(n)Hh(n)+v(n),v(n)假设为服从均值为零,方差为σ2=E[v(n)v*(n)]的高斯白噪声。
定义一个扩张数据向量D(n)=[dH(n)Qr(n) 01×r(p-1)],信道主成分状态向量 z*(n)=[zr T(n),…,zr T(n-p+1)]T,因此可以获得状态转移方程和预测方程:
z*(n)=Φz(n)z*(n-1)+N(n) (3)
r(n)=D(n)z*(n)+v(n) (4)
式中Φz(n)——代表(2)中的状态转移矩阵;
N(n)——代表过程噪声,N(n)=[ηT(n),01×r(p-1)]T,并有协方差矩阵
RN=E[N(n)NH(n)]。
通过训练序列可以估计出状态转移矩阵Φz(n)和噪声协方差矩阵RN,并将两个参数带入到卡尔曼滤波算法中就可以实现对信道的实时跟踪。
(2)AR模型的定阶及参数估计
为了寻找最佳的AR模型阶数,我们通过训练序列使用AIC准则(赤池信息量准则,Akaike information criterion),计算得到不同阶数的AIC值,来进行AR模型的定阶。p阶AR模型的AIC定义为:
AIC(p)=2p+Ntrln(S/Ntr) (5)
式中S——代表残差平方和;
Ntr——代表训练数据序列的个数。
在一般情况下,状态转移矩阵可假设为时不变,即Φz(n)=Φz,因此我们采用最小二乘法对训练序列的模型参数状态转移矩阵Φz和噪声协方差矩阵 RN进行估计,并用于后续时刻的信道跟踪过程中。
2.基于训练的自适应ASRMAE信道跟踪方法
(3)自适应卡尔曼跟踪器
卡尔曼滤波算法对信道跟踪具有很好的性能,所以在很多基于模型的信道跟踪方法中都采用卡尔曼滤波算法,但是模型的精确度对卡尔曼滤波算法的跟踪性能有很大的影响。我们在对信道主成分进行建模时,不可能完全精确地去拟合其时变性,因此卡尔曼滤波算法将有可能在滤波过程中发散,最终导致跟踪性能降低或失败。对此我们采用自适应卡尔曼跟踪器对信道主成分进行跟踪,通过可变遗忘因子α对卡尔曼增益和误差协方差进行自适应更新,自适应的校正由于模型失配带来的误差。
使用可变遗忘因子的卡尔曼滤波增益K(n)和误差协方差矩阵P(n+1|n) 的计算公式如下所示:
定义预测误差方差其中跟踪的估计误差 e(n|n-1)=r(n)-hH(n|n-1)d(n),h(n|n-1)为由h(n-1)预测出的n时刻信道状态,当有明显的增大时,则自适应的增加卡尔曼增益K(n)和误差协方差矩阵P(n+1|n),突出新的测量值对预测结果的校正作用,α(n)的值越小就表示对预测结果的作用越大。为了实现α(n)对误差的自适应校正作用,令:
α(n)=α(n-1)+μ·β (8)
式中μ——代表步长。
β用来决定可变遗忘因子的增大或减小,其值由前后两个时刻的预测误差方差的比值来决定。在跟踪过程中,若预测误差方差有较大的增加趋势,并且前后两次的比值预测误差方差达到门限值Gh时,取β=-1,否则取β=1。
即
另外,为了防止滤波发散,当α(n)超出限定范围时,需要对α(n)的值进行修正,当α(n)大于1时就修正为1,当α(n)小于最小值αmin时则修正为αmin,即
综上所述,若没有太大的变化,则α(n)始终趋于1,系统达到稳定状态;若变化较大,则通过减小α(n)来增强新的观测数据对预测结果的校正作用,这样α(n)便可以根据的的变化自适应变化。图1为自适应卡尔曼滤波器流程图。首先建立系统模型,输入发送信号和接收信号,由h(n-1)预测出n时刻的信道状态h(n|n-1),并计算信号预测误差e(n|n-1)以及误差方差然后将与上一时刻的误差方差进行比较,判断增长是否达到门限值,若达到门限值,则根据式(8)和(10)更新遗忘因子α(n) 并进行修正,然后利用新的α(n)进行校正更新得到h(n|n)和P(n+1|n),如此重复直到达到跟踪次数。
(4)跟踪器参数的设置
对于参数μ和αmin我们通过使用网格搜索法来寻优,网格搜索法是一种最基本的参数优化算法。它实质上就是将待搜索参数的各个可能的取值进行排列组合,并生成“网格”,坐标系中每一个点都代表一组参数,进而可以通过将这个给定区间内的每个点带入到跟踪器中验证其性能,使得整个跟踪器的性能最佳的那个点我们称它为最优参数。对于参数μ和αmin我们用均方误差(Mean Square Error,MSE)来评价算法的跟踪性能:
同样的,对于降秩后的维度r,也使用MSE来描述跟踪器的性能。通过训练序列即可获取性能最佳的参数设置值。
基于训练的自适应的ASRMAE信道跟踪方法的步骤如下:首先通过LMS算法获得信道的粗估计通过上述的训练过程,信道子空间Qr可以通过信道自相关矩阵的EVD(特征值分解, eigenvalue decompositio)分解获得。然后通过粗估计的可以获得AR模型阶数ptr、状态转移矩阵Φz和噪声协方差矩阵RN。这时AR 模型得以建立,通过训练确定信道降维处理后的最优维度rtr、可变遗忘因子α(n)的最优步长μtr和最小值在每一时刻,通过PASTd方法跟踪信道子空间Qr(n),并通过自适应卡尔曼跟踪器获得z*(n),再从中提取出zr(n),最后便可得到跟踪的信道
3.基于双向卡尔曼的ASRMAE信道跟踪方法
(1)双向卡尔曼跟踪器
前向卡尔曼滤波从时刻1向时刻N顺序进行跟踪,即通过时刻n的状态值预测出时刻n+1的状态值。对于后向卡尔曼,则是从时刻N向时刻1顺序进行跟踪,即通过n时刻的状态值预测出n-1的状态值。那么这里需要解决的是如何建立反向的状态转移方程,已知前向卡尔曼的状态转移方程:
z*(n)=ΦZz*(n-1)+N(n) (12)
由于ΦZ是一个非零的对角矩阵,所以可以在等式两边同乘ΦZ -1有:
ΦZ -1z*(n)=z*(n-1)+ΦZ -1N(n) (13)
将上式进行整理,可以得到后向状态转移方程:
z*(n-1)=Φbz*(n)+Nb(n) (14)
其中Φb=ΦZ -1,Nb(n)=-ΦZ -1N(n)并且其方差为ΦbRfNΦb H。
(2)信道主成分的最优估计
假设我们有两个如下的线性系统:
y1=A1w+q1 (15a)
y2=A2w+q2 (15b)
式中y1——代表系统一的观测向量;
y2——代表系统二的观测向量;
A1——代表系统一的状态转移矩阵;
A2——代表系统二的状态转移矩阵;
q1——代表系统一的高斯观测噪声;
q2——代表系统二的高斯观测噪声,与q1互不相关。
每个系统的w的线性最小均方误差(Linear Minimum Mean Square Error,LMMSE)估计由下式给出:
Rw——代表w的协方差矩阵;
Rq1——代表q1的协方差矩阵。
同理有
Rq2——代表q2的协方差矩阵。
现在我们计算同时给定y1和y2时w的估计值,有级联线性系统如下:
y=Aw+q (20)
通过求解式(20)可以得到:
在对w没有其自相关的先验信息的时候,可以去掉Rw。现在我们可以应用结论来得到n时刻信道主成分的最优估计值。本发明中定义的状态模型是线性的,即卡尔曼滤波算法是通过LMMSE实现的,因此上述的结果可以用于双向卡尔曼滤波算法中,即可以得到信道主成分z*(n)的最优估计:
M(n|n)=(Pf(n|n)-1+Pb(n|n)-1)-1 (24)
基于双向卡尔曼的ASRMAE信道跟踪方法步骤如下:在对信道主成分进行跟踪前,和基于训练的自适应ASRMAE的信道跟踪方法一样需要得到 AR模型参数和每一时刻的信道子空间Qr(n)。然后使用双向卡尔曼滤波算法分别获得前向估计值z*f(n|n)和后向估计值通过两者得到信道主成分的最优估计从中提取出zr(n),最后就可以得到信道的跟踪结果
以上内容是结合具体的优选实施方式对本发明所作的进一步详细说明,不能认定本发明的具体实施只局限于这些说明。对于本发明所属技术领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干简单推演或替换,都应当视为属于本发明的保护范围。
Claims (6)
1.一种基于训练的自适应水声信道跟踪方法,其特征在于:通过基于训练的自适应的基于子空间跟踪的降维的模基估计(ASRMAE)信道跟踪方法得到跟踪的信道,使用自适应卡尔曼跟踪器自适应的校正预测误差来改善跟踪性能,并通过训练数据获得自适应卡尔曼跟踪器的参数,通过基于双向卡尔曼ASRMAE信道跟踪方法,通过使用前向和后向时刻的测量值来改善信道主成分的跟踪。
2.根据权利要求1所述的基于训练的自适应水声信道跟踪方法,其特征在于:假设有M条路径的水声信道h(n)=[h1(n),h2(n),…,hM(n)]T,根据信道路径间的相关性,将高维度信道分解为低维度的信道主成分zr(n)和相应的信道子空间Qr(n),其中r为降秩后的维度,并有r<<M,通过分别跟踪信道主成分zr(n)和相应的信道子空间Qr(n)实现对信道的跟踪:
h(n)=Qr(n)zr(n)
对于信道主成分,首先对其建立自回归模型,简称为AR模型,然后在此基础上采用自适应卡尔曼跟踪器实现对其跟踪。
3.根据权利要求2所述的基于训练的自适应水声信道跟踪方法,其特征在于:基于训练的自适应的ASRMAE信道跟踪方法的步骤如下:首先通过LMS算法获得信道的粗估计通过训练过程,信道子空间Qr通过信道自相关矩阵的特征值分解(EVD)分解获得,然后通过粗估计的信道主成分获得AR模型阶数ptr、状态转移矩阵Φz和噪声协方差矩阵RN,这时,AR模型得以建立,通过训练确定信道降维处理后的最优维度rtr、可变遗忘因子α(n)的最优步长μtr和最小值在每一时刻,通过PASTd方法跟踪信道子空间Qr(n),并通过自适应卡尔曼跟踪器获得信道主成分状态向量z*(n),再从中提取出信道主成分zr(n),最后得到跟踪的信道
5.根据权利要求4所述的基于训练的自适应水声信道跟踪方法,其特征在于:基于训练的自适应的ASRMAE信道跟踪方法包括:
A1、信道主成分的状态转移模型;
(1)状态转移模型的建立
对于一个p阶马尔科夫离散随机过程,AR模型表示为如下形式:
式中Φ(n,l)——代表状态转移矩阵;
η(n)——代表高斯过程噪声,并有η(n)=[η1(n),η2(n),…ηr(n)]T;
为了能够得到卡尔曼的状态转移方程,将(1)改写为如下形式:
假设传输的信号向量为d(n)=[d(n),d(n-1),…,d(n-M)]T,接收信号为r(n)=d(n)Hh(n)+v(n),v(n)假设为服从均值为零,方差为σ2=E[v(n)v*(n)]的高斯白噪声;
定义一个扩张数据向量D(n)=[dH(n)Qr(n)01×r(p-1)],信道主成分状态向量z*(n)=[zr T(n),…,zr T(n-p+1)]T,因此获得状态转移方程和预测方程:
z*(n)=Φz(n)z*(n-1)+N(n) (3)
r(n)=D(n)z*(n)+v(n) (4)
式中Φz(n)——代表(2)中的状态转移矩阵;
N(n)——代表过程噪声,N(n)=[ηT(n),01×r(p-1)]T,并有协方差矩阵RN=E[N(n)NH(n)];
通过训练序列估计出状态转移矩阵Φz(n)和噪声协方差矩阵RN,并将两个参数带入到卡尔曼滤波算法中就可以实现对信道的实时跟踪;
(2)AR模型的定阶及参数估计;
为了寻找最佳的AR模型阶数,通过训练序列使用赤池信息量准则(AIC)准则,计算得到不同阶数的AIC值,来进行AR模型的定阶,p阶AR模型的AIC定义为:
AIC(p)=2p+Ntrln(S/Ntr) (5)
式中S——代表残差平方和;
Ntr——代表训练数据序列的个数;
状态转移矩阵假设为时不变,即Φz(n)=Φz,因此,采用最小二乘法对训练序列的模型参数状态转移矩阵Φz和噪声协方差矩阵RN进行估计,并用于后续时刻的信道跟踪过程中;
A2、基于训练的自适应的ASRMAE信道跟踪方法;
(1)自适应卡尔曼跟踪器;
采用自适应卡尔曼跟踪器对信道主成分进行跟踪,通过可变遗忘因子α对卡尔曼增益和误差协方差进行自适应更新,自适应的校正由于模型失配带来的误差;
使用可变遗忘因子的卡尔曼滤波增益K(n)和误差协方差矩阵P(n+1|n)的计算公式如下所示:
定义预测误差方差其中跟踪的估计误差e(n|n-1)=r(n)-hH(n|n-1)d(n),h(n|n-1)为由h(n-1)预测出的n时刻信道状态,当有明显的增大时,则自适应的增加卡尔曼增益K(n)和误差协方差矩阵P(n+1|n),突出新的测量值对预测结果的校正作用,α(n)的值越小就表示对预测结果的作用越大,为了实现α(n)对误差的自适应校正作用,令:
α(n)=α(n-1)+μ·β (8)
式中μ——代表步长;
β用来决定可变遗忘因子的增大或减小,其值由前后两个时刻的预测误差方差的比值来决定,在跟踪过程中,若预测误差方差有较大的增加趋势,并且前后两次的比值预测误差方差达到门限值Gh时,取β=-1,否则取β=1,即
另外,为了防止滤波发散,当α(n)超出限定范围时,需要对α(n)的值进行修正,当α(n)大于1时就修正为1,当α(n)小于最小值αmin时则修正为αmin,即
自适应卡尔曼滤波器流程包括:首先建立系统模型,输入发送信号和接收信号,由h(n-1)预测出n时刻的信道状态h(n|n-1),并计算信号预测误差e(n|n-1)以及误差方差然后将与上一时刻的误差方差进行比较,判断增长是否达到门限值,若达到门限值,则根据式(8)和(10)更新遗忘因子α(n)并进行修正,然后利用新的α(n)进行校正更新得到h(n|n)和P(n+1|n),如此重复直到达到跟踪次数;
(2)自适应卡尔曼跟踪器参数的设置;
对于参数μ和αmin,通过使用网格搜索法来寻优,对于参数μ和αmin,用均方误差来评价算法的跟踪性能,均方误差简称为MSE:
同样的,对于降秩后的维度r,也使用MSE来描述自适应卡尔曼跟踪器的性能,通过训练序列获取性能最佳的参数设置值。
6.根据权利要求1所述的基于训练的自适应水声信道跟踪方法,其特征在于:基于双向卡尔曼的ASRMAE信道跟踪方法包括:
B1、双向卡尔曼跟踪器;
前向卡尔曼滤波从时刻1向时刻N顺序进行跟踪,即通过时刻n的状态值预测出时刻n+1的状态值,对于后向卡尔曼,则是从时刻N向时刻1顺序进行跟踪,即通过n时刻的状态值预测出n-1的状态值,已知前向卡尔曼的状态转移方程:
z*(n)=ΦZz*(n-1)+N(n) (12)
由于ΦZ是一个非零的对角矩阵,所以在等式两边同乘ΦZ -1有:
ΦZ -1z*(n)=z*(n-1)+ΦZ -1N(n) (13)
将上式进行整理,得到后向状态转移方程:
z*(n-1)=Φbz*(n)+Nb(n) (14)
其中Φb=ΦZ -1,Nb(n)=-ΦZ -1N(n)并且其方差为ΦbRfNΦb H;
B2、信道主成分的最优估计;
假设有两个如下的线性系统:
y1=A1w+q1 (15a)
y2=A2w+q2 (15b)
式中y1——代表系统一的观测向量;
y2——代表系统二的观测向量;
A1——代表系统一的状态转移矩阵;
A2——代表系统二的状态转移矩阵;
q1——代表系统一的高斯观测噪声;
q2——代表系统二的高斯观测噪声,与q1互不相关;
线性最小均方误差简称为LMMSE,每个系统的w的LMMSE估计由下式给出:
Rw——代表w的协方差矩阵;
Rq1——代表q1的协方差矩阵;
同理有
Rq2——代表q2的协方差矩阵。
计算同时给定y1和y2时w的估计值,有级联线性系统如下:
y=Aw+q (20)
通过求解式(20)得到:
在对w没有其自相关的先验信息的时候,去掉Rw,得到信道主成分z*(n)的最优估计:
M(n|n)=(Pf(n|n)-1+Pb(n|n)-1)-1 (24)
其中Pf(n|n)和Pb(n|n)分别由前向卡尔曼滤波算法和后向卡尔曼滤波算法计算得到。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010138589.7A CN111431626A (zh) | 2020-03-03 | 2020-03-03 | 一种基于训练的自适应水声信道跟踪方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010138589.7A CN111431626A (zh) | 2020-03-03 | 2020-03-03 | 一种基于训练的自适应水声信道跟踪方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN111431626A true CN111431626A (zh) | 2020-07-17 |
Family
ID=71547645
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010138589.7A Pending CN111431626A (zh) | 2020-03-03 | 2020-03-03 | 一种基于训练的自适应水声信道跟踪方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111431626A (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115114578A (zh) * | 2021-03-19 | 2022-09-27 | Aptiv技术有限公司 | 确定对象传感器的噪声统计的方法 |
CN115460047A (zh) * | 2022-08-02 | 2022-12-09 | 北京睿信丰科技有限公司 | 复杂信道下快速跟踪频偏和距离的方法、装置及电子设备 |
CN117375731A (zh) * | 2023-10-27 | 2024-01-09 | 中国科学院声学研究所南海研究站 | 一种水声通信多途信道时变特性量化分析方法及系统 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102222911A (zh) * | 2011-04-19 | 2011-10-19 | 哈尔滨工业大学 | 基于ar模型和卡尔曼滤波的电力系统间谐波估计方法 |
CN102571666A (zh) * | 2011-08-12 | 2012-07-11 | 哈尔滨工程大学 | 基于mmse的水声ofdm判决迭代信道均衡方法 |
KR101504289B1 (ko) * | 2014-02-26 | 2015-03-30 | 연세대학교 산학협력단 | 수중 채널 응답 생성 장치 및 방법 |
CN105553898A (zh) * | 2015-12-18 | 2016-05-04 | 中国电子科技集团公司第三研究所 | 均衡器及反馈均衡方法 |
-
2020
- 2020-03-03 CN CN202010138589.7A patent/CN111431626A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102222911A (zh) * | 2011-04-19 | 2011-10-19 | 哈尔滨工业大学 | 基于ar模型和卡尔曼滤波的电力系统间谐波估计方法 |
CN102571666A (zh) * | 2011-08-12 | 2012-07-11 | 哈尔滨工程大学 | 基于mmse的水声ofdm判决迭代信道均衡方法 |
KR101504289B1 (ko) * | 2014-02-26 | 2015-03-30 | 연세대학교 산학협력단 | 수중 채널 응답 생성 장치 및 방법 |
CN105553898A (zh) * | 2015-12-18 | 2016-05-04 | 中国电子科技集团公司第三研究所 | 均衡器及反馈均衡方法 |
Non-Patent Citations (2)
Title |
---|
MURWAN BASHIR,DMITRI TRUHACHEV: "Kalman Forward-Backward Channel Tracking and Combining for OFDM in Underwater Acoustic Channels", 《2018 OCEANS-MTS/IEEE KOBE TECHNO-OCEANS(OTO)》 * |
XIAOYUAN WANG,WEI LI: "Training-based Adaptive Channel Tracking for Correlated Underwater Acoustic Channels", 《PROCEEDINGS OF THE ACM INTERNATIONAL WORKSHOP ON UNDERWATER NETWORKS & SYSTEMS》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115114578A (zh) * | 2021-03-19 | 2022-09-27 | Aptiv技术有限公司 | 确定对象传感器的噪声统计的方法 |
CN115460047A (zh) * | 2022-08-02 | 2022-12-09 | 北京睿信丰科技有限公司 | 复杂信道下快速跟踪频偏和距离的方法、装置及电子设备 |
CN115460047B (zh) * | 2022-08-02 | 2024-06-11 | 北京睿信丰科技有限公司 | 复杂信道下快速跟踪频偏和距离的方法、装置及电子设备 |
CN117375731A (zh) * | 2023-10-27 | 2024-01-09 | 中国科学院声学研究所南海研究站 | 一种水声通信多途信道时变特性量化分析方法及系统 |
CN117375731B (zh) * | 2023-10-27 | 2024-05-28 | 中国科学院声学研究所南海研究站 | 一种水声通信多途信道时变特性量化分析方法及系统 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111431626A (zh) | 一种基于训练的自适应水声信道跟踪方法 | |
US9799331B2 (en) | Feature compensation apparatus and method for speech recognition in noisy environment | |
US8976913B2 (en) | Adaptive pattern dependent noise prediction on a feed forward noise estimate | |
US20140016688A1 (en) | Equalizer and detector arrangement employing joint entropy-based calibration | |
WO2014173919A1 (en) | Estimation of reliability in speaker recognition | |
Kim | Nonstationary environment compensation based on sequential estimation | |
CN103500455B (zh) | 一种基于无偏有限冲击响应滤波器(ufir)的改进机动目标跟踪方法 | |
CN105355198A (zh) | 一种基于多重自适应的模型补偿语音识别方法 | |
CN111816200A (zh) | 一种基于时频域二值掩膜的多通道语音增强方法 | |
CN114519728A (zh) | 未知噪声下的变分贝叶斯强跟踪phd多目标跟踪方法 | |
CN108444725B (zh) | 一种针对大数据的快速噪声滤除方法 | |
CN115118251A (zh) | 一种基于新息序列的模糊修正滤波方法 | |
CN101127120B (zh) | 一种有效抑制模板漂移的目标跟踪算法 | |
CN102075295B (zh) | 基于通讯功率约束的以状态估计为目标的信号编解码方法 | |
Wang et al. | Training-based adaptive channel tracking for correlated underwater acoustic channels | |
CN113794489B (zh) | 一种通信抗强相关干扰的方法 | |
CN110110711A (zh) | 一种噪声信道下的迭代学习控制系统输入信号估计方法 | |
Kulikova | Square-root adaptive wave filtering for marine vessels | |
CN110514209A (zh) | 一种交互式多模型组合导航方法 | |
CN112636719B (zh) | 数据丢失和信道噪声干扰下的ilc系统输入信号滤波方法 | |
CN110517297B (zh) | 基于误差-模糊分解的一阶伪贝叶斯机动目标跟踪方法 | |
Gabrea | An adaptive Kalman filter for the enhancement of speech signals in colored noise | |
CN114739383B (zh) | 一种基于多目标不确定性改进的gm-phd滤波方法 | |
Gabrea | Adaptive Kalman filtering for speech signal recovery in colored noise | |
Gabrea | Robust adaptive kalman filter for speech signal recovery in colored noise |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20200717 |