matlab如何求y''-13(1-y*y)y'+y
发布网友
发布时间:2024-05-03 21:20
我来回答
共1个回答
热心网友
时间:2024-06-10 00:52
如何用matlab求微分方程y''-13(1-y*y)y'+y=0?一般用dsolve()函数来求解,但本题无法使用该函数命令其解析值。所以该微分方程应用ode()函数来求解其初值问题的数值值,其求解过程如下:
1、创建自定义微分方程函数,即
odefun=@(x,y)[y(2);13*(1-y(1)*y(1))*y(2)-y(1)];
2、确定y0初值,即 y0=[0.001,0];
3、确定x的变化范围,即xspan=[0,1];
4、使用ode45函数,得到y、y'值。即
[x,y]=ode45(odefun,xspan,y0)
5、对于数值解,也可以用曲线来表示。即
plot(x,y(:,1))%x-y(x)曲线
plot(x,y(:,2))%x-dy / dx曲线