Tổng quan nghiên cứu

Theo báo cáo thống kê từ Tổ chức Y tế Thế giới (WHO), các bệnh lý tim mạch (CVD) là nguyên nhân gây tử vong hàng đầu trên toàn cầu, cướp đi sinh mạng của khoảng 17,7 triệu người mỗi năm, chiếm tỷ lệ 31% trong tổng số các ca tử vong. Tại Việt Nam và nhiều quốc gia đang phát triển, số lượng bệnh nhân mắc các chứng rối loạn nhịp tim nguy hiểm đang gia tăng nhanh chóng do biến đổi lối sống, đái tháo đường và các bệnh lý chuyển hóa. Điện tâm đồ (ECG) là công cụ chẩn đoán không xâm lấn tiêu chuẩn để theo dõi hoạt động điện học của tim. Tuy nhiên, việc giải mã thủ công tín hiệu ECG trong các bản ghi kéo dài 24 đến 48 giờ đối mặt với thách thức lớn: biên độ điện áp rất nhỏ (thường dưới 10 mV) và thời gian xuất hiện các biến đổi bất thường chỉ dao động từ 80 đến 420 mili-giây.

Vấn đề cốt lõi của nghiên cứu là xây dựng giải pháp phân loại tự động tín hiệu điện tim đạt độ chính xác cao, khắc phục hiện tượng mất cân bằng dữ liệu giữa các nhóm bệnh và sự sai lệch tín hiệu giữa các đạo trình đo khác nhau. Đề tài tập trung vào mục tiêu ứng dụng phương pháp học chuyển đổi khoảng cách cục bộ kết hợp mạng nơ-ron nhân tạo và biến đổi sóng rời rạc để nhận dạng toàn diện 5 nhóm nhịp tim theo tiêu chuẩn quốc tế ANSI/AAMI EC57:1998.

Nghiên cứu được triển khai thực nghiệm trên cơ sở dữ liệu chuẩn quốc tế MIT-BIH Arrhythmia, bao gồm 48 bản ghi tín hiệu điện tim số hóa từ 47 bệnh nhân với xấp xỉ 110.000 nhịp tim được gắn nhãn y khoa. Kết quả của nghiên cứu đóng góp một khung thuật toán hoàn chỉnh, giúp nâng cao độ chính xác phân loại toàn diện vượt mức 98%, giảm thiểu tỷ lệ dương tính sai và tối ưu hóa thời gian tính toán nhằm phục vụ các hệ thống hỗ trợ chẩn đoán lâm sàng thời gian thực.

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

Khung lý thuyết áp dụng

Nghiên cứu được xây dựng dựa trên ba trụ cột lý thuyết nền tảng trong xử lý tín hiệu y sinh và học máy có giám sát:

Thứ nhất là lý thuyết phân tích thời gian - tần số qua biến đổi Wavelet rời rạc (DWT). Khác với biến đổi Fourier truyền thống bị giới hạn độ phân giải do chiều dài cửa sổ cố định, DWT sử dụng các hàm sóng trực giao được co giãn và tịnh tiến để phân rã tín hiệu điện tim phi tĩnh thành các dải tần số đặc trưng mà vẫn bảo toàn vị trí thời gian của các sóng P, Q, R, S, T.

Thứ hai là lý thuyết học khoảng cách Large Margin Nearest Neighbor (LMNN) trên không gian metric Mahalanobis. Thuật toán tối ưu hóa một phép biến đổi tuyến tính nhằm tạo ra biên an toàn rộng nhất: kéo các điểm dữ liệu cùng nhãn lại gần nhau trong phạm vi k lân cận mục tiêu (thiết lập k = 3), đồng thời đẩy các điểm dữ liệu khác nhãn xâm lấn ra xa khỏi ranh giới phân tách. Hàm chi phí dạng Hinge Loss được giải quyết triệt để thông qua bài toán quy hoạch bán xác định dương (SDP) với hệ số cân bằng tham số thực nghiệm được xác định ở mức 0,5.

Thứ ba là mô hình Mạng nơ-ron nhân tạo truyền thẳng nhiều lớp huấn luyện bằng thuật toán lan truyền ngược (Back-Propagation Neural Network). Mạng sử dụng các đơn vị xử lý phi tuyến với hàm kích hoạt Sigmoid và hàm mục tiêu tổng bình phương sai số (MSE) để tìm kiếm các trọng số tối ưu hóa khả năng phân tách đa lớp.

Các khái niệm chuyên ngành cốt lõi được lượng hóa bao gồm: phức hợp sóng QRS với độ rộng tiêu chuẩn 80-120 mili-giây, khoảng thời gian PR từ 120-200 mili-giây, khoảng cách Mahalanobis, ma trận hiệp phương sai và hàm chi phí lồi.

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

Nguồn dữ liệu thực nghiệm được trích xuất từ cơ sở dữ liệu MIT-BIH Arrhythmia Database do Trung tâm Y tế Beth Israel Deaconess và Viện Công nghệ Massachusetts xây dựng. Dữ liệu gồm 48 bản ghi điện tim 2 kênh, mỗi bản ghi có độ dài 30 phút, được số hóa ở tần số lấy mẫu 360 Hz (360 mẫu trên mỗi giây) với bộ chuyển đổi ADC 11-bit trên dải điện áp 10 mV. Toàn bộ cỡ mẫu gồm 224.908 nhịp tim được chia thành 5 nhóm bệnh theo khuyến nghị của Hiệp hội Phát triển Thiết bị Y tế (AAMI): nhóm N (nhịp bình thường, 181.022 mẫu), nhóm S (ngoại tâm thu trên thất, 5.552 mẫu), nhóm V (ngoại tâm thu thất, 14.464 mẫu), nhóm F (nhịp hỗn hợp, 1.604 mẫu) và nhóm Q (nhịp không xác định, 22.266 mẫu).

Phương pháp chọn mẫu và tiền xử lý tín hiệu được tiến hành bằng cách phân đoạn chính xác 200 mẫu tín hiệu cho mỗi chu kỳ tim (gồm 99 mẫu trước đỉnh R, 100 mẫu sau đỉnh R và điểm cực trị R). Tín hiệu sau đó được phân tích qua bộ lọc DWT cấp độ 4 với họ wavelet Daubechies (Db4) để trích xuất dải xấp xỉ mức 4 (tần số 0 - 11,25 Hz) và dải chi tiết mức 4 (tần số 11,25 - 22,5 Hz). Kỹ thuật giảm chiều dữ liệu phân tích thành phần chính (PCA) và phân tích phân biệt tuyến tính (LDA) được áp dụng độc lập nhằm rút gọn không gian đặc trưng, loại bỏ trên 75% dữ liệu dư thừa trước khi đưa vào mô hình học máy.

Toàn bộ quy trình tính toán, huấn luyện và đánh giá chéo được thực thi trên môi trường MATLAB R2017a kết hợp công cụ chuyên dụng WFDB Toolbox và hệ thống cơ sở dữ liệu MySQL theo tiến trình ba giai đoạn thí nghiệm độc lập kéo dài trong 12 tháng.

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 đã ghi nhận bốn phát hiện khoa học quan trọng có đối sánh số liệu định lượng cụ thể:

Thứ nhất, nghiên cứu chứng minh tác động tiêu cực rõ rệt của hiệu ứng kéo theo (Long-tail effect) do sự mất cân bằng dữ liệu tự nhiên. Trong tập dữ liệu MIT-BIH gốc, nhóm nhịp bình thường N chiếm tới 80,48% tổng số mẫu, trong khi nhóm bệnh F chỉ chiếm 0,71%. Khi sử dụng tập dữ liệu chưa cân bằng (tập DS11), mạng nơ-ron đạt độ chính xác chung cao nhưng độ nhạy đối với nhóm bệnh F và S chỉ đạt khoảng 61,5%. Sau khi áp dụng chiến lược tái cấu trúc và cân bằng mẫu (tập DS13), độ nhạy phân loại các nhóm bệnh hiếm tăng vọt lên 93,8%, khẳng định tầm quan trọng của việc chuẩn hóa tỷ lệ phân bố mẫu.

Thứ hai, phương pháp giảm chiều LDA thể hiện sự vượt trội so với PCA khi kết hợp trực tiếp với mạng nơ-ron. Cụ thể, mô hình DWT + LDA + Neural Network đạt độ chính xác trung bình 96,8%, cao hơn 4,2% so với mô hình DWT + PCA + Neural Network (chỉ đạt 92,6%). Nguyên nhân là do LDA tối ưu hóa trực tiếp tỷ số phương sai giữa các lớp và phương sai nội lớp, tạo ra các ranh giới phân tách rõ nét hơn giữa các nhóm nhịp tim.

Thứ ba, việc tích hợp phương pháp học chuyển đổi metric LMNN mang lại bước nhảy vọt về hiệu năng phân loại. Khi áp dụng LMNN sau bước giảm chiều PCA, độ chính xác của mạng nơ-ron tăng từ 92,6% lên 98,4% (mức tăng trưởng 5,8%). Đồng thời, diện tích dưới đường cong ROC (chỉ số AUC) tăng từ 0,935 lên xấp xỉ 0,991, chứng minh thuật toán đã kéo dãn thành công khoảng cách biên an toàn giữa các cụm đặc trưng phức tạp.

Thứ tư, thử nghiệm đánh giá chéo đa đạo trình chỉ ra sự suy giảm hiệu suất nghiêm trọng khi kiểm tra chéo giữa các dây điện cực. Mô hình được huấn luyện trên đạo trình MLII (tập DS13) khi kiểm tra trên đạo trình V1 (tập DS23) bị sụt giảm độ chính xác khoảng 21,3%. Tuy nhiên, khi kết hợp dữ liệu từ cả hai dây điện cực (tập DS32), độ chính xác tổng thể được khôi phục và duy trì ổn định ở mức 97,6%.

Thảo luận kết quả

Kết quả phân loại được biểu diễn trực quan thông qua ma trận nhầm lẫn (Confusion Matrix) kích thước 5x5 và hệ thống đồ thị đường cong ROC đa lớp. Ma trận nhầm lẫn ghi nhận tỷ lệ phân loại đúng (True Positive) vượt trội ở cả 5 nhóm: nhóm N đạt 99,1%, nhóm S đạt 92,4%, nhóm V đạt 96,7%, nhóm F đạt 91,8% và nhóm Q đạt 98,5%. Tỷ lệ dương tính giả (FP) và âm tính giả (FN) đối với các ca rối loạn nhịp thất nguy hiểm giảm xuống dưới ngưỡng 3,5%.

Đồ thị đường cong ROC của mô hình đề xuất có quỹ đạo áp sát góc trên cùng bên trái của hệ trục tọa độ, minh chứng cho việc mô hình duy trì độ nhạy cực cao ngay cả khi tỷ lệ báo động giả tiệm cận 0%. Cơ chế học chuyển đổi LMNN đóng vai trò quyết định trong việc biến đổi không gian khoảng cách Euclidean thông thường thành không gian Mahalanobis tối ưu, giải quyết triệt để sự chồng lấn giữa các dạng sóng ngoại tâm thu trên thất (S) và nhịp bình thường (N).

Khi so sánh với các công bố quốc tế cùng thời kỳ sử dụng máy học véc-tơ hỗ trợ (SVM) hoặc hệ mờ TSK loại 2 (thường chỉ đạt độ chính xác từ 92,0% đến 93,3% trên một vài nhóm loạn nhịp đơn lẻ), mô hình kết hợp DWT, LMNN và Back-Propagation Neural Network trong luận văn thể hiện sự vượt trội về cả độ chính xác toàn diện, độ ổn định dữ liệu và khả năng phân loại đồng thời 5 nhóm bệnh lý phức tạp.

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

Từ các phát hiện thực nghiệm, nghiên cứu đưa ra bốn nhóm giải pháp và khuyến nghị khả thi:

  1. Chuẩn hóa quy trình thu thập dữ liệu điện tim đa đạo trình: Các cơ sở y tế và nhóm nghiên cứu kỹ thuật y sinh cần áp dụng bắt buộc phương thức ghi nhận tín hiệu đồng thời trên tối thiểu hai đạo trình (như MLII kết hợp V1 hoặc V5). Mục tiêu là giảm thiểu sai số chẩn đoán chéo xuống dưới 2,5%, thực hiện áp dụng trong toàn bộ các nghiên cứu lâm sàng trong vòng 6 tháng tới.

  2. Tích hợp thuật toán LMNN và mạng nơ-ron vào thiết bị y tế nhúng: Các doanh nghiệp sản xuất thiết bị y tế cần tiến hành đóng gói mô hình thuật toán DWT-LMNN-BPNN vào các dòng vi xử lý số tín hiệu (DSP) hoặc vi mạch FPGA tiêu thụ năng lượng thấp. Hướng tới mục tiêu xử lý phân loại nhịp tim theo thời gian thực với độ trễ dưới 50 mili-giây trên các thiết bị theo dõi Holter ECG đeo tay trong lộ trình 12 tháng.

  3. Xây dựng ngân hàng dữ liệu điện tim số hóa quốc gia: Bộ Y tế cùng các bệnh viện chuyên khoa đầu ngành cần phối hợp triển khai đề án số hóa và gắn nhãn dữ liệu ECG theo tiêu chuẩn WFDB/MIT format. Đặt chỉ tiêu thu thập và chuẩn hóa hơn 500.000 bản ghi điện tim đa dạng bệnh lý vào giai đoạn 2026-2028 nhằm cung cấp nguồn dữ liệu huấn luyện máy học đạt độ tin cậy cao cho các nhà khoa học trong nước.

  4. Nâng cấp cấu trúc mô hình phân loại tự động mở rộng: Các nhóm nghiên cứu học thuật cần tiếp tục tối ưu hóa thuật toán bằng cách kết hợp cơ chế học chuyển đổi LMNN với các kiến trúc mạng học sâu tiên tiến (như CNN kết hợp Bi-LSTM hoặc Transformer). Mục tiêu nâng chỉ số F1-score của các nhóm bệnh hiếm (như nhóm F) vượt mức 96% trong khoảng thời gian nghiên cứu 18 tháng tới.

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

Nội dung và kết quả của luận văn mang lại giá trị thực tiễn và học thuật sâu sắc cho bốn nhóm đối tượng:

  1. Bác sĩ chuyên khoa tim mạch và nhân viên y tế lâm sàng: Luận văn cung cấp hiểu biết chuyên sâu về nguyên lý phân loại tự động 5 nhóm rối loạn nhịp tim theo tiêu chuẩn ANSI/AAMI. Bác sĩ có thể ứng dụng kết quả này làm cơ sở đánh giá độ tin cậy của các phần mềm đọc điện tim tự động, giúp rút ngắn 70% thời gian phân tích các bản ghi Holter ECG 24 giờ.

  2. Kỹ sư kỹ thuật y sinh (Biomedical Engineers): Tài liệu cung cấp trọn vẹn quy trình tiền xử lý tín hiệu, từ kỹ thuật lọc nhiễu, biến đổi Wavelet rời rạc 4 cấp độ đến giải thuật giảm chiều dữ liệu PCA/LDA, hỗ trợ thiết kế các mô-đun phần mềm chẩn đoán y tế có độ chính xác cao.

  3. Nhà nghiên cứu và học viên cao học chuyên ngành Trí tuệ nhân tạo và Xử lý tín hiệu: Luận văn là tài liệu tham khảo giá trị về phương pháp học chuyển đổi metric LMNN, kỹ thuật tối ưu hóa quy hoạch bán xác định (SDP) và giải pháp khắc phục hiện tượng mất cân bằng dữ liệu chuỗi thời gian sinh học.

  4. Doanh nghiệp phát triển thiết bị chăm sóc sức khỏe thông minh (IoMT): Các nhà sản xuất thiết bị theo dõi sức khỏe cá nhân có thể khai thác thuật toán để tích hợp tính năng cảnh báo sớm rối loạn nhịp tim vào đồng hồ thông minh và máy đo ECG gia đình kết nối đám mây.

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

  1. Phương pháp biến đổi Wavelet rời rạc (DWT) có ưu điểm gì vượt trội so với biến đổi Fourier trong phân tích tín hiệu ECG? Biến đổi Fourier chỉ phân tích tín hiệu trong miền tần số và đánh mất hoàn toàn mốc thời gian xuất hiện sóng, trong khi tín hiệu ECG có tính chất phi tĩnh và biến thiên liên tục. DWT giải quyết triệt để nhược điểm này bằng cách sử dụng các hàm sóng cục bộ, cho phép trích xuất chính xác năng lượng của dải tần số 0-22,5 Hz tại đúng thời điểm xuất hiện các phức hợp sóng P, QRS và T.

  2. Kỹ thuật học chuyển đổi khoảng cách LMNN giúp cải thiện độ chính xác phân loại như thế nào? LMNN tối ưu hóa một ma trận khoảng cách Mahalanobis thay vì sử dụng khoảng cách Euclidean cố định. Thuật toán thiết lập biên an toàn kéo các mẫu cùng loại bệnh lại gần 3 lân cận mục tiêu và đẩy các mẫu khác loại ra ngoài biên. Nhờ đó, ranh giới phân lớp của mạng nơ-ron trở nên rõ ràng, giúp nâng độ chính xác từ 92,6% lên 98,4%.

  3. Hiệu ứng mất cân bằng dữ liệu (Long-tail effect) ảnh hưởng thế nào đến bộ phân loại điện tim? Trong tập dữ liệu tự nhiên, số lượng nhịp tim bình thường chiếm trên 80%, trong khi các dạng loạn nhịp nguy hiểm như nhịp hỗn hợp (F) chỉ chiếm dưới 1%. Nếu không cân bằng dữ liệu, mô hình sẽ bị thiên lệch về nhóm đa số, khiến độ nhạy phát hiện các nhóm bệnh hiếm bị sụt giảm nghiêm trọng xuống mức khoảng 61,5%.

  4. Tại sao cần phải đánh giá chéo giữa các dây điện cực (đạo trình) khác nhau? Hình thái sóng điện tim thay đổi đáng kể tùy thuộc vào vị trí đặt điện cực (như MLII, V1 hay V5). Thực nghiệm cho thấy mô hình chỉ huấn luyện trên một đạo trình duy nhất sẽ giảm hơn 20% độ chính xác khi kiểm tra trên đạo trình khác. Do đó, việc đánh giá chéo giúp đảm bảo mô hình có khả năng thích ứng linh hoạt trong môi trường lâm sàng thực tế.

  5. Tiêu chuẩn quốc tế ANSI/AAMI EC57 phân loại nhịp tim thành những nhóm nào? Tiêu chuẩn ANSI/AAMI phân chia toàn bộ các dạng nhịp tim thành 5 nhóm chính: nhóm N (nhịp bình thường và block nhánh), nhóm S (ngoại tâm thu trên thất), nhóm V (ngoại tâm thu thất), nhóm F (nhịp hỗn hợp giữa thất và bình thường), và nhóm Q (nhịp không thể phân loại hoặc nhịp máy tạo nhịp). Phân loại này giúp chuẩn hóa việc so sánh hiệu năng giữa các hệ thống chẩn đoán tự động trên toàn thế giới.

Kết luận

  • Luận văn đã xây dựng thành công quy trình phân loại tự động 5 nhóm tín hiệu điện tim ECG theo chuẩn ANSI/AAMI EC57:1998 với độ chính xác tổng thể đạt 98,4%.
  • Phương pháp biến đổi Wavelet rời rạc (Db4 cấp độ 4) kết hợp kỹ thuật giảm chiều LDA chứng minh hiệu quả lọc nhiễu và nén đặc trưng vượt trội so với các phương pháp truyền thống.
  • Ứng dụng đột phá của thuật toán học chuyển đổi LMNN giúp tối ưu hóa không gian khoảng cách Mahalanobis, nâng cao diện tích dưới đường cong ROC lên mức AUC = 0,991.
  • Nghiên cứu chỉ ra quy tắc then chốt trong việc xử lý hiệu ứng mất cân bằng dữ liệu và yêu cầu bắt buộc phải kết hợp dữ liệu đa đạo trình để loại bỏ sai số phân loại chéo.
  • Trong giai đoạn tiếp theo, nhóm tác giả định hướng tối ưu hóa thuật toán trên nền tảng phần cứng nhúng thời gian thực và mở rộng mô hình mạng nơ-ron sâu lai ghép trong giai đoạn 2026-2027.

Hãy kết nối và tham khảo toàn văn công trình nghiên cứu để tiếp cận đầy đủ các thuật toán xử lý tín hiệu y sinh và mô hình học máy tiên tiến phục vụ phát triển công nghệ chẩn đoán y khoa hiện đại.