130
0
1
2
3
4
5
6
0
1
2
3
4
t
y
(t
)
Hì
nh 6-2. Nghiệm số và nghiệm giải tích của ptvp
/
cos( )
dy dt
t
=
, ode45
Lưu ý rằng số điểm chia sử dụng trong ode23 nhỏ hơn số điểm chia khi sử dụng
ode45. Hãy so sánh length(y) va length(w).
Trong ví dụ sau ta sẽ sử dụng ode45 ñể giải phương trình vi phân
5
5 ( )
dy
y
F t
dt
= −
+
, với kích động
/
( )
cos( )
t tc
F t
te
t
ω
−
=
,
với các số liệu cho trong bảng
Hằng số thời gian,
c
t
Tần số
ω
0.01 s
628 rad/s
0.1 s
6.28 rad/s
Trước hết định nghĩa phương trình vi phân trong một m-file, trong đó các thông số
c
t
và ω được khai báo ở dạng biến toàn cục (global tc w):
function ydot = eq2(t,y)
global tc w
F = t*exp(-t/tc)*cos(w*t);
ydot = 5*(F - y);
% save with file name eq2.m
Khi ghi lại phần soạn thảo trên Matlab tự động lấy tên hàm eq2 đặt tên cho m-file.
Trong của sổ lệnh, ta khai báo tc và w là các biến toàn cục
>> global tc w
Trước hết sử dụng các thông số trong dòng đầu của bảng:
>> tc = 0.01; w = 628;
