137
>> f1=cos(4.3*t);
>> plot(t,x(:,1),'k-', t,f1,'k--'), xlabel('t'),
>> axis([0 2*pi -1.2 1.9]), legend('x_1','cos(4.3t)')
0
1
2
3
4
5
6
-1
0
1
t
x
1
cos(4.3t)
Hì
nh 6-12. So sánh nghiệm x1 với cos(4.3t), điều kiện đầu [1, 1]
0
1
2
3
4
5
6
- 1
- 0.5
0
0.5
1
t
x1
cos(4.3t )
sin(4.3t )
Hì
nh 6-13. So sánh x1 với sin và cos(4.3t), khi điều kiện đầu là [0, 0]
Hai đường này bắt đầu từ hai điểm rất gần nhau, nhưng chúng xa dần khi thời gian
tăng lên. Nhưng so với trường hợp trước, khi điều kiện đầu là [0, 0], ta thấy các
nghiệm tăng dần (hình 6-13), như vậy có thể kết luận là các nghiệm không tăng lên
mãi theo thời gian. Để đánh giá được nghiệm này, ta cần phải tính toán nghiệm
trong một khoảng thời gian dài hơn, như trường hợp sau đây:
>> [t,x] = ode45('eqx2',[4*pi 20*pi],[0,0]);
>> plot(t,x(:,1)),xlabel('t')
Kết quả được đưa ra như trên hình 6-14.
