变步长的INC_MPPT算法仿真（simulink ；PV模型）

本资源为变步长的INC算法仿真的MPPT算法，采用MATLAB的simulink模块仿真，附有变步长的INC算法源代码(s函数),默认最大功率50W，其他参数均写在仿真中，供大家学习参考。 注意：需MATLABR2016b及以上版本才可运行！

function [sysx0strts] = MPPT_AD_INC（txuflag）
switch flag
  case 0
    [sysx0strts] = mdlInitializeSizes;
  case 2                                                
    sys = mdlUpdate（txu）; 
  case 3                                                
    sys = mdlOutputs（txu）;
  case {149}                                                
    sys =[];
  otherwise
    error（[‘unhandled flag = ‘num2str（flag）]）;
end

function [sysx0strts] = mdlInitializeSizes
sizes = simsizes;
sizes.NumContStates  = 0;
sizes.NumDiscStates  = 4;
sizes.NumOutputs     = 1;
sizes.NumInputs      = 2;
sizes.DirFeedthrough = 0;
sizes.NumSampleTimes = 1;
sys = simsizes（sizes）;
x0  =[0 0 0 0];
str =[];
ts  =[0.002 0];    %步长为1e-6

function sys = mdlUpdate（txu）  
% 增量电导法MPPT
dU=u（1）-x（1）;
dI=u（2）-x（2）;

slope=0;
n=0.0002;   %占空比步长

%变步长
Pnew=u（1）*u（2）;
P=x（1）*x（2）;
dP=Pnew-P;
if dU~=0
    m=dP/dU;
    k1=n*abs（m）;
else k1=0.05;
end
if k1>0.05;
    k1=0.05;
end
if dU==0
    if dI==0
    slope = 0;
    elseif dI > 0
    slope = 1;            
    elseif  dI < 0
    slope = -1;  
    end
else
    if  （dI/dU） ==（-u（2）/u（1））
        slope = 0;    
    elseif （dI/dU） >（-u（2）/u（1））
        slope = -1;   
    elseif  （dI/dU） <（-u（2）/u（1））
        slope = 1; 
    end
end
%找到最大功率点的判断
if abs（u（2）+u（1）*dI/dU）<2e-5
    x（4）=1;
else
    x（4）=0;
end

if x（4）==0
    x（3）=x（3）+k1*slope;     %更新占空比值  
else
    x（3）=x（3）;
end

x（1）=u（1）;
x（2）=u（2）;

sys=[x（1） x（2） x（3） x（4）];

function sys = mdlOutputs（txu）
if x（3）<=0
    x（3）=0.1;
end
if x（3）>=0.8
    x（3）=0.8;
end

sys=x（3）;     % 输出占空比值

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件        1718  2018-12-23 08:52  作业四\MPPT_AD_INC.m
     文件      109004  2018-12-22 22:32  作业四\my_PVtest.mdl
     目录           0  2019-01-04 10:41  作业四\