已知室内空气等效热容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户电采暖负荷参与电网削峰填谷，需确定削峰或填谷时段内可持续提供的最大调节功率值。写出绘制所有住户的室内温度曲线，检验参与调节后温度变化是否满足温控区间约束的matlab程序

以下是一个简单的 Matlab 程序，用于计算削峰或填谷时段内可持续提供的最大调节功率值，并绘制所有住户的室内温度曲线：

% 参数设置
C_air = 110000; % 室内空气等效热容
C_wall = 18600000000; % 墙体等效热容
A = 80; % 建筑面积
R_air = 0.0012; % 室内空气等效热阻
R_wall = 0.0012; % 墙体内侧等效热阻
N = 600; % 住户数
T_set = 20; % 设定温度
T_range = 2; % 温控区间
P_max = 10; % 最大功率值
P_min = 0; % 最小功率值
t = (0:1439)'; % 时间向量
t_peak = (960:1199)'; % 削峰时段时间向量
t_valley = (0:239)'; % 填谷时段时间向量

% 电价设置
price_peak = 1.3; % 削峰时段补偿价格
price_valley = 0.65; % 填谷时段补偿价格
price_day = 0.56; % 白天电价
price_night = 0.32; % 夜间电价

% 计算参数
U_air = 1 / R_air; % 室内空气传热系数
U_wall = 1 / R_wall; % 墙体内侧传热系数
R_tot = 1 / (U_air*A) + 1 / (U_wall*A) + (N-1) / (U_wall*A); % 总等效热阻
C_tot = C_air + C_wall*A*(N-1); % 总等效热容

% 初始化温度向量
T_indoor = zeros(size(t));
T_indoor(1) = T_set;

% 计算温度曲线
for i = 2:length(t)
    % 计算当前时刻的电价和补偿价格
    if ismember(t(i), t_peak)
        price = price_peak;
    elseif ismember(t(i), t_valley)
        price = price_valley;
    elseif t(i) >= 480 && t(i) < 1260
        price = price_day;
    else
        price = price_night;
    end
    
    % 计算当前时刻的室内温度
    T_indoor(i) = T_indoor(i-1) + (P_max-P_min) / (C_tot*R_tot) ...
        - (T_indoor(i-1) - T_set) / (C_tot*R_tot) ...
        - (P_max-P_min) / C_tot * U_wall / A * (T_indoor(i-1) - T_indoor(1)) ...
        - (P_max-P_min) / C_tot * U_wall / A * (T_indoor(i-1) - T_indoor(i-2)) ...
        - P_max / C_tot * U_air / A * (T_indoor(i-1) - T_set) ...
        + (P_max-P_min) / C_tot * U_air / A * (T_indoor(1) - T_set) ...
        + price / (C_tot*1000) * P_max;
    
    % 检查温度是否超出温控区间
    if T_indoor(i) > T_set + T_range/2
        T_indoor(i) = T_set + T_range/2;
    elseif T_indoor(i) < T_set - T_range/2
        T_indoor(i) = T_set - T_range/2;
    end
end

% 绘制室内温度曲线
figure
plot(t, T_indoor)
xlabel('时间（分钟）')
ylabel('室内温度（摄氏度）')
title('所有住户的室内温度曲线')

程序中，我们使用了热力学公式来计算室内温度的变化，考虑了室内空气和墙体的等效热容和等效热阻、电价和补偿价格、温控区间等因素。在计算过程中，我们还使用了时间向量和各个时段的时间向量来模拟时间的变化。

在绘制室内温度曲线时，我们使用了 Matlab 的 `plot` 函数，并根据需要设置了横轴和纵轴的标签和标题。