Nhân ma trận, DFT, và lý thuyết biểu diễn nhóm (2)

Ngô Quang Hưng | 25 tháng 11, 2005 | Bản để in Bản để in

Tiếp theo bài trước, lần này ta bàn sơ qua về DFT. Tôi học DFT lần đầu tiên vào khoảng năm 1993. Học xong thấy nó rất bí hiểm, theo kiểu: nếu lấy vector này, tính toán thế này, thì ra các hệ số thế kia, nhưng không hiểu ý tưởng nằm sau các công thức đó. Sau này tôi mới khám phá ra là DFT chẳng qua là một phép thay đổi cơ sở (basis) trong không gian tuyến tính. Từ đó thấy mọi thứ rõ ràng, dễ hiểu hẳn ra.

Một trong những thước đo độ sâu và tầm quan trọng của một ý tưởng hay một nhánh nghiên cứu là ứng dụng của nó vào nhiều nơi (kể cả những chỗ mà ta không ngờ trước được là ý tưởng có thể có ứng dụng). Harmonic analysis nói chung và Discrete Fourier Transform nói riêng là một trong những nhánh như vậy. Ứng dụng của Fourier transforms có ở hầu hết tất cả mọi ngóc ngách của toán học, vật lý, khoa học máy tính, điện/điện tử, xác suất thống kê, vân vân. Ví dụ: Hastad đã dùng Fourier analysis để chứng minh một định lý rất quan trọng trong hệ thống PCP, các phân tích tốt nhất về expansion rate của một vài phép xây dựng expander graphs đều dùng Fourier analysis, v.v. Bạn có thể xem thêm vài lecture notes của giáo sư Nathan Linial về một số ứng dụng của Fourier analysis trong toán rời rạc, và master thesis này có thảo luận vài ứng dụng trong lý thuyết tính toán. (Các links đã dẫn chỉ là chóp đỉnh của tảng băng harmonic analysis khổng lồ!)

Về căn bản, DFT không có gì ghê gớm lắm. Căn bản đại số tuyến tính các bạn có thể xem một lecture note tôi tóm tắt vài thứ cần thiết. Hãy xét một ví dụ (tương đối) đơn giản sau đây. Xét một tập X hữu hạn và tập V = \{ f \ | \ f: X \to \mathbb{C} \}. Dễ thấy rằng nếu |X| = n thì V có thể xem như \mathbb{C}^n (nói chính xác hơn thì V\mathbb{C}^n isomorphic với nhau). Giả sử ta có một orthonormal basis (hệ cơ sở trực chuẩn) \{u_1, \dots, u_n\} của không gian V, thì mọi vector f \in V đều có thể viết được dưới dạng:

f = \hat{f}_1u_1 + \hat{f}_2u_2 + \dots + \hat{f}_nu_n

(Ở đây ta ngầm hiểu các vectors đều là vectors cột.) Các hệ số \hat{f}_j cũng tạo thành một vector mà ta gọi là \hat{f}. Vector này chính là DFT của f, ký hiệu là \hat{f} = DFT(f). Cái DFT mà ta thường gặp trong xử lý tín hiệu số tương ứng với một trường hợp đặc biệt khi mà
u_j^T = [ \omega_n^{0}, \omega_n^{j}, \omega_n^{2j}, \dots, \omega_n^{(n-1)j} ]

Trong đó \omega_n = e^{-2\pi i/n}nth root of unity.

Trong miền số phức, ta thường dùng hermitian form (hiểu nôm na là tích vô hướng của hai vectors) sau đây:

\langle f, g \rangle = f_1^*g_1 + \cdots + f_n^*g_n

Vì các vectors \{u_1, \dots, u_n\} là một orthonormal basis (nghĩa là chúng vuông góc với nhau và đều có chiều dài bằng 1, dễ thấy rằng
\langle f, g \rangle = \langle \hat{f}, \hat{g} \rangle

Đẳng thức này thường được gọi là đẳng thức Parseval. Đôi khi người ta gọi trường hợp đặc biệt khi f = gđẳng thức Parseval, hoặc đẳng thức Plancherel:
\|f\|_2 = \|\hat{f}\|_2

Nói cách khác: đổi hệ cơ sở của một vector space sang một hệ cơ sở trực chuẩn khác thì không thay đổi độ dài của các vectors.

Tập X trong phân tích ở trên không nhất thiết là phải hữu hạn, miễn là bilinear form được chọn tốt là ta có thể xây dựng các hệ cơ sở trực chuẩn cho không gian tuyến tính tương ứng. Việc chọn hệ cơ sở như thế nào tùy thuộc rất nhiều vào ứng dụng đang xét. Ý tưởng chính của Fourier analysis là làm thế nào ta có thể phân tích các đặc tính của f thông qua \hat{f} và ngược lại. Ví dụ: nếu chọn X = \{0,1\}^n thì ta có thể dùng Fourier transform này để phân tích influence of variables on boolean functions (Hastad đã dùng ý tưởng này cho hệ thống PCP của ông).

Lần tới tôi sẽ nói thêm về trường hợp khi X là một nhóm (thường là nhóm Abel) và cách chọn hệ cơ sở trực chuẩn là các characters của nhóm.

Chủ đề: Thuật Toán |

1 lời bình cho bài “Nhân ma trận, DFT, và lý thuyết biểu diễn nhóm (2)”

  1. 1
    trần huy tùng viết:

    Chào bạn! Bạn có thể cho mình biết vài điều về BOMBE của Alan Turing không. Mình đang tìm hiểu ngọn nguồn của máy tính. Rất mong được giúp đỡ.Thank.

Ghi lời bình của bạn: