已知室内空气等效热容110000、墙体等效热容18600000000，建筑面积为80平方米，8：00-21：00电价为0.56，21：00-8：00电价为0.32，削峰时段16：00-20：00的补偿价格为1.3元，填谷时段0：00-4：00的补偿价格为0.65元室内空气和墙体内侧的等效热阻分别为0.0012、聚合商组织住宅区600户电采暖负荷参与电网削峰填谷，需确定削峰或填谷时段内可持续提供的最大调节功率值。问题4所解出的各时点可上调、下调功率结果是基于单纯满足温控区间约束条件的电采暖设备开关状态决定的，电采暖负荷参与功率调节将改变其原有的开关状态，进而影响后续可调节功率的时变特性，计算600户电采暖负荷在削峰时段可提供的持续最大向上调节功率值，假设上述计算所得持续最大向上、向下调节功率全部被调度中心调用，统计各时点由于参与电网调节导致开、关状态发生变化的电采暖设备数量，绘制所有住户的室内温度曲线，检验参与调节后温度变化是否满足温控区间约束的matlab程序

以下是一种可能的Matlab程序，仅供参考：

% 室内空气等效热容
C_air = 110000;
% 墙体等效热容
C_wall = 18600000000;
% 建筑面积
A = 80;
% 室内空气和墙体内侧的等效热阻
R_air = 0.0012;
R_wall = 0.0012;
% 电价
E_price = [0.56*ones(1,13),0.32*ones(1,11)];
% 削峰时段补偿价格
C_price_peak = 1.3;
% 填谷时段补偿价格
C_price_valley = 0.65;
% 削峰时段
peak_start = 16;
peak_end = 20;
% 填谷时段
valley_start = 0;
valley_end = 4;
% 电采暖负荷参与电网削峰填谷的户数
N = 600;
% 计算削峰时段可持续提供的最大向上调节功率值
P_peak_max = 0;
for t = peak_start:peak_end
    % 计算当前时刻的室内温度
    T_in = ...  % 根据电采暖负荷模型计算室内温度，略
    % 计算当前时刻的电采暖负荷
    P_load = ...  % 根据电采暖负荷模型计算负荷，略
    % 计算当前时刻的电价
    E = E_price(t);
    % 计算当前时刻的补偿价格
    C = E;
    if t>=peak_start && t<=peak_end
        C = C_price_peak;
    elseif t>=valley_start || t<=valley_end
        C = C_price_valley;
    end
    % 计算当前时刻可提供的最大向上调节功率值
    P_peak_max = max(P_peak_max,(C-E)/R_air/A-(T_in-T_set)/R_wall/A-P_load);
end
% 计算削峰时段可提供的持续最大向上调节功率值
P_peak_sustain = P_peak_max * N;
% 统计各时点由于参与电网调节导致开、关状态发生变化的电采暖设备数量
num_switch = zeros(24,1);
for t = 1:24
    % 计算当前时刻的室内温度
    T_in = ...  % 根据电采暖负荷模型计算室内温度，略
    % 计算当前时刻的电采暖负荷
    P_load = ...  % 根据电采暖负荷模型计算负荷，略
    % 计算当前时刻的电价
    E = E_price(t);
    % 计算当前时刻的补偿价格
    C = E;
    if t>=peak_start && t<=peak_end
        C = C_price_peak;
    elseif t>=valley_start || t<=valley_end
        C = C_price_valley;
    end
    % 计算当前时刻可提供的最大向上调节功率值
    P_peak_max_t = (C-E)/R_air/A-(T_in-T_set)/R_wall/A-P_load;
    % 调度中心调用持续最大向上调节功率值
    P_peak = min(P_peak_max_t,P_peak_sustain);