请给上面的代码加一段程序使得温度在一定范围内成正弦波动曲线

要让温度在一定范围内成正弦波动曲线，可以在每个时间步长内，通过计算正弦函数的值对房间内净热流进行微小的变化，从而引起温度的正弦波动。可以在原有的代码中加入如下一行代码：

Q = (P - Q1 - Q2 + 100*sin(2*pi*t(i)/86400))*dt; % 加入正弦波扰动

其中，`sin(2*pi*t(i)/86400)`表示以一天为周期的正弦波，`86400`是一天的秒数。乘以100是为了保证扰动的幅度适当。这样，每个时间步长内的净热流就会因为正弦波扰动而微小变化，从而导致温度的正弦波动。

完整代码如下所示：

C_air = 1000; % J/(m^3*K)
C_wall = 100000; % J/(m^2*K)
R1 = 0.1; % m^2*K/W
R2 = 0.2; % m^2*K/W
T_in = 20; % °C
T_wall = 20; % °C
T_out = -10; % °C
P_heat = 2000; % W
P_rat = 3000; % W
S = @(t) (mod(t, 3600) < 1800); % 一小时开一小时关

t0 = 0;
tf = 24*3600;
dt = 60;

T = T_in*ones(1, (tf-t0)/dt+1);

for i = 2:length(T)
    t = (i-1)*dt; % 当前时间
    P = S(t)*P_heat;
    C = C_air + C_wall;
    R = R1 + R2;

    Q1 = (T_in - T_wall)/R1;
    Q2 = (T_wall - T_out)/R2;

    Q = (P - Q1 - Q2 + 100*sin(2*pi*t/86400))*dt; % 加入正弦波扰动

    dT = Q/C;

    T(i) = T(i-1) + dT;
    T_wall = T_wall + (Q1 - Q2)*dt/C_wall;
end

t = t0:dt:tf;
plot(t/3600, T);
xlabel('Time (h)');
ylabel('Temperature (°C)');

注意，由于加入了正弦波扰动，每次运行程序得到的温度曲线都会有所不同。