Tổng quan nghiên cứu

Theo thống kê của Tổ chức Y tế Thế giới, các 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, cướp đi sinh mạng của khoảng 17,5 triệu người mỗi năm. Trung bình cứ 2 giây có một trường hợp tử vong do bệnh tim, 5 giây có một ca nhồi máu cơ tim và 6 giây lại xuất hiện một ca đột quỵ. Tín hiệu điện tâm đồ giữ vai trò là phương pháp chẩn đoán lâm sàng không xâm lấn, chi phí thấp và hiệu quả cao để phát hiện sớm các rối loạn dẫn truyền tim.

Tuy nhiên, trong quá trình ghi nhận thực tế, biên độ của sóng điện tim rất nhỏ, đỉnh sóng R chỉ dao động từ 1,5 mV đến 2 mV và sóng P khoảng 0,25 mV. Tín hiệu này nằm trong dải tần số từ 0,05 Hz đến 100 Hz, rất dễ bị biến dạng bởi nhiều nguồn nhiễu như nhiễu nguồn điện xoay chiều 50 Hz đến 60 Hz, nhiễu điện cơ từ 20 Hz đến 30 Hz và nhiễu do tiếp xúc điện cực.

Mục tiêu cốt lõi của nghiên cứu là xây dựng giải pháp khử nhiễu tối ưu cho tín hiệu điện tim bằng cách ứng dụng phép biến đổi Wavelet rời rạc, khắc phục hoàn toàn nhược điểm mất thông tin thời gian của phép biến đổi Fourier và hạn chế méo tín hiệu của thuật toán lọc thích nghi. Đề tài được triển khai thực hiện tại Khoa Vật lý thuộc Trường Đại học Sư phạm Thành phố Hồ Chí Minh, hoàn thành vào năm 2015 dưới sự tư vấn chuyên môn từ các bác sĩ tim mạch Bệnh viện Bình Chánh. Kết quả nghiên cứu cung cấp thuật toán xử lý tín hiệu đạt tỷ số tín hiệu trên nhiễu cao, hạ thấp sai số tương đối dưới 5% và tích hợp thành công trên giao diện phần mềm trực quan.

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 trên sự kết hợp chặt chẽ giữa cơ sở điện sinh lý học tim mạch và lý thuyết xử lý tín hiệu số y sinh hiện đại:

  • Hoạt động điện thế tim: Tim hoạt động theo chu chuyển 0,8 giây gồm ba giai đoạn: nhĩ thu 0,1 giây, thất thu 0,3 giây và tâm trương toàn bộ 0,4 giây. Sự biến thiên điện thế màng tế bào từ trạng thái nghỉ -90 mV sang trạng thái khử cực 0 mV và tái cực tạo nên các dạng sóng P, Q, R, S, T. Trục điện tim trung bình hợp với phương ngang một góc khoảng 85 độ.
  • Phép biến đổi Wavelet rời rạc: Phát triển từ thuật toán lọc đa kênh của nhà khoa học Stéphane Mallat năm 1988. Kỹ thuật này phân tách tín hiệu thành các thành phần xấp xỉ tần số thấp chứa năng lượng chính và các thành phần chi tiết tần số cao chứa thông tin nhiễu. Khác với biến đổi Fourier thời gian ngắn bị giới hạn bởi kích thước cửa sổ cố định, Wavelet giải quyết triệt để nguyên lý bất định Heisenberg nhờ khả năng co giãn tỷ lệ và dịch chuyển thời gian linh hoạt.
  • Lý thuyết lọc ngưỡng Donoho và Johnstone: Ứng dụng quy tắc đặt ngưỡng cứng và ngưỡng mềm để loại bỏ các hệ số chi tiết mang năng lượng nhiễu, tái tạo lại tín hiệu sạch nguyên bản thông qua phép biến đổi ngược.

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

  • Nguồn dữ liệu và mô hình thực nghiệm: Tín hiệu điện tim chuẩn được đưa vào mô hình thử nghiệm pha tạp nhiễu Gauss trắng với phương sai chuẩn hóa, khảo sát qua các mức tỷ lệ tín hiệu trên nhiễu từ 0 dB đến 18 dB, tập trung phân tích sâu tại 5 dB, 10 dB, 15 dB và 18 dB.
  • Cỡ mẫu và phương pháp chọn mẫu: Nghiên cứu tiến hành khảo sát thực nghiệm toàn diện trên 54 hàm Wavelet mẹ thuộc các họ trực giao và song trực giao tiêu biểu như Daubechies, Symlets, Coiflets và Biorthogonal. Quá trình phân tích thực hiện trên 4 mức phân tách điển hình gồm mức 2, mức 4, mức 6 và mức 8.
  • Phương pháp phân tích và lý do lựa chọn: Sử dụng thuật toán phân tách đa độ phân giải kết hợp các hàm xử lý ngưỡng tối ưu trong phần mềm MATLAB. Phương pháp này được lựa chọn vì tính chất tín hiệu điện tim là tín hiệu không dừng, đòi hỏi độ phân giải thời gian cực cao ở vùng tần số cao và độ phân giải tần số tốt ở vùng tần số thấp. Toàn bộ quá trình tính toán và kiểm chứng được tiến hành xuyên suốt trong giai đoạn 2014 đến 2015.

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

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

  • Phát hiện 1: Trong số 54 hàm Wavelet mẹ được khảo sát, họ Daubechies (đặc biệt là dạng db10 và db8) cùng họ Symlets mang lại hiệu suất lọc vượt trội nhất. Tỷ lệ tín hiệu trên nhiễu sau lọc tăng từ 30% đến 45% đối với các tín hiệu bị suy biến nặng ở mức 5 dB ban đầu.
  • Phát hiện 2: Mức phân tách tối ưu cho tín hiệu điện tim được xác định chuẩn xác ở mức 4 đến mức 6. Tại mức phân tách 4, chỉ số sai số trung bình bình phương giảm mạnh về mức tiệm cận tối thiểu, đồng thời sai khác căn bậc hai trung bình duy trì ổn định dưới 4,5%.
  • Phát hiện 3: Phương pháp đặt ngưỡng mềm thể hiện tính ưu việt tuyệt đối so với ngưỡng cứng. Ngưỡng mềm loại bỏ hoàn toàn các gai nhọn gián đoạn, làm mịn đường đẳng điện và nâng cao chỉ số SNR thêm 1,8 dB đến 2,5 dB ở cùng mức nhiễu đầu vào 10 dB.
  • Phát hiện 4: Xây dựng thành công giao diện tương tác người dùng trên nền tảng MATLAB GUI, cho phép nạp tín hiệu, tùy chọn hàm lọc, thiết lập tham số ngưỡng và hiển thị kết quả khử nhiễu tức thì với thời gian xử lý chưa đến 0,5 giây.

Thảo luận kết quả

Hiệu quả vượt trội của hàm Wavelet db10 xuất phát từ việc dạng sóng của hàm mẹ này có sự tương đồng hình học rất cao với phức bộ QRS tự nhiên của cơ tim người, nơi tập trung phần lớn năng lượng khử cực ở biên độ 1,5 mV đến 2 mV.

Trên các đồ thị phân bố dữ liệu biểu diễn mối tương quan giữa chỉ số SNR và MSE, đường cong đáp ứng cho thấy khi tăng số mức phân tách từ 2 lên 4, chỉ số SNR tăng vọt theo phương thẳng đứng. Tuy nhiên, nếu tiếp tục nâng mức phân tách lên mức 8, chỉ số sai số PRD bắt đầu tăng trở lại do hiện tượng làm mịn quá mức, gây suy hao biên độ của sóng P vốn chỉ có độ lớn khoảng 0,25 mV.

So với các công trình nghiên cứu sử dụng biến đổi Fourier nhanh hoặc mạng nơ-ron nhân tạo của các nhà khoa học quốc tế trong giai đoạn 2010 đến 2011, giải pháp Wavelet rời rạc đa mức kết hợp ngưỡng mềm giúp bảo toàn nguyên vẹn pha của tín hiệu và không làm trễ đỉnh sóng R, mang lại giá trị chẩn đoán lâm sàng trung thực.

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

  • Tích hợp thuật toán lọc Wavelet db10 mức 4 vào bộ vi xử lý của các máy đo điện tim cầm tay và thiết bị theo dõi bệnh nhân liên tục, hướng tới mục tiêu giảm tỷ lệ cảnh báo sai lệch lâm sàng xuống dưới 2% trong vòng 12 tháng, do các kỹ sư y sinh và doanh nghiệp sản xuất thiết bị y tế phối hợp thực hiện.
  • Chuẩn hóa quy trình xử lý số tín hiệu điện tâm đồ tại các phòng khám và bệnh viện đa khoa, loại bỏ triệt để can nhiễu từ nguồn điện lưới 50 Hz đến 60 Hz và duy trì độ sai lệch biên độ sóng dưới 3%, triển khai thực hiện trong lộ trình 6 tháng bởi đội ngũ kỹ thuật viên y tế.
  • Nâng cấp mô hình thuật toán hướng đến xử lý thời gian thực cho hệ thống điện tim 12 đạo trình đồng thời, tối ưu hóa tốc độ xử lý đạt độ trễ dưới 100 mili-giây trong giai đoạn 2024 đến 2025, do các nhóm nghiên cứu công nghệ y sinh tại các trường đại học đảm nhiệm.
  • Biên soạn tài liệu và chuyển giao công nghệ lập trình xử lý tín hiệu y sinh trên phần mềm MATLAB cho sinh viên và học viên cao học, hoàn thành trong vòng 1 học kỳ bởi Khoa Vật lý và các đơn vị đào tạo kỹ thuật liên quan.

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

  • Bác sĩ chuyên khoa tim mạch và nhân viên y tế: Giúp nắm bắt sâu sắc bản chất suy giảm chất lượng tín hiệu điện thế màng tế bào, từ đó phân tích chính xác các dạng sóng P-QRS-T đã được lọc sạch để chẩn đoán rối loạn nhịp tim chuẩn xác.
  • Kỹ sư thiết kế thiết bị y tế: Cung cấp sơ đồ khối thuật toán lọc sóng đa kênh, bộ thông số ngưỡng tối ưu và mã nguồn mẫu để nhúng trực tiếp vào các thiết bị đo điện tim thế hệ mới.
  • Giảng viên, nhà nghiên cứu và học viên cao học ngành Vật lý y sinh: Đóng vai trò là tài liệu tham khảo học thuật giá trị về phân tích miền thời gian - tần số, làm nền tảng mở rộng nghiên cứu cho các tín hiệu điện sinh học khác như điện não đồ hoặc điện cơ đồ.
  • Sinh viên ngành Kỹ thuật y sinh và Công nghệ thông tin: Hướng dẫn chi tiết phương pháp lập trình giao diện đồ họa người dùng và cách áp dụng các chỉ số SNR, MSE, PRD để đánh giá khoa học một mô hình xử lý dữ liệu.

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

  • Tại sao biến đổi Wavelet lại phù hợp hơn biến đổi Fourier trong xử lý tín hiệu điện tim? Biến đổi Fourier chỉ cung cấp thông tin tần số thuần túy mà làm mất hoàn toàn tọa độ thời gian của biến cố. Tín hiệu điện tim có tính chất không dừng, chứa các sóng nhọn xuất hiện theo từng thời điểm. Biến đổi Wavelet giải quyết triệt để vấn đề này nhờ cơ chế cửa sổ co giãn linh hoạt, vừa định vị chính xác vị trí thời gian của đỉnh sóng R, vừa giữ độ phân giải tần số ổn định trong dải 0,05 Hz đến 100 Hz.

  • Hàm Wavelet mẹ nào đạt hiệu quả khử nhiễu điện tâm đồ cao nhất? Qua quá trình thực nghiệm đối sánh 54 hàm Wavelet khác nhau, hàm db10 thuộc họ Daubechies cho kết quả tốt nhất. Dạng sóng mẹ của db10 có cấu trúc tương thích cao với phức bộ sóng QRS, giúp chỉ số SNR đạt mức cao nhất và giữ chỉ số sai lệch PRD luôn ở mức an toàn dưới 4,5%.

  • Số mức phân tách Wavelet tối ưu là bao nhiêu trong thực tế? Mức phân tách lý tưởng nhất nằm trong khoảng từ mức 4 đến mức 6. Nếu chọn mức quá thấp như mức 2, các thành phần nhiễu tần số cao vẫn chưa được bóc tách triệt để. Ngược lại, nếu chọn mức quá cao như mức 8, các chi tiết giải phẫu quan trọng như sóng P biên độ 0,25 mV sẽ bị suy giảm năng lượng dẫn đến méo dạng sóng.

  • Ngưỡng mềm có ưu điểm gì nổi bật so với ngưỡng cứng? Ngưỡng cứng triệt tiêu các giá trị nhỏ hơn ngưỡng nhưng giữ nguyên giá trị lớn hơn, tạo nên các điểm đứt gãy đột ngột trên đồ thị. Ngưỡng mềm thực hiện co các giá trị lớn về gần 0, giúp đường cong tín hiệu sau tái tạo đạt độ trơn mịn tự nhiên và cải thiện chỉ số SNR cao hơn từ 1,8 dB đến 2,5 dB so với ngưỡng cứng.

  • Nhiễu Gauss trắng mô phỏng những nguồn can nhiễu nào trên thực tế? Nhiễu Gauss trắng với dải tỷ số tín hiệu trên nhiễu từ 0 dB đến 18 dB được sử dụng làm mô hình chuẩn hóa cho hầu hết các loại nhiễu cộng tính trong môi trường lâm sàng, bao gồm nhiễu điện cơ 20 Hz đến 30 Hz do bệnh nhân run cơ, nhiễu dịch chuyển đường đẳng điện và can nhiễu từ điện cực tiếp xúc trên bề mặt da.

Kết luận

  • Đã chứng minh toàn diện tính ưu việt của phép biến đổi Wavelet rời rạc trong việc loại bỏ can nhiễu và phục hồi trung thực tín hiệu điện tâm đồ.
  • Xác lập thành công bộ thông số tối ưu gồm hàm mẹ Daubechies db10, mức phân tách 4 và phương pháp đặt ngưỡng mềm thích nghi.
  • Đạt được các tiêu chuẩn định lượng khắt khe với tỷ số tín hiệu trên nhiễu tăng cao, sai số tương đối PRD kiểm soát dưới mức 5%.
  • Hoàn thiện giao diện phần mềm MATLAB GUI trực quan, thân thiện, đáp ứng tốt cho cả mục đích nghiên cứu học thuật lẫn thực nghiệm tiền lâm sàng.
  • Đóng góp cơ sở lý thuyết và giải thuật thực nghiệm vững chắc cho chuyên ngành Vật lý kỹ thuật y sinh tại Việt Nam.

Công trình đã giải quyết trọn vẹn bài toán nâng cao chất lượng tín hiệu chẩn đoán tim mạch không xâm lấn. Trong 2 năm tới, hướng phát triển trọng tâm là mở rộng thuật toán sang xử lý tín hiệu trực tuyến đa kênh trên vi điều khiển nhúng. Hãy tham khảo và ứng dụng ngay giải pháp xử lý tín hiệu Wavelet chuyên sâu này để tối ưu hóa độ chính xác cho các hệ thống theo dõi sức khỏe tim mạch hiện đại.