70
Khi chạy chương trình chính các chương trình sẽ trao đổi dữ liệu theo sơ đồ sau:
% file_name = ptvp_vantoc.m
function dydt=ptvp_vantoc(t,y)
global m b F0
% khai bao cac bien tong the
dydt=1/m*(F0-b*y);
% file_name = Euler_method.m
function [t, z]=Euler_method(ptvp, t0, tf, y0, N)
. . . . . . . . . .
% Chuong trinh chinh
global m b F0
% khai bao cac bien tong the
% Gan gia tri cho cac tham so
m=...; b=...; ...
[t, z]=Euler_method(@ptvp_vantoc, t0, tf, y0, N);
@ptvp_vantoc, t0, tf, y0, N
t, y
t, y
dy/dt
Khi chạy chương trình chính cho ta các kết quả như trên đồ thị dưới đây.
0
1
2
3
4
5
6
7
8
9
10
0
5
10
15
20
y, z [m/s]
0
1
2
3
4
5
6
7
8
9
10
0
0.1
0.2
0.3
0.4
t [s]
|y−z|
y
z
Hình 1. Đồ thị các hàm y(t), z(t) và hiệu y(t)-z(t)
