CN107317649A - 基于wson损伤模型的最短路径的计算方法及其系统 - Google Patents

基于wson损伤模型的最短路径的计算方法及其系统 Download PDF

Info

Publication number
CN107317649A
CN107317649A CN201710440754.2A CN201710440754A CN107317649A CN 107317649 A CN107317649 A CN 107317649A CN 201710440754 A CN201710440754 A CN 201710440754A CN 107317649 A CN107317649 A CN 107317649A
Authority
CN
China
Prior art keywords
matrix
mrow
path
msub
section
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.)
Granted
Application number
CN201710440754.2A
Other languages
English (en)
Other versions
CN107317649B (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.)
Fiberhome Telecommunication Technologies Co Ltd
Original Assignee
Fiberhome Telecommunication Technologies Co Ltd
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 Fiberhome Telecommunication Technologies Co Ltd filed Critical Fiberhome Telecommunication Technologies Co Ltd
Priority to CN201710440754.2A priority Critical patent/CN107317649B/zh
Publication of CN107317649A publication Critical patent/CN107317649A/zh
Priority to PCT/CN2017/111568 priority patent/WO2018227877A1/zh
Application granted granted Critical
Publication of CN107317649B publication Critical patent/CN107317649B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04J—MULTIPLEX COMMUNICATION
    • H04J14/00—Optical multiplex systems
    • H04J14/02—Wavelength-division multiplex systems
    • H04J14/0227—Operation, administration, maintenance or provisioning [OAMP] of WDM networks, e.g. media access, routing or wavelength allocation
    • H04J14/0254—Optical medium access
    • H04J14/0267—Optical signaling or routing
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04J—MULTIPLEX COMMUNICATION
    • H04J14/00—Optical multiplex systems
    • H04J14/02—Wavelength-division multiplex systems
    • H04J14/0227—Operation, administration, maintenance or provisioning [OAMP] of WDM networks, e.g. media access, routing or wavelength allocation
    • H04J14/0254—Optical medium access
    • H04J14/0267—Optical signaling or routing
    • H04J14/0271—Impairment aware routing
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04Q—SELECTING
    • H04Q11/00—Selecting arrangements for multiplex systems
    • H04Q11/0001—Selecting arrangements for multiplex systems using optical switching
    • H04Q11/0005—Switch and router aspects
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04W—WIRELESS COMMUNICATION NETWORKS
    • H04W40/00—Communication routing or communication path finding
    • H04W40/02—Communication route or path selection, e.g. power-based or shortest path routing
    • H04W40/18—Communication route or path selection, e.g. power-based or shortest path routing based on predicted events
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04W—WIRELESS COMMUNICATION NETWORKS
    • H04W40/00—Communication routing or communication path finding
    • H04W40/02—Communication route or path selection, e.g. power-based or shortest path routing
    • H04W40/22—Communication route or path selection, e.g. power-based or shortest path routing using selective relaying for reaching a BTS [Base Transceiver Station] or an access point
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04Q—SELECTING
    • H04Q11/00—Selecting arrangements for multiplex systems
    • H04Q11/0001—Selecting arrangements for multiplex systems using optical switching
    • H04Q11/0005—Switch and router aspects
    • H04Q2011/0007—Construction
    • H04Q2011/0016—Construction using wavelength multiplexing or demultiplexing
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02D—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
    • Y02D30/00—Reducing energy consumption in communication networks
    • Y02D30/70—Reducing energy consumption in communication networks in wireless communication networks

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Data Exchanges In Wide-Area Networks (AREA)
  • Optical Communication System (AREA)

Abstract

本发明公开了一种基于WSON损伤模型的最短路径的计算方法及其系统，涉及数据和IP传输设备的技术领域。该方法包括如下步骤：根据给定的路由策略计算前K条最短路径；构造矩阵matrix[m][n]；计算数组Xsum[n]和数组rvs_Xsum[n]；构造三个迭代矩阵：矩阵trans[m][n]、矩阵sum[m][n]、矩阵rvs_sum[m][n]；从矩阵trans[m][n]的最后一列选择最小的数值，比较计算得到的K条路径的最少中继数，选择中继数最少的路径作为最短路径。该系统包括前K条最短路径计算模块、最少中继数计算模块和最短路径选择模块。本发明能够计算指定源节点到指定宿节点的最短路径，具有准确率高、实用性强的优点。

Description

基于WSON损伤模型的最短路径的计算方法及其系统
技术领域
本发明涉及数据和IP传输设备的技术领域，具体是涉及一种基于WSON损伤模型的最短路径的计算方法及其系统。
背景技术
WSON(Wavelength Switched Optical Network)也就是基于波分复用传输网的ASON，是目前IETF标准组织倡导的智能波分标准，除了传统ASON(Automatically SwitchedOptical Network)的功能外，主要解决波分网络中光纤/波长自动发现、在线波长路由选择、基于损伤模型的路由选择等问题，其核心技术依然集中在路径计算所采用的算法。
WSON是将控制平面引入到波长网络中，在路径计算时，需要考虑波长的动态变化，以及物理链路的OSNR(Optical Signal Noise Ratio)即光信噪比。通常，初始状态下，网络中每段链路可以使用的波长信号数量一样，由于硬件限制，计算得到的逻辑通道路径上各段链路必须使用统一波长的光信号，波长的动态变化，也正是因为网络拓扑中可以同时存在多条逻辑通道，经过不同路径、使用不用波长信号的通道占用了网络资源，使网络中各链路的可用波长信号发生差异。较为简单的方法是先使用常用的最短路径策略，对给定波长计算一条每段链路该波长都可用的最短路径，若算路失败，遍历每个波长，直到某个波长算路成功或遍历完所有波长。此方法实现简单，但受限于波长不可改变，造成了资源浪费。对此，可以将部分节点作为中继节点，路由经过中继节点时，允许改变后续链路传输信号的波长，同时，也允许路径经过该节点但不作为中继节点(即波长不可变)。
通常一条路径的OSNR可以由每段链路的硬件参数计算，且随着链路的增加而降低，一条符合传输要求的路径，其总OSNR必须大于给定阈值。中继节点的另一个功能就是对光信号重新放大，使整条路径的OSNR只需要满足：源节点到下一个中继节点、相邻中继节点、宿节点到上一个中继节点之间的OSNR都大于阈值。
所以，中继功能可以分为两个部分：允许变换波长、重置OSNR。如图1所示，要建立节点1到节点3的通道，由于链路1-2只有λ1波可用，链路2-3只有λ2波可用，则必须使节点2作为中继节点，才能允许通道路径的子路段信号波长不一致。又如图2所示，假设路径上波长都可用，链路1-2、2-3的OSNR都是20，而路径1-2-3的整体OSNR为15，若给定阈值为18，则必须使节点2作为中继节点，使路径1-2-3的两个子路段OSNR都大于阈值。
传统意义上的最短路径算法只需要满足节点最少、链路代价最小等路由策略，结合中继功能重新定义后的最短路径，若将路径中源节点到下一个中继节点、相邻中继节点之间、宿节点到上一个中继节点之间路径简称为子路段，其算法需求可以总结为：寻找一条指定源节点到宿节点的最短路径，要求路径中的每一条子路段选择一个都可用的统一波长，并且所有子路段OSNR不超过给定阈值。最后，由于实现中继功能的硬件代价关系，当路径经过中继节点时，该节点也可以选择不使用中继功能，作为普通节点，于是，新的算法还必须满足一个最重要的限制条件：要求路径上作为中继的节点数最少，并且此条件要更优先于路由策略。如图3所示，节点1到节点3的最短路径在节点最少的路由策略下应该是1-2-3，且节点2作为中继，但中继最少的需求更加优先时，1-4-5-6-3才是最佳路径。
发明内容
本发明的目的是为了克服上述背景技术的不足，提供一种基于WSON损伤模型的最短路径的计算方法及其系统，能够计算指定源节点到指定宿节点的最短路径，具有准确率高、实用性强的优点。
本发明提供一种基于WSON损伤模型的最短路径的计算方法，包括如下步骤：
A、根据给定的路由策略计算前K条最短路径，K为正整数；
B、计算步骤A中K条路径的最少中继数，每一条路径的最少中继数的计算步骤如下：
B1、以该路径上所有具备中继功能的节点为界，将该路径分割为若干子路段，构造矩阵matrix[m][n]，其中，n为路段数，m为波长数；
B2、将子路段的正向光纤经过的所有光放大器OA的Xi叠加值，作为子路段总光信噪比OSNR计算指标，得到数组Xsum[n]；将子路段的反向光纤经过的所有光放大器OA的Xi叠加值，作为子路段的反向OSNR计算指标，得到数组rvs_Xsum[n]；其中，Xi为OSNR的计算公式的指数部分，i为正整数，i的取值范围为1～N，N为路径单方向所经过的OA的总数，每一条子路段OSNR都大于给定阈值T；
B3、根据动态规划算法思想，输入参数matrix[m][n]、Xsum[n]、rvs_Xsum[n]，构造三个迭代矩阵：矩阵trans[m][n]、矩阵sum[m][n]、矩阵rvs_sum[m][n]；矩阵trans[m][n]随着算法迭代过程记录当前最少中继数；矩阵sum[m][n]和矩阵rvs_sum[m][n]记录当前中继最少方案下到前一个中继或源节点的OSNR计算指标叠加量；
B4、从矩阵trans[m][n]的最后一列选择最小的数值，即为整个路径最少中继数，按迭代规则回溯矩阵，得到路径的波长选择结果和中继选择结果；若矩阵trans[m][n]最后一列都是无效值，则该路径不通，动态规划失败；
C、比较步骤B中计算得到的K条路径的最少中继数，选择中继数最少的路径作为最短路径。
在上述技术方案的基础上，步骤B1中，所述矩阵matrix[m][n]的构造过程如下：以子路段按顺序映射为行，以路段对所有波长是否可用映射为列，且保证路径中的每一条子路段至少存在一个都可用的统一波长，矩阵matrix[m][n]的矩阵取值为0或1；
其中，0表示第n路段上第m波长不可用，1表示第n路段上第m波长可用。
在上述技术方案的基础上，步骤B2中，所述OSNR的计算公式为：
其中，opi表示第i个OA的输入光功率，Mi为第i个OA所在链路背景波道数，nfi为第i个OA的噪声系数；
所述Xi的计算公式为：
所述Xi叠加值满足其中，t＝10-T/10。
在上述技术方案的基础上，步骤B3中，所述矩阵trans[m][n]的构造过程如下：设置矩阵trans[m][n]内的数值初始化均为无效值，对矩阵trans[m][n]进行逐列算法迭代，每列按照从上至下的顺序计算，共计算m*n次。
在上述技术方案的基础上，步骤C中，选择中继数最少的路径作为最短路径时，若出现多条路径的中继数相同，则选择路由策略最佳的路径，作为最短路径。
本发明还提供一种基于WSON损伤模型的最短路径的计算系统，该系统包括前K条最短路径计算模块、最少中继数计算模块和最短路径选择模块；
所述前K条最短路径计算模块用于：根据给定的路由策略计算前K条最短路径，K为正整数；
所述最少中继数计算模块用于：计算K条路径的最少中继数，每一条路径的最少中继数的计算步骤如下：
以该路径上所有具备中继功能的节点为界，将该路径分割为若干子路段，构造矩阵matrix[m][n]，其中，n为路段数，m为波长数；
将子路段的正向光纤经过的所有光放大器OA的Xi叠加值作为子路段总光信噪比OSNR计算指标，得到数组Xsum[n]；将子路段的反向光纤经过的所有光放大器OA的Xi叠加值作为子路段的反向OSNR计算指标，得到数组rvs_Xsum[n]；其中，Xi为OSNR的计算公式的指数部分，i为正整数，i的取值范围为1～N，N为路径单方向所经过的OA的总数，每一条子路段OSNR都大于给定阈值T；
根据动态规划算法思想，输入参数matrix[m][n]、Xsum[n]、rvs_Xsum[n]，构造三个迭代矩阵：矩阵trans[m][n]、矩阵sum[m][n]、矩阵rvs_sum[m][n]；矩阵trans[m][n]随着算法迭代过程记录当前最少中继数；矩阵sum[m][n]和矩阵rvs_sum[m][n]记录当前中继最少方案下到前一个中继或源节点的OSNR计算指标叠加量；
从矩阵trans[m][n]的最后一列选择最小的数值，即为整个路径最少中继数，按迭代规则回溯矩阵，得到路径的波长选择结果和中继选择结果；若矩阵trans[m][n]最后一列都是无效值，则该路径不通，动态规划失败；
所述最短路径选择模块用于：比较最少中继数计算模块计算得到的K条路径的最少中继数，选择中继数最少的路径作为最短路径。
在上述技术方案的基础上，所述最少中继数计算模块具体用于：以子路段按顺序映射为行，以路段对所有波长是否可用映射为列，构造矩阵matrix[m][n]，且保证路径中的每一条子路段至少存在一个都可用的统一波长，将矩阵matrix[m][n]取值0或1；
其中，0表示第n路段上第m波长不可用，1表示第n路段上第m波长可用。
在上述技术方案的基础上，所述最少中继数计算模块具体用于：
计算OSNR的公式为：
其中，opi表示第i个OA的输入光功率，Mi为第i个OA所在链路背景波道数，nfi为第i个OA的噪声系数；
计算Xi的计算公式为：
保证Xi叠加值满足其中，t＝10-T/10。
在上述技术方案的基础上，所述最少中继数计算模块还用于：设置矩阵trans[m][n]内的数值初始化均为无效值，对矩阵trans[m][n]进行逐列算法迭代，每列按照从上至下的顺序计算，共计算m*n次。
在上述技术方案的基础上，所述最短路径选择模块还用于：选择中继数最少的路径作为最短路径时，若出现多条路径的中继数相同，最短路径选择模块选择路由策略最佳的路径，作为最短路径。
与现有技术相比，本发明的优点如下：本发明结合中继功能，提出一种针对WSON网络波长动态变化，以及OSNR损伤模型的最短路径算法，能够自动寻找一条指定源节点到指定宿节点的最短路径。本发明先根据给定的源节点、宿节点以及路由策略计算出前K条最短路径，使用动态规划(Dynamic Programming)算法思想迭代求出该条路径中哪个节点作为中继，以及其每条子路段选用哪个波长；最后，选择K条路径中使用中继最少的结果作为算法输出的最短路径。本发明的算法不仅简单，而且准确率高、实用性强。
附图说明
图1是现有的三节点的网络拓扑结构示意图，其中λ1、λ2为信号波长。
图2是现有的三节点的链路中OSNR的网络拓扑结构示意图。
图3是现有的六节点的网络拓扑结构示意图。
图4是本发明实施例基于WSON损伤模型的最短路径的计算方法的流程图。
图5是本发明实施例计算K条路径的最少中继数的流程图。
图6是本发明实施例的七节点的网络拓扑结构示意图。
图7是本发明实施例的两节点的光放大器OA的正向光纤和反向光纤的网络拓扑结构示意图。
具体实施方式
下面结合附图及具体实施例对本发明作进一步的详细描述。
参见图4所示，本发明实施例提供一种基于WSON损伤模型的最短路径的计算方法，包括如下步骤：
S1、根据给定的路由策略计算前K条最短路径，K为正整数；
S2、计算步骤S1中K条路径的最少中继数，参见图5所示，每一条路径的最少中继数的计算步骤如下：
S201、以该路径上所有具备中继功能的节点为界，将该路径分割为若干子路段，构造矩阵matrix[m][n]，其中，n为路段数，m为波长数；
S202、将子路段的正向光纤经过的所有光放大器OA的Xi叠加值作为子路段总光信噪比OSNR计算指标，得到数组Xsum[n]；将子路段的反向光纤经过的所有光放大器OA的Xi叠加值作为子路段的反向OSNR计算指标，得到数组rvs_Xsum[n]；其中，Xi为OSNR的计算公式的指数部分，i为正整数，i的取值范围为1～N，N为路径单方向所经过的OA的总数，每一条子路段OSNR都大于给定阈值T；
S203、根据动态规划算法思想，输入参数matrix[m][n]、Xsum[n]、rvs_Xsum[n]，构造三个迭代矩阵：矩阵trans[m][n]、矩阵sum[m][n]、矩阵rvs_sum[m][n]；矩阵trans[m][n]随着算法迭代过程记录当前最少中继数；矩阵sum[m][n]和矩阵rvs_sum[m][n]记录当前中继最少方案下到前一个中继或源节点的OSNR计算指标叠加量；
S204、从矩阵trans[m][n]的最后一列选择最小的数值，即为整个路径最少中继数，按迭代规则回溯矩阵，得到路径的波长选择结果和中继选择结果；若矩阵trans[m][n]最后一列都是无效值，则该路径不通，动态规划失败；
S3、比较步骤S2中计算得到的K条路径的最少中继数，选择中继数最少的路径作为最短路径。其中，步骤S1中，已有许多完善算法计算前K条最短路径，如KSP(Top K ShortestPaths)算法，K的大小可以根据实际情况调整，计算出最短路径，第二短的路径，第三短的路径，一直到第K短的路径，一共K条。先按路由策略选出前K短的K条路径，再在这K条路径中，根据波长分配和中继情况，选出一条最符合要求的路径。
其中，步骤S201中，矩阵matrix[m][n]的构造过程如下：以子路段按顺序映射为行，以路段对所有波长是否可用映射为列，且保证路径中的每一条子路段至少存在一个都可用的统一波长，矩阵matrix[m][n]的矩阵取值为0或1，其中，0表示第n路段上第m波长不可用，1表示第n路段上第m波长可用。参见图6所示，图6省略OA器件，只显示节点与链路，若其中一条路径为1-2-3-4-5-6-7，其中节点3和节点5具备中继功能，则将路径分割为1-2-3、3-4-5、5-6-7三个子路段，图中链路标注了可用波长，假如所有可用波为λ1-λ5，则分别取三段子路段的共公可用波集合，分别为λ1-λ4，λ3、λ4，λ4，以子路段按顺序映射为行，路段对所有波长是否可用映射为列(此时路段数n＝3，波长数m＝5)，构造矩阵matrix[m][n]，取值0或1，分别表示第n路段上第m波长是否可用，图6的场景可以构造如下矩阵；
步骤S202中，OSNR的计算公式为：
其中，opi表示第i个OA的输入光功率，Mi为第i个OA所在链路背景波道数，nfi为第i个OA的噪声系数；
Xi的计算公式为：
子路段单向光纤经过的所有光放大器OA的Xi叠加值，满足
其中，t＝10-T/10。
根据公式(1)、(2)可以得到：
T为给定阈值，每一条子路段OSNR都应大于给定阈值T。
由于Xi为OSNR的计算公式的指数部分，无论op、nf、M如何取值，Xi必定大于0；则无论路径怎样取，随着路径经过的OA数量N的增加而单调递增，且Xi越小越好，根据公式(1)、(2)、(4)，令t＝10-T/10，得出公式(3)。
由于每条路径可以根据中继节点分割为若干条子路段，则每一条子路段OSNR都应大于给定阈值T可以转换为：若子路段经过N个OA，则子路段必须满足公式(3)。
参见图7所示，节点1与节点2之间的子路段收发光纤各经过2个OA，其参数Xn根据公式(2)计算得到X1、X2、X3、X4，则该路段Xsum＝X1+X2，rvs_Xsum＝X3+X4。实际中，一般实际场景中路径正向传播OSNR与反向传播的OSNR会有差异。
其中，步骤S203中，设置矩阵trans[m][n]内的数值初始化均为无效值，对矩阵trans[m][n]进行逐列算法迭代，每列按照从上至下的顺序计算，共计算m*n次。在实际应用中，无效值可以根据情况指定，只要编码的时候能够判断出哪个是无效值即可。由于子路段数每增加1，中继数最多增加1，则使子路段数为n的路径中继数最少，则其子路段数为n-1的路径中继数必然也要最少，又因为当OSNR超限或者下一段子路段同波长不可用时中继数必加1，每次迭代时若第j(0<＝j<n)路段上第i(0<＝i<m)波可用，则重新计算并将trans[i][j]替换为有效值，每次迭代结果如下(0<＝i<m，0<＝j<n)：
if(sum[i][j-1]+Xsum[j]<t)&&(rvs_sum[i][j-1]+rvs_Xsum[j]<t)&&
(matrix[i][j-1]＝＝1)&&(trans[i][j-1]<(min(trans[0][j-1]to trans[m-1][j-1])+1))
trans[i][j]＝trans[i][j-1]
sum[i][j]＝sum[i][j-1]+Xsum[j]
rvssum[i][j]＝rvssum[i][j-1]+rvs_Xsum[j]
else
trans[i][j]＝min(min(trans[0][j-1]to trans[m-1][j-1])+1,trans[i][j-1])
sum[i][j]＝Xsum[j]
rvssum[i][j]＝rvs_Xsum[j]
迭代完成后，trans矩阵保存了选用每个波长在当前子路段所需要的中继次数，即隐含了该条路径中继最少的波长、中继分配方案。若某一拓扑由步骤S201得到的matrix如下：
由步骤S202得到的Xsum与rvs_Xsum都为[1，1，1]，t＝2，则完成迭代后得到的trans矩阵为：
其中“-”表示初始化的无效值则算法的9次迭代过程如下：
trans[0][0]＝0表示若子路段1分配λ1，由于为起始路段，当前需要中继次数为0；
trans[1][0]＝-表示子路段1λ2不可用，保持初始无效值；
trans[2][0]＝0同trans[0][0]；
trans[0][1]＝-同trans[1][0]；
trans[1][1]＝1表示若子路段2分配λ2，必须进行波长转换，选择前一列最小的值加1，即trans[1][1]＝trans[0][0]+1；
trans[2][1]＝0表示若子路段2分配λ3，不需要波长转换，直接使trans[2][1]＝trans[2][0]；
trans[0][2]＝-同trans[1][0]；
trans[1][2]＝1表示若子路段3分配λ2，不需要波长转换，直接使trans[1][2]＝trans[1][1]；
trans[2][2]＝1表示若子路段3分配λ3，不需要进行波长转换，但由于Xsum累加值为3>t，必须进行中继放大，选择前一列最小的值加1，即trans[2][2]＝trans[2][1]+1。
在最后一列最小值为1，可以选择给子路段3分配λ2，子路段2分配λ2，子路段1分配λ1，并且一旦发生状态转移(即矩阵元素值发生加1)，则将发生状态转移的两个子路段之间的节点作为中继节点，即子路段1与子路段2之间的节点作为中继；若子路段3选择λ3(需要的中继数一样，可以根据不同需求选择)，则子路段1与子路段2都可以分配λ3，但在子路段2与3之间发生状态转移，将子路段2与子路段3之间的节点作为中继。
在实际应用中，步骤S204中，一条路径的矩阵trans[m][n]中，若最后一列存在多个相同的最小中继数，可以根据每次算路给定的指定波长，最终选择与该指定波长较为接近的波长。
步骤S3中，选择中继数最少的路径作为最短路径时，若出现多条路径的中继数相同，则选择路由策略最佳的路径，作为最短路径。
此外，设网络拓扑中共V个节点、E条边，初始状态可用波长数量为N，一般计算单条最短路径的时间复杂度大约为O(E+VlgV)，空间复杂度为O(V)，对每条路径执行步骤S203算法时，最坏的时间复杂度为O(NNV)，空间复杂度为O(NV)，则整个算法时间复杂度为O(K(E+VlgV+NNV))，空间复杂度为O(K(V+NV))。所以，使用本方法计算路径总共时间、空间复杂度都在多项式级别内，对拓扑规模没有要求。
本发明还公开一种基于WSON损伤模型的最短路径的计算系统，该系统包括前K条最短路径计算模块、最少中继数计算模块和最短路径选择模块；
前K条最短路径计算模块用于：根据给定的路由策略计算前K条最短路径，K为正整数；
最少中继数计算模块用于：计算K条路径的最少中继数，每一条路径的最少中继数的计算步骤如下：
以该路径上所有具备中继功能的节点为界，将该路径分割为若干子路段，构造矩阵matrix[m][n]，其中，n为路段数，m为波长数；
将子路段的正向光纤经过的所有光放大器OA的Xi叠加值作为子路段总光信噪比OSNR计算指标，得到数组Xsum[n]；将子路段的反向光纤经过的所有光放大器OA的Xi叠加值作为子路段的反向OSNR计算指标，得到数组rvs_Xsum[n]；其中，Xi为OSNR的计算公式的指数部分，i为正整数，i的取值范围为1～N，N为路径单方向所经过的OA的总数，每一条子路段OSNR都大于给定阈值T；
根据动态规划算法思想，输入参数matrix[m][n]、Xsum[n]、rvs_Xsum[n]，构造三个迭代矩阵：矩阵trans[m][n]、矩阵sum[m][n]、矩阵rvs_sum[m][n]；矩阵trans[m][n]随着算法迭代过程记录当前最少中继数；矩阵sum[m][n]和矩阵rvs_sum[m][n]记录当前中继最少方案下到前一个中继或源节点的OSNR计算指标叠加量；
从矩阵trans[m][n]的最后一列选择最小的数值，即为整个路径最少中继数，按迭代规则回溯矩阵，得到路径的波长选择结果和中继选择结果；若矩阵trans[m][n]最后一列都是无效值，则该路径不通，动态规划失败；
最短路径选择模块用于：比较最少中继数计算模块计算得到的K条路径的最少中继数，选择中继数最少的路径作为最短路径。
其中，最少中继数计算模块具体用于：以子路段按顺序映射为行，以路段对所有波长是否可用映射为列，构造矩阵matrix[m][n]，且保证路径中的每一条子路段至少存在一个都可用的统一波长，将矩阵matrix[m][n]取值0或1；
其中，0表示第n路段上第m波长不可用，1表示第n路段上第m波长可用。
其中，最少中继数计算模块具体用于：
计算OSNR的公式为：
其中，opi表示第i个OA的输入光功率，Mi为第i个OA所在链路背景波道数，nfi为第i个OA的噪声系数；
计算Xi的计算公式为：
保证子路段单向光纤经过的所有光放大器OA的Xi叠加值满足其中，t＝10-T/10。
其中，最少中继数计算模块还用于：设置矩阵trans[m][n]内的数值初始化均为无效值，对矩阵trans[m][n]进行逐列算法迭代，每列按照从上至下的顺序计算，共计算m*n次。
其中，最短路径选择模块还用于：选择中继数最少的路径作为最短路径时，若出现多条路径的中继数相同，最短路径选择模块选择路由策略最佳的路径，作为最短路径。
本领域的技术人员可以对本发明实施例进行各种修改和变型，倘若这些修改和变型在本发明权利要求及其等同技术的范围之内，则这些修改和变型也在本发明的保护范围之内。

Claims (10)

1.一种基于WSON损伤模型的最短路径的计算方法，其特征在于：包括如下步骤：
A、根据给定的路由策略计算前K条最短路径，K为正整数；
B、计算步骤A中K条路径的最少中继数，每一条路径的最少中继数的计算步骤如下：
B1、以该路径上所有具备中继功能的节点为界，将该路径分割为若干子路段，构造矩阵matrix[m][n]，其中，n为路段数，m为波长数；
B2、将子路段的正向光纤经过的所有光放大器OA的Xi叠加值，作为子路段总光信噪比OSNR计算指标，得到数组Xsum[n]；将子路段的反向光纤经过的所有光放大器OA的Xi叠加值，作为子路段的反向OSNR计算指标，得到数组rvs_Xsum[n]；其中，Xi为OSNR的计算公式的指数部分，i为正整数，i的取值范围为1～N，N为路径单方向所经过的OA的总数，每一条子路段OSNR都大于给定阈值T；
B3、根据动态规划算法思想，输入参数matrix[m][n]、Xsum[n]、rvs_Xsum[n]，构造三个迭代矩阵：矩阵trans[m][n]、矩阵sum[m][n]、矩阵rvs_sum[m][n]；矩阵trans[m][n]随着算法迭代过程记录当前最少中继数；矩阵sum[m][n]和矩阵rvs_sum[m][n]记录当前中继最少方案下到前一个中继或源节点的OSNR计算指标叠加量；
B4、从矩阵trans[m][n]的最后一列选择最小的数值，即为整个路径最少中继数，按迭代规则回溯矩阵，得到路径的波长选择结果和中继选择结果；若矩阵trans[m][n]最后一列都是无效值，则该路径不通，动态规划失败；
C、比较步骤B中计算得到的K条路径的最少中继数，选择中继数最少的路径作为最短路径。
2.如权利要求1所述的基于WSON损伤模型的最短路径的计算方法，其特征在于：步骤B1中，所述矩阵matrix[m][n]的构造过程如下：以子路段按顺序映射为行，以路段对所有波长是否可用映射为列，且保证路径中的每一条子路段至少存在一个都可用的统一波长，矩阵matrix[m][n]的矩阵取值为0或1；
其中，0表示第n路段上第m波长不可用，1表示第n路段上第m波长可用。
3.如权利要求1所述的基于WSON损伤模型的最短路径的计算方法，其特征在于：步骤B2中，所述OSNR的计算公式为：
<mrow> <mi>O</mi> <mi>S</mi> <mi>N</mi> <mi>R</mi> <mo>=</mo> <mo>-</mo> <mn>10</mn> <mi>lg</mi> <mrow> <mo>(</mo> <msubsup> <mi>&amp;Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>N</mi> </msubsup> <msup> <mn>10</mn> <mrow> <mo>-</mo> <mrow> <mo>(</mo> <msub> <mi>op</mi> <mi>i</mi> </msub> <mo>-</mo> <mn>10</mn> <msub> <mi>lgM</mi> <mi>i</mi> </msub> <mo>-</mo> <msub> <mi>nf</mi> <mi>i</mi> </msub> <mo>+</mo> <mn>58</mn> <mo>)</mo> </mrow> <mo>/</mo> <mn>10</mn> </mrow> </msup> <mo>)</mo> </mrow> </mrow>
其中，opi表示第i个OA的输入光功率，Mi为第i个OA所在链路背景波道数，nfi为第i个OA的噪声系数；
所述Xi的计算公式为：
<mrow> <msub> <mi>X</mi> <mi>i</mi> </msub> <mo>=</mo> <msup> <mn>10</mn> <mrow> <mo>-</mo> <mrow> <mo>(</mo> <msub> <mi>op</mi> <mi>i</mi> </msub> <mo>-</mo> <mn>10</mn> <msub> <mi>lgM</mi> <mi>i</mi> </msub> <mo>-</mo> <msub> <mi>nf</mi> <mi>i</mi> </msub> <mo>+</mo> <mn>58</mn> <mo>)</mo> </mrow> <mo>/</mo> <mn>10</mn> </mrow> </msup> <mo>.</mo> </mrow>
所述Xi叠加值满足其中，t＝10-T/10。
4.如权利要求1所述的基于WSON损伤模型的最短路径的计算方法，其特征在于：步骤B3中，所述矩阵trans[m][n]的构造过程如下：设置矩阵trans[m][n]内的数值初始化均为无效值，对矩阵trans[m][n]进行逐列算法迭代，每列按照从上至下的顺序计算，共计算m*n次。
5.如权利要求1所述的基于WSON损伤模型的最短路径的计算方法，其特征在于：步骤C中，选择中继数最少的路径作为最短路径时，若出现多条路径的中继数相同，则选择路由策略最佳的路径，作为最短路径。
6.一种基于WSON损伤模型的最短路径的计算系统，其特征在于：包括前K条最短路径计算模块、最少中继数计算模块和最短路径选择模块；
所述前K条最短路径计算模块用于：根据给定的路由策略计算前K条最短路径，K为正整数；
所述最少中继数计算模块用于：计算K条路径的最少中继数，每一条路径的最少中继数的计算步骤如下：
以该路径上所有具备中继功能的节点为界，将该路径分割为若干子路段，构造矩阵matrix[m][n]，其中，n为路段数，m为波长数；
将子路段的正向光纤经过的所有光放大器OA的Xi叠加值作为子路段总光信噪比OSNR计算指标，得到数组Xsum[n]；将子路段的反向光纤经过的所有光放大器OA的Xi叠加值作为子路段的反向OSNR计算指标，得到数组rvs_Xsum[n]；其中，Xi为OSNR的计算公式的指数部分，i为正整数，i的取值范围为1～N，N为路径单方向所经过的OA的总数，每一条子路段OSNR都大于给定阈值T；
根据动态规划算法思想，输入参数matrix[m][n]、Xsum[n]、rvs_Xsum[n]，构造三个迭代矩阵：矩阵trans[m][n]、矩阵sum[m][n]、矩阵rvs_sum[m][n]；矩阵trans[m][n]随着算法迭代过程记录当前最少中继数；矩阵sum[m][n]和矩阵rvs_sum[m][n]记录当前中继最少方案下到前一个中继或源节点的OSNR计算指标叠加量；
从矩阵trans[m][n]的最后一列选择最小的数值，即为整个路径最少中继数，按迭代规则回溯矩阵，得到路径的波长选择结果和中继选择结果；若矩阵trans[m][n]最后一列都是无效值，则该路径不通，动态规划失败；
所述最短路径选择模块用于：比较最少中继数计算模块计算得到的K条路径的最少中继数，选择中继数最少的路径作为最短路径。
7.根据权利要求6所述的基于WSON损伤模型的最短路径的计算系统，其特征在于：所述最少中继数计算模块具体用于：以子路段按顺序映射为行，以路段对所有波长是否可用映射为列，构造矩阵matrix[m][n]，且保证路径中的每一条子路段至少存在一个都可用的统一波长，将矩阵matrix[m][n]取值0或1；
其中，0表示第n路段上第m波长不可用，1表示第n路段上第m波长可用。
8.根据权利要求6所述的基于WSON损伤模型的最短路径的计算系统，其特征在于：所述最少中继数计算模块具体用于：
计算OSNR的公式为：
<mrow> <mi>O</mi> <mi>S</mi> <mi>N</mi> <mi>R</mi> <mo>=</mo> <mo>-</mo> <mn>10</mn> <mi>lg</mi> <mrow> <mo>(</mo> <msubsup> <mi>&amp;Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>N</mi> </msubsup> <msup> <mn>10</mn> <mrow> <mo>-</mo> <mrow> <mo>(</mo> <msub> <mi>op</mi> <mi>i</mi> </msub> <mo>-</mo> <mn>10</mn> <msub> <mi>lgM</mi> <mi>i</mi> </msub> <mo>-</mo> <msub> <mi>nf</mi> <mi>i</mi> </msub> <mo>+</mo> <mn>58</mn> <mo>)</mo> </mrow> <mo>/</mo> <mn>10</mn> </mrow> </msup> <mo>)</mo> </mrow> </mrow>
其中，opi表示第i个OA的输入光功率，Mi为第i个OA所在链路背景波道数，nfi为第i个OA的噪声系数；
计算Xi的计算公式为：
<mrow> <msub> <mi>X</mi> <mi>i</mi> </msub> <mo>=</mo> <msup> <mn>10</mn> <mrow> <mo>-</mo> <mrow> <mo>(</mo> <msub> <mi>op</mi> <mi>i</mi> </msub> <mo>-</mo> <mn>10</mn> <msub> <mi>lgM</mi> <mi>i</mi> </msub> <mo>-</mo> <msub> <mi>nf</mi> <mi>i</mi> </msub> <mo>+</mo> <mn>58</mn> <mo>)</mo> </mrow> <mo>/</mo> <mn>10</mn> </mrow> </msup> <mo>.</mo> </mrow>
保证Xi叠加值满足 <mrow> <mn>0</mn> <mo>&lt;</mo> <msubsup> <mi>&amp;Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>N</mi> </msubsup> <msub> <mi>X</mi> <mi>i</mi> </msub> <mo>&lt;</mo> <mi>t</mi> <mo>,</mo> </mrow> 其中，t＝10-T/10。
9.根据权利要求6所述的基于WSON损伤模型的最短路径的计算系统，其特征在于：所述最少中继数计算模块还用于：设置矩阵trans[m][n]内的数值初始化均为无效值，对矩阵trans[m][n]进行逐列算法迭代，每列按照从上至下的顺序计算，共计算m*n次。
10.根据权利要求6所述的基于WSON损伤模型的最短路径的计算系统，其特征在于：所述最短路径选择模块还用于：选择中继数最少的路径作为最短路径时，若出现多条路径的中继数相同，最短路径选择模块选择路由策略最佳的路径，作为最短路径。
CN201710440754.2A 2017-06-13 2017-06-13 基于wson损伤模型的最短路径的计算方法及其系统 Active CN107317649B (zh)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN201710440754.2A CN107317649B (zh) 2017-06-13 2017-06-13 基于wson损伤模型的最短路径的计算方法及其系统
PCT/CN2017/111568 WO2018227877A1 (zh) 2017-06-13 2017-11-17 基于wson损伤模型的最短路径的计算方法及其系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710440754.2A CN107317649B (zh) 2017-06-13 2017-06-13 基于wson损伤模型的最短路径的计算方法及其系统

Publications (2)

Publication Number Publication Date
CN107317649A true CN107317649A (zh) 2017-11-03
CN107317649B CN107317649B (zh) 2019-01-18

Family

ID=60183241

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710440754.2A Active CN107317649B (zh) 2017-06-13 2017-06-13 基于wson损伤模型的最短路径的计算方法及其系统

Country Status (2)

Country Link
CN (1) CN107317649B (zh)
WO (1) WO2018227877A1 (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2018227877A1 (zh) * 2017-06-13 2018-12-20 烽火通信科技股份有限公司 基于wson损伤模型的最短路径的计算方法及其系统
CN113726442A (zh) * 2021-09-02 2021-11-30 苏州大学 一种基于全光波长转换的路由和波长分配方法及装置

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7391732B1 (en) * 2002-08-05 2008-06-24 At&T Corp. Scheme for randomized selection of equal cost links during restoration
CN103441929A (zh) * 2013-09-17 2013-12-11 烽火通信科技股份有限公司 一种基于wson网络减少波长连续性限制的方法
CN103812790A (zh) * 2014-01-14 2014-05-21 河北工程大学 一种灵活光网络中调制格式自适应的路由和频谱分配方法
US20150350057A1 (en) * 2014-06-03 2015-12-03 National Cheng Kung University Switchless network topology system for parallel computation and method thereof
CN105634954A (zh) * 2016-01-08 2016-06-01 烽火通信科技股份有限公司 基于wson网络考虑光损伤的最短路径计算方法

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8891382B2 (en) * 2011-07-07 2014-11-18 Futurewei Technologies, Inc. Impairment aware path computation element method and system
US9065565B2 (en) * 2013-02-11 2015-06-23 Cisco Technology, Inc. DWDM fast lightpath setup using network status information
US9166723B2 (en) * 2013-02-14 2015-10-20 Nec Laboratories America, Inc. Network fragmentation measurement in an optical wavelength division multiplexing (WDM) network
CN104486094B (zh) * 2014-12-15 2018-03-23 西安电子科技大学 基于物理拓扑光树修正的多播业务疏导方法
CN107317649B (zh) * 2017-06-13 2019-01-18 烽火通信科技股份有限公司 基于wson损伤模型的最短路径的计算方法及其系统

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7391732B1 (en) * 2002-08-05 2008-06-24 At&T Corp. Scheme for randomized selection of equal cost links during restoration
CN103441929A (zh) * 2013-09-17 2013-12-11 烽火通信科技股份有限公司 一种基于wson网络减少波长连续性限制的方法
CN103812790A (zh) * 2014-01-14 2014-05-21 河北工程大学 一种灵活光网络中调制格式自适应的路由和频谱分配方法
US20150350057A1 (en) * 2014-06-03 2015-12-03 National Cheng Kung University Switchless network topology system for parallel computation and method thereof
CN105634954A (zh) * 2016-01-08 2016-06-01 烽火通信科技股份有限公司 基于wson网络考虑光损伤的最短路径计算方法

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2018227877A1 (zh) * 2017-06-13 2018-12-20 烽火通信科技股份有限公司 基于wson损伤模型的最短路径的计算方法及其系统
CN113726442A (zh) * 2021-09-02 2021-11-30 苏州大学 一种基于全光波长转换的路由和波长分配方法及装置

Also Published As

Publication number Publication date
CN107317649B (zh) 2019-01-18
WO2018227877A1 (zh) 2018-12-20

Similar Documents

Publication Publication Date Title
Ramamurthy et al. Fixed-alternate routing and wavelength conversion in wavelength-routed optical networks
US9686176B2 (en) Constrained and viable path computation systems and methods
CN101095320B (zh) 光通信网络内的光路径路由
CN101883293A (zh) 在邻接节点间多边条件下实现k优路径算法的方法及装置
CN101699803B (zh) 一种基于分布式多层多域光网络的域间路由方法
CN102769563B (zh) 路由计算方法及装置
CN107317649A (zh) 基于wson损伤模型的最短路径的计算方法及其系统
CN107566926A (zh) 一种光传输网络规划中分配中继资源的计算方法
CN102868624B (zh) 光电复合型网络节点的控制装置、控制系统及控制方法
CN105007223A (zh) 一种基于光层次架构的光网络动态多播路由波长分配方法
Katib et al. A study on layer correlation effects through a multilayer network optimization problem
Morea et al. A critical analysis of the possible cost savings of translucent networks
CN102934491A (zh) 波分网络规划方法及设备
KR102614809B1 (ko) 네트워크 제어 방법, 장치 및 시스템
CN102457782A (zh) 一种用于波分复用光网络的路由选择鲁棒选路算法
CN105611435A (zh) 基于wson网络的跨光层的电层子波长路由的计算方法
CN101383759A (zh) 一种光网络中划分管理区域的保护方法
Asadi et al. Photonic network-on-chip: A survey
Renani et al. NLR-OP: a high-performance optical router based on North-Last turning model for multicore processors
CN108696329A (zh) 基于二维Torus架构的大规模光网络拓扑设计方法
Sangeetha et al. Wavelength assignment problem in optical WDM networks
CN114885236A (zh) 一种光网络路径规划方法和装置
Asthana Study of p-cycle based protection in optical networks and removal of its shortcomings
Wang et al. Achieving shared multicast protection in WDM networks with sparse-light-splitting constraint
Souza Filho et al. Strategies for designing translucent wide-area networks

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
GR01 Patent grant
GR01 Patent grant