باديه التقريبي

هنري بادي

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

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

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

تعريف

بفرض دالة f وعددين صحيحين m ≥ 0 و n ≥ 1 ، فإن تقريب باديه من الرتبة [ m / n ] هو الدالة الكسرية

R(x)=ج=0مأجxج1+ك=1نبكxك=أ0+أ1x+أ2x2++أمxم1+ب1x+ب2x2++بنxن،{\displaystyle R(x)={\frac {\sum _{j=0}^{m}a_{j}x^{j}}{1+\sum _{k=1}^{n}b_{k}x^{k}}}={\frac {a_{0}+a_{1}x+a_{2}x^{2}+\dots +a_{m}x^{m}}{1+b_{1}x+b_{2}x^{2}+\dots +b_{n}x^{n}}},} وهو ما يتوافق مع f ( x ) إلى أعلى رتبة ممكنة، وهو ما يعادل و(0)=R(0)،و(0)=R(0)،و"(0)=R"(0)،و(م+ن)(0)=R(م+ن)(0).{\displaystyle {\begin{aligned}f(0)&=R(0),\\f'(0)&=R'(0),\\f''(0)&=R''(0),\\&\mathrel {\;\vdots } \\f^{(m+n)}(0)&=R^{(m+n)}(0).\end{aligned}}}

وبعبارة أخرى، إذاR(x){\displaystyle R(x)}يتم توسيعها في متسلسلة ماكلورين ( متسلسلة تايلور عند 0)، وهي الأولىم+ن{\displaystyle m+n}ستكون الحدود مساوية للأولىم+ن{\displaystyle m+n}شروطو(x){\displaystyle f(x)}وبالتالي و(x)-R(x)=جم+ن+1xم+ن+1+جم+ن+2xم+ن+2+{\displaystyle f(x)-R(x)=c_{m+n+1}x^{m+n+1}+c_{m+n+2}x^{m+n+2}+\cdots }

عندما يكون موجودًا، يكون تقريب باديه فريدًا كسلسلة قوى رسمية لـ m و n المعطاة . [ 1 ]

يُشار أيضًا إلى تقريب باديه المعرّف أعلاه باسم [م/ن]و(x).{\displaystyle [m/n]_{f}(x).}

حساب

بالنسبة لقيمة f المعطاة ، يمكن حساب تقريبات باديه بواسطة خوارزمية إبسيلون لوين [ 2 ] وأيضًا بواسطة تحويلات متسلسلة أخرى [ 3 ] من المجاميع الجزئية. تيشمال(x)=ج0+ج1x+ج2x2++جشمالxشمال{\displaystyle T_{N}(x)=c_{0}+c_{1}x+c_{2}x^{2}+\cdots +c_{N}x^{N}} من متسلسلة تايلور للدالة f ، أي لدينا جك=و(ك)(0)ك!.{\displaystyle c_{k}={\frac {f^{(k)}(0)}{k!}}.}يمكن أن تكون f أيضًا سلسلة قوى رسمية ، وبالتالي، يمكن أيضًا تطبيق تقريبات باديه على مجموع المتسلسلات المتباعدة .

إحدى طرق حساب تقريب باديه هي عبر خوارزمية إقليدس الموسعة للقاسم المشترك الأكبر لكثير الحدود . [ 4 ] العلاقة R(x)=P(x)/سؤال(x)=تيم+ن(x)تعديلxم+ن+1{\displaystyle R(x)=P(x)/Q(x)=T_{m+n}(x){\bmod {x}}^{m+n+1}} يكافئ وجود عامل ماك(x){\displaystyle K(x)}بحيث P(x)=سؤال(x)تيم+ن(x)+ك(x)xم+ن+1،{\displaystyle P(x)=Q(x)T_{m+n}(x)+K(x)x^{m+n+1},} والتي يمكن تفسيرها على أنها متطابقة بيزو لخطوة واحدة في حساب القاسم المشترك الأكبر الموسع لكثيرات الحدودتيم+ن(x){\displaystyle T_{m+n}(x)}وxم+ن+1{\displaystyle x^{m+n+1}}.

تذكر أنه لحساب القاسم المشترك الأكبر لكثيرتي حدود p و q ، يتم حساب متتابعة الباقي عن طريق القسمة المطولة ر0=ص،ر1=q،رك-1=qكرك+رك+1،{\displaystyle r_{0}=p,\;r_{1}=q,\quad r_{k-1}=q_{k}r_{k}+r_{k+1},}k = 1، 2، 3، ... معدرجةرك+1<درجةرك{\displaystyle \deg r_{k+1<\deg r_{k}\,}، حتىرك+1=0{\displaystyle r_{k+1}=0}بالنسبة لهويات بيزو للقاسم المشترك الأكبر الموسع، يتم حساب متتاليتي كثيرات الحدود في آن واحد. u0=1،v0=0،u1=0،v1=1،uك+1=uك-1-qكuك،vك+1=vك-1-qكvك{\displaystyle u_{0}=1,\;v_{0}=0,\quad u_{1}=0,\;v_{1}=1,\quad u_{k+1}=u_{k-1}-q_{k}u_{k},\;v_{k+1}=v_{k-1}-q_{k}v_{k}} للحصول على هوية بيزو في كل خطوة رك(x)=uك(x)ص(x)+vك(x)q(x).{\displaystyle r_{k}(x)=u_{k}(x)p(x)+v_{k}(x)q(x).}

بالنسبة للتقريب [ m / n ] ، يتم تنفيذ خوارزمية إقليدس الموسعة لـ ر0=xم+ن+1،ر1=تيم+ن(x){\displaystyle r_{0}=x^{m+n+1},\;r_{1}=T_{m+n}(x)} ويوقفه في اللحظة الأخيرة التيvك{\displaystyle v_{k}}له درجة n أو أقل.

ثم كثيرات الحدودP=رك،سؤال=vك{\displaystyle P=r_{k},\;Q=v_{k}}أعطِ تقريب باديه [ m / n ] . إذا تم حساب جميع خطوات حساب القاسم المشترك الأكبر الموسع، فسيتم الحصول على قطر مضاد لجدول باديه .

دالة زيتا لريمان-باديه

لدراسة إعادة تجميع سلسلة متباعدة ، على سبيل المثال z=1و(z)،{\displaystyle \sum _{z=1}^{\infty }f(z),} قد يكون من المفيد تقديم دالة باديه أو ببساطة دالة زيتا النسبية كـ ζR(s)=z=1R(z)zs،{\displaystyle \zeta _{R}(s)=\sum _{z=1}^{\infty }{\frac {R(z)}{z^{s}}},} أين R(x)=[م/ن]و(x){\displaystyle R(x)=[m/n]_{f}(x)} يمثل تقريب باديه من الرتبة ( m , n ) للدالة f ( x ) . وتُعتبر قيمة تنظيم زيتا عند s = 0 هي مجموع المتسلسلة المتباعدة.

المعادلة الوظيفية لدالة زيتا باديه هذه هي ج=0نأجζR(s-ج)=ج=0مبجζ0(s-ج)،{\displaystyle \sum _{j=0}^{n}a_{j}\zeta _{R}(sj)=\sum _{j=0}^{m}b_{j}\zeta _{0}(sj),} حيث يمثل aj و bj معاملات تقريب باديه. يشير الرمز السفلي '0' إلى أن تقريب باديه من الرتبة [0/0]، وبالتالي، لدينا دالة زيتا لريمان .

طريقة DLog Padé

يمكن استخدام تقريبات باديه لاستخراج النقاط الحرجة وأسس الدوال. [ 5 ] [ 6 ] في الديناميكا الحرارية، إذا كانت الدالة f ( x ) تتصرف بطريقة غير تحليلية بالقرب من نقطة x = مثلو(x)|x-ر|ص{\displaystyle f(x)\sim |xr|^{p}}يُطلق على النقطة x = r اسم النقطة الحرجة، وعلى p اسم الأس الحرج المرتبط بها للدالة f . إذا عُرفت حدود كافية من متسلسلة تايلور للدالة f ، يُمكن استخلاص النقاط الحرجة والأسس الحرجة تقريبًا من أقطاب وبواقي تقريبات باديه، على التوالي.[ن/ن+1]ز(x){\displaystyle [n/n+1]_{g}(x)}، أينز=و/و{\displaystyle g=f'/f}.

التعميمات

يُقارب مُقارب باديه دالة في متغير واحد. ويُسمى المُقارب في متغيرين مُقارب تشيشولم (نسبةً إلى جيه إس آر تشيشولم[ 7 ] وفي عدة متغيرات يُسمى مُقارب كانتربري (نسبةً إلى غريفز-موريس في جامعة كنت ). [ 8 ]

تقريب باديه بنقطتين

يُحدد أن تقريب باديه التقليدي يُعيد إنتاج متسلسلة ماكلورين حتى رتبة معينة. لذلك، قد يكون التقريب ضعيفًا عند القيم البعيدة عن نقطة التوسع. ويتم تجنب ذلك باستخدام تقريب باديه ثنائي النقاط، وهو نوع من طرق الجمع متعددة النقاط. [ 9 ] عندx=0{\displaystyle x=0}لنفترض حالة دالةو(x){\displaystyle f(x)}والذي يتم التعبير عنه من خلال السلوك التقاربيو0(x){\displaystyle f_{0}(x)}: وو0(x)+o(و0(x))،x0،{\displaystyle f\sim f_{0}(x)+o{\big (}f_{0}(x){\big )},\quad x\to 0,} وفيx{\displaystyle x\to \infty }سلوك تقاربي إضافيو(x){\displaystyle f_{\infty }(x)}: و(x)و(x)+o(و(x))،x.{\displaystyle f(x)\sim f_{\infty }(x)+o{\big (}f_{\infty }(x){\big )},\quad x\to \infty .}

عن طريق اختيار السلوك الرئيسي لـو0(x)،و(x){\displaystyle f_{0}(x),f_{\infty }(x)}الدوال التقريبيةF(x){\displaystyle F(x)}يمكن إيجاد طريقة لإعادة إنتاج السلوك التقاربي في آن واحد من خلال تطوير تقريب باديه في حالات مختلفة. ونتيجة لذلك، عند النقطةx{\displaystyle x\to \infty }في حين أن دقة التقريب قد تكون في أدنى مستوياتها في تقريب باديه العادي، فإن دقة تقريب باديه ثنائي النقاط مضمونة. لذلك، يمكن أن يكون تقريب باديه ثنائي النقاط طريقة تعطي تقريبًا جيدًا بشكل عام لـx=0{\displaystyle x=0\sim \infty }[ 10 ] [ 11 ] [ 12 ] [ 13 ]

في الحالات التيو0(x)،و(x){\displaystyle f_{0}(x),f_{\infty }(x)}تُعبَّر عنها بواسطة كثيرات الحدود أو متسلسلات القوى السالبة، أو الدالة الأسية، أو الدالة اللوغاريتمية، أوxlnx{\displaystyle x\ln x}يمكننا تطبيق تقريب باديه ذي النقطتين علىو(x){\displaystyle f(x)}توجد طريقة لاستخدام هذا للحصول على حل تقريبي لمعادلة تفاضلية بدقة عالية. [ 9 ] كذلك، بالنسبة للأصفار غير التافهة لدالة زيتا لريمان، يمكن تقدير أول صفر غير تافه بدقة معينة من السلوك التقاربي على المحور الحقيقي. [ 9 ]

تقريب باديه متعدد النقاط

يُعد تقريب باديه متعدد النقاط امتدادًا آخر لتقريب باديه ثنائي النقاط. [ 9 ] تعالج هذه الطريقة نقاط التفرد.x=xج(ج=1،2،3،...،شمال){\displaystyle x=x_{j}(j=1,2,3,\dots ,N)}دالةو(x){\displaystyle f(x)}والتي يجب تقريبها. ضع في اعتبارك الحالات التي يتم فيها التعبير عن نقاط التفرد لدالة ما باستخدام الدليلنج{\displaystyle n_{j}}بواسطة و(x)أج(x-xج)نج،xxج.{\displaystyle f(x)\sim {\frac {A_{j}}{(x-x_{j})^{n_{j}}}},\quad x\to x_{j}.}

إلى جانب تقريب باديه ذي النقطتين، والذي يتضمن معلومات عندx=0،x{\displaystyle x=0,x\to \infty }، هذه الطريقة تقارب لتقليل خاصية التباعد عندxxج{\displaystyle x\sim x_{j}}ونتيجة لذلك، وبما أن معلومات خصوصية الدالة قد تم التقاطها، فإن تقريب الدالةو(x){\displaystyle f(x)}يمكن إجراؤها بدقة أعلى.

أمثلة

sin( x ) [ 14 ]
الخطيئة(x)126714363920x5-236318183x3+x1+44512122x2+601872784x4+12116662240x6{\displaystyle \sin(x)\approx {\frac {{\frac {12671}{4363920}}x^{5}-{\frac {2363}{18183}}x^{3}+x}{1+{\frac {445}{12122}}x^{2}+{\frac {601}{872784}}x^{4}+{\frac {121}{16662240}}x^{6}}}}
exp( x ) [ 15 ]
خبرة(x)1+12x+19x2+172x3+11008x4+130240x51-12x+19x2-172x3+11008x4-130240x5{\displaystyle \exp(x)\approx {\frac {1+{\frac {1}{2}}x+{\frac {1}{9}}x^{2}+{\frac {1}{72}}x^{3}+{\frac {1}{1008}}x^{4}+{\frac {1}{30240}}x^{5}}{1-{\frac {1}{2}}x+{\frac {1}{9}}x^{2}-{\frac {1}{72}}x^{3}+{\frac {1}{1008}}x^{4}-{\frac {1}{30240}}x^{5}}}}
ln(1+ x ) [ 16 ]
ln(1+x)x+12x21+x+16x2{\displaystyle \ln(1+x)\approx {\frac {x+{\frac {1}{2}}x^{2}}{1+x+{\frac {1}{6}}x^{2}}}}
جاكوبي sn( z |3) [ 17 ]
sن(z|3)-9851629283609260z5-5727444726821z3+z1+8594901575607z2-592203556721852z4+625315912977897230z6{\displaystyle \mathrm {sn} (z|3)\approx {\frac {-{\frac {9851629}{283609260}}z^{5}-{\frac {572744}{4726821}}z^{3}+z}{1+{\frac {859490}{1575607}}z^{2}-{\frac {5922035}{56721852}}z^{4}+{\frac {62531591}{2977897230}}z^{6}}}}
بيسيل J 5 ( x )
ج5(x)-10728416000x7+13840x51+1515550x2+14533729600x4+1339358041600x6+2767120301977600x8{\displaystyle J_{5}(x)\approx {\frac {-{\frac {107}{28416000}}x^{7}+{\frac {1}{3840}}x^{5}}{1+{\frac {151}{5550}}x^{2}+{\frac {1453}{3729600}}x^{4}+{\frac {1339}{358041600}}x^{6}+{\frac {2767}{120301977600}}x^{8}}}}
خطأ erf( x )
قطعة أرض(x)215π49140x+3570x3+739x5165x4+1330x2+3276{\displaystyle \operatorname {erf} (x)\approx {\frac {2}{15{\sqrt {\pi }}}}\cdot {\frac {49140x+3570x^{3}+739x^{5}}{165x^{4}+1330x^{2}+3276}}}
فريسنل C ( x )
ج(x)1135990791π4x9-147189744π2x5+8714684160x1749π4x8+523536π2x4+64553216{\displaystyle C(x)\approx {\frac {1}{135}}\cdot {\frac {990791\pi ^{4}x^{9}-147189744\pi ^{2}x^{5}+8714684160x}{1749\pi ^{4}x^{8}+523536\pi ^{2}x^{4}+64553216}}}

انظر أيضاً

مراجع

  1. "Padé Approximant" ، Wolfram MathWorld
  2. النظرية 1 في وين، بيتر (مارس 1966)، "حول تقارب واستقرار خوارزمية إبسيلون"، مجلة SIAM للتحليل العددي ، 3 (1): 91-122 ، Bibcode : 1966SJNA....3...91W ، doi : 10.1137/0703007 ، JSTOR 2949688 
  3. بريزنسكي، سي. (1996)، "خوارزميات الاستقراء وتقريبات باديه"، الرياضيات العددية التطبيقية ، 20 (3): 299-318 ، CiteSeerX 10.1.1.20.9528 ، doi : 10.1016/0168-9274(95)00110-7 
  4. بيني، داريو؛ بان، فيكتور (1994)، حسابات كثيرات الحدود والمصفوفات - المجلد 1. الخوارزميات الأساسية ، سلسلة التقدم في علوم الحاسوب النظرية، بيركهاوزر، المسألة 5.2ب والخوارزمية 5.2 (ص 46)، ISBN 978-0-8176-3786-6
  5. أدلر، جوان (1994)، "توسعات المتسلسلات" ، الحوسبة في الفيزياء ، 8 (3): 287، رمز Bibcode : 1994ComPh...8..287A ، doi : 10.1063/1.168493
  6. بيكر، جي إيه جونيور (2012)، "باديه التقريبي" ، سكولاربيديا ، 7 (6): 9756، رمز Bibcode : 2012SchpJ...7.9756B ، doi : 10.4249/scholarpedia.9756
  7. تشيشولم، جيه إس آر (1973)، "التقريبات النسبية المعرفة من متسلسلات القوى المزدوجة" ، رياضيات الحساب ، 27 (124): 841-848 ، doi : 10.1090/S0025-5718-1973-0382928-6 ، ISSN 0025-5718 
  8. غريفز-موريس، بي آر؛ روبرتس، دي إي (1975)، "حساب تقريبات كانتربري"، اتصالات فيزياء الحاسوب ، 10 (4): 234-244 ، رمز Bibcode : 1975CoPhC..10..234G ، doi : 10.1016/0010-4655(75)90068-5
  9. 1 2 3 4 أويوكا، يوشيكي، مقدمة في طريقة الجمع متعدد النقاط. الرياضيات التطبيقية الحديثة التي تربط هنا بما وراء اللانهاية: من متسلسلة تايلور إلى تطبيق المعادلات التفاضلية.
  10. جونز، ويليام ب.؛ نياستاد، أولاف؛ ثرون، دبليو جيه (1983)، "توسعات باديه ذات النقطتين لعائلة من الدوال التحليلية"، مجلة الرياضيات الحسابية والتطبيقية ، 9 (2): 105-123 ، doi : 10.1016/0377-0427(83)90034-1
  11. سيدي، أفرام (1980)، "بعض جوانب تقريبات باديه ذات النقطتين"، مجلة الرياضيات الحسابية والتطبيقية ، 6 (1): 9-17 ، doi : 10.1016/0771-050X(80)90012-1
  12. ماغنوس، آرني (1982)، حول بنية جدول باديه ذي النقطتين ، سلسلة محاضرات في الرياضيات، المجلد 932، ص 176، doi : 10.1007/BFb0093301 ، ISBN   978-3-540-11567-0
  13. أميندافار، هـ.؛ ريتسي، ج. أ. (1994)، "تقريبات بادي لدوال كثافة الاحتمال"، معاملات IEEE في مجال الفضاء الجوي ، 30 (2): 416-424 ، رمز Bibcode : 1994ITAES..30..416A ، doi : 10.1109/7.272264
  14. "تقريب باديه لـ sin(x)" ، موقع وولفرام ألفا ، تم الاطلاع عليه بتاريخ 16 يناير 2022
  15. "تقريب باديه للدالة الأسية exp(x)" ، موقع وولفرام ألفا ، تم الاطلاع عليه بتاريخ 3 يناير 2024
  16. "تقريب باديه لـ log(1+x)" ، موقع وولفرام ألفا ، تم الاطلاع عليه بتاريخ 16-09-2023
  17. "تقريب باديه لـ sn(x|3)" ، موقع وولفرام ألفا ، تم الاطلاع عليه بتاريخ 16 يناير 2022

الأدب