光伏发电系统仿真模型，扰动观察法，带储能控制策略，附负载突变进行验证 支持各版本MATLAB

光伏系统仿真这玩意儿说难不难，说简单也得掉几根头发。今天咱们直接撸代码实战，用MATLAB整一个带储能的光伏发电模型，重点聊聊怎么用扰动观察法追最大功率点，顺便给负载突然搞事情时系统能不能稳住。

先搞个光伏阵列的数学模型，核心就是那串I-V特性方程。这里直接上函数封装：

function [I_pv, P_pv] = PVmodule(V_pv, G, T)
    k = 1.3806e-23;
    q = 1.602e-19;
    Isc = 8.2 * (G/1000);
    Voc = 32 + (T-25)*0.08;
    Vt = (k*(T+273.15)/q)*54;
    I0 = Isc / (exp(Voc/(Vt)) - 1);
    I_pv = Isc - I0*(exp(V_pv/(Vt)) - 1);
    P_pv = V_pv .* I_pv;
end

参数里的G是光照强度，T是温度。注意指数运算那行容易数值爆炸，遇到报错就把V_pv限制在Voc的90%以内。

接着上扰动观察法（Perturb and Observe）的控制器：

classdef P_O_controller < handle
    properties
        V_step = 0.2;  % 扰动步长
        prev_P = 0;    % 上次功率
        V_ref = 30;    % 初始参考电压
    end
    
    methods
        function [V_ref] = update(obj, P_curr)
            if (P_curr - obj.prev_P) > 0
                direction = sign(obj.V_ref - 30); 
            else
                direction = -sign(obj.V_ref - 30); 
            end
            obj.V_ref = obj.V_ref + direction * obj.V_step;
            obj.prev_P = P_curr;
            V_ref = obj.V_ref;
        end
    end
end

重点在方向判断逻辑——当功率增加时保持原扰动方向，否则反向。这里有个坑：光照突变时可能误判，后面储能系统就是用来填这个坑的。

光伏发电系统仿真模型，扰动观察法，带储能控制策略，附负载突变进行验证 支持各版本MATLAB

储能部分用了个简化电池模型：

function [I_bat, soc] = BatteryControl(V_bus, soc, P_load, P_pv)
    V_nom = 48;  % 标称电压
    bat_capacity = 200;  % Ah
    
    P_diff = P_pv - P_load;
    if V_bus < V_nom*0.95 || P_diff < 0 
        I_bat = min(100, (P_diff)/V_bus);  % 放电
    elseif V_bus > V_nom*1.05 
        I_bat = -min(50, abs(P_diff)/V_bus); % 充电
    else
        I_bat = 0;
    end
    
    soc = soc - I_bat * 0.1 / bat_capacity;  % 0.1秒步长
end

当母线电压低于阈值或光伏供电不足时，电池放电；电压过高时充电。注意soc计算要考虑仿真步长，这里假设步长0.1秒。

最后来个负载突变测试场景：

% 仿真参数
t = 0:0.1:60; 
load_profile = 2000*ones(size(t)); 
load_profile(t>20 & t<=40) = 3500;  % 20秒时负载突增

% 主循环
for i = 1:length(t)
    [I_pv(i), P_pv(i)] = PVmodule(po.V_ref, 800, 25);
    V_bus(i) = po.V_ref;  % 假设理想DC/DC
    
    [I_bat(i), soc(i)] = BatteryControl(V_bus(i), soc_prev, load_profile(i), P_pv(i));
    
    % 功率平衡校验
    P_bus = P_pv(i) + V_bus(i)*I_bat(i);
    assert(abs(P_bus - load_profile(i)) < 1e-3, '功率不平衡!');
    
    % 更新MPPT
    po.update(P_pv(i));
    soc_prev = soc(i);
end

运行后会看到在20秒负载突增时，电池电流瞬间从0飙升到正方向（放电），同时MPPT控制器逐步调整光伏电压追踪新功率点。把V_step参数调大到0.5会发现系统震荡明显，这就是扰动观察法的经典trade-off。

整套模型在MATLAB R2016a到2023b都跑过，注意高版本可能提示用classdef语法需要确认，老版本如果报错就把控制器改写成函数形式。仿真结果拿plot画个三图联动（电压、功率、SOC），比论文里的波形还直观。