, 0.001 , , , , . , .
5 10.
>> hold on
>> for = 5:10
[t, xa] = ode45(g, [0:0.01:20], [0 ]);
plot(xa(:, 1), xa(:, 2))
end
>> hold off
. 4.
5, 6 7 (. 4), π, 8, 9 10 , . , 7 8.
>> hold on
>> for = 7.0:0.2:8.0
[t, ] = ode45(g, [0:0.01:20], [0 ]);
plot(xa(:, 1), xa(:, 2))
end
>> hold off
. 5.
(. 6), - 7.2 7.4. .
>> hold on
>> for = 7.2:0.05:7.4
[t, ] = ode45(g, [0:0.01:20], [0 ]);
plot(xa(:, 1), (:, 2))
end
hold off
. 6.
, 0,01 - 7.25 7.3 (. 7 8).
for a = 7.25:0.01:7.3
[t, xa] = ode45(g, [0:0.01:20], [0 a]);
plot(xa(:, 1), xa(:, 2))
end
. 7.
, .
. 8.
, - 7.29 7.3.
.
w2 = 2 c = 1
-, , (. 2013, . 117):
function dydt = vdp1(t,y)
dydt = zeros(2,1); dydt(1) = y(2);
dydt(2) = 1*(1-y(1).^2).*y(2)-2*y(1);
(0 = 2, v0= 0) - :
[t,y] = ode23(@vdp1,[0 25],[2;0]);
v
plot(t,y(:,1),t,y(:,2)), grid on
. 1.
. 1.
(. 2):
plot(y(:,1),y(:,2)), grid on
. 2.
(0 = -0.5, v0= 5) - :
|
|
[t,y] = ode23(@vdp1,[0 25],[-0.5;5]);
. 3.
. 3.
. 4.
. 4.
(0 = -0.05, v0= -0.05) - :
[t,y] = ode23(@vdp1,[0 50],[-0.05; -0.05]);
. 5.
. 5.
. 6.
. 6.