MathLabs

Giải tích

Giải tích điều hòa

Phân tích hàm và tín hiệu thành các sóng cơ bản, tổng quát hóa chuỗi Fourier sang bối cảnh rộng hơn.

Trực giácTừ sóng đến tần số

Bất kỳ âm thanh nào bạn nghe được — tiếng vĩ cầm, giọng nói, tiếng ồn giao thông — thực chất là một tín hiệu áp suất dao động theo thời gian. Nhưng tai bạn (và một bộ chỉnh âm đồ họa) có thể nhận ra nó được tạo thành từ nhiều tông thuần túy ở các cao độ khác nhau. Giải tích điều hòa là ngành toán học của phép tách đó: nó lấy một hàm ff và viết lại nó thành tổng hoặc tích phân của các dao động thuần túy e2πixξe^{2\pi i x \xi}, mỗi dao động có tần số riêng ξ\xi. Chuỗi Fourier làm điều này cho tín hiệu tuần hoàn bằng một danh sách tần số rời rạc; giải tích điều hòa mở rộng ý tưởng sang tín hiệu không tuần hoàn, nhiều chiều hơn, và thậm chí các nhóm khác (đường tròn, nhóm hữu hạn, nhóm Lie).

Widget chuỗi Fourier tương tác cho thấy một sóng được dựng từ số lượng họa ba thay đổi.
Một sóng gần vuông được dựng lại từ 24 họa ba Fourier đầu tiên; kéo n để thấy hiện tượng 'gợn sóng' Gibbs gần chỗ nhảy bậc thu hẹp về bề rộng nhưng không giảm về chiều cao.

Đại họcChuỗi Fourier trên đường tròn: câu hỏi hội tụ

Viết tổng riêng thứ NN của chuỗi Fourier một hàm tuần hoàn ff là SNf(x)=∑∣n∣≤Nf^(n) einxS_N f(x) = \sum_{|n|\le N} \hat f(n)\, e^{inx}, trong đó f^(n)=12π∫−ππf(x) e−inx dx\hat f(n) = \frac{1}{2\pi}\int_{-\pi}^{\pi} f(x)\, e^{-inx}\,dx. Trong L2L^2 mọi thứ gọn gàng: các hàm mũ einxe^{inx} tạo thành một cơ sở trực chuẩn, đẳng thức Parseval 12π∫∣f∣2=∑∣f^(n)∣2\frac{1}{2\pi}\int |f|^2 = \sum |\hat f(n)|^2 đúng, và SNf→fS_N f \to f theo chuẩn L2L^2. Câu hỏi khó là hội tụ tại từng điểm: liệu SNf(x)→f(x)S_N f(x) \to f(x) tại một điểm xx cụ thể? Sóng vuông ở trên đã hé lộ khó khăn đó — gần chỗ nhảy bậc, SNfS_N f luôn vọt lố khoảng 9% độ cao bước nhảy, dù NN lớn đến đâu (hiện tượng Gibbs, với độ vọt lố tiến về 2πSi(π)≈1.179\frac{2}{\pi}\mathrm{Si}(\pi) \approx 1.179 lần nửa bước nhảy).

Định nghĩa: Hạch Dirichlet

Tổng riêng thực chất là một phép chập: SNf(x)=(DN∗f)(x)=12π∫−ππDN(x−y) f(y) dyS_N f(x) = (D_N * f)(x) = \frac{1}{2\pi}\int_{-\pi}^{\pi} D_N(x-y)\, f(y)\, dy, trong đó hạch Dirichlet là DN(x)=∑∣n∣≤Neinx=sin⁡((N+12)x)sin⁡(x/2)D_N(x) = \sum_{|n|\le N} e^{inx} = \dfrac{\sin\left((N+\tfrac12)x\right)}{\sin(x/2)}. Mọi câu hỏi hội tụ vì thế thực chất là câu hỏi về hình dạng của DND_N.

DN(x)=∑∣n∣≤Neinx=sin⁡((N+12)x)sin⁡(x/2)D_N(x) = \sum_{|n|\le N} e^{inx} = \dfrac{\sin\left((N+\tfrac12)x\right)}{\sin(x/2)}
Biểu đồ tương tác của hạch Dirichlet, Fejér hoặc Poisson trên một chu kỳ, kèm ước lượng trực tiếp chuẩn L1.
Hạch Dirichlet DND_N với N=9N=9: khác với hạch Fejér hay Poisson, nó lõm xuống dưới 0 (các thùy san hô), và chuẩn L1L^1 của nó — hằng số Lebesgue — tăng như (4/π2)log⁡N(4/\pi^2)\log N thay vì bị chặn. Đổi hạch hoặc tăng NN để so sánh.

Hạch Fejér FNF_N (trung bình Cesàro của D0,…,DND_0,\ldots,D_N) và hạch Poisson PrP_r đều không âm, nên chúng là những xấp xỉ đơn vị đúng nghĩa: khối lượng 11, chuẩn L1L^1 bằng 11, và khối lượng dồn về 00 khi N→∞N\to\infty (hay r→1r\to1) — đủ để buộc FN∗f→fF_N * f\to f và Pr∗f→fP_r * f\to f đều với mọi ff liên tục. Hạch Dirichlet thì không: các thùy âm của nó khiến 12π∫∣DN∣\frac{1}{2\pi}\int|D_N| (hằng số Lebesgue) tăng như 4π2log⁡N\frac{4}{\pi^2}\log N, không bị chặn khi N→∞N\to\infty. Theo nguyên lý bị chặn đều Banach–Steinhaus, hằng số Lebesgue không bị chặn buộc phải tồn tại một hàm liên tục có chuỗi Fourier phân kỳ tại một điểm — chính xác là hiện tượng du Bois-Reymond dựng bằng tay năm 1873, ở đây có được miễn phí. Kolmogorov sau đó chỉ ra sự thất bại còn tệ hơn nhiều với hàm chỉ khả tích: năm 1923 ông dựng một hàm L1L^1 có chuỗi Fourier phân kỳ hầu khắp nơi, và năm 1926 một hàm phân kỳ khắp nơi.

Định lý: Định lý Fejér

Nếu ff liên tục và tuần hoàn chu kỳ 2π2\pi, thì các trung bình Cesàro σNf=1N+1∑k=0NSkf=FN∗f\sigma_N f = \frac{1}{N+1}\sum_{k=0}^N S_k f = F_N * f hội tụ đều về ff khi N→∞N\to\infty.

Vì sao đúng?

Lấy trung bình các tổng riêng triệt tiêu chính xác sự dao động khiến SNfS_N f trở nên bất trị: trong khi DND_N có các thùy âm dữ dội, trung bình chạy của nó FNF_N lại là một gò trơn, không âm, chỉ cộng các phần của ff với trọng số dương. Vì thế dù SNfS_N f có thể không hội tụ, trung bình chạy của nó luôn ổn định — cùng mẹo dùng để thuần hóa một dãy tổng riêng dao động trong giải tích bằng cách xét trung bình Cesàro của nó thay vào đó.

Chứng minh

Bước 1 (hạch Fejér là xấp xỉ đơn vị dương). Lấy trung bình D0,…,DND_0,\ldots,D_N và dùng một đẳng thức lượng giác kiểu viễn vọng cho công thức đóng FN(x)=1N+1(sin⁡(N+12x)sin⁡(x/2))2≥0F_N(x) = \frac{1}{N+1}\left(\dfrac{\sin\left(\tfrac{N+1}{2}x\right)}{\sin(x/2)}\right)^2 \ge 0. Vì FN=1N+1∑k=0NDkF_N = \frac{1}{N+1}\sum_{k=0}^N D_k và 12π∫Dk=1\frac{1}{2\pi}\int D_k = 1 với mọi kk, nên 12π∫FN=1\frac{1}{2\pi}\int F_N = 1; do không âm, 12π∫∣FN∣=1\frac{1}{2\pi}\int |F_N| = 1 cũng vậy, nên chuẩn L1L^1 không bao giờ nổ.

Bước 2 (khối lượng dồn về 00). Với δ>0\delta>0 cố định bất kỳ, trên δ≤∣x∣≤π\delta \le |x| \le \pi mẫu số sin⁡(x/2)2\sin(x/2)^2 bị chặn dưới bởi một hằng số dương, nên FN(x)=O(1/N)F_N(x) = O(1/N) đều ở đó; khối lượng ngoài [−δ,δ][-\delta,\delta] tiến về 00 khi N→∞N\to\infty.

Bước 3 (ước lượng xấp xỉ đơn vị). Viết σNf(x)−f(x)=12π∫−ππFN(y)(f(x−y)−f(x)) dy\sigma_N f(x) - f(x) = \frac{1}{2\pi}\int_{-\pi}^{\pi} F_N(y)\big(f(x-y)-f(x)\big)\,dy dùng 12π∫FN=1\frac{1}{2\pi}\int F_N = 1. Chia tích phân thành ∣y∣<δ|y|<\delta và δ≤∣y∣≤π\delta\le|y|\le\pi. Trên phần đầu, tính liên tục đều của ff cho ∣f(x−y)−f(x)∣<ε|f(x-y)-f(x)|<\varepsilon khi δ\delta nhỏ, và FN≥0F_N\ge0 với tổng khối lượng 11 chặn phần đó bởi ε\varepsilon. Trên phần sau, FN=O(1/N)F_N=O(1/N) và ff bị chặn, nên phần đó →0\to 0 khi N→∞N\to\infty. Cả hai chặn đều đều theo xx, nên σNf→f\sigma_N f\to f đều.

Nếu f∈L2(T)f\in L^2(\mathbb{T}), thì SNf(x)→f(x)S_N f(x)\to f(x) với hầu hết xx. (Hunt, 1968, mở rộng kết quả này cho mọi f∈Lp(T)f\in L^p(\mathbb{T}) với p>1p>1.)

Vì sao đúng?

Với ví dụ năm 1926 của Kolmogorov về một hàm L1L^1 phân kỳ khắp nơi, có vẻ hoàn toàn hợp lý khi nghĩ rằng L2L^2 — một điều kiện chỉ mạnh hơn chút ít — cũng sẽ thất bại theo cách tương tự, hoặc ít nhất là hầu khắp nơi. Định lý Carleson là một bất ngờ thực sự: giả thiết bình phương khả tích nhỏ bé đó loại trừ hoàn toàn sự phân kỳ, trừ trên một tập đo được bằng không. Nó khép lại một bài toán đã mở từ khi Luzin đặt giả thuyết năm 1913, và được xem là một trong những định lý sâu sắc nhất của giải tích thế kỷ 20.

Chứng minh

Chứng minh đầy đủ là một trong những lập luận khó nhất của giải tích thế kỷ 20; Fefferman đưa ra một sự đơn giản hóa nổi tiếng năm 1973, có thể phác họa chiến lược qua ba bước.

Bước 1 (kiểm soát bằng toán tử cực đại). Chỉ cần chặn toán tử cực đại Carleson S∗f(x)=sup⁡N∣SNf(x)∣S^*f(x) = \sup_N |S_N f(x)| trên L2L^2, vì một chặn kiểu yếu trên S∗S^* cùng tính trù mật của các hàm tốt (mà với chúng hội tụ là dễ) sẽ suy ra phát biểu hội tụ hầu khắp nơi đầy đủ.

Bước 2 (phân rã thời gian–tần số thành các ô). Phân rã ff bằng các gói sóng thích ứng với các ô nhị phân trên mặt phẳng thời gian–tần số — các hình chữ nhật có diện tích ∼1\sim 1, chính là các nguyên tử được bàn ở mục 'Thời gian và tần số' bên dưới. Mỗi ô mang một phần của ff khu trú cả về một khoảng thời gian lẫn một dải tần số.

Bước 3 (sắp xếp các ô thành cây và cộng lại). Vì tần số cắt NN có thể thay đổi theo xx, các ô liên quan đến S∗f(x)S^*f(x) tạo thành các 'cây' được sắp theo quan hệ chứa thời gian–tần số. Ước lượng tổ hợp then chốt của Fefferman chặn tổng năng lượng mang bởi mọi cây, kiểm soát S∗fS^*f trong L2L^2 và hoàn tất chứng minh.

Kỹ thuật cuối cùng đã phá vỡ định lý Carleson — sắp xếp các ô thời gian–tần số thành cây — hóa ra là một công cụ dùng chung trong giải tích điều hòa hiện đại (chẳng hạn, cùng kiểu lập luận đó đã giải quyết một bài toán mở của Calderón về biến đổi Hilbert song tuyến tính vào năm 1997, bởi Lacey và Thiele).

Đại họcBiến đổi Fourier trên đường thẳng thực

Định nghĩa: Biến đổi Fourier

Với một hàm khả tích ff trên R\mathbb{R}, biến đổi Fourier f^\hat f của nó ghi lại lượng mỗi tần số thuần túy ξ\xi có mặt trong ff. Nó được định nghĩa bởi f^(ξ)=∫−∞∞f(x) e−2πixξ dx\hat f(\xi) = \int_{-\infty}^{\infty} f(x)\, e^{-2\pi i x \xi}\,dx, trong đó lũy thừa e−2πixξe^{-2\pi i x \xi} là một véc-tơ quay có độ dài đơn vị; nhân f(x)f(x) với nó rồi lấy tích phân sẽ đo độ tương quan giữa ff và tốc độ quay đó.

f^(ξ)=∫−∞∞f(x) e−2πixξ dx\hat f(\xi) = \int_{-\infty}^{\infty} f(x)\, e^{-2\pi i x \xi}\,dx

Khi f^\hat f cũng khả tích, tín hiệu ban đầu có thể được dựng lại chính xác bằng cách cộng tất cả các tần số thuần túy đó trở lại: f(x)=∫−∞∞f^(ξ) e2πixξ dξf(x) = \int_{-\infty}^{\infty} \hat f(\xi)\, e^{2\pi i x \xi}\,d\xi. Công thức nghịch đảo này chính là ý nghĩa chính xác của câu 'một hàm bằng tổng các thành phần tần số của nó.'

f(x)=∫−∞∞f^(ξ) e2πixξ dξf(x) = \int_{-\infty}^{\infty} \hat f(\xi)\, e^{2\pi i x \xi}\,d\xi
Ba khía cạnh của cùng một ý tưởng
Bối cảnhMiền của tín hiệuPhía tần sốĐẳng thức chính
Chuỗi FourierĐường tròn / chu kỳSố nguyên Z\mathbb{Z}Tổng rời rạc các họa ba
Biến đổi FourierĐường thẳng thực R\mathbb{R}Đường thẳng thực R\mathbb{R}f(x)=∫−∞∞f^(ξ) e2πixξ dξf(x) = \int_{-\infty}^{\infty} \hat f(\xi)\, e^{2\pi i x \xi}\,d\xi
Tổng PoissonĐường thẳng thực, lấy mẫu tại Z\mathbb{Z}Số nguyên Z\mathbb{Z}∑n=−∞∞f(n)=∑k=−∞∞f^(k)\sum_{n=-\infty}^{\infty} f(n) = \sum_{k=-\infty}^{\infty} \hat f(k)

Nâng caoCác định lý chính

Nếu f∈L1(R)∩L2(R)f \in L^1(\mathbb{R}) \cap L^2(\mathbb{R}), thì ∫−∞∞∣f(x)∣2 dx=∫−∞∞∣f^(ξ)∣2 dξ\int_{-\infty}^{\infty} |f(x)|^2\,dx = \int_{-\infty}^{\infty} |\hat f(\xi)|^2\,d\xi. Nói cách khác: biến đổi Fourier bảo toàn tổng năng lượng, và với chuẩn hóa này không cần thêm hằng số nào.

Vì sao đúng?

Về mặt vật lý, ∫∣f∣2\int |f|^2 là năng lượng của tín hiệu (hãy nghĩ: công suất tiêu tán bởi một dạng sóng điện áp). Plancherel nói rằng bạn có thể tính năng lượng đó bằng cách quét tín hiệu theo thời gian hoặc quét phổ của nó theo tần số — một máy phân tích phổ và một máy hiện sóng phải cho cùng tổng công suất.

Chứng minh

Bước 1 (một trường hợp đặc biệt). Trước tiên kiểm tra đẳng thức cho hàm Gauss f(x)=e−πx2f(x) = e^{-\pi x^2}. Tính toán trực tiếp (bằng cách hoàn thành bình phương trong số mũ) cho ra f^(ξ)=e−πξ2\hat f(\xi) = e^{-\pi \xi^2}, nên cả hai vế của ∫−∞∞∣f(x)∣2 dx=∫−∞∞∣f^(ξ)∣2 dξ\int_{-\infty}^{\infty} |f(x)|^2\,dx = \int_{-\infty}^{\infty} |\hat f(\xi)|^2\,d\xi đều bằng cùng một tích phân Gauss ∫e−2πx2 dx\int e^{-2\pi x^2}\,dx.

Bước 2 (mở rộng bằng tuyến tính và định lý tích chập). Với các hàm 'tốt' (hàm Schwartz) ff, đặt g=f∗f~g = f * \tilde f với f~(x)=f(−x)‾\tilde f(x) = \overline{f(-x)}; khi đó g(0)=∫∣f(x)∣2 dxg(0) = \int |f(x)|^2\,dx, và theo định lý tích chập f∗g^=f^⋅g^\widehat{f * g} = \hat f \cdot \hat g cùng công thức nghịch đảo Fourier, g(0)=∫g^(ξ) dξ=∫∣f^(ξ)∣2 dξg(0) = \int \hat g(\xi)\,d\xi = \int |\hat f(\xi)|^2\,d\xi. Điều này chứng minh đẳng thức cho mọi hàm Schwartz, bao gồm cả hàm Gauss và mọi hàm trơn giảm nhanh.

Bước 3 (lập luận trù mật). Các hàm Schwartz trù mật trong L2(R)L^2(\mathbb{R}): mọi f∈L2f \in L^2 đều là giới hạn fn→ff_n \to f của các hàm Schwartz theo chuẩn L2L^2. Vì biến đổi Fourier là một phép đẳng cự trên không gian con trù mật này theo Bước 2, nó mở rộng duy nhất thành một toán tử bị chặn trên toàn bộ L2(R)L^2(\mathbb{R}) vẫn thỏa mãn ∫−∞∞∣f(x)∣2 dx=∫−∞∞∣f^(ξ)∣2 dξ\int_{-\infty}^{\infty} |f(x)|^2\,dx = \int_{-\infty}^{\infty} |\hat f(\xi)|^2\,d\xi — phần mở rộng này chính là ý nghĩa của 'biến đổi Fourier của một hàm L2L^2' khi tích phân định nghĩa f^(ξ)=∫−∞∞f(x) e−2πixξ dx\hat f(\xi) = \int_{-\infty}^{\infty} f(x)\, e^{-2\pi i x \xi}\,dx có thể không hội tụ tuyệt đối.

Với hàm Schwartz ff, ∑n=−∞∞f(n)=∑k=−∞∞f^(k)\sum_{n=-\infty}^{\infty} f(n) = \sum_{k=-\infty}^{\infty} \hat f(k): tổng của hàm trên mọi số nguyên bằng tổng biến đổi Fourier của nó trên mọi số nguyên.

Vì sao đúng?

Đây là cầu nối giữa việc lấy mẫu tín hiệu tại các điểm nguyên và tuần hoàn hóa phổ của nó: vế trái là kết quả khi cộng các mẫu f(n)f(n), vế phải là kết quả từ phổ. Đẳng thức duy nhất này là nền tảng của định lý lấy mẫu, tổng trên mạng tinh thể học, và phương trình hàm của hàm theta và hàm zeta Riemann.

Chứng minh

Bước 1 (tuần hoàn hóa). Với hàm Schwartz ff, đặt F(x)=∑n=−∞∞f(x+n)F(x) = \sum_{n=-\infty}^{\infty} f(x+n). Do ff giảm nhanh nên tổng này hội tụ tuyệt đối và đều, và FF là hàm trơn có chu kỳ 1, nên nó có chuỗi Fourier riêng trên đường tròn.

Bước 2 (tính các hệ số Fourier của hàm tuần hoàn hóa). Hệ số Fourier thứ kk của FF là ck=∫01F(x)e−2πikx dx=∑n∫01f(x+n)e−2πikx dxc_k = \int_0^1 F(x) e^{-2\pi i k x}\,dx = \sum_{n} \int_0^1 f(x+n) e^{-2\pi i k x}\,dx. Đổi biến y=x+ny = x+n trong mỗi số hạng và dùng tính tuần hoàn của e−2πikxe^{-2\pi i k x} theo nn sẽ ghép các mảnh lại thành một tích phân duy nhất trên toàn R\mathbb{R}: ck=∫−∞∞f(y)e−2πiky dy=f^(k)c_k = \int_{-\infty}^{\infty} f(y) e^{-2\pi i k y}\,dy = \hat f(k).

Bước 3 (tính tại x=0). Vì FF trơn nên chuỗi Fourier của nó hội tụ điểm về chính nó, đặc biệt tại x=0x=0: F(0)=∑kcke0F(0) = \sum_k c_k e^{0}, tức là ∑n=−∞∞f(n)=∑k=−∞∞f^(k)\sum_{n=-\infty}^{\infty} f(n) = \sum_{k=-\infty}^{\infty} \hat f(k). Vế trái đúng bằng F(0)=∑nf(n)F(0) = \sum_n f(n) theo định nghĩa, hoàn tất chứng minh.

Nâng caoBiến đổi Hilbert: tích phân kỳ dị đầu tiên

Lấy tích chập với một hạch chỉ đơn thuần bị chặn và khả tích — như hạch Poisson hay Fejér ở trên — là hiền lành. Một số toán tử quan trọng nhất trong giải tích lại lấy tích chập với một hạch có kỳ dị thực sự tại gốc, không khả tích gần y=xy=x và chỉ được cứu vãn nhờ một sự triệt tiêu tinh tế giữa phần dương và phần âm của nó. Nguyên mẫu là biến đổi Hilbert trên đường thẳng thực, Hf(x)=p.v. 1π∫−∞∞f(y)x−y dyHf(x) = \text{p.v.}\,\frac{1}{\pi}\int_{-\infty}^{\infty}\frac{f(y)}{x-y}\,dy, được hiểu theo nghĩa giá trị chính: kỳ dị tại y=xy=x được cắt bỏ đối xứng trước khi lấy giới hạn. Về phía tần số, phép chập kỳ dị này trở thành một phép nhân đơn giản đến bất ngờ: (Hf^)(ξ)=−i sgn(ξ) f^(ξ)(\widehat{Hf})(\xi) = -i\,\mathrm{sgn}(\xi)\,\hat f(\xi) — biến đổi Hilbert là một phép quay pha 90∘90^\circ thuần túy của mọi tần số, đổi dấu của ii mà không thay đổi độ lớn.

Hf(x)=p.v. 1π∫−∞∞f(y)x−y dyHf(x) = \text{p.v.}\,\frac{1}{\pi}\int_{-\infty}^{\infty}\frac{f(y)}{x-y}\,dy
(Hf^)(ξ)=−i sgn(ξ) f^(ξ)(\widehat{Hf})(\xi) = -i\,\mathrm{sgn}(\xi)\,\hat f(\xi)

Còn có một cách nhìn cổ điển thứ hai: với ff xác định trên đường thẳng thực, đặt u=Py∗fu = P_y * f là thác triển điều hòa của nó lên nửa mặt phẳng trên (tích chập với hạch Poisson ở mục trước), và đặt v=Qy∗fv = Q_y * f là hàm liên hợp điều hòa của nó, dựng từ hạch Poisson liên hợp QyQ_y sao cho F=u+ivF=u+iv là một hàm chỉnh hình theo x+iyx+iy. Khi y→0+y\to 0^+, giá trị biên của vv khôi phục chính xác biến đổi Hilbert của ff. Đây là tổ tiên giải tích phức của toàn bộ lý thuyết: một toán tử tích phân kỳ dị trên biên thực ra là cái bóng của một hàm chỉnh hình đàng hoàng sống ở một chiều cao hơn.

Biến đổi Hilbert HH bị chặn trên Lp(R)L^p(\mathbb{R}) với mọi 1<p<∞1<p<\infty: tồn tại một hằng số CpC_p sao cho ∥Hf∥Lp≤Cp∥f∥Lp\|Hf\|_{L^p} \le C_p\|f\|_{L^p} với mọi ff. (Marcel Riesz, 1927.)

Vì sao đúng?

Trường hợp L2L^2 về cơ bản là miễn phí: hệ số nhân −i sgn(ξ)-i\,\mathrm{sgn}(\xi) có mô-đun đúng bằng 11 hầu khắp nơi, nên HH đúng nghĩa là một phép đẳng cự trên L2L^2 theo Plancherel — một phép quay pha 90∘90^\circ của mọi tần số không làm thay đổi gì về tổng năng lượng. Nhưng mở rộng tính bị chặn sang các pp khác là một bài toán thực sự khác và khó hơn nhiều: không có tương tự nào của đẳng thức Parseval bên ngoài L2L^2, nên không có lối tắt đại số nào khả dụng, và chứng minh phải dựa vào các ước lượng biến thực về kích thước và hình học của tập nơi HfHf lớn.

Chứng minh

Bước 1 (L2L^2, trực tiếp). Theo Plancherel và công thức hệ số nhân, ∥Hf∥L22=∫∣Hf^(ξ)∣2 dξ=∫∣sgn(ξ)∣2∣f^(ξ)∣2 dξ=∥f∥L22\|Hf\|_{L^2}^2 = \int |\widehat{Hf}(\xi)|^2\,d\xi = \int |\mathrm{sgn}(\xi)|^2|\hat f(\xi)|^2\,d\xi = \|f\|_{L^2}^2, dùng ∣sgn(ξ)∣=1|\mathrm{sgn}(\xi)|=1 với ξ≠0\xi\ne0 (một điểm đơn lẻ không ảnh hưởng tới tích phân). Vậy HH là một phép đẳng cự trên L2(R)L^2(\mathbb{R}), đặc biệt bị chặn với C2=1C_2=1.

Bước 2 (mọi 1<p<∞1<p<\infty, tham chiếu chuyển tiếp). Toàn bộ dải số mũ không suy ra từ Bước 1 bằng bất kỳ lập luận mềm nào. Đó là một định lý thực sự của lý thuyết Calderón–Zygmund biến thực: hạch 1/(πx)1/(\pi x) của biến đổi Hilbert thỏa mãn chính xác các điều kiện trơn và triệt tiêu của một hạch Calderón–Zygmund, nên bộ máy tích phân kỳ dị tổng quát trình bày ở mục sau áp dụng được cho nó. Bộ máy đó sinh ra một chặn kiểu yếu (1,1)(1,1) từ chặn L2L^2 đã chứng minh ở Bước 1, và nội suy Marcinkiewicz giữa chặn yếu-(1,1)(1,1) với chặn L2L^2 (rồi một lập luận đối ngẫu cho p>2p>2) hoàn tất chứng minh cho mọi 1<p<∞1<p<\infty — nửa lập luận này được thực hiện đầy đủ ở mục Calderón–Zygmund ngay bên dưới, thay vì lặp lại ở đây.

Ví dụ: Biến đổi Hilbert của hàm chỉ

Đặt f=1[−1,1]f = \mathbf{1}_{[-1,1]}, hàm chỉ của [−1,1][-1,1]: một hàm bị chặn, giá compact, hoàn toàn bình thường. Tính HfHf.

Lời giải

Theo định nghĩa, Hf(x)=1π p.v.∫−11dyx−yHf(x) = \frac{1}{\pi}\,\text{p.v.}\int_{-1}^{1} \frac{dy}{x-y}. Ngoài [−1,1][-1,1] không có kỳ dị nào cần cắt bỏ, và nguyên hàm của 1/(x−y)1/(x-y) theo yy là −log⁡∣x−y∣-\log|x-y|, nên ∫−11dyx−y=[−log⁡∣x−y∣]y=−1y=1=log⁡∣x+1x−1∣\int_{-1}^{1} \frac{dy}{x-y} = \big[-\log|x-y|\big]_{y=-1}^{y=1} = \log\left|\frac{x+1}{x-1}\right| (với xx nằm trong (−1,1)(-1,1), cùng phép tính vẫn đúng theo nghĩa giá trị chính thực sự, vì hai phân kỳ tại y=x−y=x^- và y=x+y=x^+ triệt tiêu nhau). Điều này cho Hf(x)=1πlog⁡∣x+1x−1∣Hf(x) = \frac{1}{\pi}\log\left|\frac{x+1}{x-1}\right|.

Điểm gây ấn tượng là điều xảy ra tại x=±1x=\pm1: bản thân ff hoàn toàn bị chặn ở đó (nó chỉ nhảy từ 11 xuống 00), vậy mà HfHf nổ theo kiểu lôgarit chính xác tại hai điểm đó. Một hàm bị chặn, giá compact, tốt như trong sách giáo khoa bị HH biến thành một hàm có hai kỳ dị thực sự — đây chính là dấu hiệu cục bộ của mọi toán tử tích phân kỳ dị, và chính vì thế việc kiểm soát HH đòi hỏi các ước lượng tinh tế hơn nhiều so với kiểm soát một phép chập với hạch bị chặn tốt.

Nâng caoHàm cực đại Hardy–Littlewood và phân tích Calderón–Zygmund

Để hoàn tất định lý M. Riesz cho một toán tử tích phân kỳ dị tổng quát — không chỉ biến đổi Hilbert — các nhà giải tích cần một công cụ hoàn toàn biến thực, không nói gì tới biến đổi Fourier cả. Hardy và Littlewood đưa nó ra năm 1930: hàm cực đại Mf(x)=sup⁡r>012r∫∣y−x∣<r∣f(y)∣ dyMf(x) = \sup_{r>0} \frac{1}{2r}\int_{|y-x|<r} |f(y)|\,dy, trung bình lớn nhất có thể của ∣f∣|f| trên bất kỳ khoảng nào có tâm tại xx. Đây là một công cụ thô nhưng cực kỳ hiệu quả để kiểm soát mọi quá trình lấy trung bình có thể áp dụng cho ff gần xx, tất cả cùng một lúc.

MM là kiểu yếu (1,1)(1,1): ∣{x:Mf(x)>λ}∣≤Cλ∥f∥1|\{x : Mf(x) > \lambda\}| \le \frac{C}{\lambda}\|f\|_1 với một hằng số phổ dụng CC (luôn có thể lấy C=5C=5). Do đó, bằng nội suy, MM cũng bị chặn trên Lp(R)L^p(\mathbb{R}) với mọi 1<p≤∞1<p\le\infty — nhưng bản thân MM không bao giờ bị chặn trên L1(R)L^1(\mathbb{R}).

Vì sao đúng?

Mf(x)Mf(x) chặn, trong một lần duy nhất, mọi trung bình của ff có thể lấy trên một khoảng quanh xx — trung bình chạy ở mọi thang đo có thể. Kiểm soát đại lượng duy nhất đó hóa ra chính xác là điều cần để chứng minh định lý vi phân Lebesgue (trung bình của ff trên các khoảng co lại quanh xx hội tụ về f(x)f(x) với hầu hết xx): một khi biết MfMf hữu hạn hầu khắp nơi, một lập luận mềm ngắn gọn nâng cấp điều đó thành phát biểu vi phân đầy đủ. Đây là động cơ biến thực đứng sau lý thuyết Calderón–Zygmund bên dưới, đóng vai trò mà đẳng thức Plancherel đã đóng cho ước lượng L2L^2 dễ dàng.

Chứng minh

Bước 1 (một phủ bằng các quả cầu tốt). Cố định λ>0\lambda>0. Với mỗi xx sao cho Mf(x)>λMf(x)>\lambda, theo định nghĩa của cận trên đúng, tồn tại một khoảng BxB_x có tâm tại xx với 1∣Bx∣∫Bx∣f∣>λ\frac{1}{|B_x|}\int_{B_x}|f| > \lambda. Các quả cầu {Bx}\{B_x\} này phủ tập {Mf>λ}\{Mf>\lambda\}.

Bước 2 (bổ đề phủ Vitali 5r5r). Từ bất kỳ họ quả cầu nào có bán kính bị chặn, luôn có thể trích ra một họ con rời nhau đếm được {Bi}\{B_i\} sao cho các quả cầu giãn nở gấp 55 lần {5Bi}\{5B_i\} vẫn phủ hợp của toàn bộ họ ban đầu. Áp dụng điều này cho {Bx}\{B_x\} để được một họ con rời nhau {Bi}\{B_i\} với {Mf>λ}⊆⋃i5Bi\{Mf>\lambda\}\subseteq\bigcup_i 5B_i.

Bước 3 (cộng các mảnh rời nhau). Mỗi BiB_i thỏa λ∣Bi∣<∫Bi∣f∣\lambda|B_i| < \int_{B_i}|f| theo cách dựng. Vì các BiB_i đôi một rời nhau, cộng theo ii cho λ∑i∣Bi∣<∑i∫Bi∣f∣≤∫R∣f∣=∥f∥1\lambda\sum_i|B_i| < \sum_i\int_{B_i}|f| \le \int_{\mathbb{R}}|f| = \|f\|_1, nên ∑i∣Bi∣<∥f∥1/λ\sum_i|B_i| < \|f\|_1/\lambda. Cuối cùng ∣{Mf>λ}∣≤∑i∣5Bi∣=5∑i∣Bi∣<5λ∥f∥1|\{Mf>\lambda\}| \le \sum_i|5B_i| = 5\sum_i|B_i| < \dfrac{5}{\lambda}\|f\|_1, chính xác là ∣{x:Mf(x)>λ}∣≤5λ∥f∥1|\{x : Mf(x) > \lambda\}| \le \frac{5}{\lambda}\|f\|_1.

Định nghĩa: Phân tích Calderón–Zygmund

Với f∈L1(R)f\in L^1(\mathbb{R}) và một độ cao α>0\alpha>0, Calderón và Zygmund (1952) đã chỉ ra cách chia đường thẳng thành các khoảng nhị phân 'dừng' rời nhau, thích ứng với ff và α\alpha: bắt đầu từ một khoảng nhị phân lớn, rồi đệ quy chia đôi nó, dừng lại và giữ một khoảng nhị phân QjQ_j ngay lần đầu tiên trung bình của nó 1∣Qj∣∫Qj∣f∣\frac{1}{|Q_j|}\int_{Q_j}|f| vượt quá α\alpha (trung bình của khoảng cha vẫn còn ≤α\le\alpha, nên chia đôi nhiều nhất chỉ nhân đôi nó). Điều này tạo ra một họ rời nhau các khoảng dừng {Qj}\{Q_j\} với α<1∣Qj∣∫Qj∣f∣≤2α\alpha < \frac{1}{|Q_j|}\int_{Q_j}|f| \le 2\alpha, có tổng độ dài ∑∣Qj∣≤∥f∥1/α\sum|Q_j|\le \|f\|_1/\alpha, và ∣f∣≤α|f|\le\alpha hầu khắp nơi trên phần còn lại ngoài ⋃jQj\bigcup_j Q_j.

Biểu đồ tương tác của phân tích Calderón–Zygmund của một hàm cố định trên [0,1], với thanh trượt mức α có thể điều chỉnh, cho thấy các khoảng dừng, trung bình của chúng, và các đường tham chiếu tại α và 2α.
Một hàm |f| cố định, gập ghềnh trên [0,1] (có một đỉnh hẹp gần x≈0,18, một gò gần kỳ dị gần x≈0,62, và một gò rộng gần x≈0,86), với các khoảng dừng Calderón–Zygmund ở mức này được tô đỏ: mỗi khoảng tô đỏ là một khối dừng Q_j, và đoạn đỏ đậm vẽ ngang qua nó đánh dấu trung bình của nó — luôn nằm kẹp giữa hai đường tham chiếu đứt nét α và 2α. Kéo thanh trượt để đổi α và xem các khoảng dừng co lại, lớn lên, tách ra hoặc gộp lại.

Phép phân tích chia ff thành một phần 'tốt' và vô số mảnh 'xấu', f=g+∑jbjf = g + \sum_j b_j: phần tốt gg bằng ff ở ngoài mọi QjQ_j và bằng giá trị trung bình (cỡ khoảng α\alpha) của ff trên mỗi QjQ_j nơi nó được giữ, nên ∣g∣≤2α|g|\le 2\alpha ở khắp nơi và gg được xử lý bằng lý thuyết L2L^2 dễ dàng. Mỗi mảnh xấu bjb_j có giá trên chính khối QjQ_j của nó và có trung bình bằng không ở đó, ∫Qjbj dx=0\int_{Q_j} b_j\,dx = 0 — sự triệt tiêu đó chính xác là điều khiến đóng góp của bjb_j qua một hạch trơn nhỏ đi một khi đủ xa QjQ_j, vì hạch gần như hằng số trên QjQ_j, và một hạch gần hằng số lấy tích phân một hàm trung bình bằng không cho ra gần như bằng không. Kết hợp phép chia tốt–xấu này với chặn cực đại Hardy–Littlewood kiểu yếu (1,1)(1,1), thông qua nội suy Marcinkiewicz, chính xác là bộ máy biến thực mà Calderón và Zygmund đã dùng để hoàn tất định lý M. Riesz cho các toán tử tích phân kỳ dị tổng quát — không chỉ biến đổi Hilbert — mở rộng tính bị chặn LpL^p sang mọi 1<p<∞1<p<\infty, khép lại vòng tròn trở về định lý ở trên.

Đại họcỨng dụng thực tiễn và Ví dụ minh họa

Vì biến đổi Fourier chuyển 'hình dạng theo thời gian' thành 'nội dung theo tần số,' nó là công cụ tiêu chuẩn ở bất cứ đâu kỹ sư hay nhà khoa học cần cô lập, lọc, hoặc đếm dao động: bộ chỉnh âm, ảnh cộng hưởng từ MRI và kính thiên văn vô tuyến, cùng các phương trình nhiệt và sóng trong vật lý.

Ví dụ: Lọc tiếng rít khỏi một bản ghi âm

Một micrô ghi được f(t)=cos⁡(2π⋅440 t)+cos⁡(2π⋅3000 t)f(t) = \cos(2\pi \cdot 440\,t) + \cos(2\pi \cdot 3000\,t): một nốt nhạc 440 Hz (A4) trộn với tiếng rít điện tử 3000 Hz. Kỹ sư áp dụng một bộ lọc thông thấp lý tưởng đặt f^(ξ)=0\hat f(\xi) = 0 bất cứ khi nào ∣ξ∣>1000|\xi| > 1000. Tín hiệu đầu ra là gì?

Lời giải

Mỗi hàm cos là tổng của hai lũy thừa quay thuần túy ở tần số ±ν\pm\nu: trong miền tần số, f(t)=cos⁡(2π⋅440 t)+cos⁡(2π⋅3000 t)f(t) = \cos(2\pi \cdot 440\,t) + \cos(2\pi \cdot 3000\,t) có năng lượng tập trung chỉ tại ξ=±440\xi = \pm 440 và ξ=±3000\xi = \pm 3000 (như các đỉnh nhọn trong f^\hat f).

Bộ lọc giữ lại mọi thứ với ∣ξ∣≤1000|\xi| \le 1000 và loại bỏ phần còn lại. Vì 440≤1000440 \le 1000 nhưng 3000>10003000 > 1000, các đỉnh tại ±3000\pm 3000 bị loại trong khi các đỉnh tại ±440\pm 440 còn nguyên vẹn.

Áp dụng công thức nghịch đảo cho phần còn lại sẽ dựng lại chính xác các tần số sống sót: đầu ra là cos⁡(2π⋅440 t)\cos(2\pi \cdot 440\,t), nốt nhạc sạch không còn tiếng rít. Đây chính xác là cách núm lọc thông thấp của bộ chỉnh âm hoạt động.

Ví dụ: Vì sao nhiệt lan tỏa: giải phương trình nhiệt

Một thanh kim loại có phân bố nhiệt độ ban đầu u(x,0)=f(x)u(x,0) = f(x) và tuân theo phương trình nhiệt ut=uxxu_t = u_{xx}. Dùng biến đổi Fourier (theo xx) để tìm u(x,t)u(x,t) với t>0t>0.

Lời giải

Lấy biến đổi Fourier theo xx biến mỗi đạo hàm không gian ∂x\partial_x thành phép nhân với 2πiξ2\pi i \xi, nên uxxu_{xx} trở thành −4π2ξ2u^-4\pi^2\xi^2 \hat u. Phương trình ut=uxxu_t = u_{xx} trở thành phương trình vi phân thường ∂tu^(ξ,t)=−4π2ξ2 u^(ξ,t)\partial_t \hat u(\xi,t) = -4\pi^2\xi^2\, \hat u(\xi,t) theo tt với mỗi ξ\xi cố định — một sự đơn giản hóa lớn, vì một phương trình đạo hàm riêng khó đã trở thành một họ phương trình vi phân thường tách rời dễ giải.

Phương trình vi phân này có nghiệm u^(ξ,t)=f^(ξ) e−4π2ξ2t\hat u(\xi,t) = \hat f(\xi)\, e^{-4\pi^2 \xi^2 t}, một thừa số suy giảm theo hàm mũ triệt tiêu nhanh các tần số cao ξ\xi: chi tiết không gian mịn bị làm mượt nhanh chóng, khớp với thực tế hàng ngày rằng các đỉnh nhiệt độ nhọn bị làm phẳng trước tiên.

Nghịch đảo biến đổi (một tích trong miền tần số là một phép chập trong miền không gian, theo định lý tích chập áp dụng ngược lại) cho u(x,t)=(f∗Kt)(x),  Kt(x)=14πt e−x2/(4t)u(x,t) = (f * K_t)(x),\ \ K_t(x) = \dfrac{1}{\sqrt{4\pi t}}\, e^{-x^2/(4t)}: nhiệt độ tại thời điểm tt là phân bố ban đầu bị làm nhòe (chập) với một bướu Gauss lan rộng KtK_t — đúng nghĩa đen là bức tranh toán học của sự khuếch tán nhiệt.

Định nghĩa: Nguyên lý bất định

Một hàm không thể được định vị sắc nét đồng thời cả về thời gian lẫn tần số. Nếu Δx\Delta x và Δξ\Delta \xi lần lượt đo độ trải rộng của ∣f∣2|f|^2 và ∣f^∣2|\hat f|^2, thì Δx⋅Δξ≥14π\Delta x \cdot \Delta \xi \ge \dfrac{1}{4\pi}. Nén một xung theo thời gian (nhỏ Δx\Delta x) buộc phổ của nó phải trải rộng ra (lớn Δξ\Delta \xi), và ngược lại — cùng sự đánh đổi được dùng để suy ra nguyên lý bất định Heisenberg trong cơ học lượng tử.

Δx⋅Δξ≥14π\Delta x \cdot \Delta \xi \ge \dfrac{1}{4\pi}

Với quy ước f^(ξ)=∫f(x)e−2πixξ dx\hat f(\xi) = \int f(x) e^{-2\pi i x\xi}\,dx, biến đổi Fourier của hàm Gauss f(x)=e−πx2f(x) = e^{-\pi x^2} là:

Định lý Plancherel, ∫−∞∞∣f(x)∣2 dx=∫−∞∞∣f^(ξ)∣2 dξ\int_{-\infty}^{\infty} |f(x)|^2\,dx = \int_{-\infty}^{\infty} |\hat f(\xi)|^2\,d\xi, được mô tả đúng nhất như một phát biểu về:

Một máy MRI đo các mẫu của biến đổi Fourier của tín hiệu không gian ('k-space'). Nếu chỉ thu thập các mẫu tần số thấp ( ∣ξ∣|\xi| nhỏ), ảnh dựng lại sẽ:

Nguyên lý bất định Δx⋅Δξ≥14π\Delta x \cdot \Delta \xi \ge \dfrac{1}{4\pi} ngụ ý rằng:

Định lý Fejér nói rằng các trung bình Cesàro σNf\sigma_N f hội tụ đều về mọi hàm liên tục ff. Định lý Carleson là một phát biểu riêng biệt, sâu sắc hơn nhiều, về:

Theo định lý M. Riesz, biến đổi Hilbert HH bị chặn trên Lp(R)L^p(\mathbb{R}) chính xác với:

Tài liệu tham khảo

  1. Terence Tao (2003). Recent progress on the restriction conjecture · arXiv:math/0311181
  2. Elias M. Stein, Guido Weiss (1971). Introduction to Fourier Analysis on Euclidean Spaces
  3. Loukas Grafakos (2014). Classical Fourier Analysis · DOI:10.1007/978-1-4939-1194-3