帮我生成一个园区级综合能源系统简单算例matlab代码

在MATLAB中构建一个简单的园区级综合能源系统模型，我们通常会涉及电力生产、消耗以及优化控制的部分。下面是一个基础的示例，假设我们有一个包含光伏电站和储能系统的系统，目标是最大化清洁能源利用并满足园区需求：

% 定义变量
Ppv = 50; % 光伏发电功率 (千瓦)
Es = 100; % 储能系统容量 (千瓦时)
D = 80; % 园区日平均负荷 (千瓦)

% 定义光伏出力函数 (假设为连续线性)
function P_pv = pv_output(T)
    T = T; % 温度，这里可以替换为实际的光照强度等参数
    P_pv = min(Ppv, 100 - 5 * abs(T - 25)); % 随温度变化的光伏发电量
end

% 定义电池充电和放电效率
eta_charging = 0.9;
eta_discharging = 0.9;

% 功能：计算一天内最小化储能成本的操作策略
function strategy = optimal_strategy(D, Ppv, Es)
    dt = 1; % 时间步长 (小时)
    tspan = 24; % 总时间 (小时)
    T = linspace(0, 24, tspan); % 时间向量

    % 初始化状态矩阵 (初始状态为满充)
    soc = zeros(tspan, 1);
    soc(1) = 100 / Es; % 初始储能比例

    % 循环模拟每个小时
    for i = 2:tspan
        % 根据当前时刻的光伏产量和负荷情况
        P_avail = pv_output(T(i));
        P_load = D;

        % 如果光伏不足以覆盖全部负荷，则使用储能
        if P_avail < P_load
            P_from_storage = min(Es * soc(i-1) * eta_discharging, P_load - P_pv);
            P_from_grid = P_load - P_from_storage;
        else % 光伏多余，则优先充满储能
            P_from_grid = max(0, Ppv - P_load);
            P_from_storage = 0;
        end
        
        % 更新储能状态
        soc(i) = (soc(i-1) - P_from_storage / Es) / eta_charging;

        % 记录操作策略
        strategy{i} = [P_from_storage, P_from_grid];
    end
end

% 调用策略计算函数
strategy = optimal_strategy(D, Ppv, Es);

% 打印结果
disp("最优操作策略：");
for i = 1:length(strategy)
    disp(['时间 ' num2str(i*60) ': 从储能 (' num2str(strategy{i}(1)) ') 和电网 (' num2str(strategy{i}(2)) ') 提供能源']);
end

%