基于动态时间规整算法的无线链路质量预测方法
【技术领域】
[0001]本发明涉及一种无线传感器网络链路质量预测方法，特别是利用动态时间规整(Dynamic Time Warping,简称DTW)算法进行两个时间序列相似度分析从而对链路质量作出预测的方法。
【背景技术】
[0002]近年来，随着无线通信网络的快速发展，无线网络已经遍布世界各地。人们对于无线网络的需求和依赖性越来越大。拥有一个稳定可靠的网络系统对于提高工厂以及人们的工作效率有很大的帮助。无线传感网络的节点通常放置在恶劣的环境中，恶劣环境对节点间通信使用的无线电波影响较大，容易造成数据的丢失，数据丢失后一般采用重传机制，该机制虽然在一定程度上可以提高数据传输成功率但也造成能量的浪费，如果可以提前预测出链路质量较好的链路进行传输，则会提高数据传输率，并且可以节省节点的能量消耗。因此，对无线传感网络链路质量预测机制进行研宄是很有必要的。
[0003]目前国内外在无线传感网络链路质量方面进行了大量研宄，目前对无线传感器网络链路质量的预测现有的主要的技术方法有:
[0004]1、模糊神经网络方法，神经网络是模拟人脑结构的思维功能，具有较强的自学习和联想功能，人工干预少，精度较高，对专家知识的利用也较好。但缺点是它不能处理和描述模糊信息，不能很好利用已有的经验知识，特别是学习及问题的求解具有黑箱的特性，其工作不具有可解释性，同时它对样本的要求较高。
[0005]2、采用A-TEC算法，通过利用滑动窗口模型进行粗估计，但由于预测没有更为精确定量的算法，导致预测结果的精确度不够高。
[0006]高泽鹏，程良伦发表的“一种基于PRR的链路质量预测方法”(《传感器学报》，2011,24(7):1070-1074)通过对PRR值进行平滑处理，考虑了链路的不对称性，提出了A-TEC算法，在移动窗口平均模型下对PRR进行预测，虽然起到良好的效果，但是没有从定性和定量上分析算法，同时对节点的计算能力提出了很高的要求，而且执行速度慢，工作效率很低。
[0007]3、孙佩刚，赵海，罗玎玎发表的“无线传感器网络链路通信质量测量研宄”(《通信学报》，2007，28 (10):14-22)提出把几个链路性质指标用于模糊算法当中，并设计了几个对比实验来参考比较各个指标的适用情况，但是没有给出具体的预测算法和计算公式，这使得预测结果没有确定性，没能达到可靠，稳定的效果。
[0008]无线传感网络是由部署在监测区域内大量廉价微型传感器节点组成，通过无线通信方式形成的一个多跳自组织网络系统，是以应用为中心的一类网络，为人们提供与物理世界交互的手段。一方面，应用提出了苛刻的端到端可靠性要求；另一方面，底层的无线通信链路受各种因素的影响(比如电磁波的多径传播效应和背景噪声、在无线网络或其他电磁源中会产生一定量的干扰以及硬件收发器内在的噪声)，导致节点间的通信链路质量会呈现不规则的变化，因此无线链路质量会表现出一定的时变性和动态性，难以实现准确的预测。
【发明内容】
[0009]本发明的目的是要克服上述各种技术方案的缺陷，提供一种可以提高预测结果的可靠性和准确性的预测方法。
[0010]为了实现本发明的目的，所采用的技术方案是提供一种基于动态时间规整算法的无线链路质量预测方法，其无线链路质量是由无线信号的信噪比所表征，包括链路信噪比训练数据采集与处理、待预测链路信噪比序列相似度匹配和链路质量预测，其具体步骤如下:
[0011]1、一种基于动态时间规整算法的无线链路质量预测方法，其特征在于，包括链路信噪比训练数据采集与处理、待预测链路信噪比序列相似度匹配和链路质量预测，其具体步骤如下:
[0012]步骤I，使用无线传感器网络中所规定的A节点按一定周期P连续采集无线传感器网络中规定的B节点发出的无线信号的信噪比，所采集的数据组成一段长度为K的信噪比时间序列，其中，100彡K彡1000，Is彡P彡3s ;
[0013]步骤2，对步骤I所采集的信噪比时间序列进行滤波处理；
[0014]步骤3，存储滤波处理后的信噪比时间序列，标记为训练数据序列T，
[0015]T = (T1, T2, T3, , Tj，
[0016]其中10I^TpTk分别表示经过滤波处理的第1、2、3、K次采集的信噪比；
[0017]步骤4，将长度为K的训练数据序列T切分成(K-M+1)个长度为M的训练子序列(t1； t2, t3，...，tK_M+1)，其中M为小于20的正整数，各个训练子序列的表达式为:
[0018]ti= (T1,..., Tj ；
[0019]t2= [T2,..., TM+1}；
[0020]t3= {T 3,..., ΤΜ+2}；
[0021]......
[0022]tK-1+i — {Τ κ-m+D...) Τκ};
[0023]t2、t3、TK_M+A1U表不第 1、2、3、K-M+1 段训练子序列；Τ Μ、ΤΜ+1、ΤΜ+2、ΤΚ_Μ+^别表示经过滤波处理的第Μ、Μ+1、Μ+2、Κ-Μ+1次采集的信噪比；
[0024]步骤5，使用步骤I中的A节点按相同的周期P连续采集B节点当前时刻的信噪比，并组成长度为M的信噪比时间序列；
[0025]步骤6，对步骤5所采集的长度为M的信噪比时间序列进行滤波处理，并把经过滤波处理的序列标记为当前询问序列Q，
[0026]Q = (Q1, Q2, Q3, , QmI
[0027]式中，Q1, Q2, Q3, Qm分别表示经过滤波处理的当前链路第1、2、3、M次采集的信噪比；
[0028]步骤7，对步骤4中得到的每一个训练子序列，利用动态时间规整算法计算该训练子序列与步骤6中的当前询问序列Q之间的最优路径所对应的累积距离，其具体步骤包括:
[0029]步骤7.1，计算当前询问序列Q和训练子序列每个信噪比之间的相互距离，包括:列出长度均为M的两序列的信噪比值，构造成一个MXM的矩阵网格，计算每个信噪比值之间的距离，得到MXM个距离值，其中距离值的计算公式为:
[0030]Clij= (Q1-Tj)2
[0031]式中，du表示当前询问序列第i个信噪比与训练子序列第j个信噪比的距离值；
[0032]01表示当前询问序列中第i个信噪比数值；
[0033]Tj表示训练子序列中第j个信噪比数值；
[0034]步骤7.2，将步骤7.1求得距离值du按照其下标依次填入MXM矩阵网格相应的位置，其中i为行的位置，j为列的位置；
[0035]步骤7.3，按照路径约束条件，在MXM的网格中寻找出一条最优路径，并将该最优路径所经过的网格中的距离值累加求和，得出该最优路径所对应的累计距离，即该训练子序列与当前询问序列Q之间最优路径对应的累积距离；
[0036]步骤8，对步骤4中(K-M+1)个训练子序列用步骤7同样的方法计算出每个训练子序列与当前询问序列Q之间最优路径对应的累计距离，并进行比较，找出累计距离最小的训练子序列，定义该训练子序列为与当前询问序列相似度最高的训练子序列，记为R ；
[0037]步骤9，将步骤8所得到的训练子序列R按照所采集时间顺序排列的下一段训练子序列设定为预测未来一段时间的链路质量。
[0038]优选地，步骤2所述的滤波方法选用卡尔曼滤波；
[0039]优选地，步骤7.3中所述的路径约束条件包括:
[0040]边界条件:路径必须是从MXM矩阵网格的左下角出发，在右上角结束；
[0041]连续性:路径只能沿着相邻以及对角的网格这三个方向延续。
[0042]最值性:在所形成的最优路径中必须满足其路径的距离值最小。
[0043]优选地，步骤7.3中所述的最优路径是指按照步骤7.3中所述的路径约束条件所形成的一条路径。
[0044]优选地，步骤9中所述的未来一段时间为权利要求1中步骤5所述的采集M个信噪比所使用的时间，即为(MXP) S。
[0045]与现有技术相比，本发明的有益效果如下:
[0046]1.利用卡尔曼滤波算法将采集到的时间序列信息进行去噪声处理，提高预测结果的可靠性，减少了由于噪声产生的预测误差；
[0047]2.克服了传统上利用欧氏距离方法计算序列之间相似性时距离点必须对应的问题，并且对时间序列同步问题不敏感，使得计算结果的可靠性和准确性提高。
【附图说明】
[0048]图1为本发明基本方法流程示意图。
[0049]图2为本发明实施例中收集的100次链路信噪比图。
[0050]图3为本发明实施例中经过噪声处理后的100次链路信噪比图。
[0051]图4为本发明动态时间规整算法原理示意图。
[0052]图5为本发明的【具体实施方式】预测结果比较图。
【具体实施方式】
[0053]下面结合附图及实施例对本发明做进一步详细说明。
[0054]图1为本发明基本方法流程示意图，其无线链路质量是由无线信号的信噪比所表征，由该图可见，本发明包括链路信噪比训练数据采集与处理、待预测链路信噪比序列相似度匹配和链路质量预测三个部分，具体步骤如下:
[0055]步骤I，使用无线传感器网络中所规定的A节点按一定周期P连续采集无线传感器网络中规定的B节点发出的无线信号的信噪比，所采集的数据组成一段长度为K的信噪比时间序列，其中，100彡K彡1000，Is彡P彡3s。在本实施例中，K为100，P为Is即连续采集100次无线信号的信噪比，组成长度为100的信噪比时间序列。图2为本实施例中收集的100次链路信噪比图。
[0056]步骤2，对步骤I所采集的信噪比时间序列进行滤波处理。
[0057]由于所测量得到的时间序列中存在测量误差点，为了减少这些噪声对实验结果产生的影响，因此使用卡尔曼滤波的方法滤除长度为100的信噪比时间序列的噪声，具体滤波过程为:
[0058]在本实施方案中噪声测量的卡尔曼滤波模型用以下等式(I)，(2)描述:
[0059]Xk=AXk-Jwk (I)
[0060]zk= HX k+Vk (2)
[0061]式中:
[0062]Xk是第k次链路信噪比的状态变量；
[0063]A是系统状态转移矩阵，A=l；
[0064]Xlri是第k-Ι次链路信噪比的状态变量；
[0065]H是观测矩阵，H= I;
[0066]Zk是第k次链路信噪比的观测值；
[0067]Wk和V k分别为状态噪声和观测噪声，他们是相互独立的随机噪声，并服从均值为0，方差为I的正态随机分布。
[0068]步骤3，存储滤波处理后的信噪比时间序列，标记为训练数据序列T，
[0069]T = IT1, T2, T3,…T100I，
[0070]其中，10TpI^Ticitl分别表示经过滤