1.一种LTE无线资源分配方法，其特征在于，包括以下步骤：(1)、初始化功率、RB总数目及用户总数，将系统的总发射功率分均分配到每一个RB上,(2)、根据基站接收到的反馈信息对RB进行进一步分配：(21)、设置参数T，N，rinitial、rend和m，其中T为迭代次数，N为烟花弹数，m为每次爆炸新产生的烟花弹数目，rinitial为初始的烟花弹爆炸参考半径，rend为最后一次迭代时烟花弹爆炸参考半径；

(22)、在M维搜索空间内随机初始化N个烟花弹，根据基站接收到的反馈计算每一个初始化烟花弹的适应度值和本次迭代过程中烟花弹爆炸的参考半径r(t)，其中，适应度值的计算公式为：其中：

xi是一个M维的实向量，表示第i个烟花弹的位置；f(xi)表示第i个烟花弹的适应度值；k表示用户编号；K为小区内的用户总数；rk为用户k获得的传输速率；rk,min为用户k需要满足的最小速率；pk为第k个用户分得的功率；

最大半径 初始化t＝1，

其中：r(t)为迭代次数为t时的参考爆炸半径；T为最大迭代次数，rinitial和rend为常量，rinitial为初始参考爆炸半径，rend为迭代终止时的参考爆炸半径；

(23)、对N个烟花弹进行爆炸操作，得到最佳的RB分配方案；

(3)基站根据(2)得到的最佳的RB分配方案，对每个RB上的功率进行重分配：(31)设置算法参数T，N，rinitial、rend，m，其中T为迭代次数，N为烟花弹数，m为每次爆炸新产生的烟花弹数目，rinitial为初始的烟花弹爆炸参考半径，rend为最后一次迭代时烟花弹爆炸的参考半径；

(32)在M维搜索空间内随机初始化N个烟花弹，计算N个烟花弹的适应度值f(xi)和本次迭代过程中烟花弹爆炸的参考半径r(t),适应度值 其中：

xi是一个M维向量，i＝1,...m，f(xi)表示第i个烟花弹的适应度值；k表示的是用户编号；K为小区内的用户总数；pk为第k个用户分得的功率, a＝-10,b＝0；

初始化时t＝1，

其中：r(t)为迭代次数为t时的参考爆炸半径，T为最大迭代次数，rinitial和rend为常量，rinitial为是初始参考爆炸半径，rend为迭代终止时的参考爆炸半径；

(33)对N个烟花弹进行爆炸操作，

(34)产生下一代的烟花弹，保证每一代的烟花弹总数相等(341)计算新的m个烟花弹的适应度值，从所有的N+m个烟花弹中挑选出适应度值最优的 个烟花弹；

(342)从剩下的 个烟花弹中随机选取 个烟花弹；

(35)若当前迭代数t＜最大迭代数T，则返回步骤(33)，否则迭代结束，得到每个RB上的最佳功率分配方案，(4)、循环重复执行(2)～(3)直到适应度值不再变化或者小于设置的最小精度为止。

2.根据权利要求1所述的一种LTE无线资源分配方法，其特征在于，步骤(23)包括以下步骤：(231)、根据N个烟花弹的适应度值，其各自爆炸产生的新烟花的数目为：其中：xi是一个M维的实向量，表示第i个烟花弹的位置；f(xi)表示第i个烟花弹的适应度值；mi为第i个烟花弹爆炸产生的新烟花弹的数目；ymin为N个烟花弹的适应度值中的最小值；round为一个四舍五入取整函数；

(232)根据N个烟花弹的适应度值,其各自产生的烟花弹的最大爆炸半径为：其中：ri(t)为迭代次数为t时第i个烟花弹的实际爆炸半径，r(t)为迭代次数为t时的参考爆炸半径；ymax为N个烟花弹的适应度值中的最大值，f(xi)表示的是第i个烟花弹的适应度值；

(233)计算新产生的m个烟花弹位置：

其中： 为第i个烟花弹爆炸产生的第k个烟花弹，设xi＝(x1,x2,...,xM)，rand(1,M)产生一个M维向量，每个元素在0到1之间取值，(2*rand(1,M)-1)则产生一个所有元素在(-1,1)取值的M维向量；ri(t)为当前迭代的爆炸半径；

若 的第j维元素超出了最大用户序号K，则利用

重新映射回限制以内, 为 中超出限制的元素，xmax＝K，xmin＝1，％为取余符号，round表示取整函数，(234)产生下一代的烟花弹，保证每一代的烟花弹总数相等，(235)若当前迭代数t＜最大迭代数T，则返回步骤(231)，否则迭代结束，得到最佳的RB分配方案。

3.根据权利要求2所述的一种LTE无线资源分配方法，其特征在于，步骤(234)包括以下步骤(2341)计算新的m个烟花弹的适应度值

其中：xi为新产生的m个烟花弹(位置)；f(xi)表示的是第i个烟花弹的适应度值；k表示用户编号；K为小区内的用户总数；pk为第k个用户分得的功率；

其中a＝-10,b＝0；

再从所有的(N+m)个烟花弹中挑选出适应度值最优的 个烟花弹；

(2342)从剩下的 个烟花弹中随机选取 个烟花弹。

4.根据权利要求1所述的一种LTE无线资源分配方法，其特征在于，步骤(33)包括以下步骤：(331)根据当前N个烟花弹的适应度值计算各个烟花弹的各自爆炸产生的新烟花的数目为：其中：xi是一个M维的实向量，表示第i个烟花弹的位置；f(xi)表示第i个烟花弹的适应度值；mi即为第i个烟花弹爆炸产生的新烟花弹的数目；ymin为N个烟花弹的适应度值中的最小值；round函数为一个四舍五入取整函数；

(332)根据当前N个烟花弹的适应度值计算各个烟花弹各自产生的烟花弹的最大爆炸半径为：其中：

ri(t)为迭代次数为t时第i个烟花弹的实际爆炸半径，r(t)为迭代次数为t时的参考爆炸半径；ymax为N个烟花弹的适应度值中的最大值，f(xi)表示第i个烟花弹的适应度值；

(333)新产生的m个烟花弹位置：

其中：

为第i个烟花弹爆炸产生的第k个烟花弹，设xi＝(x1,x2,...,xM)，rand(1,M)产生一个M维向量，每个元素在0到1之间取值，(2*rand(1,M)-1)则产生一个所有元素在(-1,1)取值的M维向量；ri(t)为当前迭代的爆炸半径；

若 的第j维元素超出了限制每个RB所能分得的最大功率Pmax，则利用重新映射回限制以内, 为 中超出限制的元素，xmax＝Pmax，xmin＝0。