قاعدة سمبسون
'''قاعدة سمبسون''' {{إنج|Simpson's rule}} في [[تحليل عددي|التحليل العددي]] هي طريقة من طرق [[تكامل عددي|التكامل العددي]] و هي في الحقيقة حالة خاصة من [[صيغ نيوتن-كوتس]] المغلقة لتقريب تكامل الدالة f باستخدام [[متعددة الحدود|كثيرة الحدود]] التربيعية وهي طريقة محسنة لطريقة شبه المنحرف كما أنها أسرع تقارباً وأدق ويفسر ذلك من خلال أن قاعدة سمبسون تحتوي على نقطة المنتصف التي توفر توازن أفضل للتقريب. لأنه كلما زادت عدد التقسيمات في الفترة الجزئية كانت الطريقة أدق.
و تحسب صيغة [[قاعدة شبه المنحرف|شبه المنحرف]] القيمة الفعلية للتكامل عندما تكون f دالة كثيرة حدود من الدرجة الأولى على الأكثر . بينما صيغة سمبسون فإنها تحسب القيمة الفعلية للتكامل إذا كانت f دالة [[متعددة الحدود|كثيرة الحدود]] من الدرجة الثالثة أو أقل.<ref name="Richard L. Burden 2011">“Numerical Analysis”, Richard L. Burden & J. Douglas Faires, Brooks/Cole, Cengage Learning,(2011) ISBN 0-534-38216-9</ref>

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

[[ملف:simpsons method illustration.svg|تصغير|يمين| يمكن أن تستمد قاعدة سمبسون بتقريب المساحة المطلوبة تحت المنحنى الدالة f باللون الأزرق بالمساحة تحت منحنى حدودية من الدرجة الثانية لP باللون الأحمر]]

== استنتاج القاعدة ==
يمكن استنتاج قاعدة سمبسون بطرق مختلفة.
=== الاستكمال التربيعي ===
 
ولكي نحصل على صيغة سمبسون لابد بالتعويض في الصيغة الإستكمالية

<math>P_2(x) = c_\circ+c_1x +c_2x^2</math>

<math>\int_a^bf(x)dx = {\int_{x_\circ}^{x_2}}[\frac{(x-x_1)(x-x_2)}{(x_\circ-x_1)(x_\circ-x_2)}f(x_\circ)+\frac{(x-x_\circ)(x-x_2)}{(x_1-x_\circ)(x_1-x_2)}f(x_1)+\frac{(x-x_\circ)(x-x_1)}{(x_2-x_\circ)(x_2-x_1)}f(x_2)]dx + {\int_{x_\circ}^{x_2}}\frac{(x-x_\circ)(x-x_1)(x-x_2)}{6}f^{(3)}(\xi(x))dx </math>

<math>f(x)= f(x_1)+f'(x_1)(x-x_1)+\frac{f''(x_1)}{2}(x-x_1)^2+\frac{f'''(x_1)}{6}(x-x_1)^3+\frac{f^{(4)}(\xi(x)}{24}(x-x_1)^4</math>

and

<math>\int_{x_\circ}^{x_2} f(x)dx -\left[f(x_1)(x-x_1)+\frac{f'(x_1)}{2}(x-x_1)^2+\frac{f''(x_1)}{6}(x-x_1)^3+\frac{f'''(\xi(x)}{24}(x-x_1)^4\right]_{x_\circ}^{x_2} + \frac{1}{24}\int_{x_\circ}^{x_2}f^{(4)}(\xi(x))(x-x_1)^{4}dx </math>

<math>\frac{1}{24}\int_{x_\circ}^{x_2}f^{(4)}(\xi(x))(x-x_1)^{4}dx = \frac{f^{(4)}(\xi_1)}{24}\int_{x_\circ}^{x_2}(x-x_1)^4dx = \left[\frac{f^{(4)}(\xi_1)}{120}(x-x_1)^5\right]_{x_\circ}^{x_2} </math>

لبعض قيم <math>\xi_1</math> في <math>(x_\circ,x_1)</math>

أيضا

<math>h=x_2 - x_1=x_1-x_\circ</math>

بالتالي :

<math>(x_2-x_1)^2-(x_\circ-x_1)^2=(x_2-x_1)^4-(x_\circ-x_1)^4=0</math>

<math>(x_2-x_1)^3-(x_\circ-x_1)^3=2h  and (x_2-x_1)^5-(x_\circ-x_1)^5=2h^5   </math>

<math>\int_{x_\circ}^{x_2} f(x)dx = 2hf(x_1)+\frac{h^3}{3}f''(x_1)+\frac{f^{(4)}(\xi_1)}{60}h^5</math>

وبتعويض قيمة <math>f''(x_1)</math> يصبح لدينا

<math>\int_{x_\circ}^{x_2} f(x)dx = 2hf(x_1)+\frac{h^3}{3}\left\{\frac{1}{h^2}[f(x_\circ)-2f(x_1)+f(x_2)]-\frac{h^2}{12}f^{(4)}(\xi_2)\right\}+\frac{f^{(4)}(\xi_1)}{60}h^5</math>

إذاً

<math>\int_{x_\circ}^{x_2} f(x)dx=\frac{h}{3}[f(x_\circ)+4f(x_1)+f(x_2)]-\frac{h^5}{12}\left[\frac{1}{3}f^{(4)}(\xi_2)-\frac{1}{5}f^{(4)}(\xi_1)\right]</math>
<ref name="Richard L. Burden 2011"/>

== قاعدة سمبسون البسيطة ==

صيغة سمبسون البسيطة :
الفكرة التي قامت عليها صيغة سمبسون البسيطة:

تقوم على أساس تقريب [[مساحة|المساحة]] المطلوبة تحت منحنى الدالة <math>f(x)</math> بالمساحة تحت منحنى حدودية من الدرجة الثانية <math>P_2(x)</math> تمر بالنقاط الثلاث

<math>(x_\circ,f(x_\circ),(x_1,f(x_1),(x_2,f(x_2)</math>

<math>=\frac{h}{3}[f(x_0)+4f(x_1)+f(x_2)]</math> <math>\, \int_{x_0}^{x_2} f(x)\, dx</math>

بحيث تكون المساحة محصورة بين الخطين <ref name="Richard L. Burden 2011"/>

<math>a=x_\circ , b=x_2</math>

=== مثال على سمبسون البسيطة ===
استخدام قاعدة سمبسون البسيطة لإيجاد قيمة التكامل <math>\, \int_{0}^{2} x^2\, dx</math>

الحل:

<math>\, \int_{x_0}^{x_2} f(x)\, dx</math>

<math>=\frac{h}{3}[f(x_0)+4f(x_1)+f(x_2)]</math>

<math>x_0 = 0</math>

<math>x_1 = 1</math>

<math>x_2 = 2</math>

<math>h= x_1 - x_0 = 1-0 =1</math>

<math> \to \, \int_{0}^{2} x^2\, dx =\frac{1}{3} [f(0)+4f(1)+f(2)]  </math>

<math>=\frac{1}{3}[(0)^2+4(1)^2+(2)^2]</math>

<math>=\frac{1}{3}[8] =  \frac{8}{3} </math>
<ref name="Richard L. Burden 2011"/>

== قاعدة سمبسون المركبة ==

قد لا يكون استخدام [[صيغ نيوتن-كوتس]] عمليا في الكثير من الحالات خصوصا إذا كانت فترة التكامل كبيرة نسبيا ؛ حيث أننا نضطر ففي مثل هذه الحالات إلى استخدام [[متعددة الحدود|كثيرات حدود]] ذات درجات عالية من الذبذبة وهذا بدوره يترك أثر سيئا على دقة الحلول العددية، وللتغلب على هذه المشكلة وتقليل الخطأ الناتج عن تطبيق صيغة تكامل ذات رتبة منخفضة أي قيمة n صغيرة فإننا نقسم فترة التكامل [a,b] إلى فترات أصغر ونطبق صيغة التكامل على كل فترة جزئية على حدة . والصيغة الناتجة من التطبيق المتكرر لصيغة ذات رتبة منخفضة تسمى صيغة تكامل مركبة . وفيما يلي نحصل على الصيغة المركبة الناتجة من التطبيق المتكرر لصيغة سمبسون البسيطة عندما (n=2) . ولكن من شروط قاعدة سمبسون ان تكون n عدد زوجي
وكي نكرر تطبيق قاعدة سمبسون عدد r من المرات نحتاج معرفة قيم الدالة عند عدد 2r+1 من نقاط الأساس <math>x_0,x_1,x_2.....x_2r</math>
<ref>كتاب الطرق العددية والتحليل العددي - أ.د. أبو بكر أحمد السيد، [[جامعة الكويت]] ISBN9789957171353</ref>

وتكون قاعدة سمبسون المركبة كالتالي :

:<math>\int_a^b f(x) \, dx\approx\tfrac{h}{3}\bigg[f(x_0)+2\sum_{j=1}^{n/2-1}f(x_{2j})+
4\sum_{j=1}^{n/2}f(x_{2j-1})+f(x_n)
\bigg]</math>

=== مثال على قاعدة سمبسون المركبة ===

<math>\int_{-0.5}^{0.5} xln(x+1)dx\, </math>

<math>n=6</math>

الحل بستخدام قاعدة سمبسون المركبة :

<math>x_0 = -0.5</math>

<math>x_1 = -0.333</math>

<math>x_2 = -0.166</math>

<math>x_3 = 0</math>

<math>x_4 = 0.166</math>

<math>x_5 = 0.333</math>

<math>x_6 = 0.5</math>

<math>h= \frac{0.5 -(-0.5)}{6}= \frac{1}{6}</math>

<math>\int_{-0.5}^{0.5} xln(x+1)dx\, = \frac{\frac{1}{6}}{3} [f(-0.5)+4f(-0.333)+2f(-0.166)+4f(0)+2f(0.166)+4f(0.333)+f(0.5)]</math>

<math>f(-0.5) = (-0.5)ln(-0.5+1) = 0.346</math>

<math>f(-0.333) = (-0.333)ln(-0.333+1) = 0.134</math>

<math>f(-0.166) = (-0.166)ln(-0.166+1) = 0.030</math>

<math>f(0) = (0)ln(0+1) = 0</math>

<math>f(0.166) = (0.166)ln(0.166+1) = 0.025</math>

<math>f(0.333) = (0.333)ln(0.333+1) = 0.095</math>

<math>f(0.5) = (0.5)ln(0.5+1) = 0.202</math>

<math>\int_{-0.5}^{0.5} xln(x+1)dx\,</math>

<math>=\frac{1}{18} [0.346+0.536+0.06+0+0.05+0.38+0.202]</math>

<math>=\frac{1.574}{18} = 0.0874</math>
<ref name="Richard L. Burden 2011"/>

== الخطأ ==

وفي بعض التكاملات على الفترات الكبيرة كان الحل بصيغة سمبسون يحتوي نسبة خطأ أكبر من لو قام بتقسيم الفترة

نلاحظ أن الخطأ يقل بعد تقسيم الفترة ومن هذا المنطلق نتجت صيغة سمبسون المركبة ولها العديد من الصيغ

<math>\int_a^bf(x)dx = \frac{h}{3}\left[ f(x_\circ)+2\sum_{j=1}^{(\frac{n}{2})-1}f(x_{2_j})+4\sum_{j=1}^{\frac{n}{2}}f(x_{2_{j-1}})+f(x_n) \right]-\frac{h^5}{90}\sum_{j=1}^{\frac{n}{2}}f^{(4)}(\xi_j)	 </math>

<math>h=\frac{b-a}{n}</math>

والخطأ مرتبط بالتقريب :

<math>E(f)=-\frac{h^5}{90}\sum_{j=1}^{\frac{n}{2}}f^{(4)}(\xi_j)</math>
<ref name="Richard L. Burden 2011"/>

حسب نظرية القيمة المتوسطة يوجد<math>\mu\in(a,b) </math> بحيث:

<math>f^{(4)}(\mu)=\frac{n}{2}\sum_{j=1}^{\frac{n}{2}}f^{(4)}(\xi_j)</math>

بالتالي

<math>E(f)=-\frac{h^5}{90}\sum_{j=1}^{\frac{n}{2}}f^{(4)}(\xi_j)= -\frac{h^5}{180}nf^{(4)}(\mu)</math>

وبما أن

<math>h=\frac{b-a}{n}</math>

و أيضاً

<math> E(f)=-\frac{b-a}{180}h^4nf^{(4)}(\mu)</math>

== انظر ايضًا ==
* [[قاعدة شبه المنحرف]]
* [[روجر كوتس]]
* [[إسحاق نيوتن]]
* [[تكامل عددي|التكامل العددي]]
* [[تحليل عددي|التحليل العددي]]

== مراجع ==

{{مراجع}}

* كتاب الطرق العددية والتحليل العددي - أ.د. أبو بكر أحمد السيد، جامعة الكويت ISBN9789957171353
* “Numerical Analysis”, Richard L. Burden & J. Douglas Faires, Brooks/Cole, Cengage Learning,(2011) ISBN 0-534-38216-9

== ملاحظات ==
{{ملاحظات}}
{{شريط بوابات|رياضيات|تحليل رياضي}}

[[تصنيف:تحليل عددي]]
[[تصنيف:تكامل عددي]]
[[تصنيف:حساب تكاملي]]
قاعدة سمبسون  في التحليل العددي هي طريقة من طرق التكامل العددي و هي في الحقيقة حالة خاصة من صيغ نيوتن-كوتس المغلقة لتقريب تكامل الدالة f باستخدام كثيرة الحدود التربيعية وهي طريقة محسنة لطريقة شبه المنحرف كما أنها أسرع تقارباً وأدق ويفسر ذلك من خلال أن قاعدة سمبسون تحتوي على نقطة المنتصف التي توفر توازن أفضل للتقريب. لأنه كلما زادت عدد التقسيمات في الفترة الجزئية كانت الطريقة أدق.
و تحسب صيغة شبه المنحرف القيمة الفعلية للتكامل عندما تكون f دالة كثيرة حدود من الدرجة الأولى على الأكثر . بينما صيغة سمبسون فإنها تحسب القيمة الفعلية للتكامل إذا كانت f دالة كثيرة الحدود من الدرجة الثالثة أو أقل.“Numerical Analysis”, Richard L. Burden & J. Douglas Faires, Brooks/Cole, Cengage Learning,(2011) ISBN 0-534-38216-9

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

تصغير|يمين| يمكن أن تستمد قاعدة سمبسون بتقريب المساحة المطلوبة تحت المنحنى الدالة f باللون الأزرق بالمساحة تحت منحنى حدودية من الدرجة الثانية لP باللون الأحمر

 استنتاج القاعدة 
يمكن استنتاج قاعدة سمبسون بطرق مختلفة.
 الاستكمال التربيعي 
 
ولكي نحصل على صيغة سمبسون لابد بالتعويض في الصيغة الإستكمالية

and

لبعض قيم  في 

أيضا

بالتالي :

وبتعويض قيمة  يصبح لدينا

إذاً

 قاعدة سمبسون البسيطة 

صيغة سمبسون البسيطة :
الفكرة التي قامت عليها صيغة سمبسون البسيطة:

تقوم على أساس تقريب المساحة المطلوبة تحت منحنى الدالة  بالمساحة تحت منحنى حدودية من الدرجة الثانية  تمر بالنقاط الثلاث

 

بحيث تكون المساحة محصورة بين الخطين 

 مثال على سمبسون البسيطة 
استخدام قاعدة سمبسون البسيطة لإيجاد قيمة التكامل 

الحل:

 قاعدة سمبسون المركبة 

قد لا يكون استخدام صيغ نيوتن-كوتس عمليا في الكثير من الحالات خصوصا إذا كانت فترة التكامل كبيرة نسبيا ؛ حيث أننا نضطر ففي مثل هذه الحالات إلى استخدام كثيرات حدود ذات درجات عالية من الذبذبة وهذا بدوره يترك أثر سيئا على دقة الحلول العددية، وللتغلب على هذه المشكلة وتقليل الخطأ الناتج عن تطبيق صيغة تكامل ذات رتبة منخفضة أي قيمة n صغيرة فإننا نقسم فترة التكامل [a,b] إلى فترات أصغر ونطبق صيغة التكامل على كل فترة جزئية على حدة . والصيغة الناتجة من التطبيق المتكرر لصيغة ذات رتبة منخفضة تسمى صيغة تكامل مركبة . وفيما يلي نحصل على الصيغة المركبة الناتجة من التطبيق المتكرر لصيغة سمبسون البسيطة عندما (n=2) . ولكن من شروط قاعدة سمبسون ان تكون n عدد زوجي
وكي نكرر تطبيق قاعدة سمبسون عدد r من المرات نحتاج معرفة قيم الدالة عند عدد 2r+1 من نقاط الأساس 
كتاب الطرق العددية والتحليل العددي - أ.د. أبو بكر أحمد السيد، جامعة الكويت ISBN9789957171353

وتكون قاعدة سمبسون المركبة كالتالي :

 مثال على قاعدة سمبسون المركبة 

الحل بستخدام قاعدة سمبسون المركبة :

 الخطأ 

وفي بعض التكاملات على الفترات الكبيرة كان الحل بصيغة سمبسون يحتوي نسبة خطأ أكبر من لو قام بتقسيم الفترة

نلاحظ أن الخطأ يقل بعد تقسيم الفترة ومن هذا المنطلق نتجت صيغة سمبسون المركبة ولها العديد من الصيغ

والخطأ مرتبط بالتقريب :

حسب نظرية القيمة المتوسطة يوجد بحيث:

بالتالي

وبما أن

و أيضاً

 انظر ايضًا 
 قاعدة شبه المنحرف
 روجر كوتس
 إسحاق نيوتن
 التكامل العددي
 التحليل العددي

 مراجع 

 كتاب الطرق العددية والتحليل العددي - أ.د. أبو بكر أحمد السيد، جامعة الكويت ISBN9789957171353
 “Numerical Analysis”, Richard L. Burden & J. Douglas Faires, Brooks/Cole, Cengage Learning,(2011) ISBN 0-534-38216-9

 ملاحظات 

تصنيف:تحليل عددي
تصنيف:تكامل عددي
تصنيف:حساب تكاملي