CN105721086B - 基于ukfnn的无线信道场景识别方法 - Google Patents

基于ukfnn的无线信道场景识别方法 Download PDF

Info

Publication number
CN105721086B
CN105721086B CN201610141612.1A CN201610141612A CN105721086B CN 105721086 B CN105721086 B CN 105721086B CN 201610141612 A CN201610141612 A CN 201610141612A CN 105721086 B CN105721086 B CN 105721086B
Authority
CN
China
Prior art keywords
mrow
mtd
node
msub
value
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.)
Active
Application number
CN201610141612.1A
Other languages
English (en)
Other versions
CN105721086A (zh
Inventor
李太福
姚立忠
黄迪
梁晓东
周伟
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Yangzhou Yuan Electronic Technology Co Ltd
Original Assignee
Chongqing University of Science and Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Chongqing University of Science and Technology filed Critical Chongqing University of Science and Technology
Priority to CN201610141612.1A priority Critical patent/CN105721086B/zh
Publication of CN105721086A publication Critical patent/CN105721086A/zh
Application granted granted Critical
Publication of CN105721086B publication Critical patent/CN105721086B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B17/00Monitoring; Testing
    • H04B17/30Monitoring; Testing of propagation channels
    • H04B17/391Modelling the propagation channel
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B17/00Monitoring; Testing
    • H04B17/30Monitoring; Testing of propagation channels

Landscapes

  • Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Mobile Radio Communication Systems (AREA)

Abstract

本发明公开了一种基于UKFNN的无线信道场景识别方法,包括以下步骤:采集连续路段的信道数据作为训练样本;确定分段数目并得到连续信道数据分段模型;利用AP算法对每一段数据进行区域划分;将训练样本进行标记,建立连续信道数据分区模型;将待测数据带入分段模型,判断出属于哪一路段;将待测数据带入分区模型,判断出属于哪一区域。有益效果:模型建立复杂性低和计算量少,提高模型的精度;采用AP算法进行聚类,无需指定聚类数目,使信道识别方法更加简单,构建模型更加容易,能准确的识别出采集自该连续路段测试数据属于哪一路段、哪一区域。

Description

基于UKFNN的无线信道场景识别方法
技术领域
本发明涉及无线信道中模式识别领域,具体的说是一种基于UKFNN(unscentedKalman filter artificial neural network,无迹卡尔曼滤波神经网络)的无线信道场景识别方法。
背景技术
移动通信产业一直以惊人的速度迅猛发展,已成为带动全球经济发展的主要高科技产业之一,并对人类生活及社会发展产生了巨大的影响。在移动通信中,发送端和接收端之间通过电磁波来传输信号,我们可以想象两者之间有一些看不见的电磁通路,并把这些电磁通路称为无线信道。无线信道与周围的环境密切相关,不同环境下的无线信道具有一些差异化的特征。如何发现并提取这些特征并将其应用于优化无线网络,是当前的一个研究热点。
无线通信数据信息主要是以无线电波为载体通过无线信道来传输的。由于无线信道所在环境复杂多变,容易出现波的多径传播,例如出现直射、反射、散射等现象;同时,电磁波也会在各路径相互作用下发生多径衰落现象。为更好的反映实际环境中信号的传输规律,为无线网络的规划优化、无线系统的设计、测试和定型等提供重要的参考依据,提出建设无线信道模型。
但是目前无线信道建模方法可以分为统计性建模、确定性建模以及半确定性建模三种方法,但其模型复杂度高、参数获取难度大,难以满足模型精度。
人工神经网络(ANN)以其强大的非线性逼近能力,具有不依赖精确建模过程以及任意逼近非线性映射的特点。ANN在处理复杂系统的建模问题上显示出独特的优越性,被广泛用于模式识别领域。然而目前,很少有人应用ANN来对无线信道特征进行建模,且ANN只对输入输出变量进行简单的静态映射,是一种静态建模方法,对无线信道特征建模效果有限。
发明内容
针对上述问题,本发明提供了一种基于UKFNN的无线信道场景识别方法,采用无迹卡尔曼滤波神经网络方法,即:UKFNN,对无线信道状态参数估计建立动态实时滤波效果的分段模型和分区模型,使其能够反映信道的实际分布情况,从而实现对无线信道的路段、区域的识别,来解决无线信道建模复杂度高、难度大、精度低等问题,将真实信道数据进行智能分段和分区。
为达到上述目的,本发明采用的具体技术方案如下:
一种基于UKFNN的无线信道场景识别方法,包括以下步骤:
S1:采集连续路段的信道数据作为训练样本;
S2:将所得的训练样本平均分为l段,l=2,3,4,5,…,分别利用无迹卡尔曼滤波神经网络进行建模,根据建模效果确定分段数目并得到连续信道数据分段模型;
S3:利用AP算法对每一段数据进行区域划分;
S4:将所得的训练样本按照步骤S3所划分出的区域进行标记,并利用UKFNN无迹卡尔曼滤波神经网络进行区域特征建模,得到连续信道数据分区模型;
S5:将待测数据带入步骤S2所训练出的连续信道数据分段模型,判断出属于哪一路段;
S6:将待测数据带入步骤S4所训练出的连续信道数据分区模型,从而判断出属于哪一区域。
进一步描述,步骤S2中得到连续信道数据分段模型的具体步骤如下:
S21:利用霍特林变换将复数形式的信道数据转换为实数域数据;
S22:利用主成份分析法对步骤S21转换后的信道数据进行降维处理;
S23:对步骤S21和S22处理后的数据进行归一化处理,得到新数据Xm×N,其中,m为变量个数,N为样本数目;
S24:利用UKFNN无迹卡尔曼滤波神经网络对信道数据进行建模,得到无线信道模型;
S25:定义归类准则,并以无迹卡尔曼滤波神经网络的输出变量对输入样本进行归类,对信道分段场景进行识别。
再进一步描述,步骤S21中霍特林变换为:
矩阵A为复数形式的信道数据,寻求正交矩阵Q,使得QφAQT=diag(λ12,…λn),从而得到经霍特林变换后的实数矩阵B=QA,其中,λiA为矩阵A的特征值,1≤i≤n,φA为矩阵A的协方差矩阵,φA=E{(A-E(A))(A-E(A))T},E(A)为矩阵A的期望值矩阵。
再进一步描述,步骤S22中降维的具体方法为:
S221:特征中心化,即将矩阵B的每一维数据都减去该维的均值,得到矩阵BB,其中,矩阵BB的均值为0;
S222:计算矩阵BB的协方差矩阵C;
S223:按照特征值大小,选取对应的特征向量,使得矩阵B转换为主元矩阵t;
S224:确定主元个数k,从而得到新的数据集;
如何确定主元个数是降维的关键,步骤S224中采用累积贡献率CPV来确定主元个数k,即:
第i个主元的贡献率为:
前k个主元的累积贡献率为:式中,λiC为协方差矩阵C中第i个特征值,
若CPV(i)大于设定值CL,则前k个主元可替代原始数据。
优选地,设定值CL=85%。
再进一步地:步骤S23中归一化处理的具体方法为:
式中,x′i为归一化处理之前的输入变量,x″i为归一化处理之后的变量,x'min为归一化处理前输入变量x′i的最小值,x'max为归一化处理前输入变量x′i的最大值。
再进一步描述,步骤S25中定义的归类准则为:
式中,l为场景数目,y为模型输出值,Cla为所属的路段。
再进一步描述,步骤S3中AP算法具体步骤为:
S31::算法初始化,设置一个最大迭代次数K,计算点与点之间的相似度,构成矩阵s,对参考度P赋初始值;
S32:计算样本点间的吸引度值R(u,v)、归属度值A(u,v)
其中R(u,v)表示节点v对于节点u的吸引度值;s(u,v)节点u和节点v之间的聚类中心的合适度;a(u,v')表示v'对于u的归属度值;s(u,v')节点u和节点v'之间的聚类中心的合适度;R(v,v)表示节点v对于节点v的吸引度值;R(u',v)表示节点v对于节点u'的吸引度值;其中,v=!v',u=!u';S33:吸引度和归属度进行加权更新;
每一系迭代,吸引度和归属度都要与上一次迭代的吸引度和归属度进行加权更新,公式为:
其中,λ表示收敛系数,主要作用是用来调节算法的收敛速度以及迭代过程的稳定性;Ru+1(u,v)第u+1次迭代后节点v对于节点u的吸引度值;Ru(u,v)第u次迭代后节点v对于节点u的吸引度值;表示上一更新中,第u+1次迭代后节点v对于节点u的吸引度值;Au+1(u,v)表示第u+1次迭代后节点v对于节点u的归属度值;Au(u,v)表示第u次迭代后节点v对于节点u的归属度值;表示上一更新中,第u+1次迭代后节点v对于节点u的归属度值;
S34:当R(v,v)+R(v,v)>0时,为一个聚类中心,如果迭代次数超过设定值K或聚类中心在一定迭代次数中不再改变,则终止计算,以确定类中心和各类的样本点;否则,返回S32。
再进一步描述,步骤S2和步骤S4中无迹卡尔曼滤波神经网络为三层神经网络,其中,隐含层传递函数为S型函数,输出层传递函数为Purelin函数,该三层神经网络的函数表达式为:
式中,F1为输入层与隐含层之间的连接函数,F2为隐含层与输出层之间的连接函数,bn为输入层和隐含层之间的阈值,b1i为隐含层和输出层间的阈值,xk为输入样本,wij为从神经元j到神经元i的连接权值,为输入层与隐含层之间的连接权值,为隐含层与输出层之间的连接权值,m为输入神经元个数,n为隐含层神经元数目,k为0-10之间的常数。隐含层与神经元i的连接权值。
步骤S2和步骤S4中建模的具体方法为:通过无迹卡尔曼滤波对神经网络的权值、阈值进行估计,将神经网络的权值、阈值作为无迹卡尔曼滤波的状态变量,神经网络的输出作为无迹卡尔曼滤波的测量变量,从而得模型;
步骤S2中确定分段数的方法为:
分别计算训练样本平分成l(l=2,3,4,5,…)段对应的无线信道模型预测准确率ER,比较模型预测准确率ER,确定分段数目以及对应的分段模型,其中预测准确率ER计算公式如下:
其中Me表示预测准确的数据量,M表示指定数据总量。
本发明的有益效果:采用霍特林变换(K-L)变换、主成份分析法(PCA)进行降维等数据预处理,降低了无线信道建立模型结构的复杂性和计算量,提高模型的精度;采用UKFNN对连续路段的无线信道特征进行分段建模,实现对真实数据的自动分段,在此基础上,结合AP算法和UKFNN对连续路段的无线信道特征建立的分区模型,对各段自动进行区域划分,克服了人工神经网络的静态建模方法的缺陷,且AP算法,无需指定聚类数目,使方法更加简单,构建模型更加容易,当输入真实数据时,能准确的识别出采集自该连续路段测试数据属于哪一路段、哪一区域。
附图说明
图1是本发明的方法流程图;
图2是本发明AP聚类算法流程图;
图3是连续路段自动分区图;
图4是真实信道测试数据分段识别统计图;
图5是第一组测试数据区域识别统计图;
图6是第二组测试数据区域识别统计图。
具体实施方式
下面结合附图对本发明的具体实施方式以及工作原理作进一步详细说明。
实施例:
从图1可以看出,一种基于UKFNN的无线信道场景识别方法,包括以下步骤:
S1:采集连续路段的信道数据作为训练样本,;
S2:将所得的训练样本平均分为l段,l=2,3,4,5,…,分别利用无迹卡尔曼滤波神经网络进行建模,根据建模效果确定分段数目并得到连续信道数据分段模型;
步骤S2中得到连续信道数据分段模型的具体步骤如下:
S21:利用霍特林变换将复数形式的信道数据转换为实数域数据;
霍特林变换(K-L)为:
将每个信道数据看作一个二元有序实数对,利用K-L变换将其转化为一维实数。
设A=(a,b)T为复数形式的信道数据,其中,a为复数的实部,b为复数的虚部。
A的协方差矩阵其中,E(A)为矩阵A的期望值矩阵,
φA为一个正定矩阵,因此,存在正交矩阵Q=(q1,q2,qn…)T,使得QφAQT=diag(λ12,…λn),假设λ1≥λ2≥…≥λn≥0,则:φAqi=λiqi,可知qi与λiA是协方差矩阵φA的特征向量与特征值,且满足qi Tqi=1。
此时利用正交矩阵Q对样本进行K-L变换,从而得到变换后的实数矩阵B=QA,B为实数矩阵。
S22:利用主成份分析法对步骤S21转换后的信道数据进行降维处理;
S221:特征中心化,即将矩阵B的每一维数据都减去该维的均值,得到矩阵BB,其中,矩阵BB的均值为0;
S222:计算矩阵BB的协方差矩阵C;
S223:按照特征值大小,选取对应的特征向量,使得矩阵B转换为主元矩阵t,即:
S224:确定主元个数k,从而得到新的数据集;
如何确定主元个数k是降维的关键,本发明采用累积贡献率CPV来确定主元个数k,即:
第i个主元的贡献率为:
前k个主元的累积贡献率为:式中,λi为协方差矩阵C中第i个特征值,
若CPV(i)大于设定值CL,则前k个主元可替代原始数据。
作为一种优选的技术方案,本实施例中设定值CL=85%。
降维处理能够降低模型结构的复杂性和计算量,提高模型的精度。
S23:对步骤S21和S22处理后的数据进行归一化处理,得到新数据Xm×N,其中,m为变量个数,N为样本数目;
归一化处理的具体方法为:
式中,x′i为归一化处理之前的输入变量,x″i为归一化处理之后的变量,x'min为归一化处理前输入变量x′i的最小值,x'max为归一化处理前输入变量x′i的最大值。
S24:利用UKFNN无迹卡尔曼滤波神经网络对信道数据进行建模,得到无线信道模型;
无迹卡尔曼滤波神经网络为三层神经网络,其中,隐含层传递函数为S型函数,输出层传递函数为Purelin函数,该三层神经网络的函数表达式为:
式中,F1为输入层与隐含层之间的连接函数,F2为隐含层与输出层之间的连接函数,bn为输入层和隐含层之间的阈值,b1i为隐含层和输出层间的阈值,xk为输入样本,wij为从神经元j到神经元i的连接权值,为输入层与隐含层之间的连接权值,为隐含层与输出层之间的连接权值,m为输入神经元个数,n为隐含层神经元数目,k为0-10之间的常数,在本实施例中,m=26,n=13即UKFNN的结构为26-13-1。隐含层与神经元i的连接权值。
进一步地,步骤S24中建模的具体方法为:通过无迹卡尔曼滤波对神经网络的权值、阈值进行估计,将神经网络的权值、阈值作为无迹卡尔曼滤波的状态变量,神经网络的输出作为无迹卡尔曼滤波的测量变量,从而得到分段数l(l=2,3,4,5,…)的无线信道模型;
采用无迹卡尔曼滤波神经网络建立无线信道模型时,无迹卡尔曼滤波的状态初始协方差、均值,以及无迹卡尔曼滤波中UT变换中扩展因子的取值范围均在0-1的范围内随机给定。其中UT变换的参数设定为k=0;α=0.06;β=4。
再计算训练样本平分成l(l=2,3,4,5,…)段对应的无线信道模型预测准确率ER,比较模型预测准确率ER,确定分段数目以及分段模型,其中预测准确率ER计算公式如下:
其中Me表示预测准确的数据量,M表示指定数据总量。
表1训练样本不同分段数的预测准确率
表1是训练样本不同分段数的模型预测准确率,根据分段模型准确率ER以及分段数目过多或过少都会对影响区域划分的原则,本实施例中,根据模型预测准确率,将150m连续路段无线信道数据平均分成了3段。
S25:定义归类准则,并以无迹卡尔曼滤波神经网络的输出变量对输入样本进行归类,对信道分段场景进行识别。
定义的归类准则为:
式中,l为场景数目,y为模型输出值,Cla为所属的路段。
由步骤S24得,场景数目l=3,则归类准则为:
S3:利用AP算法对每一段数据进行区域划分;AP算法是一种兴起的聚类算法,它根据各个数据点之间的相似度进行聚类,且对相似度的对称性没有明确要求,同时也不需要事先指定聚类数量,而是将已知数据点都作为潜在的聚类中心处理。
步骤S3中AP算法具体步骤为:
S31:算法初始化,设置一个最大迭代次数K,计算点与点之间的相似度,构成矩阵s,对参考度P赋初始值;
S32:计算样本点间的吸引度值R(u,v)、归属度值A(u,v)
其中R(u,v)表示节点v对于节点u的吸引度值;s(u,v)节点u和节点v之间的聚类中心的合适度;a(u,v')表示v'对于u的归属度值;s(u,v')节点u和节点v'之间的聚类中心的合适度;R(v,v)表示节点v对于节点v的吸引度值;R(u',v)表示节点v对于节点u'的吸引度值;其中,v=!v',u=!u';
S33:吸引度和归属度进行加权更新;
每一系迭代,吸引度和归属度都要与上一次迭代的吸引度和归属度进行加权更新,公式为:
其中,λ表示收敛系数,主要作用是用来调节算法的收敛速度以及迭代过程的稳定性;Ru+1(u,v)第u+1次迭代后节点v对于节点u的吸引度值;Ru(u,v)第u次迭代后节点v对于节点u的吸引度值;表示上一更新中,第u+1次迭代后节点v对于节点u的吸引度值;Au+1(u,v)表示第u+1次迭代后节点v对于节点u的归属度值;Au(u,v)表示第u次迭代后节点v对于节点u的归属度值;表示上一更新中,第u+1次迭代后节点v对于节点u的归属度值;
S34:当R(v,v)+R(v,v)>0时,为一个聚类中心,如果迭代次数超过设定值K或聚类中心在一定迭代次数中不再改变,则终止计算,以确定类中心和各类的样本点;否则,返回S32。
表2样本的路段、分区情况对照表
路段1 路段2 路段3
区域1 (1,702) (5901,7154) (11801,12894)
区域2 (703,1394) (7155,8223) (12895,14771)
区域3 (1395,2056) (8224,9233) (14772,15926)
区域4 (2057,2545) (9234,11002) (15927,16937)
区域5 (2546,3168) (11002,11800) (16938,17700)
区域6 (3169,3819)
区域7 (3820,4624)
区域8 (4625,5135)
区域9 (5136,5900)
从图3和表2可以看到路段1自动分为9个区域,路段2自动分为5个区域,路段3自动分为5个区域。
S4:从表3可以看出,将所得的3段数据训练样本按照步骤S3所划分出的区域进行标记,并通过无迹卡尔曼滤波对神经网络的权值、阈值进行估计,将神经网络的权值、阈值作为无迹卡尔曼滤波的状态变量,神经网络的输出作为无迹卡尔曼滤波的测量变量,从而得分区模型;其中无迹卡尔曼滤波神经网络和步骤S24相同;
S5:将第一组测试数据和第二组测试数据带入步骤S2所训练出的连续信道数据分段模型,从图4可以看出,第一组测试数据所属第1路段的样本点最多,则第一组测试数据属于第1路段,第二组测试数据属于第1路段。
S6:将第一组测试数据和第二组测试数据带入步骤S4所训练出的连续信道数据分区模型,其中表3为连续路段信道数据按区域标记方式,通过图5可以看出,第一组测试数据属于第1路段的第3区域,从图6可以看出,第二组数据属于第1路段的第9区域。也就是说第一组测试数据是来自第1路段的第3区域,第二组数据采集自第1路段的第9区域。
表3连续路段信道数据按区域标记表
通过上述方法,采用霍特林变换、(K-L)变换主成份分析法(PCA)进行降维等数据预处理,降低了无线信道建立模型结构的复杂性和计算量,提高模型的精度;采用UKFNN对连续路段的无线信道特征进行分段建模,实现对真实数据的自动分段,在此基础上,结合AP算法和UKFNN对连续路段的无线信道特征建立的分区模型,对各段自动进行区域划分,克服了人工神经网络的静态建模方法的缺陷,且AP算法,无需指定聚类数目,使方法更加简单,构建模型更加容易,当输入真实数据时,能准确的识别出采集自该连续路段测试数据属于哪一路段、哪一区域。
应当指出的是,上述说明并非是对本发明的限制,本发明也并不仅限于上述举例,本技术领域的普通技术人员在本发明的实质范围内所做出的变化、改性、添加或替换,也应属于本发明的保护范围。

Claims (8)

1.一种基于UKFNN的无线信道场景识别方法,其特征在于包括以下步骤:
S1:采集连续路段的信道数据作为训练样本;
S2:将所得的训练样本平均分为l段,l=2,3,4,5,…,分别利用无迹卡尔曼滤波神经网络进行建模,根据建模效果确定分段数目并得到连续信道数据分段模型;
具体步骤如下:
S21:利用霍特林变换将复数形式的信道数据转换为实数域数据;
S22:利用主成份分析法对步骤S21转换后的信道数据进行降维处理;
步骤S22中降维的具体方法为:
S221:特征中心化,即将矩阵B的每一维数据都减去该维的均值,得到矩阵BB,其中,矩阵BB的均值为0;
S222:计算矩阵BB的协方差矩阵C;
S223:按照特征值大小,选取对应的特征向量,使得矩阵B转换为主元矩阵t;
S224:确定主元个数k,从而得到新的数据集;
S23:对步骤S21和S22处理后的数据进行归一化处理,得到新数据Xm×N,其中,m为变量个数,N为样本数目;
S24:利用无迹卡尔曼滤波神经网络对信道数据进行建模,得到无线信道模型;
S25:定义归类准则,并以无迹卡尔曼滤波神经网络的输出变量对输入样本进行归类,对信道分段场景进行识别;
S3:利用AP算法对每一段数据进行区域划分;
S4:将所得的训练样本按照步骤S3所划分出的区域进行标记,并利用无迹卡尔曼滤波神经网络进行区域特征建模,得到连续信道数据分区模型;
S5:将待测数据带入步骤S2所训练出的连续信道数据分段模型,判断出属于哪一路段;
S6:将待测数据带入步骤S4所训练出的连续信道数据分区模型,从而判断出属于哪一区域。
2.根据权利要求1所述的基于UKFNN的无线信道场景识别方法,其特征在于,步骤S21中霍特林变换为:
矩阵A为复数形式的信道数据,寻求正交矩阵Q,使得QφAQT=diag(λ12,…λn),从而得到经霍特林变换后的实数矩阵B=QA,其中,λiA为矩阵A的特征值,1≤i≤n,φA为矩阵A的协方差矩阵,φA=E{(A-E(A))(A-E(A))T},E(A)为矩阵A的期望值矩阵。
3.根据权利要求1所述的基于UKFNN的无线信道场景识别方法,其特征在于,步骤S224中采用累积贡献率CPV来确定主元个数k,即:
第i个主元的贡献率为:
前k个主元的累积贡献率为:式中,λiC为协方差矩阵C中第i个特征值,
若CPV(i)大于设定值CL,则前k个主元可替代原始数据。
4.根据权利要求1所述的基于UKFNN的无线信道场景识别方法,其特征在于,步骤S23中归一化处理的具体方法为:
式中,xi'为归一化处理之前的输入变量,xi”为归一化处理之后的变量,x'min为归一化处理前输入变量xi'的最小值,x'max为归一化处理前输入变量xi'的最大值。
5.根据权利要求1所述的基于UKFNN的无线信道场景识别方法,其特征在于,步骤S25中定义的归类准则为:
<mrow> <mi>C</mi> <mi>l</mi> <mi>a</mi> <mo>=</mo> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mn>1</mn> </mtd> <mtd> <mrow> <mi>y</mi> <mo>&amp;le;</mo> <mfrac> <mn>4</mn> <mn>3</mn> </mfrac> </mrow> </mtd> </mtr> <mtr> <mtd> <mn>2</mn> </mtd> <mtd> <mrow> <mfrac> <mn>4</mn> <mn>3</mn> </mfrac> <mo>&lt;</mo> <mi>y</mi> <mo>&amp;le;</mo> <mfrac> <mn>7</mn> <mn>3</mn> </mfrac> </mrow> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mo>.</mo> </mtd> <mtd> <mo>.</mo> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>l</mi> <mo>-</mo> <mn>1</mn> </mrow> </mtd> <mtd> <mrow> <mfrac> <mrow> <mn>3</mn> <mrow> <mo>(</mo> <mi>l</mi> <mo>-</mo> <mn>2</mn> <mo>)</mo> </mrow> <mo>+</mo> <mn>1</mn> </mrow> <mn>3</mn> </mfrac> <mo>&lt;</mo> <mi>y</mi> <mo>&amp;le;</mo> <mfrac> <mrow> <mn>3</mn> <mrow> <mo>(</mo> <mi>l</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <mo>+</mo> <mn>1</mn> </mrow> <mn>3</mn> </mfrac> </mrow> </mtd> </mtr> <mtr> <mtd> <mi>l</mi> </mtd> <mtd> <mrow> <mi>y</mi> <mo>&gt;</mo> <mfrac> <mrow> <mn>3</mn> <mrow> <mo>(</mo> <mi>l</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <mo>+</mo> <mn>1</mn> </mrow> <mn>3</mn> </mfrac> </mrow> </mtd> </mtr> </mtable> </mfenced> </mrow>
式中,l为场景数目,y为模型输出值,Cla为所属的路段。
6.根据权利要求1所述的基于UKFNN的无线信道场景识别方法,其特征在于,步骤S3中AP算法具体步骤为:
S31:算法初始化,设置一个最大迭代次数K,计算点与点之间的相似度,构成矩阵s,对参考度P赋初始值;
S32:计算样本点间的吸引度值R(u,v)、归属度值A(u,v)
<mrow> <mi>R</mi> <mrow> <mo>(</mo> <mi>u</mi> <mo>,</mo> <mi>v</mi> <mo>)</mo> </mrow> <mo>=</mo> <mi>s</mi> <mrow> <mo>(</mo> <mi>u</mi> <mo>,</mo> <mi>v</mi> <mo>)</mo> </mrow> <mo>-</mo> <munder> <mrow> <mi>m</mi> <mi>a</mi> <mi>x</mi> </mrow> <mrow> <msup> <mi>v</mi> <mo>&amp;prime;</mo> </msup> <mo>&amp;NotEqual;</mo> <mi>v</mi> </mrow> </munder> <mo>{</mo> <mi>a</mi> <mrow> <mo>(</mo> <mi>u</mi> <mo>,</mo> <msup> <mi>v</mi> <mo>&amp;prime;</mo> </msup> <mo>)</mo> </mrow> <mo>+</mo> <mi>s</mi> <mrow> <mo>(</mo> <mi>u</mi> <mo>,</mo> <msup> <mi>v</mi> <mo>&amp;prime;</mo> </msup> <mo>)</mo> </mrow> <mo>}</mo> </mrow>
<mrow> <mi>A</mi> <mrow> <mo>(</mo> <mi>u</mi> <mo>,</mo> <mi>v</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <mi>min</mi> <mo>{</mo> <mn>0</mn> <mo>,</mo> <mi>R</mi> <mrow> <mo>(</mo> <mi>v</mi> <mo>,</mo> <mi>v</mi> <mo>)</mo> </mrow> <mo>+</mo> <munder> <mo>&amp;Sigma;</mo> <mrow> <mi>i</mi> <mo>&amp;NotElement;</mo> <mo>{</mo> <mi>i</mi> <mo>,</mo> <mi>k</mi> <mo>}</mo> </mrow> </munder> <mi>m</mi> <mi>a</mi> <mi>x</mi> <mo>{</mo> <mn>0</mn> <mo>,</mo> <mi>R</mi> <mrow> <mo>(</mo> <msup> <mi>u</mi> <mo>&amp;prime;</mo> </msup> <mo>,</mo> <mi>v</mi> <mo>)</mo> </mrow> <mo>}</mo> <mo>}</mo> <mo>,</mo> </mrow> </mtd> <mtd> <mrow> <mi>u</mi> <mo>&amp;NotEqual;</mo> <mi>v</mi> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <munder> <mo>&amp;Sigma;</mo> <mrow> <mi>u</mi> <mo>&amp;NotEqual;</mo> <mi>v</mi> </mrow> </munder> <mi>m</mi> <mi>a</mi> <mi>x</mi> <mo>{</mo> <mn>0</mn> <mo>,</mo> <mi>R</mi> <mrow> <mo>(</mo> <msup> <mi>u</mi> <mo>&amp;prime;</mo> </msup> <mo>,</mo> <mi>v</mi> <mo>)</mo> </mrow> <mo>}</mo> <mo>,</mo> </mrow> </mtd> <mtd> <mrow> <mi>u</mi> <mo>=</mo> <mi>v</mi> </mrow> </mtd> </mtr> </mtable> </mfenced> </mrow>
其中R(u,v)表示节点v对于节点u的吸引度值;s(u,v)节点u和节点v之间的聚类中心的合适度;a(u,v')表示v'对于u的归属度值;s(u,v')节点u和节点v'之间的聚类中心的合适度;R(v,v)表示节点v对于节点v的吸引度值;R(u',v)表示节点v对于节点u'的吸引度值;其中,v=!v',u=!u';S33:吸引度和归属度进行加权更新;
每一次迭代,吸引度和归属度都要与上一次迭代的吸引度和归属度进行加权更新,公式为:
<mrow> <msub> <mi>R</mi> <mrow> <mi>u</mi> <mo>+</mo> <mn>1</mn> </mrow> </msub> <mrow> <mo>(</mo> <mi>u</mi> <mo>,</mo> <mi>v</mi> <mo>)</mo> </mrow> <mo>=</mo> <mi>&amp;lambda;</mi> <mo>&amp;CenterDot;</mo> <msub> <mi>R</mi> <mi>u</mi> </msub> <mrow> <mo>(</mo> <mi>u</mi> <mo>,</mo> <mi>v</mi> <mo>)</mo> </mrow> <mo>+</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <mi>&amp;lambda;</mi> <mo>)</mo> </mrow> <mo>&amp;CenterDot;</mo> <msubsup> <mi>R</mi> <mrow> <mi>u</mi> <mo>+</mo> <mn>1</mn> </mrow> <mrow> <mi>o</mi> <mi>l</mi> <mi>d</mi> </mrow> </msubsup> <mrow> <mo>(</mo> <mi>u</mi> <mo>,</mo> <mi>v</mi> <mo>)</mo> </mrow> <mo>,</mo> <mi>&amp;lambda;</mi> <mo>&amp;Element;</mo> <mo>&amp;lsqb;</mo> <mn>0.5</mn> <mo>,</mo> <mn>1</mn> <mo>)</mo> </mrow>
<mrow> <msub> <mi>A</mi> <mrow> <mi>u</mi> <mo>+</mo> <mn>1</mn> </mrow> </msub> <mrow> <mo>(</mo> <mi>u</mi> <mo>,</mo> <mi>v</mi> <mo>)</mo> </mrow> <mo>=</mo> <mi>&amp;lambda;</mi> <mo>&amp;CenterDot;</mo> <msub> <mi>A</mi> <mi>u</mi> </msub> <mrow> <mo>(</mo> <mi>u</mi> <mo>,</mo> <mi>v</mi> <mo>)</mo> </mrow> <mo>+</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <mi>&amp;lambda;</mi> <mo>)</mo> </mrow> <mo>&amp;CenterDot;</mo> <msubsup> <mi>A</mi> <mrow> <mi>u</mi> <mo>+</mo> <mn>1</mn> </mrow> <mrow> <mi>o</mi> <mi>l</mi> <mi>d</mi> </mrow> </msubsup> <mrow> <mo>(</mo> <mi>u</mi> <mo>,</mo> <mi>v</mi> <mo>)</mo> </mrow> <mo>,</mo> <mi>&amp;lambda;</mi> <mo>&amp;Element;</mo> <mo>&amp;lsqb;</mo> <mn>0.5</mn> <mo>,</mo> <mn>1</mn> <mo>)</mo> </mrow>
Ru+1(u,v)第u+1次迭代后节点v对于节点u的吸引度值;Ru(u,v)第u次迭代后节点v对于节点u的吸引度值;表示上一更新中,第u+1次迭代后节点v对于节点u的吸引度值;Au+1(u,v)表示第u+1次迭代后节点v对于节点u的归属度值;Au(u,v)表示第u次迭代后节点v对于节点u的归属度值;表示上一更新中,第u+1次迭代后节点v对于节点u的归属度值;λ表示收敛系数,主要作用是用来调节算法的收敛速度以及迭代过程的稳定性;S34:当R(v,v)+R(v,v)>0时,为一个聚类中心,如果迭代次数超过设定值K或聚类中心在一定迭代次数中不再改变,则终止计算,以确定类中心和各类的样本点;否则,返回S32。
7.根据权利要求1所述的基于UKFNN的无线信道场景识别方法,其特征在于,步骤S2和步骤S4中无迹卡尔曼滤波神经网络为三层神经网络,其中,隐含层传递函数为S型函数,输出层传递函数为Purelin函数,该三层神经网络的函数表达式为:
<mrow> <mi>y</mi> <mo>=</mo> <mi>h</mi> <mrow> <mo>(</mo> <msub> <mi>w</mi> <mi>k</mi> </msub> <mo>,</mo> <msub> <mi>x</mi> <mi>k</mi> </msub> <mo>)</mo> </mrow> <mo>=</mo> <msup> <mi>F</mi> <mn>2</mn> </msup> <mrow> <mo>(</mo> <msubsup> <mi>w</mi> <mi>k</mi> <mn>2</mn> </msubsup> <mo>,</mo> <msup> <mi>F</mi> <mn>1</mn> </msup> <mo>(</mo> <mrow> <msubsup> <mi>w</mi> <mi>k</mi> <mn>1</mn> </msubsup> <mo>,</mo> <msub> <mi>x</mi> <mi>k</mi> </msub> </mrow> <mo>)</mo> <mo>)</mo> </mrow> <mo>=</mo> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>n</mi> </munderover> <mfrac> <msubsup> <mi>w</mi> <mi>i</mi> <mn>2</mn> </msubsup> <mrow> <mn>1</mn> <mo>+</mo> <mi>e</mi> <mo>&amp;lsqb;</mo> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>j</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>m</mi> </munderover> <msub> <mi>w</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> <msub> <mi>x</mi> <mi>i</mi> </msub> <mo>+</mo> <msub> <mi>b</mi> <mrow> <mn>1</mn> <mi>i</mi> </mrow> </msub> <mo>&amp;rsqb;</mo> </mrow> </mfrac> <mo>+</mo> <msub> <mi>b</mi> <mi>n</mi> </msub> </mrow>
式中,F1为输入层与隐含层之间的连接函数,F2为隐含层与输出层之间的连接函数,bn为输入层和隐含层之间的阈值,b1i为隐含层和输出层间的阈值,xk为输入样本,wij为从神经元j到神经元i的连接权值,为输入层与隐含层之间的连接权值,为隐含层与输出层之间的连接权值,m为输入神经元个数,n为隐含层神经元数目,k为0-10之间的常数;隐含层与神经元i的连接权值。
8.根据权利要求1所述的基于UKFNN的无线信道场景识别方法,其特征在于,步骤S2和步骤S4中建模的具体方法为:通过无迹卡尔曼滤波对神经网络的权值、阈值进行估计,将神经网络的权值、阈值作为无迹卡尔曼滤波的状态变量,神经网络的输出作为无迹卡尔曼滤波的测量变量,从而得到分区模型或者分段模型;
步骤S2中确定分段数的方法为:
分别计算训练样本平分成l(l=2,3,4,5,…)段对应的无线信道模型预测准确率ER,比较模型预测准确率ER,确定分段数目以及分段模型,其中预测准确率ER计算公式如下:
<mrow> <mi>E</mi> <mi>R</mi> <mo>=</mo> <mfrac> <msub> <mi>M</mi> <mi>e</mi> </msub> <mi>M</mi> </mfrac> <mo>&amp;times;</mo> <mn>100</mn> <mi>%</mi> </mrow>
其中Me表示预测准确的数据量,M表示指定数据总量。
CN201610141612.1A 2016-03-11 2016-03-11 基于ukfnn的无线信道场景识别方法 Active CN105721086B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610141612.1A CN105721086B (zh) 2016-03-11 2016-03-11 基于ukfnn的无线信道场景识别方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610141612.1A CN105721086B (zh) 2016-03-11 2016-03-11 基于ukfnn的无线信道场景识别方法

Publications (2)

Publication Number Publication Date
CN105721086A CN105721086A (zh) 2016-06-29
CN105721086B true CN105721086B (zh) 2018-05-01

Family

ID=56158575

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610141612.1A Active CN105721086B (zh) 2016-03-11 2016-03-11 基于ukfnn的无线信道场景识别方法

Country Status (1)

Country Link
CN (1) CN105721086B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107734507A (zh) * 2016-08-12 2018-02-23 索尼公司 无线场景识别装置和方法以及无线通信设备和系统
CN107451605A (zh) * 2017-07-13 2017-12-08 电子科技大学 一种基于信道状态信息和支持向量机的简单目标识别方法
CN110233653B (zh) * 2019-06-10 2020-03-20 深圳大学 基于加权集成聚类算法的mimo系统的盲多径识别方法及系统
CN110263870B (zh) * 2019-06-26 2023-04-18 深圳市悦动天下科技有限公司 运动状态识别方法、装置、智能终端和存储介质
CN111835444B (zh) * 2020-06-17 2021-04-27 武汉大学 一种无线信道场景识别方法及系统

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101814160A (zh) * 2010-03-08 2010-08-25 清华大学 一种基于特征聚类的rbf神经网络建模方法
CN102682048A (zh) * 2011-10-28 2012-09-19 刘晓勇 一种自适应吸引传播聚类方法及装置
CN103177289A (zh) * 2013-03-06 2013-06-26 重庆科技学院 一种噪声不确定复杂非线性动态系统的建模方法
CN103218675A (zh) * 2013-05-06 2013-07-24 国家电网公司 一种基于聚类和滑动窗口的短期负荷预测方法
CN103345559A (zh) * 2013-07-10 2013-10-09 重庆科技学院 铝电解过程电解槽工艺能耗的动态演化建模方法
CN104200032A (zh) * 2014-09-05 2014-12-10 山东大学 广义负荷建模中基于时段性的横向时间轴聚类方法
CN104656441A (zh) * 2014-12-29 2015-05-27 重庆科技学院 基于无迹卡尔曼神经网络的天然气净化工艺建模优化方法
CN105045941A (zh) * 2015-03-13 2015-11-11 重庆科技学院 基于无迹卡尔曼滤波的抽油机参数优化方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8346712B2 (en) * 2009-11-24 2013-01-01 King Fahd University Of Petroleum And Minerals Method for identifying hammerstein models

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101814160A (zh) * 2010-03-08 2010-08-25 清华大学 一种基于特征聚类的rbf神经网络建模方法
CN102682048A (zh) * 2011-10-28 2012-09-19 刘晓勇 一种自适应吸引传播聚类方法及装置
CN103177289A (zh) * 2013-03-06 2013-06-26 重庆科技学院 一种噪声不确定复杂非线性动态系统的建模方法
CN103218675A (zh) * 2013-05-06 2013-07-24 国家电网公司 一种基于聚类和滑动窗口的短期负荷预测方法
CN103345559A (zh) * 2013-07-10 2013-10-09 重庆科技学院 铝电解过程电解槽工艺能耗的动态演化建模方法
CN104200032A (zh) * 2014-09-05 2014-12-10 山东大学 广义负荷建模中基于时段性的横向时间轴聚类方法
CN104656441A (zh) * 2014-12-29 2015-05-27 重庆科技学院 基于无迹卡尔曼神经网络的天然气净化工艺建模优化方法
CN105045941A (zh) * 2015-03-13 2015-11-11 重庆科技学院 基于无迹卡尔曼滤波的抽油机参数优化方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
一种新型的基于神经网络的无线信道模型;周银东等;《电讯技术》;20050428;正文第143页-第146页 *

Also Published As

Publication number Publication date
CN105721086A (zh) 2016-06-29

Similar Documents

Publication Publication Date Title
CN105721086B (zh) 基于ukfnn的无线信道场景识别方法
CN109814066B (zh) 基于神经网络学习的rssi室内定位测距方法、室内定位平台
CN106203478A (zh) 一种用于智能电表大数据的负荷曲线聚类方法
CN111935629B (zh) 一种基于环境特征迁移的自适配定位方法
CN112685504B (zh) 一种面向生产过程的分布式迁移图学习方法
CN106656357B (zh) 一种工频通信信道状态评估系统和方法
WO2023169589A1 (zh) 一种基于对抗网络与长短期记忆网络的预测信道建模方法
CN108182500A (zh) 基于加速寿命试验的弹药贮存可靠性预测方法
CN107333294A (zh) 一种结合AdaBoost和支持向量机的链路质量预测方法
CN107688863A (zh) 自适应迭代加强的短期风速高精度组合预测方法
CN110502989A (zh) 一种小样本高光谱人脸识别方法及系统
CN113297936A (zh) 一种基于局部图卷积网络的排球群体行为识别方法
CN112418476A (zh) 一种超短期电力负荷预测方法
CN111797979A (zh) 基于lstm模型的振动传递系统
CN112562312A (zh) 一种基于融合特征的GraphSAGE交通路网数据预测的方法
CN105653728A (zh) 基于双层模型体系的无线传感网的数据采集方法
CN113095354B (zh) 基于辐射源特征子空间知识的未知雷达目标识别方法
CN113450562B (zh) 一种基于聚类和图卷积网络的路网交通状态判别方法
CN117636183A (zh) 一种基于自监督预训练的小样本遥感图像分类方法
CN117034060A (zh) 基于ae-rcnn的洪水分级智能预报方法
CN115099385B (zh) 基于传感器布局优化和自适应Kriging模型的频谱地图构建方法
CN105792232B (zh) 基于ukfnn的无线信道“指纹”特征动态建模方法
CN113055111B (zh) 一种基于贝叶斯优化的信道建模方法及系统
CN108985563A (zh) 一种基于自组织特征映射的机电系统服役动态标记方法
Malekzadeh et al. Non-Gaussian BLE-based indoor localization via Gaussian sum filtering coupled with Wasserstein distance

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
TR01 Transfer of patent right

Effective date of registration: 20200713

Address after: 225006 No. 15 Wenchang East Road, Yangzhou City, Jiangsu Province, No. 4 Building, Dongyuan, Jiangguang Wisdom City

Patentee after: YANGZHOU YUAN ELECTRONIC TECHNOLOGY Co.,Ltd.

Address before: 401331 Shapingba University District, Chongqing City Road, No. 20

Patentee before: Chongqing University of Science & Technology

TR01 Transfer of patent right