• 基于PSO粒子群算法的MPPT最大功率跟踪Simulink仿真,PSO采用S函数实现

    目录

    1.算法描述

    2.仿真效果预览

    3.MATLAB核心程序

    4.完整MATLAB

    1.算法描述

            MPPT控制器的全称是“最大功率点跟踪”（Maximum Power Point Tracking）太阳能控制器，是传统太阳能充放电控制器的升级换代产品。MPPT控制器能够实时侦测太阳能板的发电电压，并追踪最高电压电流值(VI)，使系统以最大功率输出对蓄电池充电。应用于太阳能光伏系统中，协调太阳能电池板、蓄电池、负载的工作，是光伏系统的大脑。

           最大功率点跟踪系统是一种通过调节电气模块的工作状态，使光伏板能够输出更多电能的电气系统能够将太阳能电池板发出的直流电有效地贮存在蓄电池中，可有效地解决常规电网不能覆盖的偏远地区及旅游地区的生活和工业用电，不产生环境污染。

            目前，光伏阵列的最大功率点跟踪(MPPT)技术，国内外已有了一定的研究，发展出各种控制方法常，常用的有一下几种：恒电压跟踪法(ConstantVoltageTracking简称CVT)、干扰观察法(PerturbationAndObservationmethod简称P&O)、增量电导法(IncrementalConductancemethod简称INC)、基于梯度变步长的电导增量法等等。(这些算法只能用在无遮挡的条件下)

    1）单峰值功率输出的MPPT的算法

           目前，在无遮挡条件下，光伏阵列的最大功率点跟踪(MPPT)的控制方法常用的有以下几种：
    l恒电压跟踪法(ConstantVoltageTracking简称CVT)
    l干扰观察法(PerturbationAndObservationmethod简称P&O)
    l增量电导法(IncrementalConductancemethod简称INC)
    l基于梯度变步长的电导增量法，等等。

    2）多峰值功率输出MPPT算法

           普通的最大功率跟踪算法，如扰动观测发和电导增量法在一片云彩的遮挡下就有可能失效，不能实现真正意义的最大功率跟踪。目前，国际上也有人提出了多峰值的MPPT算法，主要包含如下三种：

    结合常规算法的复合MPPT算法
    Fibonacci法
    短路电流脉冲法
     

       PSO初始化为一群随机粒子(随机解)。然后通过迭代找到最优解。在每一次的迭代中，粒子通过跟踪两个“极值”(pbest，gbest)来更新自己。在找到这两个最优值后，粒子通过下面的公式来更新自己的速度和位置。

    在这里插入图片描述

     第①部分称为【记忆项】，表示上次速度大小和方向的影响；

    第②部分称为【自身认知项】，是从当前点指向粒子自身最好点的一个矢量，表示粒子的动作来源于自己经验的部分；

    第③部分称为【群体认知项】，是一个从当前点指向种群最好点的矢量，反映了粒子间的协同合作和知识共享。粒子就是通过自己的经验和同伴中最好的经验来决定下一步的运动。

    以上面两个公式为基础，再来看一个公式：

    在这里插入图片描述

    2.仿真效果预览

    matlab2022a仿真如下：

    3.MATLAB核心程序

     

    1. function [sys,x0,str,ts] = BFOA_PSO_pwm(t,x,u,flag)
    2. .............................
    3.  
    4.  
    5. persistent Pbest;
    6. %persistent Pbestval;
    7.  
    8. persistent best_index;
    9.  
    10. persistent c1;
    11. persistent c2;
    12. persistent r1;
    13. persistent r2;
    14.  
    15. switch flag,
    16.     case 0,
    17.         sizes = simsizes;
    18.         sizes.NumContStates = 0;
    19.         sizes.NumDiscStates = 0;
    20.         sizes.NumOutputs     = 1;
    21.         sizes.NumInputs      = 2;
    22.         sizes.DirFeedthrough = 0;
    23.         sizes.NumSampleTimes = 1;
    24.         sys = simsizes(sizes);
    25.         x0=[];
    26.         str=[];
    27.         ts=[0.004 ,0.001];
    28.  
    29.         % initialize the static variables
    30.         first = 1;
    31.         stop = 0;
    32.         i = 0;
    33.         mg = 0;
    34.         count = 3;
    35.         Nc_count = 0;
    36.         Nre_count = 0;
    37.         D_out = zeros(1, NP);                       %初始化D_out、U、fitval_current、fitval_new为全0
    38.         D_out_current = zeros(1,NP);
    39.         best_index = 1;
    40.         %U = zeros(1, NP);
    41.         %V = zeros(1, NP);
    42.         fitval_current = zeros(1, NP);
    43.         fitval_new = zeros(1, NP);
    44.         
    45.         fit_order = zeros(1, NP);
    46.         
    47.  
    48.         Pbest = zeros(1, NP);
    49.         Gbest = 0;
    50.         %Gmaxval = 0;
    51.         %X(:) = unifrnd(XL, XU, 1,NP)
    52.         D_out_current(:) = linspace(XL+0.005,XU-0.005,NP);                 %在[-1,1]间均匀取值
    53.         D_out(:) = D_out_current(:);
    54.         MoveStep = unifrnd(-MaxStep, MaxStep, 1,NP);   %初始化单个细菌的移动速度
    55.         MoveStep_PSO = zeros(1, NP);
    56.         stable_flag = 1;                               %初始化，电路未稳定
    57.         stable_count = 0;
    58.         
    59.         c1 = 0.030;
    60.         c2 = 0.030;
    61.         r1 = rand();
    62.         r2 = rand();
    63.  
    64.     case 3,
    65.         %if count == 3
    66.          %count = count + 1;
    67.             count = 1;
    68.         %迭代完成，输出最优值
    69.         if stop == 1;
    70.             [best, best_index]= max(fitval_current);
    71.             sys = D_out(best_index);
    72.             return;
    73.         end          
    74.         %判断是否迭代完成
    75.         if mg > maxgen
    76.             stop = 1;
    77.         end
    78.         
    79.         if mg == 0  % 第一代，只进行迭代，然后计算各自功率,作为初始比较功率fitval_current   
    80.             
    81.             if i == 0
    82.                 i= i + 1;
    83.                 sys = D_out_current(1);  
    84.             elseif i > 1 && i < NP
    85.                 fitval_current(i-1) = adaptfunc(D_out_current(i-1),Uin);   %记录第一代的上一个个体的功率
    86.                 Pbest(i-1) = D_out_current(i-1);
    87.                 i= i + 1;
    88.                 sys = D_out_current(i);            
    89.             elseif i == NP
    90.                 fitval_current(i-1) = adaptfunc(D_out_current(i-1),Uin);
    91.                 Pbest(i-1) = D_out_current(i-1);
    92.                 i= i + 1;
    93.                 sys = D_out_current(NP);
    94.             elseif i == NP + 1
    95.                 fitval_current(i-1) = adaptfunc(D_out_current(i-1),Uin);
    96.                 Pbest(i-1) = D_out_current(i-1);
    97.                 [best, best_index]= max(fitval_current);
    98.                 Gbest = D_out_current(best_index);
    99.                  i = 1;
    100.                  Nc_count = Nc_count + 1;
    101.                  Nre_count = Nre_count + 1;
    102.                  mg = mg + 1;
    103.                  for j = 1 : NP                       %第一代种群中个体进行游动
    104.                     MoveStep_PSO = c1*r1*(Pbest(j)-D_out_current(j)) + c2*r2*(Gbest-D_out_current(j));
    105.                     D_out(j) = D_out_current(j)+ MoveStep(j) + MoveStep_PSO;
    106.                     if D_out(j) <= XL
    107.                             D_out(j) = XL+0.001;
    108.                         elseif D_out(j) >= XU
    109.                             D_out(j) = XU-0.001;
    110.                     end
    111.                  end
    112.                  sys = D_out(1);            
    113.             else
    114.                 i= i + 1;
    115.                 sys = D_out_current(i);            
    116.             end
    117.             
    118.             return;
    119.             
    120.         else   %第一代之后      
    121.             
    122.             if i > 1 && i < NP
    123.                 fitval_new(i-1) =adaptfunc(D_out(i-1), Uin);      %计算前一次占空比下的功率
    124.                 Pbest(i-1) = D_out(i-1);
    125.                 i = i + 1;
    126.                 sys = D_out(i);
    127.             elseif i == NP
    128.                 fitval_new(i-1) =adaptfunc(D_out(i-1), Uin);
    129.                 Pbest(i-1) = D_out(i-1);
    130.                 i = i + 1;
    131.                 sys = D_out(NP);
    132.             elseif i == 1
    133.                 i = i + 1;            
    134.                 sys = D_out(i);       
    135.             elseif i == NP + 1 %种群迭代一次结束，开始进行新的一代种群繁殖                            
    136.                 %趋向性操作,判定运动方向
    137.                 fitval_new(i-1) =adaptfunc(D_out(i-1), Uin);
    138.                 Pbest(i-1) = D_out(i-1);
    139.                 r1 = rand();
    140.                 r2 = rand();
    141.                 for j = 1 : NP
    142.                 if (fitval_new(j) >= fitval_current(j))
    143.                     D_out_current(j) = D_out(j);
    144.                     fitval_current(j) =  fitval_new(j);
    145.                 else
    146.                     MoveStep(j) = - MoveStep(j);  %若该方向功率未改进，说明不适应生存，改变方向运动
    147.                 end       
    148.                 end
    149.                 
    150.                 c1 = c1/3;
    151.                 c2 = c2/3;
    152.  
    153.                  D_out(:)                              %打印最新的占空比
    154.                 [best, best_index]= max(fitval_current);  %计算最大功率点，best为最大功率，best_index为最大功率点在种群中的位置            
    155.                 mg
    156.                 D_out(best_index)                   %显示最大功率点的占空比
    157.                 %Gbest = U(best_index)               %显示最大功率点的电压
    158.                 %Gmaxval = best                      %显示最大功率点的功率
    159.                 best
    160.                 fit_order = order(fitval_current,NP)    %显示当前种群功率从大到小顺序
    161.  
    162.                 Gbest = D_out(best_index);
    163.                 % mutation
    164.                 if Nc_count == Nc %进行满Nc次趋向性操作            
    165.                     Nc_count = 0;
    166.                     Nre_count = Nre_count + 1;
    167.                     %复制操作
    168.                     MaxStep = MaxStep/2;
    169.                     MoveStep = unifrnd(-MaxStep,MaxStep,1,NP);                 %每Nc次繁殖生成一次新的随机步长
    170.                     [fitval_current,D_out_current]=Reproduction(fitval_current,D_out_current,NP);      %复制
    171.                     D_out=D_out_current;
    172.                     if Nre_count == Nre %进行满Nre次复制操作后迁移操作
    173.                         Nre_count = 0;
    174.                         for j=1:NP
    175.                         if(rand(0,1)<Ped)
    176.                             fitval_current = 0;
    177.                             %MoveStep(i) = unifrnd(-MaxStep,MaxStep);
    178.                             D_out(j) = unifrnd(XL,XU);
    179.                         end
    180.                         end
    181.                     end
    182.                 else
    183.                     Nc_count = Nc_count + 1;
    184.                     for j=1:NP
    185.                         MoveStep_PSO = c1*r1*(Pbest(j)-D_out(j)) + c2*r2*(Gbest-D_out(j));
    186.                         D_out(j) = D_out(j) + MoveStep(j) + MoveStep_PSO;
    187.                         if D_out(j) <= XL
    188.                             D_out(j) = XL+0.001;
    189.                         elseif D_out(j) >= XU
    190.                             D_out(j) = XU-0.001;
    191.                         end
    192.                     end                
    193.                 end          
    194.             %为下一次迭代做准备
    195.             sys = D_out(1);
    196.             i = 1;
    197.             mg = mg +1;
    198.             end
    199.         end
    200.         return;
    201.       
    202.     case 2,
    203.         Uin = u;
    204.     case {1,4,9},
    205.         sys = [];
    206. end
    207.  
    208. A85

    4.完整MATLAB

    V

  • 相关阅读:
    2023-06-17：说一说redis中渐进式rehash？
    4、SpringMVC获取请求参数
    sku详情接口
    基于达梦数据库开发-C#篇
    如何查看 Docker 安装的事件历史记录
    【1684. 统计一致字符串的数目】
    选择结构——字符加密
    Haproxy搭建Web集群
    ArrayList、LinkedList、Collections.singletonList、Arrays.asList与ImmutableList.of
    软件工程导论概述----软件的生命周期
  • 原文地址：https://blog.csdn.net/hlayumi1234567/article/details/128165149