1 简介

    花朵授粉算法( Flower Pollination Algorithm，FPA)是由英国剑桥大学学者Yang 于2012年提出的，其基本思想来源于对自然界花朵自花授粉、异花授粉的模拟，是一种新的元启发式群智能随机优化技术 。

     据统计，目前在自然界中被人类发现的植物大约有 37 万 种，显花植物大约占有 20 万种，而其中 80% 的植物依靠生物授粉繁衍后代。花植物已经进化了大约 1． 25 亿年，在演化过程中，花朵授粉在花植物繁衍过程中承担着举足轻重的作用，对于花植物，如果没有花，很难想象植物世界是个什么样子。人类的发展和生存与跟花植物也是息息相关的，例如人们生活中吃的苹果等水果都是花朵授粉的结果。花朵授粉是通过花粉的传播来实现，而花粉的传播主要是靠昆虫及动物来完成。在实际授粉过程中，一些花朵仅仅吸引和依靠一种特定的昆虫来成功授粉，即一些花朵与传粉者之间形成了一种非常特别的花—传粉者伙伴关系。授粉形式主要有非生物和生物两种，大概 90% 的显花植物是属于生物授粉植物，即花粉主要是通过动物和昆虫来传播而实现繁衍后代。大约 10% 的显花植物是属于非生物授粉植物，不需要任何传粉者来传播花粉，而是通过自然风或者扩散途径来完成传粉，以实现子代的繁殖。在依赖生物传粉的显花植物中大约有 85% 的植物是由蜜蜂完成传粉的，蜜蜂在实际传粉中可能只对一些特定花植物进行传粉，而忽视其他种类的花植物，这样以便以最小的代价获得最大的收益。同时对于一些花植物而言，也获得更多的传粉机会，繁衍更多的后代。根据显花植物的授粉对象不同，可分为异花授粉和自花授粉两种。在一般情况下，异花授粉是两性花，一般一朵花的雌蕊接受的花粉是另一朵花的雄蕊的花粉，这就是所谓的异花授粉。自花授粉是显花植物成熟的花粉粒传到同一朵花的柱头上或同一种显花植物的不同花之间进行传粉，并能正常地受精结实的过程。由于传粉者鸟、蜜蜂等能飞行很长的距离，故异花授粉可以发生在远距离的地方，这种方式称为全局授粉。另外，鸟和蜜蜂具有 Levy 飞行的行为，跳或飞行的步长服从 Levy飞行分布。而自花授粉称为局部授粉。

     花朵授粉算法是模拟自然界中显花植物花朵传粉的过程，该算法的理想条件假设如下:

a) 生物异花授粉是带花粉的传粉者通过 Levy 飞行进行的全局授粉过程。

b) 非生物自花授粉是局部授粉过程。

c) 花的常性可以被认为是繁衍概率，繁衍概率与参与的两朵花的相似性成比例关系。

d) 转换概率 p∈［0，1］控制全局授粉与局部授粉之间的转换，由于物理上的邻近性和风等其他因素的影响，局部授粉在整个授粉活动中是一个非常重要的部分 p。

     然而，在现实的自然界中，每一棵显花植物可以开好多朵花，每朵花产生数百万甚至数十亿的花粉配子。但是，为了把问题简单化，假设每棵显花植物仅仅只开一朵花，且每朵花仅产生一个花粉配子。因此，问题经过简化后，意味着一朵花或一个配子就对应于优化问题中的一个解。

【智能优化算法】基于花朵授粉算法求解单标目优化问题附matlab代码_无人机

2 部分代码 


          
          

           
           %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  • 1.

          
          

           
           %%%Flower Pollination Algorithm for Multimodal Optimization (MFPA)
  • 1.

          
          

           
           %Jorge G醠vez, Erik Cuevas and Omar Avalos
  • 1.

          
          

           
           %%This is the line to execute the code:
  • 1.

          
          

           
           %%[mem,bestSol,bestFit,optima,FunctionCalls]=FPA([50 0.25 500 2]);
  • 1.

          
          

           
           %FitFunc implements the function to be optimized
  • 1.

          
          

           
           %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  • 1.

          
          

           
           function [mem,bestSol,bestFit,optima,FunctionCalls]=FPA(para)
  • 1.

          
          

           
           % Default parameters
  • 1.

          
          

           
           if nargin<1,
  • 1.

          
          

           
           para=[50 0.25 500];
  • 1.

          
          

           
           end
  • 1.
  • 1.

          
          

           
           n=para(1);           % Population size
  • 1.

          
          

           
           p=para(2);           % Probabibility switch
  • 1.

          
          

           
           N_iter=para (3);  % Number of iterations
  • 1.
  • 1.

          
          

           
           phase = 1; %First state
  • 1.

          
          

           
           phaseIte= [0.5,0.9,1.01]; %State vector
  • 1.
  • 1.

          
          

           
           %Deb Function
  • 1.

          
          

           
           d = 1;
  • 1.

          
          

           
           Lb = 0;
  • 1.

          
          

           
           Ub = 1;
  • 1.

          
          

           
           optima = [.1;.3;.5;.7;.9];
  • 1.
  • 1.
  • 1.

          
          

           
           % Initialize the population
  • 1.

          
          

           
           for i=1:n,
  • 1.

          
          

           
           Sol(i,:)=Lb+(Ub-Lb).*rand(1,d);
  • 1.

          
          

           
           Fitness(i)=fitFunc(Sol(i,:));  %%Evaluate fitness function
  • 1.

          
          

           
           end
  • 1.
  • 1.

          
          

           
           % Initialice the memory
  • 1.

          
          

           
           [mem,bestSol,bestFit,worstF] = memUpdate(Sol,Fitness, [], zeros(1,d), 100000000, 0, phase,d,Ub,Lb);
  • 1.
  • 1.

          
          

           
           S = Sol;
  • 1.
  • 1.

          
          

           
           FunctionCalls = 0;
  • 1.

          
          

           
           % Main Loop
  • 1.

          
          

           
           for ite = 1 : N_iter,
  • 1.

          
          

           
           %For each pollen gamete, modify each position acoording
  • 1.

          
          

           
           %to local or global pollination
  • 1.

          
          

           
           for i = 1 : n,
  • 1.

          
          

           
           % Switch probability
  • 1.

          
          

           
           if rand>p,
  • 1.
  • 1.

          
          

           
           L=Levy(d);
  • 1.

          
          

           
           dS=L.*(Sol(i,:)-bestSol);
  • 1.

          
          

           
           S(i,:)=Sol(i,:)+dS;
  • 1.

          
          

           
           S(i,:)=simplebounds(S(i,:),Lb,Ub);
  • 1.

          
          

           
           else
  • 1.

          
          

           
           epsilon=rand;
  • 1.

          
          

           
           % Find random flowers in the neighbourhood
  • 1.

          
          

           
           JK=randperm(n);
  • 1.

          
          

           
           % As they are random, the first two entries also random
  • 1.

          
          

           
           % If the flower are the same or similar species, then
  • 1.

          
          

           
           % they can be pollenated, otherwise, no action.
  • 1.

          
          

           
           % Formula: x_i^{t+1}+epsilon*(x_j^t-x_k^t)
  • 1.

          
          

           
           S(i,:)=S(i,:)+epsilon*(Sol(JK(1),:)-Sol(JK(2),:));
  • 1.

          
          

           
           % Check if the simple limits/bounds are OK
  • 1.

          
          

           
           S(i,:)=simplebounds(S(i,:),Lb,Ub);
  • 1.

          
          

           
           end
  • 1.

          
          

           
           Fitness(i)=fitFunc(S(i,:));
  • 1.

          
          

           
           end
  • 1.

          
          

           
           %Update the memory
  • 1.

          
          

           
           [mem,bestSol,bestFit,worstF] = memUpdate(S,Fitness,mem,bestSol,bestFit,worstF,phase,d,Ub,Lb);
  • 1.
  • 1.

          
          

           
           Sol = get_best_nest(S, mem, p);
  • 1.
  • 1.

          
          

           
           FunctionCalls = FunctionCalls + n;
  • 1.
  • 1.

          
          

           
           if ite/N_iter > phaseIte(phase)
  • 1.

          
          

           
           %Next evolutionary process stage
  • 1.

          
          

           
           phase = phase + 1;
  • 1.

          
          

           
           [m,~]=size(mem);
  • 1.

          
          

           
           %Depurate the memory for each stage
  • 1.

          
          

           
           mem = cleanMemory(mem);
  • 1.

          
          

           
           FunctionCalls = FunctionCalls + m;
  • 1.

          
          

           
           end
  • 1.

          
          

           
           end
  • 1.
  • 1.

          
          

           
           %Plot the solutions (mem) founded by the multimodal framework
  • 1.

          
          

           
           x = 0:.01:1;
  • 1.

          
          

           
           y = ((sin(5.*pi.*x)).^ 6);
  • 1.

          
          

           
           plot(x,y)
  • 1.

          
          

           
           hold on
  • 1.

          
          

           
           plot(mem(:,1),-mem(:,2),'r*');
  • 1.

3 仿真结果

【智能优化算法】基于花朵授粉算法求解单标目优化问题附matlab代码_计算技术_02

4 参考文献

[1]陈西成, and 刘曙. "基于改进花朵授粉算法的防空部署优化研究." 计算技术与自动化 038.003(2019):74-78.

博主简介：擅长智能优化算法、神经网络预测、信号处理、元胞自动机、图像处理、路径规划、无人机等多种领域的Matlab仿真，相关matlab代码问题可私信交流。

部分理论引用网络文献，若有侵权联系博主删除。

【智能优化算法】基于花朵授粉算法求解单标目优化问题附matlab代码_计算技术_03