حاسبة قاعدة سيمبسون

قاعدة سيمبسون هي طريقة للتكامل العددي تُقارب قيمة التكامل المحدد باستخدام كثيرات الحدود التربيعية. فبدلاً من تقدير المساحة تحت المنحنى باستخدام مستطيلات أو أشباه منحرفات، تُركّب قطعاً مكافئاً يمر عبر كل مجموعة من ثلاث نقاط عينة متتالية، وتحسب المساحة الدقيقة تحت كل قطع مكافئ. ويمنح مجموع هذه المساحات القطعية المكافئة التقريب المطلوب.

سُميت هذه الطريقة نسبةً إلى توماس سيمبسون، وهو رياضي إنجليزي من القرن الثامن عشر، مع أن الفكرة الأساسية ظهرت قبل ذلك في أعمال آخرين. وتكمن قوتها في ملاحظة بسيطة: يمكن للقطع المكافئ أن يتبع منحنى أملس عن قرب أكبر بكثير مما يستطيعه الخط المستقيم. فعندما تكون الدالة المُكاملة أملسة، تلتصق القطع المكافئة بالمنحنى الحقيقي بإحكام، ويكون التقريب الناتج دقيقاً بشكل ملحوظ حتى مع عدد متواضع من الفترات الفرعية.

ولفهم سبب أهمية النقاط الثلاث، تذكّر أن الخط المستقيم يتحدد تماماً بنقطتين، بينما يحتاج القطع المكافئ (كثيرة حدود تربيعية) إلى ثلاث نقاط. تستغل قاعدة سيمبسون ذلك بتقسيم فترة التكامل إلى عدد زوجي من الفترات الفرعية، وتجميعها في أزواج، وتركيب قطع مكافئ واحد عبر كل زوج. يمتد كل زوج على فترتين فرعيتين وثلاث نقاط: نقطة الطرف الأيسر، ونقطة المنتصف، ونقطة الطرف الأيمن.

تُسمى الصيغة الكلاسيكية قاعدة سيمبسون 1/3 لأن العامل 1/3 يظهر في الصيغة. وهناك أيضاً قاعدة سيمبسون 3/8 التي تستخدم كثيرات حدود تكعيبية عبر مجموعات من ثلاث فترات فرعية، لكن قاعدة 1/3 هي الأداة المعيارية والنسخة المُطبقة في معظم الحاسبات، بما فيها الحاسبة الموجودة في هذه الصفحة.

تنتمي قاعدة سيمبسون إلى عائلة صيغ نيوتن-كوتس، التي تُقارب التكاملات بتقييم الدالة المُكاملة عند نقاط متساوية التباعد ودمج تلك القيم بأوزان ثابتة. وما يميز قاعدة سيمبسون هو نمط أوزانها: 1، 4، 2، 4، 2، ...، 4، 1. تحصل نقطتا الطرفين على الوزن 1، والنقاط الداخلية ذات المواضع الفردية على الوزن 4، والنقاط الداخلية ذات المواضع الزوجية على الوزن 2. وهذا النمط المتناوب هو بصمة الطريقة ومصدر دقتها.

طول قطعة الخط المستقيم سهل: قِس المسافة بين نقطتي طرفيها. أما المنحنى فمختلف؛ إذ ينحني ويلتف، فلا يمكن لأي قياس مستقيم واحد أن يلتقط امتداده الحقيقي. ولإيجاد طول منحنى ما، يقسّمه الرياضيون إلى قطع صغيرة جداً، ويُقاربون كل قطعة بأنها مستقيمة، ثم يجمعون النتائج.

تخيّل أنك تمشي على درب متعرج. إذا أخذت خطوات صغيرة جداً، تكون كل خطوة مستقيمة تقريباً، مع أن الدرب ككل منحنٍ. اجمع كل تلك الخطوات المستقيمة الصغيرة وستحصل على تقدير جيد للمسافة الكلية المقطوعة. اجعل الخطوات أصغر فأصغر، وسيتقارب التقدير نحو الطول الدقيق. وهذه هي الفكرة الجوهرية وراء طول القوس.

بالنسبة للدالة y = f(x) على الفترة [a, b]، انظر إلى خطوة أفقية صغيرة جداً dx. خلال تلك الخطوة، ترتفع الدالة بمقدار dy = f'(x) dx تقريباً، حيث f'(x) هي المشتقة (الميل). ويُعطى طول القطعة الصغيرة جداً من المنحنى بنظرية فيثاغورس: ds = √((dx)² + (dy)²) = √(1 + (f'(x))²) dx. وبجمع كل هذه القطع المتناهية الصغر عبر التكامل نحصل على صيغة طول القوس:

L = ∫_a^b √(1 + (f'(x))²) dx

لكل جزء من هذه الصيغة معنى واضح. الرمز ∫_a^b يخبرنا بتجميع المساهمات من x = a إلى x = b. والتعبير f'(x) هو المشتقة، التي تقيس مدى انحدار المنحنى صعوداً عند كل نقطة. وتربيعها ثم إضافة 1 وأخذ الجذر التربيعي يحوّل الميل إلى عامل التمدد المحلي: أي مقدار زيادة طول المنحنى عن الخط الأفقي المستوي في ذلك الموضع. فحيث يكون المنحنى مستوياً، تكون f'(x) = 0 والعامل يساوي 1. وحيث يكون المنحنى شديد الانحدار، يكبر العامل.

تتطلب هذه الصيغة أن تكون الدالة أملسة بما يكفي: يجب أن تمتلك f مشتقة متصلة على الفترة، وأن تكون الدالة المُكاملة √(1 + (f'(x))²) قابلة للتكامل. وبالنسبة لمعظم الدوال التي تُصادف في مقررات التفاضل والتكامل والأعمال الهندسية، تتحقق هذه الشروط. وتنشأ المشاكل مع الزوايا الحادة، والمماسات الرأسية، وعدم الاستمرارية، حيث تسوء سلوك المشتقة.

الصعوبة هي أن هذا التكامل غالباً لا يملك دالة أصلية أولية. فحتى بالنسبة لدوال بسيطة مثل y = x²، تتضمن الدالة الأصلية دوال زائدية عكسية. وبالنسبة للعديد من المنحنيات الواقعية، لا توجد دالة أصلية مغلقة الصيغة على الإطلاق. وهنا بالضبط تصبح الطرق العددية مثل قاعدة سيمبسون ضرورية: فهي تُقدّر التكامل مباشرة دون الحاجة إلى دالة أصلية.

يتضمن تطبيق قاعدة سيمبسون على طول القوس تسلسلاً واضحاً من الخطوات. والفكرة الأساسية هي أن قاعدة سيمبسون لا تُطبق على الدالة الأصلية f(x)، بل تُطبق على الدالة المُكاملة لطول القوس g(x) = √(1 + (f'(x))²). وهذا التمييز حاسم ومصدر شائع للأخطاء.

إليك العملية الكاملة:

1. ابدأ بالدالة. حدّد y = f(x) وتأكد من أنها قابلة للاشتقاق على فترتك.
2. احسب المشتقة. أوجد f'(x) تحليلياً أو تأكد من توفرها.
3. ابنِ الدالة المُكاملة. كوّن g(x) = √(1 + (f'(x))²). هذه هي الدالة التي ستُكاملها قاعدة سيمبسون.
4. حدّد الفترة. تأكد من الحد الأدنى a والحد الأعلى b لقطعة المنحنى التي تريد قياسها.
5. اختر التقسيمات. اختر عدداً زوجياً n من الفترات الفرعية. يجب أن يكون العدد زوجياً لأن قاعدة سيمبسون تعمل في أزواج.
6. احسب حجم الخطوة. احسب h = (b - a) / n، وهو عرض كل فترة فرعية.
7. أوجد نقاط العينة. النقاط هي x₀ = a، وx₁ = a + h، وx₂ = a + 2h، وهكذا حتى x_n = b.
8. قيّم الدالة المُكاملة. احسب g(x₀)، وg(x₁)، ...، وg(x_n) عند كل نقطة عينة.
9. طبّق الأوزان. اضرب بنمط سيمبسون: 1 لنقطتي الطرفين، و4 للنقاط الداخلية ذات المؤشرات الفردية، و2 للنقاط الداخلية ذات المؤشرات الزوجية.
10. اجمع وقِس. اجمع القيم الموزونة واضرب في h/3 للحصول على تقدير طول القوس.

صيغة قاعدة سيمبسون 1/3 المركبة لطول القوس هي:

L ≈ (h/3)[g(x₀) + 4g(x₁) + 2g(x₂) + 4g(x₃) + ... + 4g(x_n-1) + g(x_n)]

حيث h = (b - a)/n وg(x) = √(1 + (f'(x))²).

يأتي شرط أن يكون n زوجياً من بنية الازدواج. فكل قطع مكافئ يمتد على فترتين فرعيتين (ثلاث نقاط). ولو كان n فردياً، لبقيت الفترة الفرعية الأخيرة دون شريك ولتعطلت الطريقة. وعندما تحتاج لسبب ما إلى عدد فردي من الألواح، تُستخدم قاعدة شبه المنحرف أو نهج مختلط بدلاً من ذلك.

لاحظ نمط الأوزان بعناية. القيمتان الأولى والأخيرة (g(x₀) وg(x_n)) تظهران مرة واحدة بالوزن 1. والنقاط الداخلية ذات المؤشرات الفردية (x₁، وx₃، وx₅، ...) تحصل على الوزن 4 لأن كل منها يعمل كنقطة منتصف للوحه القطعي المكافئ، حيث يكون تأثير القطع المكافئ أقوى. والنقاط الداخلية ذات المؤشرات الزوجية (x₂، وx₄، ...) تحصل على الوزن 2 لأن كل منها يعمل كنقطة طرف مشتركة بين لوحين متجاورين.

تتكون صيغة قاعدة سيمبسون المركبة من عدة عناصر يؤدي كل منها دوراً محدداً. وفهمها واحداً واحداً يجعل الطريقة أقل إرهاباً بكثير.

حجم الخطوة (h). هو عرض كل فترة فرعية، ويُحسب h = (b - a)/n. كلما صغُر h زاد عدد الألواح وتحسنت الدقة عموماً، على حساب مزيد من تقييمات الدالة. وفي طول القوس، يتحكم h في مدى دقة أخذ عينات المنحنى.

حدا التكامل (a وb). يحددان قطعة المنحنى المقاسة. ويؤثر اختيار الفترة تأثيراً كبيراً في النتيجة: فتمديد b أبعد على طول المنحنى يضيف مزيداً من الطول. تحقق دائماً من أن a < b وأن كليهما ضمن مجال الدالة.

نقاط العينة. النقاط x₀ حتى x_n وعددها n+1 متساوية التباعد عبر [a, b]. تُقيّم قاعدة سيمبسون الدالة المُكاملة عند كل نقطة من هذه النقاط دون استثناء.

مساهمتا الطرفين. تدخل القيمتان g(x₀) وg(x_n) المجموع مرة واحدة بالوزن 1. وهما تُثبّتان التقريب عند الحدود.

المساهمات ذات المؤشرات الفردية. تحصل النقاط x₁، وx₃، وx₅، وهكذا على الوزن 4. وهي نقاط منتصف كل لوح قطعي مكافئ. ويعكس الوزن 4 حقيقة أن نقطة المنتصف لها أكبر تأثير على المساحة تحت القطع المكافئ.

المساهمات الداخلية ذات المؤشرات الزوجية. تحصل النقاط x₂، وx₄، وx₆، وهكذا (باستثناء نقطتي الطرفين) على الوزن 2. فكل نقطة من هذه النقاط مشتركة كنقطة طرف بين لوحين متجاورين، لذا تُضاعف مساهمتها الكلية.

المجموع الموزون. اضرب كل قيمة دالة في وزنها، واجمع كل شيء معاً، ثم اضرب في h/3. فالعامل h/3 يضبط المجموع الموزون إلى الوحدات والمقدار الصحيحين.

يوضح الجدول التالي نمط الأوزان عند n = 6:

تُقدّر الحاسبة الموجودة في هذه الصفحة طول القوس باستخدام قاعدة سيمبسون عبر عدة أوضاع إدخال. يتعامل كل وضع مع طريقة مختلفة لوصف المنحنى، وفهم الوضع المناسب لمسألتك هو الخطوة الأولى.

وضع الدالة. أدخل صيغة مثل x^2 أو sin(x) أو exp(x)، بالإضافة إلى حدي الفترة a وb. تُشتق الحاسبة عددياً أو تستخدم البنية التي تقدمها، وتبني الدالة المُكاملة لطول القوس، وتطبق قاعدة سيمبسون. وهذا هو الوضع الأكثر شيوعاً لمسائل التفاضل والتكامل.

وضع نصف القطر والزاوية. بالنسبة للأقواس الدائرية، أدخل نصف القطر والزاوية المركزية. يتبع طول القوس مباشرة من الهندسة، ولا حاجة لقاعدة سيمبسون في الحالة الدائرية نفسها، لكن الوضع مفيد للمقارنة.

الوضع البارامتري. أدخل x(t) وy(t) كدالتين للبارامتر t، بالإضافة إلى فترة البارامتر. تستخدم الحاسبة صيغة طول القوس البارامترية مع تطبيق قاعدة سيمبسون على √((dx/dt)² + (dy/dt)²).

الوضع القطبي. أدخل r(θ) والفترة الزاوية. تُكامل قاعدة سيمبسون √(r² + (dr/dθ)²) على مدى الزاوية.

وضع النقاط. أدخل قائمة من أزواج الإحداثيات. تجمع الحاسبة المسافات المستقيمة بين النقاط المتتالية. وهذا تقريب مضلّع وليس قاعدة سيمبسون، وهو مفيد عندما تكون لديك بيانات مقاسة بدلاً من صيغة.

الوضع البارامتري ثلاثي الأبعاد. أدخل x(t) وy(t) وz(t) لمنحنى في الفضاء. تُكامل قاعدة سيمبسون √((dx/dt)² + (dy/dt)² + (dz/dt)²).

لتشغيل حساب طول قوس بقاعدة سيمبسون في وضع الدالة:

1. أدخل الدالة. اكتب الصيغة بالترميز المعياري، مثل x^2 أو sin(x).
2. حدّد الحد الأدنى. أدخل قيمة x الابتدائية a لقطعة المنحنى.
3. حدّد الحد الأعلى. أدخل قيمة x النهائية b. وتأكد من أن a < b.
4. اختر التقسيمات. اختر عدداً زوجياً n. ابدأ بقيمة معتدلة مثل 20 أو 50.
5. أكّد الطريقة. هذه الصفحة تُثبّت الطريقة العددية على قاعدة سيمبسون، فلا حاجة لأي اختيار.
6. شغّل الحساب. تُقيّم الأداة الدالة المُكاملة عند n+1 نقطة وتُرجع التقدير.
7. تحقق من التقارب. أعد التشغيل بقيمة n أكبر (مثلاً ضاعفها) وقارن. إذا بالكاد تغيرت النتيجة، فهي موثوقة.

للحصول على مثال واقعي، جرّب f(x) = x² على [0, 1] مع n = 20 ثم n = 40. يجب أن تتفق النتيجتان حتى عدة منازل عشرية، مما يؤكد استقرار التقدير.

يتناول هذا القسم مثالاً كاملاً مع عرض كل خطوة وسيطة. الدالة هي y = x² على الفترة [0, 1] باستخدام n = 4 فترات فرعية. جميع القيم محسوبة بشكل مستقل ويمكن التحقق منها بأي حاسبة.

الخطوة 1: الدالة. نبدأ بـ y = x²، وهو قطع مكافئ بسيط.

الخطوة 2: المشتقة. بالاشتقاق نحصل على y' = 2x.

الخطوة 3: الدالة المُكاملة لطول القوس. بالتعويض في الصيغة: g(x) = √(1 + (2x)²) = √(1 + 4x²).

الخطوة 4: التكامل. طول القوس الدقيق هو L = ∫₀¹ √(1 + 4x²) dx. لهذا التكامل صيغة مغلقة تتضمن دوال زائدية عكسية، وقيمته تقريباً 1.4789428575. وسنرى مدى اقتراب قاعدة سيمبسون منه.

الخطوة 5: حجم الخطوة. مع n = 4 على [0, 1]، يكون h = (1 - 0)/4 = 0.25.

الخطوة 6: نقاط العينة. النقاط الخمس متساوية التباعد هي:

• x₀ = 0.00
• x₁ = 0.25
• x₂ = 0.50
• x₃ = 0.75
• x₄ = 1.00

الخطوة 7: قيم الدالة. قيّم g(x) = √(1 + 4x²) عند كل نقطة:

• g(0.00) = √(1 + 0) = 1.0000000000
• g(0.25) = √(1 + 0.25) = √1.25 ≈ 1.1180339887
• g(0.50) = √(1 + 1) = √2 ≈ 1.4142135624
• g(0.75) = √(1 + 2.25) = √3.25 ≈ 1.8027756377
• g(1.00) = √(1 + 4) = √5 ≈ 2.2360679775

الخطوة 8: طبّق الأوزان. نمط سيمبسون عند n = 4 هو 1، 4، 2، 4، 1:

• 1 × g(x₀) = 1 × 1.0000000000 = 1.0000000000
• 4 × g(x₁) = 4 × 1.1180339887 = 4.4721359550
• 2 × g(x₂) = 2 × 1.4142135624 = 2.8284271247
• 4 × g(x₃) = 4 × 1.8027756377 = 7.2111025508
• 1 × g(x₄) = 1 × 2.2360679775 = 2.2360679775

الخطوة 9: المجموع الموزون. بجمع هذه القيم نحصل على 1.0000000000 + 4.4721359550 + 2.8284271247 + 7.2111025508 + 2.2360679775 = 17.7477336082.

الخطوة 10: النتيجة النهائية. اضرب في h/3 = 0.25/3 ≈ 0.0833333333: L ≈ 0.0833333333 × 17.7477336082 ≈ 1.4789778007.

بالمقارنة مع القيمة الدقيقة 1.4789428575، يبلغ الخطأ حوالي 0.0000349، أي تقريباً 3.5 × 10^-5. ومع 4 فترات فرعية فقط (5 تقييمات للدالة)، تحقق قاعدة سيمبسون دقة حتى 4 منازل عشرية. هذه هي قوة الطريقة: تقارب سريع على الدوال المُكاملة الأملسة.

للمقارنة، المضاعفة إلى n = 6 تعطي 1.4789405084 (خطأ ≈ 2.4 × 10^-6)، وحتى التقسيم الخشن n = 2 يعطي 1.4821537045 (خطأ ≈ 0.0032). كل مضاعفة لـ n تُقلل الخطأ بشكل كبير، وهذه سمة التقارب من الرتبة الرابعة لقاعدة سيمبسون.

كل طريقة عددية تُنتج تقريباً لا إجابة دقيقة. والفرق بين التقريب والقيمة الحقيقية يُسمى الخطأ. وفهم سلوك هذا الخطأ يساعدك على اختيار الإعدادات المناسبة والوثوق بنتائجك.

قاعدة سيمبسون طريقة من الرتبة الرابعة، أي أن الخطأ يتضاءل تناسبياً مع h⁴ (أو بالمثل 1/n⁴) كلما صغُر حجم الخطوة. وعملياً، مضاعفة عدد الفترات الفرعية تُقلل الخطأ عادة بعامل يقارب 16. ولهذا تحسّن المثال المحلول أعلاه بشكل كبير من n = 2 إلى n = 4 إلى n = 6.

حد الخطأ الكلاسيكي لقاعدة سيمبسون المركبة هو:

|E_s| ≤ M(b - a)⁵ / (180n⁴)

هنا، M هو حد أعلى للقيمة المطلقة للمشتقة الرابعة للدالة المُكاملة g(x) على الفترة [a, b]، أي M ≥ |g⁽⁴⁾(x)| لكل x في [a, b]. وتقيس المشتقة الرابعة مدى سرعة تغير انحناء الدالة المُكاملة نفسها. فعندما تكون g أملسة جداً (مشتقتها الرابعة صغيرة)، يضمن الحد خطأ صغيراً. وعندما تتذبذب g بسرعة، تكون M كبيرة والحد أقل طمأنة.

تنطبق عدة تحفظات مهمة. أولاً، يتطلب هذا الحد أن تمتلك g مشتقة رابعة متصلة على الفترة. فإذا كانت الدالة المُكاملة ذات انحناء مفاجئ أو نتوء أو انقطاع في مشتقة عليا، لا ينطبق الحد. ثانياً، حساب M عملياً غالباً ما يكون صعباً، لذا فالحد أنفع للفهم النظري منه لتقدير الخطأ اليومي. ثالثاً، قد يكون الحد متشائماً جداً: فالأخطاء الفعلية غالباً أصغر بكثير مما تتنبأ به الصيغة.

وعملياً، فإن أكثر فحوصات الدقة موثوقية هو اختبار التقارب. احسب طول القوس بعدد n من التقسيمات، ثم بـ 2n، ثم بـ 4n. فإذا اتفقت النتائج حتى عدد المنازل العشرية التي تحتاجها، كان التقدير جديراً بالثقة. وإذا استمرت في التغير بشكل ملحوظ، فزِد n أكثر أو تحقق مما إذا كانت الدالة المُكاملة ذات سلوك إشكالي (انعطافات حادة، أو شبه تفردات، أو تذبذب سريع).

نعومة الدالة الأصلية f مهمة بشكل غير مباشر عبر g. فإذا كانت f ذات زاوية حادة، تقفز f' بشكل غير متصل، وترث g = √(1 + (f')²) انحناء مفاجئاً. وتظل قاعدة سيمبسون تعمل، لكن التقارب يتباطأ. وفي مثل هذه الحالات، فإن تقسيم الفترة عند النقطة الإشكالية وتطبيق قاعدة سيمبسون على كل قطعة أملسة على حدة يعيد التقارب السريع.

قاعدة شبه المنحرف هي أقرب أقارب قاعدة سيمبسون وأكثر بدائلها شيوعاً. وكلتاهما من صيغ نيوتن-كوتس التي تستخدم نقاطاً متساوية التباعد، لكنهما تختلفان جوهرياً في كيفية تقريب الدالة المُكاملة بين نقاط العينة.

تربط قاعدة شبه المنحرف نقاط العينة المتتالية بخطوط مستقيمة مكونة أشباه منحرفات. أما قاعدة سيمبسون فتُركّب قطعاً مكافئة عبر مجموعات من ثلاث نقاط. وبما أن القطع المكافئ يمكن أن ينحني ليتبع الدالة المُكاملة بينما لا يستطيع الخط المستقيم ذلك، فإن قاعدة سيمبسون تلتقط الانحناء الذي تفتقده قاعدة شبه المنحرف. وهذا الفرق الهندسي يترجم مباشرة إلى دقة: قاعدة سيمبسون من الرتبة الرابعة (خطأ يتناسب مع 1/n⁴) بينما قاعدة شبه المنحرف من الرتبة الثانية (خطأ يتناسب مع 1/n²).

وعملياً، هذا يعني أن قاعدة سيمبسون تحقق غالباً مع n = 20 ما تحتاج قاعدة شبه المنحرف n = 200 أو أكثر لمطابقته، على الأقل للدوال المُكاملة الأملسة. والثمن هو أن قاعدة سيمبسون تتطلب n زوجياً ودفتر حسابات أدق قليلاً مع الأوزان المتناوبة. وتعمل قاعدة شبه المنحرف مع أي n ولها نمط أوزان أبسط (1، 2، 2، ...، 2، 1).

ولا تتفوق أي من الطريقتين عالمياً. لقاعدة شبه المنحرف ميزة مفاجئة مع الدوال الدورية المُكاملة على فترات كاملة، حيث يمكن أن تتقارب بسرعة أسية. وقد تواجه قاعدة سيمبسون صعوبة مع الدوال غير الأملسة بما يكفي لتطبيق نظريتها من الرتبة الرابعة. وبالنسبة للبيانات الخشنة أو المليئة بالضجيج، تكون قاعدة شبه المنحرف الأبسط أحياناً أكثر متانة.

سير العمل السليم هو استخدام قاعدة سيمبسون كطريقة أساسية والتحقق المتقاطع مع قاعدة شبه المنحرف بدقة أعلى عندما يكون سلوك الدالة المُكاملة غير مؤكد. فإذا اتفقت الطريقتان، كانت الثقة عالية. وإذا اختلفتا بشكل ملحوظ، فمن المرجح أن الدالة المُكاملة تحتاج مزيداً من التقسيمات أو تقسيم الفترة.

يمكن وصف المنحنيات في عدة أنظمة إحداثيات، ولكل منها صيغة طول قوس خاصة. وتنطبق قاعدة سيمبسون على جميعها لأن كل صيغة هي تكامل محدد لدالة مُكاملة أملسة.

طول القوس الديكارتي

بالنسبة لـ y = f(x) على [a, b]، الصيغة هي L = ∫_a^b √(1 + (dy/dx)²) dx. وهذه هي الصيغة المعيارية التي تُدرّس في مقررات التفاضل والتكامل. تُكامل قاعدة سيمبسون g(x) = √(1 + (f'(x))²) كما وُصف في هذا الدليل. وهي تعمل جيداً متى كانت f قابلة للاشتقاق بمشتقة متصلة.

طول القوس البارامتري

بالنسبة لمنحنى موصوف بـ x = x(t) وy = y(t) مع t من a إلى b، الصيغة هي L = ∫_a^b √((dx/dt)² + (dy/dt)²) dt. يتعامل الشكل البارامتري مع المنحنيات التي تفشل في اختبار الخط الرأسي، مثل الدوائر والقطوع الناقصة والدويريات، التي لا يمكن كتابتها على الصورة y = f(x). تُكامل قاعدة سيمبسون دالة السرعة البارامترية √((x'(t))² + (y'(t))²)، التي تمثل سرعة تحرك نقطة على طول المنحنى مع زيادة t.

طول القوس القطبي

بالنسبة لـ r = r(θ) مع θ من α إلى β، الصيغة هي L = ∫_α^β √(r² + (dr/dθ)²) dθ. تصف الإحداثيات القطبية اللوالب والقلبيات ومنحنيات الوردة بشكل طبيعي. تجمع الدالة المُكاملة بين المسافة الشعاعية r ومعدل التغير الزاوي dr/dθ. وتنطبق قاعدة سيمبسون مباشرة، مُكاملةً بالنسبة إلى θ.

طول القوس ثلاثي الأبعاد

بالنسبة لمنحنى في الفضاء x = x(t)، وy = y(t)، وz = z(t)، تمتد الصيغة بشكل طبيعي: L = ∫_a^b √((dx/dt)² + (dy/dt)² + (dz/dt)²) dt. يساهم كل إحداثي بمربع مركبة سرعته. وهذا ضروري للمسارات الحلزونية والمسارات ثلاثية الأبعاد والمنحنيات في رسوميات الحاسوب. تتعامل قاعدة سيمبسون مع الدالة المُكاملة ثلاثية الأبعاد بنفس طريقة النسخ ثنائية الأبعاد.

في كل حالة، الاستراتيجية متطابقة: حدّد الدالة المُكاملة، وتأكد من أنها أملسة على الفترة، واختر n زوجياً، وطبّق المجموع الموزون. تدعم الحاسبة في هذه الصفحة أوضاع الدالة والبارامتري والقطبي والبارامتري ثلاثي الأبعاد، ويمر كل منها عبر قاعدة سيمبسون بشكل مناسب.

تظهر قاعدة سيمبسون لطول القوس في العديد من المجالات حيث يجب تقدير المسافات المنحنية من نماذج رياضية أو بيانات مقاسة.

الهندسة الميكانيكية. تُصمم منحنيات الكامات وأسنان التروس وقضبان التوجيه كمنحنيات رياضية أملسة. يساعد تقدير طولها في تحديد متطلبات المواد ومسافات التلامس وأنماط التآكل. وتوفر قاعدة سيمبسون فحوصات عددية سريعة أثناء دورة تكرار التصميم.

هندسة الطرق والنقل. تتبع منحنيات الطرق السريعة وتحولات السكك الحديدية وتخطيطات المدارج ملامح هندسية مصممة بعناية (اللوالب الانتقالية والأقواس الدائرية والمنحنيات الرأسية القطعية المكافئة). حساب الطول الحقيقي لهذه المحاور مهم لتقدير التكلفة ومواضع اللافتات وتخطيط الإنشاء.

التصميم بمساعدة الحاسوب (CAD). تمثل برامج التصميم المنحنيات كشرائح ومنحنيات بيزيه. يُقدّر التكامل العددي طول هذه المنحنيات المركبة لأغراض التصنيع، حيث تعتمد مسارات القطع واستخدام المواد على بيانات طول دقيقة.

الروبوتات وتخطيط الحركة. يتبع ذراع الروبوت أو المركبة ذاتية القيادة مساراً مخططاً عبر الفضاء. يؤثر طول المسار الكلي في زمن السفر واستهلاك الطاقة وتخطيط المهمة. تُكامل قاعدة سيمبسون ملف السرعة على طول المسارات البارامترية.

الحوسبة العلمية. يحتاج الباحثون الذين يصممون نماذج للظواهر الفيزيائية، من مسارات الجسيمات إلى خطوط المجال، بشكل روتيني إلى أطوال المنحنيات. وعندما لا تُنتج المعادلات الحاكمة طول قوس بصيغة مغلقة، يكون التكامل العددي هو النهج المعياري.

الهندسة الإنشائية. للجسور المقوسة والعوارض المنحنية والمنشآت القشرية أطوال ومقاييس سطح مشتقة من تعريفاتها الهندسية. ويدعم حساب طول القوس الدقيق تحليل الأحمال وحصر المواد.

النمذجة الرياضية والبحث. في الرياضيات البحتة والتطبيقية، يظهر طول القوس في الهندسة التفاضلية والمسائل التغايرية والاستمثال. وتوفر التقديرات العددية حدساً وتحققاً للنتائج النظرية.

في كل هذه السياقات، تُقدَّر قاعدة سيمبسون لأنها تقدم دقة عالية بعدد قليل نسبياً من تقييمات الدالة، مما يجعلها فعالة لكل من الأدوات التفاعلية والحوسبة واسعة النطاق.

عدة أخطاء متكررة تُوقع الطلاب والممارسين. والوعي بهذه المزالق يمنع معظم النتائج السيئة.

استخدام عدد فردي من التقسيمات. تتطلب قاعدة سيمبسون 1/3 جوهرياً n زوجياً. مع n فردي تنكسر بنية الازدواج. تحقق دائماً من أن n زوجي قبل الحساب. وإذا اقترحت مسألتك عدداً فردياً، قرّبه إلى أقرب عدد زوجي أعلى.

مشتقة خاطئة. تستخدم الدالة المُكاملة f'(x) لا f(x). ومن الأخطاء الشائعة الاشتقاق بشكل غير صحيح (مثل كتابة مشتقة x² على أنها x بدلاً من 2x) أو نسيان قاعدة السلسلة للدوال المركبة. تحقق من كل مشتقة مرتين.

تطبيق قاعدة سيمبسون على f بدلاً من g. ومما يتصل بما سبق، أن بعض المستخدمين يطبقون الأوزان مباشرة على قيم الدالة الأصلية f(x_i) بدلاً من الدالة المُكاملة لطول القوس g(x_i) = √(1 + (f'(x_i))²). وهذا يحسب تكاملاً خاطئاً تماماً.

حدود فترة غير صحيحة. إدخال a > b، أو حدود خارج مجال الدالة، يُنتج نتائج لا معنى لها. وبالنسبة للدوال ذات المجالات المقيدة (مثل √(x) التي تتطلب x ≥ 0)، تحقق من صلاحية الفترة.

تقسيمات قليلة جداً. مع أن قاعدة سيمبسون تتقارب بسرعة، قد لا تكفي n = 2 أو n = 4 للدوال المُكاملة شديدة الانحناء. أجرِ دائماً فحص تقارب بقيم n أكبر.

الخلط بين الدرجات والراديان. عندما تتضمن f دوال مثلثية، تأكد مما إذا كانت الحاسبة تتوقع الراديان أم الدرجات. فخلطهما يُفسد المشتقة والنتيجة بصمت. معظم الأدوات الرياضية تستخدم الراديان.

افتراض الدقة التامة. قاعدة سيمبسون دقيقة تماماً فقط لكثيرات الحدود حتى الدرجة 3. وبالنسبة للدوال المُكاملة العامة، النتيجة تقريب. ومعاملتها على أنها دقيقة دون دليل تقارب خطأ.

تجاهل التفردات. إذا انفجرت f' داخل [a, b] (مثلاً قرب مماس رأسي)، تصبح الدالة المُكاملة g كبيرة وقد تتقارب قاعدة سيمبسون ببطء أو بشكل غير موثوق. قسّم الفترة أو استخدم طرقاً متخصصة.

قاعدة سيمبسون خيار افتراضي ممتاز للحساب العددي لطول القوس، لكن معرفة متى تتألق ومتى تبحث عن بديل تجعلك حلال مسائل أكثر فعالية.

استخدم قاعدة سيمبسون عندما:

• تكون الدالة المُكاملة g(x) = √(1 + (f'(x))²) أملسة على الفترة، دون انحناءات مفاجئة أو تفردات.
• لا توجد دالة أصلية أولية، أو كان إيجادها جهداً أكبر مما يستحقه تقدير عددي.
• تحتاج دقة جيدة بعدد معتدل من تقييمات الدالة.
• تحسب طول القوس لأعمال دراسية أو فحوصات هندسية أو التحقق من التصميم.
• يكون المنحنى موصوفاً بارامترياً أو بإحداثيات قطبية وتكون الدالة المُكاملة الناتجة أملسة.

فكّر في بدائل أو عناية إضافية عندما:

• تكون الدالة المُكاملة ذات انقطاع أو زاوية حادة داخل الفترة. قسّم الفترة عند النقطة الإشكالية أولاً.
• تتذبذب الدالة بسرعة. قد تحتاج n كبيراً جداً، أو طريقة مصممة للتكاملات المتذبذبة.
• لديك بيانات مقاسة مليئة بالضجيج فقط بدلاً من صيغة أملسة. قد تكون قاعدة شبه المنحرف أو النهج القائمة على الشرائح أنسب.
• تُطلب دقة عالية جداً. قد تكون التربيع الغاوسي أو قاعدة سيمبسون التكيفية (التي تُحسّن التقسيمات تلقائياً حيث الحاجة) أفضل.
• تكون الفترة لانهائية أو شبه لانهائية. حوّل التكامل إلى فترة محدودة أولاً، أو استخدم تربيعاً متخصصاً.

بالنسبة لمعظم طلاب التفاضل والتكامل والمهندسين العاملين، توفر قاعدة سيمبسون مع فحص تقارب (مقارنة n و2n و4n) كل الدقة اللازمة. فهي سهلة التطبيق، ومتوفرة على نطاق واسع في البرامج، ومفهومة نظرياً بشكل جيد.