ماذا تفعل هذه الحاسبة
تتيح لك هذه الأداة إيجاد الصفر الموجب رقم s لدوال بيسل من النوع الأول \(J_{v}(x)\) ومن النوع الثاني \(Y_{v}(x)\) (المعروفة أيضًا بدالة نيومان أو فيبر) لأي رتبة حقيقية \(v\). وتُرمَز هذه الأصفار بالرمزين \(j_{v,s}\) و\(y_{v,s}\)، وهي تظهر في كثير من تطبيقات الفيزياء والهندسة، مثل اهتزاز الأغشية الدائرية (أوجه الطبول)، وانتقال الحرارة في الأسطوانات، والأدلة الموجية الكهرومغناطيسية، ومتسلسلات فورييه-بيسل. وهي أداة رياضية بحتة للدوال الخاصة، صالحة عالميًا دون أي ارتباط بمنطقة أو وحدة قياس معينة.
طريقة الاستخدام
أدخل قيمة الرتبة \(v\) (عدد حقيقي، يتراوح عادةً بين 0 و200) ثم دليل الصفر \(s\) (عدد صحيح موجب: 1، 2، 3، …). تُعيد لك الحاسبة القيمة \(j_{v,s}\)، وهي الجذر الموجب رقم s للدالة \(J_{v}(x)\)، إضافةً إلى \(y_{v,s}\)، وهو الجذر الموجب رقم s للدالة \(Y_{v}(x)\). على سبيل المثال، تُعطي القيمتان \(v = 0\) و\(s = 1\) النمط الأساسي لاهتزاز وجه الطبل.
الصيغة والطريقة الحسابية
تحقّق كلتا الدالتين معادلة بيسل: $$x^2 y'' + x y' + (x^2 - v^2) y = 0.$$ تُحسَب \(J_{v}\) من متسلسلتها الأُسّية المتقاربة، بينما تُحسَب \(Y_{v}\) عبر العلاقة $$Y_{v} = \frac{J_{v}\cos(v\pi) - J_{-v}}{\sin(v\pi)},$$ مع معالجة حالة الرتب الصحيحة كنهاية عددية. ولتحديد موقع الصفر رقم s، نبدأ بتقدير مستمد من صيغة ماكماهون التقاربية، ثم نحصر الجذر داخل مجال، ونصقله بطريقة التنصيف حتى التقارب. وبما أنّ \(Y_{v}(x)\) تؤول إلى \(-\infty\) عند \(x = 0\)، يبدأ البحث من قيمة موجبة صغيرة لتفادي النقطة الشاذة اللوغاريتمية.
مثال محلول
عند \(v = 0\) و\(s = 1\): يكون أول صفر موجب للدالة \(J_{0}(x)\) مساويًا \(2.4048255577\)، وأول صفر موجب للدالة \(Y_{0}(x)\) مساويًا \(0.8935769663\). وعند \(v = 1\) و\(s = 1\): نحصل على \(j_{1,1} = 3.8317059702\) و\(y_{1,1} = 2.1971413260\).
ماذا تعني الأصفار في التطبيقات
أصفار دالة بيسل ليست فضولاً مجردًا — بل هي القيم الذاتية المنفصلة التي تظهر كلما تم طرح مسألة موجة أو حرارة أو جهد على مجال دائري أو أسطواني. شرط الحد يفرض على الجزء الشعاعي من الحل أن يختفي (أو يختفي مشتقه) عند الحد، وهذا الشرط مرضي فقط عند الأصفار \(j_{v,s}\).
غشاء دائري متذبذب (رأس الطبل)
بالنسبة لطبل دائري مثالي بنصف قطر \(a\) مع حافة مثبتة، ينقسم الإزاحة إلى أنماط معنونة بـ \((v,s)\)، حيث يعد \(v\) أقطار العقدة الزاويّة و\(s\) يعد الدوائر العقدية الشعاعية. الترددات الذاتية المسموحة هي \(f_{v,s}=\frac{c}{2\pi a}\,j_{v,s}\)، حيث \(c\) هي سرعة الموجة. النغمة الأساسية تستخدم \(j_{0,1}=2.4048256\)؛ وكل من \(s\) و\(v\) الأعلى يرفع النبرة، وبسبب أن \(j_{v,s}\) ليست مضاعفات صحيحة لبعضها البعض فإن النغمات الزائدة للطبل غير متوافقة.
التوصيل الحراري الأسطواني
عند حل معادلة الحرارة في أسطوانة طويلة بسطح درجة حرارة ثابتة، تكون الدوال الذاتية الشعاعية \(J_0(\lambda_s r/a)\) مع \(\lambda_s=j_{0,s}\). يتحلل كل نمط مع الزمن كـ \(\exp\!\left(-\alpha (j_{0,s}/a)^2 t\right)\)، وبالتالي فإن أصغر صفر \(j_{0,1}\) يحكم الملف الشخصي لدرجة الحرارة الذي يتحلل ببطء والأطول عمراً. يعطي \(s\) الأكبر قيماً ذاتية أكبر وبالتالي يتحلل بشكل أسرع.
ترددات قطع الموجهات
في موجه معادن دائري مجوف بنصف قطر \(a\)، تقطع الأنماط المغناطيسية العرضية (TM) عند ترددات يحددها \(j_{v,s}\) والأنماط الكهربائية العرضية (TE) بأصفار المشتقة \(J_v'\). بالنسبة لأنماط TM فإن قطع هو \(f_{c}=\frac{c\,j_{v,s}}{2\pi a}\)؛ فقط فوق هذا التردد ينتشر النمط. يستخدم أقل نمط TM (TM\(_{01}\)) مرة أخرى \(j_{0,1}\).
سلسلة فورييه–بيسل
يمكن توسيع أي دالة معقولة على قرص كـ \(f(r)=\sum_{s=1}^{\infty} c_s\,J_v\!\left(j_{v,s}\,r/a\right)\). تعمل الأصفار المقاسة \(j_{v,s}/a\) تماماً مثل الأرقام الموجية \(n\pi/L\) لسلسلة جيب فورييه عادية، والتعامد للـ \(J_v(j_{v,s}r/a)\) فوق القرص (مع الوزن \(r\)) يسمح بحساب المعاملات بالتكامل، \(c_s=\frac{2}{a^2 J_{v+1}^2(j_{v,s})}\int_0^a f(r)\,J_v\!\left(j_{v,s}r/a\right) r\,dr\).
كيف تتحول الأصفار
بالنسبة لترتيب ثابت، فإن زيادة الفهرس \(s\) تزيد \(j_{v,s}\) في خطوات تقترب من \(\pi\) (التذبذب يصبح دوريًا تقريبًا بعيداً عن الأصل). بالنسبة لفهرس ثابت، فإن زيادة الترتيب \(v\) تدفع الصفر الأول للخارج تقريباً مثل \(j_{v,1}\approx v + 1.8557\,v^{1/3}+\cdots\) للـ \(v\) الكبيرة، مما يعكس كيف أن الحد المركزي \(v^2/x^2\) في معادلة بيسل يؤخر بدء التذبذب. الفائدة العملية: التعقيد الزاوي الأعلى (\(v\)) والعقد الشعاعية الأكثر (\(s\)) كلاهما يتوافق مع قيم ذاتية أكبر وبالتالي ترددات أعلى أو تحلل أسرع.
للحصول على مثال طبل محدد مع \(c=100\ \text{م/ث}\) و\(a=0.30\ \text{م}\)، التردد الأساسي هو \(f_{0,1}=\frac{100}{2\pi(0.30)}\,(2.4048256)\approx 127.6\ \text{Hz}\). النغمة الزائدة الأولى تستخدم \(j_{1,1}=3.8317060\)، مما يعطي \(\approx 203\ \text{Hz}\)، وبالتالي فإن نسبة النغمة الزائدة إلى الأساسية هي \(j_{1,1}/j_{0,1}\approx 1.593\) — يمكن سماع الفرق وليس أوكتافًا أو خمسة، وهذا هو السبب في أن الطبول تبدو غير منغمة مقارنة بالأوتار.
الأسئلة الشائعة
لماذا يقع أول صفر للدالة \(Y_{v}\) دون الواحد؟ لأنّ \(Y_{0}\) تتباعد نحو \(-\infty\) عند نقطة الأصل، ثم تعبر الصفر عند \(x \approx 0.894\) تقريبًا قبل أن تبلغ قيمتها العظمى الأولى، ولذلك يكون أول صفر لها أصغر بكثير من أول صفر للدالة \(J_{0}\).
هل يمكن أن تكون \(v\) غير صحيحة؟ نعم. فالصيغة المعرِّفة للدالة \(Y_{v}\) صالحة مباشرةً لأي رتبة \(v\) غير صحيحة، أما الرتب الصحيحة فتُعالَج كنهاية سلسة.
ما مدى دقة النتائج؟ تستخدم الحسابات الحساب مزدوج الدقة (double precision)، مما يمنح نحو 10 أرقام معنوية لقيم \(v\) و\(s\) المعتدلة.
