一种5G切片接入的配电网保护优化方法及系统与流程

一种5g切片接入的配电网保护优化方法及系统
技术领域
1.本发明涉及配电网保护技术领域，具体来说是一种5g切片接入的配电网保护优化方法及系统。


背景技术：

2.近年来，5g技术迅速发展，其所带来的超高可靠与低延迟通信业务(ultra-reliable low latency communications，urllc)具有高带宽、低时延等优点。在配电网的差动保护中应用5g技术，为故障的定位、隔离与恢复提供了技术优势。但5g通信的传输时延及抖动仍然会对差动保护判据计算的不利影响，产生保护误动和保护拒动等问题，不利于整个配电网系统的运行。而5g延时抖动变化较大的主要原因是网络资源分配不均所导致。
3.针对上述问题，申请号为cn201911175988.4公开的一种基于5g网络的配电网分布式差动保护方法及系统[zh]本发明公开了一种基于5g网络的配电网分布式差动保护方法及系统，包括：数字化保护装置之间在信号传输时，利用基于5g网络的通信信道进行通信，从而实现两个数字化保护装置数据的交互；利用5g切片网络作为通信信道，保证了通信的安全性，在5g通信本身低时延特性的基础上，通过切片网络优化数据链路，以最小化链路传输跳数，进一步降低了数据传输的时延；利用数字化保护装置对配电网线路两端的电流同时进行突变检测，以各自的突变时刻作为线路两端数字化保护装置的同步时间基准，一旦满足突变量启动判据，则数字化保护装置启动，并向对端发送启动信号，数字化保护装置启动以后，计算故障启动后的电流故障分量值，计算完成后向对端发送电流故障分量数据。该通过切片网络优化数据链路，以最小化链路传输跳数，以达到降低了数据传输的时延。但该方法对于外部扰动不具有鲁棒性。


技术实现要素：

[0004]
本发明所要解决的技术问题在于提供一种5g切片接入的配电网保护优化方法，实现外部移动设备的接入而产生的延时抖动具有一定的鲁棒性，使配电网的保护动作更加精准。
[0005]
本发明通过以下技术手段实现解决上述技术问题的：
[0006]
一种5g切片接入的配电网保护优化方法，包括以下步骤：
[0007]
步骤1，对采用5g通信网络的配电网系统进行考虑通信延时的建模：
[0008]
设配电网中有i个数据采样节点，sn(i)为第i个采样节点当前服务的切片id；在采样节点的k个可接入的接入点中，第k个接入点部署的切片集合中有切片sm(k)，m∈{1,2,
…
,m}，其中m表示切片集合大小；bm(k)为sm(k)对应的带宽；为使系统的资源利用率最高，将所有数据采样节点均匀分配资源；因此，第i个采样节点的数据传输速率如式(1)所示；
[0009][0010]
式(1)中，am(k)表示当前切片sm(k)已接入的采样节点数；d
ik
表示采样节点i与待
接入的接入点之间的距离，其传播损耗用pl(d
ik
)表示；n0为单位带宽的噪声功率；pk表示待接入的接入点的发射功率；δ表示待接入ap收到的总干扰；
[0011]
步骤2，以最小化时延为优化目标，目标函数如式(2)所示；
[0012][0013]
其中，dsi为第i个采样节点采样的数据包大小；
[0014]
约束条件如式(3)-式(6)所示：
[0015][0016]
sm(k)≥s
n(i)ꢀꢀꢀ
(4)
[0017]rik
(m)≥ri,i∈i,k∈k,m∈m
ꢀꢀꢀ
(5)
[0018][0019]
式(2)中，x
ik
(m)为需要求解的0-1变量，当x
ik
(m)＝1时表示第i个数据采样节点接入第k个ap的第m个切片，否则表示不接入；约束条件(3)表示一个数据采样节点在某一时刻只能接入一个切片；约束条件(4)表示数据采样节点所接入切片的qos服务等级应该不小于原服务切片；约束条件(5)表示数据采样节点接入切片后的数据传输速率应大于等于其最低速率要求ri；约束条件(6)表示待接入切片所能容纳的节点数应该小于等于接入前的剩余节点数um(k)；
[0020]
步骤3，用拉格朗日乘数法来求解步骤2中目标函数的最优解，
[0021]
步骤4，基于步骤3的最优解，用分支定界法得到所述目标函数的0-1求解变量x
0-1
，从而根据所述求解变量决定用户是否接入相应的切片。
[0022]
本发明提出一种5g切片接入的配电网保护优化方法。采用时延最小化的数据传输方案，得出相应的数学模型，随后根据所提方案，采用拉格朗日乘数法和分支定界法进行模型求解，使5g网络资源均匀分配。采用此种方法可降低配电网保护通信过程中产生的延时。本发明外部设备可视为干扰的一部分，融入到数据传输速率的公式中，可优化外部可移动设备接入到本5g网络时产生的抖动延时，对外部干扰具有一定的鲁棒性。本发明通过采用一种5g切片接入的配电网保护优化方法，克服由5g低延时和抖动所导致配电网保护动作误动和拒动问题，实现配电网保护动作进一步精准，使整个配电网系统更加稳定安全的运行。
[0023]
进一步的，所述步骤3具体步骤如下：根据式(7)建立所述非线性规划问题的拉格朗日函数：
[0024][0025]
式(7)中，l(x
ik
(m),λ)表示关于x
ik
(m)和拉格朗日乘数λ的拉格朗日函数，h1(x
ik
(m))，h2(x
ik
(m))，h3(x
ik
(m))，h4(x
ik
(m))，h5(x
ik
(m))分别表示5个约束函数，λ1、λ2、λ3、λ4、λ5表示对应约束函数的拉格朗日乘数，并有：
[0026][0027]
h2(x
ik
(m))＝x
ik
(m)-1
ꢀꢀꢀ
(9)
[0028]
h3(x
ik
(m))＝sn(i)-sm(k)
ꢀꢀꢀ
(10)
[0029]
h4(x
ik
(m))＝r
i-r
ik
(m)
ꢀꢀꢀ
(11)
[0030][0031]
根据式(13)-(16)建立卡罗需-库恩-塔克条件，从而通过联合卡罗需-库恩-塔克条件相关等式，求出经过松弛的非线性规划问题的最优解x
relax
：
[0032][0033]
h1(x
ik
(m))＝0，h2(x
ik
(m))≤0，h3(x
ik
(m))≤0，h4(x
ik
(m))≤0，h5(x
ik
(m))≤0
ꢀꢀ
(14)
[0034]
λ1，λ2，λ3，λ4，λ5≥0
ꢀꢀꢀ
(15)
[0035][0036]
进一步的，所述步骤4具体步骤如下：
[0037]
步骤4.1，将式(2)作为问题p-1；分支定界法的输入为：满足卡罗需-库恩-塔克条件的松弛问题最优解x
relax
，松弛问题最优目标函数值z
relax
，0-1任意值ε；初始化k＝0，l＝0，u＝z
relax
；
[0038]
步骤4.2，从最优解x
relax
中任意选择一个不符合0-1约束条件的解xj，即：xj∈(0,1)；
[0039]
步骤4.3，若0≤xj《ε成立，则将约束条件xj＝0加到问题p-1中，形成子问题ⅰ；否则，将约束条件xj＝1加到问题p-1中，形成子问题ⅱ，ε表示0到1内任意值；
[0040]
步骤4.4，k ，继续求出子问题ⅰ或子问题ⅱ的松弛问题解，记为xk，并将对应的最优目标函数值记为zk；
[0041]
步骤4.5,找出最优目标函数最大值u作为新的上界，即：
[0042]
u＝max{zk′
|k
′
＝1,2,
…
,k},xk′
∈[0,1]；
[0043]
步骤4.6,再从符合0-1条件的分支中，找出目标函数最大值l作为新的下界，即：
[0044]
l＝max{zk′
|k
′
＝1,2,
…
,k},xk′
∈{0,1}；
[0045]
步骤4.7,若zk′
《l，剪掉相应分支,否则zk′
＞l，且不符合0-1条件,返回步骤4.3；
[0046]
步骤4.8,否则就意味着所有分支的最优目标函数值等于下界：zk′
＝l,将zk′
赋值z
0-1
，将xk′
赋值给x
0-1
，并作为问题p-1的最优解，其中，x
0-1
表示符合0-1约束条件的最优解。
[0047]
与上述方法对应的，本发明还提供一种5g切片接入的配电网保护优化系统，包括：
[0048]
建模模块，对采用5g通信网络的配电网系统进行考虑通信延时的建模：
[0049]
设配电网中有i个数据采样节点，sn(i)为第i个采样节点当前服务的切片id；在采样节点的k个可接入的接入点中，第k个接入点部署的切片集合中有切片sm(k)，m∈{1,
2,
…
,m}，其中m表示切片集合大小；bm(k)为sm(k)对应的带宽；为使系统的资源利用率最高，将所有数据采样节点均匀分配资源；因此，第i个采样节点的数据传输速率如式(1)所示；
[0050][0051]
式(1)中，am(k)表示当前切片sm(k)已接入的采样节点数；d
ik
表示采样节点i与待接入的接入点之间的距离，其传播损耗用pl(d
ik
)表示；n0为单位带宽的噪声功率；pk表示待接入的接入点的发射功率；δ表示待接入ap收到的总干扰；
[0052]
目标函数建立模块，以最小化时延为优化目标，目标函数如式(2)所示；
[0053][0054]
其中，dsi为第i个采样节点采样的数据包大小；
[0055]
约束条件如式(3)-式(6)所示：
[0056][0057]
sm(k)≥s
n(i)ꢀꢀꢀ
(4)
[0058]rik
(m)≥ri,i∈i,k∈k,m∈m
ꢀꢀꢀ
(5)
[0059][0060]
式(2)中，x
ik
(m)为需要求解的0-1变量，当x
ik
(m)＝1时表示第i个数据采样节点接入第k个ap的第m个切片，否则表示不接入；约束条件(3)表示一个数据采样节点在某一时刻只能接入一个切片；约束条件(4)表示数据采样节点所接入切片的qos服务等级应该不小于原服务切片；约束条件(5)表示数据采样节点接入切片后的数据传输速率应大于等于其最低速率要求ri；约束条件(6)表示待接入切片所能容纳的节点数应该小于等于接入前的剩余节点数um(k)；
[0061]
目标函数求解模块，用拉格朗日乘数法来求解步骤2中的目标函数，
[0062]
判断模块，目标函数的解，用分支定界法得到所述目标函数的0-1求解变量x
0-1
，从而根据所述求解变量决定用户是否接入相应的切片。
[0063]
进一步的，所述目标函数求解模块具体步骤如下：根据式(7)建立所述非线性规划问题的拉格朗日函数：
[0064][0065]
式(7)中，l(x
ik
(m),λ)表示关于x
ik
(m)和拉格朗日乘数λ的拉格朗日函数，h1(x
ik
(m))，h2(x
ik
(m))，h3(x
ik
(m))，h4(x
ik
(m))，h5(x
ik
(m))分别表示5个约束函数，λ1、λ2、λ3、λ4、λ5表示对应约束函数的拉格朗日乘数，并有：
[0066]
[0067]
h2(x
ik
(m))＝x
ik
(m)-1
ꢀꢀꢀ
(9)
[0068]
h3(x
ik
(m))＝sn(i)-sm(k)
ꢀꢀꢀ
(10)
[0069]
h4(x
ik
(m))＝r
i-r
ik
(m)
ꢀꢀꢀ
(11)
[0070][0071]
根据式(13)-(16)建立卡罗需-库恩-塔克条件，从而通过联合卡罗需-库恩-塔克条件相关等式，求出经过松弛的非线性规划问题的最优解x
relax
：
[0072][0073]
h1(x
ik
(m))＝0，h2(x
ik
(m))≤0，h3(x
ik
(m))≤0，h4(x
ik
(m))≤0，h5(x
ik
(m))≤0
ꢀꢀꢀ
(14)
[0074]
λ1，λ2，λ3，λ4，λ5≥0
ꢀꢀꢀ
(15)
[0075][0076]
进一步的，所述判断模块具体步骤如下：
[0077]
步骤4.1，将式(2)作为问题p-1；分支定界法的输入为：满足卡罗需-库恩-塔克条件的松弛问题最优解x
relax
，松弛问题最优目标函数值z
relax
，0-1任意值ε；初始化k＝0，l＝0，u＝z
relax
；
[0078]
步骤4.2，从最优解x
relax
中任意选择一个不符合0-1约束条件的解xj，即：xj∈(0,1)；
[0079]
步骤4.3，若0≤xj《ε成立，则将约束条件xj＝0加到问题p-1中，形成子问题ⅰ；否则，将约束条件xj＝1加到问题p-1中，形成子问题ⅱ，ε表示0到1内任意值；
[0080]
步骤4.4，k ，继续求出子问题ⅰ或子问题ⅱ的松弛问题解，记为xk，并将对应的最优目标函数值记为zk；
[0081]
步骤4.5,找出最优目标函数最大值u作为新的上界，即：
[0082]
u＝max{zk′
|k
′
＝1,2,
…
,k},xk′
∈[0,1]；
[0083]
步骤4.6,再从符合0-1条件的分支中，找出目标函数最大值l作为新的下界，即：
[0084]
l＝max{zk′
|k
′
＝1,2,
…
,k},xk′
∈{0,1}；
[0085]
步骤4.7,若zk′
《l，剪掉相应分支,否则zk′
＞l，且不符合0-1条件,返回步骤4.3；
[0086]
步骤4.8,否则就意味着所有分支的最优目标函数值等于下界：zk′
＝l,将zk′
赋值z
0-1
，将xk′
赋值给x
0-1
，并作为问题p-1的最优解，其中，x
0-1
表示符合0-1约束条件的最优解。
[0087]
本发明还提供一种处理设备，包括至少一个处理器，以及与所述处理器通信连接的至少一个存储器，其中：所述存储器存储有可被处理器执行的程序指令，所述处理器调用所述程序指令能够执行上述的方法。
[0088]
本发明还提供一种计算机可读存储介质，所述计算机可读存储介质存储计算机指令，所述计算机指令使所述计算机执行上述的方法。
[0089]
本发明的优点在于：
[0090]
本发明提出一种5g切片接入的配电网保护优化方法。采用时延最小化的数据传输方案，得出相应的数学模型，随后根据所提方案，采用拉格朗日乘数法和分支定界法进行模型求解，使5g网络资源均匀分配。采用此种方法可降低配电网保护通信过程中产生的延时。本发明外部设备可视为干扰的一部分，融入到数据传输速率的公式中，可优化外部可移动设备接入到本5g网络时产生的抖动延时，对外部干扰具有一定的鲁棒性。本发明通过采用一种5g切片接入的配电网保护优化方法，克服由5g低延时和抖动所导致配电网保护动作误动和拒动问题，实现配电网保护动作进一步精准，使整个配电网系统更加稳定安全的运行。
附图说明
[0091]
图1为本发明实施例中5g切片接入的配电网保护优化方法的流程框图；
[0092]
图2为本发明实施例中5g切片接入的配电网保护优化方法应用场景；
[0093]
图3为本发明实施例中实例1未采用优化方法和采用本发明优化方法的效果对比；
[0094]
图4为本发明实施例中实例25g传输时延在有无用户接入时的效果对比。
具体实施方式
[0095]
为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
[0096]
本实施例提供一种5g切片接入的配电网保护优化方法，包括以下步骤：
[0097]
步骤1，对采用5g通信网络的配电网系统进行考虑通信延时的建模：
[0098]
本实施例所考虑使用5g通信网络进行配电网差动保护的场景如图1所示。5g配电网保护包括1个5g基站和多个接入该基站的保护设备。每个保护设备都具有电流采样的能力，安装在线路的两侧，进行差动保护电流数据的采样。5g基站提供5g网络，为保护设备之间提供通信服务，同时也为其它用户接入的移动端设备提供网络服务。可通过5g通信的方式实现保护设备之间线路电流数据的对比，从而判断该线路是否要进行差动保护动作的判断。
[0099]
设配电网中有i个数据采集节点，dsi为第i个节点采样的数据包大小，sn(i)为第i个采样节点当前服务的切片id；在采集节点的k个可接入的接入点(access point,ap)中，第k个ap部署的切片集合中有任意的一类切片sm(k)，m∈{1,2,
…
,m}，其中m表示切片集合大小；bm(k)为sm(k)对应的带宽。为使系统的资源利用率最高，将所有数据采集节点均匀分配资源。因此，第i个节点的数据传输速率如式(1)所示。
[0100][0101]
式(1)中，am(k)表示当前切片sm(k)已接入的节点数；d
ik
表示节点i与待接入ap之间的距离，其传播损耗用pl(d
ik
)表示；n0为单位带宽的噪声功率；pk表示待接入ap的发射功率；δ表示待接入ap收到的总干扰。
[0102]
步骤2，为了使数据采集节点在数据传输时，最大化地利用切片资源，满足时延和
吞吐量要求，以最小化时延为优化目标，目标函数如式(2)所示。
[0103][0104]
约束条件如式(3)-式(6)所示：
[0105][0106]
sm(k)≥s
n(i)ꢀꢀꢀ
(4)
[0107]rik
(m)≥ri,i∈i,k∈k,m∈m
ꢀꢀꢀ
(5)
[0108][0109]
式(2)中，x
ik
(m)为需要求解的0-1变量，当x
ik
(m)＝1时表示第i个数据采集节点接入第k个ap的第m个切片，否则表示不接入；约束条件(3)表示一个数据采集节点在某一时刻只能接入一个切片；约束条件(4)表示数据采集节点所接入切片的qos服务等级应该不小于原服务切片；约束条件(5)表示数据采集节点接入切片后的数据传输速率应大于等于其最低速率要求ri；约束条件(6)表示待接入切片所能容纳的节点数应该小于等于接入前的剩余节点数um(k)。
[0110]
步骤3，用拉格朗日乘数法来求解上述模型，根据式(7)建立所述非线性规划问题的拉格朗日函数：
[0111][0112]
式(7)中，l(x
ik
(m),λ)表示关于x
ik
(m)和拉格朗日乘数λ的拉格朗日函数，h1(x
ik
(m))，h2(x
ik
(m))，h3(x
ik
(m))，h4(x
ik
(m))，h5(x
ik
(m))分别表示5个约束函数，λ1、λ2、λ3、λ4、λ5表示对应约束函数的拉格朗日乘数，并有：
[0113][0114]
h2(x
ik
(m))＝x
ik
(m)-1
ꢀꢀꢀ
(9)
[0115]
h3(x
ik
(m))＝sn(i)-sm(k)
ꢀꢀꢀ
(10)
[0116]
h4(x
ik
(m))＝r
i-r
ik
(m)
ꢀꢀꢀ
(11)
[0117][0118]
根据式(13)-(16)建立kkt条件(karush-kuhn-tucker conditions，卡罗需-库恩-塔克条件)，从而通过联合kkt条件相关等式，求出经过松弛的非线性规划问题的最优解x
relax
：
[0119][0120]
h1(x
ik
(m))＝0，h2(x
ik
(m))≤0，h3(x
ik
(m))≤0，h4(x
ik
(m))≤0，h5(x
ik
(m))≤0
ꢀꢀꢀ
(14)
[0121]
λ1，λ2，λ3，λ4，λ5≥0
ꢀꢀꢀ
(15)
[0122][0123]
步骤4，用分支定界法得到所述目标函数的0-1求解变量x
0-1
，从而根据所述求解变量决定用户是否接入相应的切片。具体步骤如下：
[0124]
步骤4.1，将式(2)作为问题p-1。分支定界法的输入为：满足kkt条件的松弛问题最优解x
relax
，松弛问题最优目标函数值z
relax
，0-1任意值ε。初始化k＝0，l＝0，u＝z
relax
。
[0125]
步骤4.2，从最优解x
relax
中任意选择一个不符合0-1约束条件的解xj，即：xj∈(0,1)；
[0126]
步骤4.3，若0≤xj《ε成立，则将约束条件xj＝0加到问题p-1中，形成子问题ⅰ；否则，将约束条件xj＝1加到问题p-1中，形成子问题ⅱ，ε表示0到1内任意值；
[0127]
步骤4.4，k ，继续求出子问题ⅰ或子问题ⅱ的松弛问题解，记为xk，并将对应的最优目标函数值记为zk；
[0128]
步骤4.5,找出最优目标函数最大值u作为新的上界，即：
[0129]
u＝max{zk′
|k
′
＝1,2,
…
,k},xk′
∈[0,1]；
[0130]
步骤4.6,再从符合0-1条件的分支中，找出目标函数最大值l作为新的下界，即：
[0131]
l＝max{zk′
|k
′
＝1,2,
…
,k},xk′
∈{0,1}；
[0132]
步骤4.7,若zk′
《l，剪掉相应分支,否则zk′
＞l，且不符合0-1条件,返回步骤4.3；
[0133]
步骤4.8,否则就意味着所有分支的最优目标函数值等于下界：zk′
＝l,将zk′
赋值z
0-1
，将xk′
赋值给x
0-1
，并作为问题p-1的最优解，其中，x
0-1
表示符合0-1约束条件的最优解。
[0134]
本实施例公开了一种5g切片接入的配电网保护优化方法，针对于配电网保护采集电流数据通过5g传输时产生的时延和抖动的优化问题，首先对5g通信在配电网中使用场景进行分析，并建立5g通信接入延时的数学模型；然后根据该模型，研究差动信号时延最小化的优化目标函数，并依据所提优化目标函数采用拉格朗日乘数法和分支定界法求解最大化利用通信信道资源的5g切片选择，实现差动信号传输的时延最小，从而提升电网系统的安全性，保证电网系统的稳定运行。
[0135]
以下结合具体的实施例对以上步骤进行说明。
[0136]
实例1：对采用本发明优化方法与未采用优化方法产生的时延曲线进行对比，两种情况的比较如图2所示，在没有采用本文最小化时延方法的情况下，通过5g传输产生的时延会随时间的变化而随机变化，不利于差动保护动作的判定。采用本文使用的方法后，通过5g传输产生的时延会随时间的变化逐步减小，最终趋于一个恒定值。由上述的对比，可以体现本文所提出最小化优化方法的优越性。
[0137]
实例2：在采用本发明所提的优化方法下，在1s和7s分别有两次100个移动端设备
接入，5g基站会分配一定的网络资源供其使用，这一过程会使差动保护传输的时延略有增大，之后由于采用算法的优化作用，会使产生的时延趋于平稳，与没有移动端设备接入的情况近乎相同，表现出所提方法的鲁棒性。具体的结果如图3所示。
[0138]
本实施例还提供一种5g切片接入的配电网保护优化系统，与上述方法对应，包括：
[0139]
建模模块，对采用5g通信网络的配电网系统进行考虑通信延时的建模：
[0140]
设配电网中有i个数据采样节点，sn(i)为第i个采样节点当前服务的切片id；在采样节点的k个可接入的接入点中，第k个接入点部署的切片集合中有切片sm(k)，m∈{1,2,
…
,m}，其中m表示切片集合大小；bm(k)为sm(k)对应的带宽；为使系统的资源利用率最高，将所有数据采样节点均匀分配资源；因此，第i个采样节点的数据传输速率如式(1)所示；
[0141][0142]
式(1)中，am(k)表示当前切片sm(k)已接入的采样节点数；d
ik
表示采样节点i与待接入的接入点之间的距离，其传播损耗用pl(d
ik
)表示；n0为单位带宽的噪声功率；pk表示待接入的接入点的发射功率；δ表示待接入ap收到的总干扰；
[0143]
目标函数建立模块，以最小化时延为优化目标，目标函数如式(2)所示；
[0144][0145]
其中，dsi为第i个采样节点采样的数据包大小；
[0146]
约束条件如式(3)-式(6)所示：
[0147][0148]
sm(k)≥s
n(i)ꢀꢀꢀ
(4)
[0149]rik
(m)≥ri,i∈i,k∈k,m∈m
ꢀꢀꢀ
(5)
[0150][0151]
式(2)中，x
ik
(m)为需要求解的0-1变量，当x
ik
(m)＝1时表示第i个数据采样节点接入第k个ap的第m个切片，否则表示不接入；约束条件(3)表示一个数据采样节点在某一时刻只能接入一个切片；约束条件(4)表示数据采样节点所接入切片的qos服务等级应该不小于原服务切片；约束条件(5)表示数据采样节点接入切片后的数据传输速率应大于等于其最低速率要求ri；约束条件(6)表示待接入切片所能容纳的节点数应该小于等于接入前的剩余节点数um(k)；
[0152]
目标函数求解模块，用拉格朗日乘数法来求解步骤2中的目标函数，
[0153]
判断模块，目标函数的解，用分支定界法得到所述目标函数的0-1求解变量x
0-1
，从而根据所述求解变量决定用户是否接入相应的切片。
[0154]
所述目标函数求解模块具体步骤如下：根据式(7)建立所述非线性规划问题的拉格朗日函数：
[0155][0156]
式(7)中，l(x
ik
(m),λ)表示关于x
ik
(m)和拉格朗日乘数λ的拉格朗日函数，h1(x
ik
(m))，h2(x
ik
(m))，h3(x
ik
(m))，h4(x
ik
(m))，h5(x
ik
(m))分别表示5个约束函数，λ1、λ2、λ3、λ4、λ5表示对应约束函数的拉格朗日乘数，并有：
[0157][0158]
h2(x
ik
(m))＝x
ik
(m)-1
ꢀꢀꢀ
(9)
[0159]
h3(x
ik
(m))＝sn(i)-sm(k)
ꢀꢀꢀ
(10)
[0160]
h4(x
ik
(m))＝r
i-r
ik
(m)
ꢀꢀꢀ
(11)
[0161][0162]
根据式(13)-(16)建立卡罗需-库恩-塔克条件，从而通过联合卡罗需-库恩-塔克条件相关等式，求出经过松弛的非线性规划问题的最优解x
relax
：
[0163][0164]
h1(x
ik
(m))＝0，h2(x
ik
(m))≤0，h3(x
ik
(m))≤0，h4(x
ik
(m))≤0，h5(x
ik
(m))≤0
ꢀꢀꢀ
(14)
[0165]
λ1，λ2，λ3，λ4，λ5≥0
ꢀꢀꢀ
(15)
[0166][0167]
进一步的，所述判断模块具体步骤如下：
[0168]
步骤4.1，将式(2)作为问题p-1；分支定界法的输入为：满足卡罗需-库恩-塔克条件的松弛问题最优解x
relax
，松弛问题最优目标函数值z
relax
，0-1任意值ε；初始化k＝0，l＝0，u＝z
relax
；
[0169]
步骤4.2，从最优解x
relax
中任意选择一个不符合0-1约束条件的解xj，即：xj∈(0,1)；
[0170]
步骤4.3，若0≤xj《ε成立，则将约束条件xj＝0加到问题p-1中，形成子问题ⅰ；否则，将约束条件xj＝1加到问题p-1中，形成子问题ⅱ，ε表示0到1内任意值；
[0171]
步骤4.4，k ，继续求出子问题ⅰ或子问题ⅱ的松弛问题解，记为xk，并将对应的最优目标函数值记为zk；
[0172]
步骤4.5,找出最优目标函数最大值u作为新的上界，即：
[0173]
u＝max{zk′
|k
′
＝1,2,
…
,k},xk′
∈[0,1]；
[0174]
步骤4.6,再从符合0-1条件的分支中，找出目标函数最大值l作为新的下界，即：
[0175]
l＝max{zk′
|k
′
＝1,2,
…
,k},xk′
∈{0,1}；
[0176]
步骤4.7,若zk′
《l，剪掉相应分支,否则zk′
＞l，且不符合0-1条件,返回步骤4.3；
[0177]
步骤4.8,否则就意味着所有分支的最优目标函数值等于下界：zk′
＝l,将zk′
赋值z0-1
，将xk′
赋值给x
0-1
，并作为问题p-1的最优解，其中，x
0-1
表示符合0-1约束条件的最优解。
[0178]
本发明还提供一种处理设备，包括至少一个处理器，以及与所述处理器通信连接的至少一个存储器，其中：所述存储器存储有可被处理器执行的程序指令，所述处理器调用所述程序指令能够执行上述的方法。
[0179]
本发明还提供一种计算机可读存储介质，所述计算机可读存储介质存储计算机指令，所述计算机指令使所述计算机执行上述的方法。
[0180]
以上实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。