Sóng Morlet khác gì so với biến đổi Fourier truyền thống?

Biến đổi Fourier truyền thống chỉ phân tích tín hiệu trên miền tần số thuần túy và làm mất hoàn toàn thông tin về thời điểm xảy ra biến cố. Ngược lại, sóng Morlet sử dụng cửa sổ Gauss co giãn theo tần số, cho phép xác định chính xác thời điểm và dải tần số của các hiện tượng phi dừng.

Tại sao tần số trung tâm của sóng Morlet thường được chọn bằng 6?

Khi tần số trung tâm bằng 6, giá trị phổ Fourier tại điểm 0 của sóng Morlet chỉ xấp xỉ 1.5e-8, đủ nhỏ để triệt tiêu trên thực tế tính toán số và thỏa mãn điều kiện chấp nhận wavelet mà không cần bổ sung số hạng hiệu chỉnh phức tạp.

Sóng Morlet là sóng thực hay sóng phức?

Sóng Morlet chuẩn là một hàm phức bao gồm phần thực cosin và phần ảo sin nhân với bao Gauss. Cấu trúc phức này cho phép người phân tích tách biệt và đo lường trực tiếp cả biên độ phổ tức thời lẫn góc pha dao động cục bộ của tín hiệu.

Khi nào không nên sử dụng sóng Morlet?

Sóng Morlet không phù hợp cho các bài toán yêu cầu nén dữ liệu hoặc biến đổi trực giao không dư thừa (như chuẩn nén ảnh JPEG 2000), do phép biến đổi liên tục CWT của sóng Morlet tạo ra lượng dữ liệu hệ số lớn và có tính dư thừa thông tin cao.

Sóng Morlet (Morlet wavelet) là gì?

Sóng Morlet là một hàm sóng con liên tục (continuous wavelet) thuộc lớp hàm phức, được xây dựng từ một hàm sóng mang điều hòa nhân với một hàm bao hình Gauss nhằm phân tích các tín hiệu phi dừng trong không gian thời gian - tần số. Nhờ cấu trúc giải tích đối xứng, sóng Morlet mang lại sự cân bằng tối ưu giữa độ phân giải thời gian và độ phân giải tần số, cho phép trích xuất đồng thời thông tin về biên độ dao động tức thời và pha của các thành phần dao động cục bộ. Mục từ này trình bày định nghĩa toán học, nguồn gốc lịch sử, nguyên lý biến đổi liên tục, tính chất định vị theo nguyên lý bất định Heisenberg, cùng các ứng dụng tiêu biểu và giới hạn kỹ thuật của sóng Morlet.

Khái niệm và nguồn gốc lịch sử

Trong phân tích chuỗi thời gian và tín hiệu vật lý, các công cụ cổ điển như biến đổi Fourier toàn cục bộc lộ hạn chế căn bản khi xử lý các tín hiệu có đặc tính biến đổi theo thời gian (tín hiệu phi dừng). Biến đổi Fourier trải rộng các hàm cơ sở hình sin và cosin trên toàn bộ trục thời gian vô hạn, do đó làm mất hoàn toàn thông tin về thời điểm xuất hiện của các biến cố hoặc đột biến tần số. Nhằm khắc phục hạn chế này, biến đổi Fourier thời gian ngắn (Short-Time Fourier Transform, STFT) do Dennis Gabor đề xuất đã áp dụng một cửa sổ thời gian cố định trượt dọc theo tín hiệu. Tuy nhiên, chiều rộng cửa sổ cố định của STFT tạo ra một độ phân giải bất biến, không thể đồng thời phân tích tốt các thành phần tần số cao xảy ra chớp nhoáng và các thành phần tần số thấp biến đổi chậm.

Vào năm 1982, kỹ sư địa vật lý người Pháp Jean Morlet công tác tại tập đoàn dầu khí Elf Aquitaine đã đề xuất ý tưởng sử dụng các hàm dao động có độ rộng co giãn tự nhiên tỷ lệ nghịch với tần số, gọi là các sóng con có hình dạng không đổi (ondelettes de forme constante), phục vụ xử lý sóng địa chấn phản xạ trong thăm dò địa tầng. Ngay sau đó, Jean Morlet hợp tác với nhà vật lý lý thuyết Alex Grossmann để hoàn thiện nền tảng toán học giải tích hàm trên không gian Hilbert. Đến năm 1984, Grossmann và Morlet chính thức công bố công trình chuyên khảo đặt nền móng giải tích cho biến đổi wavelet liên tục. Trong cùng năm 1984, Goupillaud, Grossmann và Morlet tiếp tục công bố các thuật toán chu kỳ bát độ (cycle-octave) ứng dụng trực tiếp sóng Morlet vào phân tích dữ liệu địa chấn phản xạ, đánh dấu sự ra đời của một nhánh toán học ứng dụng hiện đại.

Về mặt phân loại học thuật, sóng Morlet nguyên bản là một wavelet dạng phức (complex Morlet wavelet). Trong một số bài toán kỹ thuật chỉ yêu cầu phân tích thành phần đối xứng, phần thực của hàm (dạng sóng cosin nhân bao Gauss) được trích xuất riêng biệt và thường được gọi là sóng Morlet thực (real Morlet wavelet) hoặc sóng Gabor thực. Tuy nhiên, dạng phức vẫn là cấu trúc tiêu chuẩn và mang lại giá trị phân tích vượt trội nhờ tính chất giải tích tự nhiên.

Cấu trúc toán học và biểu diễn giải tích

Về mặt toán học, hàm sóng con mẹ Morlet dạng phức được định nghĩa trên trục thời gian liên tục theo biểu thức sau:

\psi(t) = \pi^{-1/4} e^{i \omega_0 t} e^{-t^2/2}

Trong biểu thức trên, các đại lượng được định nghĩa cụ thể như sau:

Biểu thức $t$ đại diện cho biến số thời gian không thứ nguyên.
Hệ số $\pi^{-1/4}$ là thừa số chuẩn hóa năng lượng, đảm bảo năng lượng toàn phần của hàm sóng con trên toàn trục số thực đạt giá trị đơn vị, tức tích phân bình phương độ lớn của hàm bằng 1.
Thành phần $e^{i \omega_0 t}$ là hàm sóng mang điều hòa phức, dao động với tần số trung tâm không thứ nguyên $\omega_0$ . Theo công thức Euler, thành phần này được khai triển thành $\cos(\omega_0 t) + i \sin(\omega_0 t)$ .
Hàm $e^{-t^2/2}$ là hàm bao hình Gauss có độ lệch chuẩn đơn vị, đóng vai trò cửa sổ suy giảm trơn tru giúp cô lập năng lượng dao động trong một khoảng thời gian hữu hạn xung quanh gốc tọa độ.

Khai triển chi tiết theo dạng đại số phức tách biệt, sóng Morlet được phân rã thành hai phần trực giao:

\text{Re}\{\psi(t)\} = \pi^{-1/4} \cos(\omega_0 t) e^{-t^2/2}

\text{Im}\{\psi(t)\} = \pi^{-1/4} \sin(\omega_0 t) e^{-t^2/2}

Trong đó, phần thực $\text{Re}\{\psi(t)\}$ là một hàm chẵn đối xứng qua gốc tọa độ, đóng vai trò bộ lọc đồng pha với các dao động đỉnh; phần ảo $\text{Im}\{\psi(t)\}$ là một hàm lẻ phản đối xứng qua gốc tọa độ, đóng vai trò bộ lọc lệch pha một góc vuông. Sự kết hợp trực giao này tạo nên tính chất giải tích (analytic signal), giúp hàm wavelet không bị triệt tiêu biên độ tại các điểm cắt không của tín hiệu khảo sát.

Phổ tần số của sóng Morlet được xác định thông qua phép biến đổi Fourier thuận:

\hat{\psi}(\omega) = \int_{-\infty}^{\infty} \psi(t) e^{-i \omega t} dt = \pi^{-1/4} \sqrt{2\pi} e^{-(\omega - \omega_0)^2 / 2}

Phổ Fourier biểu diễn một đường cong phân phối Gauss tập trung đối xứng quanh tần số trung tâm $\omega = \omega_0$ . Phổ này có tính chất đơn cực (chỉ tập trung chủ yếu ở miền tần số dương khi $\omega_0$ đủ lớn), giúp loại bỏ hiện tượng giao thoa giữa các thành phần tần số âm và tần số dương khi tái cấu trúc tín hiệu.

Nguyên lý biến đổi wavelet liên tục và điều kiện chấp nhận

Trong khuôn khổ phép biến đổi wavelet liên tục (Continuous Wavelet Transform, CWT), từ hàm sóng con mẹ $\psi(t)$ , một họ các sóng con con được sinh ra thông qua hai phép biến đổi hình học cơ bản là co giãn tỷ lệ (scaling) và dịch chuyển vị trí thời gian (translation):

\psi_{a,b}(t) = \frac{1}{\sqrt{|a|}} \psi\left(\frac{t - b}{a}\right)

Trong đó, tham số $a$ là hệ số co giãn biểu thị tỷ lệ thời gian (scale), tỷ lệ nghịch với tần số vật lý; tham số $b$ là hệ số dịch chuyển vị trí xác định tọa độ thời gian cục bộ cần khảo sát. Hệ số chuẩn hóa $\frac{1}{\sqrt{|a|}}$ bảo toàn mức năng lượng không đổi của sóng con qua mọi cấp tỷ lệ.

Hệ số biến đổi wavelet liên tục của một tín hiệu thực tế $x(t)$ được tính bằng tích phân xen phủ:

W(a, b) = \int_{-\infty}^{\infty} x(t) \psi_{a,b}^*(t) dt = \frac{1}{\sqrt{|a|}} \int_{-\infty}^{\infty} x(t) \psi^*\left(\frac{t - b}{a}\right) dt

Với $\psi^*$ là liên hợp phức của hàm sóng con. Vì sóng Morlet là hàm phức, hệ số wavelet thu được $W(a, b)$ cũng là một đại lượng phức. Từ đó, người phân tích dễ dàng tách biệt hai thông số vật lý quan trọng:

Biên độ năng lượng phổ thời gian - tần số: $|W(a, b)| = \sqrt{\text{Re}\{W(a, b)\}^2 + \text{Im}\{W(a, b)\}^2}$ , phản ánh cường độ tức thời của dao động tại vị trí thời gian $b$ và tỷ lệ $a$ .
Góc pha tức thời: $\theta(a, b) = \arctan\left(\frac{\text{Im}\{W(a, b)\}}{\text{Re}\{W(a, b)\}}\right)$ , phản ánh trạng thái dao động cục bộ và là cơ sở tính toán tính đồng bộ pha giữa các kênh tín hiệu.

Điều kiện chấp nhận toán học (Admissibility Condition)

Để một hàm sóng có thể dùng làm cơ sở cho biến đổi liên tục và đảm bảo khả năng khôi phục hoàn toàn tín hiệu gốc thông qua phép biến đổi ngược, hàm sóng phải thỏa mãn điều kiện chấp nhận giải tích:

C_\psi = \int_{-\infty}^{\infty} \frac{|\hat{\psi}(\omega)|^2}{|\omega|} d\omega < \infty

Điều kiện tích phân hữu hạn này đòi hỏi phổ Fourier của hàm sóng con tại điểm tần số bằng 0 phải triệt tiêu hoàn toàn, tức $\hat{\psi}(0) = 0$ , đồng nghĩa với việc diện tích đại số của hàm sóng con trên toàn bộ trục thời gian phải bằng 0:

\int_{-\infty}^{\infty} \psi(t) dt = 0

Xét theo định nghĩa giải tích của hàm Morlet chuẩn, giá trị phổ Fourier tại gốc tần số là:

\hat{\psi}(0) = \pi^{-1/4} \sqrt{2\pi} e^{-\omega_0^2 / 2}

Rõ ràng, đại lượng này không hoàn toàn bằng 0 tuyệt đối đối với một giá trị $\omega_0$ hữu hạn. Tuy nhiên, trong nghiên cứu thực hành chuẩn mực do Torrence và Compo công bố năm 1998, các tác giả đã chỉ ra rằng khi lựa chọn tham số tần số trung tâm bằng 6, giá trị $e^{-\omega_0^2 / 2} = e^{-18}$ chỉ đạt xấp xỉ 1.5e-8. Đây là một giá trị vô cùng nhỏ, nằm sâu dưới ngưỡng sai số dấu phẩy động của máy tính kỹ thuật số, do đó điều kiện chấp nhận được coi là thỏa mãn gần như tuyệt đối trên thực tế tính toán.

Đối với các bài toán yêu cầu tính toán giải tích chặt chẽ tuyệt đối, người ta sử dụng dạng sóng Morlet hiệu chỉnh (corrected Morlet wavelet) bằng cách trừ đi số hạng điều chỉnh dịch chuyển thành phần một chiều:

\psi_{\text{adj}}(t) = \pi^{-1/4} \left( e^{i \omega_0 t} - e^{-\omega_0^2/2} \right) e^{-t^2/2}

Khi đó, tích phân toàn phần của hàm trên trục số thực đảm bảo triệt tiêu hoàn toàn về 0 với mọi giá trị của tham số $\omega_0$ .

Mối quan hệ giữa tỷ lệ co giãn và chu kỳ Fourier tương đương

Trong biến đổi wavelet, trục tỷ lệ $a$ không đồng nhất trực tiếp với chu kỳ dao động hình sin truyền thống. Để phiên giải kết quả sang miền tần số vật lý quen thuộc, cần xác định bước sóng Fourier tương đương $\lambda$ . Theo công thức dẫn xuất của Torrence và Compo (1998), chu kỳ Fourier tương đương của sóng Morlet được liên hệ với tỷ lệ co giãn $a$ theo hệ thức:

\lambda = \frac{4\pi a}{\omega_0 + \sqrt{2 + \omega_0^2}}

Khi thiết lập tham số chuẩn $\omega_0 = 6$ , tỷ số giữa chu kỳ Fourier và tỷ lệ co giãn xấp xỉ bằng 1.03. Tỷ số gần bằng 1 này mang lại ưu điểm trực quan đặc biệt lớn cho sóng Morlet: giá trị tỷ lệ $a$ có thể xem gần như tương đương trực tiếp với chu kỳ dao động Fourier của hiện tượng vật lý đang khảo sát mà không đòi hỏi các phép hiệu chỉnh phức tạp.

Tính chất định vị và nguyên lý bất định Heisenberg-Gabor

Mọi công cụ biến đổi thời gian - tần số đều chịu sự ràng buộc cơ bản của nguyên lý bất định Heisenberg-Gabor. Nguyên lý này phát biểu rằng không có một phép biến đổi tuyến tính nào có thể đồng thời định vị chính xác vô hạn cả thời điểm xuất hiện và thành phần tần số của một tín hiệu. Độ bất định thời gian $\Delta t$ và độ bất định tần số góc $\Delta \omega$ luôn bị chặn dưới bởi một hằng số vật lý:

\Delta t \cdot \Delta \omega \ge 0.5

Trong công trình tổng quan kinh điển năm 1990 trên tạp chí IEEE Transactions on Information Theory, nhà toán học Ingrid Daubechies đã chứng minh rằng sóng Morlet là một trong số rất ít các họ hàm cơ sở đạt được giới hạn tối ưu này. Do được bao phủ bởi hàm mũ Gauss (hàm tối ưu duy nhất trong giải tích điều hòa đạt dấu bằng của bất đẳng thức Heisenberg), sóng Morlet thỏa mãn chính xác:

\Delta t \cdot \Delta \omega = 0.5

Điều này có nghĩa là diện tích hộp bất định (uncertainty box) trong mặt phẳng thời gian - tần số của sóng Morlet đạt mức nhỏ nhất có thể về mặt lý thuyết. Khi dịch chuyển sang các dải tỷ lệ khác nhau:

Ở vùng tỷ lệ $a$ nhỏ (tương ứng vùng tần số cao), độ rộng cửa sổ thời gian co hẹp lại, cung cấp khả năng xác định thời điểm xảy ra biến cố cực kỳ chuẩn xác, chấp nhận độ bất định tần số mở rộng.
Ở vùng tỷ lệ $a$ lớn (tương ứng vùng tần số thấp), độ rộng cửa sổ thời gian giãn rộng ra, cung cấp khả năng nhận biết tần số dao động cực kỳ tinh tế, chấp nhận độ bất định thời gian tăng lên.

Cơ chế co giãn tự thích nghi này phản ánh đúng bản chất vật lý của các quá trình tự nhiên: dao động tần số cao thường diễn ra nhanh và mang tính chuyển tiếp, trong khi dao động tần số thấp cần nhiều chu kỳ kéo dài để thiết lập trạng thái ổn định.

So sánh sóng Morlet với các họ wavelet phổ biến

Để lựa chọn công cụ phù hợp trong thực tế nghiên cứu, bảng đối chiếu dưới đây so sánh các đặc tính kỹ thuật cốt lõi giữa sóng Morlet và các họ sóng con thông dụng khác:

Họ sóng con	Dạng đại số	Tính chất trực giao	Khả năng tách pha	Lĩnh vực ứng dụng ưu tiên
Morlet	Phức (sóng mang điều hòa nhân bao Gauss)	Không trực giao	Trực tiếp qua hàm giải tích	Phân tích phổ thời gian - tần số, trích xuất pha, tính đồng pha trong chuỗi thời gian
Mexican Hat	Thực (đạo hàm bậc hai của hàm Gauss)	Không trực giao	Không trích xuất được góc pha	Phát hiện biên cạnh, nhận diện điểm uốn và cực trị cục bộ trong tín hiệu địa chấn
Haar	Thực (hàm bước nhảy gián đoạn)	Trực giao hoàn toàn	Không áp dụng	Phát hiện đột biến bước nhảy, biến đổi rời rạc tốc độ cao, xử lý dữ liệu logic
Daubechies	Thực (hàm trơn compact có số mômen triệt tiêu xác định)	Trực giao trực tiếp	Không áp dụng	Nén dữ liệu số, khử nhiễu tín hiệu, biến đổi wavelet rời rạc đa cấp độ phân giải

Qua bảng so sánh trên, có thể thấy sóng Morlet vượt trội trong các bài toán đo lường và khảo sát phổ liên tục nhờ cấu trúc trơn và khả năng tách pha độc lập. Ngược lại, đối với các ứng dụng yêu cầu nén dữ liệu hoặc hoàn nguyên tín hiệu không dư thừa qua thuật toán phân tích đa phân giải tháp lọc Mallat, các họ wavelet trực giao như Daubechies lại chiếm ưu thế tuyệt đối.

Ứng dụng thực tế trong khoa học và kỹ thuật

Nhờ đặc tính cân bằng thời gian - tần số và biểu diễn pha giải tích, sóng Morlet đã trở thành tiêu chuẩn vàng trong nhiều nhánh khoa học thực nghiệm:

1. Xử lý tín hiệu địa chấn và thăm dò địa vật lý

Kể từ các công trình khởi xướng năm 1984 của Goupillaud, Grossmann và Morlet, biến đổi sóng con Morlet đã được áp dụng rộng rãi để bóc tách các sóng địa chấn phản xạ từ các mặt ranh giới địa tầng ngầm dưới lòng đất. Trong môi trường đất đá có tính hấp thụ và tán xạ tần số, các xung sóng truyền qua bị suy giảm phi tuyến. Biểu đồ thời gian - tần số của sóng Morlet cho phép các nhà địa chấn học định lượng chính xác hệ số suy giảm hấp thụ tầng lọc, phát hiện dị thường túi khí hoặc tích tụ dầu mỏ trong các tầng cát kết sâu.

2. Điện não đồ và khoa học thần kinh nhận thức

Trong sinh lý học thần kinh, não bộ hoạt động thông qua sự phối hợp nhịp nhàng của các mạng lưới neuron phức tạp, biểu hiện qua các dải sóng não đặc trưng. Năm 1996, nghiên cứu bước ngoặt của Tallon-Baudry và các cộng sự công bố trên The Journal of Neuroscience đã áp dụng sóng Morlet phức để phân tách hoạt động dao động trong dải tần gamma 40 Hz từ các bản ghi điện não đồ (EEG) của người. Nhờ khả năng phân tách biên độ và pha độc lập, nhóm tác giả đã chứng minh sự tồn tại của hai cơ chế thần kinh riêng biệt: dao động khóa pha với kích thích thị giác và dao động không khóa pha phản ánh quá trình liên kết nhận thức đối tượng trong vỏ não. Ngày nay, sóng Morlet là công cụ trung tâm trong các phần mềm xử lý tín hiệu thần kinh mã nguồn mở như EEGLAB và Brainstorm để tính toán hệ số đồng bộ pha liên vùng não.

3. Khí tượng thủy văn và nghiên cứu biến đổi khí hậu

Các chuỗi thời gian khí hậu, chẳng hạn như chỉ số dao động phương Nam (SOI) hay chỉ số dao động Bắc Đại Tây Dương (NAO), thường thể hiện các chu kỳ dao động không liên tục và chịu ảnh hưởng của các hiện tượng cực đoan ngẫu nhiên. Công trình kinh điển năm 1998 của Torrence và Compo đã sử dụng sóng Morlet để phân tích hiện tượng El Niño - Dao động phương Nam (ENSO), chứng minh rằng chu kỳ thống trị của ENSO dịch chuyển biến thiên từ khoảng chu kỳ ngắn hạn đến dài hạn, đồng thời định lượng hóa sự tương tác pha giữa các vùng biển xích đạo.

4. Phân tích chuỗi thời gian kinh tế và tài chính tại Việt Nam

Tại Việt Nam, các mô hình phân tích liên kết sóng con (Wavelet Coherency) và liên kết từng phần dựa trên hàm cơ sở Morlet đang được các nhà nghiên cứu tài chính ứng dụng mạnh mẽ trên các ấn phẩm học thuật. Phương pháp này giúp đo lường mức độ đồng chuyển động (co-movement) và hiện tượng lan truyền rủi ro giữa chỉ số thị trường chứng khoán VN-Index, biến động tỷ giá hối đoái tiền tệ và giá vàng quốc tế theo từng khung thời gian ngắn hạn, trung hạn và dài hạn.

Hạn chế và lưu ý kỹ thuật khi áp dụng

Bên cạnh các ưu thế vượt trội, việc triển khai sóng Morlet trong thực nghiệm đòi hỏi người nghiên cứu phải kiểm soát chặt chẽ các giới hạn kỹ thuật sau:

Vùng ảnh hưởng của biên (Cone of Influence, COI): Khi tính toán biến đổi wavelet liên tục trên một chuỗi tín hiệu có độ dài thời gian hữu hạn, các phép tính tích phân xen phủ tại điểm đầu và điểm cuối của dãy số liệu bị thiếu dữ liệu quá khứ hoặc tương lai, dẫn tới hiện tượng méo biên (edge effects). Vùng hình nón ảnh hưởng COI xác định phạm vi mà tại đó biên độ sóng con suy giảm xuống $e^{-2}$ so với đỉnh cực đại; mọi hiện tượng hoặc công suất phổ quan sát được nằm ngoài vùng COI đều không đảm bảo ý nghĩa thống kê và cần bị loại bỏ khi diễn giải dữ liệu.
Tính dư thừa thông tin (Information Redundancy): Biến đổi wavelet liên tục ánh xạ một tín hiệu một chiều thành một ma trận hai chiều liên tục theo thời gian và tỷ lệ. Do sóng Morlet không tạo thành một hệ cơ sở trực giao chặt chẽ, các hệ số wavelet lân cận có mức độ tương quan cao, dẫn tới dung lượng dữ liệu phình to và đòi hỏi chi phí lưu trữ lớn.
Độ trễ tính toán: Việc tính toán tích phân xen phủ qua hàng trăm cấp tỷ lệ co giãn đòi hỏi tài nguyên xử lý lớn hơn nhiều so với thuật toán biến đổi Fourier nhanh (FFT) hay biến đổi wavelet rời rạc (DWT). Trong thực hành tính toán, phép biến đổi CWT của sóng Morlet thường được thực thi gián tiếp trong miền tần số thông qua tích của phép biến đổi Fourier nhanh để tối ưu hóa tốc độ xử lý.