صيغة أويلر-ماكلورين

من ويكيبيديا، الموسوعة الحرة
اذهب إلى: تصفح، ‏ ابحث

في الرياضيات, تعطي صيغة أويلر-ماكلورين ارتباطا وثيقا بين التكامل (انظر التفاضل والتكامل) والمجموع.

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

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

الصيغة[عدل]

إذا كان n عددا طبيعيا, وكانت f(x)\, دالة أكثر سلاسة(بمعنى أنها قابلة للاشتقاق بما يكفي) معرفة لجميع الأعداد الحقيقية x بين 0 وn, فإن التكامل

I=\int_0^n f(x)\,dx

يمكن تقريبه بالمجموع (والعكس صحيح)


S=\frac{1}{2}f(0)+f\left(1\right) +\cdots+f\left(n-1\right) +\frac{1}{2}f(n)

تعطي صيغة أويلر-ماكلورين تعابير للفرق بين المجموع والتكامل بدلالة المشتقات العليا f^{(k)}\, عند اطراف النقاط على الفترة 0 وn. وبوضوح, لأي عدد طبيعي p, يكون لدينا (انظر قاعدة شبه المنحرف).

 S-I = \sum_{k=2}^p {\frac{B_k}{k!}\left(f^{(k-1)}(n)-f^{(k-1)}(0)\right)} +R

حيث B1 = −1/2، B2 = 1/6، B3 = 0، B4 = −1/30, B5 = 0، B6 = 1/42، B7 = 0، B8 = −1/30,... هي أعداد بيرنولي, وR هو حد الخطأ وعادة مايكون صغيرا عند قيم مناسبة لـp. (غالبا يتم كتابة الصيغة بالقيم الزوجية, بما أن أعداد بيرنولي الفردية أصفار عدا B1.)

لاحظ أن

-B_1(f(n)+f(0)) =\frac{1}{2}(f(n)+f(0)).

وعليه يمكن أيضا كتابة الصيغة كما يلي:


\begin{align}
& \quad \sum_{i=0}^n f(i) = \int^n_0f(x)\,dx-B_1(f(n)+f(0))+\sum_{k=2}^p\frac{B_k}{k!}\left(f^{(k-1)}(n)-f^{(k-1)}(0)\right)+R.
\end{align}

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

الحد المتبقي[عدل]

الحد المتبقي R يعبر عنه غالبا باستخدام كثيرات حدود بيرنولي الدورية (Pn(x. تعرف كثيرة حدود بيرنولي Bn(xn = 0, 1, 2, ... بشكل تعاودي على أنها

B_0(x) = 1, \,
 B_n'(x) = nB_{n-1}(x)\mbox{ and }\int_0^1 B_n(x)\,dx = 0\mbox{ for }n \ge 1.

حينئذ فإن دوال بيرنولي الدورية Pn تعرف بالشكل

 P_n(x) = B_n(x - \lfloor x\rfloor)\mbox{ for }x> 0, \,

حيث ترمز\scriptstyle \lfloor x\rfloor إلى العدد الصحيح الأكبر والذي لا يكون أكبر من x. وعليه, بدلالة Pn(x), فإن الحد المتبقي R يمكن كتابته بالصورة

 R = (-1)^{p+1} \int_0^n f^{(p)}(x) {P_p(x) \over p!}\,dx,

بما أن |P_p(x)|\, عادة أقل من 2p!/(2\pi)^p\,, فيمكن توقع حجم الحد المتبقي باستخدام

\left|R\right|\leq\frac{2}{(2\pi)^p}\int_0^n\left|f^{(p)}(x)\right|\,dx.

تطبيقات[عدل]

مسألة بازل[عدل]

كانت مسألة بازل تتسأل عن إيجاد المجموع

 1 + \frac14 + \frac19 + \frac1{16} + \frac1{25} + \cdots = \sum_{n=1}^\infty \frac{1}{n^2}.

قام اويلر بحساب هذا المجموع بدقة عشرين خانة عشرية باستعمال حدود معدودة من صيغة أويلر ماكلورين سنة 1735. وربما أقنعه هذا بأن المجوموع مكافئا ل π2 / 6, والتي أثبتها في نفس العام.[1]

مجاميع ذات كثيرة حدود[عدل]

إذا كانت f دالة كثيرة حدود وp كبيرة بشكل كاف, فإن الحد المتبقي يتلاشى. على سبيل المثال، إذا كانت f(x) = x3,, فيمكننا اختيار p = 2 لتوضيح الاتي بعد التبسيط

\sum_{i=0}^n i^3=\left(\frac{n(n+1)}{2}\right)^2

(انظر صيغة فاولابر).

التكامل العددي[عدل]

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

نشر المحاميع المقارب[عدل]

في سياق حساب المتسلسلات المقاربة للمجاميع والمتسلسلات, فإن الشكل التالي يكون عالبا أفضل صيغة في أويلر ماكلورين

\sum_{n=a}^b f(n) \sim \int_a^b f(x)\,dx + \frac{f(a)+f(b)}{2} + \sum_{k=1}^\infty \,\frac{B_{2k}}{(2k)!}\left(f^{(2k-1)}(b)-f^{(2k-1)}(a)\right), \,

حيثa وb أعداد صحيحة.[2] وفي الغالب يبقى النشر مشروعا حتى بعد أحذ النهايات {\scriptstyle a\to -\infty}\, أو {\scriptstyle b\to +\infty}\,, أو كلاهما. في حالات عديدة يمكن تقدير التكامل في الجانب الأيمن في شكل مغلق بدلالة الدوال الأساسية حتى ولو أن الجمع على الجانب الأيسر غير قادر على ذلك. حينئذ يمكن التعبير عن جميع الحدود في السلسلة المقاربة بدلالة دوال أساسية.

\sum_{k=0}^\infty \frac{1}{(z+k)^2} \sim \underbrace{\int_0^\infty\frac{1}{(z+k)^2}\,dk}_{=1/z}+\frac{1}{2z^2}
+\sum_{t=1}^\infty \frac{B_{2t}}{z^{2t+1}}.\,

هنا الجانب الأيسر يساوي {\scriptstyle \psi^{(1)}(z)}\,, بالاسم دالة متعدد غاما من الرتبة الأولى المعرفة من خلال {\scriptstyle \psi^{(1)}(z)=\frac{d^2}{dz^2}\ln \Gamma(z)}\,; دالة غاما {\scriptstyle \Gamma(z)}\, تساوي {\scriptstyle (z-1)!}\, إذا كان {\scriptstyle z}\, عدد صحيح. ينتج عن هذا نشر مقارب لـ {\scriptstyle \psi^{(1)}(z)}\,. وهذا النشر بدوره, يخدمنا كنقطة بداية لواحدة من اشتقاقات دقيقة لتوقع الخطأ في تقريب ستيرلينغ لدالة المضروب.

البراهين[عدل]

الاشتقاق بالاستقراء الرياضي[عدل]

باتباع النقاش المعطى في (Apostol) [3].

يمكن تعريف كثيرات حدود بيرنولي Bn(x), n = 0, 1, 2, ... بشكل معاود كما يلي:

B_0(x) = 1, \,
 B_n'(x) = nB_{n-1}(x)\mbox{ and }\int_0^1 B_n(x)\,dx = 0\mbox{ for }n \ge 1.

البعض الأولى من هذه هي

 B_1(x)=x-1/2, \quad B_2(x)=x^2-x+1/6,
 B_3(x) = x^3-\frac{3}{2}x^2+\frac{1}{2}x, \quad B_4(x)=x^4-2x^3+x^2-\frac{1}{30}, \dots

القيم Bn(1) هي أعداد بيرنولي. لاحظ أن لأجل n ≥ 2 يكون لدينا

B_n(0) = B_n(1) = B_n\quad(:n\text{th Bernoulli number}).

نعرف دوال بيرنولي الدورية Pn بالصورة

 P_n(x) = B_n(x - \lfloor x\rfloor)\mbox{ for }0 <x <1, \,

حيث  \lfloor x\rfloor ترمز لأكبر عدد صحيح ليس أكبر من x. وعليه فإن Pn تتفق مع كثيرات حدود بيرنولي على الفترة (0, 1) كما أنها دورية إذا كان الدور = 1. إذن,

 P_n(0) = P_n(1)= B_n\quad \text{for } n>1.

لأجلn = 1,

 P_1(0) = - P_1(1)=B_1.

والان, باعتبار التكامل

 \int_k^{k+1} f(x)\,dx = \int u\,dv,

حيث

\begin{align}
u &{}= f(x), \\
du &{}= f'(x)\,dx, \\
v &{}= P_1(x),\\
dv &{}= P_0(x)\,dx \quad (\mbox{since }P_0(x)=1). \\
\end{align}

التكامل بالتجزيء, نحصل على

\begin{align}
\int_k^{k+1} f(x)\,dx &= uv - \int v\,du &{}\\
&= \Big[f(x)P_1(x) \Big]_k^{k+1} - \int_k^{k+1} f'(x)P_1(x)\,dx \\[8pt]
&=-B_1(f(k) + f(k+1)) - \int_k^{k+1} f'(x)P_1(x)\,dx.
\end{align}

وبجمع الصورة السابقة k = 0 to k = n − 1, نحصل على


\begin{align}
&\int_0^1 f(x)\,dx+\dotsb+\int_{n-1}^n f(x)\,dx \\
&= \int_0^n f(x)\, dx  \\
&= \frac{f(0)}{2}+ f(1) + \dotsb + f(n-1) + {f(n) \over 2} - \int_1^n f'(x) P_1(x)\,dx. 
\end{align}

بإضافة (ƒ(0) + ƒ(n))/2 إلى كلا الطرفين واعادة الترتيب نحصل على

 \sum_{k=0}^n f(k) = \int_0^n f(x)\,dx + {f(0) + f(n) \over 2} + \int_0^n f'(x) P_1(x)\,dx.\qquad (1)

يعطي الحدين الأخيرين بالتالي الخطأ عند تقريب التكامل بالمجموع.

فيما يلي لنعتبر

 \int_k^{k+1} f'(x)P_1(x)\,dx = \int u\,dv,

حيث

\begin{align}
u &{}= f'(x), \\
du &{}= f''(x)\,dx, \\
v &{}= P_2(x)/2\\
dv &{}= P_1(x)\,dx.
\end{align}

وبالتكامل بالتجزيء مرة أخرى, نحصل على

\begin{align}
uv - \int v\,du &{}= \left[ {f'(x)P_2(x) \over 2} \right]_k^{k+1} - {1 \over 2}\int_k^{k+1} f''(x)P_2(x)\,dx \\  \\
&{}= {f'(k+1) - f'(k) \over 12} -{1 \over 2}\int_k^{k+1} f''(x)P_2(x)\,dx.
\end{align}

الجمع من k = 0 إلى k = n − 1, ومن ثم بتعويض التكامل الأخير في (1) بالذي أثبتنا للتو أنه مساو له, نجد أنe

 \sum_{k=0}^n f(k) = \int_0^n f(x)\,dx + {f(0) + f(n) \over 2} + \frac{B_2}{2}(f'(n) - f'(0)) - {1 \over 2}\int_0^n f''(x)P_2(x)\,dx.

والآن, لابد أن القارئ قد حزر بأن هذه العملية تكرارية. استطعنا الحصول على إثبات لصيغة مجموع أويلر-ماكلورين بالاستقراء الرياضي, حيث اعتمدت خطوة الاستقراء على التكامل بالتجزيء ومطابقات دوال بيرنولي الدورية.

الاشتقاق بالتحليل الدالي[عدل]

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

في البداية يتم حصر المسألة في مجال فترة الوحدة [0,1]. لتكن B_n(x)\, هي كثيرات حدود بيرنولي. تعطى زمرة دول ثنائية لكثيرات حدود بيرنولي بالعلاقة:

\tilde{B}_n(x)=\frac{(-1)^{n+1}}{n!} \left[ 
\delta^{(n-1)}(x-1) - \delta^{(n-1)}(x) \right]

حيث δ هي دالة دايراك دلتا. العلاقة السابقة هي رمز لفكرة أخذ المشتقات عند نقطة; بالتالي يكون لدينا

\int_0^1 \tilde{B}_n(x) f(x)\, dx = \frac{1}{n!} \left[ 
f^{(n-1)}(1) - f^{(n-1)}(0) \right]

من أجل n > 0 ودالة اعتباطية ولكنها تفاضلية ƒ(x) على دالة الوحدة. في الحالة n = 0, يمكن تعريف \tilde{B}_0(x)=1. كثيرات حدود بيرنولي, على امتداد ثنائياتها, من مجموعة حالات متعامدة على فترة الوحدة: لدينا

\int_0^1 \tilde{B}_m(x) B_n(x)\, dx = \delta_{mn}

و

\sum_{n=0}^\infty B_n(x) \tilde{B}_n(y) = \delta (x-y).

حينئذ تتبع صيغة مجموع أويلر ماكلورين كتكامل على الآخر. لدينا

f(x)=\int_0^1 \sum_{n=0}^\infty B_n(x) \tilde{B}_n(y) f(y)\, dy
=\int_0^1 f(y)\,dy + 
\sum_{n=1}^N B_n(x) \frac{1}{n!} 
\left[ f^{(n-1)}(1) - f^{(n-1)}(0) \right] 
- \frac{1}{(N+1)!} \int_0^1 B_{N+1}(x-y) f^{(N)}(y)\, dy.

وبوضع x = 0 وترتيب الحدود, يمكن الحصول على تعبير لـ ƒ(0). لاحظ أن أعداد بيرنولي تعرف بأنها Bn = Bn(0), وأنها تتلاشى للقيم الفردية n أكبر من 1.

وبالتالي, باستخدام دالة بيرولي الدورية Pn المعرفة أعلاه وبإعادة النقاش حول الفترة [1,2], يمكننا الحصول على تعبير ƒ(1). وعلى هذا المنوال يمكن الحصول على تعبير لـ ƒ(n), n = 0, 1, 2, ..., N, وبإضافتها فوق بعض نحصل على صيغة أويلر ماكلورين. لاحظ أن هذا الاشتقاق لايفترض أن ƒ(x) قابلة للتفاضل بشكل كاف وذات سلوك صحيح ; خاصة أن ƒ يمكن تقريبها بكثيرات حدود; وبشكل مكافئ بأن ƒ دالة تحليلية حقيقية.

يمكن إذن النظر إلى صيغة مجموع أويلر-ماكلورين على أنها حصيلة تمثيل دوال على فترة الوحدة بالضرب المباشر لكثيرات حدود بيرنولي وثنائياتها. ومع ذلك لاحظ أن التمثيل ليس مكتملا على مجموعة دوال مربعة قابلة للتكامل. النشر بدلالة كثيرات حدود بيرنولي ليس له نواة عادية. بشكل خاص، sin(2πnx) تقع على النواة; تكامل sin(2πnx) يتلاشى على فترة الوحدة, كما هو الفرق بين مشتقاته على النقاط الطرفية.

ملاحظات[عدل]

  1. ^ David J. Pengelley, "Dances between continuous and discrete: Euler's summation formula", in: Robert Bradley and Ed Sandifer (Eds), Proceedings, Euler 2K+2 Conference (Rumford, Maine, 2002)، Euler Society, 2003.
  2. ^ Abramowitz & Stegun (1972), 23.1.30
  3. ^ Apostol، T. M. (1 May 1999). "An Elementary View of Euler's Summation Formula". The American Mathematical Monthly (Mathematical Association of America) 106 (5): 409–418. doi:10.2307/2589145. ISSN 00029890. JSTOR 2589145.  edit

المصادر[عدل]

انظر أيضا[عدل]