ما الذي تقوم به هذه الحاسبة
تنشئ هذه الأداة جدولاً لقيم دالة بيسل من النوع الأول، التي تُكتب \(J_{v}(x)\)، عند رتبة ثابتة \(v\) مع تغيير الوسيط \(x\) تدريجياً. تختار قيمة \(x\) الابتدائية، ومقدار الزيادة، وعدد الصفوف المطلوب توليدها، فتُعيد لك الحاسبة جدولاً أنيقاً من عمودين يقابل بين \(x\) والقيمة المناظرة لها \(J_{v}(x)\). وتظهر دوال بيسل من النوع الأول في شتى مجالات الفيزياء والهندسة: اهتزازات الطبلة الدائرية، وتوصيل الحرارة في الأسطوانات، والموجات الكهرومغناطيسية في أدلة الموجات، ومعالجة الإشارات (الفصوص الجانبية في تضمين التردد FM).
طريقة الاستخدام
أدخل الرتبة \(v\) (أي عدد حقيقي — 0 أو 1 أو 2، أو قيمة كسرية مثل 0.5، أو قيمة سالبة). ثم حدد القيمة الابتدائية لـ \(x\)، ومقدار الزيادة (المسافة بين كل قيمة وما يليها من قيم \(x\)؛ ويمكن أن تكون سالبة لمسح تنازلي أو صفراً لتكرار النقطة نفسها)، وعدد التكرارات (عدد الصفوف، من 1 وحتى 10000). يستخدم الصف رقم \(i\) القيمة \(x = \text{startX} + i \times \text{stepX}\).
شرح الصيغة
تُعرَّف الدالة بمتسلسلة القوى $$J_{v}(x) = \sum_{k=0}^{\infty} \frac{(-1)^{k}}{k!\,\Gamma(k+v+1)} \left(\frac{x}{2}\right)^{v+2k}$$ حيث \(\Gamma\) هي دالة جاما. تحسب الحاسبة هذه المتسلسلة حداً بحد باستخدام علاقة تكرارية مستقرة: يُشتق كل حد من الحد السابق بضربه في \(-\dfrac{x^{2}/4}{(k+1)(k+v+1)}\)، ما يتجنّب طفحان المضروب (factorial). وتُحسب دالة جاما عبر تقريب لانكزوس (Lanczos) بحيث تعمل الأداة مع الرتب غير الصحيحة والسالبة. أما للرتب الصحيحة السالبة فتُستخدم المتطابقة \(J_{-n}(x) = (-1)^{n} J_{n}(x)\).
مثال محلول
عند \(v = 0\) وstartX = 0 وstepX = 0.2 وloopCount = 6، يعطي الجدول \(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\)، و\(J_{0}(1.0) \approx 0.765198\) — وهو ما يطابق القيمة المجدولة القياسية \(J_{0}(1) = 0.7651976866\).
القيم المرجعية لـ J_v(x)
الجدول أدناه يُدرج دالة بيسل من النوع الأول \(J_v(x)\) للرتب \(v=0,1,2\) عند عدة متغيّرات معيارية. تُقرّب القيم إلى ستة منازل عشرية وتنتج من السلسلة \(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 |
كفحص يدوي، قيّم \(J_0(1)\): الحدود الرائدة تعطي \(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.
الأصفار البارزة (الجذور)
الأصفار الموجبة هي قيم \(x\) حيث \(J_v(x)=0\)؛ وهي تحدد أنماط الطبلة وحدود قطع الموجات وشروط حدية مشابهة.
| فهرس الجذر \(s\) | الصفر الـ \(s\) لـ \(J_0\) | الصفر الـ \(s\) لـ \(J_1\) |
|---|---|---|
| 1 | 2.404826 | 3.831706 |
| 2 | 5.520078 | 7.015587 |
| 3 | 8.653728 | 10.173468 |
| 4 | 11.791534 | 13.323692 |
لاحظ أن \(x=0\) هو صفر لـ \(J_v\) لكل رتبة \(v>0\)، لكنه لا يُعتبر من بين الجذور الموجبة المذكورة أعلاه.
التعاريف والمسرد
- الرتبة \(v\)
- المعامل (هنا حقل النموذج order) الذي يختار أي عضو من عائلة بيسل يتم حسابه. يمكن أن يكون أي رقم حقيقي — تنشأ الرتب الصحيحة في المسائل الأسطوانية، والرتب نصف الصحيحة \(v=n+\tfrac12\) تعطي دوال بيسل الكروية.
- المتغير \(x\)
- المتغير المستقل الذي يتم عنده تقييم \(J_v\). في هذا الجدول يبدأ من startX ويتقدم بمقدار stepX لـ loopCount صفوف.
- دالة جاما \(\Gamma\)
- التوسع المستمر للعاملي، حيث \(\Gamma(n+1)=n!\) للأعداد الصحيحة غير السالبة. تظهر في المقام \(\Gamma(v+k+1)\) في السلسلة بحيث يتم تعريف الرتب غير الصحيحة بشكل جيد.
- دالة بيسل من النوع الأول \(J_v(x)\)
- حل معادلة بيسل التفاضلية \(x^2 y''+x y'+(x^2-v^2)y=0\) الذي يبقى محدودًا عند الأصل (لـ \(v\ge 0\)). يُعطى بواسطة سلسلة القوى في الصيغة أعلاه.
- الأصفار / الجذور
- قيم \(x\) التي \(J_v(x)=0\). كل رتبة لها عدد لا نهائي من الأصفار الموجبة، بفواصل متزايدة موحدة وغير متساوية بشكل تقاربي بمقدار \(\pi\).
- الرتبة نصف الصحيحة (الكروية)
- عندما \(v=n+\tfrac12\)، \(J_v\) تتعلق بدوال بيسل الكروية \(j_n(x)=\sqrt{\tfrac{\pi}{2x}}\,J_{n+1/2}(x)\)، التي تصف الأجزاء الشعاعية من معادلات الموجة في الإحداثيات الكروية.
- نسبة حد التكرار
- الحدود المتتالية من السلسلة تحقق \(\frac{a_{k+1}}{a_k}=\frac{-(x/2)^2}{(k+1)(v+k+1)}\)، والتي تُستخدم داخليًا لتوليد كل حد من السابق وتقييم التقارب.
تفسير جدولك
تساعد بعض الحقائق في قراءة الأعمدة التي ينتجها الفحص:
- القيم الابتدائية. \(J_0(0)=1\)، بينما \(J_v(0)=0\) لكل رتبة \(v>0\). لذا يبدأ الجدول عند \(x=0\) بـ 1 فقط للرتبة الصفرية.
- التذبذب مع التضاؤل. بالنسبة للقيم الكبيرة من \(x\)، \(J_v(x)\approx\sqrt{\tfrac{2}{\pi x}}\cos\!\left(x-\tfrac{v\pi}{2}-\tfrac{\pi}{4}\right)\). الدالة تتذبذب مثل جيب التمام المزاح بينما تتناقص سعتها كـ \(1/\sqrt{x}\). لذا تتقلص القيم العظمى المتتالية ببطء مع نمو \(x\).
- تغييرات الإشارة تحدد الأصفار. حيثما يتغير عمود بين صفين، يوجد جذر لـ \(J_v\) في تلك الفترة (على سبيل المثال \(J_0\) يتغير بين \(x=2\) و \(x=3\)، مما يحصر جذره الأول \(\approx 2.4048\)). بالنسبة للمتغيرات الكبيرة، تتباعد الأصفار المتتالية بحوالي \(\pi\).
- العقد الفيزيائية. تلك الأصفار تقابل شروط الحدود الفيزيائية: الأنماط الشعاعية لرأس طبلة دائرية مهتزة، وترددات القطع في أنابيب الموجات الأسطوانية، وأنماط المجال في الألياف البصرية جميعها مفهرسة بواسطة أصفار \(J_v\).
- الحجم. بالنسبة لـ \(x\) ثابت، تبدأ الرتب الأعلى \(v\) بالقرب من الصفر وترتفع ببطء أكثر؛ بالنسبة للقيم الصغيرة من \(x\) السلوك الرائد هو \(J_v(x)\sim \frac{1}{\Gamma(v+1)}\left(\frac{x}{2}\right)^{v}\)، لذا تبقى القيم الأكبر \(v\) أصغر حتى يصبح \(x\) قابلًا للمقارنة مع \(v\).
تنتج هذه الملاحظات من السلسلة المعروفة والأشكال التقاربية أعلاه وتنطبق على أي رتبة تدخلها.
الأسئلة الشائعة
هل يمكن أن تكون الرتبة كسرية أو سالبة؟ نعم. تدعم المتسلسلة المعتمدة على دالة جاما أي رتبة حقيقية، بما في ذلك أنصاف الأعداد الصحيحة (التي تعطي صيغ بيسل الكروية) والقيم السالبة.
ماذا يحدث عند \(x = 0\)؟ تكون \(J_{0}(0) = 1\) بينما \(J_{v}(0) = 0\) لكل \(v > 0\)، لأن العامل الرائد \((x/2)^{v}\) يتلاشى.
ما مدى دقتها عند قيم \(x\) الكبيرة؟ المتسلسلة بدقة مزدوجة (double precision) دقيقة في النطاقات المعتادة (\(x\) حتى نحو 20–30). أما عند قيم \(x\) الكبيرة جداً فقد يؤدي الإلغاء الكارثي إلى تراجع الدقة؛ ويُفضَّل في هذه الحالة استخدام الصيغة التقاربية $$J_{v}(x) \approx \sqrt{\frac{2}{\pi x}} \cos\left(x - \frac{v\pi}{2} - \frac{\pi}{4}\right).$$
