Bối cảnh và vấn đề nghiên cứu

Sự phát triển của công nghệ thông tin trong các hoạt động kinh tế, xã hội và giáo dục đã dẫn đến sự hình thành của các kho dữ liệu có quy mô lớn. Các kỹ thuật quản trị và phân tích cơ sở dữ liệu truyền thống bộc lộ hạn chế trong việc trích xuất các tri thức tiềm ẩn phục vụ quá trình hỗ trợ ra quyết định. Bối cảnh này thúc đẩy sự phát triển của lĩnh vực Phát hiện tri thức và Khai phá dữ liệu (KDD – Knowledge Discovery and Data Mining), một hướng tiếp cận ra đời từ cuối thập niên 1980 (được định nghĩa bởi Fayyad, Piatetsky-Shapiro và Smyth vào năm 1989).

Khác với cơ sở dữ liệu tĩnh, cơ sở dữ liệu chuỗi thời gian (Time Series) chứa các tập bản ghi gắn liền với nhãn thời gian, trong đó yếu tố thời gian là một thuộc tính không thể tách rời trong quá trình phân tích. Xuất phát từ nhu cầu theo dõi sự biến thiên của các đại lượng theo thời gian và ứng dụng dự báo trong giáo dục, khóa luận tập trung nghiên cứu đề tài "Khai phá dữ liệu chuỗi thời gian".

Khóa luận xác định các mục tiêu nghiên cứu cụ thể như sau:

  1. Tìm hiểu về phương pháp phân tích chuỗi tuần tự theo thời gian và dự báo.
  2. Tìm hiểu về cơ sở dữ liệu chuỗi thời gian.
  3. Tìm hiểu về kỹ thuật khám phá tri thức trên cơ sở dữ liệu chuỗi thời gian.
  4. Áp dụng vào bài toán cụ thể: khám phá tri thức trên cơ sở dữ liệu điểm học sinh, sau đó thực hiện dự báo kết quả học tập của học sinh trong các giai đoạn tiếp theo.

Nhiệm vụ nghiên cứu tương ứng được phân định theo 4 nội dung:

  1. Nắm vững phương pháp phân tích chuỗi tuần tự theo thời gian và cách thức dự báo.
  2. Hệ thống hóa các khái niệm và đặc trưng của cơ sở dữ liệu chuỗi thời gian.
  3. Phân tích các phương pháp khám phá tri thức và kỹ thuật khai phá trên cơ sở dữ liệu chuỗi thời gian.
  4. Thực nghiệm mô hình dự báo kết quả học tập của học sinh dựa trên dữ liệu điểm số thực tế.

Đối tượng nghiên cứu của đề tài gồm hai thành phần chính: cơ sở dữ liệu chuỗi thời gian nói chung và việc khám phá tri thức trên cơ sở dữ liệu điểm học sinh nhằm dự báo kết quả học tập tiếp theo. Phạm vi nghiên cứu về mặt dữ liệu thực nghiệm được giới hạn trong bảng điểm tổng kết 4 học kỳ (lớp 10 và lớp 11) của 30 học sinh thuộc hai lớp 10A1 và 11A1 tại Trường THPT Nguyễn Thị Minh Khai (huyện Đức Thọ, tỉnh Hà Tĩnh), khóa học 2007 – 2010.

Cơ sở lý thuyết và phương pháp

Khóa luận xây dựng nội dung dựa trên các khung lý thuyết, mô hình toán học và phương pháp khai phá dữ liệu sau:

  • Khung khái niệm chuỗi thời gian: Chuỗi thời gian $X = (x_1, x_2, \dots, x_n)$ là tập hợp $n$ giá trị của một biến thực được ghi nhận tuần tự trên các khoảng thời gian bằng nhau. Các khái niệm liên quan bao gồm dãy con (subsequence) chiều dài $m$, hàm đại diện $x(t)$, hàm xấp xỉ tuyến tính $\beta(t) = \alpha t + \beta$ với hệ số góc $\alpha$ đóng vai trò là đặc trưng khuynh hướng, và khái niệm "sự kiện" (event) khi có sự biến đổi đặc trưng giữa hai đoạn liên tiếp.
  • Mô hình cấu trúc chuỗi thời gian: Phân rã chuỗi thời gian thành 4 thành phần: xu hướng dài hạn ($T_t$), mùa vụ ($S_t$), chu kỳ ($C_t$), và bất thường ($I_t$). Hai mô hình toán học cơ bản được sử dụng gồm mô hình nhân ($X_t = T_t \times S_t \times C_t \times I_t$) và mô hình cộng ($X_t = T_t + S_t + C_t + I_t$).
  • Độ đo tương đồng chuỗi thời gian: Đánh giá khoảng cách giữa các chuỗi thông qua ba độ đo:
    1. Khoảng cách Euclide: Tính toán khoảng cách hình học trực tiếp giữa hai vector cùng chiều dài $n$.
    2. Chuẩn hóa chuỗi thời gian (Zero-mean Normalization): Khử sai khác về biên độ bằng cách trừ giá trị trung bình $E_x$, sau đó tính khoảng cách Euclide trên chuỗi đã chuẩn hóa.
    3. Dãy con chung dài nhất (Longest Common Subsequence - LCS): Phương pháp do Vlachos và cộng sự đề xuất năm 2004, cho phép so khớp các chuỗi con tương đồng và bỏ qua các điểm dị biệt hoặc nhiễu.
  • Phương pháp gom cụm (Clustering): Sử dụng hàm chi phí gom cụm ($m_C$) và hai thuật toán phân hoạch chính:
    1. Thuật toán K-Means: Phân hoạch tập dữ liệu thành $k$ cụm dựa trên khoảng cách tới các điểm trọng tâm (centroid), với độ phức tạp thuật toán là $O(knl)$ (với $k$ là số cụm, $n$ là số phần tử, $l$ là số vòng lặp).
    2. Thuật toán Greedy: Gom cụm dựa trên số cụm $k$ và ngưỡng bán kính nhóm $d$ xác định trước.
  • Mô hình hồi quy theo thời gian (Linear Regression): Xây dựng phương trình xu thế tuyến tính $\hat{Y} = a + bt$. Khóa luận áp dụng kỹ thuật đổi biến thời gian đối xứng sao cho tổng chỉ số thời gian triệt tiêu ($\sum t = 0$), từ đó tính nhanh các hệ số hồi quy theo công thức: $$a = \frac{\sum y}{n}, \quad b = \frac{\sum yt}{\sum t^2}$$
  • Nguồn dữ liệu nghiên cứu: Dữ liệu thứ cấp thu thập từ sổ điểm học tập thực tế của 30 học sinh tại Trường THPT Nguyễn Thị Minh Khai, Hà Tĩnh.

Nội dung chính theo từng chương

Chương I: Tìm hiểu về phân tích chuỗi tuần tự theo thời gian và dự báo (Time series Analysis and Forecoasting)

Chương I trình bày tổng quan lý thuyết về chuỗi thời gian và các kỹ thuật dự báo thống kê cổ điển. Tác giả phân tích chi tiết bốn thành phần cấu thành chuỗi thời gian: xu hướng dài hạn (thể hiện bằng đường thẳng hoặc đường cong trơn), mùa vụ (biến thiên lặp lại trong vòng một năm), chu kỳ (biến động lặp lại với chu kỳ trên một năm) và bất thường (các biến cố ngẫu nhiên không thể dự đoán trước).

Các phương pháp làm trơn số liệu (Smoothing methods) được phân tích nhằm giảm thiểu tác động của yếu tố mùa vụ và bất thường:

  • Phương pháp trung bình dịch chuyển (Moving Average): Thay thế giá trị $X_t$ bằng trung bình của $(2m+1)$ điểm lân cận.
  • Phương pháp làm trơn bằng hàm số mũ đơn giản (Simple Exponential Smoothing): Do Brown đề xuất năm 1954 theo phương trình $F_{t+1} = \alpha X_t + (1-\alpha) F_t$, với hằng số làm trơn $0 \le \alpha \le 1$. Tác giả đưa ra kinh nghiệm chọn $\alpha = \frac{2}{n+1}$ và thực nghiệm tính toán trên bộ dữ liệu doanh số (CA) của công ty Z trong 22 quý (5,5 năm) để so sánh giữa hai tham số $\alpha = 0,3$ và $\alpha = 0,6$, sử dụng chỉ số sai số độ lệch tuyệt đối trung bình (MAD).
Tham số làm trơn Đặc điểm ứng xử với dữ liệu Đánh giá sai số qua chỉ số MAD
$\alpha = 0,3$ Trọng số phân bổ đều hơn về quá khứ, phù hợp chuỗi có tính ổn định cao Mức độ điều chỉnh theo biến động thực tế chậm hơn
$\alpha = 0,6$ Chú trọng các quan sát gần nhất, thích hợp chuỗi biến động không ổn định Bám sát các đỉnh biến động của chuỗi doanh số qua 22 quý

Chương này cũng phân loại dự báo theo thời gian (ngắn hạn dưới 3 năm, dài hạn từ 5 năm trở lên), theo dạng kết quả (dự báo điểm, dự báo khoảng) và xác lập quy trình dự báo 4 bước: Thu thập số liệu $\to$ Xử lý sơ bộ (chia tập dữ liệu thành nhóm đầu và nhóm kiểm tra) $\to$ Lựa chọn phương pháp và xây dựng mô hình $\to$ Thực hiện dự báo và phân tích sai số.

Hệ thống các phương pháp dự báo thông dụng được tác giả khảo sát gồm:

  1. Phương pháp giản đơn (Naive): Gồm 4 biến thể (gán giá trị gần nhất $F_{t+1} = X_t$, bổ sung thành phần xu hướng $F_{t+1} = X_t + (X_t - X_{t-1})$, tính theo tỷ lệ $F_{t+1} = X_t \times (X_t / X_{t-1})$, và tích hợp cả yếu tố mùa vụ lẫn xu hướng dựa trên dữ liệu cùng kỳ năm trước).
  2. Phương pháp trung bình: Trung bình giản đơn và trung bình dịch chuyển $n$ quan sát.
  3. Phương pháp tự hồi quy (Autoregressive - AR): Mô hình AR bậc 1, bậc 2 và bậc $p$.

Chương II: Cơ sở dữ liệu chuỗi thời gian

Chương II tập trung vào các kỹ thuật biểu diễn dữ liệu và khai phá tri thức trên cơ sở dữ liệu chuỗi thời gian. Quá trình KDD trên chuỗi thời gian được chia thành hai giai đoạn chính: Tiền xử lý dữ liệu và Khai phá dữ liệu.

Giai đoạn tiền xử lý thực hiện phân đoạn chuỗi thời gian theo các cửa sổ trượt có độ rộng cố định, thực hiện gom cụm các dãy con để chuyển đổi cơ sở dữ liệu chuỗi thời gian gốc thành cơ sở dữ liệu các mẫu đặc trưng (hoặc cơ sở dữ liệu các sự kiện). Giai đoạn khai phá tiếp theo sẽ trích xuất các luật kết hợp theo thời gian, tiêu biểu là dạng luật $A \xrightarrow{T} B$ ("Nếu mẫu A xuất hiện thì mẫu B sẽ xuất hiện sau khoảng thời gian $T$") hoặc dạng luật định vị thời điểm cụ thể ("Nếu A xuất hiện tại $T_1$ thì B xuất hiện tại $T_2$").

Nội dung trọng tâm của chương là phân tích các độ đo tương đồng và thuật toán gom cụm:

  • Phân tích ưu nhược điểm của khoảng cách Euclide, phương pháp chuẩn hóa chuỗi và độ đo LCS.
  • Trình bày chi tiết thuật toán gom cụm K-Means: khởi tạo $k$ tâm ngẫu nhiên, gán đối tượng vào cụm gần nhất, cập nhật lại tọa độ trọng tâm và lặp lại cho đến khi các cụm không thay đổi.
Thuật toán gom cụm Cơ chế phân hoạch Tham số đầu vào Ưu điểm và hạn chế nêu trong văn bản
K-Means Phân hoạch tập dữ liệu thành $k$ nhóm xoay quanh các điểm trọng tâm đại diện Số cụm $k$ (nguyên dương cố định) Đơn giản, dễ cài đặt; tuy nhiên phụ thuộc vào việc chọn $k$ và chi phí tính toán lớn khi dữ liệu tăng
Greedy Tìm kiếm tuần tự và gom các điểm vào nhóm nếu khoảng cách nhỏ hơn ngưỡng bán kính Số nhóm $k$ và bán kính nhóm $d$ Phù hợp khi đã biết trước bán kính cụm $d$; khoảng cách giữa hai trọng tâm tối thiểu là $d$

Chương II cung cấp bài toán thực nghiệm từng bước của K-Means trên tập dữ liệu mẫu gồm 4 đối tượng A(1,1), B(2,1), C(4,3), D(5,4) với $k=2$. Qua 3 bước lặp tính toán các ma trận khoảng cách $D_0, D_1, D_2$ và cập nhật tâm $c_1, c_2$, thuật toán hội tụ và phân hoạch thành 2 nhóm: Nhóm 1 gồm {A, B} và Nhóm 2 gồm {C, D}.

Chương III: Bài toán vận dụng – Sử dụng mô hình hồi quy theo thời gian để dự báo kết quả học tập của học sinh trong tương lai

Chương III vận dụng lý thuyết đã nghiên cứu vào việc giải quyết bài toán thực tế trong lĩnh vực giáo dục. Tác giả chỉ ra rằng công tác đánh giá học sinh truyền thống thường mang tính cảm tính, cào bằng qua điểm trung bình tĩnh mà chưa phản ánh được tiến trình phát triển năng lực của học sinh qua thời gian.

Để mô hình hóa, điểm số học tập của từng học sinh qua các học kỳ được xem như một chuỗi thời gian độc lập. Bộ dữ liệu gồm 30 học sinh với điểm tổng kết của 4 học kỳ liên tiếp: HK I, HK II (lớp 10) và HK III, HK IV (lớp 11). Với $n=4$ mốc thời gian, tác giả gán các chỉ số thời gian đối xứng là $t \in {-3, -1, 1, 3}$ để thỏa mãn điều kiện $\sum t = 0$, giúp đơn giản hóa việc tính toán các tham số $a$ và $b$ trong phương trình hồi quy tuyến tính $\hat{Y} = a + bt$.

Mã học sinh HK I ($t=-3$) HK II ($t=-1$) HK III ($t=1$) HK IV ($t=3$) Phương trình hồi quy xu thế $\hat{Y}$ Điểm xấp xỉ $\hat{Y}_1$ Điểm xấp xỉ $\hat{Y}_2$ Điểm xấp xỉ $\hat{Y}_3$ Điểm xấp xỉ $\hat{Y}_4$
A101 8,2 8,3 8,4 8,1 $\hat{Y} = 8,25 - 0,01t$ 8,28 8,26 8,24 8,22
A102 8,5 8,7 8,9 8,6 $\hat{Y} = 8,675 + 0,025t$ 8,60 8,65 8,70 8,75
A103 7,9 8,0 8,1 8,0 $\hat{Y} = 8,00 + 0,02t$ 7,94 7,98 8,02 8,06
A104 8,5 8,7 8,8 8,8 $\hat{Y} = 8,70 + 0,05t$ 8,55 8,65 8,75 8,85
A105 8,0 8,2 8,3 8,5 $\hat{Y} = 8,25 + 0,08t$ 8,01 8,17 8,33 8,49
A106 6,5 6,8 7,0 6,6 $\hat{Y} = 6,60 + 0,025t$ 6,525 6,575 6,625 6,675
A107 8,8 9,1 9,1 9,2 $\hat{Y} = 9,05 + 0,06t$ 8,87 8,99 9,11 9,23
A108 8,0 8,2 8,2 8,2 $\hat{Y} = 8,15 + 0,03t$ 8,06 8,12 8,18 8,24

Phương trình hồi quy cho phép xác định rõ xu hướng học tập của từng cá nhân: học sinh A101 có hệ số $b = -0,01$ thể hiện chiều hướng giảm sút nhẹ, trong khi học sinh A105 có hệ số $b = +0,08$ cho thấy mức độ tiến bộ nhanh nhất trong nhóm khảo sát.

Kết quả và đóng góp

Khóa luận đạt được các kết quả cụ thể sau:

  • Hệ thống hóa lý thuyết: Tổng hợp toàn diện các phương pháp xử lý chuỗi thời gian, từ các mô hình thống kê dự báo (làm trơn số mũ Brown, trung bình dịch chuyển, tự hồi quy AR) đến các kỹ thuật khai phá dữ liệu hiện đại (đo khoảng cách Euclide, LCS, chuẩn hóa chuỗi và thuật toán gom cụm K-Means, Greedy).
  • Thực nghiệm tính toán: Xây dựng thành công phương trình hồi quy xu thế cá thể hóa cho 30 học sinh thuộc Trường THPT Nguyễn Thị Minh Khai. Kết quả tính toán chỉ ra sự khác biệt về độ dốc biến thiên (hệ số $b$ dao động từ $-0,01$ đến $+0,08$) giữa các học sinh có cùng mức điểm trung bình tĩnh.
  • Đề xuất giải pháp: Đưa ra phương thức tiếp cận mới cho công tác quản lý giáo dục: chuyển đổi sổ điểm định kỳ thành cơ sở dữ liệu chuỗi thời gian, áp dụng mô hình toán học để ngoại suy kết quả học tập trong tương lai nhằm giảm thiểu tính cảm tính trong đánh giá học lực.
  • Đóng góp của đề tài: Khóa luận đã bước đầu kết nối lý thuyết khai phá chuỗi thời gian với bài toán đánh giá giáo dục phổ thông, chỉ ra tính khả thi của việc sử dụng các mô hình chuỗi thời gian thay thế cho các phương pháp thống kê tĩnh truyền thống.

Hạn chế và hướng nghiên cứu tiếp

Dựa trên phạm vi và nội dung văn bản, đề tài có các hạn chế và hướng mở sau:

  • Hạn chế:
    • Quy mô dữ liệu thực nghiệm còn hẹp, mới thực hiện trên 30 học sinh thuộc một trường THPT trong phạm vi 4 học kỳ.
    • Phần bài toán vận dụng ở Chương III mới triển khai mô hình hồi quy tuyến tính đơn biến theo thời gian, chưa kết hợp thử nghiệm thuật toán gom cụm K-Means hoặc trích xuất luật kết hợp $A \xrightarrow{T} B$ trên tập dữ liệu điểm số như khung lý thuyết đã xây dựng ở Chương II.
  • Hướng nghiên cứu tiếp:
    • Mở rộng tập dữ liệu điểm số trên quy mô nhiều khối lớp, nhiều trường học với chuỗi thời gian dài hơn.
    • Ứng dụng đầy đủ các kỹ thuật gom cụm chuỗi thời gian và khai phá luật liên kết thời gian để phân nhóm các mẫu học tập đặc trưng và dự báo đa chiều hơn.

Giá trị tham khảo

Khóa luận là tài liệu tham khảo phù hợp cho sinh viên, học viên ngành Sư phạm Tin học, Khoa học Máy tính, Hệ thống Thông tin và Công nghệ Thông tin đang nghiên cứu về khai phá dữ liệu chuỗi thời gian. Tài liệu cũng có giá trị tham khảo đối với cán bộ quản lý giáo dục muốn tìm hiểu phương pháp định lượng trong theo dõi tiến trình học tập của học sinh.

Các nội dung có giá trị tham khảo cao trong khóa luận gồm:

  • Phần tổng quan các công thức dự báo giản đơn, làm trơn số mũ và tính toán chỉ số MAD trong Chương I.
  • Quy trình từng bước thực thi thuật toán K-Means với ma trận khoảng cách cụ thể trong Chương II.
  • Phương pháp đặt biến thời gian triệt tiêu ($\sum t = 0$) để thiết lập phương trình hồi quy dự báo xu thế trong Chương III.

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

1. Kỹ thuật đặt biến thời gian đối xứng $\sum t = 0$ trong Chương III được thiết lập như thế nào với chuỗi 4 học kỳ?
Với $n=4$ học kỳ (số quan sát chẵn), tác giả gán thứ tự thời gian $t$ lần lượt là các giá trị đối xứng qua gốc tọa độ: $t_1 = -3$, $t_2 = -1$, $t_3 = 1$, $t_4 = 3$. Tổng các giá trị thời gian bằng 0 ($\sum t = -3 - 1 + 1 + 3 = 0$) và tổng bình phương $\sum t^2 = (-3)^2 + (-1)^2 + 1^2 + 3^2 = 20$, giúp loại bỏ thành phần liên đới trong công thức tính và xác định trực tiếp hệ số chặn $a = \frac{\sum y}{4}$ và hệ số góc $b = \frac{\sum yt}{20}$.

2. Thuật toán K-Means trong Chương II đạt điều kiện dừng sau bao nhiêu bước lặp trên ví dụ 4 loại thuốc?
Thuật toán K-Means trong ví dụ minh họa thực hiện qua 8 bước tính toán (3 lần phân nhóm). Tại bước 8, ma trận phân nhóm $G_2$ trùng khớp hoàn toàn với ma trận $G_1$ từ vòng lặp trước đó ($G_2 = G_1$), cho thấy không còn sự dịch chuyển nhóm của bất kỳ đối tượng nào, thuật toán dừng lại và phân chia thành 2 nhóm ổn định: Nhóm 1 {A, B} và Nhóm 2 {C, D}.

3. Khóa luận nêu ưu điểm nổi bật nào của độ đo LCS so với khoảng cách Euclide thông thường?
Độ đo dãy con chung dài nhất (LCS) do Vlachos và cộng sự đề xuất (2004) có ưu điểm trực quan và cho phép bỏ qua những điểm bất thường (nhiễu hoặc đột biến cục bộ) khi so sánh hai chuỗi thời gian. Trong khi đó, khoảng cách Euclide yêu cầu hai chuỗi phải có độ dài bằng nhau và rất nhạy cảm với sự dịch chuyển biên độ hoặc độ trễ thời gian.

4. Dữ liệu thực nghiệm của công ty Z trong Chương I được sử dụng nhằm mục đích gì?
Dữ liệu doanh số 22 quý (5,5 năm) của công ty Z được tác giả sử dụng để minh họa phương pháp làm trơn hàm mũ đơn giản của Brown (1954), so sánh độ thích ứng của chuỗi dự báo khi thay đổi hằng số làm trơn $\alpha$ giữa hai giá trị $\alpha = 0,3$ và $\alpha = 0,6$.

5. Kết quả phương trình hồi quy của học sinh A101 và A105 phản ánh xu hướng học tập khác nhau ra sao?
Học sinh A101 có phương trình $\hat{Y} = 8,25 - 0,01t$ với hệ số góc $b = -0,01 < 0$, biểu thị xu hướng điểm số suy giảm nhẹ qua các kỳ. Ngược lại, học sinh A105 có phương trình $\hat{Y} = 8,25 + 0,08t$ với hệ số góc $b = +0,08 > 0$, biểu thị xu thế tiến bộ rõ rệt qua 4 học kỳ.

Kết luận

Khóa luận "Khai phá dữ liệu chuỗi thời gian" của tác giả Trần Thị Tân đã hệ thống hóa cơ sở lý thuyết về phân tích, gom cụm và khám phá tri thức trên dữ liệu chuỗi thời gian. Đề tài đã vận dụng mô hình hồi quy xu thế tuyến tính để xử lý bảng điểm 4 học kỳ của 30 học sinh Trường THPT Nguyễn Thị Minh Khai (Hà Tĩnh), chứng minh tính khả thi của việc mô hình hóa kết quả học tập thành chuỗi thời gian. Công trình cung cấp tài liệu tham khảo có cấu trúc rõ ràng về mặt phương pháp luận thống kê và khai phá dữ liệu ứng dụng trong giáo dục.