Tổng quan luận án

Tính cấp thiết và khoảng trống nghiên cứu

Bệnh lý tim mạch đang gia tăng và là nguyên nhân gây tử vong hàng đầu trên toàn cầu, đặc biệt tập trung tại các quốc gia đang phát triển theo thống kê của Tổ chức Y tế Thế giới (WHO). Mỗi năm trên thế giới ghi nhận hàng chục triệu người mắc bệnh tim mạch [43]; nhiều trường hợp rơi vào tình trạng nguy kịch như nhồi máu cơ tim hoặc đột quỵ do việc chẩn đoán và điều trị diễn ra muộn. Tại Việt Nam, thống kê cho thấy cứ 3 người trưởng thành thì có 1 người có nguy cơ mắc bệnh tim mạch. Dự báo từ Hội Tim mạch học Việt Nam chỉ ra rằng đến năm 2017, tỷ lệ dân số mắc các bệnh tim mạch và tăng huyết áp sẽ chiếm khoảng 1/5. Xu hướng bệnh lý tim mạch diễn biến phức tạp với nhiều biến chứng, đồng thời độ tuổi trung bình của bệnh nhân có xu hướng giảm dần về phía lực lượng lao động, dẫn đến tình trạng quá tải nghiêm trọng tại các bệnh viện tuyến trung ương trong bối cảnh cơ sở vật chất và nguồn nhân lực y tế chuyên khoa chưa đáp ứng kịp thời.

Khám bệnh tim mạch thông qua ghi đo tín hiệu điện tim (Electrocardiography - ECG) là quy trình cơ bản nhằm phát hiện các bệnh lý không thường xuyên như thiếu máu cục bộ cơ tim hay loạn nhịp tim. Do đó, việc ứng dụng dịch vụ y tế từ xa (telemedicine) để theo dõi và chăm sóc sức khỏe bệnh nhân tại nhà hoặc di động đang trở thành hướng đi cần thiết, giúp giảm chi phí đi lại, hạn chế ách tắc giao thông và giảm tải áp lực cho hệ thống bệnh viện. Tuy nhiên, tại Việt Nam chưa có công trình nghiên cứu tổng thể và toàn diện nào về việc hỗ trợ theo dõi, chăm sóc bệnh nhân tim mạch tại các bệnh viện có chuyên khoa tim mạch từ tuyến Trung ương đến tuyến Tỉnh/Thành phố sử dụng công nghệ vô tuyến.

Khi truyền tín hiệu ECG thời gian thực qua môi trường vô tuyến từ thiết bị số trợ giúp cá nhân (Personal Digital Assistant - PDA) dùng pin, hệ thống đối mặt với hai thách thức kỹ thuật lớn:

  • Tín hiệu ECG thu nhận từ cơ thể người bị can nhiễu nghiêm trọng, đặc biệt là nhiễu nguồn điện xoay chiều (50Hz/60Hz) và nhiễu trôi dạt đường cơ sở (baseline wander) do hoạt động hô hấp và thay đổi trở kháng tiếp xúc điện cực - da.
  • Dữ liệu ECG đa đạo trình có dung lượng lớn (ví dụ tín hiệu 12 đạo trình, tần số lấy mẫu 360Hz, độ phân giải 11-bit ghi trong 1 ngày tạo ra dung lượng khoảng 746MB; nếu tần số lấy mẫu 1000Hz thì dung lượng lên tới 2072MB/ngày). Khi kích thước gói tin lớn, việc truyền qua các mạng vô tuyến băng thông hẹp hoặc không ổn định dễ gây mất gói, tăng thời gian trễ, tăng năng lượng tiêu thụ của khối phát PDA và tốn kém tài nguyên lưu trữ lâu dài [4], [7], [24], [44]. Do đó, việc xây dựng giải pháp lọc nhiễu kết hợp nén dữ liệu điện tim trước khi truyền vô tuyến là vấn đề cấp thiết cần giải quyết.

Mục tiêu nghiên cứu của luận án

Luận án xác định mục tiêu tổng quát là xây dựng gói giải pháp liên hoàn cho việc truyền tín hiệu ECG bằng công nghệ không dây, bao gồm 3 mục tiêu cụ thể:

  1. Nâng cao chất lượng lọc nhiễu trôi dạt đường cơ sở và nhiễu nguồn xoay chiều có trong tín hiệu ECG (hai loại nhiễu ảnh hưởng nhiều nhất đến chất lượng thu nhận tín hiệu ECG).
  2. Đề xuất phương pháp nén tín hiệu điện tim nhằm nâng cao hiệu quả và chất lượng truyền tín hiệu, giúp việc theo dõi và chẩn đoán bệnh tim mạch được chính xác và hiệu quả.
  3. Đánh giá hiệu quả của thuật toán nén hai trạng thái tín hiệu điện tim dựa trên hệ thống truyền - nhận trong môi trường thực ứng dụng công nghệ WiFi.

Để hoàn thành các mục tiêu trên, luận án tập trung giải quyết 3 vấn đề chính:

  1. Nghiên cứu phương pháp lọc hiệu quả nhiễu nguồn xoay chiều và nhiễu trôi dạt đường cơ sở có trong tín hiệu điện tim.
  2. Nghiên cứu phương pháp nén dữ liệu điện tim trước khi truyền và giải nén sau khi thu nhận nhằm giảm lỗi mất gói trên đường truyền và giảm năng lượng phát của PDA, đảm bảo phía thu nhận khôi phục chính xác tín hiệu.
  3. Xây dựng và thử nghiệm mô hình hệ thống truyền - nhận tín hiệu ECG sử dụng công nghệ WiFi khi truyền tín hiệu ECG được nén và giải nén bằng thuật toán đề xuất; đánh giá thuật toán nén hai trạng thái trong môi trường thời gian thực và điều kiện thực tế.

Đối tượng và phạm vi nghiên cứu

  • Đối tượng nghiên cứu: Tín hiệu điện tim (ECG), các dạng nhiễu can thiệp vào tín hiệu ECG (nhiễu nguồn xoay chiều 50Hz/60Hz, nhiễu trôi dạt đường cơ sở, nhiễu trắng), các giải pháp xử lý lọc số, thuật toán nén dữ liệu ECG và các giao thức truyền thông vô tuyến (WiFi, Zigbee, Bluetooth, GPRS/3G).
  • Phạm vi không gian và công nghệ: Xử lý tiền xử lý tín hiệu tại thiết bị PDA đeo trên người; truyền nhận tín hiệu ECG qua mạng không dây cục bộ (WiFi IEEE 802.11) và kết nối Internet tới máy chủ trung tâm hoặc trạm giám sát của bác sĩ.
  • Phạm vi dữ liệu: Dữ liệu mô phỏng toán học trên phần mềm LabVIEW và các bộ cơ sở dữ liệu điện tim chuẩn quốc tế từ PhysioNet/MIT-BIH, bao gồm: MIT-BIH Arrhythmia Database (với các bản ghi chuẩn như 103, 105, 106, 220 và toàn bộ 48 bản ghi loạn nhịp), MIT-BIH ST Change Database (các bản ghi 301, 302) cùng 9 bản ghi nhịp nhanh tâm thất (Malignant Ventricular Ectopy Database).

Tổng quan tài liệu và vị trí của luận án

Các hướng nghiên cứu trước được luận án tổng quan

Tác giả đã hệ thống hóa các nghiên cứu trong nước và quốc tế theo các nhóm nội dung:

  1. Công nghệ vô tuyến trong truyền dẫn tín hiệu y sinh và ECG:

    • Chuẩn Zigbee (IEEE 802.15.4): Hoạt động ở các dải tần 868MHz (Châu Âu, 20kbps), 915MHz (Mỹ/Nhật, 40kbps) và 2.4GHz (toàn cầu, 250kbps) [69]. Các nghiên cứu [69], [77] sử dụng phần cứng Xbee Pro kết hợp LabVIEW và Visual Basic để truyền tín hiệu ECG mô phỏng. Luận án chỉ ra Zigbee có ưu điểm chi phí thấp, công suất tiêu thụ thấp nhưng nhược điểm là tốc độ dữ liệu thấp, trễ đường truyền cao, dễ bị can nhiễu bởi các mạng IEEE 802.11 và thiết bị vi sóng [12].
    • Chuẩn Bluetooth (IEEE 802.15.1): Băng thông từ 1-3Mbps (phiên bản 2.0), cự ly ngắn 2-10m. Rasid và cộng sự [52] cùng các nghiên cứu [12], [20] đã áp dụng mô hình Master-Slave dùng Bluetooth truyền từ PDA đến điện thoại thông minh, sau đó kết nối GPRS về máy chủ bệnh viện. Nhược điểm là băng thông hẹp, cự ly ngắn và dễ giao thoa sóng vô tuyến.
    • Chuẩn WiFi (IEEE 802.11b/g/n): Tốc độ từ 11Mbps đến 600Mbps trên dải tần ISM 2.4GHz và 5GHz [30]. CaiKen và cộng sự [15] xây dựng hệ thống truyền - nhận ECG qua WiFi nhưng mới dừng ở mức thu nhận và hiển thị, chưa đánh giá trong môi trường thực. Kyungtae Kang và cộng sự [45], [46] can thiệp vào lớp điều khiển truy nhập môi trường (MAC) để nâng cao chất lượng dịch vụ (QoS) và khắc phục trượt bit. Tài liệu [44] đề xuất mô hình mạng thích nghi để truyền ECG thời gian thực nhưng nhấn mạnh thách thức dung lượng dữ liệu 500MB/ngày/bệnh nhân gây biến động băng thông và tăng tỷ lệ lỗi kênh.
    • Mạng di động GPRS/3G: Tốc độ GPRS đạt 40-55kbps, 3G đạt 144kbps - 2Mbps [27], khi di chuyển giảm xuống 128kbps [60]. Các công trình [52], [59], [81] ứng dụng GPRS/3G cho y tế từ xa nhưng chi phí đường truyền tính theo dung lượng tải khá cao và dễ xảy ra nghẽn mạng ở vùng sóng yếu.

  2. Kỹ thuật lọc nhiễu nguồn xoay chiều (Power-line Interference):

    • Các thuật toán lọc thích nghi: Bình phương trung bình nhỏ nhất (LMS), bình phương trung bình nhỏ nhất chuẩn hóa (NLMS), bình phương nhỏ nhất đệ quy (RLS) được nghiên cứu bởi I. Romeo [34], Sachin Singh và cộng sự [68], Syed A. Rehman [75] và các tác giả [2], [6], [17], [28], [74]. RLS có tốc độ hội tụ nhanh nhưng độ phức tạp tính toán cao và kém ổn định so với LMS/NLMS. Syed A. Rehman [75] thử nghiệm lọc nhiễu 60Hz (1mV) trên các bản ghi MIT-BIH 103, 105, 106 cho thấy tỷ số tín hiệu trên nhiễu (SNR) trung bình đạt 9.42 dB (LMS/NLMS) và 11.23 dB (DLMS).
    • Bộ lọc triệt tần thích nghi: Nghiên cứu [2] kết hợp bộ lọc triệt tần với tiếp cận sóng nhỏ (wavelet) và biến đổi Fourier để làm trơn phổ tín hiệu, đánh giá qua sai số bình phương trung bình (MSE). Tuy nhiên, khâu dò tìm tần số nhiễu trong [2] phức tạp và hiệu quả lọc suy giảm khi biên độ nhiễu lớn (chiếm 30%-40% biên độ tín hiệu ECG).

  3. Kỹ thuật lọc nhiễu trôi dạt đường cơ sở (Baseline Wander):

    • Bộ lọc thông cao cổ điển với tần số cắt khoảng 0.5Hz làm méo dạng sóng ECG, đặc biệt là đoạn ST [39].
    • Lọc thích nghi [47], [64], [68] gặp khó khăn trong việc thiết lập tín hiệu tham chiếu phù hợp.
    • Phân tích biến đổi Wavelet [10], [18], [33], [47], [63]: Tác giả luận án đã phân tích trong công bố tại Hội nghị IEEE ICUFN-2013 [công bố số 2], chỉ ra rằng độ chính xác lọc của wavelet phụ thuộc chặt chẽ vào tần số ngưỡng và biên độ nhiễu; do nhiễu trôi dạt không ổn định nên việc chọn ngưỡng cố định rất khó đạt hiệu quả tối ưu trong thời gian thực.
    • Xử lý trong miền thời gian: Jungkuk Kim và cộng sự [39] loại bỏ độ lệch sóng (epoch) tần số cao bằng 4 sóng đánh dấu độ dốc nhưng không đưa ra tiêu chí định lượng đánh giá chất lượng; A. Fasano và cộng sự [8] dùng phương pháp giảm biến bậc 2; Na Pan và cộng sự [54] dùng phương pháp phân tách theo kinh nghiệm (Empirical Mode Decomposition - EMD) trên cơ sở dữ liệu MIT-BIH, đạt hệ số tương quan trung bình 0.9746 trên bản ghi 301.

  4. Tác động của kích thước gói tin đến truyền dẫn không dây:

    • Jari Korhonen và Ye Wang [35] chứng minh kích thước gói tin lớn trong môi trường không dây tốc độ thấp làm tăng đáng kể trễ và tỷ lệ mất gói do phải truyền lại nhiều lần. Yaakob và cộng sự [55] khảo sát kích thước gói từ 30 đến 130 byte, chỉ ra sự cần thiết phải giảm kích thước gói khi tỷ lệ lỗi bit (BER) cao. Liang và Balasingham [53] lưu ý nếu kích thước gói quá nhỏ sẽ làm tăng tần suất phát của các nút mạng, gây lãng phí tài nguyên kênh.

Khoảng trống luận án lựa chọn giải quyết

Luận án tập trung giải quyết khoảng trống về một giải pháp tích hợp hoàn chỉnh: (1) Nâng cao hiệu quả lọc nhiễu nguồn 50Hz/60Hz bằng cách đơn giản hóa việc dò tìm tần số qua biến đổi Fourier nhanh và áp dụng lọc nhiều vòng lặp với ngưỡng dừng biên độ; (2) Lọc triệt để nhiễu trôi dạt đường cơ sở trong miền thời gian dựa trên các điểm đẳng thế sinh học kết hợp nội suy Cubic Spline; (3) Đề xuất phương pháp nén dữ liệu hai trạng thái phù hợp để giảm kích thước gói tin, giảm lỗi mất gói và giảm trễ khi truyền tín hiệu ECG thời gian thực qua hạ tầng mạng WiFi.

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

Khung lý thuyết và khái niệm cốt lõi

  1. Đặc tính tín hiệu điện tim (ECG): Tín hiệu ECG ghi nhận hoạt động điện sinh học của tim, bao gồm 5 sóng thành phần cơ bản: P, Q, R, S, T và đôi khi có sóng U (xuất hiện khi nhịp tim < 60 nhịp/phút ghi ở đạo trình trước tim bên phải) [23], [61].
    • Dải tần số chẩn đoán y khoa: 0.05Hz $\div$ 100Hz.
    • Dải tần số giám sát monitor: 0.5Hz $\div$ 40Hz.
Sóng / Đoạn thành phần Biên độ chuẩn (mV) Khoảng thời gian chuẩn (s)
Sóng P < 0.25 0.16
Phức bộ QRS 1.60 0.08 $\div$ 0.10
Mức ST 0.00 0.05 $\div$ 0.15
Đoạn PR - 0.12 $\div$ 0.20
Đoạn QT - 0.35 $\div$ 0.44
Sóng T 0.10 $\div$ 0.50 -

  1. Mô hình toán học lọc nhiễu nguồn xoay chiều thích nghi: Tín hiệu đo được $S(n) = U(n) + N(n)$, trong đó $U(n)$ là tín hiệu ECG gốc và $N(n)$ là nhiễu xoay chiều 50Hz/60Hz. Biến đổi Fourier rời rạc của tín hiệu $X(n)$ được xác định bởi: $$X(\omega) = \sum_{n=-\infty}^{\infty} X(n) e^{-i\omega n}$$ Hàm truyền đạt của bộ lọc triệt tần thích nghi bậc 2 được xây dựng dựa trên công trình của B. Widrow [14]: $$H(z) = \frac{1 + b_1 z^{-1} + z^{-2}}{1 + a_1 z^{-1} + 0.97 z^{-2}}$$ với các hệ số được tính toán theo tần số nhiễu $\omega_0$: $a_0 = 1$; $b_0 = 1$; $a_1 = 1.97\cos\omega_0$; $b_1 = -2\cos\omega_0$.

  2. Lý thuyết điểm đẳng thế (Isoelectric Points) và nội suy Cubic Spline: Điểm đẳng thế bao gồm: khởi đầu sóng P ($P_{onset}$), kết thúc sóng P ($P_{offset}$), khởi đầu phức bộ QRS ($QRS_{onset}$), kết thúc phức bộ QRS ($QRS_{offset}$) và kết thúc sóng T ($T_{offset}$). Tại các thời điểm này, hoạt động điện tim nghỉ, điện áp chênh lệch bằng 0. Khi có nhiễu trôi dạt đường cơ sở, biên độ các điểm này lệch khỏi 0.

    • Nhiễu trôi dạt đường cơ sở được định nghĩa là đường nội suy Cubic Spline đi qua tập hợp các điểm đẳng thế tìm được.
    • Các tính chất toán học xác định điểm đẳng thế:
      • Tính chất 1: $QRS_{onset}$ và $QRS_{offset}$ là điểm uốn, tương ứng với điểm vượt 0 của đạo hàm bậc 2 ($QRS_{onset}$ vượt 0 từ âm lên dương; $QRS_{offset}$ vượt 0 từ dương xuống âm).
      • Tính chất 2: $T_{offset}$ là điểm cực tiểu sau sóng T, tương ứng điểm vượt 0 từ âm lên dương của đạo hàm bậc 1; $P_{onset}$ là điểm vượt 0 từ dương của đạo hàm bậc 1.
      • Tính chất 3: $T_{offset}$ là điểm vượt 0 từ âm lên dương của đạo hàm bậc 2 (sau đỉnh T); $P_{onset}$ là điểm vượt 0 từ dương của đạo hàm bậc 2 (trước đỉnh P).
    • Đoạn nội suy $Y$ trên khoảng $[X_i, X_{i+1}]$ được tính bằng: $$y = A y_i + B y_{i+1} + C y''i + D y''{i+1}$$ trong đó $A = \frac{X_{i+1} - X}{X_{i+1} - X_i}$; $B = 1 - A$; $C = \frac{1}{6}(A^3 - A)(X_{i+1} - X_i)^2$; $D = \frac{1}{6}(B^3 - B)(X_{i+1} - X_i)^2$.

Phương pháp nghiên cứu thực tế

  • Phương pháp tiếp cận: Kết hợp nghiên cứu lý thuyết giải tích, mô phỏng số học thuật toán và thử nghiệm đo kiểm thực nghiệm.
  • Công cụ phần mềm và xử lý: Môi trường lập trình LabVIEW (sử dụng các công cụ Spline Interpolation.VI, FFT, bộ lọc trung bình làm mịn), ngôn ngữ C và Visual C để lập trình module truyền thông và giao diện hiển thị.
  • Tiêu chí đánh giá định lượng:
    • Sai số bình phương trung bình (MSE): $$MSE = \frac{1}{N} \sum_{i=1}^{N} (x_i - \hat{x}_i)^2$$
    • Tỷ số tín hiệu trên nhiễu (SNR, tính theo dB): $$SNR = 10 \log_{10} \left( \frac{\sum_{n} U^2(n)}{\sum_{n} [S(n) - U(n)]^2} \right)$$
    • Hệ số tương quan chéo ($r$): $$r = \frac{\sum_{n=1}^N x(n)y(n)}{\sqrt{\sum_{n=1}^N x^2(n) \cdot \sum_{n=1}^N y^2(n)}}$$

Nguồn dữ liệu thực nghiệm

  1. Tín hiệu mô phỏng LabVIEW: Nhịp tim 72 nhịp/phút (1.2Hz), tần số lấy mẫu $f_s = 360$Hz, độ dài 10 giây, biên độ đỉnh 1.6mV, cộng nhiễu trôi dạt tần số $\le 1.2$Hz biên độ 0.5mV và nhiễu trắng $0.05 \div 0.2$mV.
  2. Dữ liệu chuẩn quốc tế MIT-BIH: Các bản ghi 103, 105, 106, 220 (thuộc MIT-BIH Arrhythmia Database); các bản ghi 301, 302 (thuộc MIT-BIH ST Change Database); 48 bản ghi loạn nhịp và 9 bản ghi nhịp nhanh tâm thất.

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

Chương 1: Tổng quan về thu nhận, xử lý và truyền tín hiệu điện tim trong hệ thống theo dõi và hỗ trợ chẩn đoán bệnh nhân tim mạch

Chương 1 hệ thống hóa cơ sở sinh lý học của tín hiệu điện tim và các đạo trình đo chuẩn của Willem Einthoven [61]. Tác giả phân tích các dạng méo tín hiệu điện tim tương ứng với các bệnh lý thường gặp như nhịp nhanh thất, block nhĩ thất độ 3, nhồi máu cơ tim và cơn đau thắt ngực (đoạn ST chênh lên hoặc chênh xuống). Chương này phân loại chi tiết các nguồn nhiễu ảnh hưởng đến ECG: nhiễu nguồn xoay chiều (biên độ chiếm 30% - 50% biên độ đỉnh - đỉnh của ECG), nhiễu trôi dạt đường cơ sở (tần số 0.05Hz đến dưới 0.8Hz do nhịp thở), nhiễu điện cơ và nhiễu tiếp xúc.

Công nghệ vô tuyến Chuẩn IEEE Băng thông / Tốc độ Dải tần hoạt động Khoảng cách truyền
Zigbee IEEE 802.15.4 20 kbps $\div$ 250 kbps 868MHz, 915MHz, 2.4GHz 10 $\div$ 100 m
Bluetooth IEEE 802.15.1 1 Mbps $\div$ 3 Mbps 2.4 GHz 2 $\div$ 10 m
WiFi IEEE 802.11b/g/n 11 Mbps $\div$ 600 Mbps 2.4 GHz, 5 GHz 30 $\div$ 100 m
GPRS / 3G Mạng tế bào 40 kbps $\div$ 2 Mbps Băng tần viễn thông Toàn mạng di động

Chương 1 phân tích cấu trúc hệ thống theo dõi bệnh nhân có dây (monitor đầu giường nối tập trung) và không dây (hệ thống monitor trung tâm vô tuyến 4-16 máy trạm; hệ thống theo dõi di động kết hợp PDA, WiFi, GPRS/3G và Internet). Trên cơ sở tính toán dung lượng dữ liệu ECG phát sinh và các hạn chế về băng thông, trễ gói và năng lượng nguồn pin, chương 1 xác lập định hướng nghiên cứu trọng tâm vào hai khối: xử lý lọc sạch hai loại nhiễu chính và nén dữ liệu ECG trước khi phát sóng vô tuyến.

Chương 2: Đề xuất giải pháp nâng cao hiệu quả lọc nhiễu trong tín hiệu điện tim

Chương 2 trình bày hai giải pháp kỹ thuật cụ thể:

1. Lọc nhiễu nguồn xoay chiều 50Hz/60Hz bằng bộ lọc triệt tần thích nghi FFT nhiều vòng lặp

  • Quy trình thực hiện: Sử dụng biến đổi Fourier nhanh (FFT) trên tín hiệu $S(n)$ để dò tìm chính xác tần số nhiễu có biên độ cực đại trong dải hẹp (47-53Hz đối với nguồn 50Hz; 57-63Hz đối với nguồn 60Hz). Tính toán hàm truyền $H(n)$ tương ứng và thực hiện lọc $S(n) = S(n) * H(n)$. Quá trình lọc lặp lại nhiều vòng cho đến khi biên độ nhiễu sau lọc nhỏ hơn ngưỡng thiết lập.
  • Khảo sát ngưỡng dừng lọc: Khảo sát ngưỡng từ 0.5 $\mu$V đến 30 $\mu$V trên 3 bản ghi 103, 105, 106 với nhiễu cộng 1mV cho thấy:
    • Với nhiễu 50Hz: Ngưỡng tối ưu đạt 15 $\mu$V $\div$ 18.5 $\mu$V, cho SNR trung bình cực đại đạt 15.35 dB.
    • Với nhiễu 60Hz: Ngưỡng tối ưu đạt 20 $\mu$V, cho SNR trung bình cực đại đạt 16.27 dB.
Ngưỡng sau lọc ($\mu$V) SNR trung bình lọc 50Hz (dB) SNR trung bình lọc 60Hz (dB)
0.5 13.82 12.63
5.0 14.30 14.20
10.0 14.85 15.15
15.0 15.35 15.80
18.5 15.35 -
20.0 15.10 16.27
25.0 14.45 15.40
30.0 13.80 14.60
  • So sánh với các nghiên cứu trước:
    • So sánh với Syed A. Rehman [75] trên 3 bản ghi 103, 105, 106 (nhiễu 60Hz biên độ 1mV): LMS đạt 9.42 dB, NLMS đạt 9.42 dB, DLMS đạt 11.23 dB, trong khi phương pháp đề xuất đạt 16.27 dB.
    • So sánh với nghiên cứu [2] trên bản ghi 220mitdb (nhiễu 50Hz): Khi biên độ nhiễu tăng cao ($\ge 0.4$mV), phương pháp đề xuất duy trì MSE thấp vượt trội nhờ cơ chế lọc nhiều vòng lặp.
Biên độ nhiễu 50Hz (mV) MSE nghiên cứu [2] Thiết lập ngưỡng ($\mu$V) Số vòng lặp đề xuất MSE phương pháp đề xuất
0.1 $11 \times 10^{-5}$ 15 1 $11 \times 10^{-5}$
0.2 $12 \times 10^{-5}$ 15 1 $12 \times 10^{-5}$
0.3 $15 \times 10^{-5}$ 15 1 $15 \times 10^{-5}$
0.4 $49 \times 10^{-3}$ 15 2 $40 \times 10^{-4}$
0.5 $77 \times 10^{-3}$ 15 3 $48 \times 10^{-4}$

2. Lọc nhiễu trôi dạt đường cơ sở trong miền thời gian dựa trên các điểm đẳng thế

  • Thuật toán 5 bước:
    • Bước 1: Giảm mẫu tín hiệu $S(t)$, tính đạo hàm cấp 1 và cấp 2 để triệt giảm ảnh hưởng của nhiễu trắng $W(t)$.
    • Bước 2: Tìm sườn dốc cực trị trước/sau đỉnh R từ đạo hàm cấp 1; xác định các điểm vượt 0 của đạo hàm cấp 1 và cấp 2 để định vị sườn âm - dương của phức bộ QRS.
    • Bước 3: Tìm các đỉnh P, Q, S, T từ các điểm vượt 0 của đạo hàm cấp 1.
    • Bước 4: Tìm các điểm $P_{onset}$, $T_{offset}$ (điểm uốn đạo hàm cấp 1 và 2), $QRS_{onset}$ (vượt 0 âm của đạo hàm cấp 2), $QRS_{offset}$ (vượt 0 dương của đạo hàm cấp 2).
    • Bước 5: Kiểm tra điều kiện sinh học loại trừ sai lệch do nhiễu ($P_{onset}$ nằm ở nửa sau chu kỳ T-T; $T_{offset}$ nằm ở 2/3 phần trước chu kỳ T-T; $QRS_{onset}$ nằm sau P và trước T).
  • Nội suy và trừ nhiễu: Áp dụng công cụ Spline Interpolation.VI nội suy đường cong $B'(t)$ từ các điểm đẳng thế và lấy $P(t) = S(t) - B'(t)$.
  • Kết quả thực nghiệm:
    • Độ chính xác nhận dạng tại mức nhiễu trắng 0.05mV: MSE xác định đỉnh R là 1, đỉnh S là 7, đỉnh Q là 9, đỉnh P là 77, đỉnh T là 84; MSE xác định $QRS_{offset}$ là 8, $QRS_{onset}$ là 11, $P_{onset}$ là 53, $T_{offset}$ là 59.
    • Tín hiệu sau khi lọc nhiễu trôi dạt đạt sai số MSE = 0.0039 $\text{mV}^2$.
    • So sánh với phương pháp EMD của Na Pan và cộng sự [54] trên bản ghi 301 MIT-BIH: Hệ số tương quan trung bình của phương pháp đề xuất đạt 0.9984 (so với 0.9746 của Na Pan). Thử nghiệm trên bản ghi 301 và 302 đạt MSE = 0.0018 $\text{mV}^2$ và SNR = 26.54 dB.

Chương 3: Đề xuất phương pháp nén dữ liệu điện tim (ECG) ứng dụng trong hệ thống theo dõi và chẩn đoán bệnh tim mạch từ xa

Chương 3 phân tích cơ sở lý thuyết về mối quan hệ giữa kích thước gói tin, tỷ lệ lỗi bit (BER), tỷ lệ lỗi gói (PER) và độ trễ truyền dẫn không dây dựa trên các công trình của Jari Korhonen & Ye Wang [35], Yaakob [55] và Liang & Balasingham [53]. Tác giả đề xuất thuật toán nén hai trạng thái mới dành cho tín hiệu ECG, tối ưu hóa kích thước gói tin để truyền qua mạng vô tuyến. Thuật toán được thử nghiệm và đối sánh trên hai bộ cơ sở dữ liệu chuẩn gồm 48 bản ghi loạn nhịp và 9 bản ghi nhịp nhanh tâm thất từ MIT-BIH.

Chương 4: Thiết lập hệ thống thử nghiệm trong môi trường thực tế để truyền-nhận tín hiệu ECG sử dụng công nghệ WiFi

Chương 4 mô tả thiết lập hệ thống thử nghiệm truyền - nhận thực tế qua mạng WiFi kết nối Internet. Tác giả tiến hành đo kiểm, so sánh hai luồng dữ liệu: dữ liệu điện tim không nén và dữ liệu điện tim được nén bằng thuật toán hai trạng thái đề xuất; đánh giá chất lượng tín hiệu ECG thu nhận thông qua các tham số kỹ thuật về tỷ lệ mất gói, độ trễ thời gian thực và độ trung thực của dạng sóng sau giải nén.

Kết quả và những đóng góp mới

Các đóng góp mới của luận án

  1. Về mặt lý luận và thuật toán:
    • Đề xuất phương pháp lọc nhiễu nguồn xoay chiều kết hợp biến đổi Fourier nhanh (FFT) để quét tần số tự động trong dải hẹp (47-53Hz hoặc 57-63Hz) và áp dụng cơ chế lọc triệt tần thích nghi nhiều vòng lặp theo ngưỡng biên độ dư (15 $\mu$V cho 50Hz; 20 $\mu$V cho 60Hz), giải quyết triệt để tình trạng méo tín hiệu khi biên độ nhiễu lớn.
    • Đề xuất phương pháp lọc nhiễu trôi dạt đường cơ sở trong miền thời gian dựa trên các điểm đẳng thế sinh học ($P_{onset}$, $P_{offset}$, $QRS_{onset}$, $QRS_{offset}$, $T_{offset}$) kết hợp nội suy Cubic Spline, giúp giữ nguyên vẹn đoạn sóng ST mà không làm thay đổi các đặc trưng bệnh lý tim mạch.
    • Đề xuất phương pháp nén hai trạng thái tín hiệu ECG, giúp giảm nhỏ kích thước gói tin trước khi đưa lên kênh truyền không dây.
  2. Về mặt thực tiễn và ứng dụng:
    • Xây dựng hoàn chỉnh bộ công cụ phần mềm xử lý lọc nhiễu và nén tín hiệu ECG trên nền tảng LabVIEW và C/Visual C.
    • Triển khai thử nghiệm thành công mô hình truyền - nhận tín hiệu ECG thời gian thực qua mạng WiFi trong môi trường thực tế, chứng minh tính khả thi của việc ứng dụng hệ thống theo dõi và hỗ trợ chẩn đoán bệnh nhân tim mạch tại các cơ sở y tế và gia đình tại Việt Nam.

Kiến nghị và giải pháp đề xuất

  • Ứng dụng giải pháp lọc nhiễu và nén hai trạng thái vào khối tiền xử lý của các thiết bị PDA di động gắn trên người bệnh nhân nhằm tiết kiệm năng lượng pin, giảm lưu lượng đường truyền và tăng độ tin cậy kết nối.
  • Sử dụng thuật toán nén để tối ưu hóa hệ thống lưu trữ dữ liệu hồ sơ bệnh án điện tử ECG lâu dài tại các bệnh viện, giúp tiết kiệm đáng kể dung lượng ổ đĩa cứng.

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

  • Hạn chế: Phạm vi đánh giá thực nghiệm của hệ thống truyền nhận vô tuyến trong luận án tập trung chủ yếu vào cơ sở hạ tầng mạng WiFi cục bộ và Internet; các thử nghiệm về lọc nhiễu và nén dữ liệu được thực hiện trên các tập dữ liệu mô phỏng LabVIEW và cơ sở dữ liệu chuẩn MIT-BIH (48 bản ghi loạn nhịp, 9 bản ghi nhịp nhanh thất, các bản ghi 103, 105, 106, 220, 301, 302).
  • Hướng nghiên cứu tiếp:
    • Mở rộng thử nghiệm hệ thống truyền nhận trên các môi trường di động diện rộng đa dạng hơn (kết hợp chuyển vùng giữa WiFi, 3G/4G).
    • Tích hợp sâu các thuật toán xử lý vào chip phần cứng nhúng (như ARM hoặc vi điều khiển chuyên dụng trong PDA) để tối ưu hóa thời gian thực và giảm thiểu tối đa kích thước vật lý của thiết bị đeo.

Giá trị tham khảo

  • Đối với nghiên cứu sinh, học viên và giảng viên ngành Kỹ thuật Y sinh, Điện tử - Viễn thông và Công nghệ Thông tin: Luận án cung cấp tài liệu tham khảo chi tiết về đặc tính giải phẫu điện tim, mô hình toán học của các bộ lọc thích nghi, kỹ thuật nhận dạng điểm đẳng thế dựa trên đạo hàm bậc 1 và bậc 2, cùng các phân tích định lượng về hiệu năng mạng truyền thông vô tuyến trong y tế.
  • Đối với các kỹ sư phát triển thiết bị y tế và hệ thống Telemedicine: Các sơ đồ khối thiết kế PDA, lưu đồ thuật toán lọc nhiễu 50Hz/60Hz nhiều vòng lặp và mã nguồn xử lý tín hiệu trên LabVIEW/C tại Chương 2, Chương 3 và Chương 4 là tài liệu kỹ thuật có giá trị ứng dụng trực tiếp.
  • Đối với các nhà quản lý y tế và bệnh viện: Cung cấp cơ sở khoa học để xây dựng mô hình mạng giám sát bệnh nhân tim mạch không dây, hỗ trợ giảm tải cho các bệnh viện tuyến trung ương.

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

1. Tại sao nhiễu nguồn xoay chiều và nhiễu trôi dạt đường cơ sở lại là hai loại nhiễu ảnh hưởng nghiêm trọng nhất đến tín hiệu ECG?

Tần số của hai loại nhiễu này nằm xen kẽ trực tiếp trong dải tần hoạt động của tín hiệu ECG (dải chẩn đoán từ 0.05Hz đến 100Hz). Nhiễu nguồn xoay chiều (50Hz hoặc 60Hz) có biên độ lớn (chiếm 30% - 50% biên độ đỉnh - đỉnh của ECG), dễ làm biến dạng sóng P, sóng T, phức bộ QRS và đoạn QT. Nhiễu trôi dạt đường cơ sở (tần số 0.05Hz đến dưới 0.8Hz sinh ra do hô hấp) làm trôi đường biên, gây khó khăn cho việc nhận dạng đỉnh QRS và phân tích đoạn S-T, dẫn đến nguy cơ chẩn đoán nhầm các bệnh lý như nhồi máu cơ tim hay cơn đau thắt ngực.

2. Thuật toán lọc nhiễu 50Hz/60Hz đề xuất có điểm gì khác biệt so với các nghiên cứu trước đây?

Thuật toán đề xuất có hai điểm cải tiến cốt lõi: (1) Sử dụng biến đổi Fourier nhanh (FFT) để dò tìm tự động và đơn giản tần số nhiễu trong dải hẹp (47-53Hz cho lưới điện 50Hz; 57-63Hz cho lưới điện 60Hz); (2) Thiết lập ngưỡng biên độ dư sau lọc (15 $\mu$V cho 50Hz; 20 $\mu$V cho 60Hz) để điều khiển vòng lặp lọc nhiều lần. Khi biên độ nhiễu ban đầu lớn ($\ge 0.4$mV), thuật toán tự động tăng số vòng lặp (2-3 vòng), giúp duy trì sai số MSE thấp ($40 \times 10^{-4}$ đến $48 \times 10^{-4}$) trong khi phương pháp trong tài liệu [2] bị suy giảm hiệu quả nghiêm trọng (MSE tăng lên $49 \times 10^{-3}$ và $77 \times 10^{-3}$).

3. Phương pháp xác định các điểm đẳng thế trong Chương 2 dựa trên những tính chất toán học nào?

Phương pháp dựa trên 3 tính chất toán học:

  • $QRS_{onset}$ và $QRS_{offset}$ là các điểm uốn, tương ứng với điểm vượt 0 của đạo hàm bậc 2 ($QRS_{onset}$ vượt 0 từ âm lên dương; $QRS_{offset}$ vượt 0 từ dương xuống âm).
  • $T_{offset}$ là điểm cực tiểu sau sóng T (điểm vượt 0 từ âm lên dương của đạo hàm bậc 1); $P_{onset}$ là điểm vượt 0 từ dương của đạo hàm bậc 1.
  • $T_{offset}$ là điểm vượt 0 từ âm lên dương của đạo hàm bậc 2 (sau đỉnh T); $P_{onset}$ là điểm vượt 0 từ dương của đạo hàm bậc 2 (trước đỉnh P).

4. Kết quả lọc nhiễu trôi dạt đường cơ sở của luận án vượt trội hơn công trình của Na Pan và cộng sự [54] như thế nào?

Khi thử nghiệm trên cùng bản ghi chuẩn 301 từ cơ sở dữ liệu MIT-BIH ST Change Database, phương pháp đề xuất dựa trên điểm đẳng thế và nội suy Cubic Spline đạt hệ số tương quan trung bình $r = 0.9984$, cao hơn đáng kể so với phương pháp phân tách theo kinh nghiệm (EMD) của Na Pan ($r = 0.9746$). Khi thử nghiệm trên các bản ghi 301 và 302, phương pháp đạt sai số MSE rất nhỏ (0.0018 $\text{mV}^2$) và SNR đạt 26.54 dB.

Kết luận

Luận án đã giải quyết bài toán truyền dẫn tín hiệu điện tim thời gian thực phục vụ theo dõi và hỗ trợ chẩn đoán bệnh nhân tim mạch từ xa qua công nghệ vô tuyến WiFi. Tác giả đã đề xuất thành công hai giải pháp tiền xử lý: bộ lọc triệt tần thích nghi FFT nhiều vòng lặp loại bỏ nhiễu nguồn 50Hz/60Hz và giải pháp nội suy Cubic Spline qua các điểm đẳng thế trong miền thời gian để khử nhiễu trôi dạt đường cơ sở. Cùng với phương pháp nén dữ liệu hai trạng thái, công trình đã hình thành gói giải pháp kỹ thuật hoàn chỉnh, giúp nâng cao chất lượng tín hiệu ECG thu nhận, giảm tỷ lệ mất gói và tối ưu hóa tài nguyên truyền dẫn cũng như lưu trữ y tế.