[发明专利]面向负载均衡的IP网络弹性路由层优化方法

权利要求书

1.一种面向负载均衡的IP网络弹性路由层优化方法，其特征在于包括以下步骤：

步骤一、建立层次化RRL技术体系，最底层为RRL技术平台；第二层为路由子层生成算法，用于得到确定的路由子层分配方案；第三层为转发流量的分发方式；第四层为RRL的应用模式；

步骤二、RRL定义的矩阵表示；

网络IP层拓扑是由节点和链路构成的无向图，记为G(V,E)，其中，V表示顶点的集合，E表示边的集合；若图的顶点个数为n，则表示为一个n×n的矩阵，用D=(d_ij)_n×n表示，其中

式中，(i,j)表示连接顶点i和j的一条边；

对于带权值的图

式中，w(i,j)表示链路(i,j)的权值，矩阵D沿主对角线对称；

定义1：若拓扑矩阵D^F与D¹、D²、…、D^l满足如下关系，

DF=1l-1Σk=1lDk---(3)]]>

且满足矩阵中的元素全部为非零或非∞元素，该条件表示各拓扑子层必须是连通的，n为矩阵的阶数，则矩阵D¹、D²、…、D^l所表示的网络拓扑组合为矩阵D^F所表示的网络拓扑生成的一组RRL；其中，l为拓扑子层的层数；

定义2：最短路径矩阵P=(p_ij)_n×n表示一个n阶图中各顶点间的最短路径，对于无权图，元素p_ij为顶点i与顶点j之间的最少边数，对于有权图，元素p_ij为顶点i与顶点j之间最短通路的权值和；

步骤三、面向负载均衡的IP网络RRL优化问题描述；

（1）给定已知常量；

a.网络拓扑信息，包括：节点数n=|V|，链路数|E|，拓扑矩阵

b.流量需求矩阵D^T；

c.链路容量矩阵C；

d.故障状态集合F，无故障状态表示为f⁰，第e个链路故障表示为f^e；

e.平均最短重路由路径及负载均衡调节权重ω_sp、ω_lb；

f.链路重要度权重表示网络中第e个链路的重要度；

由已知常量，根据Dijkstra算法计算得到以下量，将用于目标函数的计算中：

a.变量正常状态f⁰时，链路(i,j)承载节点I到节点J的流量，则否则(i,j)∈E，I、J∈V；

b.变量从节点I到节点J的无故障传输路径上，若包含故障状态f^e的故障链路，则从节点I到故障链路上游节点I^e的路径上，链路(i,j)包含其中，则否则，Xij,IIe0=0;]]>

c.变量从节点I到节点J的无故障传输路径上，若包含故障状态f^e的故障链路，该故障链路下游节点J_e到节点J的路径上，链路(i,j)包含其中，则否则，Xij,JeJ0=0;]]>

（2）给定决策变量；

a.拓扑子层层数l；

b.拓扑子层矩阵D^k(F_p)，1≤k≤l，F_p表示该拓扑子层保护的故障状态集合；

c.变量故障状态f^e时，在保护该故障状态的弹性路由拓扑子层上，链路(i,j)承载节点I到节点J的重路由流量，则否则

由决策变量得到：

d.故障状态为f^e时，平均重路由路径增加值为

Δse=ΣI=1nΣJ=1npgt(Σi=1nΣj=1n(Xij,IJe·dijk(Fp)-(Xij,IJ0-Xij,IIe0)·dijF))n×n---(4)]]>

其中，函数pgt定义如下

e.故障状态为f^e时，最大链路利用率为

ηe=maxi,j(ΣI=1nΣJ=1n(Xij,IJ0+Xij,IJe-Xij,JeJ0)·dIJTcij)---(6)]]>

（3）建立面向负载均衡的RRL生成优化模型

优化问题描述如下：

minωsp·Σe=1|E|ωle·Δse+ρ·ωlb·Σe=1|E|ωle·ηe---(7)]]>

s.t.  2≤l≤ξ·|E|，0≤ξ≤1     (8)

DF=1l-1Σk=1lDk---(9)]]>

Σi=1n-1(Dk)(i)=(dijk)n×n∩Πi=1nΠj=1ndijk≠0---(10)]]>

DF-Dk=(dij(F-k))n×n∩Σi=1nΣj=1ndij(F-k)≤ceil(Σi=1nΣj=1ndijFl)(k=1,2,···,l)---(11)]]>

Xij,IJe≤dijk(Fp),∀fe∈Fp---(12)]]>

XIi1,IJe·Xi1i2,IJe·Xi2i3,IJe···XinIJJ,IJe+(1-XIi1IJe)·(1-Xi1i2IJe)·(1-Xi2i3,IJe)···(1-XinIJJ,IJe)=1---(13)]]>

Σi=1nΣj=1nXij,IJe·dijk(Fp)≤Σi=1nΣj=1n(Xij,IXe·dijk(Fp)+Xij,XJe·dijk(Fp)),∀fe∈Fp,X∈V---(14)]]>

ΣI=1nΣJ=1n(Xij,IJ0+Xij,IJe-Xij,JeJ0)·dIJT≤cij---(15)]]>

该优化问题中，目标函数为最小化平均重路由路径增加值和最大链路利用率的加权和，当ω_sp较大时，优化目标更侧重于最短重路由路径，保证重路由的传输效率；当ω_lb较大时，优化目标更侧重于链路传输的负载均衡问题，以有效避免链路拥塞；其中，ρ为调节因子，尽可能使重路由路径增加值和最大链路利用率在同一个数量级上，以保证权重因子的调节作用；公式（8）约束了弹性路由拓扑子层的层数l，l最小取值为2，最大取值由ξ进行调节，使得最大取值介于2和链路数|E|之间，ξ根据实际网络拓扑情况决定；公式（9）和公式（10）根据定义1给出，确保了生成的路由子层为弹性路由子层；公式（11）的给出使得每一个弹性路由子层保护的链路数相当，有利于生成恢复路径短的RRL结果；

公式（12）-（15）对网络故障状态下的数据进行了约束；公式（12）限定了的取值，对于任意的故障状态f^e，只有对该故障状态保护的拓扑子层中，节点i与j之间存在链路时的值有可能取值为1，当节点i与j之间不存在链路时只取值为0；公式（13）进一步限定了的取值，表示对于任意的故障状态f^e∈F_p，在拓扑子层D^k(F_p)上，节点I与J之间的路径对应的一系列值同时取1或同时取0，这样限定将使得算法得到的节点I与节点J之间的链路是连通的；公式（14）表示，对于任意故障状态f^e∈F_p，节点I与J在拓扑子层D^k(F_p)上的重路由路径是最短路径；公式（15）表示在任意故障状态情况下，链路(i,j)上承载的流量均小于该链路的容量；

步骤四、单亲遗传算法求解RRL生成优化模型的具体步骤；

Step1:初始编码；令进化代数g=0，按照编码规则给出初始群体，群体中个体数目为n_o；具体的编码规则为：将初始拓扑中的每条链路用符号①、②、③、...进行编号，不同的RRL生成结果看成是这些编码的不同分组；

Step2:适应度计算；由公式（16）计算每个个体的适应度

f=1ωsp·Σe=1|E|ωle·Δse+ρ·ωlb·Σe=1|E|ωle·ηe---(16)]]>

Step3:单亲繁殖；采用单亲遗传方式完成个体的繁殖操作，采用单点基因换位和两点基因换位两种算子；

Step4:淘汰个体；将得到的每一个个体由公式（10）进行连通性验证，若不满足连通性验证则直接淘汰，对于每一个个体，根据公式（12）-（14）的法则计算函数的值，再根据公式（15）判断是否会避免拥塞的发生，若不满足则直接淘汰；重复繁殖过程，直到繁殖的个数达到子代个体数目总值n_p；

Step5:竞争选择；采用选择算子选择出新一代群体，采取父子竞争选择模式，经过家庭竞争和社会竞争两轮竞争来完成选择；竞争选择完成后包含个体数目为n_o；令g=g+1；

Step6:结束判断；如果终止条件满足g大于终止代数G，则算法结束；否则，转到Step2。