Công cụ lập bảng hàm Bessel Y là gì?
Công cụ này lập bảng hàm Bessel loại hai — còn gọi là hàm Weber hay hàm Neumann, ký hiệu \(Y_{v}(x)\). Đây là nghiệm thứ hai độc lập tuyến tính của phương trình vi phân Bessel. Với một bậc thực \(v\) cố định, công cụ tính giá trị \(Y_{v}(x)\) tại một dãy các điểm \(x\) được xác định bởi giá trị khởi đầu, bước nhảy và số điểm, từ đó tạo ra một bảng số đầy đủ.
Cách sử dụng
Bạn nhập bậc \(v\) (có thể là số không nguyên hoặc số âm), giá trị \(x\) ban đầu, bước nhảy (khoảng cách) giữa các điểm, và số lần lặp (số dòng). Công cụ sẽ dựng dãy \(x_i = \text{startX} + i \cdot \text{stepX}\) với \(i\) chạy từ 0 đến \(\text{pointCount} - 1\), rồi liệt kê \(Y_{v}(x)\) cho từng giá trị. Lưu ý rằng \(Y_{v}(x)\) phân kỳ tới âm vô cực tại \(x = 0\) và chỉ nhận giá trị thực khi \(x > 0\), nên mọi dòng có \(x \le 0\) đều được đánh dấu là không xác định.
Công thức
Với bậc không nguyên:
$$Y_{\nu}(x) = \frac{J_{\nu}(x)\cos(\nu\pi) - J_{-\nu}(x)}{\sin(\nu\pi)}$$Với bậc nguyên \(n\), lấy giới hạn ta được công thức đóng gồm một số hạng logarit chứa \(J_{n}(x)\cdot\ln(x/2)\), một phần hiệu chỉnh chuỗi lũy thừa hữu hạn, và một chuỗi digamma. Hàm Bessel loại một \(J_{v}(x)\) được tính bằng cách lấy tổng chuỗi lũy thừa của nó, còn hàm Gamma được tính bằng phép xấp xỉ Lanczos.
Ví dụ minh họa
Với \(v = 0\), \(\text{startX} = 0\), \(\text{stepX} = 0{,}2\), \(\text{pointCount} = 51\), các dòng sẽ chạy từ \(x = 0{,}0\) đến \(10{,}0\). \(Y_{0}(0)\) không xác định (\(-\infty\)), \(Y_{0}(0{,}2) \approx -1{,}0811\), \(Y_{0}(1{,}0) \approx 0{,}0883\), \(Y_{0}(2{,}0) \approx 0{,}5104\), và \(Y_{0}(10{,}0) \approx 0{,}0557\). Dòng kết quả nổi bật "giá trị hữu hạn đầu tiên" sẽ hiển thị \(-1{,}0811\).
Định nghĩa & Bảng Thuật ngữ
- Bậc \(\nu\)
-
Tham số (trường
order) lập chỉ số cho họ các hàm Bessel. Nó có thể là bất kỳ số thực nào. Các bậc nguyên (0, 1, 2, …) là phổ biến nhất trong các bài toán vật lý có đối xứng trụ; các bậc nửa nguyên cho hàm Bessel hình cầu. - Hàm Bessel loại thứ hai \(Y_\nu(x)\)
- Còn gọi là hàm Weber hoặc Neumann (đôi khi viết là \(N_\nu\)). Nó là một nghiệm của phương trình Bessel không bị chặn (kỳ dị) tại gốc. Được định nghĩa cho \(\nu\) không nguyên bởi \(Y_\nu(x) = \dfrac{J_\nu(x)\cos(\nu\pi) - J_{-\nu}(x)}{\sin(\nu\pi)}\), với trường hợp nguyên được lấy làm giới hạn.
- \(J_\nu\) so với \(Y_\nu\)
- \(J_\nu(x)\) (loại thứ nhất) hữu hạn tại \(x=0\); \(Y_\nu(x)\) (loại thứ hai) phân kỳ đến \(-\infty\) khi \(x\to 0^+\). Cùng nhau chúng tạo thành cặp đầy đủ các nghiệm độc lập của phương trình Bessel.
- Phương trình vi phân Bessel
- Phương trình ODE tuyến tính \(x^2 y'' + x y' + (x^2 - \nu^2) y = 0\). Nghiệm tổng quát của nó là \(y = c_1 J_\nu(x) + c_2 Y_\nu(x)\).
- Hàm Gamma \(\Gamma(z)\)
- Sự mở rộng liên tục của giai thừa, \(\Gamma(n+1) = n!\), xuất hiện trong các hệ số chuỗi của \(J_\nu\) và \(Y_\nu\).
- Hàm Digamma \(\psi(z)\)
- Đạo hàm logarit \(\psi(z) = \Gamma'(z)/\Gamma(z)\). Nó xuất hiện một cách rõ ràng trong chuỗi cho \(Y_n(x)\) bậc nguyên, chứa một số hạng logarit \(\tfrac{2}{\pi}\ln(x/2)J_n(x)\) cộng với các hệ số có trọng số digamma.
- Xấp xỉ Lanczos
- Một phương pháp số rất chính xác để đánh giá hàm Gamma \(\Gamma(z)\) cho đối số phức hoặc thực, thường được sử dụng bên trong các quy trình hàm Bessel để tính các hệ số chuỗi.
- Nghiệm độc lập tuyến tính
- Một nghiệm thứ hai không thể được biểu thị dưới dạng bội số hằng của nghiệm thứ nhất. Bởi vì \(J_\nu\) một mình không thể biểu diễn các nghiệm kỳ dị tại gốc, \(Y_\nu\) cung cấp người bạn độc lập cần thiết cho nghiệm tổng quát.
Câu hỏi thường gặp
Vì sao dòng đầu tiên lại không xác định? \(Y_{v}(x)\) có một điểm kỳ dị tại \(x = 0\), phân kỳ tới \(-\infty\), nên tại đó nó không có giá trị hữu hạn.
Bậc \(v\) có thể là số âm không? Có. Với bậc nguyên âm, ta áp dụng tính đối xứng \(Y_{-n}(x) = (-1)^{n}Y_{n}(x)\); còn với bậc không nguyên âm thì dùng trực tiếp công thức tổng quát.
Độ chính xác ra sao? Các chuỗi được cộng dồn cho đến khi số hạng nhỏ hơn sai số máy, cho khoảng 6–7 chữ số có nghĩa với các giá trị \(x\) vừa phải.
