AUTOMATIC SEGMENTATION OF MEDICAL ULTRASONIC IMAGE USING SELF-CREATING AND ORGANIZING NEURAL NETWORK

1.   引言

突如其来的新冠疫情给人们的生活和工作带来了翻天覆地的变化。一夜之间，许多工作及娱乐都需要在网络上展开，催生出类似在线游戏、远程会议、网络直播等业务，这类业务的带宽大多具有随时间变化的特点^[1]。在大多数情况下，为了确保服务质量，网络运营商按照时变业务的峰值带宽需求为其分配频谱资源，从而忽略了时变业务带宽随时间变化的特点，当时变业务需求带宽降低时，网络分配的频谱资源得不到充分利用，导致严重的资源浪费^[2,3]。因此，根据时变业务特点定制资源分配方法是十分必要的。

弹性光网络(Elastic Optical Networks, EONs)能够有效地克服传统波分复用光网络粗糙粒度带宽资源划分和固定不变调制格式的缺点，以灵活的带宽分配和可变的调制格式提高资源的使用率，被广泛地认为是一个极具发展潜力的支持迅速发展的互联网数据和时变业务需求的智能光网络^[4]。在EONs中，路由和频谱分配(Routing and Spectrum Assignment, RSA)问题成为一个基本问题^[5]。研究者提出的大多数RSA方法没有考虑到时变业务特点^[6,7]。

对于EONs中时变业务的资源分配问题，Kumar等人^[8]提出了灵活频谱重分配和频谱扩展/压缩(Spectrum Expansion/Contraction, SEC)方式，其中的SEC可以在不中断业务传输的同时完成时变业务带宽的更新，极大地提高了频谱使用率。文献[9]提出了一种“频谱交易”的时变业务资源分配策略，该策略根据业务请求的累计信用实现频谱资源的共享，但是，当业务累计信用值小于设定阈值时，业务将被阻止使用网络中的频谱资源。为了进一步降低时变业务的阻塞率，Pathak等人^[10]提出，通过将相邻时变业务聚合，达到降低阻塞率的目的。针对时变业务带宽变化的方向问题，文献[11]提出了动态高扩展低收缩和动态交替方向(Dynamic Alternate Direction, DAD)策略缓解时变业务被阻塞的概率，但是，其交替扩展收缩的频谱分配策略对网络高负载时的资源分配显得较为“僵硬”，使资源使用率下降。

另外，生存性RSA是解决EONs的节点或者链路出现故障情况下保障业务服务传输质量的重要手段^[12-15]。Ding^[16]研究了生存性EONs中时变业务的路由和频谱分配问题，提出了基于专有路径保护策略的动态可生存性路由算法(Dynamic Survivable Path Routing Algorithm based on Dedicated Path Protection, DSPRA-DPP)，选择跳数少和频谱碎片小的工作和专有保护路径；为了降低时变业务的保护开销，进一步提出基于共享备份路径保护(Shared Backup Path Protection, SBPP)的动态生存性路由算法(Dynamic Survivable Path Routing Algorithm-SBPP, DSPRA-SBPP)。上述2种保护策略在频谱分配阶段采用首次命中(First Fit, FF) 频谱分配策略和中心频谱压缩扩展的SEC策略，导致业务SEC后频谱分配失败的风险较大。Paira等人^[17]设计了一种基于碎片感知的距离自适应生存性多路径的启发式路由频谱分配(Fragmentation-aware survivable Multipath Distance-Adaptive RSA, FMDA-RSA)策略，该策略采用距离自适应调制格式的多路径共享保护方法选择业务所需路由，为业务选择跳数少、频谱碎片化较小、调制格式较高的路由和频谱块，提高了业务生存性和频谱使用率。但是，FMDA-RSA侧重减少频谱碎片，忽略了频谱块的业务承载能力和SEC特性，对时变业务带宽扩展的适应能力较差。

针对生存性EONs中时变业务的路由和频谱分配问题，本文提出一种基于频谱窗滑动的时变业务共享保护(Time-varying Traffic Sharing Protection based on Spectrum Window Sliding, TTSP-SWS)算法。TTSP-SWS算法的贡献在于：(1)设计考虑可用频谱块承载权重和保护频谱块共享度的保护路径代价函数，用于选择时变业务的保护路径；(2)设计基于频谱窗滑动的保护路径的频谱分配策略为时变业务分配所需频谱块；(3)根据时变业务带宽变化情况，设计基于频谱窗滑动的频谱扩展/压缩的调整策略。

2.   问题描述

EONs拓扑抽象为G(V, E, F)，其中，V表示EONs节点集合，E表示光纤链路集合，F表示每条链路提供的频隙资源集合。时变业务表示为r(s, d, B, q)，其中，s,d分别表示时变业务的源、目的节点，B表示时变业务所需带宽，q表示时变业务需要的保护等级。

当时变业务r(s, d, B, q)到达EONs时，根据其源、目的节点以及所需带宽，采用最短路径算法确定两条链路分离的传输路径。

然后，为业务选择工作路径和保护路径物理传输距离约束的信号最高调制等级^[2,3,17]。

计算时变业务在保护路径上所需频隙(Frequency Slots, FSs)数目如式(1)

其中，C_f表示单位频隙的带宽，ρ_m表示信号调制等级，GB为业务之间的保护频隙。

最后，在频谱分配过程中，业务在工作和保护路径上分配的频谱应满足约束条件式(2)—式(5)

式(2)和式(3)保证了时变业务r(s, d, B, q)在工作路径和保护路径频谱分配都需要满足频谱一致性原则，其中，$P_{{\text{sd}}}^w$表示工作路径，$P_{{\text{sd}}}^p$为保护路径，$E(P_{{\text{sd}}}^w)$为工作路径的链路集合，$E(P_{{\text{sd}}}^p)$为保护路径的链路集合；${r_{e,{f_w}}}(P_{{\text{sd}}}^w)$表示工作路径的链路e的频隙f_w是否被占用，若占用，${r_{e,{f_w}}}(P_{{\text{sd}}}^w)$=1；若保护路径的链路e上分配的频隙f_p被占用，则${r_{e,{f_p}}} (P_{{\text{sd}}}^p)$=1；$F(P_{{\text{sd}}}^w)$和$F(P_{{\text{sd}}}^p)$分别表示工作路径和保护路径的频隙集合，式(4)和式(5)则保证了时变业务r(s, d, B, q)在工作路径和保护路径上分配的频谱需满足频谱连续性原则。

若不同业务的工作路径不相交，在单链路故障发生时，他们相交的保护路径中的频谱资源可以共享，以进一步减少频谱资源占用。

相较于带宽固定不变的业务，时变业务在带宽变化后，可以围绕着固定频率(Anchor Frequency, AF)进行动态扩展/压缩频谱。设$I({l_r})$,$I({u_r})$分别表示时变业务r(s, d, B, q)在其链路e的AF左、右两边的已经分配频谱块的索引值，则AF对应的已分配频谱块带宽为：${A_r} = I({u_r}) - I({l_r})$频隙；${L_r}$,${R_r}$分别表示时变业务r(s, d, B, q)在链路 e上其AF左、右两边潜在可以分配的空闲频隙数目，其计算公式为

式(6)中，R₁表示在链路 e上的时变业务r 的AF左边频隙索引位置的业务集合；式(7)中，R₂表示在链路 e上的时变业务r 的AF右边频隙索引位置的业务集合，|F|为每光纤频隙数目。

另外，根据时变业务的保护等级要求，在时变业务的所有保护路径上应至少提供$q \times B$的保护带宽资源，保证业务的生存性传输需求，即满足约束条件式(8)

其中，|P|表示时变业务r(s, d, B, q)的工作和保护路径数目，$B_r^k$表示时变业务r(s, d, B, q)第k条路径分配的带宽，$B_r^w$表示时变业务r(s, d, B, q)在工作路径中分配的带宽值。

EONs为了服务更多时变业务，在分配AF频谱的左、右两端预留一些空闲频隙支持时变业务的频谱扩展，导致频谱使用率降低，为满足动态时变业务频谱分配的实时性要求，如何在提高频谱使用率和减少生存性时变业务L_r和R_r频隙闲置和碎片化问题是时间复杂度极高的工作^[18]。为此，本文设计TTSP-SWS启发式算法求解上述动态时变业务在EONs中生存性传输时的RSA问题。

3.   基于频谱窗滑动的时变业务共享保护(TTSP-SWS)方法

TTSP-SWS包括：执行最短路径算法确定时变业务的工作路径，采用中心命中的FF方法在工作路径上分配业务的频隙资源；其次，设计基于频谱窗滑动的共享频谱生存性RSA策略，选择共享保护代价值最小的保护路径，选择频谱共享度值最大的频谱块；最后，当时变业务带宽需求发生变化时，根据当前路径中频谱资源使用情况以及时变业务占用频隙位置信息，利用DAD策略滑动频谱窗，实现工作路径和保护路径上频谱分配的压缩或扩展的调整。

3.1   基于频谱窗滑动的共享频谱生存性RSA策略

采用K最短路径算法确定时变业务的工作/保护边分离的保护路径集合。在保护路径的可用频谱块中，可用频谱块包括共享FS和空闲FS，从减少网络频谱资源使用角度，考虑尽可能使时变业务占用可共享的保护FS。为此，本文设计共享保护代价函数为

其中，$\left| {B_{\text{A}}^i} \right|$,$\left| {B_{\text{P}}^j} \right|$分别表示第i个可用频谱块和第j个保护频谱块的FS数目；N,M分别表示保护路径中可用频谱块和共享频谱块的数目；W_g表示第g个FS的共享度，共享度为该FS可保护业务的数量^[16]。

选择$C(P_{{\text{sd}}}^p)$值最小且时变业务在保护路径所需要的FS数目的候选路径为时变业务的保护路径$P_{{\text{sd}}}^p$。然后，根据所选保护路径$P_{{\text{sd}}}^p$中共享频谱块的FS共享状态，创建频谱窗(Spectrum Window, SW)，SW由一组连续的FSs构成，计算各SW频谱窗共享度，选择频谱窗共享度最大的频谱块分配给时变业务，其中，第m个SW频谱窗共享度为

利用频谱窗滑动方法，搜寻候选保护路径中的频谱窗共享度值最大频谱块，能提高保护路径中的频谱块利用率，降低生存性业务的阻塞率。基于频谱窗滑动的频谱共享生存性RSA策略见算法1。

算法1　基于频谱窗滑动的频谱共享生存性RSA策略
输入：弹性光网络拓扑G(V, E, F)，时变业务r(s, d, B, q)，工作路径$P_{{\text{sd}}}^w$以及频谱分配情况；
输出：时变业务r(s, d, B, q)的保护路径$P_{{\text{sd}}}^p$及其频谱分配结果；
(1) 将时变业务的工作路径$P_{{\text{sd}}}^w$从网络拓扑G(V, E, F)中移除；
(2) 执行K最短路径算法，根据式(9)，计算每条候选保护路径的共享保护代价值$C(P_{{\text{sd}}}^p)$，并升序排序候选保护路径在集合 　　 $P_{{\text{sd}}}^{{\text{ps}}} = \{ P_{{\text{sd}}}^j,j = 1,2,\cdots,K\}$中；
(3) 根据各候选保护路径的长度，计算确定业务所需最高调制等级，并由式(1)计算业务在各候选保护路径中所需FS数目f；
(4) 依次统计各候选保护路径中满足时变业务所需FS的可用频谱块，并保存在业务的可用频谱块集合Bj, A中；
(5) 若Bj, A集合非空，转步骤(6)；否则，阻塞时变业务r(s, d, B, q)，结束算法；
(6) 在非空的Bj, A集合中，选择路径序号j最小的可用频谱块的候选保护路径$P_{{\text{sd}}}^j$标识为$P_{{\text{sd}}}^p$，将$P_{{\text{sd}}}^p$中可用频谱块的共享FS保存在集合 　　 BP中，若BP中存在满足时变业务r(s, d, B, q)所需FS数目的频谱块；转至步骤(7)，否则，转至步骤(11)；
(7) 从BP中选择包含共享FS数目最多的共享频谱块，标识为BP,m，将共享频谱块BP,m的起始索引值记为index(FP)，末位索引值记为 　　 index(EP)；
(8) 在该保护路径上创建业务的滑动频谱窗SW, SW的大小由时变业务r(s, d, B, q)所需FS数目f决定；
(9) 根据式(10)，依次计算${\text{S} }{ {\text{W} }_{ {\text{index} }({F_{\text{P} } })} }{\text{～}}{\text{S} }{ {\text{W} }_{ {\text{index} }({E_{\text{P} } } - {f_{} } + 1)} }$的频谱窗共享度$W_m^s$值，选取$W_m^s$值最大的SW，并将时变业务r(s, d, B, q)所需 　　 保护频隙数目分配在频隙索引值区间$[m,m + f - 1]$，输出结果，结束算法1；如果存在$W_z^s = W_n^s$, $z,n \in m$，并且$z \ne n$，执行下一 　　 步骤；
(10) 计算$d_{\text{L}}^{\text{P}} = \varepsilon - {\text{index}}({F_{\text{P}}}),d_{\text{R}}^{\text{P}} = {\text{index}}({E_{\text{P}}}) - (\varepsilon + f - 1)$，${d_\varepsilon } = |d_{\text{R}}^{\text{P}} - d_{\text{L}}^{\text{P}}|$，其中$\varepsilon \in (z,n)$；选取${d_\varepsilon }$值最小的频谱窗，并将时变业务 　　 r(s, d, B, q)所需保护频隙数目分配在频隙索引值区间$[\varepsilon ,\varepsilon + f - 1]$，输出结果，结束算法1；
(11) 从保护路径$P_{{\text{sd}}}^p$的Bj, A中选择包含共享FS最多的可用频谱块BA, m作为预分配频谱块，该频谱块包含频隙数目为$|{\text{B}}{{\text{P}}_{{\text{max}}}}|$，并确定BA, m 　　 的起始索引值index(F)和末位索引值index(E)，并确定该可用频谱块内的最大共享频谱块的起始索引值index(FP)和末位索引值 　　 index(EP)；
(12) 计算$d_{\text{L}}^{\text{A}} = {\text{index}}({F_{\text{P}}}) - {\text{index}}(F)$, $d_{\text{R}}^{\text{A}} = {\text{index}}(E) - {\text{index}}({E_{\text{P}}})$，并计算$\varDelta = {f_{} } - |{\text{B} }{ {\text{P} }_{ {\text{max} } } }|$；
(13) 若$d_{\text{L}}^{\text{A}} > d_{\text{R}}^{\text{A}}$，则更新${\text{index}}({F_{\text{P}}}) = {\text{index}}({F_{\text{P}}}) - 1$；
(14) 若$d_{\text{L} }^{\text{A} } \le d_{\text{R} }^{\text{A} }$，则更新${\text{index}}({E_{\text{P}}}) = {\text{index}}({E_{\text{P}}}) + 1$；
(15) 更新$\varDelta = \varDelta - 1$；
(16) 若$\varDelta \ne 0$，返回步骤(12)；否则，将时变业务r(s, d, B, q)所需频隙数目分配在频隙索引值区间$[{\text{index}}({F_{\text{P}}}),{\text{index}}({E_{\text{P}}})]$，输出结果，结 　　 束算法1。

下载: 导出CSV 

| 显示表格

3.2   基于频谱窗滑动的扩展/压缩频谱分配策略

若时变业务的带宽B值变大，导致业务所需FS数目增加，即$f_{}^n > {f_{}}$，其中$f_{}^n$表示时变业务带宽发生变化后所需的FS数目，则工作路径和保护路径的频谱分配均需要做出改变；如果工作路径和保护路径当前占用的频谱块能够满足频谱扩展的要求，则工作路径的频谱扩展采用DAD策略即可^[11]。对于保护路径中频谱扩展，若时变业务当前占用的共享频谱块可以满足业务带宽扩展的需求，则根据生存性保护需所需FS数目，建立滑动频谱窗SW，从当前共享频谱块中选择频谱窗共享度最高的SW分配给业务；若时变业务当前占用的共享频谱块不能满足频谱扩展需求，则需要在保护路径上寻找可用频谱块中对生存性业务所需FS数目进行频谱窗的调整，在频谱窗扩展调整中，为了提高FS的共享度，尽量使用时变业务当前占用的共享频谱块，并将时变业务的共享FS滑动在可用频谱块的中间位置。

基于频谱窗滑动的生存性扩展频谱分配策略详细过程见算法2。

算法2　基于频谱窗滑动的生存性扩展频谱分配策略
输入：时变业务r(s, d, B, q)，带宽扩展后的时变业务r(s, d, Bn, q)，当前保护路径的FS状态；
输出：带宽变化后时变业务r(s, d, Bn, q)的频谱分配结果；
(1) 根据式(1)，计算带宽发生变化后时变业务r(s, d, Bn, q)在保护路径上所需保护FS的数目${f^n}$；
(2) 若时变业务r(s, d, Bn, q)当前占用保护路径的可用频谱块满足${f^n}$需求，即${B_{{\text{A}},{\text{ }}m}} \ge f_{}^n$，其中，BA, m表示当前时变业务r(s, d, Bn, q) 　　 在保护路径上占用可用频谱块的大小，并将该频谱块保存在集合EA, m中；否则，从业务的保护路径集合$P_{{\text{sd}}}^{{\text{ps}}}$的可用频谱块集合Bj, A中 　　 寻找满足${f^n}$需求的其他可用频谱块，并确定对应的保护路径$P_{{\text{sd}}}^j$，可用频谱块保存在集合EA, m中；若EA, m非空，转步骤(3)，否则，阻 　　 塞业务；
(3) 若EA, m的保护频谱块能满足${f^n}$需求，即${E_{{\text{P}},{\text{ }}m}} \ge f_{}^n$，则将该保护频谱块保存在集合EP, m中，转步骤(4)；转步骤(7)；
(4) 从EP, m中选择包含共享FS数目最多的共享频谱块，记录该共享频谱块的起始索引值记为index(FP)，末位索引值记为index(EP)，根据 　　 时变业务r(s, d, Bn, q)带宽变化所需FS数目${f^n}$，创建滑动频谱窗SW；
(5) 根据式(10)，依次计算${\text{S} }{ {\text{W} }_{ {\text{index} }({F_{\text{P} } })} }{\text{～}}{\text{S} }{ {\text{W} }_{ {\text{index} }({E_{\text{P} } } - {f^n} + 1)} }$的频谱窗共享度$W_m^s$值，选取$W_m^s$值最大的SW，并将时变业务r(s, d, Bn, q)所 　　 需保护频隙数目分配在频隙索引值区间$[m,m + {f^n} - 1]$，输出结果，结束算法2；如果存在$W_z^s = W_n^s$，$z,n \in m$，并且$z \ne n$，转步 　　 骤(6)；
(6) 计算$d_{\text{L}}^{\text{P}} = \varepsilon - {\text{index}}({F_{\text{P}}}){,_{}}d_{\text{R}}^{\text{P}} = {\text{index}}({E_{\text{P}}}) - (\varepsilon + {f^n} - 1)$, ${d_\varepsilon } = |d_{\text{R}}^{\text{P}} - d_{\text{L}}^{\text{P}}|$，其中$\varepsilon \in (z,n)$；选取${d_\varepsilon }$值最小的频谱窗，并将时变业务 　　 r(s, d, Bn, q)所需保护频隙数目分配在频隙索引值区间$[\varepsilon ,\varepsilon + {f_{}} - 1]$，输出结果，结束算法2；
(7) 从保护路径的可用频谱块保存在集合EA, m中选择包含共享FS最多的可用频谱块作为预分配频谱块，其大小记为BE，确定该可用频谱 　　 块的起始索引值index(F)和末位索引值index(E)，并确定该可用频谱块内的最大共享频谱块的起始索引值index(FP)和末位索引值 　　 index(EP)；
(8) 计算$d_{\text{L}}^{\text{E}} = {\text{index}}({F_{\text{P}}}) - {\text{index}}(F)$，$d_{\text{R}}^{\text{E}} = {\text{index}}(E) - {\text{index}}({E_{\text{P}}})$，并计算${\varDelta _{\text{E} } } = f_{}^n - {B_{\text{E} } }$；
(9) 如果$d_{\text{L}}^{\text{E}} > d_{\text{R}}^{\text{E}}$，则更新${\text{index}}({F_{\text{P}}}) = {\text{index}}({F_{\text{P}}}) + 1$；
(10) 如果$d_{\text{L}}^{\text{E}} \le d_{\text{R}}^{\text{E}}$，则更新${\text{index}}({E_{\text{P}}}) = {\text{index}}({E_{\text{P}}}) + 1$；
(11) 更新${\varDelta _{\text{E} } } = {\varDelta _{\text{E} } } - 1$；
(12) 如果${\varDelta _{\text{E} } } \ne 0$，转至步骤(8)；否则，将时变业务r(s, d, Bn, q)所需保护频隙数目$f_{}^n$分配在频隙索引值$[{\text{index}}({F_{\text{P}}}),{\text{index}}({E_{\text{P}}})]$区间，输 　　 出结果，结束算法2。

下载: 导出CSV 

| 显示表格

如果时变业务r(s, d, B, q)的带宽B值变小，则在工作路径和保护路径上分配的频谱块比业务实际需要的频谱块大，在工作路径上采用DAD策略进行频谱压缩^[11]；在保护路径上，在保证频谱一致性条件下，采用基于频谱窗滑动的生存性频谱压缩分配策略，滑动频谱窗，选择出共享度最高的频谱块分配给时变业务r(s, d, B, q)。基于频谱窗滑动的生存性压缩频谱分配策略详细过程见算法3。

算法3　基于频谱窗滑动的生存性压缩频谱分配策略
输入：时变业务r(s, d, B, q)，带宽压缩后的时变业务r(s, d, Bw, q)，当前保护路径资源分配结果；
输出：带宽压缩后时变业务r(s, d, Bw, q)的频谱分配结果；
(1) 根据式(1)，计算带宽压缩后时变业务r(s, d, Bw, q)在保护路径上所需保护FS的数目${f^c}$，将时变业务在当前保护路径分配频谱块的起始 　　 索引值index(F)和末位索引值index(E)；
(2) 计算时变业务带宽压缩后，在保护路径上需要减少的FS压缩${\varDelta _{\text{C}}} = {f_{}} - f_{}^c$；
(3) 计算时变业务r(s, d, B, q)在当前保护路径占用频谱块的起始FS的共享度${W_{{\text{index}}(F)}}$；若${W_{{\text{index}}(F)}} = 1$，转步骤(4)；否则，转步骤(7)；
(4) 更新频谱块的起始频隙索引值，${\text{index}}(F) = {\text{index}}(F) + 1$；
(5) 更新${\varDelta _{\text{C}}} = {\varDelta _{\text{C}}} - 1$；
(6) 若${\varDelta _{\text{C}}} = 0$，将时变业务带宽压缩后所需保护频隙数目${f^c}$分配在频隙索引值区间$[{\text{index}}(F),{\text{index}}(E)]$的频谱块上；若${\varDelta _{\text{C}}} \ne 0$，转至步 　　 骤(3)；
(7) 计算时变业务在当前保护路径占用频谱块的截止FS的共享度${W_{{\text{index}}(E)}}$；若${W_{{\text{index}}(E)}} = 1$，转步骤(8)；否则，转步骤(11)；
(8) 更新${\text{index}}(E) = {\text{index}}(E) - 1$；
(9) 更新${\varDelta _{\text{C}}} = {\varDelta _{\text{C}}} - 1$；
(10) 若${\varDelta _{\text{C}}} = 0$，将时变业务带宽压缩后所需保护频隙数目${f^c}$分配在频隙索引值区间$[{\text{index}}(F),{\text{index}}(E)]$的频谱块上，输出结果，结束 　　 表3算法；若${\varDelta _{\text{C}}} \ne 0$，转至步骤(7)；
(11) 根据时变业务r(s, d, Bw, q)带宽压缩后所需FS数目${f^c}$创建滑动频谱窗SW；
(12) 根据式(10)，依次计算${\text{S} }{ {\text{W} }_{ {\text{index} }(F)} }{\text{～}}{\text{S} }{ {\text{W} }_{ {\text{index} }(E - f_{}^c + 1)} }$的频谱共享度$W_m^s$值，选取$W_m^s$值最大的SW，将时变业务带宽压缩后所需保护 　　 频隙数目${f^c}$分配在频隙索引值区间$[m,m + f_{}^c - 1]$的频谱块上，输出结果，结束算法3；如果存在$W_z^s = W_n^s$, $z,n \in m$并且$z \ne n$， 　　 转至步骤(13)；
(13) 计算$d_{\text{L}}^{\text{C}} = \varepsilon - {\text{index}}(F){,_{}}d_{\text{R}}^{\text{C}} = {\text{index}}(E) - (\varepsilon + f_{}^c - 1)$, $d_\varepsilon ^{\text{C}} = |d_{\text{R}}^{\text{C}} - d_{\text{L}}^{\text{C}}|$，其中$\varepsilon \in (z,n)$；选取最小值$d_\varepsilon ^{\text{C}}$对应的ε值，并将时变业务 　　 带宽压缩后所需保护频隙数目${f^c}$分配在频隙索引值区间$[\varepsilon ,\varepsilon + f_{}^c - 1]$的频谱块上，输出结果，结束算法3。

下载: 导出CSV 

| 显示表格

4.   仿真分析

4.1   仿真环境及评价指标

为了验证本文所提TTSP-SWS算法的性能，本文分别对DSPRA-DPP^[16], DSPRA-SBPP^[16]和FMDA-RSA^[17]生存性策略在图1所示的国家科学基金网(National Science Foundation Network, NSFNET)和美国网络(United States of America Network, USNET)拓扑^[17]中的业务阻塞率、保护冗余度和频谱使用率性能进行仿真，其中，NSFNET拓扑具有14个节点，21条链路，USNET拓扑包含24个节点，43条链路^[17]，设K=3，业务间保护频隙GB=1 fs，图1的链路旁边数字表示节点之间的物理长度，单位km。其他默认仿真主要参数如表1。

图  1  仿真网络拓扑图

下载: 全尺寸图片 幻灯片

表  1  默认仿真参数

参数	数值
时变业务请求数目	105 个
每条链路提供的频隙数目F	320 fs
单位频隙提供带宽Cf	12.5 GHz
业务保护等级q	(0, 1)
业务大小	[25, 125] Gb/s
业务请求到达率	服从泊松分布
业务请求持续时间	服从负指数分布

下载: 导出CSV 

| 显示表格

4.2   仿真结果及分析

图2显示了不同负载情况下的业务阻塞率性能。随着业务负载的增多，图2中所有算法的业务阻塞率都增加，这是因为网络中的可用频谱资源随着负载的增加而减少，当负载较大时，网络中的资源竞争加剧，业务在生存性路由选择和频谱分配中失败的概率增加。相比其他算法，本文所提TTSP-SWS算法在相同负载下的业务阻塞率最低，这是因为：TTSP-SWS算法在保护路径选择时，优先考虑可用或者保护频谱资源更集中的路径，选择共享保护代价值更小的保护路径；在频谱分配时，利用滑动频谱窗，将时变业务所需的FS分配在保护路径中频隙共享度较高的位置，提高了频隙的共享度和频谱窗的利用率。而对比的DSPRA-DPP算法的业务阻塞率最高，原因是DSPRA-DPP算法采用专有路径保护，频谱使用率较低；与之对应的是，DSPRA-SBPP算法，FMDA-RSA方法和TTSP-SWS算法采用共享路径保护，相较于DSPRA-DPP专有路径保护，共享路径保护可以减少保护资源的使用。相比DSPRA-SBPP算法，采用多路径生存性策略和最小频谱碎片的FMDS-RSA算法不能适应时变业务的SEC，当业务带宽变化时，FMDS-RSA算法使频谱碎片率和业务阻塞率增加。

图  2  业务阻塞率随负载的变化

下载: 全尺寸图片 幻灯片

图2(b)是4种算法在USNET拓扑中的业务阻塞率性能，各算法在USNET中的业务阻塞率曲线呈现与图2(a)相同的变化趋势；不同的是，在相同负载情况下，NSFNET中各算法的业务阻塞率高于USNET，原因是USNET具有更多的节点和链路数，网络的连通度更高，时变业务的路由成功概率更高。

图3展示了4种算法在不同网络负载下的频谱使用率性能。随着业务负载增加，EONs的频谱资源被更多业务使用，提高了频谱使用率。在负载较大时，本文所提TTSP-SWS算法在NSFNET和USNET拓扑中均获得最大的频谱使用率，这是因为：TTSP-SWS算法不仅在业务路由选择阶段考虑网络中保护频谱块的频隙共享情况，在频谱分配阶段以及时变业务的带宽扩展或带宽压缩中，采用基于频谱窗滑动的生存性频谱扩展或频谱压缩策略，能将带宽变化后的业务分配在保护路径的频谱窗共享度高的频谱块位置，提高频谱资源的利用率。相较于DSPRA-SBPP算法，虽然DSPRA-DPP的冗余保护需要消耗更多的频谱资源，但是其业务阻塞率较大，所以其频谱使用率低于DSPRA-SBPP算法；而FMDA-RSA算法虽然在低负载时的业务阻塞率略低于DSPRA-SBPP算法，由于其频谱分配采用碎片感知的最小频偏率频谱块方法，减少了碎片频谱产生概率，因而其频谱使用率略低。相比于采用共享保护的DSPRA-SBPP算法和FMDS-RSA算法，TTSP-SWS算法在频谱扩展或压缩中考虑了时变业务带宽调整方向，增加了相邻时变业务之间的频谱窗共享度，进一步提高了共享频谱资源的利用率。

图  3  频谱使用率随负载的变化

下载: 全尺寸图片 幻灯片

图4是DSPRA-DPP, DSPRA-SBPP, FMDS-RSA和TTSP-SWS算法的保护冗余度性能对比。由于DSPRA-DPP算法采用专有路径保护，需要为业务预留更多的保护频谱资源开销，所以DSPRA-DPP的保护冗余度在4种算法中最高。相较于DSPRA-SBPP和FMDS-RSA共享保护算法，TTSP-SWS算法考虑了保护路径中频隙的共享度，并根据可用频谱块以及保护频谱块的大小，为时变业务选择共享保护代价值较小的保护路径，提高了保护路径上的频谱使用率，降低了生存性保护的资源开销；同时，在频谱分配和带宽调整中，进一步比较保护路径上频隙的共享度，选择频谱窗共享度大的频谱块，进一步降低了业务使用频谱资源的数目。另外，随着网络中负载的增多，DSPRA-SBPP, FMDS-RSA和TTSP-SWS算法的备份路径频谱保护冗余度呈现出缓慢下降的趋势，这是因为随着请求数目的增加，被业务预留的保护频谱资源增多，增加了保护频隙被其他业务共享概率，从而降低了生存性业务预留的保护频谱资源冗余度。

图  4  保护冗余度随负载的变化

下载: 全尺寸图片 幻灯片

5.   结论

针对在EONs中时变业务高阻塞率和频谱使用率低的问题，本文提出一种基于频谱窗滑动的时变业务共享保护(TTSP-SWS)方法，为业务选择频谱使用率高和资源共享度高的保护路径；同时，通过滑动频谱窗，为业务选择频谱窗共享度高的频谱块以适应时变业务带宽扩展和压缩情况，提高了频谱资源的使用率。随着互联网中时变业务发展和生存性传输需求的增长，通过共享频谱的方式提高时变生存性业务可靠性和降低频谱资源的保护开销，对进一步增强EONs网络服务能力和质量保证能力具有积极意义。