Tổng quan nghiên cứu

Bệnh lý tim mạch hiện là nguyên nhân gây tử vong hàng đầu trên toàn cầu, chiếm khoảng 31% đến 33% tổng số ca tử vong mỗi năm theo thống kê y tế quốc tế. Tại Việt Nam, tần suất ghi nhận các ca rối loạn nhịp tim và nhồi máu cơ tim có xu hướng gia tăng liên tục qua từng năm. Điện tim đồ (ECG - Electrocardiogram) với 12 chuyển đạo tiêu chuẩn là công cụ chẩn đoán không xâm lấn thiết yếu trong lâm sàng. Tuy nhiên, tín hiệu ECG thực tế có đặc tính phi dừng, biên độ điện thế rất nhỏ (khoảng 0.12 mV đến 2 mV) và dải tần số thấp (0.05 Hz đến 100 Hz), khiến tín hiệu dễ bị suy giảm chất lượng bởi nhiều nguồn tạp âm phức tạp như nhiễu co cơ EMG, nhiễu mạng điện 50 Hz và nhiễu trôi đường đồng điện dưới 0.5 Hz.

Luận văn thạc sĩ kỹ thuật điện tử - viễn thông của tác giả Phạm Thị Quỳnh Trang, dưới sự hướng dẫn của Tiến sĩ Trịnh Anh Vũ tại Trường Đại học Công nghệ – Đại học Quốc gia Hà Nội, tập trung giải quyết bài toán cốt lõi: Ứng dụng cơ sở toán học của phép biến đổi Wavelet (Wavelet Transform) để xây dựng hệ thống thuật toán tự động tách lọc, nhận dạng các điểm đặc trưng của tín hiệu điện tim. Mục tiêu cụ thể của nghiên cứu là định vị chính xác phức bộ QRS, xác định điểm bắt đầu và kết thúc của sóng P, sóng T, cũng như tự động đo đạc khoảng cách QT trên các bản ghi ECG chuẩn.

Nghiên cứu được triển khai thử nghiệm trên môi trường mô phỏng Matlab với tín hiệu số hóa tại tần số lấy mẫu 250 Hz. Kết quả thực nghiệm chứng minh giải thuật Wavelet nâng cao độ nhạy phát hiện sóng tim đạt trên 99.2%, khắc phục hạn chế của các bộ lọc miền thời gian truyền thống (vốn chỉ đạt hiệu suất 96% đến 98%), đồng thời kiểm soát sai lệch trễ thời gian dưới 5 mẫu tín hiệu (tương đương 20 ms), mang lại giá trị thực tiễn cao cho việc phát triển thiết bị y tế thông minh.

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 thiết lập hệ thống cơ sở khoa học dựa trên hai trụ cột lý thuyết chính:

  1. Lý thuyết phân tích thời gian - tần số và biến đổi Wavelet: Phép biến đổi Fourier cổ điển (FT) chuyển tín hiệu từ miền thời gian sang miền tần số nhưng làm mất hoàn toàn mốc thời gian diễn ra sự kiện. Biến đổi Fourier cửa sổ trượt (STFT) do Dennis Gabor đề xuất năm 1946 đã cải thiện việc định vị thời gian nhưng bị ràng buộc bởi kích thước cửa sổ cố định theo nguyên lý bất định Heisenberg. Phép biến đổi Wavelet giải quyết triệt để vấn đề này bằng việc sử dụng các hàm cơ sở co giãn linh hoạt theo tỷ lệ (scale $a$) và dịch chuyển (translation $b$). Luận văn áp dụng cấu trúc đa phân giải MRA (Multiresolution Analysis) của Stephane Mallat trên lưới nhị phân Dyadic ($a_0 = 2, b_0 = 1$), cho phép phân tích tín hiệu ECG thành các thành phần xấp xỉ tần số thấp và chi tiết tần số cao.
  2. Lý thuyết điện sinh lý tim và hình thái học 12 chuyển đạo: Hoạt động co bóp của cơ tim tạo ra điện trường lan truyền khắp cơ thể. Tín hiệu thu thập gồm 3 chuyển đạo lưỡng cực chi (D1, D2, D3), 3 chuyển đạo đơn cực chi (VR, VL, VF) và 6 chuyển đạo trước tim (V1 đến V6). Một chu kỳ tim hoàn chỉnh gồm các thành phần: sóng P (khử cực nhĩ, độ dài 0.08 s đến 0.11 s, biên độ 0.12 mV đến 0.2 mV), khoảng PQ (dẫn truyền nhĩ - thất, dài 0.15 s đến 0.20 s), phức bộ QRS (khử cực thất, độ dốc lớn nhất), đoạn ST và sóng T (tái cực thất).

Khung nghiên cứu vận dụng 4 khái niệm trung tâm:

  • Hàm Wavelet Spline bậc hai: Là đạo hàm bậc nhất của hàm làm trơn, có tính chất đối xứng và pha tuyến tính, giúp duy trì ổn định dạng hình học của các đỉnh sóng.
  • Cực trị Modulus (Modulus Maxima): Cặp giá trị cực đại - cực tiểu của hệ số Wavelet phản ánh điểm uốn và vị trí biến đổi đột ngột của tín hiệu.
  • Số mũ Lipschitz: Tham số định lượng mức độ kỳ dị của tín hiệu, giúp phân biệt đỉnh R thực sự với các xung nhiễu nhọn tần số cao.
  • Khoảng cách QT: Chỉ số sinh học quan trọng phản ánh toàn bộ thời gian khử cực và tái cực tâm thất, dao động bình thường từ 250 ms đến 550 ms.

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

  • Nguồn dữ liệu và kích thước mẫu: Nghiên cứu sử dụng tập dữ liệu tín hiệu điện tim số hóa với tần số lấy mẫu chuẩn 250 Hz (khoảng cách lấy mẫu 4 ms). Tập mẫu bao gồm hàng nghìn chu kỳ tim được trích xuất từ 12 chuyển đạo độc lập, mô phỏng các hình thái điện tim chuẩn và điện tim bệnh lý có kèm nhiễu cơ học, nhiễu trôi nền.
  • Phương pháp chọn mẫu: Lựa chọn mẫu phân tầng theo các phân nhóm hình thái sóng, tập trung vào chuyển đạo D2 (nơi sóng P và QRS có biên độ rõ rệt nhất) và chuyển đạo V2, V4 (tối ưu cho việc phân tích sóng T và khoảng cách QT).
  • Phương pháp phân tích: Tín hiệu được đưa qua dàn lọc phân rã đa tầng (Decimator và Un-decimator) từ mức tỷ lệ $2^1$ đến $2^4$. Tác giả xây dựng thuật toán lọc thông thấp FIR có tần số cắt 16 Hz tại điểm suy giảm 3 dB (độ khuếch đại 25 tại 0 Hz, độ trễ 20 ms) kết hợp hàm biến đổi chiều dài đường cong $L(n)$ với kích thước cửa sổ trượt $w = 160\text{ ms}$ (tương ứng 40 mẫu).
  • Lý do lựa chọn phương pháp: Phép biến đổi Wavelet Spline bậc hai có tính trực giao và đối xứng pha, loại bỏ hiện tượng méo dạng hình học. Phương pháp tách lọc đa tầng tại tỷ lệ $2^1 - 2^4$ cho phép cô lập năng lượng QRS ở tỷ lệ $2^3$ mà không bị ảnh hưởng bởi nhiễu tần số cao tại $2^1$ hay nhiễu trôi đường đồng điện tại $2^5$.
  • Thời gian thực hiện: Toàn bộ quá trình nghiên cứu lý thuyết, phát triển thuật toán và kiểm thử mô phỏng trên nền tảng Matlab được thực hiện trong 12 tháng.

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

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

  1. Phân bố năng lượng phổ của phức bộ QRS: Thực nghiệm chứng minh năng lượng của phức bộ QRS tập trung chủ yếu tại hai mức tỷ lệ đặc trưng là $2^3$ và $2^4$, trong đó mật độ năng lượng đạt cực đại tại tỷ lệ $2^3$ (chiếm hơn 72% tổng năng lượng của phức bộ). Tại tỷ lệ $2^1$, thành phần tín hiệu chủ yếu là nhiễu trắng và dao động tần số cao; ngược lại, tại tỷ lệ lớn hơn $2^5$, ảnh hưởng của hiện tượng trôi đường đẳng điện tăng lên đáng kể.
  2. Độ chính xác vượt trội trong phát hiện đỉnh R: Thuật toán dò tìm đường Modulus cực trị từ tỷ lệ lớn ($2^4$) xuống tỷ lệ nhỏ ($2^1$) kết hợp quy tắc khử cực trị dư thừa (dựa trên tỷ số biên độ - khoảng cách $A/L > 1.2$) giúp định vị đỉnh R chính xác tại điểm đi qua mức 0 của cặp cực trị âm - dương ở tỷ lệ $2^1$. Độ chính xác phát hiện đỉnh R đạt 99.4%, cao hơn 2.4% so với phương pháp lọc thông dải miền thời gian cổ điển (đạt 97.0%).
  3. Hiệu quả phân tách ranh giới sóng P, sóng T và khoảng cách QT: Tại tỷ lệ $2^4$, sóng T tạo ra cặp cực trị rõ nét với độ trễ được hiệu chỉnh chính xác bằng $2^{j-1}-1$ mẫu. Thuật toán kết hợp hàm độ dài $L(n)$ đã xác định chính xác điểm kết thúc của sóng T trong phạm vi độ dài QT chuẩn từ 250 ms đến 550 ms, giảm thiểu sai số nhận dạng sai lệch đỉnh sóng xuống dưới 0.9%.

Thảo luận kết quả

Khả năng phân tách tín hiệu chính xác của mô hình bắt nguồn từ đặc tính toán học của hàm Wavelet Spline. Nhờ đạo hàm bậc nhất của hàm làm trơn, phép biến đổi chuyển các điểm uốn có độ dốc lớn của phức bộ QRS thành các cặp cực trị đối xứng qua điểm 0. Việc phân tích đồng thời trên miền thời gian - tần số khắc phục triệt để nhược điểm của phép biến đổi Fourier vốn làm mất thông tin tức thời.

Khi so sánh với các nghiên cứu sử dụng mạng Nơ-ron nhân tạo (ANN), thuật toán Wavelet thể hiện ưu thế vượt trội về chi phí tính toán: thời gian xử lý một chu kỳ tim giảm khoảng 45% đến 60%, không đòi hỏi pha huấn luyện phức tạp và không phụ thuộc vào dữ liệu mồi. So với biến đổi Fourier thời gian ngắn (STFT) với cửa sổ Gauss cố định, Wavelet tự động thu hẹp cửa sổ thời gian ở vùng tần số cao, cho phép nhận dạng các biến đổi đột ngột của sóng tim với độ phân giải thời gian dưới 4 ms.

Các kết quả định lượng có thể được trực quan hóa qua hai dạng biểu diễn kỹ thuật:

  • Bảng so sánh hiệu suất lọc: Tổng hợp tỷ số tín hiệu trên nhiễu (SNR) và tỷ lệ phát hiện đúng giữa 4 mức tỷ lệ phân rã ($2^1, 2^2, 2^3, 2^4$), làm nổi bật sự vượt trội của tỷ lệ $2^3$ với độ nhạy đạt 99.4%.
  • Đồ thị đáp ứng tần số - biên độ: Biểu diễn đường cong suy giảm 3 dB của bộ lọc FIR tại tần số cắt 16 Hz và đồ thị dịch chuyển pha tuyến tính với độ trễ cố định 5 mẫu, minh chứng khả năng khôi phục nguyên vẹn dạng sóng gốc sau khi khử nhiễu.

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

  1. Tích hợp thuật toán Wavelet Spline vào vi xử lý nhúng của thiết bị Holter ECG 24 giờ:

    • Hành động: Tối ưu hóa mã nguồn thuật toán từ môi trường Matlab sang ngôn ngữ C/C++ nhúng và nạp vào chip DSP hoặc vi điều khiển ARM Cortex-M.
    • Chỉ số mục tiêu: Đạt tốc độ xử lý thời gian thực với độ trễ tính toán dưới 15 ms cho mỗi chu kỳ tim, nâng thời lượng pin thiết bị lên trên 48 giờ hoạt động liên tục.
    • Thời gian thực hiện: Giai đoạn 6 đến 12 tháng tới.
    • Chủ thể thực hiện: Kỹ sư R&D tại các doanh nghiệp sản xuất thiết bị y tế công nghệ cao.

  2. Hoàn thiện mô-đun lọc thích nghi triệt tiêu nhiễu trôi đường đồng điện:

    • Hành động: Xây dựng thuật toán ước lượng đường nền động dựa trên hệ số xấp xỉ ở cấp phân rã thứ 5 ($2^5$) kết hợp bộ lọc thông thấp cắt 0.5 Hz để loại bỏ hoàn toàn nhiễu hô hấp.
    • Chỉ số mục tiêu: Giảm tỷ lệ dương tính giả khi phát hiện biến đổi đoạn ST xuống dưới 0.5%.
    • Thời gian thực hiện: Hoàn thành trong vòng 18 tháng.
    • Chủ thể thực hiện: Nhóm nghiên cứu xử lý tín hiệu sinh học tại các viện nghiên cứu và trường đại học.

  3. Mở rộng cơ sở dữ liệu kiểm chuẩn điện tim 12 chuyển đạo tại các bệnh viện:

    • Hành động: Thiết lập kho dữ liệu ECG số hóa đa dạng bệnh lý (nhồi máu cơ tim, rung nhĩ, ngoại tâm thu thất) với tối thiểu 10.000 bản ghi được gắn nhãn lâm sàng bởi các bác sĩ chuyên khoa.
    • Chỉ số mục tiêu: Độ tin cậy chẩn đoán tự động trên tập dữ liệu thử nghiệm độc lập đạt trên 99.5%.
    • Thời gian thực hiện: Lộ trình 24 tháng (2027 - 2028).
    • Chủ thể thực hiện: Bộ Y tế phối hợp cùng các bệnh viện tim mạch tuyến trung ương.

  4. Phát triển phần mềm hỗ trợ chẩn đoán tự động trên nền tảng y tế từ xa (Telemedicine):

    • Hành động: Đóng gói thuật toán thành các dịch vụ API đám mây, cho phép tự động phân tích điện tim gửi về từ các trạm y tế cơ sở và cảnh báo khẩn cấp các ca loạn nhịp nguy hiểm.
    • Chỉ số mục tiêu: Rút ngắn thời gian trả kết quả phân tích sơ bộ xuống dưới 3 giây/bản ghi.
    • Thời gian thực hiện: Giai đoạn 2028 - 2030.
    • Chủ thể thực hiện: Các công ty công nghệ y tế (MedTech) và cơ sở khám chữa bệnh từ xa.

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

  1. Học viên cao học và nghiên cứu sinh chuyên ngành Kỹ thuật Điện tử, Viễn thông và Kỹ thuật Y sinh:
    • Lợi ích: Nắm vững cơ sở toán học của biến đổi Wavelet liên tục (CWT), biến đổi Wavelet rời rạc (DWT), thuật toán phân tích đa phân giải Mallat và kỹ thuật xử lý tín hiệu y sinh phi dừng.
    • Tình huống ứng dụng: Xây dựng các đề tài nghiên cứu về xử lý ảnh y tế, nén tín hiệu số hoặc thiết kế các bộ lọc thích nghi.
  2. Kỹ sư phát triển phần mềm và phần cứng thiết bị y tế (Biomedical Engineers):
    • Lợi ích: Tiếp cận sơ đồ khối chi tiết của dàn lọc phân rã đa tầng và quy tắc loại bỏ cực trị dư thừa để lập trình trực tiếp lên thiết bị đo ECG cầm tay.
    • Tình huống ứng dụng: Thiết kế firmware cho máy theo dõi bệnh nhân (Patient Monitor) và máy ghi điện tim gắng sức.
  3. Bác sĩ chuyên khoa tim mạch và kỹ thuật viên thăm dò chức năng tim:
    • Lợi ích: Hiểu rõ cơ chế tính toán tự động các thông số hình thái học của 12 chuyển đạo, nguyên nhân gây sai số đường đồng điện và ý nghĩa của các khoảng thời gian sinh học PQ, QT, ST.
    • Tình huống ứng dụng: Phối hợp cùng chuyên gia kỹ thuật để hiệu chỉnh ngưỡng cảnh báo lâm sàng và đánh giá độ tin cậy của phần mềm máy đo ECG.
  4. Giảng viên đại học giảng dạy chuyên đề Xử lý tín hiệu số và Toán ứng dụng:
    • Lợi ích: Khai thác hệ thống công thức giải tích chuẩn xác, các minh họa so sánh giữa Fourier, STFT và Wavelet cùng các đoạn mã giả định vị điểm đặc trưng.
    • Tình huống ứng dụng: Sử dụng làm tài liệu bài giảng mẫu và bài tập lớn mô phỏng Matlab cho sinh viên năm cuối.

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

Phép biến đổi Wavelet vượt trội hơn phép biến đổi Fourier như thế nào trong phân tích ECG?

Phép biến đổi Fourier chỉ phân tích tín hiệu trên miền tần số thuần túy và làm mất hoàn toàn tọa độ thời gian, không thể phát hiện các biến đổi tức thời. Biến đổi Wavelet cho phép khảo sát tín hiệu đồng thời trên cả miền thời gian và miền tần số nhờ khả năng co giãn kích thước cửa sổ linh hoạt, giúp định vị chính xác các đỉnh sóng nhọn có độ dài chỉ vài mili giây với sai số thời gian dưới 4 ms.

Tại sao hàm Wavelet Spline bậc hai được chọn làm hàm cơ sở trong nghiên cứu này?

Hàm Wavelet Spline bậc hai có tính chất đối xứng và độ dịch pha tuyến tính, giúp giữ nguyên vẹn dạng hình học của các đỉnh sóng mà không gây méo tín hiệu sau lọc. Ngoài ra, hàm này là đạo hàm bậc nhất của hàm làm trơn, tạo ra các cặp cực trị âm - dương đối xứng qua mức 0 tại các điểm có độ dốc lớn nhất, rất phù hợp với đặc tính dốc đứng của phức bộ QRS.

Thuật toán phân biệt đỉnh sóng R thực sự với các xung nhiễu tần số cao bằng cách nào?

Thuật toán tính toán số mũ suy giảm Lipschitz $\alpha$ tại các tỷ lệ từ $2^1$ đến $2^4$. Sóng R có năng lượng tập trung tại $2^3$ nên tổng giá trị độ dốc $\alpha_1 + \alpha_2 > 0$. Ngược lại, các xung nhiễu co cơ hoặc nhiễu tần số cao suy giảm rất nhanh ở tỷ lệ lớn, tạo ra $\alpha_1 < 0, \alpha_2 < 0$, giúp loại bỏ hoàn toàn các cực trị giả với tỷ lệ chính xác trên 99.4%.

Tần số lấy mẫu 250 Hz đáp ứng yêu cầu phân tích tín hiệu điện tim như thế nào?

Theo định lý lấy mẫu Nyquist, tần số 250 Hz (chu kỳ lấy mẫu 4 ms) hoàn toàn bao phủ dải tần hữu ích của tín hiệu ECG (từ 0.05 Hz đến 100 Hz). Tần số này đảm bảo độ phân giải thời gian đủ mịn để xác định bề rộng phức bộ QRS (chuẩn 80 ms đến 120 ms) và độ trễ bộ lọc chỉ 5 mẫu (20 ms), đồng thời tối ưu hóa khối lượng tính toán trên hệ thống nhúng.

Khoảng cách QT được xác định tự động bằng phương pháp nào trong luận văn?

Khoảng cách QT được đo tự động qua 3 bước: phát hiện điểm đầu phức bộ QRS tại chuyển đạo V2, chọn chuyển đạo tối ưu có biên độ sóng T nổi trội nhất, và tính toán độ dài đường cong biến đổi $L(n)$ với cửa sổ trượt 160 ms (40 mẫu). Điểm xa nhất từ tín hiệu $L(n)$ đến đường nối chuẩn được xác định là điểm kết thúc sóng T, cho phép đo chính xác khoảng QT trong phạm vi từ 250 ms đến 550 ms.

Kết luận

  • Khẳng định tính ưu việt của công cụ Wavelet: Luận văn đã chứng minh phép biến đổi Wavelet rời rạc (DWT) và cấu trúc đa phân giải Mallat là giải pháp tối ưu để phân tích các tín hiệu y sinh phi dừng, vượt trội hoàn toàn so với biến đổi Fourier cổ điển.
  • Xây dựng thành công hệ thuật toán nhận dạng sóng tim: Nghiên cứu đã phát triển hoàn chỉnh chuỗi thuật toán phát hiện phức bộ QRS, sóng P, sóng T và đo khoảng cách QT đạt độ chính xác trên 99.2% trên tập dữ liệu chuẩn 250 Hz.
  • Tối ưu hóa khả năng kháng nhiễu: Giải quyết triệt để bài toán loại bỏ nhiễu tần số cao tại tỷ lệ $2^1$ và nhiễu trôi đường đồng điện tại tỷ lệ $2^5$ nhờ việc phân lập năng lượng sóng tim tại các tỷ lệ đặc trưng $2^3$ và $2^4$.
  • Cung cấp mô hình toán học khả thi cho công nghệ nhúng: Thuật toán lọc FIR pha tuyến tính kết hợp hàm Wavelet Spline bậc hai có cấu trúc tính toán tinh gọn, phù hợp để chuyển đổi lên các dòng vi xử lý DSP thời gian thực.
  • Lộ trình phát triển tiếp theo: Trong giai đoạn 2027 - 2030, các nhóm nghiên cứu cần đẩy mạnh việc chuyển giao giải thuật vào các thiết bị Holter ECG thương mại và tích hợp trí tuệ nhân tạo để mở rộng khả năng chẩn đoán tự động trên hệ thống y tế số.

Quý độc giả, nhà nghiên cứu và các kỹ sư y sinh quan tâm có thể khai thác toàn văn tài liệu luận văn để ứng dụng trực tiếp các mô hình toán học và thuật toán mô phỏng vào các dự án phát triển thiết bị chăm sóc sức khỏe tim mạch hiện đại.