Math Labs

Giải tích

Hàm đặc biệt

Hàm Gamma, Beta, zeta và Bessel mở rộng các hàm sơ cấp, xuất hiện khắp toán học và vật lý.

Trực giácVượt ra ngoài giai thừa: lấp đầy khoảng trống giữa các số nguyên

Giai thừa n!=1⋅2⋯nn! = 1\cdot 2 \cdots n chỉ có nghĩa với các số nguyên, vậy mà nó tăng trưởng mượt mà đến mức — 1,2,6,24,120,…1, 2, 6, 24, 120, \dots — dường như đòi hỏi một đường cong đi qua từng điểm đó, xác định tại mọi đầu vào thực hay thậm chí phức. Hàm đặc biệt chính là điều đó: một họ nhỏ các hàm, sinh ra từ những tích phân hay tổng vô hạn tự nhiên, nội suy và tổng quát hóa các quy luật số học sơ cấp (giai thừa, lũy thừa, tổng) thành một cỗ máy mượt mà, phổ quát, được xác suất, lý thuyết số và vật lý cùng chia sẻ.

Đồ tô màu miền của hàm mũ phức, với các dải màu lặp lại theo chiều dọc và các vòng độ sáng lan ra từ bên trái.
Tô màu miền của hàm mũ eze^z: màu sắc mã hóa argument, các vòng độ sáng mã hóa module. Ánh xạ này nguyên hàm và không bao giờ bằng 0 — một ví dụ mẫu về thác triển giải tích không có cực điểm, khác với hàm Gamma mà nó giúp định nghĩa (Γ(z)Γ(1−z)=π/sin⁡(πz)\Gamma(z)\Gamma(1-z)=\pi/\sin(\pi z)).

Đại họcHàm Gamma, Beta và zeta

Định nghĩa: Hàm Gamma

Với Re⁡(z)>0\operatorname{Re}(z) > 0, định nghĩa Γ(z)=∫0∞tz−1e−t dt\Gamma(z)=\displaystyle\int_0^\infty t^{z-1}e^{-t}\,dt. Lấy tích phân từng phần cho thấy Γ(z+1)=z Γ(z)\Gamma(z+1) = z\,\Gamma(z), và vì Γ(1)=∫0∞e−t dt=1\Gamma(1) = \int_0^\infty e^{-t}\,dt = 1, quy nạp cho Γ(n+1)=n!\Gamma(n+1)=n! với mọi số nguyên không âm nn: hàm Gamma nội suy giai thừa. Công thức truy hồi cũng mở rộng Γ\Gamma bằng thác triển phân hình tới toàn bộ zz trừ z=0,−1,−2,…z = 0, -1, -2, \dots, nơi nó có các cực điểm đơn.

Γ(z)=∫0∞tz−1e−t dt\Gamma(z)=\int_0^\infty t^{z-1}e^{-t}\,dt

Hàm Beta B(x,y)=∫01tx−1(1−t)y−1 dtB(x,y) = \displaystyle\int_0^1 t^{x-1}(1-t)^{y-1}\,dt gói gọn một tỉ số các hàm Gamma, B(x,y)=Γ(x)Γ(y)Γ(x+y)B(x,y) = \dfrac{\Gamma(x)\Gamma(y)}{\Gamma(x+y)}, và là hằng số chuẩn hóa tự nhiên đứng sau phân phối xác suất Beta. Một đẳng thức sâu hơn liên kết hai hàm này một cách đối xứng: công thức phản xạ Euler Γ(z)Γ(1−z)=πsin⁡(πz)\Gamma(z)\Gamma(1-z) = \dfrac{\pi}{\sin(\pi z)}, đúng với mọi zz không nguyên, gắn các cực điểm của Γ(z)\Gamma(z) với các không điểm của sin⁡(πz)\sin(\pi z).

Γ(z)Γ(1−z)=πsin⁡(πz)\Gamma(z)\Gamma(1-z) = \dfrac{\pi}{\sin(\pi z)}

Định nghĩa: Hàm zeta Riemann và thác triển giải tích của nó

Với Re⁡(s)>1\operatorname{Re}(s) > 1, định nghĩa ζ(s)=∑n=1∞1ns\zeta(s) = \displaystyle\sum_{n=1}^{\infty} \dfrac{1}{n^s}. Chuỗi này phân kỳ ở nơi khác, nhưng ζ(s)\zeta(s) mở rộng thành một hàm phân hình trên toàn bộ C\mathbb{C} với đúng một cực điểm đơn tại s=1s=1, thông qua phương trình hàm ζ(s)=2sπs−1sin⁡ ⁣(πs2)Γ(1−s) ζ(1−s)\zeta(s) = 2^s \pi^{s-1} \sin\!\left(\dfrac{\pi s}{2}\right) \Gamma(1-s)\, \zeta(1-s), biểu diễn ζ(s)\zeta(s) theo ζ(1−s)\zeta(1-s) và một thừa số Gamma. Đây là điều cho phép hiểu chặt chẽ những giá trị như ζ(−1)=−112\zeta(-1) = -\dfrac{1}{12}, vốn không liên quan gì tới việc cộng theo nghĩa đen 1+2+3+⋯1+2+3+\cdots.

ζ(s)=2sπs−1sin⁡ ⁣(πs2)Γ(1−s) ζ(1−s)\zeta(s) = 2^s \pi^{s-1} \sin\!\left(\dfrac{\pi s}{2}\right) \Gamma(1-s)\, \zeta(1-s)
Bốn hàm đặc biệt nhìn tổng quan
HàmCông thức định nghĩaTính chất chính
Gamma, Γ(z)\Gamma(z)Γ(z)=∫0∞tz−1e−t dt\Gamma(z)=\displaystyle\int_0^\infty t^{z-1}e^{-t}\,dtNội suy n!n!; cực điểm đơn tại 0,−1,−2,…0,-1,-2,\dots
Beta, B(x,y)B(x,y)B(x,y)=∫01tx−1(1−t)y−1 dtB(x,y) = \displaystyle\int_0^1 t^{x-1}(1-t)^{y-1}\,dtB(x,y)=Γ(x)Γ(y)Γ(x+y)B(x,y) = \dfrac{\Gamma(x)\Gamma(y)}{\Gamma(x+y)}
Zeta, ζ(s)\zeta(s)ζ(s)=∑n=1∞1ns\zeta(s) = \displaystyle\sum_{n=1}^{\infty} \dfrac{1}{n^s}Thác triển phân hình; một cực điểm tại s=1s=1
Bessel, Jν(x)J_\nu(x)Jν(x)=∑k=0∞(−1)kk! Γ(k+ν+1)(x2)2k+νJ_\nu(x) = \displaystyle\sum_{k=0}^{\infty} \dfrac{(-1)^k}{k!\,\Gamma(k+\nu+1)} \left(\dfrac{x}{2}\right)^{2k+\nu}Nghiệm của x2y′′+xy′+(x2−ν2)y=0x^2y''+xy'+(x^2-\nu^2)y=0

Nâng caoCác định lý then chốt

Với mọi số nguyên không âm nn, Γ(n+1)=n!\Gamma(n+1)=n!, trong đó Γ\Gamma được định nghĩa bởi Γ(z)=∫0∞tz−1e−t dt\Gamma(z)=\displaystyle\int_0^\infty t^{z-1}e^{-t}\,dt với Re⁡(z)>0\operatorname{Re}(z) > 0.

Vì sao đúng?

Điều này cho thấy tích phân Gamma không phải là một tổng quát hóa tùy tiện, mà là sự tiếp nối mượt mà duy nhất tái tạo đúng cấu trúc nhân của giai thừa, đó là lý do nó là đối tượng đúng để thay vào các số không nguyên.

Chứng minh

Bước 1 (Trường hợp cơ sở). Γ(1)=∫0∞t0e−t dt=∫0∞e−t dt=[−e−t]0∞=0−(−1)=1=0!\Gamma(1) = \int_0^\infty t^0 e^{-t}\,dt = \int_0^\infty e^{-t}\,dt = \big[-e^{-t}\big]_0^\infty = 0-(-1) = 1 = 0!, nên công thức đúng với n=0n=0.

Bước 2 (Công thức truy hồi bằng tích phân từng phần). Với Re⁡(z)>0\operatorname{Re}(z) > 0, lấy tích phân từng phần Γ(z)=∫0∞tz−1e−t dt\Gamma(z)=\displaystyle\int_0^\infty t^{z-1}e^{-t}\,dt với u=tzu=t^z, dv=e−t dtdv=e^{-t}\,dt: Γ(z+1)=∫0∞tze−t dt=[−tze−t]0∞+z∫0∞tz−1e−t dt=0+z Γ(z)\Gamma(z+1) = \int_0^\infty t^z e^{-t}\,dt = \big[-t^z e^{-t}\big]_0^\infty + z\int_0^\infty t^{z-1}e^{-t}\,dt = 0 + z\,\Gamma(z), dùng việc số hạng biên triệt tiêu ở cả hai đầu. Điều này chứng minh Γ(z+1)=z Γ(z)\Gamma(z+1) = z\,\Gamma(z).

Bước 3 (Quy nạp). Giả sử Γ(k+1)=k!\Gamma(k+1) = k! với một số nguyên không âm kk nào đó. Theo Bước 2 với z=k+1z=k+1, Γ(k+2)=(k+1) Γ(k+1)=(k+1)⋅k!=(k+1)!\Gamma(k+2) = (k+1)\,\Gamma(k+1) = (k+1)\cdot k! = (k+1)!.

Bước 4 (Kết luận). Vì trường hợp cơ sở n=0n=0 đúng (Bước 1) và bước quy nạp chuyển đẳng thức từ n=kn=k sang n=k+1n=k+1 (Bước 3), Γ(n+1)=n!\Gamma(n+1)=n! đúng với mọi số nguyên không âm nn theo nguyên lý quy nạp toán học.

Với mọi zz thỏa 0<Re⁡(z)<10 < \operatorname{Re}(z) < 1, Γ(z)Γ(1−z)=πsin⁡(πz)\Gamma(z)\Gamma(1-z) = \dfrac{\pi}{\sin(\pi z)}.

Vì sao đúng?

Nó liên kết một đối tượng trông thuần đại số (tỉ số hai giá trị Gamma) với lượng giác, và chính là đẳng thức khiến các cực điểm của Γ\Gamma tại các số nguyên không dương khớp với các không điểm của sin⁡(πz)\sin(\pi z) tại mọi số nguyên — phép tính thặng dư từ tích phân đường chính là thứ tạo ra hàm sin ở vế phải.

Chứng minh

Bước 1 (Quy về hàm Beta). Dùng B(x,y)=Γ(x)Γ(y)Γ(x+y)B(x,y) = \dfrac{\Gamma(x)\Gamma(y)}{\Gamma(x+y)} với y=1−zy=1-z, Γ(z)Γ(1−z)=B(z,1−z)=∫01tz−1(1−t)−z dt\Gamma(z)\Gamma(1-z) = B(z,1-z) = \int_0^1 t^{z-1}(1-t)^{-z}\,dt. Thế u=t/(1−t)u = t/(1-t) biến đổi thành ∫0∞uz−11+u du\displaystyle\int_0^\infty \dfrac{u^{z-1}}{1+u}\,du với 0<Re⁡(z)<10<\operatorname{Re}(z)<1.

Bước 2 (Thiết lập đường cong lỗ khóa). Để tính I=∫0∞uz−11+u duI=\int_0^\infty \dfrac{u^{z-1}}{1+u}\,du, dùng nhánh của wz−1w^{z-1} cắt dọc trục thực dương và lấy tích phân wz−11+w\dfrac{w^{z-1}}{1+w} quanh một đường cong lỗ khóa: đi ra ngay phía trên vết cắt, một đường tròn lớn bán kính RR, quay lại ngay phía dưới vết cắt, và một đường tròn nhỏ bán kính ε\varepsilon quanh 00.

Bước 3 (Thặng dư tại cực điểm duy nhất). Điểm kỳ dị duy nhất của wz−11+w\dfrac{w^{z-1}}{1+w} bên trong đường cong là cực điểm đơn tại w=−1w=-1, tại đó (dùng nhánh với argument π\pi ở đó) Res=(−1)z−1=eiπ(z−1)\mathrm{Res} = (-1)^{z-1} = e^{i\pi(z-1)}, nên định lý thặng dư cho tích phân lỗ khóa bằng 2πi eiπ(z−1)2\pi i\, e^{i\pi(z-1)}.

Bước 4 (So sánh hai phía vết cắt và cho R→∞R\to\infty, ε→0\varepsilon\to 0). Đường tròn lớn và nhỏ triệt tiêu khi 0<Re⁡(z)<10<\operatorname{Re}(z)<1, và hai đoạn thẳng khác nhau bởi pha mà wz−1w^{z-1} nhận được khi đi qua vết cắt, cho I(1−e2πi(z−1))=2πi eiπ(z−1)I\big(1 - e^{2\pi i(z-1)}\big) = 2\pi i\, e^{i\pi(z-1)}. Giải ra II và rút gọn các hàm mũ dùng sin⁡θ=eiθ−e−iθ2i\sin\theta = \frac{e^{i\theta}-e^{-i\theta}}{2i} cho I=πsin⁡(πz)I = \dfrac{\pi}{\sin(\pi z)}, nên kết hợp với Bước 1, Γ(z)Γ(1−z)=πsin⁡(πz)\Gamma(z)\Gamma(1-z) = \dfrac{\pi}{\sin(\pi z)}.

Nâng caoỨng dụng thực tiễn và Ví dụ minh họa

Hàm Gamma chuẩn hóa các phân phối xác suất dùng để mô hình hóa thời gian chờ và tiên nghiệm Bayes (Gamma, Beta, chi bình phương, Student's t); hàm Bessel là nghiệm tự nhiên của phương trình sóng và nhiệt trong hình học trụ, nên chúng xuất hiện trong âm học của mặt trống, mẫu nhiễu xạ quang học, và mẫu bức xạ ăng-ten; còn hàm zeta neo giữ lý thuyết số giải tích (sự phân bố số nguyên tố), đồng thời cũng xuất hiện trong vật lý qua phép chính quy hóa bằng hàm zeta, dùng để gán giá trị hữu hạn cho các tổng phân kỳ trong hiệu ứng Casimir và lý thuyết dây.

Ví dụ: Chuẩn hóa phân phối chi bình phương

Một nhà thống kê mô hình hóa tổng bình phương của k=5k=5 biến chuẩn tắc độc lập bằng mật độ chi bình phương f(x)=C xk/2−1e−x/2f(x) = C\, x^{k/2-1}e^{-x/2} với x>0x>0. Tìm hằng số chuẩn hóa CC sao cho ∫0∞f(x) dx=1\int_0^\infty f(x)\,dx = 1.

Lời giải

Bước 1 (Đổi biến để khớp tích phân Gamma). Đặt t=x/2t = x/2, nên x=2tx=2t, dx=2 dtdx = 2\,dt: ∫0∞xk/2−1e−x/2 dx=∫0∞(2t)k/2−1e−t 2 dt=2k/2∫0∞tk/2−1e−t dt\int_0^\infty x^{k/2-1}e^{-x/2}\,dx = \int_0^\infty (2t)^{k/2-1}e^{-t}\, 2\,dt = 2^{k/2}\int_0^\infty t^{k/2-1}e^{-t}\,dt.

Bước 2 (Nhận ra hàm Gamma). Tích phân còn lại chính xác là Γ(k/2)\Gamma(k/2) theo định nghĩa, nên ∫0∞xk/2−1e−x/2 dx=2k/2 Γ(k/2)\int_0^\infty x^{k/2-1}e^{-x/2}\,dx = 2^{k/2}\,\Gamma(k/2).

Bước 3 (Thay k=5k=5). Ở đây k/2=5/2k/2 = 5/2, và dùng công thức truy hồi Γ(5/2)=(3/2)(1/2)Γ(1/2)=34π\Gamma(5/2) = (3/2)(1/2)\Gamma(1/2) = \dfrac{3}{4}\sqrt{\pi}, nên tích phân bằng 42⋅34π=32π4\sqrt{2}\cdot\frac{3}{4}\sqrt{\pi}=3\sqrt{2\pi} sau khi rút gọn 25/2=422^{5/2}=4\sqrt{2}.

Bước 4 (Giải tìm CC). Vì C⋅32π=1C \cdot 3\sqrt{2\pi} = 1, hằng số chuẩn hóa là C=132πC = \dfrac{1}{3\sqrt{2\pi}}, khớp với công thức chi bình phương tổng quát C=12k/2Γ(k/2)C = \dfrac{1}{2^{k/2}\Gamma(k/2)}.

Ví dụ: Tần số cơ bản của mặt trống hình tròn

Một mặt trống tròn bán kính aa có một mode dao động đối xứng xuyên tâm u(r,t)=J0(kr)cos⁡(ωt)u(r,t) = J_0(kr)\cos(\omega t) phải triệt tiêu tại vành cứng, u(a,t)=0u(a,t)=0 với mọi tt. Biết không điểm dương đầu tiên của J0J_0 là j0,1≈2.405j_{0,1}\approx 2.405, hãy biểu diễn tần số góc cơ bản ω1\omega_1 theo aa và tốc độ sóng cc (với ω=ck\omega = ck).

Lời giải

Bước 1 (Áp dụng điều kiện biên). Vì u(a,t)=J0(ka)cos⁡(ωt)u(a,t)=J_0(ka)\cos(\omega t) phải triệt tiêu với mọi tt, ta cần J0(ka)=0J_0(ka) = 0, nên kaka phải là một không điểm của hàm Bessel J0J_0.

Bước 2 (Chọn mode cơ bản). Mode cơ bản ứng với kk dương nhỏ nhất, tức không điểm dương nhỏ nhất j0,1≈2.405j_{0,1}\approx 2.405 của J0J_0; do đó k1a=j0,1k_1 a = j_{0,1}, nên k1=j0,1/ak_1 = j_{0,1}/a.

Bước 3 (Chuyển sang tần số góc). Dùng hệ thức tán sắc ω=ck\omega = ck của phương trình sóng màng, ω1=ck1=c j0,1a\omega_1 = c k_1 = \dfrac{c\, j_{0,1}}{a}.

Bước 4 (Diễn giải). Thay giá trị số cho ω1≈2.405 ca\omega_1 \approx \dfrac{2.405\,c}{a}: tăng gấp đôi bán kính aa làm giảm một nửa tần số cơ bản, đúng như dự đoán vật lý rằng trống lớn hơn tạo ra âm trầm hơn.

Γ(6)\Gamma(6) bằng bao nhiêu?

Γ ⁣(12)\Gamma\!\left(\dfrac{1}{2}\right) bằng bao nhiêu?

Thác triển phân hình của ζ(s)\zeta(s) có cực điểm duy nhất ở đâu, và thuộc loại gì?

Một kỹ sư âm thanh thiết kế một mặt trống hình tròn và cần các hình dạng mode giữ hữu hạn tại tâm và triệt tiêu tại vành cứng. Họ cần dùng họ hàm đặc biệt nào để có các hình dạng mode xuyên tâm này?

Tài liệu tham khảo

  1. Frank W. J. Olver, Ronald F. Boisvert, Daniel W. Lozier, Charles W. Clark (eds.) (2023). NIST Digital Library of Mathematical Functions, Chapter 5: Gamma Function
  2. Frank W. J. Olver, Ronald F. Boisvert, Daniel W. Lozier, Charles W. Clark (eds.) (2023). NIST Digital Library of Mathematical Functions, Chapter 25: Zeta and Related Functions
  3. David J. Platt, Timothy S. Trudgian (2021). The Riemann Hypothesis Is True Up to 3×10^12 · arXiv:2004.09765