level 1
子狐艺灵
楼主
试了很多方法都不行，这是一个非自治的微分方程，一直显示有奇点或者是刚性方程，有没有办法弄过去，吧里的方法都试过了也不行
s = NDSolve[{x'[t] == 0.3 y[t],
y'[t] == -0.1 x[t] - 0.01 y[t] + 0.2 y[t]^2 + 0.5 Cos[2 Pi t] + 4,
x[0] == 1, y[0] == 1}, {x[t], y[t]}, {t, 0, 50 Pi}]
ParametricPlot[Evaluate[{x[t], y[t]} /. s], {t, 0, 50 Pi},
PlotStyle -> Red]
2019年01月17日 01点01分
1
s = NDSolve[{x'[t] == 0.3 y[t],
y'[t] == -0.1 x[t] - 0.01 y[t] + 0.2 y[t]^2 + 0.5 Cos[2 Pi t] + 4,
x[0] == 1, y[0] == 1}, {x[t], y[t]}, {t, 0, 50 Pi}]
ParametricPlot[Evaluate[{x[t], y[t]} /. s], {t, 0, 50 Pi},
PlotStyle -> Red]
