Tách nguồn mù và trích đặc trưng bằng giải thuật FastICA tại HCMUTE

Chuyên khảo phân tích Hcmute ứng dụng giải thuật fastica trong tách nguồn mù và trích đặc trưng, đánh giá các khía cạnh quan trọng, đề xuất hướng nghiên cứu tiếp theo.

• Tác giả: Ngô Quốc Cường
• Trường học: Trường Đại Học Sư Phạm Kỹ Thuật Thành Phố Hồ Chí Minh
• Chuyên ngành: Khoa Điện – Điện Tử
• Đề tài: Ứng Dụng Giải Thuật FastICA Trong Tách Nguồn Mù Và Trích Đặc Trưng
• Loại tài liệu: báo cáo tổng kết đề tài KH&CN cấp trường
• Năm xuất bản: 2013
• Địa điểm: Thành Phố Hồ Chí Minh

Mục lục chi tiết

• MỞ ĐẦU
• 1. CHƯƠNG 1: CƠ SỞ LÝ THUYẾT
• 1.1. Lý Thuyết Chung Về ICA
• 1.2. Giả thiết trong ICA
• 1.3. Ước lượng Trong ICA
• 1.4. Tiền Xử Lý trong ICA
• 1.5. ICA dùng phương pháp cực đại tính phi Gauss
• 1.6. Phương Pháp Tái Hiện Hình Ảnh Và Tín Hiệu Não Người
• 1.6.1. Phương pháp EEG
• 1.6.2. Phương pháp MRI
• 1.6.3. Phương pháp Quang Phổ Cận Hồng Ngoại fNIRS
• 2. CHƯƠNG 2: ỨNG DỤNG FAST-ICA TRONG TÁCH NGUỒN MÙ
• 2.1. Nguồn Tín Hiệu Sin và Răng Cưa
• 2.2. Nguồn Tín Hiệu Âm Thanh
• 3. ỨNG DỤNG FASTICA TRONG TRÍCH ĐẶC TRƯNG TÍN HIỆU NIRS
• 3.1. Các Vùng Chức Năng Não Người
• 3.2. Thiết Bị và Tín Hiệu fNIRS
• 3.3. Thiết Lập Thí Nghiệm Thu Dữ Liệu Cho Hoạt Động Gõ Tay
• 3.4. Tiền xử lý dữ liệu
• 3.5. Mô hình trích đặc trưng áp dụng fastICA
• 4. CHƯƠNG 4
• 4.1. Hướng phát triển
• DANH SÁCH HÌNH
• DANH SÁCH TỪ VIẾT TẮT
• THÔNG TIN KẾT QUẢ NGHIÊN CỨU

Tóm tắt

I. Giới thiệu về giải thuật FastICA và ứng dụng trong xử lý tín hiệu

Bài báo nghiên cứu ứng dụng giải thuật FastICA trong tách nguồn mù và trích đặc trưng, tập trung vào phân tích thành phần độc lập (ICA). FastICA là một thuật toán nhanh chóng và hiệu quả để thực hiện ICA, một kỹ thuật xử lý tín hiệu nhằm tách các nguồn tín hiệu độc lập từ hỗn hợp tín hiệu quan sát được. Giải thuật tách nguồn mù đóng vai trò quan trọng trong nhiều ứng dụng, đặc biệt là khi thông tin về nguồn tín hiệu không đầy đủ hoặc chưa biết trước. Ứng dụng FastICA trong xử lý tín hiệu mang lại nhiều lợi ích, bao gồm khả năng tách biệt các nguồn tín hiệu chồng chéo, cải thiện chất lượng tín hiệu và trích xuất thông tin hữu ích.

1.1 Cơ sở lý thuyết ICA

Mô hình ICA cơ bản giả định rằng tín hiệu quan sát được là sự kết hợp tuyến tính của các nguồn tín hiệu độc lập. Mục tiêu của ICA là ước lượng các nguồn tín hiệu độc lập này và ma trận hỗn hợp. Các giả định chính của ICA bao gồm tính độc lập thống kê giữa các nguồn, tính phi Gauss của các nguồn, và ma trận hỗn hợp là ma trận vuông. Các phương pháp ước lượng trong ICA thường tập trung vào việc tối đa hóa tính phi Gauss của các nguồn, ví dụ sử dụng kurtosis hoặc negentropy. Giải thuật ICA được áp dụng rộng rãi trong nhiều lĩnh vực, từ xử lý âm thanh, hình ảnh đến xử lý tín hiệu sinh học. Thuật toán FastICA là một trong những thuật toán ICA hiệu quả nhất, nổi bật với tốc độ hội tụ nhanh và độ tin cậy cao. Giải thuật tách nguồn mù dựa trên ICA đã được chứng minh hiệu quả trong nhiều bài toán thực tế.

1.2 Thuật toán FastICA và các bước thực hiện

Thuật toán FastICA dựa trên nguyên lý cực đại hóa tính phi Gauss của các nguồn tín hiệu. Thuật toán sử dụng phương pháp lặp điểm cố định để tìm kiếm các vectơ trọng số tối ưu. Các bước thực hiện FastICA bao gồm: tiền xử lý dữ liệu (quy tâm và trắng hóa), cập nhật lặp vectơ trọng số, và chuẩn hóa kết quả. FastICA sử dụng kurtosis làm hàm đo lường tính phi Gauss. Sự lựa chọn hàm đo lường này ảnh hưởng đến hiệu quả của thuật toán. Việc lựa chọn tham số thích hợp trong FastICA cũng rất quan trọng để đảm bảo hiệu quả và độ chính xác của kết quả. Ứng dụng FastICA trong các hệ thống thực tế đòi hỏi việc hiểu rõ các tham số và tối ưu hóa thuật toán cho từng ứng dụng cụ thể. Hiệu quả của FastICA phụ thuộc vào tính chất của dữ liệu đầu vào và các giả định của mô hình ICA.

II. Ứng dụng FastICA trong tách nguồn mù tại HCMUTE

Nghiên cứu ứng dụng FastICA tại Đại học Công nghệ Thành phố Hồ Chí Minh (HCMUTE) tập trung vào tách nguồn mù trong các tín hiệu khác nhau, bao gồm tín hiệu âm thanh và tín hiệu sinh học. Kết quả thực nghiệm cho thấy FastICA có khả năng tách biệt hiệu quả các nguồn tín hiệu gốc từ hỗn hợp tín hiệu quan sát được, ngay cả trong trường hợp có nhiễu. Việc so sánh các giải thuật tách nguồn mù khác cũng được thực hiện để đánh giá hiệu quả của FastICA. Nghiên cứu này đã chứng minh tính khả thi và hiệu quả của FastICA trong tách tín hiệu.

2.1 Tách nguồn mù trong tín hiệu âm thanh

Nghiên cứu đã áp dụng FastICA để tách nguồn mù trong các tín hiệu âm thanh được trộn lẫn. Kết quả cho thấy FastICA có khả năng tách biệt các nguồn âm thanh gốc với độ chính xác cao, ngay cả khi các nguồn tín hiệu có tần số gần nhau. Hiệu quả của FastICA được đánh giá thông qua các chỉ số chất lượng tín hiệu như tỷ số tín hiệu trên nhiễu (SNR) và độ méo hài hòa (THD). Các thí nghiệm được thực hiện bằng cách sử dụng các nguồn tín hiệu mô phỏng và các nguồn tín hiệu thu thập từ thực tế. Việc so sánh các giải thuật tách nguồn mù khác như PCA cho thấy FastICA mang lại hiệu quả tốt hơn trong việc tách các nguồn tín hiệu độc lập, đặc biệt là trong trường hợp tín hiệu có độ chồng chéo cao. Phân tích tần số được sử dụng để đánh giá chất lượng tín hiệu sau khi tách. Python FastICA và Matlab FastICA là hai công cụ phổ biến được sử dụng trong nghiên cứu này.

2.2 Tách nguồn mù trong tín hiệu sinh học NIRS

Nghiên cứu cũng ứng dụng FastICA trong phân tích dữ liệu đa chiều thu được từ thiết bị quang phổ cận hồng ngoại (fNIRS). FastICA được sử dụng để trích xuất đặc trưng từ tín hiệu NIRS nhằm nhận dạng các hoạt động của não. Trích xuất đặc trưng hình ảnh và trích xuất đặc trưng âm thanh là hai hướng ứng dụng quan trọng. Trích xuất đặc trưng NIRS giúp phân tích hoạt động não bộ liên quan đến các nhiệm vụ nhận thức như gõ tay. Mô hình trích đặc trưng dựa trên FastICA cho thấy khả năng phân biệt giữa các hoạt động khác nhau. Tiền xử lý dữ liệu như lọc nhiễu và chuẩn hóa là những bước cần thiết để đảm bảo độ chính xác của kết quả. Thực nghiệm được thiết kế để thu thập dữ liệu từ các đối tượng tham gia và phân tích kết quả. Kết quả thực nghiệm cho thấy FastICA hiệu quả trong giảm chiều dữ liệu và trích xuất đặc trưng từ tín hiệu NIRS, mở ra hướng nghiên cứu mới trong lĩnh vực xử lý tín hiệu sinh học. Hạn chế của FastICA trong ứng dụng này cần được nghiên cứu thêm.

III. Kết luận và hướng phát triển

Nghiên cứu đã chứng minh hiệu quả của giải thuật FastICA trong tách nguồn mù và trích đặc trưng trên nhiều loại tín hiệu. FastICA là một công cụ mạnh mẽ cho xử lý tín hiệu và xử lý dữ liệu. Hiệu quả của FastICA phụ thuộc vào việc lựa chọn các tham số và tiền xử lý dữ liệu. Các thư viện FastICA cung cấp hỗ trợ mạnh mẽ cho việc triển khai thuật toán. Nghiên cứu mở ra nhiều hướng phát triển, bao gồm cải tiến thuật toán, ứng dụng trong các lĩnh vực khác và tích hợp với các kỹ thuật xử lý tín hiệu khác.

3.1 Đánh giá tổng quan về hiệu quả của FastICA

Nghiên cứu cho thấy FastICA đạt hiệu quả cao trong việc tách nguồn mù và trích xuất đặc trưng, đặc biệt trong các trường hợp tín hiệu có nhiễu và chồng chéo. Hiệu quả của FastICA được thể hiện qua kết quả thực nghiệm với độ chính xác cao. Tuy nhiên, hạn chế của FastICA cũng cần được xem xét, ví dụ như độ nhạy cảm với các tham số đầu vào và yêu cầu về tính phi Gauss của nguồn tín hiệu. Cải tiến FastICA là một hướng nghiên cứu quan trọng để khắc phục các hạn chế này. Đánh giá hiệu năng của FastICA nên được tiến hành trên nhiều tập dữ liệu khác nhau để đảm bảo tính tổng quát của kết quả. So sánh FastICA với các thuật toán khác cũng là một cách để đánh giá hiệu quả của nó. Mô hình toán học của FastICA được xây dựng trên cơ sở lý thuyết vững chắc của ICA.

3.2 Hướng phát triển trong tương lai

Các hướng phát triển trong tương lai bao gồm: cải tiến FastICA để xử lý dữ liệu có nhiễu cao hơn và độ chồng chéo lớn hơn; ứng dụng FastICA trong các lĩnh vực mới như xử lý video, phân tích dữ liệu lớn; tích hợp FastICA với các kỹ thuật học máy khác để nâng cao hiệu quả; nghiên cứu hạn chế của FastICA và tìm kiếm các giải pháp khắc phục; phát triển các thư viện FastICA mở rộng hỗ trợ nhiều loại dữ liệu và nền tảng; nghiên cứu khoa học về FastICA cần được tiếp tục để thúc đẩy sự phát triển của công nghệ xử lý tín hiệu. Bài báo khoa học về ứng dụng FastICA có giá trị thực tiễn cao và đóng góp quan trọng vào lĩnh vực xử lý tín hiệu.

Nội dung chính

ỨNG DỤNG GIẢI THUẬT FASTICA TRONG TÁCH NGUỒN MÙ VÀ TRÍCH ĐẶC TRƯNG TÍN HIỆU NÃO fNIRS

---

Tóm tắt nghiên cứu (200-250 từ)

Nghiên cứu khoa học "Ứng dụng giải thuật FastICA trong tách nguồn mù và trích đặc trưng" (Mã số đề tài: T2013-13), do ThS. Ngô Quốc Cường chủ nhiệm cùng sự tham gia của ThS. Nguyễn Ngô Lâm tại Trường Đại học Sư phạm Kỹ thuật TP.HCM, tập trung giải quyết bài toán phân tách nguồn mù (Blind Source Separation - BSS) và trích xuất đặc trưng sinh học thần kinh phức tạp.

Câu hỏi nghiên cứu cốt lõi đặt ra là: Làm thế nào để ứng dụng hiệu quả giải thuật phân tích thành phần độc lập tốc độ cao (FastICA) nhằm khôi phục các tín hiệu gốc từ hỗn hợp trộn chưa biết thông tin kênh truyền, đồng thời trích xuất đặc trưng vận động từ tín hiệu quang phổ cận hồng ngoại chức năng (fNIRS)?

Phương pháp luận của đề tài kết hợp giữa mô phỏng toán học trên các dạng sóng tổng hợp (sin, răng cưa, xung vuông, tín hiệu âm thanh) và thực nghiệm lâm sàng thu thập dữ liệu huyết động học não bộ bằng hệ thống FOIRE-3000 (Shimadzu) với giao thức vận động gõ ngón tay.

Kết quả nổi bật cho thấy giải thuật FastICA dựa trên tối đa hóa tính phi Gauss qua đại lượng Kurtosis có khả năng tách chính xác các hài tần số gốc ngay cả trong điều kiện có nhiễu (SNR = 16 dB). Đặc biệt, khi mô hình hóa bài toán fNIRS thành sơ đồ tương đương "tiệc cocktail", FastICA đã trích xuất thành công đặc trưng phổ trong dải tần số 6 – 7 Hz của tín hiệu Oxy-Hb, giúp phân biệt rõ rệt giữa vận động tay trái (đạt độ chính xác nhận định 27/32 mẫu) và tay phải (17/32 mẫu). Nghiên cứu mở ra hướng tiếp cận triển vọng cho phát triển hệ thống giao diện não - máy tính (BCI).

---

Bối cảnh và tầm quan trọng (300-350 từ)

Trong kỷ nguyên xử lý tín hiệu số hiện đại, bài toán tách nguồn mù (BSS) – điển hình là hiện tượng "tiệc cocktail" (Cocktail Party Problem) – luôn là một thách thức lớn. Vấn đề cốt lõi nằm ở việc khôi phục các nguồn tín hiệu thành phần ban đầu từ các kênh cảm biến quan sát được mà hoàn toàn không biết trước thông tin về ma trận trộn lẫn cũng như đặc tính vật lý của môi trường truyền dẫn. Mặc dù các kỹ thuật cổ điển như PCA hay lọc thích nghi tuyến tính mang lại một số kết quả nhất định, chúng chủ yếu dựa vào các thống kê bậc hai (trung bình, hiệp phương sai), vốn không đủ khả năng phân ly các biến ngẫu nhiên phi Gauss độc lập.

Sự xuất hiện của lý thuyết Phân tích thành phần độc lập (ICA - Independent Component Analysis) và giải thuật điểm cố định FastICA đã tạo ra bước đột phá nhờ khai thác các thống kê bậc cao (Kurtosis, Negentropy). Tuy nhiên, việc đưa FastICA vào xử lý tín hiệu y sinh phi tuyến, đặc biệt là tín hiệu não bộ, vẫn tồn tại nhiều khoảng trống nghiên cứu.

Tại Việt Nam vào thời điểm nghiên cứu (năm 2013), phần lớn các nghiên cứu tín hiệu não tập trung vào Điện não đồ (EEG) hoặc Cộng hưởng từ chức năng (fMRI). Kỹ thuật quang phổ cận hồng ngoại chức năng (fNIRS - functional Near-Infrared Spectroscopy) là một hướng tiếp cận hoàn toàn mới, sở hữu ưu thế vượt trội: không xâm lấn, chi phí hợp lý, độ phân giải không gian tốt hơn EEG và tính cơ động cao hơn fMRI. Dù vậy, tín hiệu fNIRS rất dễ bị suy biến bởi nhiễu sinh lý (nhịp tim, huyết áp, hô hấp) và nhiễu do dịch chuyển đầu đo.

Đề tài mang tính cấp thiết cao khi kết hợp giải thuật xử lý tín hiệu tiên tiến FastICA với công nghệ thu nhận quang học fNIRS để giải mã tín hiệu vùng vỏ não vận động (motor cortex). Nghiên cứu không chỉ đặt nền móng lý thuyết và thực nghiệm cho chuyên ngành Điện tử - Viễn thông và Kỹ thuật Y sinh trong nước, mà còn đóng góp trực tiếp vào xu hướng phát triển công nghệ tương tác người - máy hiện đại.

---

Methodology và approach (350-400 từ)

Nghiên cứu được thiết kế theo quy trình khoa học chặt chẽ, kết hợp tuần tự từ cơ sở lý thuyết, mô phỏng toán học chuyên sâu đến thực nghiệm đo đạc sinh học:

[Mô hình toán ICA & FastICA] 
         │
         ▼
[Tiền xử lý: Quy tâm + Trắng hóa (EVD)]
         │
         ▼
[Kiểm thử mô phỏng: Sóng tuần hoàn & Âm thanh]
         │
         ▼
[Thực nghiệm fNIRS: Hệ thống FOIRE-3000]
         │
         ▼
[Lọc Savitzky-Golay (N=11, Bậc 3)]
         │
         ▼
[Mô hình hóa Cocktail Party & Tách FastICA]
         │
         ▼
[Phân tích phổ & Đánh giá phân biệt vận động]

1. Thiết kế mô hình toán học và giải thuật

Mô hình toán học xuất phát từ phương trình tuyến tính vector: $$\mathbf{x} = \mathbf{A}\mathbf{s}$$ Trong đó:

• $\mathbf{x}$ là vector tín hiệu trộn quan sát được.
• $\mathbf{A}$ là ma trận trộn chưa biết.
• $\mathbf{s}$ là vector các nguồn độc lập thống kê phi Gauss.

Thuật toán thực thi qua các bước:

1. Tiền xử lý dữ liệu:
   • Quy tâm (Centering): Đưa kỳ vọng toán học của tín hiệu về 0: $$\mathbf{x} \leftarrow \mathbf{x} - E{\mathbf{x}}$$
   • Trắng hóa (Whitening): Sử dụng phương pháp phân giải trị riêng (EVD - Eigenvalue Decomposition) của ma trận hiệp phương sai để loại bỏ tương quan bậc hai: $$\mathbf{z} = \mathbf{V}\mathbf{x} = \mathbf{E}\mathbf{D}^{-1/2}\mathbf{E}^T\mathbf{x}$$
2. Ước lượng FastICA: Tối đa hóa giá trị tuyệt đối của Kurtosis bậc bốn ($kurt(y) = E{y^4} - 3(E{y^2})^2$) thông qua quy tắc lặp điểm cố định có nhân tử Lagrange: $$\mathbf{w} \leftarrow E{\mathbf{z}(\mathbf{w}^T\mathbf{z})^3} - 3\mathbf{w}$$ Sau mỗi vòng lặp, vector trọng số $\mathbf{w}$ được chuẩn hóa ($\mathbf{w} \leftarrow \mathbf{w} / |\mathbf{w}|$) cho đến khi đạt điều kiện hội tụ.

2. Thiết lập thực nghiệm và thu thập dữ liệu

• Thiết bị đo: Hệ thống FOIRE-3000 (Shimadzu, Nhật Bản) với 3 bước sóng cận hồng ngoại: $780\text{ nm}$, $805\text{ nm}$, $830\text{ nm}$, tính toán nồng độ Hemoglobin dựa trên định luật Beer-Lambert mở rộng: $$\Delta \text{OD}(\lambda) = [\varepsilon_{\text{Hb}}(\lambda)[\text{HbR}] + \varepsilon_{\text{HbO}_2}(\lambda)[\text{HbO}_2]] \cdot L \cdot \text{DPF}$$
• Đối tượng & Bố trí: Chủ thể nam 25 tuổi, thuận tay phải. Bố trí 20 kênh đo (10 kênh bán cầu não trái, 10 kênh não phải) phủ lên vùng vỏ não vận động (motor control - thùy đỉnh). Khoảng cách giữa các đầu phát (optode emitter) và đầu thu (detector) cố định ở mức $3\text{ cm}$.
• Giao thức kích thích (Paradigm): Chu kỳ 60 giây gồm $20\text{ s}$ nghỉ (Rest) $\rightarrow$ $20\text{ s}$ thực hiện gõ tay liên tục khoảng 10 lần (Task) $\rightarrow$ $20\text{ s}$ nghỉ (Rest).

3. Khử nhiễu và trích xuất đặc trưng

Dữ liệu Oxy-Hb thô được làm phẳng bằng bộ lọc Savitzky-Golay (FIR) với kích thước cửa sổ $N = 11$ và đa thức làm phẳng bậc 3, giúp triệt tiêu artifacts dịch chuyển cơ học mà vẫn bảo toàn nguyên vẹn đỉnh sóng biến động huyết động học. Tiếp đó, hệ thống các kênh đo (kênh 2, 9 từ đầu thu 1; kênh 5, 6 từ đầu thu 4) được ánh xạ tương đương thành hệ đa kênh thu - phát để đưa vào giải thuật FastICA nhằm ước lượng nguồn kích hoạt từ các optode phát.

---

Phát hiện chính (400-450 từ)

Quá trình khảo sát thực nghiệm và mô phỏng của đề tài đã mang lại 4 phát hiện khoa học quan trọng:

┌────────────────────────────────────────────────────────────────────────┐
│                        CÁC PHÁT HIỆN CHÍNH CỦA ĐỀ TÀI                  │
├───────────────────────────────────┬────────────────────────────────────┤
│ 1. Độ tin cậy giải thuật FastICA  │ 2. Giới hạn xử lý trong môi trường │
│    - Tách chính xác tần số cơ bản │    nhiễu (Noise Robustness)        │
│      (100Hz, 300Hz, 500Hz, 1000Hz)│    - Tách tốt nguồn ở SNR = 16 dB  │
│    - Giữ nguyên thành phần hài phổ│    - Cần kết hợp tiền lọc tín hiệu │
├───────────────────────────────────┼────────────────────────────────────┤
│ 3. Tối ưu tiền lọc Savitzky-Golay │ 4. Đặc trưng phổ fNIRS (6 - 7 Hz)  │
│    - Cửa sổ N=11, đa thức bậc 3   │    - Ngưỡng biên độ 0.05           │
│    - Triệt tiêu gai nhiễu artifacts│   - Phân biệt vận động 2 bán cầu   │
└───────────────────────────────────┴────────────────────────────────────┘

1. Tính hiệu quả và sự bảo toàn tần số trong tách nguồn mô phỏng

Thử nghiệm trên các cặp tín hiệu có phân bố phi Gauss đa dạng khẳng định FastICA có tốc độ hội tụ cực nhanh (chỉ sau vài bước lặp):

• Cặp Sin (300 Hz) và Răng cưa (100 Hz, duty cycle 70%): Tín hiệu ước lượng giữ chính xác thành phần hài cơ bản $100\text{ Hz}$ và $300\text{ Hz}$ trên miền tần số.
• Cặp Sin (3000 Hz) và Răng cưa (1000 Hz) cùng Sin (500 Hz) và Vuông (200 Hz): Phổ tần số sau khi tách phản ánh sự tương thích hoàn toàn với tín hiệu gốc.
• Đặc tính cố hữu: Do bản chất toán học của ICA, các tín hiệu phục hồi tồn tại sự không đồng định về dấu/pha (sign/phase ambiguity) và tỷ lệ biên độ (scale indeterminacy), nhưng thông tin dạng sóng và cấu trúc phổ được bảo toàn trọn vẹn.

2. Hành vi của FastICA trong môi trường tín hiệu có nhiễu

Khi kiểm thử với tín hiệu Sin ($800\text{ Hz}$) và Vuông ($1200\text{ Hz}$) bị tác động bởi tạp âm Gauss trắng ở mức tỷ số tín hiệu trên nhiễu $\text{SNR} = 16\text{ dB}$:

• FastICA vẫn tách thành công các thành phần tần số độc lập cốt lõi ($800\text{ Hz}$ và $1200\text{ Hz}$).
• Hạn chế phát hiện: Giải thuật FastICA thuần túy không có chức năng triệt tiêu thành phần nhiễu cộng; nhiễu vẫn tồn tại trên nền tín hiệu sau khi phân ly. Điều này chứng minh tính tất yếu của việc tích hợp các bộ lọc chuyên dụng trước khi đưa tín hiệu vào mô hình ICA.

3. Hiệu quả làm phẳng của bộ lọc Savitzky-Golay trên tín hiệu fNIRS

Tín hiệu huyết động học Oxy-Hb thu được từ máy FOIRE-3000 chịu ảnh hưởng lớn bởi dao động mạch máu cục bộ. Việc áp dụng bộ lọc Savitzky-Golay với tham số $N = 11$, bậc 3 cho thấy khả năng loại bỏ triệt để các xung gai đột biến (motion artifacts) mà không làm suy giảm biên độ hay làm trễ pha các biến đổi nồng độ Oxyhemoglobin sinh lý.

4. Phát hiện dấu hiệu phân biệt vận động não trong dải tần 6 – 7 Hz

Đây là phát hiện thực nghiệm quan trọng nhất của đề tài:

• Sau khi chuyển đổi hệ thống thu nhận quang học vùng motor sang mô hình "tiệc cocktail" và phân ly bằng FastICA, phân tích phổ Fourier cho thấy sự khác biệt rõ nét ở dải tần số $6\text{ Hz} - 7\text{ Hz}$ với ngưỡng biên độ chọn lọc là $0.05$.
• Khi đối tượng thực hiện gõ tay trái, các thành phần phổ ước lượng tại bán cầu não tương ứng có biên độ tiệm cận hoặc vượt ngưỡng $0.05$. Ngược lại, khi gõ tay phải, biên độ trong dải tần này nằm hoàn toàn dưới ngưỡng $0.05$.
• Khảo sát thống kê trên tập $64$ mẫu thực nghiệm mang lại kết quả:

$$\begin{array}{|l|c|c|} \hline \textbf{Hành động vận động} & \textbf{Độ phân biệt đúng theo nhận định} & \textbf{Tỷ lệ chính xác} \ \hline \text{Gõ tay trái} & 27 / 32 & 84.375% \ \text{Gõ tay phải} & 17 / 32 & 53.125% \ \hline \end{array}$$

Kết quả khẳng định tính khả thi vượt trội trong việc sử dụng phổ thành phần độc lập FastICA làm vector đặc trưng cho bài toán nhận dạng trạng thái não bộ.

---

Đóng góp khoa học (250-300 từ)

Đề tài nghiên cứu của ThS. Ngô Quốc Cường và ThS. Nguyễn Ngô Lâm đã mang lại những giá trị đóng góp thiết thực trên nhiều phương diện:

• Về mặt phương pháp luận (Methodological Innovations): Đề tài đã đề xuất thành công mô hình chuyển đổi không gian đo quang học fNIRS sang sơ đồ tương đương của bài toán BSS ("Cocktail Party"). Việc kết hợp chuỗi thuật toán: Tiền xử lý Savitzky-Golay $\rightarrow$ Chuẩn hóa/Trắng hóa EVD $\rightarrow$ Phân tích thành phần độc lập FastICA $\rightarrow$ Phân tích miền tần số đã thiết lập một quy trình chuẩn hóa, tối ưu cho việc khai phá dữ liệu huyết động học não.

• Về mặt ứng dụng thực tiễn (Practical Applications): Nghiên cứu chứng minh rằng dữ liệu quang phổ cận hồng ngoại chức năng hoàn toàn có thể trích xuất được các đặc trưng vận động cục bộ với độ tin cậy cao. Đây là tiền đề công nghệ đặc biệt quan trọng để xây dựng các hệ thống Giao diện Não - Máy tính (Brain-Computer Interface - BCI) không xâm lấn, hỗ trợ điều khiển thiết bị trợ vận động, xe lăn điện hoặc phục hồi chức năng sau tai biến mạch máu não.

• Về mặt giáo dục và chuyển giao học thuật (Academic & Educational Value): Đề tài cung cấp tài liệu nghiên cứu và thực nghiệm hệ thống về lý thuyết thống kê bậc cao, phép phân tích ICA và xử lý tín hiệu y sinh. Sản phẩm nghiên cứu được ứng dụng trực tiếp vào công tác giảng dạy các môn học chuyên ngành như Thực tập mô phỏng, Xử lý tín hiệu số, và công bố học thuật trong môi trường đại học.

---

Đối tượng quan tâm (200-250 từ)

Công trình nghiên cứu này mang lại giá trị tham khảo và ứng dụng chuyên sâu cho nhiều nhóm đối tượng:

• Nhà nghiên cứu Khoa học Thần kinh & Kỹ thuật Y sinh (Biomedical Engineers): Cung cấp tài liệu tham khảo thực nghiệm chi tiết về cách thức triển khai thiết bị fNIRS (đặc biệt là dòng máy FOIRE-3000), phương pháp xử lý dữ liệu hấp thụ đa bước sóng ($780\text{ nm}, 805\text{ nm}, 830\text{ nm}$) và cơ chế tính toán nồng độ Oxy-Hb/Deoxy-Hb/Total-Hb theo định luật Beer-Lambert.

• Kỹ sư Xử lý Tín hiệu số & AI/Machine Learning: Cung cấp mô hình tham chiếu thực tế trong việc giải quyết bài toán BSS bằng FastICA, phương pháp tối ưu hàm phi Gauss (Kurtosis/Negentropy), kỹ thuật làm phẳng dữ liệu sinh học và phương pháp trích chọn đặc trưng (feature extraction) phục vụ các mạng phân loại nơ-ron hay SVM.

• Nhà phát triển thiết bị Y tế & Hệ thống BCI (Brain-Computer Interface): Nắm bắt giải pháp thiết kế phần cứng/phần mềm thu nhận tín hiệu não không xâm lấn với độ trễ thấp, phục vụ các ứng dụng tương tác người - máy phục hồi chức năng y khoa.

• Giảng viên và Sinh viên chuyên ngành Điện tử - Viễn thông: Tài liệu học thuật trực quan phục vụ các học phần Xử lý tín hiệu nâng cao, Thực hành mô phỏng MATLAB và Đồ án tốt nghiệp.

---

FAQ - Các câu hỏi thường gặp (250-300 từ)

1. Phát hiện quan trọng nhất của đề tài nghiên cứu này là gì?

Phát hiện cốt lõi là giải thuật FastICA có thể áp dụng hiệu quả trên mô hình kênh đo fNIRS tương đương để trích xuất đặc trưng phổ Oxy-Hb trong dải tần $6 - 7\text{ Hz}$, tạo ra dấu hiệu định lượng (ngưỡng biên độ $0.05$) giúp phân biệt rõ rệt hoạt động vận động giữa tay trái và tay phải.

2. Điểm đặc biệt trong phương pháp luận (Methodology) của nghiên cứu là gì?

Điểm sáng tạo nằm ở việc mô hình hóa các cặp đầu phát quang (emitter) và đầu thu (detector) trên vỏ não thành hệ thống đa nguồn phát - đa kênh thu trong bài toán "tiệc cocktail", kết hợp bộ lọc Savitzky-Golay bậc 3 để khử nhiễu chuyển động trước khi áp dụng FastICA lặp điểm cố định.

3. Kết quả nghiên cứu có khả năng tổng quát hóa (Generalize) không?

Mô hình toán học FastICA có tính tổng quát tuyệt đối cho mọi bài toán BSS tuyến tính. Đối với tín hiệu fNIRS, nguyên lý trích đặc trưng phổ có thể áp dụng rộng rãi, tuy nhiên để đạt độ tin cậy cao trên diện rộng cần mở rộng tập dữ liệu thực nghiệm trên nhiều nhóm chủ thể với độ tuổi và thể trạng khác nhau.

4. Bước phát triển tiếp theo (Next Steps) của nghiên cứu là gì?

Tích hợp các thuật toán học máy phân loại tự động (như Support Vector Machine - SVM, Mạng Nơ-ron nhân tạo) trên vector đặc trưng FastICA thu được, đồng thời mở rộng thí nghiệm đa kênh trên cả hai bán cầu não và tăng cỡ mẫu thử nghiệm.

5. Những ứng dụng thực tế khả thi của nghiên cứu này là gì?

Ứng dụng trực tiếp trong xây dựng hệ thống BCI điều khiển chi giả sinh học, thiết bị hỗ trợ người tàn tật, đánh giá vùng vỏ não tổn thương trong y học phục hồi chức năng và giám sát trạng thái nhận thức/tập trung của con người.

---

Kết luận (150 từ)

Công trình nghiên cứu khoa học "Ứng dụng giải thuật FastICA trong tách nguồn mù và trích đặc trưng" đã chứng minh tính hiệu quả vượt trội của thuật toán phân tích thành phần độc lập điểm cố định trong cả bài toán mô phỏng tách nguồn lẫn bài toán thực nghiệm y sinh phức tạp. Việc khai phá thành công tín hiệu quang phổ cận hồng ngoại chức năng (fNIRS) mở ra hướng đi đầy tiềm năng cho lĩnh vực kỹ thuật thần kinh tại Việt Nam.

Trong tương lai, việc kết hợp FastICA với các mô hình học sâu (Deep Learning) và mở rộng quy mô mẫu thử nghiệm hứa hẹn sẽ hoàn thiện các hệ thống nhận dạng hành vi não bộ theo thời gian thực với độ chính xác cao.

Trích đoạn nội dung tài liệu

Mở Đầu  Kết quả của đề tài có thể đƣợc sử dụng vào việc nghiên cứu việc giảng dạy cho sinh viên chuyên ngành điện tử viễn thông. Mục tiêu đề tài  Tách nguồn mù từ các tín hiệu trộn.  Trích đặc trƣng của đối tƣợng. Cách tiếp cận, phương pháp nghiên cứu, phạm vi nghiên cứu Cách tiếp cận  Lý thuyết đến thực nghiệm Phƣơng pháp nghiên cứu  Nghiên cứu lý thuyết  Mô phỏng  Thực nghiệm Phạm vi nghiên cứu  Ứng dụng giải thuật fastICA trong tách nguồn mù (nguồn tín hiệu mô phỏng) và trích đặc trƣng của đối tƣợng (tín hiệu NIRS).

Nội dung nghiên cứu  Cơ sở lý thuyết phân tích thành phần độc lập ICA  Thuật toán fast ICA  Ứng dụng thuật toán trong tách nguồn mù  Ứng dụng thuật toán trong trích đặc trƣng trên tín hiệu NIRS 2 Luan van 1. Cơ Sở Lý Thuyết Chƣơng 1 CƠ SỞ LÝ THUYẾT 1.1 Lý Thuyết Chung Về ICA 1.1 Mô Hình ICA Cơ Bản Để định nghĩa ICA [3], ta có thể sử dụng một mô hình thống kê „biến ẩn‟. Giả sử ta quan sát n biến ngẫu nhiên x1,x2,. xn, các biến này đƣợc biểu diễn theo sự kết hợp tuyến tính với của n biến ngẫu nhiên s1,s2.sn: xi = ai1s1 + ai2s2 +…+ ainsn (1.,n là các hệ số thực.

Theo định nghĩa, sj là độc lập thống kê với nhau. Đây là mô hình ICA cơ bản. Nó mô tả dữ liệu thu đƣợc bằng cách trộn các hành phần sj với nhau. Thành phần độc lập sj (viết tắt: ICs) là biến ẩn, không thể quan sát trực tiếp đƣợc.

Ngoài ra các hệ số trộn aij cũng đƣợc giả sử là không biết. Tất cả những gì ta quan sát đƣợc là các biến ngẫu nhiên xi, và ta phải ƣớc lƣợng các hệ số trộn aij và các thành phần sj bằng cách sử dụng xi. Ở đây biến thời gian t đã đƣợc bỏ qua. Bởi vì trong mô hình ICA cơ bản này, ta giả định rằng xi cũng nhƣ mỗi thành phần độc lập sj là biến ngẫu nhiên, thay vì biểu diễn theo thời gian.

Nhƣ vậy, ta có thể bỏ qua thời gian trễ trong quá trình kết hợp các thành phần độc lập chƣa biết. ICA có thể đƣợc miêu tả nhƣ bài toán tách nguồn mù. Nguồn ở đây chính là các tín hiệu gốc – thành phần độc lập chƣa biết. Do thông tin về sự kết hợp này cũng chƣa biết trƣớc nên ta có thể so sánh mô hình trong công thức (1.1) tƣơng ứng với bài toán tách nguồn mù và đƣợc sử dụng làm cơ sở cho các ứng dụng về sau.

Sử dụng biểu diễn theo dạng ma trận cho công thức (1.2) Trong đó, x đại diện cho x1,x2,. và s là các thành phần s1.sn, ma trận A với các phần tử aij.2 Giả thiết trong ICA Mô hình ICA cơ bản đặt ra các giả thiết sau cho việc phân ly (ƣớc lƣợng) các thành phần độc lập:  Các nguồn s độc lập thống kê nhau, nghĩa là biết đƣợc một nguồn không thể suy ra các nguồn còn lại.  Các hàm phân bố xác suất của các nguồn có trị trung bình bằng không.  Không có nguồn (thành phần độc lập) nào có phân bố Gauss.

 Ma trận trộn A là ma trận vuông tức số lƣợng nguồn và số lƣợng trộn bằng nhau.2 Ƣớc lƣợng Trong ICA Mô hình ICA đặt ra một hạn chế là các thành phần độc lập phải có tính phi Gauss (non-gaussianity), tức không có phân bố (hàm mật độ xác suất) là Gauss. Lý do tính phi Gauss nằm ở chỗ là các biến ngẫu nhiên Gauss đƣợc xác định hoàn toàn bởi các thống kê bậc một (trị trung bình) và bậc hai (phƣơng sai), các thống kê bậc cao hơn bằng không. Trong lúc đó, mô hình ICA cần các thống kê bậc cao hơn của các thành phần độc lập để thực hiện sự phân ly (ƣớc lƣợng các thành phần độc lập). Nhƣ vậy, sự phi tuyến, tính phi Gauss dẫn đến sự độc lập thống kê.

Nhiều giải thuật khác nhau đã đƣợc đề xuất để giải quyết bài toán ƣớc lƣợng các thành phần độc lập trong ICA:  Cực đại hóa tính phi Gauss (non-gaussianity)  Ƣớc lƣợng khả năng cực đại (maximum likelihood)  Cực tiểu hoá thông tin tƣơng hỗ (mutual information) 1.1 Đo Tính Phi Gauss Bằng Kurtosis Đầu tiên là phép đo dựa trên kurtosis của một biến ngẫu nhiên y có trung bình bằng không là cumulant bậc bốn: kurt( y) = E{ y 4 } - 3 ( E{ y 2 }) 2 (1.3) Thật ra vì ta giả sử y có phƣơng sai đơn vị, nên kurtosis là kurt(y)=E{y4}–3, tức kurtosis là phiên bản chuẩn hóa của momen thứ tƣ E{y4}. Khi y có phân bố Gauss momen thứ tƣ bằng 3(E{y2})2 nên kurtosis bằng không đối với các biến ngẫu nhiên Gauss. Hầu hết các biến ngẫu nhiên không phải Gauss thì kurtosis khác không. Nếu kurtosis là dƣơng thì biến ngẫu nhiên có phân bố siêu Gauss (supergaussian), còn nếu kurotsis là âm thì biến ngẫu nhiên có phân bố dƣới Gauss (subgaussian).

Phân bố siêu Gauss không còn dạng hình chuông nhƣ Gauss mà tăng nhanh ở trung tâm tƣơng tự nhƣ phân bố Laplace, còn phân bố dƣới Gauss không nhô lên ở phần giữa nhƣ Gauss mà tiến đến phân bố đều với biên độ rất nhỏ ở xa trung tâm.2 Đo Tính Phi Gauss Bằng Negentropy Một phép đo tính phi Gauss khác là negentropy. Negentropy là đại lƣợng dựa trên lý thuyết thông tin gọi là entropy vi sai. Entropy của một biến ngẫu nhiên là số đo lƣợng thông tin trung bình của nó.4) Đặc tính quan trọng của entropy là biến ngẫu nhiên Gauss có entropy lớn nhất trong các biến ngẫu nhiên có cùng phƣơng sai. Nhƣ vậy entropy, và negentropy định nghĩa theo entropy, có thể dùng để đo tính phi Gauss của một biến ngẫu nhiên.

Để có đƣợc một số đo tính phi Gauss sao cho bằng không đối với biến Gauss và luôn không âm, ta định nghĩa negentropy của vectơ ngẫu nhiên y: 4 Luan van 1.5) Trong đó yGauss là một vectơ ngẫu nhiên Gauss cùng ma trận hiệp phƣơng sai (hay ma trận tƣơng quan vì các dữ liệu đƣợc giả sử có trung bình là không). Tuy nhiên tính toán negentropy lại khó khăn. Một số tính toán xấp xỉ đã đƣợc phát triển.3 Tiền Xử Lý trong ICA Tiền xử lý là quá trình xử lý tín hiệu thu đƣợc x trƣớc khi thực hiện ƣớc lƣợng tín hiệu gốc. Tiền xử lý đƣợc sử dụng để việc giải bài toán ICA trở nên đơn giản hơn, gồm có quy tâm và trắng hóa tín hiệu.1 Qui Tâm Tín Hiệu Nhƣ đã nêu ở trƣớc, các tín hiệu trộn quan sát đƣợc x phải có trị trung bình m=E[x] bằng không, điều này cũng có nghĩa là các tín hiệu nguồn s cũng có trị trung bình là không.

Nếu các tín hiệu chƣa có trị trung bình là không ta thực hiện phép qui tâm tức trừ phân bố của các biến ngẫu nhiên với các trị trung bình của chúng: x = x - E{x} (1.6) Trong đó, x′ là vectơ ngẫu nhiên chƣa có trung bình là không. Sau khi đã ƣớc lƣợng ma trận A và các thành phần s ta có thể thêm trở lại các trị trung bình của chúng. Khi vectơ ngẫu nhiên x (hoặc s) có trị trung bình là không thì hiệp phƣơng sai và tƣơng quan của nó giống nhau.2 Trắng Hoá Sau khi đã quy tâm, ta thực hiện một biến đổi trên x để đƣợc một vector mới gọi là vector trắng. Mục đích của việc làm trắng này là làm cho dữ liệu không tƣơng quan.

x đƣợc gọi là trắng khi E{xxT }  I (1.7) Việc làm trắng là một biến đổi tuyến tính thông qua ma trận V z = Vx (1.8) Trong đó x là dữ liệu cần làm trắng và V là ma trận làm trắng, z là dữ liệu đã trắng hóa. Một trong những phƣơng pháp phổ biến để trắng hóa là sử dụng phƣơng pháp phân giải trị riêng (EVD) của ma trận hiệp phƣơng sai.9) Với E là ma trận trực giao của các vectơ riêng của E{xxT} và D là đƣờng chéo ma trận của các giá trị riêng của nó, D = diag(d1,…,d2). Ma trận làm trắng đƣợc tính: V  ED1/2 ET (1. Cơ Sở Lý Thuyết Ma trận D-1/2 đƣợc tính bằng một thao tác đơn giản D1/ 2  diag( d11/2 ,.

Ma trận trắng hóa tính toán theo cách này đƣợc ký hiệu là E{xxT}-1/2 hoặc C-1/2. Trắng hóa biến đổi ma trận trộn thành ma trận mới. Ta có ~ z  VAs  As (1.11) ~ Sau khi đã ƣớc lƣợng đƣợc ma trận trắng hóa A , các thành phần độc lập đƣợc tính s  Wz (1. Từ đây, ma trận A đƣợc ƣớc lƣợng ~ A  VA  D1 / 2ET A (1.13) ~ A  ED1 / 2 A  ED1 / 2 W 1  ED1 / 2 WT (1.14) Vì W trực giao nên W1  WT .4 ICA dùng phƣơng pháp cực đại tính phi Gauss Để thực hiện ƣớc lƣợng tính phi Gauss trong ICA, kurtosis đƣợc sử dụng.

Nhƣ đã đề cập, kurtosis thể hiện cumulant bậc 4 của một biến ngẫu nhiên. Kurtosis của y, kurt(y) đƣợc định nghĩa kurt( y )  E{ y 4 }  3( E{ y 2 }) 2 (1.15) Trong thực tế, để cực đại trị tuyệt đối của kurtosis, ta sẽ bắt đầu từ vector w, tính toán hƣớng mà trị tuyệt đối kurtosis của y  wT z tăng nhanh trên các mẫu z(1), …, z(T), và sau đó chuyển vector w theo hƣớng đó. Ý tƣởng này đƣợc thực thi theo phƣơng pháp gradient. Gradient của trị tuyệt đối kurtosis của wT z có thể đƣợc tính  kurt(w T z )  4sign(kurt(w T z ))[ E{z (w T z )3}  3w w ] 2 (1.16) w Ta có đƣợc giải thuật gradient w  sign(kurt(wT z ))[ E{z (wT z )3} (1.18) Để tăng tốc độ hội tụ, giải thuật lặp điểm cố định có thể đƣợc sử dụng trong việc tìm w.

Sử dụng các nhân tử Lagrange, từ (1.19) Từ đây, ta tính giá trị mới cho w w  E{z (wT z )3}  3w (1.20) Sau mỗi vòng lặp, w đƣợc chuẩn hóa bằng cách chia cho độ lớn của nó (vì vậy w  1 và có thể bỏ qua trong công thức 1. Cơ Sở Lý Thuyết Sự hội tụ của giải thuật lặp điểm cố định đó là giá trị w cũ và mới đều có cùng hƣớng, tức là nhân chấm của chúng bằng 1. Giải thuật này đƣợc chứng minh là hội tụ nhanh và đáng tin cậy nên đƣợc gọi là FastICA.5 Phƣơng Pháp Tái Hiện Hình Ảnh Và Tín Hiệu Não Ngƣời 1.1 Phƣơng pháp EEG EEG (Electroencephalography) – điện não đồ là một phƣơng pháp ghi lại các hoạt động về điện trên da đầu của ngƣời. Phƣơng pháp này đo sự thay đổi điện áp do việc di chuyển của các hạt điện tích bên trong các nơ- ron của não.