基于多目标决策的LEO卫星网络多业务路由优化方法，包括：卫星s接受到数据传输要求，目的节点为d；获取拓扑结构时间片，业务需求时延为a1，带宽为a2，误包率为a3；筛选可行链路集合V，筛选后得到V′，根据业务分类建立链路评价矩阵Α；判断理想链路v*是否存在，如果存在求解结束；如果理想链路v*不存在，求最接近理想链路的路径vi，作为下一个通信链路。针对当前业务和实时的链路状态为业务选择合适的路径，保证卫星网络资源整体利用率。
1.一种基于多目标决策的LEO卫星网络多业务路由优化方法，其特征在于，包括：首先根据卫星网络的拓扑结构进行时间片划分，针对不同的业务分类计算权值系数，对可行链路V＝{v1,v2,…vp}运用优选法对链路集合进行筛选，淘汰一些处于劣势的链路方案，得到筛选后的链路集合V'＝{v1,v2,…vq}，q≤p，路由过程中，对可行链路进行针对不同业务的矩阵评估，然后求出目标函数的理想链路，在路径集中找出最接近理想链路的路径从而得到问题的帕累托最优解；
可行链路域V′中，链路v*的各项属性满足Fi(v*)≤Fi(v)或Fi(v*)≥Fi(v)，其中 则称v*为理想链路：对于A类业务，有链路 分别都满足时延 剩余带宽 和误包
率 的要求，则称 为多目标最优链路，最优链路构成链路集合，链路集存在一个理想链路 满足同理，B类和C类业务分别有 和 所得到的这个理想值可能是存在着的某条
链路，也可能是在某几条链路之间的虚拟路径，在约束条件下，采用动态规划方法得到所需理想链路。
2.根据权利要求1所述基于多目标决策的LEO卫星网络多业务路由优化方法，其特征在于，本申请把业务划分成三种：A类是实时业务；B类是允许一定时延业务，C类是可靠性敏感业务，在此基础上对业务的优先级进一步规定，使A类业务高于B类业务，B类业务高于C类业务。
3.根据权利要求1所述基于多目标决策的LEO卫星网络多业务路由优化方法，其特征在于，将卫星运行周期T分成n个时间片,[t0＝0,t1],[t1,t2],[t2,t3]…[tn-1,tn＝T]；在每个时间片内，拓扑结构不变，且链路的切换和网络拓扑的变化只在时间点t0,t1,…,tn时刻发生，用G(V,E)表示卫星网络拓扑结构模型；其中，V＝M×N表示在星座中共分布于M条卫星轨道，每条轨道有N颗卫星；E代表卫星之间的星间链路。
4.根据权利要求1所述基于多目标决策的LEO卫星网络多业务路由优化方法，其特征在于，用 来分别表示链路Ek的n个属性值,即时延、剩余带宽和误包率；链路的总时延为：其中，用ws,d表示网络中所有可能的源、目的节点对，则Q(ws,d)＝{E1,E2,...,Ek,...,EK}代表一条路径序列，是SD节点之间通过K条链路连接；为了表示单个链路上的流量，令链路包含函数为：式中，若路径经过链路Ek则 取1，反之取0；用C表示每个节点向其他节点发送的数据包数量，则某个路径上的流量计算公式：剩余带宽为：
5.根据权利要求1所述基于多目标决策的LEO卫星网络多业务路由优化方法，其特征在于，获得v'是通过：其中，属性函数集F＝{f1,f2,...,fn}中fj用于评价每条链路的第j个链路属性，即aij＝fj(vi)，i＝1,2,…,m，j＝1,2,…,n； 表示第i类业务第j个属性的理想值，appr表示求逼近；不同QoS业务对于链路属性的要求强度不同，即wk＝(w1,w2,...,wn)表示第k类业务对于n个属性的要求程度，满足 wj表示某类业务的第j个链路属性权值，Vi属于链路集。
6.根据权利要求1所述基于多目标决策的LEO卫星网络多业务路由优化方法，其特征在于，链路评价矩阵
基于业务特征的权重是指把第1个属性时延的权w1和第2个属性剩余带宽的权w2之比记为α12，第2个属性剩余带宽的权w2和第3个属性误包率的权w3之比记为α23，以此类推，构成决策矩阵。
7.根据权利要求6所述基于多目标决策的LEO卫星网络多业务路由优化方法，其特征在于，(A-nI)w＝0，式中I是单位矩阵，若矩阵A中的值估计准确，上式等于0，求出权值系数w，若估计不够准确，则A中元素的摄动代表本征值的摄动，于是有Aw＝λmaxw
式中λmax是矩阵A的最大本征值，根据该式求得本征向量即权向量w＝[w1,w2,…,wn]T。
8.根据权利要求7所述基于多目标决策的LEO卫星网络多业务路由优化方法，其特征在于，判定矩阵A是否被接受： 若比率CI>0.1，各元素的估计一致性太差，应重新估计权向量w；若CI<0.1，各元素的估计基本一致，可用上式求得权向量w。
9.根据权利要求1所述基于多目标决策的LEO卫星网络多业务路由优化方法，其特征在于，实际通信链路的性能与理想链路性能的差异是指两条链路属性向量的加权欧氏距离；
实际通信链路和理想通信链路的性能分别用n维向量vi＝(vi1,vi2,…,vin)和表示，则它们之间的差距为：其中w1,w2,…,wn即为n个目标函数的权值，满足 且wi∈[0,1](i＝1,2,…,n)。
基于多目标决策的LEO卫星网络多业务路由优化方法
技术领域
[0001] 本发明属于LEO卫星网络领域，具体说是一种基于多目标决策的LEO卫星网络多业务路由优化方法。
背景技术
[0002] 随着卫星网络技术的不断发展，低轨卫星系统已能较好地实现全球移动通信。具有星间链路的LEO卫星网络能实现全球覆盖，与GEO卫星网络相比，它能够有效降低传输时延，减少卫星对地面节点的依赖，并能够更好地支持地面移动终端。但LEO卫星网络不同于一般的地面网络，它具有高误码率、长时延等空间通信的特点，与此同时，卫星网络业务类型不同，其对端到端传输时延、传输带宽等服务的需求也有所不同，因此卫星网络中不仅要满足不同业务传输的QoS参数要求，而且还需尽可能地提高网络传输效率，充分利用网络资源。而针对QoS所提出的路由，无论是按需路由还是流量分配路由，大多都是考虑某一种或两种链路属性来决定，从而忽略其他约束条件，这样容易导致网络的局部负载过大。因此，路由算法需要在兼顾多约束条件的情况下,尽量平衡地利用网络资源。
[0003] 目前，有关LEO卫星网络的路由算法中考虑链路状态特征的算法有：一种应用于节点的精确负载均衡explicit load balancing,ELB策略，它根据下一跳链路的时延，当节点出现链路数据拥塞时，发送信号给邻居节点，邻节点选择次优路径，从而减少网络拥塞。一种受限最短路径优先constraints shortest path first,CSPF算法,这是一种改进的最短路径优先算法，它为了避免网络或节点拥塞，将链路带宽的反比定义链路权重，根据业务的特定要求，在链路状态数据库的基础上，得到最终最短路径。多路径QoS路由算法multi-path QoS routing，MPQR是当卫星收到传输请求时，计算同时满足时延和带宽限制的最优路径。
[0004] 此外，考虑QoS业务分类的算法有：根据业务分类，一种多服务按需路由协议multiservice on-demand routing,MOR，它单独对各类服务流量进行路由。一种多业务类QoS路由算法multi-class QoS routing,MQoSR，该算法根据时延和带宽将业务分为两类，利用相对空闲链路来减少链路拥塞。由于这些算法有的只考虑链路状态信息，有的只考虑到业务分类，有针对当前业务和实时的链路状态为业务选择合适的路径，这样很难保证卫星网络资源整体利用率。
发明内容
[0005] 针对现有技术存在的上述缺点和不足，本发明提供了一种基于多目标决策的LEO卫星网络多业务路由优化方法，针对当前业务和实时的链路状态为业务选择合适的路径，保证卫星网络资源整体利用率。
[0006] 为实现上述目的，本发明提供了一种基于多目标决策的LEO卫星网络多业务路由优化方法，包括：
[0007] S1：卫星s接受到数据传输要求，目的节点为d；
[0008] S2：获取拓扑结构时间片，业务需求时延为a1，带宽为a2，误包率为a3；
[0009] S3：筛选可行链路集合V，筛选后得到V′，根据业务分类建立链路评价矩阵Α；
[0010] S4：判断理想链路v*是否存在，如果存在求解结束；
[0011] S5：如果理想链路v*不存在，求最接近理想链路的路径vi，作为下一个通信链路。
[0012] 进一步的，在进行步骤S1之前把业务划分成三种：A类是实时业务；B类是允许一定时延业务；C类是可靠性敏感业务；在此基础上对业务的优先级进一步规定，使A类业务高于B类业务，B类业务高于C类业务。
[0013] 进一步的，将卫星运行周期T分成n个时间片,[t0＝0,t1],[t1,t2],[t2,t3]…[tn-1,tn＝T]；在每个时间片内，拓扑结构不变，且链路的切换和网络拓扑的变化只在时间点t0,t1,…,tn时刻发生，用G(V,E)表示卫星网络拓扑结构模型；其中，V＝M×N表示在星座中共分布于M条卫星轨道，每条轨道有N颗卫星；E代表卫星之间的星间链路。
[0014] 进一步的，用 来分别表示链路Ek的n个属性值,即时延、剩余带宽和误包率；链路的总时延为：
[0015]
[0016] 其中，用ws,d表示网络中所有可能的源、目的节点对，则Q(ws,d)＝{E1,E2,…,Ek,…,EK}代表一条路径序列，是SD节点之间通过K条链路连接；为了表示单个链路上的流量，令链路包含函数为：
[0017]
[0018] 式中，若路径经过链路Ek则 取1，反之取0；用C表示每个节点向其他节点发送的数据包数量，则某个路径上的流量计算公式：
[0019]
[0020] 剩余带宽为：
[0021]
[0022] 误包率为：
[0023]
[0024] 进一步的，获得v'是通过：
[0025]
[0026] 其中，属性函数集F＝{f1,f2,...,fn}中fj用于评价每条链路的第j个链路属性，即aij＝fj(vi)，i＝1,2,…,m，j＝1,2,…,n； 表示第i类业务第j个属性的理想值，appr表示求逼近；不同QoS业务对于链路属性的要求强度不同，即wk＝(w1,w2,...,wn)表示第k类业务对于n个属性的要求程度，满足 wj表示某类业务的第j个链路属性权值，Vi属于链路集。
[0027] 进一步的，
[0028]
[0029] 基于业务特征的权重是指把第1个属性时延的权w1和第2个属性剩余带宽的权w2之比记为α12，第2个属性剩余带宽的权w2和第3个属性误包率的权w3之比记为α23，以此类推，构成决策矩阵。
[0030] 更进一步的：(Α-nI)w＝0，式中I是单位矩阵，若矩阵Α中的值估计准确，上式等于0，求出权值系数w，若估计不够准确，则Α中元素的摄动代表本征值的摄动，于是有[0031] Αw＝λmaxw
[0032] 式中λmax是矩阵Α的最大本征值，根据该式求得本征向量即权向量w＝[w1,w2,…,wn]T。
[0033] 更进一步的，判定矩阵Α是否被接受： 若比率CI>0.1，各元素的估计一致性太差，应重新估计；若CI<0.1，各元素的估计基本一致，可用上式求得w。
[0034] 作为更进一步的，可行链路域V′中，链路v*的各项属性满足Fi(v*)≤Fi(v)或Fi(v*)≥Fi(v)，其中 则称v*为理想链路：
[0035] 对于A类业务，有链路 分别都满足时延 剩余带宽 和误包率 的要求，则称 为多目标最优链路，最优链路构成链路集合，链路集存在一个理想链路 满足
[0036]
[0037] 同理，B类和C类业务分别有 和
[0038] 作为更进一步的，实际通信链路的性能与理想链路性能的差异是指两条链路属性向量的加权欧氏距离；
[0039] 实际通信链路和理想通信链路的性能分别用n维向量vi＝(vi1,vi2,…,vin)和表示，则它们之间的差距为：
[0040]
[0041] 其中w1,w2,…,wn即为n个目标函数的权值，满足 且wi∈[0,1](i＝1,2,…,n)。
[0042] 本发明由于采用以上技术方案，能够取得如下的技术效果：不仅保证了不同业务的不同要求，而且在吞吐量和负载分布指数等方面性能均有较明显的提升，能够有效地均衡卫星网络负载。
附图说明
[0043] 本发明共有附图1幅：
[0044] 图1为卫星网络拓扑结构模型示意图。
具体实施方式
[0045] 下面通过实施例，并结合附图，对本发明的技术方案作进一步的具体说明。
[0046] 实施例1
[0047] 一种基于多目标决策的LEO卫星网络多业务路由优化方法，根据数据语音视频等多样化的业务的QoS要求，把业务划分成三种：A类是实时业务，即时延敏感业务，如对时间要求高的指令语音等；B类是允许一定时延业务，即带宽敏感业务，例如对地观测业务；C类是可靠性敏感业务，主要体现在对误包率要求较为苛刻的业务上；在此基础上对业务的优先级进一步规定，使A类业务高于B类业务，B类业务高于C类业务；具体包括如下步骤：
[0048] S1：卫星s接受到数据传输要求，目的节点为d；
[0049] S2：获取拓扑结构时间片，业务需求时延为a1，带宽为a2，误包率为a3；将卫星运行周期T分成n个时间片,[t0＝0,t1],[t1,t2],[t2,t3]…[tn-1,tn＝T]；在每个时间片内，拓扑结构不变，且链路的切换和网络拓扑的变化只在时间点t0,t1,…,tn时刻发生，用G(V,E)表示卫星网络拓扑结构模型；其中，V＝M×N表示在星座中共分布于M条卫星轨道，每条轨道有N颗卫星；E代表卫星之间的星间链路(ISLs)；
[0050] S3：筛选可行链路集合V，筛选后得到V′，根据业务分类建立链路评价矩阵Α；
[0051] S4：判断理想链路v*是否存在，如果存在求解结束；
[0052] S5：如果理想链路v*不存在，求最接近理想链路的路径vi，作为下一个通信链路。
[0053] 实施例2
[0054] 与实施例1具有相同的技术方案，更为具体的是，用 来分别表示链路Ek的n个属性值,即时延、剩余带宽和误包率；链路的总时延为：
[0055]
[0056] 其中，用ws,d表示网络中所有可能的源、目的节点SD对，则Q(ws,d)＝{E1,E2,…,Ek,…,EK}代表一条路径序列，是SD节点之间通过K条链路连接；为了表示单个链路上的流量，令链路包含函数为：
[0057]
[0058] 式中，若路径经过链路Ek则 取1，反之取0；用C表示每个节点向其他节点发送的数据包数量，则某个路径上的流量计算公式：
[0059]
[0060] 剩余带宽为：
[0061]
[0062] 误包率为：
[0063]
[0064] 实施例3
[0065] 作为实施例1或2的补充，多业务路由的多目标决策问题，是指在可行链路域中，根据当前链路状态以及业务QoS要求，如何获得最优可行链路域v'是通过：
[0066]
[0067] 其中，属性函数集F＝{f1,f2,...,fn}中fj用于评价每条链路的第j个链路属性，即aij＝fj(vi)，i＝1,2,…,m，j＝1,2,…,n； 表示第i类业务第j个属性的理想值，appr表示求逼近；不同QoS业务对于链路属性的要求强度不同，即wk＝(w1,w2,...,wn)表示第k类业务对于n个属性的要求程度，满足 wj表示某类业务的第j个链路属性权值，Vi属于链路集。
[0068] 实施例4
[0069] 作为实施例3的补充说明，基于业务特征的权值计算，把n个属性的重要性成对比较，把第p个属性对第q个属性的相对重要性记为αpq，并认为，这就是属性p的权wp和属性q的权wq之比的近似值，αpq＝wp/wq，n个目标成对比较的结果形成矩阵A。
[0070]
[0071] 基于业务特征的权重是指把第1个属性时延的权w1和第2个属性剩余带宽的权w2之比记为α12，第2个属性剩余带宽的权w2和第3个属性误包率的权w3之比记为α23，以此类推，构成决策矩阵。
[0072] 作为进一步补充的：(Α-nI)w＝0，式中I是单位矩阵，若矩阵Α中的值估计准确，上式等于0，求出权值系数w，若估计不够准确，则Α中元素的摄动代表本征值的摄动，于是有
[0073] Αw＝λmaxw
[0074] 式中λmax是矩阵Α的最大本征值，根据该式求得本征向量即权向量w＝[w1,w2,…,wn]T。
[0075] 作为进一步补充的，为了判定矩阵Α在此方法中的科学性，引入一致性比率consistence rate,CR的概念，用一致性指标consistence index,CI与随机指标random index,RI的比值来表示，判定矩阵Α是否被接受： 对于阶数为n的矩阵对应的RI值如下：
[0076]
[0077] 若比率CR>0.1，说明各元素αpq的估计一致性太差，应重新估计。若CR<0.1，可认为αpq的估计基本一致，可用上式求得w。
[0078] 实施例5
[0079] 对上述实施例进行进一步补充：在路由选择的开始阶段，首先根据卫星网络的拓扑结构进行时间片划分，针对不同的业务分类计算权值系数，对可行链路V＝{v1,v2,…vp}运用优选法对链路集合进行筛选，淘汰一些处于劣势的链路方案，得到筛选后的链路集合V'＝{v1,v2,…vq}，q≤p，路由过程中，对可行链路进行针对不同业务的矩阵评估，然后求出目标函数的理想链路，在路径集中找出最接近理想链路的路径从而得到问题的帕累托最优解。
[0080] 对于LEO卫星网络多业务路由决策问题的求解基于以下定义：可行链路域V′中，链* * * *路v的各项属性满足Fi(v)≤Fi(v)或Fi(v)≥Fi(v)，其中 则称v为理想链路：
[0081] 对于A类业务，有链路 分别都满足时延 剩余带宽 和误包率 的要求，则称 为多目标最优链路，最优链路构成链路集合，链路集存在一个理想链路 满足
[0082]
[0083] 同理，B类和C类业务分别有 和 所得到的这个理想值可能是存在着的某条链路，也可能是在某几条链路之间的虚拟路径，在约束条件下，可采用动态规划等方法得到所需理想链路。
[0084] 作为更进一步补充的，实际通信链路的性能与理想链路性能的差异是指两条链路属性向量的加权欧氏距离；
[0085] 实际通信链路和理想通信链路的性能分别用n维向量vi＝(vi1,vi2,…,vin)和表示，则它们之间的差距为：
[0086] 其中w1,w2,…,wn即为n个目标函数的权值，满足 且wi∈[0,1](i＝1,2,…,n)。
[0087] 上述方法中约束条件是确保某两个节点之间有且只有一条路径，表示为：
[0088] 本申请的时间复杂度分析如下：由于针对每个节点需计算三个属性值，时间复杂度O(kn)，k＝3；同时还要计算链路的加权欧氏距离，即O(n2)；因而时间复杂度为O(n2+kn)，其中k是需要评价的节点属性数量，最终算法时间复杂度为O(n2)。
[0089] 以上所述，仅为本发明较佳的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明披露的技术范围内，根据本发明的技术方案及其发明构思加以等同替换或改变，都应涵盖在本发明的保护范围之内。
