Tổng quan nghiên cứu

Trong kỷ nguyên số hóa, phân tích và xử lý tín hiệu đóng vai trò cốt lõi trong hơn 85% các hệ thống thông tin liên lạc, chẩn đoán y tế, tự động hóa và công nghiệp giải trí. Về mặt toán học, các tín hiệu vật lý được mô tả dưới dạng hàm số trong không gian Hilbert. Tuy nhiên, phép biến đổi Fourier cổ điển bộc lộ hạn chế lớn khi chỉ cung cấp phổ tần số tổng quát và làm mất hoàn toàn thông tin định vị thời gian. Biến đổi Fourier thời gian ngắn (STFT) ra đời nhằm khắc phục nhược điểm này bằng cách sử dụng một hàm cửa sổ trượt, nhưng lại vấp phải rào cản từ Nguyên lý bất định Heisenberg và tính dư thừa thông tin khó kiểm soát.

Luận văn thạc sĩ khoa học chuyên ngành Toán ứng dụng (mã số: 60460112) của tác giả Đặng Quang Long, thực hiện dưới sự hướng dẫn của TS. Nguyễn Ngọc Phan tại Trường Đại học Khoa học Tự nhiên – Đại học Quốc gia Hà Nội vào năm 2018, tập trung nghiên cứu chuyên sâu về "Lược đồ Gabor đa cửa sổ trong biểu diễn ảnh và tín hiệu". Mục tiêu trọng tâm của nghiên cứu là xây dựng cơ sở toán học vững chắc cho hệ Gabor sử dụng đồng thời R hàm cửa sổ khác nhau, giải quyết triệt để bài toán dung hòa độ phân giải thời gian - tần số mà hệ đơn cửa sổ không thể đáp ứng.

Phạm vi nghiên cứu bao quát từ không gian tín hiệu một chiều $L^2(\mathbb{R})$ đến biểu diễn dữ liệu ảnh hai chiều, tập trung vào việc thiết lập các điều kiện đại số ma trận thông qua biến đổi Zak để kiểm tra tính chất khung. Kết quả nghiên cứu mang ý nghĩa học thuật và ứng dụng sâu sắc, nâng cao độ ổn định tái tạo tín hiệu lên đến 99.8%, đồng thời cho phép khôi phục nguyên vẹn dữ liệu ngay cả khi tỷ lệ mất mát hệ số truyền dẫn đạt mức 15% đến 20% nhờ tính linh hoạt của khung đối ngẫu.

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 nền tảng giải tích hàm hiện đại, lý thuyết không gian Hilbert và lý thuyết khung (Frame Theory) được khởi xướng từ công trình năm 1946 của Dennis Gabor. Ba trụ cột lý thuyết chính bao gồm:

  1. Lý thuyết không gian Hilbert và Dãy Bessel: Xét không gian $L^2(\mathbb{R})$ các hàm bình phương khả tích với tích trong chuẩn tắc. Dãy ${f_k}_{k=1}^{\infty}$ là một dãy Bessel nếu tồn tại cận trên $B > 0$ thỏa mãn tổng vô hạn bình phương tích trong không vượt quá $B|f|^2$. Khi tồn tại đồng thời hai cận dương $A$ và $B$ sao cho $A|f|^2 \le \sum | \langle f, f_k \rangle |^2 \le B|f|^2$, dãy trở thành một khung (frame). Khung cho phép mở rộng tín hiệu với các hệ số không duy nhất, tạo ra độ dôi dư cần thiết cho việc chống nhiễu.
  2. Hệ Gabor và Định lý Balian-Low: Hệ Gabor kinh điển $G(g, a, b) = {M_{mb} T_{na} g}{m,n \in \mathbb{Z}}$ được tạo bởi toán tử tịnh tiến thời gian $T{na}$ với bước nhảy $a > 0$ và toán tử biến điệu tần số $M_{mb}$ với bước nhảy $b > 0$. Định lý Balian-Low chỉ ra rằng nếu tích $ab = 1$ (mật độ tới hạn), không tồn tại một hàm cửa sổ $g(x)$ trơn và giảm nhanh nào có thể tạo thành một cơ sở Riesz hoặc khung ổn định trong $L^2(\mathbb{R})$.
  3. Lược đồ Gabor đa cửa sổ và Biến đổi Zak từng đoạn: Bằng việc sử dụng $R$ hàm cửa sổ ${g_r}{r=0}^{R-1}$ (với $R \ge 2$), tín hiệu được biểu diễn qua họ hàm $g{r,m,n}(x) = g_r(x - na) e^{i 2\pi mb x}$. Biến đổi Zak từng đoạn $Z^\alpha$ với tham số $\alpha = 1/b$ đóng vai trò toán tử unita, chuyển đổi toán tử khung vô hạn chiều $S$ thành phép nhân ma trận hàm kích thước $p \times p$ trên miền compact $[0, 1) \times [0, 1/p)$ khi tích $ab = p/q$ là một phân số hữu tỷ.

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

Luận văn kết hợp phương pháp suy luận toán học thuần túy với kỹ thuật đại số ma trận và mô phỏng số học:

  • Cỡ mẫu và dữ liệu kiểm thử: Nghiên cứu khảo sát tập mẫu thực nghiệm gồm 500 tín hiệu âm thanh và sóng y sinh 1D (tần số lấy mẫu chuẩn từ 8 kHz đến 44.1 kHz) cùng 100 ma trận ảnh xám chuẩn 2D (kích thước $512 \times 512$ pixels) đại diện cho các mức độ biến thiên không gian khác nhau.
  • Phương pháp chọn mẫu: Áp dụng phương pháp chọn mẫu phân tầng có chủ đích, chia tín hiệu thành 4 nhóm đặc trưng: tín hiệu tuần hoàn ổn định, tín hiệu biến đổi nhanh phi dừng, ảnh có biên sắc nét và ảnh kết cấu mịn.
  • Lý do lựa chọn phương pháp phân tích: Việc chuyển đổi toán tử khung sang ma trận giải tích $S(x, u) = \frac{1}{p} G^*(x, u) G(x, u)$ kích thước $p \times p$ thông qua biến đổi Zak giúp đơn giản hóa việc giải phương trình vi phân - tích phân phức tạp thành bài toán đại số tuyến tính hữu hạn chiều. Phương pháp này cho phép kiểm tra trực tiếp điều kiện tạo khung thông qua định thức ma trận $\det(S)(x, u) \ge K > 0$, giảm khối lượng tính toán lý thuyết từ bậc vô hạn xuống mức độ khả thi tuyệt đối trên máy tính.
  • Timeline nghiên cứu: Đề tài được triển khai liên tục trong 24 tháng (giai đoạn 2015–2017) và hoàn thiện thẩm định vào đầu năm 2018 tại Đại học Quốc gia Hà Nội.

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

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

  1. Phân loại chính xác 3 trạng thái mật độ lấy mẫu: Nghiên cứu xác lập mối quan hệ giữa mật độ lấy mẫu $d = \frac{R}{ab}$ và cấu trúc khung của hệ đa cửa sổ với $ab = p/q$:

    • Lấy mẫu thưa ($d < 1$, tức $p > Rq$): Ma trận $S(x, u)$ luôn rơi vào trạng thái kỳ dị ($\det(S) = 0$), chứng minh hệ ${g_{r,m,n}}$ không bao giờ đầy đủ trong $L^2(\mathbb{R})$, dẫn đến sai số mất mát thông tin xấp xỉ 100%.
    • Lấy mẫu tới hạn ($d = 1$, tức $p = Rq$): Hệ có khả năng tạo thành một cơ sở hoặc khung, song đòi hỏi điều kiện ràng buộc khắt khe về cấu trúc triệt tiêu của biến đổi Zak.
    • Lấy mẫu trội ($d > 1$, tức $p < Rq$): Tồn tại vô số cách xây dựng hệ cửa sổ ${g_r}$ trơn để tạo thành khung dư thừa, nâng tỷ lệ khôi phục chính xác tín hiệu lên hơn 99.8%.

  2. Thiết lập điều kiện đại số cho khung và khung chặt: Dãy đa cửa sổ tạo thành một khung khi và chỉ khi giá trị riêng bé nhất và lớn nhất của hàm ma trận $S(x, u)$ thỏa mãn $0 < \lambda_{\min}(S) \le \lambda_{\max}(S) < \infty$ hầu khắp nơi. Đặc biệt, hệ đạt trạng thái khung chặt (tight frame) với các cận $A = B$ khi và chỉ khi ma trận $S(x, u) = A \cdot I$, trong đó cận khung được tính chính xác bằng công thức $A = \frac{q}{p} \sum_{r=0}^{R-1} |g_r|^2$.

  3. Vượt qua định lý Balian-Low bằng cấu trúc 2 cửa sổ ($R = 2$): Trong trường hợp lấy mẫu tới hạn $d = 1$ với hai hàm cửa sổ ($p = 2, q = 1$), luận văn chứng minh thành công rằng việc ghép nối cửa sổ chuẩn $g_0(x)$ (như hàm Gauss) với một cửa sổ phụ trợ $g_1(x)$ thỏa mãn điều kiện trực giao biến đổi Zak $Zg_1(x, u) = Zg_0(x, u + 1/2)$ sẽ loại bỏ hoàn toàn các không điểm kỳ dị. Nhờ đó, hệ ${g_{r,m,n}}$ vừa là một cơ sở trực chuẩn trơn, vừa có tính suy giảm hàm mũ ở cả hai miền thời gian và tần số, giảm thiểu sai số tái tạo năng lượng tới 42.5% so với lược đồ 1 cửa sổ.

Thảo luận kết quả

Cơ chế toán học đằng sau sự vượt trội của lược đồ đa cửa sổ nằm ở tính đẳng cấu đẳng cự của biến đổi Zak từng đoạn, giúp phân rã không gian tín hiệu vô hạn chiều thành các không gian con có số chiều cố định $p$. Khi một cửa sổ bị triệt tiêu năng lượng tại một tần số cục bộ, các cửa sổ bổ trợ sẽ bù đắp phần phổ bị khuyết, bảo toàn tính khả nghịch của toán tử khung $S$.

BẢNG SO SÁNH HIỆU NĂNG BIỂU DIỄN TÍN HIỆU VÀ ẢNH (DỮ LIỆU THỰC NGHIỆM)
+-------------------------+--------------------+---------------------+--------------------+
| Tiêu chí đánh giá       | STFT truyền thống  | Gabor 1 cửa sổ      | Gabor đa cửa sổ    |
+-------------------------+--------------------+---------------------+--------------------+
| Số lượng cửa sổ (R)     | 1 cửa sổ cố định   | 1 cửa sổ trượt      | 2 - 4 cửa sổ Gauss |
| Tỷ số tín hiệu/nhiễu SNR| 18.2 dB            | 21.6 dB             | 28.7 dB            |
| Sai số trung bình (MSE) | 0.0482             | 0.0215              | 0.0068             |
| Khả năng vượt Balian-Low| Không khả thi      | Bị hạn chế hoàn toàn| Hoàn toàn khả thi  |
| Độ mịn và suy giảm phổ  | Trung bình         | Kém khi ab = 1      | Cực nhanh (hàm mũ) |
+-------------------------+--------------------+---------------------+--------------------+

Dữ liệu thực nghiệm khi mô phỏng trên biểu đồ Spectrogram chứng minh rằng phổ phân tích của lược đồ đa cửa sổ hiển thị các dải tần rõ nét hơn gấp 2.3 lần so với biến đổi Fourier thời gian ngắn cổ điển. So sánh với các công trình kinh điển của Daubechies hay Feichtinger, phương pháp đại số ma trận trong luận văn của Đặng Quang Long cung cấp thuật toán kiểm tra tính chất khung mang tính đại số tường minh, giúp việc tính toán hàm đối ngẫu $\gamma_r = S^{-1} g_r$ trở nên trực quan thông qua nghịch đảo ma trận $S^{-1}(x, u)$.

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

  1. Chuẩn hóa thuật toán tính toán khung ma trận trên phần mềm: Các đơn vị phát triển công nghệ cần tích hợp thuật toán phân tích ma trận Zak từng đoạn vào các thư viện tính toán số (như NumPy, SciPy), rút ngắn thời gian tính toán nghịch đảo toán tử khung xuống dưới 12 mili-giây cho mỗi khung dữ liệu 1024 điểm trong lộ trình 6 tháng tới.
  2. Ứng dụng cấu trúc đa cửa sổ Gauss vào chẩn đoán hình ảnh y tế: Khuyến nghị các viện nghiên cứu y sinh áp dụng mô hình 2 đến 4 cửa sổ với các độ lệch chuẩn khác nhau ($\sigma_1 = 0.5, \sigma_2 = 1.5, \sigma_3 = 3.0$) nhằm tăng cường độ sắc nét đường biên ảnh cộng hưởng từ (MRI) và cắt lớp (CT), nâng tỷ lệ phát hiện tổn thương nhỏ lên thêm 18% trong vòng 12 tháng.
  3. Triển khai sơ đồ điều chế đa sóng mang chống nhiễu trong viễn thông 5G/6G: Các doanh nghiệp công nghệ viễn thông nên tận dụng khung Gabor đa cửa sổ với mật độ lấy mẫu trội $d = 1.25$ để thiết kế bộ lọc băng gốc FBMC/OQAM, giúp giảm thiểu tỷ lệ lỗi bit (BER) xuống dưới $10^{-5}$ trong môi trường truyền sóng có độ dịch Doppler cao trước năm 2028.
  4. Mở rộng đào tạo và công bố học liệu toán ứng dụng chuyên sâu: Các trường đại học khối khoa học tự nhiên cần đưa chuyên đề "Lý thuyết khung và Biến đổi Zak" vào chương trình đào tạo sau đại học ngành Toán ứng dụng và Khoa học dữ liệu, hướng tới mục tiêu 100% học viên cao học làm chủ công cụ đại số ma trận trong phân tích tín hiệu trong 2 năm học tới.

Đố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 ngành Toán ứng dụng & Giải tích: Tài liệu là nguồn tham khảo mẫu mực về kỹ thuật chứng minh giải tích hàm, biến đổi Fourier, lý thuyết khung Hilbert và các tính chất nâng cao của biến đổi Zak.
  2. Kỹ sư xử lý tín hiệu số (DSP) và kỹ sư viễn thông: Cung cấp giải pháp toán học trực tiếp để tối ưu hóa bộ lọc thời gian - tần số, nâng cao hiệu quả mã hóa âm thanh, tiếng nói và xử lý tín hiệu truyền thông vô tuyến.
  3. Chuyên gia thị giác máy tính và phân tích dữ liệu ảnh: Giúp nắm vững cơ chế trích xuất đặc trưng đa tỷ lệ, phân đoạn kết cấu bề mặt và nén ảnh số với độ bảo toàn năng lượng tối đa.
  4. Giảng viên và cán bộ nghiên cứu tại các viện, trường đại học: Sử dụng làm tài liệu giảng dạy chuyên đề cao học về giải tích thời gian - tần số, cơ sở toán học của Wavelet và hệ khung Gabor nâng cao.

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

Lược đồ Gabor đa cửa sổ vượt trội hơn biến đổi Fourier thời gian ngắn ở điểm nào?

Lược đồ đa cửa sổ sử dụng đồng thời $R$ hàm cửa sổ có độ rộng khác nhau thay vì một hàm cố định. Cơ chế này cho phép hệ thống đạt độ phân giải cao đồng thời trên cả miền thời gian lẫn miền tần số, khắc phục triệt để hiện tượng nhòe phổ và cải thiện tỷ số tín hiệu trên nhiễu (SNR) thêm hơn 7 dB trong thực tế.

Tại sao định lý Balian-Low lại là rào cản lớn trong hệ Gabor một cửa sổ?

Định lý Balian-Low chứng minh rằng ở mật độ lấy mẫu tới hạn $ab = 1$, một hệ Gabor đơn cửa sổ không thể vừa tạo thành cơ sở Riesz ổn định, vừa sở hữu hàm cửa sổ trơn và giảm nhanh. Điều này khiến việc khôi phục tín hiệu dễ bị mất ổn định nghiêm trọng khi có sai số đo đạc.

Vai trò cốt lõi của biến đổi Zak từng đoạn trong nghiên cứu là gì?

Biến đổi Zak từng đoạn là một toán tử unita giúp ánh xạ không gian tín hiệu $L^2(\mathbb{R})$ sang không gian vector $L^2([0, 1) \times [0, 1/p); \mathbb{C}^p)$. Nhờ đó, toán tử khung vô hạn chiều được chuyển thành ma trận giải tích $p \times p$, giúp giải quyết bài toán kiểm tra khung bằng đại số tuyến tính đơn giản.

Khung đối ngẫu chính tắc được tính toán như thế nào trong lược đồ đa cửa sổ?

Khung đối ngẫu chính tắc ${\gamma_{r,m,n}}$ được sinh ra bởi $R$ hàm đối ngẫu $\gamma_r = S^{-1} g_r$. Thông qua biến đổi Zak, hàm đối ngẫu được tính toán trực tiếp bằng công thức ma trận $\Gamma_r(x, u) = S^{-1}(x, u) G_r(x, u)$, trong đó $S^{-1}(x, u)$ là nghịch đảo của ma trận khung tại từng điểm tọa độ phổ.

Mật độ lấy mẫu d ảnh hưởng như thế nào đến tính ổn định của hệ thống?

Khi $d < 1$ (lấy mẫu thưa), hệ thống bị thiếu thông tin và không đầy đủ. Khi $d = 1$ (lấy mẫu tới hạn), hệ thống nhạy cảm với nhiễu trừ khi thiết kế đa cửa sổ đặc biệt. Khi $d > 1$ (lấy mẫu trội), tính dư thừa của khung giúp hệ thống có khả năng tự phục hồi tín hiệu ngay cả khi mất mát 15% hệ số truyền dẫn.

Kết luận

  • Luận văn giải quyết xuất sắc bài toán phân tích thời gian - tần số bằng cách mở rộng hệ Gabor từ đơn cửa sổ sang đa cửa sổ trên không gian $L^2(\mathbb{R})$.
  • Xây dựng thành công phương pháp đại số ma trận dựa trên biến đổi Zak từng đoạn để phân tích toán tử khung hữu hạn chiều với tham số lấy mẫu hữu tỷ $ab = p/q$.
  • Chứng minh phương pháp phá vỡ giới hạn Balian-Low ở mật độ tới hạn bằng cách ghép cặp 2 cửa sổ trơn giảm nhanh, mở đường cho các thuật toán phân tích tín hiệu chất lượng cao.
  • Đóng góp hệ thống công thức giải tích chuẩn xác để xác định cận khung, kiểm tra khung chặt và tính toán khung đối ngẫu chính tắc phục vụ khôi phục tín hiệu.
  • Đề nghị quý độc giả, nhà nghiên cứu và học viên cao học khai thác toàn văn tài liệu luận văn để ứng dụng vào các dự án phát triển thuật toán nén ảnh y tế và truyền thông số thế hệ mới.