132
>> grid on
>> subplot(2,1,2); plot(t,y,'k-'), xlabel('t'), ylabel('y(t)'),
>> grid on
Và nhận được đồ thị như trên hình 6-4.
Từ hình trên ta thấy lực kích động lớn hơn ñến cả 100 lần đáp ứng, nhưng chúng có
dạng khá giống nhau. Cả hai ñều có dạng dao động tắt dần.
Bây giờ ta thực hiện giải phương trình vi phân đã cho với các số liệu:
>> tc = 0.1; w = 6.28;
Bởi vì hằng số thời gian tăng 10 lần nên, ta cho thời gian tích phân cũng tăng lên
10 lần:
>> final_time = 1.0;
Gọi lại lệnh
>> [t,y] = ode45('eq2',[0 fi nal_time],y0);
>> plot(t,y,'k-'), xlabel('t'), ylabel('y(t)'), grid on
Kết quả trong trường hợp này được thể hiện trên hình 6-5 và 4-6.
0
0.2
0.4
0.6
0.8
1
- 0.01
0
0.01
0.02
t
y
(t
)
Hì
nh 6-5. Nghiệm số với tc = 0.1 và w = 6.28
0
0.2
0.4
0.6
0.8
1
- 0.05
0
0.05
t
F
(t
)
0
0.2
0.4
0.6
0.8
1
- 0.02
0
0.02
t
y
(t
)
Hì
nh 6-6. So sánh kích động F(t) và đáp ứng y(t) với tc = 0.1 và w = 6.28
