已知室内空气等效热容110000、墙体等效热容18600000000，建筑面积为80平方米，假设供暖期为180天，室外平均温度及持续天数如表2 所示，试估算各室外温度下该住宅区电采暖负荷参与削峰填谷的上调、下调功率，根据辅助服务补偿价格，用matlab写出绘制所有住户的室内温度曲线的程序

对于估算各室外温度下该住宅区电采暖负荷参与削峰填谷的上调、下调功率，可以使用热传导方程：

Q = U*A*ΔT

其中，Q为传热量，U为传热系数，A为传热面积，ΔT为温度差。

将上式代入室内空气和墙体的等效热容公式中，得到室内温度随时间变化的微分方程：

C1*dT/dt = Uw*A*(Tw-T) + Qh

C2*dT/dt = Qh - Qc

其中，C1和C2分别为室内空气和墙体的等效热容，Tw为室外温度，T为室内温度，Qh为供暖功率，Qc为散热功率。

根据削峰填谷的要求，需要计算在室外温度变化时，需要上调或下调的功率，即：

P_up = max(Qh) - Qh

P_down = Qh - min(Qh)

为了绘制所有住户的室内温度曲线，可以使用ode45函数求解微分方程，并将结果绘制在图表上。

以下是用Matlab实现的程序：

% 室内空气等效热容
C1 = 110000;
% 墙体等效热容
C2 = 18600000000;
% 建筑面积
A = 80;
% 供暖期天数
days = 180;
% 室外温度及持续天数
Tw = [-15, -10, -5, 0, 5, 10, 15, 20, 25, 30];
Tdays = [30, 30, 30, 30, 30, 30, 30, 30, 30, 30];
% 室内初始温度
T0 = 15;
% 供暖功率
Qh = zeros(1, days);
% 散热功率
Qc = zeros(1, days);
% 室内温度
T = zeros(1, days);
T(1) = T0;
for i = 1:length(Tw)
    % 计算传热系数
    if Tw(i) < 0
        Uw = 3.6;
    elseif Tw(i) < 10
        Uw = 3.2;
    else
        Uw = 2.8;
    end
    % 计算供暖功率和散热功率
    Qh(Tdays(i)+1:Tdays(i)+30) = Uw*A*(Tw(i)-T0);
    Qc(Tdays(i)+1:Tdays(i)+30) = Uw*A*(T0-Tw(i));
    % 求解微分方程
    [T1, Y] = ode45(@(t, y) [(Qh(Tdays(i)+round(t))-Uw*A*(Tw(i)-y(1)))/C1; (Uw*A*(Tw(i)-y(1))-Uw*A*(y(1)-y(2)))/C2], [0 30*24*3600], [T0 T0]);
    T(Tdays(i)+1:Tdays(i)+30) = Y(:, 1)';
    % 更新初始温度
    T0 = Y(end, 1);
end
% 计算上调和下调功率
P_up = max(Qh) - Qh;
P_down = Qh - min(Qh);
% 绘制曲线图
figure;
plot(T);
xlabel('时间（天）');
ylabel('温度（℃）');
title('所有住户的室内温度曲线');

注意，该程序只是一个简单的估算模型，实际情况会受到很多因素的影响，如建筑结构、采暖设备、室内人数等，需要根据实际情况进行调整。同时，在使用该程序时需要注意保护个人隐私，不要泄露住户的信息。