Question

matlab如何求解二阶微分方程数值解的图象

我想解这个微分方程组，但是解析解无法解，求那位大佬教教小白如何敲代码用数值解解出他的图像

Answer 1

题主的问题可以使用ode45函数求解，得到【x，y，y'】的数值解。其求解方法：

1、确定初值，y0=[pi,0];

2、确定x的范围，xspan=[0,10];

3、自定义函数，odefun=@(x,y)[y(2);sin(y(1))*cos(y(1))*((2/9*pi)/(2+0.4*(sin(y(1)))^2))^2-5*sin(y(1))];

4、求解，[x,y]=ode45(odefun,xspan,y0,ff);

5、绘图，plot(x,y,'.-')

6、标注图例，legend('y(x)','dy/dx','Location','northwest')

7、标注坐标轴名称，xlabel('x');ylabel({'y(x)';'dy/dx'});

8、输出结果

追问

万分谢谢楼主啊，就是我还有个问题：问题因为我是刚刚接触这个matlab，因为我目前只用他来解微分方程，不知道楼主有没有相关数值解微分方程内容的链接让小白我学习一下。因为这个方程只是我需要解的一个。我还有几个参数改变的方程我也要解。

追答

这个你网上查找应该能搜索到的。具体可以通过私信，这里不能贴有链接的。

追问

还有一个问题楼主，这个ff如何自定义的啊，因为我把这个复制进去它显示：未定义函数或变量 'ff'

追答

ff=odeset;ff.AbsTol=1e-6;ff.RelTol=1e-6;

Answer 2

别用图片,直接用代码,好用来解下.