134
>> axis([0 10 –1.12 1.12])
Ta nhận được đồ thị hai hàm x1 và x2 theo thời gian như trên hình 6-6.
0
2
4
6
8
10
-1
- 0.5
0
0.5
1
t
Hì
nh 6-7. Hàm x tức x1 nét liền, và hàm y tức x2 nét đứt
Để vẽ quĩ đạo pha, sử dụng lệnh sau
>> plot(x(:,1),x(:,2), 'k-'), axis([-0.65 0.65 -0.4 1.1])
và ta nhận được hình ảnh pha như trên hình 6-8.
- 0.6
- 0.4
- 0.2
0
0.2
0.4
0.6
- 0.4
- 0.2
0
0.2
0.4
0.6
0.8
1
Hì
nh 6-8. Quỹ đạo pha (x1,x2)
Ví dụ dưới đây trình bày việc giải phương trình vi phân bậc hai. Hãy tìm nghiệm
của phương trình vi phân sau
16
sin(4.3 )
y
y
t
+
=
ɺɺ
với điều kiện đầu
(0)
0, (0)
0
y
y
=
=
ɺ
.
Bằng cách đặt
1
2
,
x
y
x
y
=
= ɺ
, ta nhận được hệ hai phương trình vi phân cấp một:
