قاعدة شبه المنحرف

ملف:Trapezoidal rule illustration.png
الدالة f(x) (باللون الأزرق) تم تقريبها بدالة خطية (باللون الأحمر).
ملف:Composite trapezoidal rule illustration.png
توضيح لقاعدة شبه المنحرف المركبة (بتشبيك غير منتظم).
ملف:Trapezoidal rule illustration small.svg
شكل توضيحي لقاعدة شبه المنحرف (بتشبيك منتظم).

في الرياضيات، تعتبر قاعدة شبة المنحرف إحدى طرق الحساب التقريبي للتكامل المحدد.

abf(x)dx.

تعمل قاعدة شبه المنحرف بتقريب المنطقة تحت منحنى الدالة f(x) بشبه منحرف وحساب مساحته. ينجم عن ذلك

abf(x)dx(ba)f(a)+f(b)2.

لحساب التكامل بدقة أفضل, يمكن فصل فترة التكامل [a,b] أولا إلىn فترات أصغر, ومن ثم تطبيق قاعدة شبه المنحرف على كل فترة. يمكن تحصيل قاعدة شبه المنحرف المركب:

abf(x)dxban[f(a)+f(b)2+k=1n1f(a+kban)].

ويمكن صياغة هذا بشكل اخر:

abf(x)dxba2n(f(x0)+2f(x1)+2f(x2)++2f(xn1)+f(xn))

حيث

xk=a+kban, for k=0,1,,n

تحليل الخطأ

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

error=abf(x)dxban[f(a)+f(b)2+k=1n1f(a+kban)].

يمكن كتابة هذا الخطأ بالشكل

error=(ba)312n2f(ξ),

حيثξ عدد ما بين a وb.[١]

يعطى تخمين الخطأ المقارب لـ n → ∞ بالعلاقة

error=(ba)212n2(f(b)f(a))+O(n3). [٢]

الحدود الأخرى لهذا الخطأ يمكن إيجادها من صيغة مجموع أويلر-ماكلورين.

البرمجة

مثال على قاعدة شبه المنحرف مكتوب بلغة البايثون <source lang="python">

  1. !/usr/bin/env python

def trapezoidal_rule(f, a, b, N):

   """Approximate the definite integral of f from a to b by the
   composite trapezoidal rule, using N subintervals"""
   return (b-a) * (f(a)/2 + f(b)/2 + sum([f(a + (b-a)*k/N) for k in range(1,N)])) / N
  1. test

print trapezoidal_rule(lambda x:x**9, 0.0, 10.0, 100000) </source>

إنظر أيضا

ملاحظات

  1. ^ Atkinson (1989), equation (5.1.7)
  2. ^ Atkinson (1989), equation (5.1.9)

مراجع

ca:Mètode trapezial de:Trapezregel Trapezoidal rule]] es:Regla del trapecio fi:Puolisuunnikassääntö fr:Méthode des trapèzes hu:Trapézszabály it:Regola del trapezio ja:台形公式 ko:사다리꼴 공식 nl:Trapeziumregel pl:Wzór trapezów ru:Метод трапеций sr:Трапезоидно правило sv:Trapetsmetoden uk:Метод трапецій (математика) zh:梯形公式