136
- 0.6
- 0.4
-0.2
0
0.2
0.4
0.6
- 3
- 2
- 1
0
1
2
3
x
1
x
2
Hì
nh 6-10. Hình ảnh quĩ đạo pha với điều kiện đầu [0, 0]
Để so sánh, ta sẽ giải hệ trên với điều kiện đầu
(0)
(0)
1
y
y
=
=
ɺ
. Lặp lại các dòng
lệnh trên, ta nhận được đồ thị như trên các hình 6-11,
>> [t,x] = ode45('eqx2',[0 2*pi],[1,1]);
>> plot(t,x(:,1), 'k-',t,x(:,2),'k--'), xlabel('t'),
>> axis([0 2*pi -4.2 4.2]), legend('x_1','x_2')
0
1
2
3
4
5
6
- 4
- 2
0
2
4
t
x
1
x
2
Hì
nh 6-11. Ứng xử của hệ với điều kiện đầu [1, 1],
Như thấy trên đồ thị, nghiệm của hệ đã thay đổi rất nhiều khi điều kiện đầu thay
đổi từ [0, 0] sang [1, 1]. Các nghiệm tăng sớm hơn ngay trong những giây đầu tiên.
Để so sánh x1 với cos(4.3t) ta vẽ hai đại lượng này trên cùng một hệ trục tọa độ,
như trên hình 6-12.
