Công cụ này làm gì
Công cụ này lập bảng giá trị hàm Bessel biến đổi loại một, ký hiệu \(I_{v}(x)\), với một bậc thực \(v\) cố định trên một dãy giá trị \(x\). Bạn chỉ cần nhập bậc, giá trị x ban đầu, bước nhảy (số gia) và số dòng cần tạo; công cụ sẽ dựng dãy \(x_{i} = \text{đầu} + i\cdot\text{bước}\) rồi tính \(I_{v}(x_{i})\) tại từng điểm, trả về cả bảng số liệu lẫn đồ thị. Đây là công cụ thuần toán học về hàm đặc biệt, áp dụng được ở mọi nơi (không phụ thuộc quy định vùng miền hay đơn vị nào).
Công thức
Hàm Bessel biến đổi \(I_{v}(x)\) là nghiệm của phương trình Bessel biến đổi \(x^{2}y'' + xy' - (x^{2} + v^{2})y = 0\). Ở đây nó được tính từ chuỗi lũy thừa:
$$I_{v}(x) = \sum_{k=0}^{\infty} \frac{1}{k!\;\Gamma(v+k+1)}\left(\frac{x}{2}\right)^{v+2k}$$Nhờ giai thừa và hàm Gamma, \(v\) có thể là số thực bất kỳ. Để đảm bảo ổn định số học, mỗi số hạng được tính trong miền logarit bằng xấp xỉ Lanczos của \(\ln\Gamma\), sau đó cộng dồn cho đến khi các số hạng trở nên không đáng kể.
Cách sử dụng
Nhập Bậc v (ví dụ 0, 1 hoặc 2.5), Giá trị x ban đầu, Số gia cộng vào \(x\) ở mỗi dòng và Số lần lặp (số dòng). Bấm tính để nhận bảng hai cột gồm \(x\) và \(I_{v}(x)\) cùng đồ thị trên cùng khoảng giá trị đó.
Ví dụ minh họa
Với \(v = 0\), đầu = 0, bước = 0.5, số dòng = 5 bạn được \(x = 0, 0.5, 1, 1.5, 2\) và:
\(I_{0}(0) = 1\), \(I_{0}(0.5) \approx 1.0634834\), \(I_{0}(1) \approx 1.2660658\), \(I_{0}(1.5) \approx 1.6467232\), \(I_{0}(2) \approx 2.2795853\). Các kết quả này khớp với bảng tra chuẩn.
Câu hỏi thường gặp
Bậc có thể âm hoặc không nguyên không? Có. Với bậc nguyên âm, công cụ dùng đẳng thức \(I_{-n}(x) = I_{n}(x)\). Bậc \(v\) không nguyên được hỗ trợ khi \(x \geq 0\); còn với \(x < 0\) mà \(v\) không nguyên thì giá trị là số phức, nên kết quả trả về là NaN.
Vì sao \(I_{v}(x)\) tăng nhanh đến vậy? Khác với hàm Bessel thường \(J_{v}\) dao động, hàm biến đổi tăng xấp xỉ theo \(e^{x}/\sqrt{2\pi x}\) khi \(x\) lớn, nên với \(x\) lớn giá trị có thể tràn số thành vô cực.
\(I_{v}(0)\) bằng bao nhiêu? \(I_{0}(0) = 1\), và \(I_{v}(0) = 0\) với \(v > 0\).