ich benutzeode45 um eine einfache ODE zu lösen:

function dCdt=u_vent(t,C)
if t>   600 &&  t<= 720
    Q=Q2;               
elseif t>   1320    &&  t<= 1440
    Q=Q2;           
elseif t>   2040    &&  t<= 2160
    Q=Q2;           
elseif t>   2760    &&  t<= 2880
    Q=Q2;               
elseif t>   3480    &&  t<= 3600
    Q=Q2;
else
    Q=Q1;
end

V=100;
C_i=400;
S=100;

dCdt=Q/V*C_i+S/V-Q/V*C(1);
return

Ich benutze dann zu lösen:

[t,C]=ode45(@u_vent, [0 1*3600], 400);

Ich möchte eine periodische Funktion wie die auf dem Bild für erstellenQ, 0<t<3600, ohne diese zu verwendenif Aussagen ... irgendwelche Gedanken?