已知室内空气等效热容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户电采暖负荷参与电网削峰填谷（削峰时段、填谷时段及补偿价格见附件B），需确定削峰或填谷时段内可持续提供的最大调节功率值。问题4所解出的各时点可上调、下调功率结果是基于单纯满足温控区间约束条件的电采暖设备开关状态决定的，电采暖负荷参与功率调节将改变其原有的开关状态，进而影响后续可调节功率的时变特性，计算600户电采暖负荷在削峰时段可提供的持续最大向上调节功率值得到，假设上述计算所得持续最大向上、向下调节功率全部被调度中心调用，统计各时点由于参与电网调节导致开、关状态发生变化的电采暖设备数量，绘制所有住户的室内温度曲线，检验参与调节后温度变化是否满足温控区间约束的matlab程序

由于该问题涉及到电力系统、建筑热力学和Matlab编程等多个领域，以下是一些思路和参考代码，但需要根据具体情况进行修改和完善。

1. 计算削峰时段可提供的持续最大向上调节功率值

可以先按照问题4的方法，计算出每个时刻可调节的向上和向下功率，然后在削峰时段内，统计所有可调节的向上功率的总和，即为可持续提供的最大向上调节功率值。

2. 统计开关状态发生变化的电采暖设备数量

可以先根据问题4的方法，计算出每个时刻所有电采暖设备的开关状态，然后在削峰时段内，统计与前一时刻开关状态不同的设备数量。

3. 绘制所有住户的室内温度曲线

可以利用建筑热力学模型，根据室内外温度、电采暖设备功率等参数，计算出每个时刻室内温度的变化，并绘制出温度曲线图。

4. 检验参与调节后温度变化是否满足温控区间约束

可以根据设定的温控区间，检验每个时刻的室内温度是否在该区间内，并计算出室内温度不满足温控区间的比例，作为检验指标。

以下是一些参考代码，供参考：

% 读入参数
C_air = 110000; % 室内空气等效热容
C_wall = 18600000000; % 墙体等效热容
A = 80; % 建筑面积
T_out = [2 -2 -3 -5 -6 -6 -6 -6 -5 -4 -2 1 2 3 4 4 3 2 1 0 -1 -2 -2]; % 外界温度数据
T_set = [20 20 20 20 20 20 20 20 22 22 22 22 22 22 22 20 20 20 20 20 20 20 20]; % 温控区间设定值
P_up = ... % 每个时刻可上调的功率值，根据问题4计算得到
P_down = ... % 每个时刻可下调的功率值，根据问题4计算得到

% 计算削峰时段可提供的持续最大向上调节功率值
P_up_max = sum(P_up(17:20)); % 削峰时段为第17~20个时刻

% 统计开关状态发生变化的电采暖设备数量
N_change = 0;
for t = 2:length(P_up)
    for i = 1:600
        % 根据当前时刻的功率值和上一时刻的开关状态，计算出当前时刻的开关状态
        % 如果与上一时刻不同，则开关状态发生了变化
        % 这里省略了具体代码实现
        if state_cur ~= state_prev
            N_change = N_change + 1;
        end
    end
end

% 计算室内温度变化
T_in = zeros(1, length(T_out));
T_in(1) = T_out(1); % 假设初始室内温度等于外界温度
for t = 2:length(T_out)
    % 根据建筑热力学模型，计算出当前时刻的室内温度
    % 这里省略了具体代码实现
    T_in(t) = ...;
end

% 绘制室内温度曲线
figure;
plot(T_in);
title('Indoor Temperature');

% 检验温度变化是否满足温控区间约束
N_violate = 0;
for t = 1:length(T_in)
    if T_in(t) < T_set(t)-0.5 || T_in(t) > T_set(t)+0.5 % 这里将温控区间的容忍度设为0.5度
        N_violate = N_violate + 1;
    end
end
P_violate = N_violate / length(T_in); % 计算出温度不满足温控区间的比例