带有阶跃函数的微分方程，一般只能用数值的方法来求解。可以用ode45函数求解。
1、调用格式，[t,y] = ode45(@myodefun,tspan,y0);
式中：myodefun——微分方程的自定义函数，tspan——时间区间，y0——初始条件
2、微分方程的自定义函数，主要内容为
dy(1)=y(2);
dy(2)=-4*y(2)-4*y(1)+df+3*f;
3、部分计算结果及图形
