虽然你的系统是高度非线性的,
数值
可以评估
sinc ()
函数。然而,你选择的初始值
导致轨迹发散
,奇点(除零)发生在这一项
。
导致轨迹发散
,奇点(除零)发生在这一项
。
(x, y) = meshgrid (-4.9:0.35:4.9);
z1 = 1。/ (5 ^ 2 - x ^ 2);
z2 = 1。/ (5 ^ 2 - y ^ 2);
w = - * (z1 (x。/ z2 + 1) + 2 * z1。* (x ^ 2)。* sinc (y)。/ z2 + y。y ^ 2 +);
u = y + x。* sin (y);
v = y。^ 2 + x + w;
l = streamslice (x, y, u, v);
轴紧
例如,选择另一组初始值
导致轨迹收敛到原点。
导致轨迹收敛到原点。
tspan = linspace (0, 101);
x0 = (-4.5 - 4.5);
(t, x) =数值(@midterm2 tspan, x0);
情节(t x,“线宽”1.5),网格在
包含(“t”),ylabel (“\ bf {x}”),
传奇(”间的{1}”,”间的{2}”)
函数dstatesdt = midterm2 (t,状态)
x1 =状态(1);
x2 =状态(2);
z1 = 1 / (5 ^ 2 - x1 ^ 2);
z2 = 1 / (5 ^ 2 - x2 ^ 2);
% u = (x1 * (z1 / z2 + 1) + z1 * x1 ^ 2 * sinc (x2) / z2 + x2 ^ 2 + x2); x
u = - * (z1 / z2 + 1) (x1 + 2 * z1 * x1 ^ 2 * sinc (x2) / z2 + x2 ^ 2 + x2);
% u = (x1 * (z1 / z2 + 1) + z1 * x1 ^ 2 * sinc (x2) / z2 + x2 ^ 2 + x2 + x2 / z1);
x1dot = x2 + x1 * sin (x2);
x2dot x2 = ^ 2 + x1 + u;
dstatesdt = [x1dot;x2dot];
结束
