131
Đưa vào thời gian cuối final_time, cùng với điều kiện đầu y(0):
>> final_time = 0.1; y0 = 0;
Sau đó gọi lệnh ode45 ñể giải theo mẫu cú pháp:
>> [t,y] = ode45('eq2',[0:0.0005:final_time],y0);
>> plot(t,y,'k-'), xlabel('t'), ylabel('y(t)'), grid on
Kết quả được thể hiện như trên hình 6-3.
0
0.02
0.04
0.06
0.08
0.1
- 4
- 2
0
2
x 10
-5
t
y
(t
)
Hì
nh 6-3. Nghiệm số với tc = 0.01 và w = 628
0
0.02
0.04
0.06
0.08
0.1
-5
0
5
x 10
-3
t
F
(t
)
0
0.02
0.04
0.06
0.08
0.1
-4
-2
0
2
x 10
-5
t
y
(t
)
Hì
nh 6-4. So sánh kích động F(t) và đáp ứng y(t)
Để so sánh với hàm kích động
( )
F t
, ta viết vào dòng lệnh:
>> F = t.*exp(-t./tc).*cos(w*t);
>> subplot(2,1,1); plot(t,F,'k-'), xlabel('t'), ylabel('F(t)'),
