Question

用matlab求微分方程 dy/dt+y tany=cosy;y0=1的数值解,并绘制曲线.

Answer 1

如何用matlab求微分方程 dy/dt+y tany=cosy;y(0)=1的数值解,并绘制曲线.

这个问题我们可以这样来考虑：

第一步，考虑用ode函数求解的自定义函数，即

function dy=odefun(t,y)

dy=cos(y)-y*tan(y);

第二步，确定y的初值，即y0=1；

第三步，确定t的变化范围，如tspan=【0，5】；

第四步，使用ode45函数得到其数值解，即

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

第五步，使用plot函数，绘制t—y的曲线图。即

plot(t,y,'k-')

第六步，标注坐标轴名称，即

xlabel('t'),ylabel('y(t)')

第七步，执行上述代码，得到如下结果