基于链路拥塞强度分布的网络拥塞链路诊断方法及系统与流程

本发明涉及网络拥塞控制领域，更具体地说，涉及基于链路拥塞强度分布的网络拥塞链路诊断方法及系统。

背景技术：

ip网络链路性能推断技术对广大网络用户，特别是网络运营、管理商非常重要。随着ip网络规模的迅速扩大，网络结构日益多样，人工定期检查已不能适应大规模网络的需要。网络层析成像能够在网络内部节点不提供测量协作的情况下，根据端到端的测量结果，间接地估计网络内部链路性能参数，是一种重要的网络测量手段，能直接指导网络管理和网络优化。

在当前类traceroute工具所测数据包含有低转发优先级的icmp报文测量部分以及internet中大部分玩过不提供测量协助的背景下，基于端到端测量的网络层析成像技术被提出来了。网络层析成像将目标网络视为一个黑盒子。自从1996年vardi首次提出在ip网络性能推断中使用类似医学层析扫描(tomograph)技术以来，借助网络tomography技术推断ip网络内部链路性能的方法主要包括两类：第一类方法使用多播方式或多簇单播模拟多播的方式，通过构建ip网络线性方程组求解各链路丢包率，并假设链路性能服从特定分布或具有时空独立性和平稳性等。出于安全等因素，当前ip网络中的大部分路由器对单播的支持度高于多播，时间相关性难以保证，且tomography技术以尽可能少的e2e路径探测推断网络内部链路性能，常导致构建的系统线性方程组因系数矩阵维度过大，求逆计算复杂，甚至导致算法失效。第二类方法将路径及链路性能借助布尔(boolean)二进制代数值0、1进行表示。

h.x.nguyen，v.n.padmanabhan，n.g.duffield等人提出借助不相关的e2e路径探测，推断ip网络最有可能发生拥塞的链路集合，简化了链路性能推断过程。其中，huangx.等人提出测clink算法较不使用先验概率的mcmc(montecarlomarkovchain)算法及使用一致先验概率的scfs(smallestconsistentfailuresetalgorithm)算法在推断性能上有较大程度的提升。

技术实现要素：

本发明要解决的技术问题在于，针对上述目前传统借助主动网络层析成像技术推断ip网络拥塞链路的方法均不能对动态路由ip网络进行准确建模，且忽略网络拥塞强度的客观存在，从而使得推断拥塞链路产生一定误差的技术问题，提出了基于链路拥塞强度分布的网络拥塞链路诊断方法及系统来解决上述问题。

基于链路拥塞强度分布的网络拥塞链路诊断方法，包括：

s1、对待测ip网络中端到端路径进行t时刻快照，获得t时刻待测ip网络中各端到端路径性能及拓扑结构的探测结果，使用快照获得待测ip网络链路的初始概率，初始概率包括链路的初始拥塞概率以及链路的初始拥塞强度概率；

s2、根据t时刻快照所得到的初始概率使用em算法模型计算得到待测ip网络的先验概率，先验概率包括先验拥塞概率和链路的先验拥塞强度概率；

s3、建立待测ip网络中各链路的先验拥塞概率求解的线性方程组；

s4、根据步骤s2得到待测ip网络的先验概率，结合步骤s3得到的线性方程组，计算推断待测ip网络中最有可能发生拥塞的链路集合。

进一步的，所述步骤s2计算先验概率具体包括：

a、根据所述探测结果，获得某一时间片t对应的链路的拥塞概率向量p和链路的拥塞强度概率分布d^t；

b、根据t时刻得到的初始概率以及ip网络路径状态z，计算对数似然ll(d|z，x)关于x的期望，得到最可能发生的链路状态x；

c、根据链路状态x，重新计算每条链路的拥塞强度概率分布，得到新的d^t+1；

d、当迭代到一定次数，或者算法收敛到一定精度，结束算法，得到所述待测网络中的先验链路拥塞强度概率分布；否则，回到步骤b。

进一步的，建立所述线性方程组具体包括：

a、在所述待测ip网络中，将各端到端与其路径的各链路之间的关系使用布尔代数模型表示为：

nc为各端到端路径经过的链路总和，zi为第i条路径的状态变量，xj为第j条链路的状态变量。当zi≥1时，表示第i条路径为拥塞状态，当zi＝0时，表示第i条路径为正常状态；当xj＝1时，表示第j条链路位拥塞状态，当xj＝0时，表示第j条链路为正常状态；dij为路径和链路的关系变量，当dij＝1时，表示第i条路径经过第j条链路，当dij＝0时，表示第i条路径不经过第j条链路；

b、在所述待测ip网络中，将各端到端与其路径的各链路之间的拥塞强度的关系使用线性代数模型表示为：

nc为各端到端路径经过的链路总和，mi为第i条路径的状态变量，yj为第j条链路的状态变量。当mi≥1时，表示第i条路径为拥塞状态，当mi＝0时，表示第i条路径为正常状态；当yj≥1时，表示第j条链路位拥塞状态，当yj＝0时，表示第j条链路为正常状态；dij为路径和链路的关系变量，当dij＝1时，表示第i条路径经过第j条链路，当dij＝0时，表示第i条路径不经过第j条链路；

c、对所述a步骤中的布尔代数模型，找到向量x＝[x1x2…xn]，使得如下条件概率公式得到最大化：

其中z＝[z1z2…zm]表示快照得到的待测ip网络的状态，n表示待测ip网络中链路的个数，m表示待测ip网络中路径的个数，

pg为测量标记为正常链路的集合，pc为测量标记为拥塞链路的集合，pi表示路径i。

进一步的，所述步骤c中最大化概率公式可以化简为：

其中，εr表示去除所有经过正常路径的链路以及确定为正常链路之后所得到的路由矩阵，χ表示判断的拥塞链路集合，pk表示链路k的拥塞概率。

进一步的，在当前推断时刻拥塞路径链路集合中查找权值最高的链路的判断公式如下：

其中，dk表示链路k经过的路径的集合，|dk|即为链路k经过的路径的个数。

基于链路拥塞强度分布的网络拥塞链路诊断系统，包括：处理器及存储设备；所述处理器加载并执行所述存储设备中的指令及数据用于实现上述的任意基于链路拥塞强度分布的网络拥塞链路诊断方法。

与现有技术相比，本发明的优势在于：本发明针对ip网络路由的动态特性，根据多路端到端探测获取的路径性能及拓扑结构的探测结果，对ip网络拥塞链路进行推断，建立em算法模型得到待测ip网络的各链路的拥塞强度概率分布，在此基础上基于该模型对ip网络拥塞链路进行推断，实现了动态路由下ip网络拥塞链路的推断。通过对不同类型及规模ip网络进行模拟实验和仿真实验，在相同场景下，分别利用传统clink和scfs算法及本发明的方法进行拥塞链路推断，实验证明了本发明的方法较clink算法有更高的推断准确性。

附图说明

下面将结合附图及实施例对本发明作进一步说明，附图中：

图1为本发明基于链路拥塞强度分布的网络拥塞链路诊断方法的流程图；

图2为本发明中em算法建模流程图；

图3是本发明中不同拥塞链路比例下clink、scfs、dacli-clink算法的dr对比图；

图4是本发明中不同拥塞链路比例下clink、scfs、dacli-clink算法的fpr对比图；

图5是本发明中不同拥塞链路数目下clink、scfs、dacli-clink算法的dr对比图；

图6是本发明中不同拥塞链路数目下clink、scfs、dacli-clink算法的fpr对比图；

图7是本发明中判断拥塞链路的伪代码示意图。

具体实施方式

为使本发明的目的、技术方案和优点更加清楚，下面将结合附图对本发明实施方式作进一步地描述。

如图1所示，基于链路拥塞强度分布的网络拥塞链路诊断方法，包括：

s1、对待测ip网络中端到端路径进行t时刻快照，获得t时刻待测ip网络中各端到端路径性能及拓扑结构的探测结果，使用快照获得待测ip网络链路的初始概率，初始概率包括链路的初始拥塞概率以及链路的初始拥塞强度概率；

s2、根据t时刻快照所得到的初始概率使用em算法模型计算得到待测ip网络的先验概率，先验概率包括先验拥塞概率和链路的先验拥塞强度概率，em算法建模流程如图2所示；

s3、建立待测ip网络中各链路的先验拥塞概率求解的线性方程组；

s4、根据步骤s2得到待测ip网络的先验概率，结合步骤s3得到的线性方程组，计算推断待测ip网络中最有可能发生拥塞的链路集合。

步骤s2计算先验概率具体包括：

a、根据所述探测结果，获得某一时间片t对应的链路的拥塞概率向量p和链路的拥塞强度概率分布d^t；

b、根据t时刻得到的初始概率以及ip网络路径状态z，计算对数似然ll(d|z，x)关于x的期望，得到最可能发生的链路状态x；

c、根据链路状态x，重新计算每条链路的拥塞强度概率分布，得到新的d^t+1；

d、当迭代到一定次数，或者算法收敛到一定精度，结束算法，得+到所述待测网络中的先验链路拥塞强度概率分布；否则，回到步骤b。

步骤s3建立所述线性方程组具体包括：

a、在所述待测ip网络中，将各端到端与其路径的各链路之间的关系使用布尔代数模

型表示为：

nc为各端到端路径经过的链路总和，zi为第i条路径的状态变量，xj为第j条链路的状态变量。当zi≥1时，表示第i条路径为拥塞状态，当zi＝0时，表示第i条路径为正常状态；当xj＝1时，表示第j条链路位拥塞状态，当xj＝0时，表示第j条链路为正常状态；dij为路径和链路的关系变量，当dij＝1时，表示第i条路径经过第j条链路，当dij＝0时，表示第i条路径不经过第j条链路；

b、在所述待测ip网络中，将各端到端与其路径的各链路之间的拥塞强度的关系使用线性代数模型表示为：

nc为各端到端路径经过的链路总和，mi为第i条路径的状态变量，yj为第j条链路的状态变量。当mi≥1时，表示第i条路径为拥塞状态，当mi＝0时，表示第i条路径为正常状态；当yj≥1时，表示第j条链路位拥塞状态，当yj＝0时，表示第j条链路为正常状态；dij为路径和链路的关系变量，当dij＝1时，表示第i条路径经过第j条链路，当dij＝0时，表示第i条路径不经过第j条链路；

c、对所述a步骤中的布尔代数模型，找到向量x＝[x1x2…xn]，使得如下条件概率公式得到最大化：

其中z＝[z1z2…zm]表示快照得到的待测ip网络的状态，n表示待测ip网络中链路的个数，m表示待测ip网络中路径的个数；

根据步骤d-h得到p：

d、对所述c步骤中条件概率公式，根据贝叶斯方法，变换得到：

e、对所述d步骤中的条件概率公式，可以得到如下的等效最大化问题：

f、对所述e步骤中的等效最大化问题公式，因为链路状态xk是独立于随机变量，所以得到如下的公式：

g、根据所述a步骤中的布尔代数模型，得到如下的公式：

h、根据所述g步骤中的公式，得到如下的公式：

其中，pg为测量标记为正常链路的集合，pc为测量标记为拥塞链路的集合，pi表示路径i。

所述步骤c中最大化概率公式可以化简为：

其中，εr表示去除所有经过正常路径的链路以及确定为正常链路之后所得到的路由矩阵，χ表示判断的拥塞链路集合，pk表示链路k的拥塞概率。

在当前推断时刻拥塞路径链路集合中查找权值最高的链路的判断公式如下：

其中，dk表示链路k经过的路径的集合，|dk|即为链路k经过的路径的个数。

本实施例验证本发明的优越性的具体过程如下：

引入准确率dr、误诊率fpr，dr和fpr的表达式如下所示：

χ为算法推断出的拥塞链路集合，为推断时刻实际发生拥塞的链路。

(1)不同拥塞链路比例对算法的影响

为了验证ip网络在拥塞链路比例不同的情况下对算法性能带来的影响。我们引入拥塞链路比例参数，进行dacli-clink、scfs、clink算法之间性能的比较，模拟实验中设置的拥塞链路比例参数在10％到60％之间，得到dacli-clink、scfs、clink三种算法的dr和fpr，实验结果如图3、图4所示。

在相同的ip网路环境下，随着拥塞链路的比例的逐渐提高，dacli-clink、scfs、clink三种算法的dr均呈下降趋势，fpr均呈现上升趋势，dacli-clink算法的dr和fpr明显高于scfs和clink两种算法。当拥塞链路比例超过50％时，scfs和clink算法的dr均低于0.1，dr随着拥塞链路比例大幅下降，dacli-clink算法的dr在拥塞链路比例超过50％的时候，仍保持在0.7左右。

(2)不同拥塞链路数目对算法的影响

为了验证ip网络在拥塞链路数目不同的情况下对算法性能带来的影响，进行dacli-clink、scfs、clink算法之间性能的比较，模拟试验中设置拥塞链路数目在1到15条之间，得到dacli-clink、scfs、clink三种算法的dr和fpr，实验结果如图5、图6所示。

在相同的ip网络环境下，随着拥塞链路数目的逐渐增加，在三种算法的fpr并无较大差距的情况下，dacli-clink算法的推断准确率均优于scfs算法和clink算法。

通过对不同类型及规模ip网络进行模拟实验和仿真实验，在相同场景下，分别利用传统clink和scfs算法及本发明的方法进行拥塞链路推断，实验证明了本发明的方法较clink算法和scfs算法有更高的推断准确性。

该算法伪代码如图7所示，伪代码中各符号表示为：

mk：端到端测量结果集合；

∈a：未判断的链路集合；

domain(ek)：链路k经过的路径数；

ρc：拥塞路径集合；

yk：链路的拥塞强度；

χ：算法推断出的拥塞链路集合。

上面结合附图对本发明的实施例进行了描述，但是本发明并不局限于上述的具体实施方式，上述的具体实施方式仅仅是示意性的，而不是限制性的，本领域的普通技术人员在本发明的启示下，在不脱离本发明宗旨和权利要求所保护的范围情况下，还可做出很多形式，这些均属于本发明的保护之内。