Cin、Cwall分别为室内空气等效热容、墙体等效热容；R1、R2分别为室内空气和墙体内侧的等效热阻、墙体外侧和室外空气的等效热阻；自变量室内温度、墙体温度、室外温度；室外温度为0，-5，-10，-15，-20，-25，Pheat(t)为电采暖设备制热功率，Pheat(t)=S(t)PN，其中PN为电采暖设备的额定功率，S(t)为电采暖设备的开关状态，关闭时取0，开启时取1。用matlab描述房间温变过程的集总参数常微分方程，室内初始温度为20℃，在表1给定的室外温度下，计算并绘制一日24h的室内温度变化和相应的电采暖设备开关状态曲线，统计相关特征量填入表1，并分析室外温度对电采暖设备运行特性及耗电量的影响程序

为了描述房间温变过程的集总参数常微分方程，我们可以采用以下的数学模型：

Cin * dTin/dt = (Twall - Tin) / R1 + Pheat(t)
Cwall * dTwall/dt = (Tin - Twall) / R1 + (Tout - Twall) / R2

其中，Tin是室内温度，Twall是墙体温度，Tout是室外温度，Pheat(t)是电采暖设备制热功率，Cin和Cwall分别为室内空气等效热容和墙体等效热容，R1和R2分别为室内空气和墙体内侧的等效热阻以及墙体外侧和室外空气的等效热阻。

为了求解上述方程组，我们可以采用ode45函数进行求解。具体的matlab程序如下：

% 定义常数
Cin = 1.5 * 1000 * 2.8 * 3;
Cwall = 200 * 1000 * 0.15 * 3;
R1 = 0.1;
R2 = 0.04;
PN = 10;
Tin0 = 20;

% 定义时间范围和初始条件
tspan = [0, 24*3600];
y0 = [Tin0, Tin0, Tin0];

% 定义室外温度数组
Tout_arr = [0, -5, -10, -15, -20, -25];

% 定义电采暖设备开关状态函数
S = @(t) (mod(floor(t/3600), 24) >= 6 && mod(floor(t/3600), 24) < 22) * 1;

% 循环计算所有室外温度下的结果
for i = 1:length(Tout_arr)
    Tout = Tout_arr(i);
    [t, y] = ode45(@(t, y) myode(t, y, R1, R2, Cin, Cwall, S, PN, Tout), tspan, y0);
    Tin = y(:, 1);
    Pheat = S(t) * PN;
    % 统计相关特征量
    T_in_avg = mean(Tin);
    T_in_max = max(Tin);
    T_in_min = min(Tin);
    Pheat_sum = trapz(t, Pheat);
    % 绘制图像
    figure();
    plot(t/3600, Tin, 'r');
    hold on;
    plot(t/3600, S(t)*Tin_max, 'b', 'LineWidth', 2);
    xlabel('时间/h');
    ylabel('温度/℃');
    title(['室外温度为', num2str(Tout), '℃时的温度变化']);
    legend('室内温度', '电采暖设备开关状态');
    % 显示结果
    disp(['室外温度为', num2str(Tout), '℃时的统计特征量：']);
    disp(['平均室内温度：', num2str(T_in_avg), '℃']);
    disp(['最高室内温度：', num2str(T_in_max), '℃']);
    disp(['最低室内温度：', num2str(T_in_min), '℃']);
    disp(['电采暖设备功率总和：', num2str(Pheat_sum), 'W']);
end

function dydt = myode(t, y, R1, R2, Cin, Cwall, S, PN, Tout)
    Tin = y(1);
    T1 = y(2);
    T2 = y(3);
    Pheat = S(t) * PN;
    dTin_dt = (T1 - Tin) / R1 + Pheat / Cin;
    dT1_dt = (Tin - T1) / R1 + (T2 - T1) / R2;
    dT2_dt = (Tout - T2) / R2 + (T1 - T2) / R2;
    dydt = [dTin_dt; dT1_dt; dT2_dt];
end

在上面的程序中，我们采用了ode45函数对常微分方程组进行了求解。其中，myode函数定义了常微分方程组的形式。在循环中，我们依次计算了所有室外温度下的结果，并统计了相关特征量，最后绘制了图像。

注意，在上述程序中，我们假设电采暖设备在每天的6:00至22:00之间开启，其余时间关闭。这是一个简化的假设，实际情况中可能会有所不同。此外，我们还假设墙体和室内空气的温度是相同的，这也是一个简化的假设。如果需要更精细的模拟，可以考虑将墙体温度分为内侧和外侧两部分，并对其进行不同的建模。