طاولة بادي

هنري بادي

في التحليل المركب ، جدول باديه هو مصفوفة، قد تكون ذات امتداد لانهائي، من تقريبات باديه النسبية

R m , n

إلى متسلسلة قوى رسمية معقدة معينة . يمكن في كثير من الأحيان إثبات أن بعض متواليات التقريبات التي تقع ضمن جدول باديه تتوافق مع التقاربات المتتالية لتمثيل الكسر المستمر لدالة هولومورفية أو ميرومورفية .

تاريخ

على الرغم من أن علماء الرياضيات السابقين قد توصلوا إلى نتائج متفرقة تتضمن متواليات من التقريبات النسبية للدوال المتسامية ، إلا أن فروبينيوس (عام 1881) كان على ما يبدو أول من نظم هذه التقريبات في شكل جدول. وقد وسّع هنري بادي هذا المفهوم في أطروحته للدكتوراه بعنوان "حول التمثيل التقريبي لدالة بواسطة الكسور النسبية" عام 1892. وعلى مدى السنوات الست عشرة اللاحقة، نشر بادي 28 بحثًا إضافيًا استكشف فيها خصائص جدوله، وربط الجدول بالكسور المستمرة التحليلية. [ 1 ]

أُعيد إحياء الاهتمام الحديث بجداول باديه على يد إتش إس وول وأوسكار بيرون ، اللذين انصبّ اهتمامهما بالدرجة الأولى على الروابط بين هذه الجداول وأنواع معينة من الكسور المستمرة. ونشر دانيال شانكس وبيتر وين أبحاثًا مؤثرة حوالي عام 1955، وحقق دبليو بي غراغ نتائج تقارب واسعة النطاق خلال سبعينيات القرن العشرين. وفي الآونة الأخيرة، حفّز الاستخدام الواسع النطاق للحواسيب الإلكترونية اهتمامًا إضافيًا كبيرًا بهذا الموضوع. [ 2 ]

الترميز

يتم تمثيل الدالة f ( z ) بواسطة متسلسلة قوى رسمية :

و(z)=ج0+ج1z+ج2z2+=ل=0جلzل،{\displaystyle f(z)=c_{0}+c_{1}z+c_{2}z^{2}+\cdots =\sum _{l=0}^{\infty }c_{l}z^{l},}

حيث c ≠ 0 ، حسب الاصطلاح. يُعطى المدخل ( m , n ) [ 3 ] R m, n في جدول باديه للدالة f ( z ) بالصيغة التالية:

Rم،ن(z)=Pم(z)سؤالن(z)=أ0+أ1z+أ2z2++أمzمب0+ب1z+ب2z2++بنzن{\displaystyle R_{m,n}(z)={\frac {P_{m}(z)}{Q_{n}(z)}}={\frac {a_{0}+a_{1}z+a_{2}z^{2}+\cdots +a_{m}z^{m}}{b_{0}+b_{1}z+b_{2}z^{2}+\cdots +b_{n}z^{n}}}}

حيث أن P <sub>m</sub> ( z ) و Q <sub>n</sub> ( z ) كثيرتا حدود من الدرجة m و n على التوالي. ويمكن دائمًا إيجاد المعاملات { a <sub>i</sub> } و { b<sub> i</sub> } من خلال النظر في التعبير

و(z)ل=0م+نجلzل=:وأصصرox(z){\displaystyle f(z)\approx \sum _{l=0}^{m+n}c_{l}z^{l}=:f_{\mathrm {approx} }(z)}
سؤالن(z)وأصصرox(z)=Pم(z){\displaystyle Q_{n}(z)f_{\mathrm {approx} }(z)=P_{m}(z)}
سؤالن(z)(ج0+ج1z+ج2z2++جم+نzم+ن)=Pم(z){\displaystyle Q_{n}(z)\left(c_{0}+c_{1}z+c_{2}z^{2}+\cdots +c_{m+n}z^{m+n}\right)=P_{m}(z)}

and equating coefficients of like powers of z up through m + n. For the coefficients of powers m + 1 to m + n, the right hand side is 0 and the resulting system of linear equations contains a homogeneous system of n equations in the n + 1 unknowns bi, and so admits of infinitely many solutions each of which determines a possible Qn. Pm is then easily found by equating the first m coefficients of the equation above. However, it can be shown that, due to cancellation, the generated rational functions Rm, n are all the same, so that the (m, n)th entry in the Padé table is unique.[2] Alternatively, we may require that b0 = 1, thus putting the table in a standard form.

Although the entries in the Padé table can always be generated by solving this system of equations, that approach is computationally expensive. Usage of the Padé table has been extended to meromorphic functions by newer, timesaving methods such as the epsilon algorithm.[4]

The block theorem and normal approximants

Because of the way the (m, n)th approximant is constructed, the difference

Qn(z)f(z)  Pm(z)

is a power series whose first term is of degree no less than

m + n + 1.

If the first term of that difference is of degree

m + n + r + 1, r > 0,

then the rational functionRm, n occupies

(r + 1)2

cells in the Padé table, from position (m, n) through position (m+r, n+r), inclusive. In other words, if the same rational function appears more than once in the table, that rational function occupies a square block of cells within the table. This result is known as the block theorem.

إذا ظهرت دالة كسرية معينة مرة واحدة فقط في جدول باديه، تُسمى تقريبًا طبيعيًا للدالة f ( z ). إذا كان كل عنصر في جدول باديه الكامل طبيعيًا، يُقال إن الجدول نفسه طبيعي. يمكن توصيف تقريبات باديه الطبيعية باستخدام محددات المعاملات c <sub>n</sub> في متسلسلة تايلور للدالة f ( z )، كما يلي: يُعرَّف المحدد ( m, n) كما يلي : 

دم،ن=|جمجم-1...جم-ن+2جم-ن+1جم+1جم...جم-ن+3جم-ن+2جم+ن-2جم+ن-3...جمجم-1جم+ن-1جم+ن-2...جم+1جم|{\displaystyle D_{m,n}=\left|{\begin{matrix}c_{m}&c_{m-1}&\ldots &c_{m-n+2}&c_{m-n+1}\\c_{m+1}&c_{m}&\ldots &c_{m-n+3}&c_{m-n+2}\\\vdots &\vdots &&\vdots &\vdots \\c_{m+n-2}&c_{m+n-3}&\ldots &c_{m}&c_{m-1}\\c_{m+n-1}&c_{m+n-2}&\ldots &c_{m+1}&c_{m}\\\end{matrix}}\right|}

مع D m ,0 = 1، و D m ,1 = c m ، و c k  =  0 لـ k  <  0. إذن

  • يكون التقريب ( m , n ) للدالة f ( z ) طبيعيًا إذا وفقط إذا لم تكن أي من المحددات الأربعة D<sub> m , n - 1</sub> ، D <sub> m,n </sub> ، D <sub> m +1, n </sub> ، و D <sub> m +1, n +1 </sub> معدومة؛
  • يكون جدول باديه طبيعياً إذا وفقط إذا لم يكن أي من المحددات D m,n يساوي صفرًا (لاحظ على وجه الخصوص أن هذا يعني أنه لا يمكن أن يكون أي من المعاملات c k في تمثيل المتسلسلة لـ f ( z ) صفرًا). [ 5 ]

العلاقة بالكسور المستمرة

من أهم الأشكال التي يمكن أن يظهر بها الكسر المستمر التحليلي هو الكسر المستمر العادي ، وهو كسر مستمر على الصورة التالية:

و(z)=ب0+أ1z1-أ2z1-أ3z1-أ4z1-.$

حيث أن a i 0 هي ثوابت مركبة، و z هو متغير مركب.

ثمة ارتباط وثيق بين الكسور المستمرة المنتظمة وجداول باديه ذات التقريبات الطبيعية على طول القطر الرئيسي : فتسلسل "الدرج" لتقريبات باديه R₀ , ، R₁ , ₀ ، R₁ , ، R₂ , ، ... يكون طبيعيًا إذا وفقط إذا تطابق هذا التسلسل مع التقاربات المتتالية لكسر مستمر منتظم. بعبارة أخرى، إذا كان جدول باديه طبيعيًا على طول القطر الرئيسي، فيمكن استخدامه لإنشاء كسر مستمر منتظم، وإذا وُجد تمثيل كسر مستمر منتظم للدالة f ( z )، فإن القطر الرئيسي لجدول باديه الذي يمثل f ( z ) يكون طبيعيًا. [ 2 ]

مثال على ذلك – الدالة الأسية

فيما يلي مثال على جدول باديه، للدالة الأسية .

جزء من جدول باديه للدالة الأسية e z
ن
م
01234
011{\displaystyle {\frac {1}{1}}}11-z{\displaystyle {\frac {1}{1-z}}}11-z+12z2{\displaystyle {\frac {1}{1-z+{\scriptstyle {\frac {1}{2}}}z^{2}}}}11-z+12z2-16z3{\displaystyle {\frac {1}{1-z+{\scriptstyle {\frac {1}{2}}}z^{2}-{\scriptstyle {\frac {1}{6}}}z^{3}}}}11-z+12z2-16z3+124z4{\displaystyle {\frac {1}{1-z+{\frac {1}{2}}z^{2}-{\frac {1}{6}}z^{3}+{\frac {1}{24}}z^{4}}}}
11+z1{\displaystyle {\frac {1+z}{1}}}1+12z1-12z{\displaystyle {\frac {1+{\scriptstyle {\frac {1}{2}}}z}{1-{\scriptstyle {\frac {1}{2}}}z}}}1+13z1-23z+16z2{\displaystyle {\frac {1+{\scriptstyle {\frac {1}{3}}}z}{1-{\scriptstyle {\frac {2}{3}}}z+{\scriptstyle {\frac {1}{6}}}z^{2}}}}1+14z1-34z+14z2-124z3{\displaystyle {\frac {1+{\scriptstyle {\frac {1}{4}}}z}{1-{\scriptstyle {\frac {3}{4}}}z+{\scriptstyle {\frac {1}{4}}}z^{2}-{\scriptstyle {\frac {1}{24}}}z^{3}}}}1+15z1-45z+310z2-115z3+1120z4{\displaystyle {\frac {1+{\frac {1}{5}}z}{1-{\frac {4}{5}}z+{\frac {3}{10}}z^{2}-{\frac {1}{15}}z^{3}+{\frac {1}{120}}z^{4}}}}
21+z+12z21{\displaystyle {\frac {1+z+{\scriptstyle {\frac {1}{2}}}z^{2}}{1}}}1+23z+16z21-13z{\displaystyle {\frac {1+{\scriptstyle {\frac {2}{3}}}z+{\scriptstyle {\frac {1}{6}}}z^{2}}{1-{\scriptstyle {\frac {1}{3}}}z}}}1+12z+112z21-12z+112z2{\displaystyle {\frac {1+{\scriptstyle {\frac {1}{2}}}z+{\scriptstyle {\frac {1}{12}}}z^{2}}{1-{\scriptstyle {\frac {1}{2}}}z+{\scriptstyle {\frac {1}{12}}}z^{2}}}}1+25z+120z21-35z+320z2-160z3{\displaystyle {\frac {1+{\scriptstyle {\frac {2}{5}}}z+{\scriptstyle {\frac {1}{20}}}z^{2}}{1-{\scriptstyle {\frac {3}{5}}}z+{\scriptstyle {\frac {3}{20}}}z^{2}-{\scriptstyle {\frac {1}{60}}}z^{3}}}}1+13z+130z21-23z+15z2-130z3+1360z4{\displaystyle {\frac {1+{\frac {1}{3}}z+{\frac {1}{30}}z^{2}}{1-{\frac {2}{3}}z+{\frac {1}{5}}z^{2}-{\frac {1}{30}}z^{3}+{\frac {1}{360}}z^{4}}}}
31+z+12z2+16z31{\displaystyle {\frac {1+z+{\scriptstyle {\frac {1}{2}}}z^{2}+{\scriptstyle {\frac {1}{6}}}z^{3}}{1}}}1+34z+14z2+124z31-14z{\displaystyle {\frac {1+{\scriptstyle {\frac {3}{4}}}z+{\scriptstyle {\frac {1}{4}}}z^{2}+{\scriptstyle {\frac {1}{24}}}z^{3}}{1-{\scriptstyle {\frac {1}{4}}}z}}}1+35z+320z2+160z31-25z+120z2{\displaystyle {\frac {1+{\scriptstyle {\frac {3}{5}}}z+{\scriptstyle {\frac {3}{20}}}z^{2}+{\scriptstyle {\frac {1}{60}}}z^{3}}{1-{\scriptstyle {\frac {2}{5}}}z+{\scriptstyle {\frac {1}{20}}}z^{2}}}}1+12z+110z2+1120z31-12z+110z2-1120z3{\displaystyle {\frac {1+{\scriptstyle {\frac {1}{2}}}z+{\scriptstyle {\frac {1}{10}}}z^{2}+{\scriptstyle {\frac {1}{120}}}z^{3}}{1-{\scriptstyle {\frac {1}{2}}}z+{\scriptstyle {\frac {1}{10}}}z^{2}-{\scriptstyle {\frac {1}{120}}}z^{3}}}}1+37z+114z2+1210z31-47z+17z2-2105z3+1840z4{\displaystyle {\frac {1+{\frac {3}{7}}z+{\frac {1}{14}}z^{2}+{\frac {1}{210}}z^{3}}{1-{\frac {4}{7}}z+{\frac {1}{7}}z^{2}-{\frac {2}{105}}z^{3}+{\frac {1}{840}}z^{4}}}}
41+z+12z2+16z3+124z41{\displaystyle {\frac {1+z+{\scriptstyle {\frac {1}{2}}}z^{2}+{\scriptstyle {\frac {1}{6}}}z^{3}+{\scriptstyle {\frac {1}{24}}}z^{4}}{1}}}1+45z+310z2+115z3+1120z41-15z{\displaystyle {\frac {1+{\scriptstyle {\frac {4}{5}}}z+{\scriptstyle {\frac {3}{10}}}z^{2}+{\scriptstyle {\frac {1}{15}}}z^{3}+{\scriptstyle {\frac {1}{120}}}z^{4}}{1-{\scriptstyle {\frac {1}{5}}}z}}}1+23z+15z2+130z3+1360z41-13z+130z2{\displaystyle {\frac {1+{\scriptstyle {\frac {2}{3}}}z+{\scriptstyle {\frac {1}{5}}}z^{2}+{\scriptstyle {\frac {1}{30}}}z^{3}+{\scriptstyle {\frac {1}{360}}}z^{4}}{1-{\scriptstyle {\frac {1}{3}}}z+{\scriptstyle {\frac {1}{30}}}z^{2}}}}1+47z+17z2+2105z3+1840z41-37z+114z2-1210z3{\displaystyle {\frac {1+{\scriptstyle {\frac {4}{7}}}z+{\scriptstyle {\frac {1}{7}}}z^{2}+{\scriptstyle {\frac {2}{105}}}z^{3}+{\scriptstyle {\frac {1}{840}}}z^{4}}{1-{\scriptstyle {\frac {3}{7}}}z+{\scriptstyle {\frac {1}{14}}}z^{2}-{\scriptstyle {\frac {1}{210}}}z^{3}}}}1+12z+328z2+184z3+11680z41-12z+328z2-184z3+11680z4{\displaystyle {\frac {1+{\frac {1}{2}}z+{\frac {3}{28}}z^{2}+{\frac {1}{84}}z^{3}+{\frac {1}{1680}}z^{4}}{1-{\frac {1}{2}}z+{\frac {3}{28}}z^{2}-{\frac {1}{84}}z^{3}+{\frac {1}{1680}}z^{4}}}}

تظهر عدة سمات على الفور.

  • يتكون العمود الأول من الجدول من عمليات القطع المتتالية لسلسلة تايلور لـ e z .
  • وبالمثل، يحتوي الصف الأول على مقلوبات عمليات القطع المتتالية لتوسيع السلسلة لـ e z .
  • التقريبان R<sub> m,n</sub> و R <sub>n,m </sub> متناظران تمامًا - حيث يتبادل البسط والمقام، وتختلف أنماط الإشارات الموجبة والسالبة، لكن المعاملات نفسها تظهر في كليهما. ويمكن التعبير عنهما بدلالة دوال خاصة كما يلي:
Rم،ن=1F1(-م؛-م-ن؛z)1F1(-ن؛-م-ن؛-z)=ن!2مθم(z2؛ن-م+2،2)م!2نθن(-z2؛م-ن+2،2){\displaystyle R_{m,n}={\frac {{}_{1}F_{1}(-m;-m-n;z)}{{}_{1}F_{1}(-n;-m-n;-z)}}={\frac {n!\,2^{m}\theta _{m}\left({\frac {z}{2}};n-m+2,2\right)}{m!\,2^{n}\theta _{n}\left(-{\frac {z}{2}};m-n+2,2\right)}}}،
أين1F1(أ؛ب؛z){\displaystyle {}_{1}F_{1}(a;b;z)}هي متسلسلة هندسية فائقة معممة وθن(x؛α،β){\displaystyle \theta _{n}(x;\alpha ,\beta )}هي متعددة حدود بيسل معكوسة معممة . [ 6 ]
تُختزل التعبيرات الموجودة على القطر الرئيسي إلىRن،ن=θن(z/2)/θن(-z/2){\displaystyle R_{n,n}=\theta _{n}(z/2)/\theta _{n}(-z/2)}، أينθن(x){\displaystyle \theta _{n}(x)}هي متعددة حدود بيسل معكوسة . [ 7 ]
  • يمكن إجراء العمليات الحسابية التي تتضمن R n , n (على القطر الرئيسي) بكفاءة عالية. على سبيل المثال، يُعيد R 3,3 إنتاج متسلسلة القوى للدالة الأسية بدقة تامة حتى 1 / 720 z 6 ، ولكن نظرًا لتناظر كثيرتي الحدود التكعيبيتين، يمكن ابتكار خوارزمية تقييم سريعة جدًا.

يمكن تطبيق الإجراء المستخدم لاشتقاق الكسر المستمر لغوس على سلسلة هندسية فائقة متقاربة معينة لاشتقاق توسيع الكسر C التالي للدالة الأسية، وهو صالح في جميع أنحاء المستوى المركب بأكمله:

هـz=1+z1-12z1+16z1-16z1+110z1-110z1+-.{\displaystyle e^{z}=1+{\cfrac {z}{1-{\cfrac {{\frac {1}{2}}z}{1+{\cfrac {{\frac {1}{6}}z}{1-{\cfrac {{\frac {1}{6}}z}{1+{\cfrac {{\frac {1}{10}}z}{1-{\cfrac {{\frac {1}{10}}z}{1+-\ddots }}}}}}}}}}}}.}

بتطبيق صيغ التكرار الأساسية، يمكن التحقق بسهولة من أن التقاربات المتتالية لهذا الكسر C هي متتالية متدرجة من تقريبات باديه R 0,0 ، R 1,0 ، R 1,1 ، ... في هذه الحالة تحديدًا، يمكن الحصول على كسر مستمر وثيق الصلة من المتطابقة

هـz=1هـ-z؛{\displaystyle e^{z}={\frac {1}{e^{-z}}};}

يبدو هذا الكسر المستمر على النحو التالي:

هـz=11-z1+12z1-16z1+16z1-110z1+110z1-+.{\displaystyle e^{z}={\cfrac {1}{1-{\cfrac {z}{1+{\cfrac {{\frac {1}{2}}z}{1-{\cfrac {{\frac {1}{6}}z}{1+{\cfrac {{\frac {1}{6}}z}{1-{\cfrac {{\frac {1}{10}}z}{1+{\cfrac {{\frac {1}{10}}z}{1-+\ddots }}}}}}}}}}}}}}.}

تظهر التقاربات المتتالية لهذا الكسر أيضًا في جدول باديه، وتشكل المتتالية R 0,0 ، R 0,1 ، R 1,1 ، R 1,2 ، R 2,2 ، ...

التعميمات

تكون متسلسلة نيوتن الرسمية L على الشكل التالي

L(z)=ج0+ن=1جنك=1ن(z-βك){\displaystyle L(z)=c_{0}+\sum _{n=1}^{\infty }c_{n}\prod _{k=1}^{n}(z-\beta _{k})}

حيث تُعرف متتالية النقاط { βk } في المستوى المركب بمجموعة نقاط الاستيفاء . ويمكن تكوين متتالية من التقريبات النسبية Rm ,n لمثل هذه المتسلسلة L بطريقة مماثلة تمامًا للإجراء الموصوف أعلاه، ويمكن ترتيب هذه التقريبات في جدول نيوتن-باديه . وقد بُيّن [ 8 ] أن بعض متتاليات "الدرج" في جدول نيوتن-باديه تتوافق مع التقاربات المتتالية لكسر مستمر من نوع ثيل، والذي يكون على الشكل التالي:

أ0+أ1(z-β1)1-أ2(z-β2)1-أ3(z-β3)1-.{\displaystyle a_{0}+{\cfrac {a_{1}(z-\beta _{1})}{1-{\cfrac {a_{2}(z-\beta _{2})}{1-{\cfrac {a_{3}(z-\beta _{3})}{1-\ddots }}}}}}.}

قام علماء الرياضيات أيضًا بإنشاء جداول باديه ذات النقطتين من خلال النظر في سلسلتين، إحداهما في قوى z ، والأخرى في قوى 1/ z ، والتي تمثل بالتناوب الدالة f ( z ) في جوار الصفر وفي جوار اللانهاية. [ 2 ]

انظر أيضاً

ملحوظات

  1. أوكونور، جون جيه؛ روبرتسون، إدموند إف ، "جدول باديه" ، أرشيف ماك تيوتور لتاريخ الرياضيات ، جامعة سانت أندروز
  2. 1 2 3 4 جونز وثورن، 1980.
  3. يعتبر المدخل( m ، n ) موجودًا في الصف m والعمود n ، ويبدأ ترقيم الصفوف والأعمدة من (0، 0).
  4. وين، بيتر (أبريل 1956). "حول جهاز لحساب تحويل e^ m ( Sn ) ". الجداول الرياضية وغيرها من الوسائل المساعدة في الحساب . 10 (54). الجمعية الرياضية الأمريكية: 91-96 . doi : 10.2307/2002183 . JSTOR 2002183 . 
  5. غراغ، دبليو بي (يناير 1972). "جدول باديه وعلاقته ببعض خوارزميات التحليل العددي" . مجلة SIAM . 14 (1): 1-62 . doi : 10.1137/1014001 . ISSN 0036-1445 . JSTOR 2028911 .  
    • أندرهيل، سي. (1986). "بعض الخصائص التقاربية لتقريبات باديه لـهـ-x{\displaystyle e^{-x}}". رياضيات الحوسبة . 47 (175): 253– 263. JSTOR 2008092 . 
  6. ^ ثييل، تينيسي (1909). الاستيفاءات rechnung . لايبزيغ: تيوبنر. رقم ISBN 1-4297-0249-4.{{cite book}}عدم توافق رقم ISBN / التاريخ ( مساعدة )

مراجع