Question

matlab中dsolve求微分方程的问题

y’’(t)+ 4y’(t)+4y(t)=f’(t)+3f(t) f(t)= exp(-t)u(t) ;起始条件都为零
syms t y
f=exp(-t).*heaviside(t);
dsolve(['D2y+4*Dy+4*y=','Df+3*f'],'y(0)=3','Dy(0)=0')
plot(t,y)

警告: Unable to display result due to a previously interrupted computation or out of memory. Execute
'reset(symengine)' and repeat the commands to regenerate the result.
> In sym/disp (line 44)
In sym/display>displayVariable (line 85)
In sym/display (line 47)
In savintest (line 3)
错误使用 plot
数据必须为可转换为双精度值的数值、日期时间、持续时间或数组。

出错 savintest (line 4)
plot(t,y)

Answer 1

带有阶跃函数的微分方程，一般只能用数值的方法来求解。可以用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、部分计算结果及图形