Cin、Cwall分别为室内空气等效热容、墙体等效热容；R1、R2分别为室内空气和墙体内侧的等效热阻、墙体外侧和室外空气的等效热阻；自变量室内温度、墙体温度、室外温度；室外温度为0，-5，-10，-15，-20，-25，Pheat(t)为电采暖设备制热功率，Pheat(t)=S(t)PN，其中PN为电采暖设备的额定功率，S(t)为电采暖设备的开关状态，关闭时取0，开启时取1。用matlab描述房间温变过程的集总参数常微分方程，室内初始温度为20℃，峰时（8:00-21:00）电价为0.56，谷时（21:00-次日8:00）电价为0.32，在表1给定的室外温度下，计算并绘制一日24h的室内温度变化和相应的电采暖设备开关状态曲线,计算平均升温时长，平均降温时长，周期，平均占空比，日用电量，日平均用电功率，日用电成本

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

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;
peak_price = 0.56;
offpeak_price = 0.32;

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

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

% 定义电采暖设备开关状态函数
S = @(t) (mod(floor(t/3600), 24) >= 8 && mod(floor(t/3600), 24) < 21) * 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);
    % 计算电费
    daily_energy = Pheat_sum / 3600 / 1000;
    peak_hours = sum(S(t)) * (21-8);
    offpeak_hours = 24 - peak_hours;
    daily_cost = peak_price * peak_hours + offpeak_price * offpeak_hours;
    % 计算平均升温时长、平均降温时长、周期、平均占空比
    [t_rise_avg, t_fall_avg, period, duty_cycle] = calc_features(Tin, S(t));
    % 绘制图像
    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']);
    disp(['平均升温时长：', num2str(t_rise_avg), 'h']);
    disp(['平均降温时长：', num2str(t_fall_avg), 'h']);
    disp(['周期：', num2str(period), 'h']);
    disp(['平均占空比：', num2str(duty_cycle), '%']);
    disp(['日用电量：', num2str(daily_energy), 'kWh']);
    disp(['日平均用电功率：', num2str(daily_energy/24), 'kW']);
    disp(['日用电成本：', num2str(daily_cost), '元']);
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

function [t_rise_avg, t_fall_avg, period, duty_cycle] = calc_features(Tin, S)
    % 计算平均升温时长
    t_rise_arr = [];
    for i = 1:length(Tin)-1
        if S(i) == 1 && S(i+1) == 1 && Tin(i+1) > Tin(i)
            t_rise_arr = [t_rise_arr, 1];
        end
    end
    t_rise_avg = mean(t_rise_arr);
    % 计算平均降温时长
    t_fall_arr = [];
    for i = 1:length(Tin)-1
        if S(i) == 0 && S(i+1) == 0 && Tin(i+1) < Tin(i)
            t_fall_arr = [t_fall_arr, 1];
        end
    end
    t_fall_avg = mean(t_fall_arr);
    % 计算周期和平均占空比
    period_arr = [];
    duty_cycle_arr = [];
    for i = 1:length(Tin)-1
        if S(i) == 1 && S(i+1) == 0
            period_arr = [period_arr, 1];
            duty_cycle_arr = [duty_cycle_arr, (i+1-mod(i, 3600))/3600];
        end
    end
    period = mean(period_arr);
    duty_cycle = mean(duty_cycle_arr) / 24 * 100;
end

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

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