Tổng quan nghiên cứu

Trong thời đại bùng nổ của trí tuệ nhân tạo và khoa học dữ liệu, ước tính có hơn 80% khối lượng dữ liệu phát sinh tại các tổ chức kinh tế, tài chính và khí tượng là dữ liệu chuỗi thời gian. Dự báo chuỗi thời gian giữ vai trò then chốt trong việc kiểm soát rủi ro và tối ưu hóa quyết định đầu tư. Tuy nhiên, các chuỗi dữ liệu thực tế thường mang đặc tính phi tuyến tính cao, chuyển động phức tạp và chứa nhiều yếu tố nhiễu. Các mô hình thống kê cổ điển như ARIMA hay làm trơn hàm mũ thường bộc lộ nhiều hạn chế khi xử lý các không gian dữ liệu đa chiều phức tạp.

Mạng nơ-ron học sâu DBN (Deep Belief Networks), do nhà khoa học Geoffrey Hinton đề xuất vào năm 2006, đã tạo nên bước ngoặt lớn khi giải quyết triệt để hiện tượng triệt tiêu độ dốc (vanishing gradient). Dù vậy, hiệu quả dự báo của DBN lại phụ thuộc mật thiết vào 4 siêu tham số cốt lõi: hệ số học RBM, hệ số học lan truyền ngược, số nơ-ron đầu vào và số nơ-ron lớp ẩn. Việc thiết lập các tham số này theo cách thủ công hoặc ngẫu nhiên thường tiêu tốn hàng trăm giờ tính toán mà không đảm bảo đạt được tối ưu toàn cục.

Đề tài luận văn thạc sĩ chuyên ngành Khoa học Máy tính do học viên Đỗ Ngọc Lưu thực hiện dưới sự hướng dẫn của PGS.TS Dương Tuấn Anh tại Trường Đại học Bách Khoa – Đại học Quốc Gia TP. Hồ Chí Minh (hoàn thành vào tháng 12 năm 2019) tập trung giải quyết bài toán trên. Mục tiêu nghiên cứu là xây dựng mô hình dự báo chuỗi thời gian kết hợp mạng DBN và mạng Perceptron đa lớp (MLP), tích hợp giải thuật siêu phỏng đoán Harmony Search (đặc biệt là biến thể Global-best Harmony Search) để tự động ước lượng bộ siêu tham số tối ưu. Công trình được kiểm chứng trên 4 bộ dữ liệu với quy mô hàng trăm chu kỳ thời gian, giúp cải thiện độ chính xác dự báo và giảm thiểu sai số toàn phương trung bình (MSE) từ 15% đến 35% so với các phương pháp truyền thống.

Cơ sở lý thuyết và phương pháp nghiên cứu

Khung lý thuyết áp dụng

Mô hình nghiên cứu kết hợp giữa lý thuyết mạng nơ-ron học sâu và giải thuật tối ưu hóa mô phỏng sinh học. Trong đó, 2 nền tảng lý thuyết chủ đạo được tích hợp gồm:

  1. Lý thuyết mạng Deep Belief Networks (DBN) và máy Boltzmann giới hạn (RBM): RBM là mô hình sinh dựa trên năng lượng bao gồm 2 lớp: lớp quan sát (visible units) và lớp ẩn (hidden units) không có liên kết nội lớp. DBN được xây dựng bằng cách xếp chồng các lớp RBM lên nhau. Quá trình huấn luyện không giám sát (pre-training) áp dụng giải thuật Phân kỳ tương phản (Contrastive Divergence - CD) kết hợp lấy mẫu Gibbs để trích xuất đặc trưng sâu từ dữ liệu thô. Sau đó, mạng Perceptron đa lớp (MLP) ở tầng trên cùng được tinh chỉnh có giám sát bằng thuật toán lan truyền ngược (Back-Propagation - BP) để xuất giá trị dự báo.
  2. Lý thuyết giải thuật Harmony Search (HS) và Global-best Harmony Search (GHS): Được Zong Woo Geem giới thiệu năm 2001, Harmony Search mô phỏng quá trình người nhạc sĩ sáng tạo hợp âm hoàn hảo. Giải thuật vận hành dựa trên các tham số chính: kích thước bộ nhớ hòa âm (Harmony Memory Size - HMS thiết lập từ 50 đến 100), tỷ lệ xem xét bộ nhớ (HMCR) và tỷ lệ điều chỉnh cao độ (Pitch Adjusting Rate - PAR). Biến thể GHS cập nhật động tham số PAR qua từng vòng lặp, tận dụng thông tin từ giải pháp tốt nhất toàn cục để tăng tốc độ hội tụ.

Các khái niệm chuyên ngành cốt lõi được áp dụng xuyên suốt gồm: hàm năng lượng RBM, xấp xỉ phân kỳ tương phản CD-k, cơ chế nhớ hòa âm, và chỉ số sai số toàn phương trung bình (MSE).

Phương pháp nghiên cứu

  • Nguồn dữ liệu và cỡ mẫu: Nghiên cứu thu thập và phân tích 4 bộ dữ liệu mẫu đa dạng:
    • Bộ dữ liệu thực tế 1: Chuỗi đo lường trung bình vệt đen mặt trời hàng tháng từ tháng 1 năm 1749 đến tháng 3 năm 1977 với cỡ mẫu 305 mốc thời gian (288 mẫu dùng để huấn luyện và 17 mẫu kiểm thử).
    • Bộ dữ liệu thực tế 2: Chỉ số giá tiêu dùng (CPI) tại Tây Ban Nha hàng tháng từ tháng 1 năm 1960 đến tháng 6 năm 2005 với quy mô 546 mẫu quan sát.
    • Bộ dữ liệu thực tế 3: Tỷ giá trung bình giữa USD và Bảng Anh hàng tháng từ tháng 1 năm 1981 đến tháng 7 năm 2005 gồm 295 mẫu dữ liệu tài chính.
    • Bộ dữ liệu tổng hợp 4: Chuỗi thời gian phi tuyến hỗn loạn sinh ra từ hệ phương trình Lorenz với 1000 điểm dữ liệu mô phỏng.
  • Phương pháp chọn mẫu: Dữ liệu được trích xuất tuần tự theo phương pháp lấy mẫu chuỗi thời gian liên tục và cấu trúc lại qua kỹ thuật cửa sổ trượt (sliding window) với các bước trễ không gian nhúng, đảm bảo giữ nguyên tính phụ thuộc thời gian và quan hệ nhân quả của chuỗi dữ liệu.
  • Phương pháp phân tích và lý do lựa chọn: Toàn bộ hệ thống được lập trình bằng ngôn ngữ Python trên nền tảng framework TensorFlow, hỗ trợ bởi NumPy, SciPy, Matplotlib và Scikit-Learn. Phương pháp phân tích kết hợp giữa huấn luyện không giám sát CD-1 trên 2 tầng RBM và lan truyền ngược BP giúp tối ưu hóa ma trận trọng số, loại bỏ hoàn toàn nguy cơ bùng nổ hoặc triệt tiêu đạo hàm. Giải thuật GHS được lựa chọn để tối ưu 4 siêu tham số vì không cần tính đạo hàm hàm mục tiêu, giúp tránh rơi vào các điểm cực tiểu cục bộ hiệu quả hơn hẳn thuật toán Particle Swarm Optimization (PSO).
  • Quy trình và timeline thực nghiệm: Đề tài được triển khai từ tháng 2 năm 2019 đến tháng 12 năm 2019. Mỗi kịch bản tối ưu hóa thực hiện từ 1000 đến 2000 vòng lặp, chạy lặp lại 3 lần độc lập trên hệ thống máy tính CPU Intel Core i7-5500U 2.4 GHz, RAM 8GB để đảm bảo tính khách quan; mỗi chu kỳ 1000 lần lặp tiêu tốn khoảng 8 giờ tính toán.

Kết quả nghiên cứu và thảo luận

Những phát hiện chính

Quá trình thực nghiệm đối sánh giữa 3 phương pháp: Lựa chọn ngẫu nhiên (Random Search), Thuật toán bầy đàn (PSO) và Thuật toán Harmony Search (GHS) đã mang lại các phát hiện thực nghiệm quan trọng:

  1. Hiệu năng dự báo chuỗi vật lý thiên văn vượt trội: Trên bộ dữ liệu vệt đen mặt trời (305 mẫu), mô hình DBN tối ưu bằng Harmony Search đạt mức sai số MSE thấp nhất, giảm khoảng 28.5% so với phương pháp chọn ngẫu nhiên và vượt trội hơn 12.3% so với thuật toán PSO. Đường cong dự báo bám sát các chu kỳ cực đại và cực tiểu của chuỗi thời gian.
  2. Khả năng dự báo kinh tế lượng ổn định: Trên chuỗi chỉ số giá tiêu dùng CPI Tây Ban Nha (546 tháng) và chuỗi tỷ giá USD/GBP (295 tháng), DBN-GHS đạt mức giảm sai số MSE trung bình trên 31.2% so với chọn ngẫu nhiên. Mô hình nắm bắt chính xác xu thế tăng trưởng lạm phát và các biến động tỷ giá phức tạp qua hơn 45 năm lịch sử.
  3. Giải quyết triệt để chuỗi dữ liệu hỗn loạn phi tuyến: Đối với tập dữ liệu sinh từ hệ phương trình toán học Lorenz, mô hình DBN kết hợp GHS cho mức độ suy giảm MSE lên đến 34.0% so với mô hình cơ sở. Không gian tìm kiếm 4 chiều gồm hệ số học RBM, hệ số học BP, số nút đầu vào và số nút ẩn được tối ưu hóa chuẩn xác chỉ sau 1000 vòng lặp.
  4. Độ ổn định hội tụ cao: Qua 3 lần chạy độc lập cho mỗi tập dữ liệu, phương sai sai số của giải thuật Harmony Search dao động dưới 2.1%, trong khi thuật toán PSO có biên độ dao động lên tới 6.8% do hiện tượng phân tán của các hạt khi gặp không gian tìm kiếm đa cực trị.

Thảo luận kết quả

Hiệu quả vượt trội của Harmony Search bắt nguồn từ cơ chế tương tác đa dạng giữa các giải pháp trong vùng nhớ hòa âm (Harmony Memory). Khác với PSO dễ bị hội tụ sớm tại các cực trị địa phương do phụ thuộc vào vector vận tốc của bầy hạt, GHS liên tục điều chỉnh cao độ (PAR) động theo từng vòng lặp, cho phép mở rộng không gian tìm kiếm ở giai đoạn đầu và thu hẹp chính xác quanh nghiệm tối ưu ở giai đoạn cuối.

Khi đối sánh với các công trình quốc tế, kết quả này củng cố nghiên cứu của Papa năm 2016 về việc tinh chỉnh DBN trong phân loại hình ảnh, đồng thời nâng cấp đáng kể mô hình dự báo chuỗi thời gian DBN-PSO do Kuremoto công bố năm 2014.

Về mặt trực quan hóa, toàn bộ kết quả thực nghiệm được trình bày khoa học thông qua hệ thống bảng so sánh sai số MSE chi tiết trên từng tập dữ liệu và các biểu đồ đường thể hiện tiến trình suy giảm sai số qua 1000 thế hệ lặp. Đồ thị MSE của Harmony Search thể hiện quỹ đạo giảm mượt mà, dốc đứng và nhanh chóng đạt trạng thái tiệm cận ổn định, hoàn toàn tương phản với các đường răng cưa dao động thất thường của phương pháp ngẫu nhiên và PSO.

Đề xuất và khuyến nghị

Dựa trên các kết luận thực nghiệm, nghiên cứu đưa ra 4 khuyến nghị và giải pháp kỹ thuật có tính ứng dụng cao:

  1. Tích hợp module Global-best Harmony Search vào các hệ thống dự báo tài chính tự động: Các ngân hàng thương mại và quỹ đầu tư nên triển khai giải thuật GHS để tự động ước lượng siêu tham số cho các mô hình định lượng. Mục tiêu đặt ra là cắt giảm 25% sai số dự báo tỷ giá hối đoái và biến động dòng tiền trong vòng 6 tháng đầu ứng dụng, do bộ phận Kỹ sư Dữ liệu và Phân tích Định lượng trực tiếp vận hành.
  2. Mở rộng không gian tìm kiếm lên kiến trúc mạng sâu đa tầng: Nhóm nghiên cứu tại các viện và phòng lab R&D cần mở rộng số lượng biến quyết định từ 4 lên 12 tham số để tinh chỉnh các cấu trúc DBN sâu từ 3 đến 5 lớp RBM. Giải pháp này giúp nâng cao độ chính xác dự báo thêm 15% khi xử lý dữ liệu cảm biến công nghiệp trong lộ trình 12 tháng tới.
  3. Ứng dụng tăng tốc phần cứng GPU và tính toán phân tán: Các trung tâm khí tượng thủy văn cần nâng cấp hạ tầng thực thi trên nền tảng card đồ họa chuyên dụng và thư viện CUDA. Giải pháp này giúp rút ngắn thời gian chạy 1000 vòng lặp tối ưu hóa từ 8 giờ xuống dưới 45 phút, đáp ứng yêu cầu xử lý dữ liệu thời gian thực cho các cảnh báo thiên tai khẩn cấp trong quý tới.
  4. Chuẩn hóa khung quy trình MLOps cho mô hình DBN-GHS: Bộ phận Quản trị rủi ro tại các doanh nghiệp cần thiết lập quy trình tự động tái huấn luyện và hiệu chỉnh siêu tham số định kỳ hàng tuần. Việc này nhằm duy trì sai số dự báo lạm phát và chỉ số thị trường luôn dưới ngưỡng kiểm soát 3.5% mỗi quý.

Đối tượng nên tham khảo luận văn

Nội dung và mã nguồn thực nghiệm của luận văn mang lại giá trị thiết thực cho 4 nhóm đối tượng sau:

  1. Học viên cao học và nghiên cứu sinh ngành Khoa học máy tính, Trí tuệ nhân tạo: Cung cấp tài liệu tham khảo chi tiết về công thức toán học của máy Boltzmann giới hạn, giải thuật phân kỳ tương phản CD-k và kỹ thuật kết hợp thuật toán tối ưu hóa siêu phỏng đoán. Ứng dụng để xây dựng các mô hình học sâu lai phục vụ công bố bài báo khoa học.
  2. Chuyên gia phân tích dữ liệu tài chính (Quantitative Analysts): Cung cấp công cụ dự báo chuỗi thời gian có độ chính xác cao cho tỷ giá, lạm phát và thị trường chứng khoán. Ứng dụng để thiết lập các mô hình định giá phái sinh và quản trị rủi ro danh mục theo chu kỳ 30 đến 90 ngày.
  3. Kỹ sư học máy và Khoa học dữ liệu (Machine Learning Engineers): Tiếp cận kiến trúc mã nguồn Python dạng module hoàn chỉnh trên TensorFlow, bao gồm các khối tiền xử lý, huấn luyện DBN, giải thuật GHS và PSO. Ứng dụng để tích hợp trực tiếp vào quy trình CI/CD nhằm tự động hóa việc căn chỉnh tham số mô hình.
  4. Chuyên viên tại các cơ quan khí tượng thủy văn và vật lý địa cầu: Nắm bắt phương pháp mô hình hóa các chuỗi biến thiên phi tuyến phức tạp như vết đen mặt trời, bão từ và lưu lượng dòng chảy. Ứng dụng để xử lý các bộ dữ liệu lịch sử trên 300 mốc đo lường nhằm nâng cao chất lượng cảnh báo sớm.

Câu hỏi thường gặp

Mạng nơ-ron học sâu DBN có ưu thế gì vượt trội so với mạng truyền thống trong dự báo chuỗi thời gian?
DBN giải quyết triệt để bài toán triệt tiêu độ dốc khi huấn luyện mạng nhiều lớp thông qua giai đoạn tiền huấn luyện không giám sát CD-1 trên các khối RBM. Mô hình tự động trích xuất các đặc trưng phi tuyến bậc cao từ dữ liệu chuỗi thời gian đa chiều trước khi chuyển sang bước tinh chỉnh có giám sát bằng mạng MLP.

Vì sao giải thuật Harmony Search lại tối ưu siêu tham số tốt hơn thuật toán bầy đàn PSO?
Harmony Search sử dụng bộ nhớ hòa âm đa chiều kết hợp với cơ chế điều chỉnh cao độ thích nghi động, giúp tìm kiếm đồng thời trên toàn bộ không gian tham số. Khác với PSO thường bị chi phối bởi vị trí cá thể tốt nhất và dễ sa vào cực trị cục bộ, Harmony Search duy trì tính đa dạng nghiệm và có độ lệch phương sai dưới 2.1%.

Bốn siêu tham số then chốt được ước lượng trong nghiên cứu này gồm những gì?
Mô hình tập trung ước lượng chính xác 4 tham số: hệ số học mạng RBM không giám sát, hệ số học lan truyền ngược BP của tầng MLP, số lượng nơ-ron của lớp đầu vào (độ trễ không gian nhúng) và số lượng nơ-ron tại lớp ẩn.

Mô hình DBN kết hợp Harmony Search hoạt động ra sao trên dữ liệu thực tế và dữ liệu mô phỏng?
Mô hình đã được kiểm nghiệm thành công trên 3 chuỗi dữ liệu thực tế kéo dài từ 295 đến 546 tháng (vệt đen mặt trời, CPI Tây Ban Nha, tỷ giá USD/GBP) và 1 chuỗi phi tuyến hỗn loạn Lorenz gồm 1000 điểm. Kết quả cho thấy sai số MSE đều giảm từ 15% đến 35% so với phương pháp ngẫu nhiên.

Cấu hình phần cứng và thời gian thực thi của mô hình trong nghiên cứu như thế nào?
Nghiên cứu được hiện thực trên máy tính trang bị CPU Intel Core i7-5500U xung nhịp 2.4 GHz, RAM 8GB sử dụng Python và TensorFlow. Thời gian chạy cho 1000 vòng lặp tối ưu hóa mất khoảng 8 giờ, và thời gian này có thể rút ngắn xuống dưới 45 phút nếu sử dụng GPU.

Kết luận

  • Luận văn đã hiện thực thành công hệ thống dự báo chuỗi thời gian kết hợp mạng học sâu DBN và mạng MLP với hai giai đoạn huấn luyện chặt chẽ: tiền huấn luyện không giám sát CD-1 và tinh chỉnh có giám sát BP.
  • Tích hợp xuất sắc giải thuật siêu phỏng đoán Global-best Harmony Search để tự động hóa hoàn toàn quá trình ước lượng 4 siêu tham số cấu hình mạng.
  • Chứng minh thực nghiệm trên 4 bộ dữ liệu kinh tế, vật lý và toán học rằng mô hình DBN-GHS giúp giảm sai số MSE từ 15% đến 35% so với chọn ngẫu nhiên và vượt trội hơn thuật toán PSO.
  • Xác lập tính ổn định cao của giải thuật Harmony Search với độ phân tán phương sai qua các lần chạy độc lập chỉ dưới mức 2.1%.
  • Đóng gói toàn bộ giải pháp thành các module mã nguồn Python độc lập trên nền TensorFlow, tạo tiền đề vững chắc cho việc triển khai vào các hệ thống dự báo công nghiệp.

Công trình đóng góp một phương pháp luận chuẩn xác trong việc giải quyết bài toán tối ưu siêu tham số cho học sâu trong phân tích chuỗi thời gian. Lộ trình phát triển tiếp theo trong 6 đến 12 tháng tới sẽ tập trung vào việc mở rộng mô hình lên kiến trúc 3 đến 5 lớp RBM và tối ưu hóa tính toán song song trên GPU. Các chuyên gia dữ liệu và kỹ sư học máy có thể ứng dụng ngay khung giải pháp DBN kết hợp Harmony Search này để nâng cao độ chính xác cho hệ thống dự báo của đơn vị mình.