Tiểu mục này trình bày một số kết quả thực nghiệm thu được từ mô hình
FTS2NT-CMPSO trên hai tập dữ liệu khác nhau: dữ liệu về nhiệt độ trung bình hàng ngày tại Đài Bắc, Đài Loan trong Bảng 2.7 và thị trường chứng khoán Đài Loan
TAIFEX [54]. Các kết quả thu được từ mô hình FTS2NT-CMPSO được đưa ra để đánh giá, so sánh với các mô hình dự báo khác. Ngoài ra, để khẳng định tính hiệu quả của mô hình dựa trên quan hệ mờ bậc cao hai nhân tố, mô hình FTS2NT- CMPSO cũng được xem xét để so sánh với mô hình một nhân tố FTS1NT- CMPSO đã được đề cập trong Mục 3.1.1 trên cùng tập dữ liệu. Trước khi thực
nghiệm mô hình dự báo được đề xuất, các tham số của mô hình sử dụng tối ưu PSO như: số lượng cá thể, số lần lặp tối đa, giới hạn vận tốc và vị trí của mỗi cá thể cần được thiết lập để đảm bảo các rằng buộc của hàm tối ưu cũng như miền giới hạn tìm kiếm trong PSO. Tùy thuộc vào mỗi tập dữ liệu thì các tham số này có thể được thiết lập với các giá trị khác nhau, như được chỉ ra trong các Bảng 3.24 và Bảng 3.25.
Bảng 3.24: Các tham số của mô hình FTS2NT-CMPSO được thiết lập cho dự báo nhiệt độ
Các tham số của mô hình Giá trị thiết lập cho dự báo nhiệt độ
Nhiệt độ Mật độ của mây Số lượng cá thể trong quần thể: pn = 30
Số lần lặp tối đa (số thế hệ): iter_max = 150
Phạm vi của bậc cần tối ưu m = 1-8
Trọng số quán tính ω (Giảm tuyến tính) 𝒜𝒜𝒜𝒜= 1.4 to 𝒜𝒜 𝒜𝒜 = 0.4
Hệ số tự tin cậy và hệ số xã hội C1 = C2 2.05
Phạm vi vị trí của mỗi cá thể: X = [23, 32] [0, 100]
Phạm vi vận tốc của mỗi cá thể: V = [-5, 5] [-50, 50] Bảng 3.25: Các tham số của mô hình FTS2NT-CMPSO được thiết lập cho dự báo
thị trường chứng khoán
Các tham số của mô hình Giá trị thiết lập cho dự thị trường chứng khoán
TAIFEX TAIEX
Số lượng cá thể trong quần thể: pn = 30
Số lần lặp tối đa (số thế hệ): iter_max = 150
Phạm vi của bậc cần tối ưu m = 1-8
Trọng số quán tính ω (Giảm tuyến tính) 𝒜𝒜𝒜𝒜= 1.4 to 𝒜𝒜𝒜𝒜 = 0.4
Hệ số tự tin cậy và hệ số xã hội C1 = C2 2.05
Miền giới hạn vị trí của mỗi cá thể: X = [6200, 7560] [6251, 7599]
Miền giới hạn vận tốc của mỗi cá thể: V = [-750, 750] [-750, 750] 3.2.2.1 Áp dụng dự báo trên tập dữ liệu nhiệt độ
Từ các tham số trong Bảng 3.24, thực hiện một số lần chạy mô hình
FTS2NT- CMPSO trên từng tháng của tập dữ liệu nhiệt độ hàng ngày. Kết quả của
lần chạy tương ứng với khoảng tốt nhất cho tất cả các bậc của mô hình và chọn bậc có kết quả dự báo với sai số nhỏ nhất làm kết quả dự báo cuối cùng. Sau khi kết thúc của mỗi
lần chạy, mô hình FTS2NT-CMPSO sẽ hiển thị độ dài khoảng tối ưu và sai số dự báo của mỗi bậc (từ bậc 1 đến bậc 8). Trong phần thực nghiệm này, luận án sử dụng hàm MAPE (1.10) làm tiêu chí đánh giá về sai số dự báo để tiện so sánh kết quả thực nghiệm với các mô hình khác cùng sử dụng tiêu chí này.
Kết quả dự báo nhiệt độ hàng ngày từ tháng 6/1996 đến tháng 9/ 1996 của mô hình FTS2NT-CMPSO dựa trên bậc tốt nhất và độ dài khoảng tối ưu với 9 khoảng chia cho nhân tố chính “nhiệt độ trung bình” và 7 khoảng chia cho nhân tố thứ hai “mật độ của mây” được thể hiện trong Bảng 3.26. Thêm nữa sai số dự báo MAPE khi thực hiện tối ưu các bậc của mô hình cũng được ghi lại trong Bảng 3.27.
Bảng 3.26: Kết quả dự báo nhiệt độ trung bình hàng ngày tại Đài Bắc, Đài Loan theo từng tháng từ ngày 01/06/1996 đến 30/09/1996. Ngày, Tháng (T) Nhiệt độ_T6 Dự báo_T6 Nhiệt độ_T7 Dự báo_T7 Nhiệt độ_T8 Dự báo_T8 Nhiệt độ_T9 Dự báo_T9 1 26.1 -- 29.9 -- 27.1 -- 27.5 -- 2 27.6 -- 28.4 -- 28.9 -- 26.8 -- 3 29 -- 29.2 -- 28.9 -- 26.4 -- 4 30.5 -- 29.4 -- 29.3 29.3 27.5 -- 5 30 -- 29.9 29.9 28.8 28.8 26.6 -- 6 29.5 -- 29.6 29.6 28.7 28.7 28.2 28.21 7 29.7 -- 30.1 30.1 29 29 29.2 29.2 8 29.4 29.4 29.3 29.31 28.2 28.21 29 29 9 28.8 28.8 28.1 28.1 27 27 30.3 30.3 10 29.4 29.4 28.9 28.9 28.3 28.3 29.9 29.9 11 29.3 29.32 28.4 28.39 28.9 28.9 29.9 29.9 12 28.5 28.49 29.6 29.6 28.1 28.11 30.5 30.47 13 28.7 28.7 27.8 27.8 29.9 29.89 30.2 30.22 14 27.5 27.5 29.1 29.09 27.6 27.59 30.3 30.3 15 29.5 29.5 27.7 27.7 26.8 26.8 29.5 29.48 16 28.8 28.8 28.1 28.1 27.6 27.59 28.3 28.3 17 29 29 28.7 28.71 27.9 27.9 28.6 28.6 18 30.3 30.31 29.9 29.9 29 29 28.1 28.11 19 30.2 30.19 30.8 30.8 29.2 29.2 28.4 28.4 20 30.9 30.9 31.6 31.57 29.8 29.81 28.3 28.3 21 30.8 30.81 31.4 31.39 29.6 29.58 26.4 26.4 22 28.7 28.7 31.3 31.3 29.3 29.3 25.7 25.79 23 27.8 27.81 31.3 31.3 28 27.99 25 25.01 24 27.4 27.4 31.3 31.3 28.3 28.3 27 26.94 25 27.7 27.69 28.9 28.9 28.6 28.59 25.8 25.79 26 27.1 27.1 28 28 28.7 28.7 26.4 26.4 27 28.4 28.38 28.6 28.61 29 29 25.6 25.6 28 27.8 27.81 28 28 27.7 27.71 24.2 24.17
29 29 29 29.3 29.31 26.2 26.22 23.3 23.34
30 30.2 30.19 27.9 27.9 26 26.07 23.5 23.52
31 26.9 26.94 27.7 27.71
MAPE 0.018% 0.018% 0.027% 0.053%
Bảng 3.26 thể hiện kết quả dự báo trên từng tháng của tập dữ liệu nhiệt độ từ tháng 6 đến tháng 9 năm 1996. Quan sát bảng này cho thấy, các giá trị dự báo thu được trên từng tháng của tập dữ liệu nhiệt độ bám rất sát so với dữ liệu thực thế dựa trên các bậc tối ưu đạt được từ mô hình. Từ những kết quả này có thể nhận thấy việc kết hợp giữa kỹ thuật tối ưu PSO và phương pháp phân khoảng FCM trong mô hình chuỗi thời gian mờ hai nhân tố FTS2NT-CMPSO là phù hợp và đáng tin cậy.
Bảng 3.27: Sai số dự báo MAPE (%) trên từng bậc của mô hình FTS2NT-CMPSO Tháng Bậc 1 Bậ 2 Bậc 3 Bậc 4 Bậc 5 Bậc 6 Bậc 7 Bậc 8
T6 1.485 0.117 0.031 0.025 0.021 0.021 0.018 0.019
T7 1.64 0.023 0.02 0.018 0.019 0.019 0.02 0.02
T8 0.92 0.038 0.027 0.028 0.029 0.031 0.032 0.032
T9 1.302 0.258 0.262 0.061 0.053 0.054 0.057 0.059 Ngoài ra Bảng 3.27 cũng đưa ra độ chính xác trên từng bậc của mô hình và lựa chọn được bậc với độ chính xác tốt nhất khi sử dụng kỹ thuật tối ưu PSO. Điều này có thể thấy việc đề xuất cải tiến thuật toán sử dụng PSO để tối ưu đồng thời bậc của mô hình và độ dài khoảng là rất khả quan và hữu hiệu trong mô hình chuỗi thời mờ hai nhân tố. Mô tả trực quan hơn, có thể thấy đường cong biểu diễn giữa giá trị dự báo với giá trị thực tế trong tháng 6 như được minh họa trên Hình 3.9.
Hình 3.9: Đường cong biểu diễn giữa giá trị thực tế và giá trị dự báo nhiệt độ trung bình hàng ngày trong tháng 6, dựa trên mô hình FTS2NT-CMPSO bậc 7
Tiếp theo, để chứng tỏ tính ưu việt của mô hình FTS2NT-CMPSO dựa trên chuỗi thời gian mờ bậc cao hai nhân tố với số lượng khoảng chia bằng 9 cho nhân tố chính “nhiệt độ” và bằng 7 cho nhân tố thứ hai “mật độ của mây”, các mô hình trong các
bài báo [17, 34, 40, 54, 56, 57] được lựa chọn cho việc so sánh. Các kết quả về sai số dự báo MAPE trên từng tháng (từ tháng 6 đến tháng 9 năm 1996) giữa mô hình
FTS2NT-CMPSO với các mô hình so sánh dựa trên các bậc khác nhau (từ bậc 2
đến bậc 8) được đưa ra từ Bảng 3.28 đến Bảng 3.31 như sau:
Bảng 3.28: So sánh sai số dự báo MAPE (%) trong tháng 6 giữa mô hình FTS2NT-
CMPSO và các mô hình khác dựa trên các bậc khác nhau
Mô hình Bậc của mô hình Bậc 2 Bậc 3 Bậc 4 Bậc 5 Bậc 6 Bậc 7 Bậc 8 [56] 0.28 0.29 0.30 0.29 0.29 0.28 0.29 [54] 0.8 0.76 0.79 0.76 0.79 0.79 0.81 α=0.25 0.44 0.42 0.42 0.42 0.44 0.40 0.40 [34] 𝒜 =0.5 0.50 0.45 0.42 0.38 0.43 0.39 0.46 𝒜 =0.9 0.46 0.42 0.44 0.42 0.41 0.46 0.39 [40] 0.36 0.34 0.32 0.31 0.31 0.28 0.29 [17] 0.25 0.26 0.27 0.27 0.27 0.28 0.29 [57] 0.27 0.28 0.27 0.25 0.25 0.18 0.18 FTS2NT-CMPSO 0.117 0.031 0.025 0.021 0.021 0.018 0.019
Bảng 3.29: So sánh sai số dự báo MAPE (%) trong tháng 7 giữa mô hình FTS2NT-
CMPSO và các mô hình khác dựa trên các bậc khác nhau
Mô hình Bậc của mô hình
Bậc 2 Bậc 3 Bậc 4 Bậc 5 Bậc 6 Bậc 7 Bậc 8 [56] 0.34 0.35 0.34 0.34 0.35 0.33 0.32 [54] 0.96 0.96 0.98 0.97 1.00 0.98 0.99 α=0.25 0.42 0.42 0.44 0.40 0.40 0.40 0.40 [34] 𝒜 =0.5 0.42 0.38 0.43 0.39 0.46 0.46 0.46 𝒜 =0.9 0.44 0.42 0.41 0.46 0.39 0.39 0.39 [40] 0.34 0.33 0.33 0.32 0.32 0.34 0.33 [17] 0.28 0.29 0.28 0.29 0.30 0.29 0.30 [57] 0.23 0.23 0.23 0.23 0.23 0.22 0.22 FTS2NT-CMPSO 0.023 0.02 0.018 0.019 0.019 0.02 0.02
Bảng 3.30: So sánh sai số dự báo MAPE (%) trong tháng 8 giữa mô hình FTS2NT-
CMPSO và các mô hình khác dựa trên các bậc khác nhau.
Mô hình Bậc của mô hình
Bậc 2 Bậc 3 Bậc 4 Bậc 5 Bậc 6 Bậc 7 Bậc 8 [56] 0.23 0.22 0.22 0.22 0.23 0.23 0.22 [54] 1.07 1.06 1.08 1.08 1.09 1.07 1.07 α=0.25 0.42 0.42 0.44 0.40 0.40 0.40 0.40 [34] 𝒜 =0.5 0.42 0.38 0.43 0.39 0.46 0.46 0.46 𝒜 =0.9 0.44 0.42 0.41 0.46 0.39 0.39 0.39 [40] 0.31 0.34 0.33 0.33 0.33 0.34 0.35
[17] 0.35 0.35 0.35 0.36 0.37 0.35 0.36
[57] 0.12 0.12 0.13 0.13 0.17 0.13 0.13
FTS2NT-CMPSO 0.038 0.027 0.028 0.029 0.031 0.032 0.032
Bảng 3.31: So sánh sai số dự báo MAPE (%) trong tháng 9 giữa mô hình FTS2NT-
CMPSO và các mô hình khác dựa trên các bậc khác nhau
Mô hình Bậc của mô hình
Bậc 2 Bậc 3 Bậc 4 Bậc 5 Bậc 6 Bậc 7 Bậc 8 [56] 0.51 0.49 0.51 0.51 0.53 0.50 0.51 [54] 1.01 0.90 0.94 0.96 0.95 0.95 0.95 α=0.25 0.42 0.42 0.44 0.40 0.40 0.40 0.40 [34] 𝒜 =0.5 0.42 0.38 0.43 0.39 0.46 0.46 0.46 𝒜 =0.9 0.44 0.42 0.41 0.46 0.39 0.39 0.39 [40] 0.54 0.56 0.54 0.50 0.51 0.52 0.41 [17] 0.36 0.35 0.35 0.37 0.38 0.37 0.38 [57] 0.27 0.28 0.28 0.29 0.30 0.30 0.29 FTS2NT-CMPSO 0.258 0.262 0.061 0.053 0.054 0.057 0.059
Quan sát sai số dự báo MAPE trên từng bậc của các mô hình dự báo từ Bảng 3.28 đến Bảng 3.31, cho thấy mô hình dự báo FTS2NT-CMPSO đưa ra sai số nhỏ hơn so với các mô hình so sánh tương ứng với từng tháng trong tập dữ liệu nhiệt độ trung bình hàng ngày tại Đài Bắc, Đài Loan. Cụ thể mô hình FTS2NT-CMPSO đạt
giá trị MAPE nhỏ hơn nhiều so với các mô hình dự báo so sánh trên tất cả các bậc của mô hình. Đặc biệt, đưa ra sai số dự báo tốt nhất tương ứng với từng tháng dựa trên các bậc tối ưu lần lượt là: kết quả dự báo tháng 6 đạt sai số với giá trị MAPE =
0.018% nhỏ nhất dựa trên bậc 7, tháng 7 đưa ra sai số MAPE = 0.018% nhỏ nhất ở
bậc 4, tháng 8 đưa ra sai số MAPE = 0.027% nhỏ nhất ở bậc 3 và tháng 9 đạt giá trị
MAPE = 0.053% nhỏ nhất trong trường hợp bậc 5. Từ kết quả này cho thấy mô
hình đề xuất FTS2NT-CMPSO ưu việt hơn một số mô hình trước đây trên tập dữ liệu nhiệt độ trung bình hàng ngày tại Đài Bắc, Đài Loan.
3.2.2.2 Áp dụng dự báo trên tập dữ liệu thị trường chứng khoán.
Trong tiểu mục này, mô hình dự báo FTS2NT-CMPSO được áp dụng để dự báo thị trường chứng khoán Đài Loan [54] từ ngày 8/3/1998 đến 9/30/1998 bao gồm hai nhân tố có các chỉ số tương ứng là TAIFEX và chỉ số TAIEX.
Để chứng tỏ tính hiệu quả dự báo của mô hình FTS2NT-CMPSO dựa trên chuỗi thời gian bậc cao hai nhân tố, các mô hình có tên như: L06 [54], L08 [34], MPTSO [40], THPSO [35] và các mô hình trong [17, 56, 91] được lựa chọn để so sánh và đánh giá về sai số dự báo MSE (1.8). Từ các tham số được lựa chọn trong Bảng 3.32, mô hình FTS2NT-CMPSO được thử nghiệm với số khoảng là 16 cho cả
hai nhân tố và đưa ra kết quả dự báo với bậc có giá trị MSE nhỏ nhất.
Ngoài ra, khi thực hiện tối ưu bậc và độ dài khoảng một cách đồng thời, mô hình đề xuất cũng đưa ra thời gian tính toán và tốc độ hội tụ theo giá trị MAPE với số khoảng cố định bằng 16. Đồ thị biểu diễn tốc độ hội tụ theo giá trị MAPE với 150 lần lặp trên dữ liệu TAIFEX được minh họa trong Hình 3.10. Trong thử nghiệm này, mô hình hội tụ tới một giá trị MAPE thấp nhất có thể chấp nhận được là (MAPE =
0.23% hay MSE = 4.8) tại 90 lần lặp.
Kết quả và sai số dự báo tương ứng với bậc tối ưu của mô hình được đề xuất
FTS2NT-CMPSO và các mô hình so sánh được đưa ra trong Bảng 3.32.
Hình 3.10: Đồ thị thể hiện thời gian tính toán và tốc độ hội tụ của mô hình FTS2NT-CMPSO với 150 lần lặp
Bảng 3.32: So sánh kết quả và sai số dự báo MSE giữa mô hình FTS2NT-CMPSO với các mô hình khác dựa trên QHM bậc cao với 16 khoảng chia cho mỗi nhân tố
Ngày tháng DL
thực L06[54] [17] [56] L08[34] MTPSO[40] [91] THPSO[35] FTS2NT-CMPSO
8/3/1998 7552 -- -- -- -- -- -- -- -- 8/4/1998 7560 -- -- -- -- -- -- -- -- 8/5/1998 7487 -- -- -- -- -- -- -- -- 8/6/1998 7462 7450 7456.63 7474.5 -- -- -- -- -- 8/7/1998 7515 7550 7515.13 7515.8 -- -- -- -- 7514.84 8/10/1998 7365 7350 7381.13 7413.5 -- -- -- -- 7365.88 8/11/1998 7360 7350 7381.13 7362.5 -- -- -- -- 7356.89 8/12/1998 7330 7350 7325 7345 7329 7325.28 7325.5 7325 7329.92 8/13/1998 7291 7250 7280.33 7295.5 7289.5 7287.48 7292.5 7287.5 7292.64 9/29/1998 6806 6850 6781.93 6833.5 6796 6781.01 6796.6 6794.3 6805.86 9/30/1998 6787 6750 6781.93 6788.5 6796 6781.01 6796.6 6794.3 6786.9 10/01/1998 N/A 6750.5 MSE 1364.56 287.32 252.47 105.02 92.17 69.98 55.96 4.8
Từ kết quả trong Bảng 3.32 cho thấy, mô hình FTS2NT-CMPSO đưa ra giá trị dự báo trên từng ngày gần như sát với dữ thực tế và đạt sai số dự báo với giá trị (MSE = 4.8) nhỏ hơn nhiều so với các mô hình so sánh dựa trên QHM bậc 4.
Chi tiết đánh giá dựa trên ba điểm khác chính là cách NQHM, phương pháp phân khoảng tập nền và quy tắc giải mờ đầu ra. Kết quả tính toán cho thấy ưu thế của mô hình FTS2NT-CMPSO được thể hiện trong việc sử dụng dụng nhóm quan hệ mờ phụ thuộc thời gian, trong khi các mô hình còn lại đều sử dụng NQHM của Chen [10]. Với công việc phân khoảng tập nền, thì hai mô hình trong [17, 57] sử dụng phân cụm tự động để tìm khoảng tối ưu, mô hình [91] sử dụng tìm kiếm Tuba với hệ thống suy diễn mờ để điều chỉnh độ dài khoảng chia tập nền. Còn lại các mô hình L06 [54], L08 [34], MTPSO [40], THPSO [35] và mô hình FTS2NT-CMPSO đều áp dụng thuật toán tối ưu để điều chỉnh khoảng ban đầu và tìm khoảng chia tối ưu trên mỗi tập nền của từng nhân tố, trong đó hai mô hình trước là sử dụng GA còn 3 mô hình còn lại sử dụng PSO.
Bên cạnh đó, để thấy được hiệu quả của việc mở rộng mô hình chuỗi thời gian mờ một nhân tố thành hai nhân tố, mô hình FTS2NT-CMPSO được so sánh
với mô hình chuỗi thời gian mờ một nhân tố (FTS1NT-CMPSO) trên tập dữ liệu thị trường chứng khoán TAIFEX. Quan sát Bảng 3.32 cho thấy sai số dự báo của mô hình hai nhân tố FTS2NT-CMPSO với giá trị MSE = 4.8 nhỏ hơn so với mô