如何用matlab求二阶微分方程的图像解

Question

我需要用matlab解y"=siny/-0.02-0.04t这个二阶微分方程，其中y0=0.2,y'(0)=2.6，matlab无法给出数值解，求助如何得到图像解（比如ode45等函数），代码请写详细些，谢谢。

Answer 1

如何用matlab求二阶微分方程的图像解。可以按下列方法来求解。

1、自定义二阶微分方程函数，ode_fun(t,y)，其格式

function dy=ode_fun(t,y)

dy=[y(2);sin(y(1))/(-0.02)-0.04*t];

2、利用ode45函数求解[t,y]值

y0=[0.2,2.6];

tspan=[0,10];

[t,y]=ode45(@ode_fun,tspan,y0)

3、利用plot函数绘出二阶微分方程的图像解

polt(t,y)

xlabel('t'),ylabel('y(t)','dy/dt');

4、执行结果