粒子群算法代码.docVIP

粒子群算法C++代码 一、问题重述 某一灾区有N名受灾群众，现有一批救灾物资要发放给这些受灾者。物资共有M种，每种物资的数量有限；各受灾者的灾情不同，对每种物资的急需程度和需求量不同。 (1)你作为一名物资分配者，请制定分配原则并给出合理的分配方法。 (2)试给出一个符合题意的数值算例。 二、模型假设 1、对于每一种物资，它的数量一定，且在分配过程中是以份为单位的； 2、灾民对获得一急需度可用某区间内的整型值表示； 3、灾民对得到一份某一种物资的满意度等价于获得该份该种物资的急需度。 4、灾民对某一物资的需求量用一整数来表示； 5、对于每种物资的分配是相对独立的，即只要每种物资的分配达到最优，则总体便最优。 三、符号定义 论文中使用的符号 对应在程序中用的符号 表示意义 N Dim 待分配人数 n[ij] X[ij] 第i个灾民实际分得的第j种物资的量 m[ij] Xup[ij] 第i个灾民对j种物资的需求量 X[ij] Value 第i个灾民对j种物资的急需度 fj（n[ij]） GetFit[j] 所有灾民对所获得的第j种物资的满意度函数 f（n[ij]） 所有灾民对所获得的所有物资的满意度函数 N[j] GdsCnt 第j种物资的数量 X 表示第i个人得到n[ij]份第j种物资的满意度 M 物资种数 T? ?? ??? PNum 在粒子群算法中粒子的个数 Pbest pso算法中，粒子的个体极值 Gbest pso算法中，粒子群的全局极值 X[ij][t]? ? 第ｉ个人分到第ｔ份第ｊ种物资的的满意度 #include stdafx.h #include math.h #include time.h #include iostream #include fstream using namespace std; int c1=2;? ?? ?? ?? ?? ?? ?? ?? ?? ?? ?? ?//加速因子 int c2=2;? ?? ?? ?? ?? ?? ?? ?? ?? ?? ?? ?//加速因子 double w=1;? ?? ?? ?? ?? ?? ?? ?? ?? ?? ?? ?//惯性权重 double Wmax=1;? ?? ?? ?? ?? ?? ?? ? //最大惯性权重 double Wmin=0.6;? ?? ?? ?? ?? ?? ? //最小惯性权重 int Kmax=110;? ?? ?? ?? ?? ?? ?? ?? ? //迭代次数 int GdsCnt;? ?? ?? ?? ?? ?? ?? ?? ?? ?? ?//物资总数 int const Dim=10;? ?? ?? ?? ?? ?? ???//粒子维数 int const PNum=50;? ?? ?? ?? ?? ???//粒子个数 int GBIndex=0;? ?? ?? ?? ?? ?? ?? ?? ?//最优粒子索引 double a=0.6;? ?? ?? ?? ?? ?? ?? ?? ???//适应度调整因子 double b=0.5;? ?? ?? ?? ?? ?? ?? ?? ???//适应度调整因子 int Xup[Dim];? ?? ?? ?? ?? ?? ?? ?? ???//粒子位置上界数组 int Xdown[Dim]={0};? ?? ?? ?? ?? ?//粒子位置下界数组 int Value[Dim];? ?? ?? ?? ?? ?? ?? ?? ?//初始急需度数组 int Vmax[Dim];? ?? ?? ?? ?? ?? ?? ?? ?//最大速度数组 class PARTICLE;? ?? ?? ?? ?? ?? ?? ???//申明粒子节点 void Check(PARTICLE,int);? ?? ?//约束函数 void Input(ifstream);? ?? ?? ?? ???//输入变量 void Initial();? ?? ?? ?? ?? ?? ?? ?? ?? ?//初始化相关变量 double GetFit(PARTICLE);? ?? ? //计算适应度 void CalculateFit();? ?? ?? ?? ?? ?? ?? ?? ?//计算适应度 void BirdsFly();? ?? ?? ?? ?? ?? ?? ?? ?//粒子飞翔 void Run(ofstream,int=2000); //运行函数 //微粒类 class PARTICLE { public: ? ???int X[Dim];? ?? ?? ?? ?? ?? ?? ? //微粒的坐标数组 ? ???int XBest[Dim];? ?? ?? ?? ?? ?? ?