矩阵实验室与微分方程矩阵实验室能否求解下面这个微分方程?x''+[(x+2)^2]x'-5=0x'=dx/dt,x''=dx'/dtx(0)=0,x'(0)=0若不能给出解析解,能否用数值计算的办法画出x(t)的曲线?怎么没人知道？

Question

矩阵实验室与微分方程
矩阵实验室能否求解下面这个微分方程?
x''+[(x+2)^2]x'-5=0
x'=dx/dt,x''=dx'/dt
x(0)=0,x'(0)=0
若不能给出解析解,能否用数值计算的办法画出x(t)的曲线?
怎么没人知道？

Answer 1

先做个函数文件,保存在eqx.m中,文件内容如下——
function xdot=eqx(t,x)
xdot=zeros(2,1);xdot(1)=x(2);xdot(2)=5-(x(1)+2)^2*x(2);
再在命令窗口中输入——
[t,x]=ode45('eqx',[0,2],[0,0]);
plot(t,x(:,1),t,x(:,2),':')
注意[0,2]是t的取值范围,可按需要选取.