Công cụ này làm gì
Công cụ này lập bảng hàm Bessel loại một, ký hiệu \(J_{v}(x)\), với bậc v cố định trong khi quét giá trị của biến x. Bạn chọn giá trị x ban đầu, bước nhảy và số dòng cần tạo, sau đó công cụ trả về một bảng hai cột gọn gàng thể hiện x và \(J_{v}(x)\) tương ứng. Hàm Bessel loại một xuất hiện ở khắp nơi trong vật lý và kỹ thuật: dao động của màng trống tròn, dẫn nhiệt trong hình trụ, sóng điện từ trong ống dẫn sóng và xử lý tín hiệu (các dải biên trong điều chế tần số FM).
Cách sử dụng
Nhập Bậc v (số thực bất kỳ — 0, 1, 2, phân số như 0,5 hay số âm). Đặt Giá trị x ban đầu, Bước nhảy (khoảng cách giữa các giá trị x liên tiếp; có thể âm để quét giảm dần hoặc bằng 0 để lặp lại một điểm), và Số lần lặp (số dòng, từ 1 đến 10000). Dòng thứ i dùng \(x = \text{startX} + i \times \text{stepX}\).
Giải thích công thức
Hàm được định nghĩa bằng chuỗi lũy thừa $$J_{v}(x) = \sum_{k=0}^{\infty} \frac{(-1)^{k}}{k!\,\Gamma(k+v+1)} \left(\frac{x}{2}\right)^{v+2k}$$ trong đó \(\Gamma\) là hàm gamma. Công cụ tính chuỗi này theo từng số hạng bằng một công thức truy hồi ổn định: mỗi số hạng được suy ra từ số hạng trước bằng cách nhân với \(-\dfrac{x^{2}/4}{(k+1)(k+v+1)}\), nhờ đó tránh được hiện tượng tràn số do giai thừa. Hàm gamma được tính bằng xấp xỉ Lanczos để hỗ trợ cả bậc không nguyên và bậc âm. Với bậc nguyên âm, công cụ sử dụng đẳng thức \(J_{-n}(x) = (-1)^{n} J_{n}(x)\).
Ví dụ minh họa
Với \(v = 0\), \(\text{startX} = 0\), \(\text{stepX} = 0{,}2\), \(\text{loopCount} = 6\), bảng cho ra \(J_{0}(0) = 1\), \(J_{0}(0{,}2) \approx 0{,}990025\), \(J_{0}(0{,}4) \approx 0{,}960398\), \(J_{0}(0{,}6) \approx 0{,}912005\), \(J_{0}(0{,}8) \approx 0{,}846287\), và \(J_{0}(1{,}0) \approx 0{,}765198\) — khớp với giá trị tra bảng chuẩn \(J_{0}(1) = 0{,}7651976866\).
Giá Trị Tham Chiếu của J_v(x)
Bảng dưới đây liệt kê hàm Bessel loại thứ nhất \(J_v(x)\) cho các bậc \(v=0,1,2\) tại nhiều đối số tiêu chuẩn. Các giá trị được làm tròn đến sáu chữ số thập phân và tuân theo chuỗi \(J_{v}(x)=\sum_{k=0}^{\infty}\frac{(-1)^{k}}{k!\,\Gamma(v+k+1)}\left(\frac{x}{2}\right)^{2k+v}\).
| \(x\) | \(J_0(x)\) | \(J_1(x)\) | \(J_2(x)\) |
|---|---|---|---|
| 0 | 1.000000 | 0.000000 | 0.000000 |
| 0.5 | 0.938470 | 0.242268 | 0.030604 |
| 1 | 0.765198 | 0.440051 | 0.114903 |
| 2 | 0.223891 | 0.576725 | 0.352834 |
| 3 | −0.260052 | 0.339059 | 0.486091 |
| 5 | −0.177597 | −0.327579 | 0.046565 |
| 10 | −0.245936 | 0.043473 | 0.254630 |
Để kiểm tra bằng một ví dụ cụ thể, tính \(J_0(1)\): các số hạng đầu tiên cho \(1-\tfrac{(0.5)^2}{1}+\tfrac{(0.5)^4}{4}-\tfrac{(0.5)^6}{36}+\dots = 1-0.25+0.015625-0.000434+\dots\approx\) 0.765198.
Các Không Điểm Đáng Chú Ý (Căn Số)
Các không điểm dương là các giá trị của \(x\) khi \(J_v(x)=0\); chúng xác định các chế độ dao động của trống, tần số cắt của dẫn sóng và các điều kiện biên tương tự.
| Chỉ số căn \(s\) | Không điểm thứ \(s\) của \(J_0\) | Không điểm thứ \(s\) của \(J_1\) |
|---|---|---|
| 1 | 2.404826 | 3.831706 |
| 2 | 5.520078 | 7.015587 |
| 3 | 8.653728 | 10.173468 |
| 4 | 11.791534 | 13.323692 |
Lưu ý rằng \(x=0\) là một không điểm của \(J_v\) cho mọi bậc \(v>0\), nhưng nó không được tính trong các căn dương ở trên.
Định Nghĩa & Bảng Thuật Ngữ
- Bậc \(v\)
- Tham số (ở đây trường mẫu bậc) xác định thành viên nào của họ Bessel được tính toán. Nó có thể là bất kỳ số thực nào — các bậc nguyên phát sinh trong các bài toán hình trụ, các bậc nửa nguyên \(v=n+\tfrac12\) cho các hàm Bessel hình cầu.
- Đối số \(x\)
- Biến độc lập tại đó \(J_v\) được tính. Trong bảng này nó bắt đầu từ startX và tăng theo stepX cho loopCount hàng.
- Hàm Gamma \(\Gamma\)
- Sự mở rộng liên tục của giai thừa, với \(\Gamma(n+1)=n!\) cho các số nguyên không âm. Nó xuất hiện trong mẫu \(\Gamma(v+k+1)\) của chuỗi để các bậc không nguyên được xác định tốt.
- Hàm Bessel loại thứ nhất \(J_v(x)\)
- Nghiệm của phương trình vi phân Bessel \(x^2 y''+x y'+(x^2-v^2)y=0\) tại điểm gốc (với \(v\ge 0\)). Nó được cho bởi chuỗi luỹ thừa trong công thức ở trên.
- Không điểm / căn số
- Các giá trị của \(x\) khi \(J_v(x)=0\). Mỗi bậc có vô số không điểm dương, ngày càng cách đều nhau và tiệm cận được tách bởi \(\pi\).
- Bậc nửa nguyên (hình cầu)
- Khi \(v=n+\tfrac12\), \(J_v\) liên hệ với các hàm Bessel hình cầu \(j_n(x)=\sqrt{\tfrac{\pi}{2x}}\,J_{n+1/2}(x)\), mô tả các phần bán kính của các phương trình sóng trong tọa độ hình cầu.
- Tỉ lệ số hạng truy hồi
- Các số hạng liên tiếp của chuỗi thỏa mãn \(\frac{a_{k+1}}{a_k}=\frac{-(x/2)^2}{(k+1)(v+k+1)}\), được sử dụng nội bộ để sinh ra mỗi số hạng từ số hạng trước đó và để đánh giá sự hội tụ.
Diễn Giải Bảng Của Bạn
Một số dữ kiện giúp đọc các cột mà quét của bạn tạo ra:
- Các giá trị khởi đầu. \(J_0(0)=1\), trong khi \(J_v(0)=0\) cho mọi bậc \(v>0\). Vì vậy một bảng bắt đầu từ \(x=0\) chỉ bắt đầu ở 1 đối với bậc không.
- Dao động với giảm dần. Đối với \(x\) lớn, \(J_v(x)\approx\sqrt{\tfrac{2}{\pi x}}\cos\!\left(x-\tfrac{v\pi}{2}-\tfrac{\pi}{4}\right)\). Hàm dao động như một cosin dịch pha trong khi biên độ của nó giảm như \(1/\sqrt{x}\). Do đó các cực đại liên tiếp co lại chậm khi \(x\) tăng.
- Đổi dấu đánh dấu các không điểm. Bất cứ nơi nào một cột thay đổi dấu giữa hai hàng, một căn của \(J_v\) nằm trong khoảng đó (ví dụ: \(J_0\) chuyển đổi dấu giữa \(x=2\) và \(x=3\), khoanh vùng căn thứ nhất của nó \(\approx 2.4048\)). Đối với các đối số lớn, các không điểm liên tiếp được cách nhau khoảng \(\pi\).
- Nút vật lý. Các không điểm đó tương ứng với các điều kiện biên vật lý: các chế độ bán kính của một trống tròn dao động, tần số cắt của các dẫn sóng hình trụ, và các mẫu trường trong các sợi quang học đều được lập chỉ mục bởi các không điểm của \(J_v\).
- Độ lớn. Đối với \(x\) cố định, các bậc cao hơn \(v\) bắt đầu gần không và tăng chậm hơn; đối với \(x\) nhỏ, hành vi hàng đầu là \(J_v(x)\sim \frac{1}{\Gamma(v+1)}\left(\frac{x}{2}\right)^{v}\), vì vậy \(v\) lớn hơn vẫn nhỏ hơn cho đến khi \(x\) trở nên so sánh được với \(v\).
Những quan sát này tuân theo các chuỗi và các dạng tiệm cận đã được thiết lập ở trên và áp dụng cho bất kỳ bậc nào bạn nhập.
Câu hỏi thường gặp
Bậc có thể là phân số hay số âm không? Có. Chuỗi dựa trên hàm gamma hỗ trợ mọi bậc thực, kể cả bậc bán nguyên (cho ra dạng Bessel cầu) và bậc âm.
Tại x = 0 thì sao? \(J_{0}(0) = 1\) và \(J_{v}(0) = 0\) với \(v > 0\), vì thừa số dẫn đầu \((x/2)^{v}\) triệt tiêu.
Độ chính xác với x lớn ra sao? Chuỗi tính bằng độ chính xác kép cho kết quả tốt trong các dải thông thường (x đến khoảng 20–30). Với x rất lớn, hiện tượng triệt tiêu nghiêm trọng có thể làm giảm độ chính xác; trong trường hợp đó nên dùng dạng tiệm cận $$J_{v}(x) \approx \sqrt{\frac{2}{\pi x}} \cos\left(x - \frac{v\pi}{2} - \frac{\pi}{4}\right)$$
