CN107864506A - 基于萤火虫算法的Ka多波束卫星通信系统资源分配方法 - Google Patents

Abstract

本发明公开了一种基于萤火虫算法的Ka多波束卫星通信系统资源分配方法。该方法为：首先确定每个波束内的链接申请的载波分配状态信息；然后采用萤火虫算法动态调整链接申请的载波编码模式：收集链接申请的编码的参数并确定萤火虫的目标函数值；确定群体中萤火虫的相对亮度和吸引度，并确定目标函数值最大的萤火虫位置；更新萤火虫的控件位置向量；检验新产生的萤火虫位置是否满足限制条件，如果满足进入下一步，否则返回上一步；重复多次，得到链接申请编码的最佳编码模式；最后通过时隙封装将链接申请分配到载波上，完成该波束内的链路申请分配。本发明从时隙和载波两个维度，按MF‑TDMA卫星通信系统的约束条件，实现了多波束卫星系统资源高效分配。

Description

基于萤火虫算法的Ka多波束卫星通信系统资源分配方法

技术领域

本发明属于卫星通信技术领域，特别是一种基于萤火虫算法的Ka多波束卫星通信系统资源分配方法。

背景技术

在MF-TDMA卫星通信系统中，待分配的系统资源是二维的频域和时域资源。NCC 将频域和时域资源以载波和时隙的形式进行划分，把不同的载波上的不同时隙封装给不同的用户以满足其业务需求。与SCPC系统中的一维资源分配问题相比，MF-TDMA二维的资源分配策略更加灵活，但是也正是这种灵活性，增加了其资源分配的难度。此外，与SCPC技术相比，MF-TDMA技术需要多个载波在时间上都实现同步，增加了网络管理的复杂度。

MF-TDMA卫星通信系统中，当地面站所支持的链接有业务需要时，地面站会向 NCC发送链接申请。在一个帧时间内，NCC综合考虑所有接收到的链接申请，根据资源分配算法动态地将不同的链接申请，分配到不同的载波上的不同时隙中，然后将分配结果以突发时间计划(Burst Time Plan,BTP)的形式，周期性的下发给所有的地面站。地面站收到BTP后，对其进行解析，寻找属于自己的载波和时隙，然后在相应的载波和时隙上发送业务数据。与FDMA卫星通信系统不同，在MF-TDMA卫星通信系统中，资源分配的对象不再是每个地面站(用户)，而是每个地面站所支持的链接。因为在 MF-TDMA卫星通信系统中，同一个地面站下链接可以在一个帧周期内与不同的地面站的链接进行通信。而FDMA卫星通信系统无法做到这点，只能做到同一个地面站下的链接与另外一个地面站下的链接进行通信。

发明内容

本发明的目的在于提供一种基于萤火虫算法的Ka多波束卫星通信系统资源分配方法，通过固定链接申请编码模式，在充分利用时频资源和满足链接申请业务量的基础上，尽可能的减小链路申请发送端地面站的发射功率。

实现本发明目的的技术解决方案为：一种基于萤火虫算法的Ka多波束卫星通信系统资源分配方法，包括以下步骤：

步骤1，确定每个波束内的链接申请的载波分配状态信息；

步骤2，采用萤火虫算法，对链接申请的调制编码模式进行动态调整；

步骤3，通过时隙封装将链接申请分配到载波上，完成该波束内的链路申请分配。

进一步地，步骤1所述的确定每个波束内的链接申请的载波分配状态信息，包括：波束载波的总数目N，链接申请的总数目W，载波速率总和S_total，链接申请的总业务量D_total，所有载波的集合Φ，Φ＝{F₁,F₂,…,F_N}，其中F_j为第j条载波，并且按照载波速率的大小升序排列，即S₁≤S₂≤……S_N，S_j为第j条载波的载波速率；X为所有链接申请的集合，X＝{C₁,C₂,……,C_W}，其中C_i为第i个链接申请，并且按照链接申请业务量的大小升序排列，即D₁≤D₂≤……D_W，D_i为第i个链接申请的业务量；Y_j为第j条载波的时隙数目。

进一步地，步骤2所述的采用萤火虫算法，对链接申请的调制编码模式进行动态调整，具体步骤如下：

步骤2.1，根据链接申请的载波分配情况，初始化萤火虫算法的参数和种群中每个萤火虫的位置向量，并确定萤火虫的目标函数值；

步骤2.2，计算群体中萤火虫的相对亮度和吸引度，根据相对亮度决定萤火虫的转移方向，确定目标函数最大的萤火虫位置；

步骤2.3，更新萤火虫的空间位置向量；

步骤2.4，检验新产生的萤火虫位置是否满足限制条件，如果满足，进入步骤2.5，否则，返回步骤2.3；

步骤2.5，重复步骤2.2～步骤2.4N_g次，输出目标函数值的最优解，其中N_g为萤火虫方法的最大迭代次数。

进一步地，步骤2.1所述的根据链接申请的载波分配情况，初始化萤火虫算法的参数和种群中每个萤火虫的位置向量，并确定萤火虫的目标函数值，具体为：

(1)初始化萤火虫种群大小NF_p、萤火虫算法最大迭代次数NF_g、光强吸收系数γ、最大吸引度β₀、步长因子ε，其中NF_p∈[100,130]、NF_g∈[500,600]、γ∈[0.9,1.1]、β₀∈[0.9,1.1]、ε∈[0.3,0.4]；令gen＝0，初始化每个链接申请的最高阶调制编码模式

(2)初始化第gen次迭代时种群中萤火虫位置向量x(gen)＝[c(gen),p(gen)]，其中c(gen)为第gen次迭代时种群中萤火虫链接申请向量，p(gen)为第gen次迭代时种群中萤火虫编码分配向量，并且K为多波束系统中波束数，N 为多波束卫星通信系统中单个波束内的载波数目，W为多波束卫星通信系统中单个波束内的链接申请的数目，为时刻t波束k中链接申请w在载波n上分配载波指示变量，且满足公式(1)：

为时刻t波束k中链接申请w在载波n上调制编码分配指示变量，且满足公式(2)：

同时链接申请编码模式的选择必须满足预算方程(3)和(4)：

其中，[]表示一种运算，[x]＝10lg(x)；Mⁱ为第i个链路申请C_i的链路余量，Con 为大于0的常量；D_i为第i个链路申请C_i的业务量；b为波尔兹曼常数；

(3)确定第gen次迭代时萤火虫i的目标函数值f(x_i(gen))，如公式(5)：

其中，为t时刻波束k中链接申请w在载波n上发射功率，JFI(t)为时刻t用户的Jain公平指数，α为罚函数因子且α∈[10,12]，令萤火虫的最大荧光亮度 LF₀＝f(x_i(gen))。

进一步地，步骤2.2所述的计算群体中萤火虫的相对亮度和吸引度，根据相对亮度决定萤火虫的转移方向，确定目标函数最大的萤火虫位置计，具体为：

根据公式(6)计算群体中萤火虫的相对亮度LF：

LF＝LF₀×exp(-γD_i,j) (6)

其中，LF₀为最大相对亮度；

根据公式(7)计算吸引度β：

其中，β₀为最大吸引度

根据相对亮度决定萤火虫的移动方向，确定目标函数值最大的萤火虫位置 x_gbest(gen)，D_i,j为萤火虫i和j的空间距离。

进一步地，步骤2.3所述的更新萤火虫的空间位置向量，具体为：根据公式(8) 更新萤火虫的空间位置，对处于最佳位置的萤火虫进行随机扰动；

设定为第gen+1次迭代时萤火虫i的空间位置，则：

其中，为第gen次迭代时萤火虫i的空间位置；为第gen次迭代时萤火虫j的空间位置；ε为步长因子；ξ为[0,1]上服从均与分布的随机因子。

进一步地，步骤2.4中所述的检验新产生的萤火虫位置是否满足限制条件，具体为：检验新产生的萤火虫位置是否满足公式(1)、(2)，如果满足则进入步骤2.5，否则返回步骤2.3。

进一步地，步骤2.5中所述的输出目标函数的最优解，具体为：

令gen←gen+1，重复步骤2.2～步骤2.4NF_g次，输出x_gbest(gen)作为最优解，并且采用公式(9)确定波束内的链接申请模式F(x_gbest(gen))

其中，为t时刻波束k中链接申请w在载波n上分配载波指示变量；为t时刻波束k中链接申请w在载波n上分配发射功率指示变量。

进一步地，步骤3中所述的通过时隙封装将链接申请分配到载波上，具体为：在完成链接申请载波的分配和链接申请的编码模式分配之后，通过给定的启发式算法来实现实时隙封装，步骤如下：

步骤3.1，选择波束内空余时间最大的载波，如果存在空余时间相同的两个载波，则选择速率小的载波，为此载波上未分配时隙的链接申请分配时隙；

步骤3.2，若该载波上没有足够的时隙来容纳待分配的链接申请，则选择待分配队列的下一个链接申请；

步骤3.3，若选择的链接申请与载波上已分配的链接申请没有产生MF-TDMA卫星限制条件冲突，则分配成功，否则分配失败；

步骤3.4，重复步骤3.1-3.3，直至波束内的载波可用时隙均分配完成，结束该波束的时隙封装过程。

本发明与现有技术相比，其显著优点在于：(1)在卫星通信系统调制编码开始阶段，将所有链接申请编码模式定为最低阶的调制编码模式，通过动态调整编码，节省了地面站发射功率；(2)选择余量最大的链接申请可以充分利用现有的时频资源，尽量减小未利用的时频资源‘(3)载波和功率联合分配大大提高了多波束卫星系统的资源利用率。

附图说明

图1为本发明基于萤火虫算法的Ka多波束卫星通信系统资源分配方法的流程图。

图2为本发明中基于萤火虫方法动态调整分配流程图。

图3为本发明中链路申请的调制编码模式分配算法流程图。

图4为本发明中链路申请的时隙封装流程图。

具体实施方式

下面结合附图及具体实施例对本发明作进一步详细描述。

结合图1，本发明基于基于萤火虫算法的Ka多波束卫星通信资源分配方法，在充分利用在Ka频段下的多波束卫星通信系统时隙资源和满足链接申请业务申请量的基础上，在最开始将每个波束内的所有链接申请的编码模式设置为最低阶，为了尽可能的减少链路申请发送端地面站的发射功率，通过采用萤火虫算法通过动态调整链接申请的载波编码模式：收集链接申请的编码的参数并确定萤火虫的目标函数值；确定群体中萤火虫的相对亮度和吸引度，并确定目标函数值最大的萤火虫位置；更新萤火虫的控件位置向量；检验新产生的萤火虫位置是否满足限制条件，如果满足进入下一步，否则返回上一步；重复多次，得到链接申请编码的最佳编码模式；最后通过时隙封装将链接申请分配到载波上，完成该波束内的链接申请分配。具体包括以下步骤：

步骤1，确定每个波束内的链接申请的载波分配状态信息；

所述的确定每个波束内的链接申请的载波分配状态信息，包括：波束载波的总数目 N，链接申请的总数目W，载波速率总和S_total，链接申请的总业务量D_total，所有载波的集合Φ，Φ＝{F₁,F₂,…,F_N}，其中F_j为第j条载波，并且按照载波速率的大小升序排列，即S₁≤S₂≤……S_N，S_j为第j条载波的载波速率；X为所有链接申请的集合，X＝{C₁,C₂,……,C_W}，其中C_i为第i个链接申请，并且按照链接申请业务量的大小升序排列，即D₁≤D₂≤……D_W，D_i为第i个链接申请的业务量；Y_j为第j条载波的时隙数目。

步骤2，采用萤火虫算法，对链接申请的调制编码模式进行动态调整；

步骤2.1，根据链接申请的载波分配情况，初始化萤火虫算法的参数和种群中每个萤火虫的位置向量，并确定萤火虫的目标函数值，算法流程如图3所示；

(1)初始化萤火虫种群大小NF_p、萤火虫算法最大迭代次数NF_g、光强吸收系数γ、最大吸引度β₀、步长因子ε，其中NF_p∈[100,130]、NF_g∈[500,600]、γ∈[0.9,1.1]、β₀∈[0.9,1.1]、ε∈[0.3,0.4]；令gen＝0，初始化每个链接申请的最高阶调制编码模式

(2)结合图2，初始化第gen次迭代时种群中萤火虫位置向量 x(gen)＝[c(gen),p(gen)]，其中c(gen)为第gen次迭代时种群中萤火虫链接申请向量， p(gen)为第gen次迭代时种群中萤火虫编码分配向量，并且 K为多波束系统中波束数，N为多波束卫星通信系统中单个波束内的载波数目，W为多波束卫星通信系统中单个波束内的链接申请的数目，为时刻t波束k中链接申请 w在载波n上分配载波指示变量，且满足公式(1)：

为时刻t波束k中链接申请w在载波n上调制编码分配指示变量，且满足公式(2)：

同时链接申请编码模式的选择必须满足预算方程(3)和(4)：

其中，[]表示一种运算，[x]＝10lg(x)；Mⁱ为第i个链路申请C_i的链路余量，Con 为大于0的常量；D_i为第i个链路申请C_i的业务量；b为波尔兹曼常数；

(3)确定第gen次迭代时萤火虫i的目标函数值f(x_i(gen))，如公式(5)：

其中，为t时刻波束k中链接申请w在载波n上发射功率，JFI(t)为时刻t用户的Jain公平指数，α为罚函数因子且α∈[10,12]，令萤火虫的最大荧光亮度 LF₀＝f(x_i(gen))。

步骤2.2，计算群体中萤火虫的相对亮度和吸引度，根据相对亮度决定萤火虫的转移方向，确定目标函数最大的萤火虫位置；

根据公式(6)计算群体中萤火虫的相对亮度LF：

LF＝LF₀×exp(-γD_i,j) (6)

其中，LF₀为最大相对亮度；

根据公式(7)计算吸引度β：

其中，β₀为最大吸引度

根据相对亮度LF决定萤火虫的移动方向，确定目标函数值最大的萤火虫位置x_gbest(gen)，D_i,j为萤火虫i和萤火虫j的空间距离；

步骤2.3，更新萤火虫的空间位置向量，具体为：根据公式(8)更新萤火虫的空间位置，对处于最佳位置的萤火虫进行随机扰动；

设定为第gen+1次迭代时萤火虫i的空间位置，则：

其中，为第gen次迭代时萤火虫i的空间位置；为第gen次迭代时萤火虫j的空间位置；ε为步长因子；ξ为[0,1]上服从均与分布的随机因子；

步骤2.4，检验新产生的萤火虫位置是否满足限制条件，如果满足，进入步骤2.5，否则，返回步骤2.3；

所述检验新产生的萤火虫位置是否满足限制条件，包括检验新产生的萤火虫位置是否满足公式(10)～(11)，如果满足，进入步骤2.5，否则，返回步骤2.3；

步骤2.5，重复步骤2.2～步骤2.4N_g次，输出编码模式和地面站发射功率即目标函数值的最优解，N_g为返回萤火虫方法的最大迭代次数；具体为：

令gen←gen+1，重复步骤2.2～步骤2.4NF_g次，输出x_gbest(gen)作为最优解，并且采用公式(9)确定波束内的链接申请最佳编码模式F(x_gbest(gen))：

其中，为t时刻波束k中链接申请w在载波n上分配载波指示变量；为t时刻波束k中链接申请w在载波n上分配发射功率指示变量。

步骤3，通过时隙封装将链接申请分配到载波上，完成该波束内的链路申请分配。

如图4所示，所述的通过时隙封装将链接申请分配到载波上，具体为：在完成链接申请载波的分配和链接申请的编码模式分配之后，通过给定的启发式算法来实现实时隙封装，步骤如下：

步骤3.1，选择波束内空余时间最大的载波，如果存在空余时间相同的两个载波，则选择速率小的载波，为此载波上未分配时隙的链接申请分配时隙；

步骤3.2，若该载波上没有足够的时隙来容纳待分配的链接申请，则选择待分配队列的下一个链接申请；

步骤3.3，若选择的链接申请与载波上已分配的链接申请没有产生MF-TDMA卫星限制条件冲突，则分配成功，否则分配失败；

步骤3.4，重复步骤3.1-3.3，直至波束内的载波可用时隙均分配完成，结束该波束的时隙封装过程。

下面结合具体实施例对本发明作进一步说明。

实施实例1

结合图1，本发明基于萤火虫算法的Ka多波束卫星通信系统资源分配优化方法，各步骤中参数如下：

步骤1，确定每个波束内的链接申请的载波分配状态信息。

帧长T_frame为110ms，载波数目N为4，载波符号速率S_n为 64K,2×64K,4×64K,8×128K，每个载波时隙长度L_n为10,7,5.5,3，时隙头长度H_n为 128符号，时隙尾长度T_n为0.1，系统支持的调制编码模式为系统允许的最大误码率为e^-6，调制编码模式的解调时的门限比特信噪比为3,4.2,5.4,6.5,7.8,9.8，每种地面站的EIRP值为 76,68,64,60,48,43，链接申请的业务量服从均值为60K的指数分布，系统总带宽W＝500 MHz，波束数K＝10。

步骤2，采用萤火虫算法，对链接申请的调制编码模式进行动态调整。

基于萤火虫算法的Ka多波束系统中有K个波束，W个链接申请，N＝4个载波；结合图2中基于萤火虫方法的动态调整分配流程：首先，系统初始化萤火虫方法的参数、初始化种群中的每个萤火虫位置向量和确定萤火虫的目标函数值，NF_p＝110、 NF_g＝560、γ＝1、β₀＝1和ε＝0.35，令gen＝1，初始化初始化 x(gen)＝[c(gen),p(gen)]，满足公式(1)，满足公式(2)，采用公式(5)确定 f(x_i(gen))，令LF₀＝f(x_i(gen))，α＝11，萤火虫的最大荧光亮度LF₀＝1。然后，确定群体中萤火虫的相对亮度和吸引度，根据相对亮度决定萤火虫的移动方向，确定目标函数值最大的萤火虫位置，用公式(6)和公式(7)确定群体中萤火虫的相对亮度LF和吸引度β；根据相对亮度决定萤火虫的移动方向，确定目标函数值最大的萤火虫位置x_gbest(gen)；其次，更新萤火虫的空间位置向量，根据公式(8)更新萤火虫的空间位置。对处于最佳位置的萤火虫进行随机扰动；再次，检验新产生的萤火虫位置是否满足公式(1)、公式(2)和公式(5)，如果满足，进入下一步，否则返回上一步重新更新萤火虫的空间位置向量；最后令gen←gen+1，重复上述步骤NF_g次，输出x_gbest(gen)作为最优解，并且采用公式(9)确定基于编码调整Ka多波束系统中波束内的链接申请最佳编码和发射功率。

步骤3，通过时隙封装将链接申请分配到载波上，完成该波束内的链路申请分配。

综上所述，本发明基于萤火虫算法的Ka多波束卫星通信系统资源分配方法，充分利用了载波的时隙资源，优化了链接申请的编码模式，使链接申请占用更少的时隙，可以容纳更多的链接申请。

Claims (9)

1.一种基于萤火虫算法的Ka多波束卫星通信系统资源分配方法，其特征在于，包括以下步骤：

步骤1，确定每个波束内的链接申请的载波分配状态信息；

步骤2，采用萤火虫算法，对链接申请的调制编码模式进行动态调整；

步骤3，通过时隙封装将链接申请分配到载波上，完成该波束内的链路申请分配。

2.根据权利要求1所述的基于萤火虫算法的Ka多波束卫星通信系统资源分配方法，其特征在于，步骤1所述的确定每个波束内的链接申请的载波分配状态信息，包括：波束载波的总数目N，链接申请的总数目W，载波速率总和S_total，链接申请的总业务量D_total，所有载波的集合Φ，Φ＝{F₁,F₂,…,F_N}，其中F_j为第j条载波，并且按照载波速率的大小升序排列，即S₁≤S₂≤……S_N，S_j为第j条载波的载波速率；X为所有链接申请的集合，X＝{C₁,C₂,……,C_W}，其中C_i为第i个链接申请，并且按照链接申请业务量的大小升序排列，即D₁≤D₂≤……D_W，D_i为第i个链接申请的业务量；Y_j为第j条载波的时隙数目。

3.根据权利要求1所述的基于萤火虫算法的Ka多波束卫星通信系统资源分配方法，其特征在于，步骤2所述的采用萤火虫算法，对链接申请的调制编码模式进行动态调整，具体步骤如下：

步骤2.1，根据链接申请的载波分配情况，初始化萤火虫算法的参数和种群中每个萤火虫的位置向量，并确定萤火虫的目标函数值；

步骤2.2，计算群体中萤火虫的相对亮度和吸引度，根据相对亮度决定萤火虫的转移方向，确定目标函数最大的萤火虫位置；

步骤2.3，更新萤火虫的空间位置向量；

步骤2.4，检验新产生的萤火虫位置是否满足限制条件，如果满足，进入步骤2.5，否则，返回步骤2.3；

步骤2.5，重复步骤2.2～步骤2.4N_g次，输出目标函数值的最优解，其中N_g为萤火虫方法的最大迭代次数。

4.根据权利要求3所述的基于萤火虫算法的Ka多波束卫星通信系统资源分配方法，其特征在于，步骤2.1所述的根据链接申请的载波分配情况，初始化萤火虫算法的参数和种群中每个萤火虫的位置向量，并确定萤火虫的目标函数值，具体为：

(1)初始化萤火虫种群大小NF_p、萤火虫算法最大迭代次数NF_g、光强吸收系数γ、最大吸引度β₀、步长因子ε，其中NF_p∈[100,130]、NF_g∈[500,600]、γ∈[0.9,1.1]、β₀∈[0.9,1.1]、ε∈[0.3,0.4]；令gen＝0，初始化每个链接申请的最高阶调制编码模式

(2)初始化第gen次迭代时种群中萤火虫位置向量x(gen)＝[c(gen),p(gen)]，其中c(gen)为第gen次迭代时种群中萤火虫链接申请向量，p(gen)为第gen次迭代时种群中萤火虫编码分配向量，并且K为多波束系统中波束数，N为多波束卫星通信系统中单个波束内的载波数目，W为多波束卫星通信系统中单个波束内的链接申请的数目，为时刻t波束k中链接申请w在载波n上分配载波指示变量，且满足公式(1)：

为时刻t波束k中链接申请w在载波n上调制编码分配指示变量，且满足公式(2)：

同时链接申请编码模式的选择必须满足预算方程(3)和(4)：

其中，[]表示一种运算，[x]＝10lg(x)；Mⁱ为第i个链路申请C_i的链路余量，Con为大于0的常量；D_i为第i个链路申请C_i的业务量；b为波尔兹曼常数；

(3)确定第gen次迭代时萤火虫i的目标函数值f(x_i(gen))，如公式(5)：

其中，为t时刻波束k中链接申请w在载波n上发射功率，JFI(t)为时刻t用户的Jain公平指数，α为罚函数因子且α∈[10,12]，令萤火虫的最大荧光亮度LF₀＝f(x_i(gen))。

5.根据权利要求3所述的基于萤火虫算法的Ka多波束卫星通信系统资源分配方法，其特征在于，步骤2.2所述的计算群体中萤火虫的相对亮度和吸引度，根据相对亮度决定萤火虫的转移方向，确定目标函数最大的萤火虫位置计，具体为：

根据公式(6)计算群体中萤火虫的相对亮度LF：

LF＝LF₀×exp(-γD_i,j) (6)

其中，LF₀为最大相对亮度；

根据公式(7)计算吸引度β：

其中，β₀为最大吸引度

根据相对亮度决定萤火虫的移动方向，确定目标函数值最大的萤火虫位置x_gbest(gen)，D_i,j为萤火虫i和j的空间距离。

6.根据权利要求3所述的基于萤火虫算法的Ka多波束卫星通信系统资源分配方法，其特征在于，步骤2.3所述的更新萤火虫的空间位置向量，具体为：根据公式(8)更新萤火虫的空间位置，对处于最佳位置的萤火虫进行随机扰动；

设定为第gen+1次迭代时萤火虫i的空间位置，则：

其中，为第gen次迭代时萤火虫i的空间位置；为第gen次迭代时萤火虫j的空间位置；ε为步长因子；ξ为[0,1]上服从均与分布的随机因子。

7.根据权利要求3所述的基于萤火虫算法的Ka多波束卫星通信系统资源分配方法，其特征在于，步骤2.4中所述的检验新产生的萤火虫位置是否满足限制条件，具体为：检验新产生的萤火虫位置是否满足公式(1)、(2)，如果满足则进入步骤2.5，否则返回步骤2.3。

8.根据权利要求3所述的基于萤火虫算法的Ka多波束卫星通信系统资源分配方法，其特征在于，步骤2.5中所述的输出目标函数的最优解，具体为：

令gen←gen+1，重复步骤2.2～步骤2.4NF_g次，输出x_gbest(gen)作为最优解，并且采用公式(9)确定波束内的链接申请模式F(x_gbest(gen))

其中，为t时刻波束k中链接申请w在载波n上分配载波指示变量；为t时刻波束k中链接申请w在载波n上分配发射功率指示变量。

9.根据权利要求1或3所述的基于萤火虫算法的Ka多波束卫星通信系统资源分配方法，其特征在于，步骤3中所述的通过时隙封装将链接申请分配到载波上，具体为：在完成链接申请载波的分配和链接申请的编码模式分配之后，通过给定的启发式算法来实现实时隙封装，步骤如下：

步骤3.1，选择波束内空余时间最大的载波，如果存在空余时间相同的两个载波，则选择速率小的载波，为此载波上未分配时隙的链接申请分配时隙；

步骤3.2，若该载波上没有足够的时隙来容纳待分配的链接申请，则选择待分配队列的下一个链接申请；

步骤3.3，若选择的链接申请与载波上已分配的链接申请没有产生MF-TDMA卫星限制条件冲突，则分配成功，否则分配失败；

步骤3.4，重复步骤3.1-3.3，直至波束内的载波可用时隙均分配完成，结束该波束的时隙封装过程。