Chương 8: PHÂN TÍCH CHUỖI TUẦN TỰ THEO THỜI
GIAN VÀ DỰ BÁO
( Time series Analysis and Forecasting)
8.1 Chuỗi tuần tự theo thời gian
( Time series)
8.1.1 Đònh nghóa:
- Chuỗi tuần tự thao thời gian là một chuỗi các gía trò của một đại lượng nào
đó được ghi nhận tuần tự theo thời gian.
Ví dụ:
• Số lượng hàng bán được trong 12 tháng của một công ty.
- Các gía trò của chuỗi tuần tự theo thời gian của đại lượng X được ký hiệu X
1,
X
2
, ………, X
t
, …. X
n
, với X
t
, là gía trò quan sát của X ở thời điểm t.
8.1.2 Các thành phần của chuỗi tuần tự theo thời gian: (Components of
time series)
Các thành thống kê thường chia chuỗi tuần tự theo thời gian ra làm 4 thành
phần:
- Thành phần xu hướng dài hạn (long-term teend component)
- Thành phần mùa (Seasonal component)
- Thành phần chu kỳ ( Cyclical component)
8.1.2.3 Thành phần chu kỳ:
Thành phần này chỉ thay đổi của đại lượng X theo chu kỳ. Sự khác biệt của
thành phần này so với thành phần mùa là chu kỳ của nó dài hơn một năm. Để đánh gía
thành phần chu kỳ các gía trò của chuỗi tuần tự theo thời gian sẽ được quan sát hằng
năm
Ví dụ:
Lượng dòng chảy đến hồ chứa Trò An từ năm 1959 đến 1985 t
Chất
đốt
Xu hướng tăng theo thời gian
Thay đổi theo mùa
6
12 6 12 6 12
t (năm)
Xt: gía trò của đại lượng X ở thời điểm t
Tt, St, Ct, It: Gía trò của thành phần xu hướng, mùa, chu kỳ và bất thường ở thời
điểm t.
Trong thực tế, để dự báo gía trò của đại lượng X ta có thểa phối hợp cả 2 loại
mô hình trên.
8.2. Các phương pháp làm trơn: (Smoothing methods)
Trong một số chõi tuần tự theo thời gian thành phần mùa và thành phần bất
thường thay đổi quá lớn làm cho việc xác đònh thành phần xu hướng và thành phần chu
kỳ gặp nhiều khó khăn. Sự thay đổi lớn này có thể được giảm nhỏ bằng các phương
pháp làm trơn. Các phương pháp làm trơn được đề cập trong chương này gồm phương
pháp trung bình dòch chuyển và phương pháp làm trơn bằng hàm số mũ. (Moving
average and exponential smoothing methods)
8.2.1 Phương pháp trung bình dòch chuyển: (Trung bình trượt – Moving
average)
X
t
= T
t
* S
t
* C
t
* I
t
Xt = Tt * St *Ct
Xt = Tt + St + Ct + It
Nội dung của phương pháp này là thay thế gía trò quan sát
Xt bằng gíatrò trung bình của chính nó với m gía trò trước nó và m gía trò sau nó.
Nghóa là thay Xt bằng X
+
−
Nếu t = 3
X
*
3
=
5
XXXXX
5
4321
+
++
+
8.2.2 Phương pháp làm trơn bằng hàm số mũ đơn giản:
(Simple exponential smoothing method)
Phương pháp làm trơn hàm số mũ tiến hành dựa trên việc xem xét một cách
liên tục các gía trò của quá khứ, dựa trên trung bình có trọng số của chuỗi dữ liệu.
Trong phương pháp này trọng số có gía trò càng nhỏ khi nó càng cách xa thời điểm dự
báo. Với ý nghóa trên ta có:
t
X
= α X
t
+ (1-2) α X
t-1
3
2t
2
1t1t
X)1( X)1(X)1(X)1(X)1( αα−++αα−+αα−+α−=α−
−
−−−−
Từ các phương trình trên ta có:
Trong tóm tắt ta lấy
X
t
*
=
∑
−=
+
+
m
mj
jt
X
1m2
1
XX =
Và các trò số X
t
được lấy từ số liệu quan sát
Hàm 2
Cho (1-α)
122
XXX +α=
233
XXX +α=
. . . . . . . . . . . . .
1ttt
X)1(XX
−
α−+α=
. . . . . . . . . . . . . .
1nnn
X)1(XX
−
α−+α=
Ví dụ:
Bảng số liệu sau đây cho thấy số lượng máy tính được bán ra của cửa hàng
trong 24 tháng vừa qua.
a) Tìm chuỗi thời gian trung bình trượt 5 tháng
b) Tìm các gía trò trung bình bằng phương pháp làm trơn hàm mũ với hệ số làm trơn
STT
t
Tháng Số liệu quan sát
19
20
21
22
23
24
2 năm trước
Một
Hai
Ba
Tư
Năm
Sáu
Bảy
Tám
Chín
Mười
Mười Một
Mười Hai
1 năm trước
Môt
Hai
Ba
Tư
Năm
Sáu
Bảy
Tám
18.4
17.2
17.6
19.4
20.8
23.0
25.0
24.6
24.0
24.0
23.8
24.4
27.6
29.0
29.0
28.8
27.6
25.0
23.0
23.4
21.0
20.0
19.8
Đoán trước sự vận động của đối tượng nghiên cứu trong tương lai. Đó là kết quả
nhận thức chủ quan của con người dựa trên một số cơ sở nhất đònh. Có thể nêu mức độ
tiên đoán ở 3 khía cạnh.
a) Tiên đoán không tưởng:
Đó là những tiên đoán không có cơ sở khoa học, chỉ dựa trên những mối liên hệ
không tưởng thiếu căn cứ.
b) Tiên đoán kinh nghiệm:
Đó là những tiên đoán dựa trên chuỗi thông tin lòch sử. Mức độ ít nhiều có cơ sở
khách quan, tuy nhiên có nhược điểm là loại tiên đoán này không giải thích được xu
thế vận động của đối tượng nghiên cứu và đa số dừng lại ở bước đònh tính.
c) Tiên đoán khoa học:
Đó là tiên đoán dựa trên phân tích mối liên hệ qua lại giữa các đối tượng
nghiên cứu và phương pháp xử lý thông tin khoa học nhằm phát hiện tính quy luật của
đối tượng .
8.3.1.1.2 Dự báo (Forecasting)
Dự báo là tiên đoán khoa học mang tính xác suất và tính phương án trong
khoảng thời gian hữu hạn về tương lai phát triển của đối tượng nghiên cứu.
a) Tính xác suất :
Do dự báo dựa trên việc xử lý chuỗi thông tin bao gồm cả 2 yếu tố xu thế phát
triển và yếu tố ngẫu nhiên, vì vậy kết quả dự báo so với thực tế có sự chên lệch mang
tính xác suất.
b) Tính chất phương án:
Dự báo được thể hiện bằng nhiều dạng kết quả có thể xảy ra trong tương lai
(dạng đònh tính, dạng đònh lượng, dạng khoảng, dạng điểm, …)
c) Tính chất thời gian hữu hạn:
Sự chênh lệch giữa thời điểm dự báo và thời điểm hiện tại được gọi là khoảng
cách dự báo (tầm xa dự báo l), khoảng cách này không thể tùy tiện mà nó phụ thuộc
vào mức độ ổn đònh của đối tượng nghiên cứu trong quá trình phát triển . Vì vậy dự
báo được tiến hành với khoảng cách dự báo thích hợp tương ứng một khoảng thời gian
hữu hạn nào đó.
Yêu cầu phải có số liệu
- Chính xác
- Đúng mục đích dự báo
F
t+l
= A
F
t+l
= A
ε
±
Đây là phần khó khăn và tốn thời gian.
8.3.2.2 Xử lý sơ bộ số liệu:
- Bỏ những số liệu không cần thiết, không chính xác.
- Bổ xung những số liệu còn thiếu.
- Chia tập số liệu thành 2 nhóm : nhóm đầu và nhóm kiểm tra.
8.3.2.3 Lựa chọn phương pháp và xây dựng mô hình dự báo
- Phương pháp dự báo được chọn sao cho phù hợp với số kiệu thuộc nhóm
đầu và với đối tượng nghiên cứu.
- Lập mô hình dự báo sao cho sai số dự báo là nhỏ nhất. Sai số dự báo được
kiểm đònh bởi nhóm số liệu kiểm tra.
8.3.2.4 Dự báo:
- Từ mô hình dự báo xác đònh giá trò dự báo
- Phân tích kết quả nhận được.
8.3.3 Các phương pháp dự báo thông dụng:
Các phương pháp dự báo thường dùng có thể chia làm 3 loại:
- Phương pháp giản đơn (naive method)
– F
t+1
Ví dụ:
Lượng hàng bán ra theo qúy của một công ty từ năm 1979 đến 1985 được cho
trong bảng số liệu.
Dựa vào bảng số liệu lập mô hình dự báo lượng hàng bán ra theo quý trong tương
lai.
Năm Quý T Lương sp
bán được
1979
1980
1981
1982
1
2
3
4
1
2
3
4
1
2
3
4
1
2
3
4
5
6
7
8
9
10
11
12
200
150
400
550
350
250
550
550
400
350
600
750
500
400
650
850
600
450
700
Giải:
Nếu chúng ta sử dụng số liệu từ năm 1979 đến năm 1984 làm phần đầu của
nhóm số liệu và sử dụng số liệu năm 1985 là phần kiểm tra
Giá trò dự báo chko q I năm 1985 (thời điểm 25)
F
25
Trong phương pháp này, người ta thêm vào ác số hạng để dánh giá ảnh
hưởng của thành phần xu hướng trong chuỗi thời gian. Mô hình có thể có dạng:
a)
Đánh giá ảnh hưởng của thành phần xu hướng
Ví dụ:
Giá trò dự báo ở thời điểm 25
F
25
= X
24
+ ( X
24
- X
23
)
= 650 + ( 650 – 400)
F25 = 900
Sai số dự báo ở thời điểm 25
E25 = X25 - F25
= 850 - 900
e25 = -50
b)
F
t+1
= X
t
+ ( X
t
F
25
= 1050
e
25
= X
25
- F
25
= 850 -1056
e
25
= -206
C. Phương pháp 3:
Trong phương pháp này, người ta có chú ý đến ảnh hưởng của thành
phần mùa trong chuỗi thời gian. Với công thức này lượng sản phẩm bán ra ở quý dự báo bằng lượng sản
phẩm bán ra trong quý tương ứng ở năm trước đó. Nhược điểm của phương pháp này là
không chú ý tới tác động khác nhau giữa các năm và thành phần xu hướng.
Ví dụ:
F
25
= X
21
= 750
21
+
4
)XX()XX)XX()XX(
2021212222232324
−
+
−
+
−
+
−
= 750 +
4
)600750()750500()500400()400650( −+−+−+−
= 750 +12,5
F
25
= 762,5
e
25
= X
25
- F
25
= 850 - 762.5
e
F
t+1
=
n
1
∑
=
n
1t
t
X
F
25
=
24
1
∑
=
24
1t
t
X
=
)9800(*
24
1
F
t
25
1t
t
XX
25
1
X
25
1
=
[]
8509800
25
1
+
F
26
=
25
1
* 10650 = 426
e
26
= X
26
- F
26
= 600 - 426
• n = 1⇒ F
t+1
= X
t
: Phương pháp nạve
• n = 1⇒ F
t+1
=
∑
=
n
1t
t
X
n
1
:Phương pháp trung bình giản đơn.
- Tổng quát, phương pháp trung bình dich chuyển tốt hơn phương pháp trung bình
giản đơn.
Ví dụ:
Giá trò dự báo cho q I/1985 với n = 4
F
25
=
4
1
(X
24
+X
F
t+1
:giá trò dự báo ở thời điểm t+1
F
t
:giá trò dự báo ở thời điểm t
X
t
:giá trò quan sát ở thời điểm t
α :hằng số làm tròn , 10
≤
α
≤
(toàn bộ số
liệu)
F
t+1
= α X
t
+ (1-
α
) F
t
t
= β +Φ
1
X
t-1
+ Φ
2
X
t-2
+ a
t
X
t
=
β
+ Φ
1
X
t-1
+
Φ
2
X
t-2
+ … +
Φ
p
X
t-p
+a